[发明专利]一种基于JC算法的热障涂层界面氧化失效可靠性评估方法

摘要

一种基于JC算法的热障涂层界面氧化失效可靠性评估方法，属于热障涂层材料的可靠性分析技术领域。建立热障涂层的失效准则确立其极限状态方程，分析状态方程中各参量的随机统计特性；设定各参量的初始验算点，将功能函数在验算点进行一阶泰勒展开，迭代计算出最优验算点，输出该失效模式的一阶失效概率；当功能函数非线性程度较高时，将功能函数在验算点处进行二阶泰勒展开，计算二阶失效概率；拟合失效概率和各参量的二次函数，计算各参量的敏感性因子。本发明借鉴工程可靠性分析的JC算法，能够简单、快速、定量地评估热障涂层的可靠性，还可根据敏感性因子来分析各参量对热障涂层失效的影响程度，对热障涂层的可靠性评估具有重大意义。

主权项

1.一种基于JC算法的热障涂层界面氧化失效可靠性评估方法，其特征在于，该方法步骤如下：步骤一：根据热障涂层的主要失效模式，建立失效准则；步骤二：根据失效准则建立极限状态方程：Z＝g_X(X₁,X₂,…,X_n)＝0    (1)其中，Z＝g_X(X₁,X₂,…,X_n)为功能函数，X_i(i＝1,2,…,n)为影响热障涂层失效的参量；步骤三：确定各参量X_i(i＝1,2,...n)的随机统计特性，包括其分布函数类型，平均值和标准差；步骤四：假设各参量X_i(i＝1,2,…,n)的设计验算点p^*的初始值为x*=x1*x2*...xn*T,]]>一般可取其平均值：x*=μ1*μ2*...μn*T---(2)]]>其中，为各参量X_i(i＝1,2,…,n)的平均值；步骤五：对非正态分布参量X_i(i＝1,2,…,n)，需要进行当量正态化处理，并用当量正态化X_i'的平均值和标准差替换对应X_i的平均值和标准差：μXi′=xi*-Φ-1[FXi(xi*)]σXi′---(3)]]>其中，为非正态分布参量X_i(i=1,2,…n)的分布函数，为其逆函数，为非正态分布参量X_i(i=1,2,…n)的概率密度函数；步骤六：求出各参量X_i(i＝1,2,…,n)的灵敏度系数cosθ_xi：cosθXi=-∂gX(x*)∂Xiσxi′Σi=1n[∂gX(x*)∂Xiσxi′]2---(5)]]>其中，为功能函数Z＝g_X(X₁,X₂,…,X_n)在验算点x^*处对X_i的一阶偏导数，为当量正态分布参量X_i'的标准差；步骤七：将功能函数Z＝g_X(X₁,X₂,…,X_n)在验算点x^*处进行一阶泰勒级数展开并取至一次项，求出可靠指标β：β=μzLσzL=gX(x1*,x2*,...,xn*)+Σi=1n∂gX(x*)∂Xi(μXi′-xi*)Σi=1n[∂gX(x*)∂Xiσxi′]2---(6)]]>其中，可靠指标β表示在标准正态空间中，原点到极限状态曲线的最短距离，为功能函数在验算点x^*处的一阶泰勒展开式的平均值，为其标准差；步骤八：由上述β和cosθ_xi可得到新的验算点x^*，其坐标为：xi*=μXi+βσXicosθ.---(7)]]>其中，为参量X_i的平均值，β为可靠指标，为参量X_i的标准差，为参量X_i的灵敏度系数；步骤九：以新的重复步骤四～步骤八，直至前后两次的||x^*||＜ε，ε为设置的允许误差精度，则热障涂层的一阶失效概率：p_f＝Φ(-β)    (8)其中，β为可靠指标，Φ(-β)为-β的标准正态分布函数；步骤十：当功能函数的非线性程度较高时，在上述步骤的基础上，计算单位向量α_Y：αY=αY1αY2...αYnT=-▿gY(y*)||▿gY(y*)||---(9)]]>其中，▽g_Y(y^*)为在标准正态空间变量Y空间下，功能函数Z＝g_Y(Y₁,Y₂,...Y_n)在验算点y^*处的一阶偏导数，变量X空间和标准正态变量Y空间的转化关系为αYi=-∂gY(y*)∂Yi/||▿gY(y*)||(i=1,2,...n),]]>||▿gY(y*)||=Σi=1n(∂gY(y*)∂Yi)2;]]>步骤十一：以α_Y作为第n列向量用正交规范化处理技术构造一个正交矩阵H；步骤十二：利用式(10)计算矩阵Q：Q=-▿2gY(y*)||▿gY(y*)||---(10)]]>其中，▽²g_Y(y^*)为在标准正态变量Y空间下，功能函数Z=g_Y(Y₁,Y₂,…,Y_n)在验算点y^*处的二阶偏导数，||▿gY(y*)||=Σi=1n(∂gY(y*)∂Yi)2;]]>步骤十三：利用式(11)计算热障涂层的二阶失效概率p_fQ：pfQ≈Φ(-β)det[I-β(HTQH)n-1]---(11)]]>其中，Φ(-β)为-β的标准正态分布函数，det[Ι-β(H^TQH)_n-1]为矩阵[Ι-β(H^TQH)_n-1]的行列式；步骤十四：假定计算出的失效概率p_F表示为某一参量x的二次函数：p_F＝ax²+bx+c   (12)其中，p_F为步骤九中的p_f或步骤十三中的p_fQ；步骤十五：对(x_i,p_Fi),i＝1,2,3,4,5进行加权最小平方回归，拟合二次函数式(13)，相应的加权因子为(∂1,∂2,∂3,∂4,∂5)=(1,2,4,2,1);]]>其中，(x_i),i＝1,2,3,4,5为该参量取5个不同的值，p_Fi为对应各x_i的失效概率，x₃为敏感性分析参考点；步骤十六：计算该参量x的敏感性因子：wpx=|∂pF∂x|x=xr|×xrpFr×pxr≈|2axr+b|×xrpFr×p---(13)]]>其中，x_r对应敏感性分析参考点x₃，a为二次函数式(12)的二次项系数，b为二次函数式(12)的一次项系数，为x偏离x_r的概率。