结构可靠度计算方法(一次二阶矩)课件

《结构可靠度计算方法(一次二阶矩)课件》由会员分享，可在线阅读，更多相关《结构可靠度计算方法(一次二阶矩)课件（56页珍藏版）》请在装配图网上搜索。

1、 s o u t h w e s t j I a o t o n g u n I v e r s I t y西南交通大学西南交通大学Southwest Jiaotong University隧道与地下结构可靠度隧道与地下结构可靠度课程课程第第 三三 讲讲结构可靠度计算方法结构可靠度计算方法主要内容主要内容 基本概念基本概念 一一次二阶矩理论的次二阶矩理论的中心点法中心点法 一一次二阶矩理论的次二阶矩理论的验算点法验算点法(JC法法)映射变换法映射变换法 实用分析法实用分析法 s o u t h w e s t j I a o t o n g w n I v e r s I t y一、基本概念一

2、、基本概念西南交通大学西南交通大学Southwest Jiaotong University现现代的结构可靠度理论是以代的结构可靠度理论是以概率论概率论和和数理统计学数理统计学为基础发展起来的，要解为基础发展起来的，要解决的决的中心问题中心问题是围绕着怎样描述和分析是围绕着怎样描述和分析可靠度可靠度，以及研究影响可靠度各基本变，以及研究影响可靠度各基本变量的量的概率模型概率模型。1 1、解决的问题、解决的问题结构可靠度结构可靠度计算方法计算方法分分精确法精确法和和近似法近似法两种。两种。精确法：求解结构的失效概率精确法：求解结构的失效概率 pf 的方法，的方法，通常称为全概率法；通常称为全概率

3、法；近似法：一次二阶矩计算方法等，虽然近似法：一次二阶矩计算方法等，虽然是近似的，但仍属概率法。是近似的，但仍属概率法。2 2、计算方法、计算方法结构功能函数大多是结构功能函数大多是非线性非线性函数，且非线函数，且非线性不是很强的条件下，但又性不是很强的条件下，但又不能不能直接精确积分直接精确积分计算得到结构的可靠度，而计算得到结构的可靠度，而通过通过计算结构计算结构可靠可靠指标指标，近似得到结构可靠度的计算方法。，近似得到结构可靠度的计算方法。在通常情况下，结构功能函数的在通常情况下，结构功能函数的一阶矩一阶矩(均值均值)和和二阶矩二阶矩(方差方差)较容易得到较容易得到，故称之为故称之为一次

4、二阶矩法一次二阶矩法。3 3、一次二阶矩法、一次二阶矩法一次二阶矩法一次二阶矩法是一种在随机变量的分布尚是一种在随机变量的分布尚不清楚不清楚时，采用时，采用均值均值和和标准差标准差的数学模型，求的数学模型，求解结构的可靠指标、结构可靠度的方法。解结构的可靠指标、结构可靠度的方法。该法将功能函数该法将功能函数 在在某点用泰勒级数展开，使之线性化某点用泰勒级数展开，使之线性化，然后求，然后求解结构的可靠度，因此称为一次二阶矩。解结构的可靠度，因此称为一次二阶矩。),(21nxxxgZ s o u t h w e s t j I a o t o n g w n I v e r s I t y二、一次

5、二阶矩理论的中心点法二、一次二阶矩理论的中心点法西南交通大学西南交通大学Southwest Jiaotong University中心点法是结构可靠度研究初期提出的中心点法是结构可靠度研究初期提出的一种方法一种方法。其基本思想其基本思想：首先，将非线性功能函数首先，将非线性功能函数在随机变量的在随机变量的平均值平均值(中心点中心点)处作处作泰勒级数展泰勒级数展开开，并保留至并保留至一次项一次项；然后，近似计算功能函；然后，近似计算功能函数的数的平均值平均值和和标准差标准差。1 1、一次二阶矩中心点法、一次二阶矩中心点法设设X1,X2,Xn是结构中是结构中n个相互独立的随机个相互独立的随机变量，

6、其变量，其平均值平均值为为 ,标准差标准差为为 ,功能函数功能函数将功能函数将功能函数Z在平均值在平均值P*(X1,X2,Xn)处处展开且保留至一次项，即展开且保留至一次项，即 (3-1),(21nxxxgZiXiXinXiniiXXXXXggZ1),(212 2、推导过程、推导过程ZL平均值和方差为平均值和方差为：(3-2)2212211),()(),()(2121iLninLXniiLLZXXXniXiniiXXXLZXgZEZEgXEXggZE结构可靠指标为结构可靠指标为 (3-3)212,21inLLXniiXXXZZXgg可靠指标可靠指标的几何意义是什么？证明如下的几何意义是什么？证

7、明如下功能函数泰勒级数展开至一次项，即功能函数泰勒级数展开至一次项，即 (3-4)假定正态变换，即：假定正态变换，即：(3-5)inXiniiXXXXXggZ1),(21iiXXiiXX3 3、几何意义、几何意义将将(3-5)式代入式代入(3-4)式，得式，得 (3-6)inXniiiXXXXXggZ1),(21(3-6)式式为一个超平面方程，点为一个超平面方程，点P*(X1,X2,Xn)到平面的距离为：到平面的距离为：(3-7)niXiXXXinXggd122),(21中心点法中心点法验算点法验算点法极限方程曲面极限方程曲面可靠区可靠区均值点均值点显然，点显然，点P*(X1,X2,Xn)到平

8、面的距离到平面的距离d，就是所求的可靠指标值就是所求的可靠指标值，两者是相等的。，两者是相等的。P*优点：优点：计算简便。计算简便。缺点：缺点：对于非线性功能函数，对于非线性功能函数，均值点均值点一般在可靠区一般在可靠区内内，而，而不在极限边界上；不在极限边界上；选择选择不同极限状态方程不同极限状态方程（数学表达式不同，（数学表达式不同，同样物理含义），得到的同样物理含义），得到的可靠指标不同可靠指标不同。例例如：如：p30例例3-1。适用条件：适用条件：结果比较粗糙，适用于可靠度要求结果比较粗糙，适用于可靠度要求不高的情况，如钢筋混凝土结构正常使用极限不高的情况，如钢筋混凝土结构正常使用极限

9、状态的可靠度分析。状态的可靠度分析。4 4、优缺点、优缺点例题例题1设设X1,X2,Xn是结构中是结构中n个相互个相互独立的随机变量，其平均值为独立的随机变量，其平均值为xi（i=1,2,n)，标准差为标准差为xi(i=1,2,n)，功能函数，功能函数Z=g(X1,X2,Xn)。求结构可靠指标。求结构可靠指标？解解 将功能函数将功能函数Z在随机变量的平均值处在随机变量的平均值处泰勒级数展开，且保留一次项，即泰勒级数展开，且保留一次项，即iXiniinXXggZ121),(5 5、举例、举例 ZL的平均值和方差为的平均值和方差为：结构可靠指标为结构可靠指标为：),()(21nXXXZgZE221

10、22122)()(iLXniiiiiniiZXgpXEXXgZEZEniXiXXXZZinXgg122),(21例题例题2 某结构构件正截面强度的功能函数某结构构件正截面强度的功能函数为为Zg(R,S)=R-S，其中抗力，其中抗力R服从对数正态分服从对数正态分布，布，R=100kNm，R=0.12；荷载效应；荷载效应S服从极服从极值值I型分布，型分布，S=50kNm，S=0.15。试用中心点。试用中心点法求结构失效概率法求结构失效概率Pf?解解：1212.0100RRR5.715.050SSS 结构可靠指标结构可靠指标 结构失效概率结构失效概率 533.35.712501002222SRSR4

11、10078.2)533.3()(fp s o u t h w e s t j I a o t o n g w n I v e r s I t y三、一次二阶矩理论的验算点法三、一次二阶矩理论的验算点法西南交通大学西南交通大学Southwest Jiaotong UniversityJC法是法是Hasofer,Lind,Rackwitz和和Fiessler,Paloheimo和和Hannus等人提出的验算点法。等人提出的验算点法。适用于随机变量为适用于随机变量为非正态分布非正态分布的结构可靠指标的结构可靠指标的计算。的计算。通俗易懂，计算精度又能满足工程实际需要。通俗易懂，计算精度又能满足工程实

12、际需要。国际结构安全度联合委员会国际结构安全度联合委员会(JCSS)推荐使用，推荐使用，故称为故称为JC法法。我国我国建筑结构设计统一标准建筑结构设计统一标准(GBJ68-84)和和铁路工程结构设计统一标准铁路工程结构设计统一标准(GB50216-94)中都中都规定采用规定采用JC法进行结构可靠度计算法进行结构可靠度计算。1 1、验算点法（、验算点法（JCJC法）法）将将P*（X*1,X*2,X*n)定义为验算点（设计点），故定义为验算点（设计点），故称之为称之为验算点法验算点法。又因为是在中心点法的基础上改进的，。又因为是在中心点法的基础上改进的，故故称为称为一次二阶矩的改进方法一次二阶矩的

13、改进方法。数学推导过程如下：数学推导过程如下：设设X1,X2,Xn(i=1,2,n)为基本变量，且相互独立，则为基本变量，且相互独立，则极限极限状态功能函数状态功能函数方方程为：程为：(3-8)将极限方程用泰勒级数在将极限方程用泰勒级数在P*（X*1,X*2,X*n)点上展点上展开，取一次项，可得极限方程为：开，取一次项，可得极限方程为：(3-9),(21nXXXgZ0),(*1*2*1*iipniinXXXgXXXgZ2 2、推导过程、推导过程设设 (3-10)有有 (3-11)将将(3-11)代入代入(3-9)，得，得 (3-12)iiXXiiXX*0),(*1*2*1*iXiipniin

14、XXXgXXXgiXiiiiiXgXXXgXg*Z的平均值为：的平均值为：(3-13)验算点在极限边界上，即验算点在极限边界上，即又又 (3-14)将将(3-14)代入代入(3-13)，得，得 (3-15)*1*2*1)(),()(*iiXpniinZXXEXgXXXgZEi0),(*2*1nXXXgiXiXE)(iiXiXpniiZXXg*1*2.1 2.1 按定义推导按定义推导Z的标准差的标准差Z为：为：(3-16)则可靠指标则可靠指标为：为：(3-17)2112*niXPiZiXg21121*niXPiniXXipiiiiXgXXg随机变量满足正态分布，即随机变量满足正态分布，即 (3-

15、18)其中：其中：(3-19)*coscos*iiiiiiXXXiXiXXiXXXniXPiXPiXiiiXgXg12*cos由由(3-12)，得得 (3-19)此为超平面方程，均值点此为超平面方程，均值点P(X1,X2,Xn)到超平面的距离到超平面的距离d为：为：(3-20)0*1*iXiipniiXXXg21121*niXPiniXXipiiiiXgXXgd2.2 2.2 按几何意义推导按几何意义推导各变量的方向余弦为：各变量的方向余弦为：(3-21)显然，两种方法得到的结果是一致的。显然，两种方法得到的结果是一致的。niXPiXPiXiiiXgXg12*cos将将(3-8)与与(3-18

16、)联立，求得联立，求得和各变量值，和各变量值，再代入到再代入到(3-8)和和(3-18)，且联立求解，得到新，且联立求解，得到新的一组的一组和各变量值。和各变量值。直到满足下式为止，即直到满足下式为止，即 (3-22)迭代结束，计算完成。迭代结束，计算完成。1nn2.3 2.3 迭代过程迭代过程两个随机变量为正态分布时，其极限方程为两个随机变量为正态分布时，其极限方程为 0),(SRSRgZRRSSRRSS0SRSRSR标准化变换标准化变换极限状态方程变为极限状态方程变为(3-23)(3-24)(3-25)3 3、正态分布时的推导过程、正态分布时的推导过程0coscosRSRS222222co

17、scosSRSRSRRRSRsS式中式中将将(3-25)变为标准法线式直线方程变为标准法线式直线方程(3-26)(3-27)是坐标系是坐标系 中原点中原点 到极限状态直到极限状态直线的距离线的距离 (其中其中P*为垂足为垂足)。在验算点法中，在验算点法中，的计算就转化为求的计算就转化为求 的长度。的长度。RSOO*PO*PO两个正态随机变量的极限状态方程和设计验算点*)*,(*RSPOSRSSSRRRO 非正态分布时，可采取以下三种方法：非正态分布时，可采取以下三种方法：当量正态化法（当量正态化法（JC法）法）映射变换法映射变换法实用分析法实用分析法 JC法为法为当量正态化法当量正态化法，将原

18、来非正态分布随，将原来非正态分布随机变量机变量Xi用等效正态分布代替，用等效正态分布代替，要求满足要求满足以下以下2个条件个条件：原函数值原函数值F(xi*)与当量正态函数值与当量正态函数值F(xi*)相等相等原概率密度值原概率密度值f(xi*)与当量正态分布概率密度值与当量正态分布概率密度值f(xi*)相等相等iX4 4、非正态分布时、非正态分布时JCJC法的等效正态分布图法的等效正态分布图*ix)(iXxfi原分布原分布FXi(xi*),fXi(xi*)等效正态分布等效正态分布FXi(xi*),fXi(xi*)ixO条件条件(1)和和(2)的数学表达式为的数学表达式为 (3-20)(3-2

19、1)由由(3-20)，得，得 (3-22)由由(3-21)，得，得 (3-23)iiXXiixxF*)()()()()(*iXiXiXiXxfxfxFxFiiiiiiiiiiiXXXiiXXiiiXiXxdXxddXxdFxf1)()(*4.1 4.1 当量正态化法当量正态化法-JC-JC法法由由(3-22)，得，得 (3-24)将将(3-24)代入代入(3-23)，得，得 (3-25)由由(3-24)和和(3-25)，得，得 (3-26)*1*iXXXixFxiiiiiiXiXiXxFxf1)()(*1*)(*1*1*iXiXXXiXiXxfxFxFxiiiiii将将(3-26)代入代入(3

20、-18)、(3-19)和和(3-20)进行迭进行迭代计算，就可求解随机变量非正态分布的可靠代计算，就可求解随机变量非正态分布的可靠度问题。度问题。显然，显然，JC法通过法通过当量变换当量变换，使得非正态分，使得非正态分布的随机变量满足正态分布要求，进而应用满布的随机变量满足正态分布要求，进而应用满足正态分布的方法进行迭代计算，求解非正态足正态分布的方法进行迭代计算，求解非正态分布随机变量的可靠度问题。分布随机变量的可靠度问题。李云贵李云贵(1993)提出映射变换法。具体数学过程如下：提出映射变换法。具体数学过程如下：设结构中的设结构中的n个相互独立的随机变量为个相互独立的随机变量为X1,X2,

21、Xn，其概率分布函数为其概率分布函数为Fi(xi)(i=1,2,n)，概率密度函数为，概率密度函数为fi(xi)(i=1,2,n)，极限状态方程为，极限状态方程为 Zx=g(X1,X2,Xn)=0 (3-27)映射变换映射变换 (3-28)则则 (3-29),2,1)()(niYXFiii iiiiiiXFYYFX114.2 4.2 映射变换法映射变换法将将(3-29)代入代入(3-27)，得，得 (3-30)由于由于Yi是一个标准正态随机变量是一个标准正态随机变量，则，则 (3-31)于是于是 (3-32)(3-33)(3-34),(,211212111nnnYYYYGYFYFYFgZ1,0

22、iiYYiYiycos*niPiPiYYGYGi12*cos0),(*2*1nyyyG其中：其中：(3-35)对于常用的几种概率对于常用的几种概率(1)Xi服从正态分布服从正态分布 (3-36)(2)Xi服从对数正态分布服从对数正态分布 (3-37),2,1(niYXXGYGiiiiiiiXiiXiXiYXYXiiiXiiiXiXiXYXYXlnlnln)exp(3)Xi服从极值服从极值I型分布型分布 (3-38)式中，式中，)(ln)()(lnlniiiiiiiYYYYXYXiiiXXX16,45.0Paloheimo和和Hannus(1972)在赫尔辛基工程力学学术讨论在赫尔辛基工程力学学

23、术讨论会提出了会提出了分位值法分位值法。有的著作中是中国铁科院姚明初有的著作中是中国铁科院姚明初(1993)提出的分位值法。提出的分位值法。本文认为仍然是本文认为仍然是Paloheimo和和Hannus(1972)提出的。提出的。所谓分位值法就是所谓分位值法就是映射变换法映射变换法的另一种表述。的另一种表述。n个独立随机变量个独立随机变量X1,X2,Xn，结构极限状态方程为，结构极限状态方程为 (3-39)映射变换映射变换 (3-40)则有则有 (3-41)0),(21nXXXGZ),2,1)()(niYXFiii iiiiiiXFYYFX114.3 4.3 分位值法分位值法将将(3-39)在

24、在 点进行泰勒级数展开且保留一次点进行泰勒级数展开且保留一次项，即：项，即：(3-42)于是有于是有 (3-43)其中其中niiPiiPiniPiPiYXXGXXGYGYGi12212*cosiiiiiiXXGYXXGYG),(*2*1*nYYYP0),(*121iiniinYYYGXXXG (3-44)是基本变量是基本变量 对应分位概率为对应分位概率为 的分位值。的分位值。(3-45)各变量各变量 的的“分项可靠指标分项可靠指标”可用可用(3-45)求解得到。求解得到。于是得到于是得到“设计值设计值”，即，即 (3-46)*iXniYiniiYiniiYGYYGXXXG121*2*1*),(

25、iX)(*iYiYiYcos*iYiX*1*iiiYFX*iX变量变量Xi的的“设计值设计值”与变量与变量Xi的标准值的标准值Xik之比定义为分项之比定义为分项系数，即系数，即 (3-47)*iXikiiXX*JCJC法法-举例举例例题例题 已知极限状态方程：已知极限状态方程：(1)(2)随机变量随机变量f，W均服从正态分布，均服从正态分布，f=38，f=0.1；W=54，W=0.05，求：两个极限状态方程条件下的，求：两个极限状态方程条件下的及及f和和W的的验算点之值验算点之值f*，W*。解解 1.首先求首先求(1)式条件下的式条件下的及及f和和W的验算点之值的验算点之值f*，W*。0114

26、0),(fWWfgZ06.80),(WfWfgZ;8.31.038fff7.205.054WWW由由 (3)得得niXpiXpiXXgXgi12*11cos5792.08.37.27.2cos8152.08.37.28.3cos222222WWfff由由 (4)得得 (4-1)将上式代入将上式代入(1)式，得式，得 (4-2)将将(4-2)代入代入(4-1)，得，得由于由于 为常数，所以不需为常数，所以不需要迭代就求解出要迭代就求解出，f*，W*。5638.154cos*0978.3388152.08.338cos*WWWfffWfiiiXXXiXcos*5742.23106.495638.1

27、54*7106.289987.20978.3380978.338*Wf5792.0cos,8152.0cosWf对于极限状态方程对于极限状态方程(1)，若采用中心点求解可靠指标，得，若采用中心点求解可靠指标，得当极限状态方程为线性方程时，采用验算点法与中心点法当极限状态方程为线性方程时，采用验算点法与中心点法求解得到的结果是一致的。求解得到的结果是一致的。5742.27.28.380543880),(22221221WfWfniXinZZiXgXXXg2.按极限状态方程按极限状态方程(2)求解求解及及f和和W的验算点之值的验算点之值f*，W*。由由(3)式，得式，得由由(4-1)，得，得222

28、2222244.1429.77.2cos44.1429.78.3cosWffWfWWfWWWfffWWWWffffWfcos7.254cos*cos8.338cos*将上式代入将上式代入(2)式，得式，得 (5)显然，显然，(5)式的求解需要采用迭代法求解。式的求解需要采用迭代法求解。2.1 第一次迭代第一次迭代取取将上式代入将上式代入(5)，得，得09.88)cos10cos20(coscos2WfWf54,38*WfWf4472.0cos,8944.0cosWf02.2229.552307.4第二次迭代第二次迭代得得显然，需要再迭代求解。显然，需要再迭代求解。8.48*36.23)8944

29、.0(307.48.338cos8.338*Wff322.0cos9467.0cosWf001.00452.0307.42618.4001.02618.41212第三次迭代第三次迭代得得满足迭代精度要求，迭代求解结束。满足迭代精度要求，迭代求解结束。29.50*67.22*Wf3049.0cos9524.0cosWf001.00009.02618.42627.4001.02627.42312将将 代入代入(4-1)，得，得若极限状态方程若极限状态方程(2)应用中心点法求解可靠指标，有：应用中心点法求解可靠指标，有：29.50*67.22*Wf2627.4912114054381140WfZ42.229)7.238()8.354(2222WfZWgfg98.342.229912ZZ将中心点法与验算点法计算结果进行比较将中心点法与验算点法计算结果进行比较经比较发现，中心点法在极限状态方程是非线性方程，求经比较发现，中心点法在极限状态方程是非线性方程，求解得到的结果与验算点法有较大误差。但若极限状态方程为线解得到的结果与验算点法有较大误差。但若极限状态方程为线性方程，则两者求解结果一致。性方程，则两者求解结果一致。2627.4JC934.02627.498.3JCCen98.3Cen