第三章结构可靠性分析方法.ppt

单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,3.1,一次可靠度分析法,一次可靠度分析法,(First Order Reliability Method,FORM),计算结构构件可靠度的基本思路是：首先将结构构件功能函数,Z=g(X,l,，,X,2,，,，,X,n,),展开成,Taylor,级数，忽略高阶项，仅保留线性项，再利用基本随机变量,X=(X,l,X,2,X,n,),的一阶矩、二阶矩求取,Z,的均值,z,与标准差,z,，从而确定结构构件可靠指标。根据功能函数线性化点的取法不同以及是否考虑基本随机变量的分布类型，一次可靠度分析法分为：,均值一次二阶矩法,(,中心点法,)改进的一次二阶矩法,(,验算点法,),和,JC,法,。,3.1,一次可靠度分析法,泰勒,(Taylor),中值定理,(,一元,):,如果函数,f,(,x,),在含有,x,0,的某个开区间,(a,，,b),内具有直到,(n+1),阶导数，则当,x,在,(a,，,b),时，,f,(,x,),可表示为,(,x,-,x,0,),的一个,n,次多项式与一个余项,R,n,(,x,),之和：,一元函数,3.1,一次可靠度分析法,泰勒公式,(,二元,):,设,z=,f,(,x,y,),在点,(,x,0,y,0,),的某一邻域内连续且有直到,(n+1),阶导数，有：,一次可靠度分析常取前面两项。即线性项。,可以推广至有,n,元情况,2.2.1,均值一次二阶矩法,2.2,一次可靠度分析法,1,、均值一次二阶矩法,(,中心点法,),当功能函数包含有多个相互独立的,正态随机变量,X=(X,l,X,2,X,n,),，状态函数为：,Z=g(X,l,，,X,2,，,，,X,n,),。,随机变量标准差与其函数标准差的近似表达。,3,.1.1,均值一次二阶矩法,3.1,一次可靠度分析法,1,、均值一次二阶矩法,(,中心点法,),计算步骤：,当功能函数包含有多个相互独立的,正态随机变量,X=(X,l,X,2,X,n,),，状态函数为：,Z=g(X,l,，,X,2,，,，,X,n,),。,(1),用各随机变量的均值代入功能函数，得出功能函数的均值,Z,；,(2),求功能函数,的标准差,Z,；,(3),求,和,P,f,。,将随机变量的均值代入,例题：可靠度分析,解：结构基本变量为,f,y,和,d,，荷载极限状态方程：,右图所示圆截面直杆，承受拉力,P=100KN,，已知材料的强度设计值,f,y,的均值,fy,=290MPa,，标准差,fy,=25MPa,，杆直径,d,的均值,d,=30mm,，标准差,d,=3mm,，求此杆的可靠指标。,例题：可靠度分析,解：结构基本变量为,f,y,和,d,，应力极限状态方程：,右图所示圆截面直杆，承受拉力,P=100KN,，已知材料的强度设计值,f,y,的均值,fy,=290MPa,，标准差,fy,=25MPa,，杆直径,d,的均值,d,=30mm,，标准差,d,=3mm,，求此杆的可靠指标。,均值一次二阶矩法简单，使用方便。但存在缺陷：其一是对同样条件下的同一结构，若采用不同的功能函数来描述结构的同一功能要求，将得出不同的,值；其二是选取基本随机变量均值点作为功能函数的线性化点来求取,将产生较大误差。,3,.1.1,均值一次二阶矩法,3.1,一次可靠度分析法,2,、可靠指标,的几何意义,设,Z=g(X,l,，,X,2,，,，,X,n,),是线性函数，极限状态方程为：,D(AX+BY)=A,2,D(X)+B,2,D(Y),3,.1.1,均值一次二阶矩法,3.1一次可靠度分析法,2,、可靠指标,的几何意义,将,X,l,，,X,2,，,，,X,n,作,标准化,变换：,U,i,在,空间的均值为零，标准差为,1,。有：,原结构极限状态方程：,在,空间极限状态方程：,3,.1.1,均值一次二阶矩法,3.1,一次可靠度分析法,2,、可靠指标,的几何意义,该方程表示,U,空间中的一个,超平面,。由解析几何知识可知，在,U,空间中坐标原点,(,即中心点,M),到此极限状态超平面的距离为：,在,空间极限状态方程：,两维情况：,Z=R-S,3,.1.1,均值一次二阶矩法,3.1,一次可靠度分析法,2,、可靠指标,的几何意义,上式说明了可靠指标,的几何意义：指在经过标准化变换得到的,U,空间中，从坐标原点,(,即中心点,O),到相应的极限状态超平面的距离。当结构的功能函数为非线性函数时，可以得出相同的结论。中心点在可靠区内，它离开极限状态超平面越远，表明结果越可靠。,3,.2.1,改进的一次二阶矩法,3.2,一次可靠度分析法,-(,改进一次二阶矩法,),针对均值一次二阶矩法将功能函数线性化点取作基本随机变量均值点带来的问题，改进的一次二阶矩法将功能函数线性化点取在设计验算点，从而提高了计算,的精度，并保证了对同一结构问题,的唯一性。改进的一次二阶矩法也称为验算点法。,当极限状态方程中包含有多个,相互独立,的,正态随机变量,X=(X,l,X,2,X,n,),假设方程为：,Z=g(X,l,，,X,2,，,，,X,n,)=0,，则此超曲面,Z=0,上距离中心点,M=(,X1,，,X2,，,，,Xn,),最近的点,P,=(,x,1,，,x,2,，,，,x,n,),为设计验算点，简称验算点。显然，,x,i,(I=1,，,2,，,，,n),满足极限状态方程：,3,.2.1,改进的一次二阶矩法,3.2,一次可靠度分析法,偏导数在验算点的值,3,.2.1,改进的一次二阶矩法,3.2,一次可靠度分析法,结构可靠指标为：,引入灵敏系数,i,灵敏系数,i,实际上是各基本变量的不定性对可靠度影响的“权”，有：,于是有：,上式展开：,3,.2.1,改进的一次二阶矩法,3.2,一次可靠度分析法,上式展开：,g(,x,i,*,)0,，必有：,有,x,1,，,x,2,，,，,x,n,和,的方程组,(n+1),：,解方程组可得验算点,P*=(,x,1,，,x,2,，,，,x,n,),和,值。但这仅是理论上的结论，实际上求解此方程组是相当困难的。通常在给定,X,i,的统计参数,Xi,，,Xi,后，进行迭代运算，求出,x,i,*,和,的近似值。,3,.2.1,改进的一次二阶矩法,3.2,一次可靠度分析法,改进的一次二阶矩法迭代步骤：,(,第,1,步：,假定,=,2.0),；,第,2,步：,设验算点为,x,i,*,，,i=1,2,n,，第一步取基本变量平均值,x,i,*Xi,；,第,3,步：,计算非正态随机变量的当量正态分布均值和标准差；,第,4,步：,计算偏导数 在验算点,x,i,*,的值；,第,5,步：,计算灵敏系数,i,在验算点,x,i,*,的值；,第,6,步：,计算验算点,x,i,*,新值,重复步骤,36,，直到灵敏系数,i,收敛,(,前后两次的之差的绝对值小于,0.005),。一旦,i,收敛，就把,作为未知参数。如果不检验,的收敛性，则删除第,1,步。,第,7,步：,将,视为未知数，把,x,i,*,值代入极限状态方程，解方程得出,值，并计算设计验算点,x,i,*,的新值。,第,8,步：,重复步骤,27,，直至前后两次算出的,值之差小于允许误差一般,0.01),代入当量标准差。,代入当量均值。,3,.2.2 JC,法,3.2,一次可靠度分析法,-JC,法,1,、当量正态化,一次二阶矩法适用于结构功能函数所含随机变量为,独立,、,正态变量,情况。但在可靠性分析中，极限状态方程常常包含非正态分布的随机变量。,JC,法的基本思路是：,对非正态变量当量正态化，将其转换为等效正态随机变量，即可利用一次二阶矩法求结构可靠指标。,“当量正态化”的条件是：,(1),在设计验算点,x,i,处，非正态变量,X,i,(,其均值为,Xi,，标准差为,Xi,),的分布函数值,F,Xi,(,x,i,),与当量正态变量,X,i,(,其均值为,Xi,，标准差为,Xi,),的分布函数值,F,xi,(,x,i,),相等；,(2),在设计验算点,x,i,处，非正态变量,X,i,的概率密度函数值,f,Xi,(,x,i,),与当量正态变量,X,i,的概率密度函数值,f,Xi,(,x,i,),相等。,根据以上两个条件可以计算出当量正态分布的均值和标准差。,3,.3.1 JC,法,3.3,一次可靠度分析法,1,、当量正态化,非正态分布函数,当量正态分布函数,第一个条件：分布函数相等。,3,.3.1 JC,法,3.3,一次可靠度分析法,1,、当量正态化,非正态分布密度函数,当量正态分布密度函数,第二个条件：分布密度函数相等。,3,.3.1 JC,法,3.3,一次可靠度分析法,1,、当量正态化,非正态分布密度函数,当量正态分布密度函数,第二个条件：分布密度函数相等。,3,.3.1 JC,法,3.3,一次可靠度分析法,1,、当量正态化,当量正态化分析步骤：,(1),根据分布密度函数相等得出当量正态分布的标准差。,(2),根据分布函数相等得出当量正态分布的均值。,将验算点,x,i,代入非正态分布函数，得一计算值，反查正态分布表，将该值代入标准正态分布密度函数得分子值。将验算点,x,i,代入非正态分布密度函数，得分母值。,(2),在极限状态方程中，非正态随机变量的当量正态的均值和标准差求得后，即将问题化为正态变量的情形。,前提条件是必须得出非正态分布的分布函数和分布密度函数。,3,.3.1 JC,法,3.3,一次可靠度分析法,1,、当量正态化,对于对数正态分布变量,X,i,，其分布函数和分布密度函数为：,分布密度函数：,在正态分布公式中令,z=(t-)/,，可将随机变量,X,标准化，标准化后的随机变量,z,服从标准正态分布。,对数正态分布的密度函数为：,正态分布的密度函数为：,3,.3.1 JC,法,3.3,一次可靠度分析法,1,、当量正态化,对于对数正态分布的分布函数和分布密度函数：,由公式，当量标准差为：,由公式，当量均值为：,正态与对数正态分布转换,P29,页。,关键的两个公式：,例题：当量正态化,解：从前面的分析中得到如下公式：,设某随机变量服从对数正态分布，其样本均值和方差分别为：,试求其当量正态分布的均值,Xi,和标准差,Xi,。,验算点,x,i,即为均值点。题目已知的是样本的均值和方差，但公式中为样本取对数后的均值,X,和方差,X,。由,P29,页：,解方程组，可得均值,X,和标准差,X,。,对数当量正态化公式,Xi,为样本的均值，,X,为样本的变异系数,x,i,为验算点，,Xi,为当量正态分布的均值，,Xi,为当量正态分布的标准差。,例题：改进一次二阶矩法,解,：,1,第一次迭代,已知极限状态方程,Z=g(B,W)=BW-1140=0,，随机变量,B,服从对数正态分布，平均值,B,=38,，标准差,B,=3.8,；随机变量,W,服从正态分布，平均值,W,=54,，标准差,W,=2.7,。求可靠指标,值。,(2),取验算点,x,1,=b,=38,，,x,2,=W,=54(,为均值点,),。,(3),对于,B,随机变量当量正态化，,B,=3.8/38=0.1,。,对数正态分布当量正态化公式,(1),设,=3,。,(4),例题：改进一次二阶矩法,(6),计算设计验算点,x,i,的新值：,(5),计算灵敏系数,i,在设计验算点的值,x,i,。,重复步骤,(3)(6),，直到灵敏系数,i,收敛,(,前后两次的,i,之差的绝对值小于允许误差,(,一般,0.005),经过三次迭代，,B,和,W,的值收敛为,0.8800674,和,0.4748486,。,此时计算的,B,N,=36.698103,，标准差,B,N,=2.8942793,；,例题：改进一次二阶矩法,(7),将,视为未知数，把 代入极限状态方程，解方程计算,值，并计算设计验算点,x,i,的新值。,(8),重复,(2)(7),步，直到前后两次计算得到的,相差小于,0.01,。,3.,26,2,-184.36,+841.7=0,两个解：一个为,5.002,，另一个为,51.543,，取第一个解,例题：,对于这种比较复杂的迭代计算，一般习惯于将计算结果列成一个表格。,3,.3.2 JC,法,3.3,一次可靠度分析法,2,、相关随机变量的处理,对于随机变量相关的情形，需将它们先变换为相互独立的变量，然后再运用一次二阶矩法或改进的一次二阶矩法求可靠指标和失效概率。,将相关变量变换为不相关变量的方法是：通过正交变换。,(1),相关正态随机变量的情况,假设相关正态基本随机变量,X,i,=(X,1,X,2,X,n,),，其协方差矩阵为：,3,.3.2 JC,法,3.3,一次可靠度分析法,2,、相关随机变量的处理,相关正态基本随机变量,X,i,=(X,1,X,2,X,n,),，其协方差矩阵为：,如果,X,i,=(X,1,X,2,X,n,),是不相关正态基本随机变量，其协方差矩阵为：,3,.3.2 JC,法,3.3,一次可靠度分析法,2,、相关随机变量的处理,若,X,i,=(X,1,X,2,X,n,),转换为不相关随机变量,Y,i,=(Y,1,Y,2,Y,n,),，令：,将,X,值代入极限状态方程,Z=g(X,1,X,2,X,n,)=0,，得,:,Z=,f,(Y,1,Y,2,Y,n,)=0,式中：,A,是正交矩阵，表示,D=A,T,C,X,A,D,矩阵为特征向量构成的矩阵，其对角上的值为对应不相关随机变量的方差。,D,即为不相关随机变量的协方差矩阵。,Matlab,求解式为：,V,D=eig(C),C=3 6 1;6 4 0;1 0 2,A,D=eig(C),3,.3.2 JC,法,3.3,一次可靠度分析法,2,、相关随机变量的处理,由,X,iX,1,X,2,X,n,),相关随机变量转换为不相关随机变量,Y,i,=(Y,1,Y,2,Y,n,),的方法与步骤：,(1),根据相关随机变量,X,i,的协方差矩阵,C,X,，求得特征值矩阵,D(,即不相关随机变量办方差矩阵,),和正交矩阵,A,，由,D,可求得不相关随机变量的方差。,由 可得不相关随机变量的均值。,(2),由,X=(A,T,),-1,Y,，可得出,X,i,，将其代入极限状态方程,Z=g(X,1,X,2,X,n,)=0,，得,:Z=,f,(Y,1,Y,2,Y,n,)=0,Z,即为由不相关随机变量构成的极限状态方程。,例题：,JC,法,有一极限状态方程为：,Z=50+5X,1,-4X,2,-2X,3,+X,4,，四个变量均为正态分布，,X,1,，,X,2,，,X,3,为相关变量且协方差矩阵为：,均值 标准差,X,1,15 3.5,X,2,10 3.0,X,3,15 2.5,X,4,5 1.5,协方差矩阵：,X,1,X,2,X,3,X,1,12.25 -5 2,X,2,-5 9 -4,X,3,2 -4 6.25,求系统可靠性指标：,(1),假设没有相关影响；,(2),考虑相关影响。,解：,(1),假设没有相关影响：用,y,i,=(X,i,-,Xi,)/,Xi,进行标准化处理。,Z=50+5(,y,1,1,+,1,)-4(,y,2,2,+,2,)-2(,y,3,3,3,)+(,y,4,4,+,4,),=60+17.5,y,1,-12,y,2,-5,y,3,+1.5,y,4,可靠指标,的迭代过程为：,y,1,y,2,y,3,y,4,Z,1 0 0 0 0 0 60,2 -2.20 1.51 0.63 -0.19 2.75 0.0,3 -2.2 1.51 0.63 -0.19 2.75 0.0,例题：,JC,法,解：,(2),考虑相关影响,由协方差矩阵可得,A,和,D,协方差矩阵：,X,1,X,2,X,3,X,1,12.25 -5 2,X,2,-5 9 -4,X,3,2 -4 6.25,C=12.25 5 2;,-5 9 4;,2 -4 6.25,A,Deig(C),D=3.1233 0 0,0 7.0280 0,0 0 17.3487,A=0.2068 0.6596 -0.7226,0.6646 0.4473 0.5985,0.7180 -0.6040 -0.3459,EX=15;10;15,EY=A*EX,EY=20.5178;5.3071;-10.0427,将,Y,标准化,y,i,=(Y,i,-,Yi,)/,Yi,:,例题：,JC,法,将,X,代入至极限状态方程：,得,Z=50+5X,1,-4X,2,-2X,3,+X,4,=60+21.96y,1,+7.072y,2,-5.33y,3,+1.5y,4,Z,极限状态方程中,y,i,各变量不相关，且得出了其均值与方差，可用均值一次二阶矩法求解可靠指标。,可靠指标,的迭代过程为：,y,1,y,2,y,3,y,4,Z,1 0 0 0 0 0 60,2 -2.34 0.75 0.57 -0.16 2.53 0.0,3 -2.34 0.75 0.57 -0.16 2.53 0.0,可以看出，考虑相关影响时，可靠指标为,2.53,，不考虑相关影响时，可靠指标增加，为,2.75,3.4,混沌优化可靠度分析法,3.4.1,混沌优化,混沌优化是利用,Logistic,迭代方程,X,n+1,=uX,n,(1-X,n,),遍历,(0,1),区间的混沌特性，用来搜索最优解的方法。,Logistic,迭代方程中,u=4,时，系统进入完全混沌状态。,Matlab,模拟程序,x1=0.8,u=4,for i=1:100,x1=u*x1*(1-x1),y(i)=x1,end,plot(1:100,y,*),hold on,plot(1:100,y,-),3.4,混沌优化可靠度分析法,3.4.1,混沌优化,迭代方程,X,n+1,=4X,n,(1-X,n,),遍历,(0,1),区间的混沌特性，用来搜索最优解的方法。,如求方程：,f1=x,2,-32X+10,的极小值。,步骤：,1),赋初值,wmin=10000;2),迭代初值,x1=0.82,；,3),混沌搜索区间,10,20,，,xx1=a1+(a2-a1)*x1,；,4),搜索最优解,200,次；当计算,f1wmin,侧放弃记录；直到搜索达到,200,次。,%,搜索最优解,wmin=1000000,x1=0.82,a1=10,a2=20,for ii=1:200,x1=x1*4*(1-x1),xx1=a1+(a2-a1)*x1,f1=xx1.2-32*xx1+10,if f1wmin,rx=xx1,wmin=f1,end,end,%,画图,for i=1:160,xx2=i*0.2,ff1(i)=xx2.2-32*xx2+10,j(i)=xx2,end,plot(j,ff1,-),设功能函数,Z=R-S,，,R,和,S,均为正态分布R,均值为,3,，方差为,1,，,S,均值为,2,，方差为,1,。,值的几何意义是中心点,M,至,P*,点的最近距离。,2,2,3,S,R,1,P,*,3.4,混沌优化可靠度分析法,3.4.2,混沌优化计算结构可靠度,%,搜索最优解,bt=10000,x1=0.512,a1=0,a2=4,for ii=1:200,x1=x1*4*(1-x1),xx1=a1+(a2-a1)*x1,f1=sqrt(xx1-2)2+(xx1-3)2),if f1bt,rx=xx1,bt=f1,end,end,M,计算方法与步骤：,(1),取初始值,wmin,，迭代,x0,值,(,有,m,随机变量取,m-1,个初始值,),，搜索范围，迭代次数,N,；,(2),设,Z=g(X,1,X,2,X,m,)=0,，有,m,个随机变量，将方程变为,X,m,=g(X,1,X,2,X,m-1,),，确保搜索在超平面上进行。,(3),进行混沌迭代，,V,n+1,i,=4*V,n,i,(1-V,n,i,),，,i=1,2,(m-1),(4),将搜索值映射到搜索范围，,VV,i,=a1+(a2-a1)*V,n+1,i,；,i=1,2,(m-1),(5),根据,VV,i,值，计算,X,m,=g(X,1,X,2,X,m-1,),，得到曲平面上的搜索点,P*(VV,1,，,VV,1,，,，,VV,m-1,，,X,m,),(6),计算均值点,M,点至,P*,点的距离，得,wz,，如果,wzwmin,，放弃,wz,，转步,(3),，继续迭代；,(7),如果迭代达到规定的次数,N,，,wmin,即为,M,点至,P*,点的最小距离，即,=wmin,。,3.4,混沌优化可靠度分析法,3.4.2,混沌优化计算结构可靠度,例题：混沌优化可靠度分析方法,有一极限状态方程为：,Z=50+5X,1,-4X,2,-2X,3,+X,4,，四个变量均为正态分布，,X,1,，,X,2,，,X,3,为不相关随机变量，求系统可靠性指标。,均值 标准差,X,1,15 3.5,X,2,10 3.0,X,3,15 2.5,X,4,5 1.5,解：用,y,i,=(X,i,-,Xi,)/,Xi,进行标准化处理。,Z=50+5(,y,1,1,+,1,)-4(,y,2,2,+,2,)-2(,y,3,3,+,3,)+(,y,4,4,+,4,),=60+17.5,y,1,-12,y,2,-5,y,3,+1.5,y,4,则中心点座标为,(0,0,0,0),，搜索范围,(-3,2),，计算程序如右侧：,可靠指标,的迭代计算结果为：,y,1,y,2,y,3,y,4,Z,1 0 0 0 0 0 60,2 -2.20 1.51 0.63 -0.19 2.75 0.0,3 -2.2 1.51 0.63 -0.19 2.75 0.0,wmin=100,x1=0.512x2=0.485;x3=0.31,a1=-3;b1=2;,for ii=1:1000,x1=x1*4*(1-x1),xx1=a1+(b1-a1)*x1,x2=4*x2*(1-x2),xx2=a1+(b1-a1)*x2,x3=4*x3*(1-x3),xx3=a1+(b1-a1)*x3,xx4=-(60+17.5*xx1-12*xx2-5*xx3)/1.5,f1=sqrt(xx12+xx22+xx32+xx42),if f1wmin,rx=xx1 xx2 xx3 xx4,wmin=f1,end,end,z=60+17.5*xx1-12*xx2-5*xx3+1.5*xx4,蒙特卡洛,(Monte Carl),法又称随机模拟法或统计试验法，它源于第二次世界大战期间，,Von Neuman,和,Ulam,在对裂变物质的中子随机扩散进行直接模拟中，以世界闻名的赌城蒙特卡洛作为该研究的秘密称呼而得名。,设功能函数,Z=g(X,1,X,2,X,n,),，式中：,X,i,(i=1,2,n),为具有任意分布的随机变量。对,X,i,进行,N,次抽样，得,N,组,X,i,j,(i=1,2,n,；,j=1,2,N),随机变量样本。将第,j,组的,X,i,j,(i=1,2,n),的值代入功能函数,Z,j,=g(X,1,X,2,X,n,),，如果,N,个,Z,j,值中有,N,f,个,Z,j,0,，则结构失效的概率为：,3.5,蒙特卡洛模拟,3.5.1,蒙特卡洛模拟的基本原理,当,N,时，根据伯努利大数定理，,P,f,=N,f,/N,。,与其它可靠性分析方法相比，蒙特卡罗模拟应重点解决,两个基本问题,：一是模拟的精度和效率，这主要涉及随机抽样数,N,，二是任意分布随机变量的随机抽样方法。,由于随机抽样试验是概率为,p,的伯努利试验，所以估计值,P,f,的期望为：,3.5,蒙特卡洛模拟,3.5.2,蒙特卡洛模拟的精度和效率,K,/2,为标准正态分布的百分位值，,K,/2,与,的关系为：,K,/2,0.01 2.5758,0.02 2.3263,0.05 1.96,P,f,的方差为：,P,f,的标准差为是伯努利试验的方差。,根据概率统计检验，显著水平为,的置信区间为,3.5,蒙特卡洛模拟,3.5.2,蒙特卡洛模拟的精度和效率,由,误差,为：,需要的试验次数为,当,p,较小时，所需的试验次数,N,较大，随着,N,的增大，误差,减小并趋于收敛。当,p=0.001,，,=20,，,K,/2,=1.96,2,时,当,p=0.1,，,=20,，,K,/2,=1.96,2,时,3.5,蒙特卡洛模拟,3.5.3,蒙特卡洛模拟的随机变量抽样方法,随机变量的抽样是通过产生随机数的方法来完成。通常要分两步进行：,(1),产生在区间,0,1,上的均匀分布随机数；,(2),变换成给定分布的随机数。,人们已提出了各种数学方法产生均匀分布随机数，例如取中法、加同余法、乘同余法、混合同余法和组合同余法。,1,、混合同余法产生伪随机数的方法为：,3.5,蒙特卡洛模拟,3.5.3,蒙特卡洛模拟的随机变量抽样方法,1,、混合同余法产生伪随机数的方法为：,如果,M=16,a=5,c=2,x,0,=7,i x,i,ax,i,+c mod(ax,i,+c,16)u,i,1 7 38 6 6/16,2 6 33 1 1/16,3 1 8 8 8/16,4 8 43 11 11/16,.,要求,c,与,M,互为质数；若,M,能被,4,整除，则,a-1,也能被,4,整除。,x1=7;a=5;c=3;m=16,for i=140,x1=mod(x1*a+c),m),xx(i)=x1,jj(i)=i,end,plot(jj,xx,-),hold on,plot(jj,xx,*),3.5,蒙特卡洛模拟,3.5.3,蒙特卡洛模拟的随机变量抽样方法,2,、正态分布随机数产生,如果,M=16,a=5,c=3,x,0,=7,x1=7;a=5;c=3;m=16,for i=1:40,x1=mod(x1*a+c),m),if x11,x1=7,end,xx(i)=x1,jj(i)=i,end,u=xx/16,for i=1:2:30,x(i)=(sqrt(-2.0*log(u(i)*cos(2*pi*u(i+1),x(i+1)sqrt(-2.0*log(u(i)*sin(2*pi*u(i+1),j(i)=i,j(i+1)=i+1,end,plot(j,x,-),hold on,plot(j,x,*),产生,40,个标准正态随机数,3.5,蒙特卡洛模拟,3.5.3,蒙特卡洛模拟的随机变量抽样方法,2,、正态分布随机数产生,如果,M=16,a=5,c=3,x,0,=7,正态分布均值为,100,标准差为,20,x1=7;a=5;c=3;m=16,for i=1:40,x1=mod(x1*a+c),m),if x11,x1=7,end,xx(i)=x1,jj(i)=i,end,u=xx/16,for i=1:2:30,x(i)=(sqrt(-2.0*log(u(i)*cos(2*pi*u(i+1)*20+100,x(i+1)=(sqrt(-2.0*log(u(i)*sin(2*pi*u(i+1)*20+100,j(i)=i,j(i+1)=i+1,end,plot(j,x,-),hold on,plot(j,x,*),产生,30,个一般正态分布随机数,在正态分布公式中令,z=(t-)/,，可将随机变,量,X,标准化，标准化后的,随机变量,z,服从标准正态分布。,则：,t=+z,3.5,蒙特卡洛模拟,3.5.3,蒙特卡洛模拟的随机变量抽样方法,3,、对数正态分布随机数产生,对数正态分布随机数可以由正态分布随机数转化而来。假设随机变量,X,服从对数正态分布，则随机变量,Y=lnX,服从正态分布。,Y,的标准差和变异系数可以由,X,的标准差和变异系数求得：,然后根据正态分布随机数产生的方法求得,Y,的随机数,y,i,，则,x,i,服从对数正态分布的随机数为,x,i,=exp(y,i,),。,3.5,蒙特卡洛模拟,3.5.3,蒙特卡洛模拟的随机变量抽样方法,3,、对数正态分布随机数产生,对数正态分布随机数产生随机数方法如下：,(1),首先求相应正态分布的均值和标准差。,(3)x,i,服从对数正态分布的随机数为,x,i,=exp(y,i,),。,(2),产生正态分布随机变量,一般正,态分布,3.5,蒙特卡洛模拟,3.5.3,蒙特卡洛模拟的随机变量抽样方法,2,、对数正态分布随机数与正态分布随机数产生的区别,如果,M=16,a=5,c=3,x,0,=7,，样本均值为,100,，标准差为,20,x1=7;a=5;c=3;m=16,for i=1:40,x1=mod(x1*a+c),m),if x11,x1=7,end,xx(i)=x1,jj(i)=i,end,u=xx/16,for i=1:2:30,x(i)=(sqrt(-2.0*log(u(i)*cos(2*pi*u(i+1)*20+100,x(i+1)=(sqrt(-2.0*log(u(i)*sin(2*pi*u(i+1)*20+100,j(i)=i,j(i+1)=i+1,end,plot(j,x,-),hold on,plot(j,x,*),产生,30,个一般正态分布随机数,x1=7;a=5;c=3;m=16,for i=1:40,x1=mod(x1*a+c),m),if x11,x1=7,end,xx(i)=x1,jj(i)=i,end,u=xx/16,d=0.1;mu=log(100/sqrt(1+0.22);m=log(1+0.22),for i=1:2:30,x(i)=(sqrt(-2.0*log(u(i)*cos(2*pi*u(i+1)*m+mu,x(i)=exp(x(i),x(i+1)=(sqrt(-2.0*log(ui)*sin(2*pi*u(i+1)*m+mu,x(i+1)=exp(x(i+1),j(i)=i,j(i+1)=i+1,end,plot(j,x,-),hold on,plot(j,x,*),产生,30,个对数正态分布随机数,例题,设结构裂缝宽度服从对数正态分布，已知样本的均值为,2.67,，方差为,2.90,，试产生两个对数正态分布随机数。,解：首先求正态分布的均值和标准差。,=2.9,1/2,/2.67=0.6389,根据下列公式产生两个随机数：,由,M=16,a=5,c=3,x,0,=7,，产生两个,(0,1),分布随机数：,6/16,和,1/16,3.5,蒙特卡洛模拟,3.5.3,蒙特卡洛模拟的随机变量抽样方法,4,、极值,I,型分布随机数产生,极值,I,型分布函数为：,若产生均匀分布随机数为,u,i,，极值,I,型分布随机数产生采用下式：,例题：蒙特卡洛模拟可靠度分析方法,设某构件正截面强度计算的极限状态方程为,Z=R-S=0,。其中,R,和,S,分别为正态和极值,I,型分布的随机变量，其统计量为,R(100,20),和,S(80,24),，,20,和,24,为标准差。试用蒙特卡罗模拟求解构件失效概率。,解：,(1),产生,(0,，,1),伪随机数,(2),产生变量,R,和,S,的随机数：,N=800,次,失效概率为,24.25%,，不同方法产生随机数，计算结果稍有差别x1=7;a=5;c=3;m=32,for i=1:1100,x1=mod(x1*a+c),m),if x11,x1=7,end,xx(i)=x1,jj(i)=i,end,u=xx/16,for i=1:2:1000,x(i)=(sqrt(-2.0*log(u(i)*cos(2*pi*u(i+1)*20+100,x(i+1)=(sqrt(-2.0*log(u(i)*sin(2*pi*u(i+1)*20+100,end,for i=1:1000,xx(i)=80-0.45*24-0.7797*24*log(-log(u(i),end,nm=0,N=800,for i=1:N,z=x(i)-xx(i),if z=0,nm=nm+1,end,end,pf=nm/N,3.5,蒙特卡洛模拟,3.5.3,蒙特卡洛模拟的随机变量抽样方法,蒙特卡洛法模拟步骤：,设功能函数,Z=g(X,1,X,2,X,n,),，式中：,X,i,(i=1,2,n),为具有任意分布的随机变更。,(1),产生在区间,0,1,上均匀分布的,N,组随机数；,(2),将产生的均匀分布随机数变换成给定分布的随机数，得,N,组,X,i,j,(i=1,2,n,；,j=1,2,N),随机变量样本。,(3),将第,j,组的,X,i,j,(i=1,2,n),的值代入功能函数,Z,j,=g(X,1,X,2,X,n,),，如果,N,个,Z,j,值中有,N,f,个,Z,j,0,，则结构失效的概率为：,