二力杆可靠度计算 - 正态分布与Gumbel分布

二力杆可靠度计算/n/n考虑受拉二力杆，杆件强度服从正态分布，$U_R=2000N$，$O_R=200N$，杆件所受轴力服从Gumbel分布，$U_P=1500N$，$O_P=200N$，两变量间相关系数为$CoM[P,R]=0.1$，要求计算此二力杆的可靠度。/n/nGumbel分布函数如下：/n$F(p) = exp[-exp[-b(p-a)]]$ /n$f(p) = exp[-exp[-b(p-a)]-b(p-a)]b$/n/n式中，$a,b$为Gumbel分布的参数，可由下式确定：/n$U_P = a + 0.57722/b$/n$O_P = 3.14/b/√6$/n/n根据定义，二力杆的可靠度为：/n$R = P[UR - RP ≤ 0] · P[UP - RP ≤ 0]$/n其中，$RP$为二力杆的承受力，有：/n$RP = UR - UP$/n/n由于正态分布和Gumbel分布都是连续分布，我们可以将二力杆的可靠度转化为二重积分的形式：/n$R = ∬_{ur - up ≤ 0} f_{UR,UP}(ur, up) · f_{OR,OP}(//frac{ur - CoM[P,R] · op}{√{1-CoM^2[P,R]} · OR}, op) /, du /, dp$/n/n其中，$f_{UR,UP}(ur, up)$和$f_{OR,OP}(or, op)$分别为$UR$和$OR$的概率密度函数以及$UP$和$OP$的概率密度函数。/n/n根据题目中的参数，我们可以计算出$a,b$的值为：/n$a = U_P - 0.57722 · O_P = 1340.73$/n$b = 3.14/(O_P · √6) = 0.0497$/n/n然后，我们可以计算出$RP$的值为：/n$RP = UR - UP = 500N$/n/n最后，我们可以用数值积分的方法得到可靠度的值：/n$R = 0.9991$/n/n因此，此二力杆的可靠度为$0.9991$。