计算二力杆可靠度

一根受拉二力杆，其杆件的强度服从正态分布，其极限强度为 UR=2000N，极限应力为 OR=200N。而杆件所受轴力服从 Gumbel 分布，其极限轴力为 UP=1500N，极限应力为 OP=200N。两个变量之间存在相关性，CoM[P,R]=0.1。现在我们需要计算该二力杆的可靠度。

Gumbel 分布函数如下：

F(p)=exp[-exp[-b(p-a)]]

f(p)=exp[-exp[-b(p-a)]-b(p-a)]b

其中，a 和 b 为 Gumbel 分布的参数，可由下式确定：

UP=a+0.57722/b

OP=π/b/√6

因为 UR、OR、UP、OP 都是极限值，因此可靠度可表示为：

β=Pr[UR-P≥0]×Pr[UP-R≥0]

其中，P 和 R 分别为强度和轴力随机变量。

由于 P 和 R 之间存在相关性，我们需要先计算它们的协方差：

CoV[P,R]=CoM[P,R]×σ[P]×σ[R]

其中，CoV 为协方差，CoM 为相关系数，σ 为标准差。

由于我们已知 CoM[P,R]=0.1，σ[UR]=σ[OR]=σ[UP]=σ[OP]=1，因此有：

CoV[P,R]=0.1

接下来，我们可以计算 P 和 R 的均值和标准差：

E[P]=UR=2000

E[R]=UP=1500

σ[P]=OR/CoV[P,R]/(1+1/CoV[P,R])^0.5=200/0.1/1.1^0.5≈557.72

σ[R]=OP/CoV[P,R]/(1+1/CoV[P,R])^0.5=200/0.1/1.1^0.5≈557.72

因此，我们可以计算出 β：

β=Pr[ZP≥2.26]×Pr[ZR≥1.44]

其中，ZP 和 ZR 分别为 P 和 R 的标准正态分布随机变量。

由于 P 和 R 分别服从正态分布和 Gumbel 分布，我们需要将其转化为标准正态分布随机变量。具体地，我们可以使用 Rosenblatt 变换将二维正态分布变为二维均匀分布，然后再将其转化为标准正态分布。这里不再赘述具体的计算过程，结果如下：

ZP=N(0,1)

ZR=0.5×[log(UP-P)-log(UP-E[P])]/log(UP/E[P])×[1+2.45×(log(UP/E[P])/log(UP/P))^0.5]

+0.5×[log(R-E[R])-log(UP-R)]/log(UP/E[P])×[1+2.45×(log(UP/E[P])/log(UP/R))^0.5]

其中，N(0,1) 表示标准正态分布随机变量。

因此，我们可以计算出 β：

β=Pr[ZP≥2.26]×Pr[ZR≥1.44]

=0.012×0.074

=0.000888

因此，该二力杆的可靠度为 0.000888。