一种基于配点型算法的随机参数结构可靠性评估方法技术

本发明专利技术提供一种基于配点型算法的随机参数结构可靠性评估方法，涉及工程结构可靠性分析技术领域。步骤包括：对于具有多个独立随机变量的工程结构设计问题，建立包含随机参数的结构功能函数；建立结构功能函数积分的等效弱形式，得到真实结构功能函数的一元替代算式；确定结构功能函数统计量，包括均值、方差及高阶原点矩和中心矩的近似表达式；在随机变量定义域内进行配点，确定配点方案；确定结构失效概率，用来评估结构可靠性。本发明专利技术综合考虑了工程结构可靠性问题在计算精度和计算效率上的平衡，基于结构功能函数的一元等效积分弱形式和配点型方法给出结构的失效概率，从而减缓了随机响应的正态性假设所带来的模型误差。

【技术实现步骤摘要】
一种基于配点型算法的随机参数结构可靠性评估方法
本专利技术涉及工程结构可靠性分析
，尤其涉及一种基于配点型算法的随机参数结构可靠性评估方法。
技术介绍
工程结构设计需要保证结构在规定的时间内和给定的载荷下完成预定的功能。然而，在设计过程和实际操作环境中往往存在大量的不确定性，大致可分为物理不确定性、统计不确定性和模型不确定性。物理不确定性由载荷、材料、几何和边界条件等参数的变异性所致。统计不确定性是由样本数量的不足而对分布参数产生的估计偏差所造成的，也经常是由于真实分布形式与预估分布形式存在偏差所造成的。而模型不确定性是指所使用的数学、物理模型不能完全精确地反映问题的本质，只是对客观情况的一种近似描述。为保证结构的安全性能，需要对这些不确定性带来的负面效果进行合理地量化和管理。在具有充足的样本信息时，对工程结构设计和操作环境中的不确定性使用概率模型进行描述是一种合理的方法，可使用可靠度或失效概率的概念量化结构完成预定功能的概率。在工程结构的可靠性分析中，由于各随机变量的联合概率分布函数通常极难获得，所以精确的全概率积分方法很难应用于实际问题之中。作为工程中一种现实的选择，人们转向研究只使用随机变量前两阶矩的一次二阶矩方法，其中包括中心点法、验算点法等。这类方法将非线性功能函数在随机变量域内的某点处展开为线性函数得到，计算过程简单但具有较低的计算精度。无论是中心点法还是验算点法，都没有正面解决响应函数Z的分布问题，而只满足于在求得可靠度指标β值后，通过假定Z的分布(通常为正态)来解决β与失效概率的转换问题。这在一定程度上引入了模型不确定性。为改善一次二阶矩方法的计算精度，国内外广泛开展了二次二阶矩、二次四阶矩方法的研究。虽然随着展开的阶数和所使用统计量的增多，计算精度稍显增加，但这些方法都是基于Taylor展开的，很重要的环节是需要计算响应函数关于各随机变量的灵敏度，工程问题中往往以差商运算代替微商运算，而这个过程也将引入难以估量的误差。鉴于基于Taylor展开方法的精度限制，统计类可靠性分析方法也得到了长足的发展，尤其是在处理需要高可靠度指标的工程问题中得以应用。MonteCarlo仿真方法的理论基础是概率论中的大数定理，具有不受应用范围限制的优点，但其缺点是需要消耗极大的计算代价，虽然已出现了改进的多种抽样方法，但其在实际工程计算中仍较少采用，一般多用于检验一些新提出方法的精度。所以，在目前工程结构可靠度分析方法中，缺少一类具有较高计算精度，但计算量又可以承受的方法。
技术实现思路
针对现有技术的缺陷，本专利技术提供一种基于配点型算法的随机参数结构可靠性评估方法，综合考虑了工程结构可靠性问题在计算精度和计算效率上的平衡，基于结构功能函数的一元等效积分弱形式和配点型方法给出结构的失效概率，从而减缓了随机响应的正态性假设所带来的模型误差。一种基于配点型算法的随机参数结构可靠性评估方法，包括以下步骤：步骤1：对于包含多个独立随机变量的工程结构设计问题，建立参数化结构功能函数，定义其结构功能函数为X＝g(h)＝g(h1，…hN)(1)其中，随机变量i＝1，…N，N为正整数，和分别为第i个随机变量hi的下界值和上界值，Ω为一超长方体形成的凸域，并记域Ω的中心点为半径为步骤2：确定结构功能函数的等效积分弱形式；根据结构功能函数g(h)在中心点hc处具有收敛的Taylor展开式近似计算函数g(h)在凸域Ω上的积分，并记为I[g(h)]；取g(h)的一元近似函数为式(2)；将g(h)的一元近似函数同样进行Taylor展开并积分，记为使用作为I[g(h)]的近似函数，得到具有4阶小量的残差估计；步骤3：确定结构功能函数统计量；基于积分运算和求和运算的可交换性，得到结构功能函数g(h)的期望μg和方差一元表示形式分别如式(3)和式(4)所示；其中，为一元分解函数的均值，表示第i个一元分解函数，为第i个一元分解函数的方差，gc为结构功能函数g(h)在点h＝hc处的值；结构功能函数g(h)的m阶原点矩和中心矩分别如式(5)和式(6)所示；其中，表示第i个一元分解函数的第m阶原点矩，表示组合方法，C表示数学排列组合中的组合符号，表示第i个一元分解函数的第k阶中心矩；步骤4：在随机变量定义域内进行配点，得到式(7)所示的均值的数值积分格式；其中，f(hi)为随机变量hi在域Γ内的概率密度分布函数，为第i个一元分解函数在随机变量hi所在区间内的第j个配点，j＝1，…，NP，NP为配点的个数，为相应的权重因子；权重因子与配点的数目和位置是相关的，需根据具体的分布进行确定；步骤5：为避免使用功能函数的正态性假设所带来的模型误差，考虑基于非高斯随机变量X的统计量，根据一维概率密度函数的渐进展式以近似得到结构功能函数的概率密度分布函数：其中，z所对应的随机变量Z由随机变量X标准化得到，记Z＝(X-μX)/σX，这里μX为X的均值，而σX为X的标准差，p(z)为高斯概率密度函数，cn为待定系数；使用所得到的原点矩和中心距计算概率密度函数p*(z)，得到结构的失效概率Pf的估计式如式(9)所示；根据所得到的结构失效概率Pf和工程结构的失效标准来评估结构的可靠性。由上述技术方案可知，本专利技术的有益效果在于：本专利技术提供的基于配点型算法的随机参数结构可靠性评估方法，可达到4阶计算精度，并且计算量随着随机变量数目的增加只是线性增长，有效提高了计算精度和计算效率。根据非高斯随机变量概率密度函数的Edgeworth渐进展式，得到结构响应的逼近概率密度函数，对其直接进行积分得到结构的失效概率，从而减缓了正态性假设所带来的模型不确定性。附图说明图1为本专利技术实施例提供的基于配点型算法的随机参数结构可靠性评估方法流程图。具体实施方式下面结合附图和实施例，对本专利技术的具体实施方式作进一步详细描述。以下实施例用于说明本专利技术，但不用来限制本专利技术的范围。一种基于配点型算法的随机参数结构可靠性评估方法，如图1所示，本实施例的方法如下所述。步骤1：建立参数化结构功能函数。工程结构中的随机变量包括材料属性、环境参数、模型误差、认知能力等，对于包含多个独立随机变量的工程结构设计问题，其结构功能函数定义为X＝g(h)＝g(h1，…hN)(1)式中，i＝1，…N，N为正整数，Ω为一超长方体形成的凸域。并记域Ω的中点为半径为步骤2：确定结构功能函数的等效积分弱形式。假设功能函数g(h)在中心点hc处具有收敛的Taylor展开式式中，算子根据式(10)，可近似计算函数g(h)在凸域Ω上的积分引入下面的关系式式中，mi为一非负整数。值得注意的是，当且仅当集合{mi，i＝1，2，…，N}中至少有一个元素为奇数时，式(12)成立。将式(10)代入式(11)，并使用式(12)进行简化，可得式(13)。另一方面，取函数将函数式(2)进行Taylor展开并进行积分，利用关系式(12)，得到使用作为I[g(h)]的近似函数，可得残差估计式(15)由式(15)可知，所得到的残差为4阶小量。工程问题一般将响应函数展开至1-2阶，即可获得较为满意的数值精度，所以，函数式(2)可在上述的积分形式下近似等效于结构功能函数g(h)。为方便，记则式(2)可重写为其中，称为第i个一元分解函数。结构功能函数积分本文档来自技高网...

【技术保护点】
一种基于配点型算法的随机参数结构可靠性评估方法，其特征在于：该方法包括以下步骤：步骤1：对于包含多个独立随机变量的工程结构设计问题，建立参数化结构功能函数，定义其结构功能函数为X＝g(h)＝g(h

【技术特征摘要】
1.一种基于配点型算法的随机参数结构可靠性评估方法，其特征在于：该方法包括以下步骤：步骤1：对于包含多个独立随机变量的工程结构设计问题，建立参数化结构功能函数，定义其结构功能函数为X＝g(h)＝g(h1,…hN)(1)其中，随机变量N为正整数，和分别为第i个随机变量hi的下界值和上界值，Ω为一超长方体形成的凸域，并记域Ω的中心点为半径为步骤2：确定结构功能函数的等效积分弱形式；根据结构功能函数g(h)在中心点hc处具有收敛的Taylor展开式近似计算函数g(h)在凸域Ω上的积分，并记为I[g(h)]；取g(h)的一元近似函数为式(2)；将g(h)的一元近似函数同样进行Taylor展开并积分，记为使用作为I[g(h)]的近似函数，得到具有4阶小量的残差估计；步骤3：确定结构功能函数统计量；基于积分运算和求和运算的可交换性，得到结构功能函数g(h)的期望μg和方差一元表示形式分别如式(3)和式(4)所示；其中，为一元分解函数的均值，表示第i个一元分解函数，为第i个一元分解函数的方差，gc为结构功能函数g(h)在点h＝hc处的值；结构功能函数g(h)的m阶原点矩和中心矩分别如式(5)和式(6)所示；

【专利技术属性】
技术研发人员：祁武超，王晓军，
申请(专利权)人：沈阳航空航天大学，北京航空航天大学，
类型：发明
国别省市：辽宁,21