《一次二阶矩法ppt课件.ppt》由会员分享,可在线阅读,更多相关《一次二阶矩法ppt课件.ppt(39页珍藏版)》请在三一办公上搜索。
1、ch4一次二阶矩法,引言一次二阶矩法是求解非时变荷载作用下结构可靠度问题的行之有效的近似方法。它既有较高的精度,又有较高的计算效率 。一次二阶矩法是在基本变量xi(i= 1,2n),的概率分布尚不清楚中时,采用只有均值(又称为一阶原点矩)和标准差(又称为二阶中心矩)的数学模型去求解结构可靠度的方法 。,Ch4 一次二阶矩法,4.1 一次二阶矩中心点法4.2 一次二阶矩验算点法4.3 JC法,4.1一次二阶矩中心点法,4.1.1.简单问题可靠指标的求法4.1.2.中心点法基本原理 4.1.3.中心点法举例 4.1.4. 对中心点法的评述,在介绍中心点法之前,先回顾功能函数为线性函数,随机变量为正
2、态分布情况下的可靠指标问题。,I随机变量为正态分布,且功能函数为线性函数,假定抗力R和荷载效应S均服从正态分布,对于功能函数Z=R-S,,由于Z是R、S的线性函数,根据正态随机变量的特性,Z也服从正态分布,其平均值及标准差分别为:,则定义可靠指标:,为什么可以定义来衡量结构的可靠程度?,4.1.1简单问题可靠指标(或失效概率)的求法,正态分布功能函数Z,其失效概率与可靠指标之间的精确关系,这种情况下可靠指标是唯一的,且与失效概率之间有精确的对应关系。,II随机变量为对数正态分布,功能函数为线性函数,对数正态分布的定义:若lnR、lnS 服从正态分布,则称R、S服从对数正态分布。假定抗力R和荷载
3、效应S均服从对数正态分布,且结构功能函数可表示为Z=lnR-lnS,由于lnR、lnS均服从正态分布,Z也服从正态分布,是lnR、lnS的线性函数,根据正态随机变量的特性,其平均值及标准差也可以精确得到:,可靠指标为:,根据x的统计参数,求 lnx的统计参数,结论:随机变量为正态分布(或对数正态分布),且功能函数为线性函数的条件下,结构的可靠指标或失效概率可根据随机变量的统计参数方便的求出。,随机变量为正态分布,由其形成的线性函数也服从正态分布。,4.1.2中心点法引言,假定随机变量x1、x2xn服从任意分布,功能函数 不是线性函数,这时, 精确求解Z的平均值和标准差是非常困难的,即便能够求得
4、,Z也不服从正态分布,也不能用上面方法来计算结构的可靠指标。 若将非线性功能函数作为泰勒级数展开,并取其一次展开式(前两项)。但一次展开式已不是原来的功能函数,所计算可靠指标与结构失效概率之间不再存在精确的对应关系。在这种情况下如何选择展开点,从而使近似计算结果与精确失效概率的误差最小,成为一次二阶矩法要研究的问题。,4.1.2中心点法 基本原理,设结构的极限状态方程为 将极限状态函数在中心点M= ( ) 处展开为泰勒级数,并作线性化处理,得根据概率论中随机变量参数估计 ,Z*的统计参数为: 结构的可靠指标,在中心点处展开为泰勒级数:,4.1.3中心点法举例,钢梁承受确定性弯矩 ,截面塑性抵抗
5、矩W和屈服强度f都是随机变量,已知统计参数为:W:f:试用中心点法计算该钢梁抗弯的可靠指标,4.1.3中心点法举例,解法1取功能函数 极限状态方程为,4.1.3中心点法举例,解法2取功能函数 极限状态方程为,4.1.4对中心点法的评述,与随机变量为正态分布、功能函数为线性函数的情况下求可靠指标,形式一样( ),却有本质的区别:(1)中心点法不论R、S的概率分布如何,直接用其平均值、标准差求可靠指标。事实上,平均值、标准差相同的随机变量有无穷多个,每一种情况的可靠度都是不同的,而这里是同一个结果。(2)前者的可靠指标是唯一的,并与失效概率有精确的对应关系,而中心点法可靠指标与建立的功能函数表达式
6、的形式有关。,4.1.3对中心点法的评述,中心点法的主要弱点没有考虑基本变量的概率分布 均值、方差及可靠指标的计算式是误差传递公式同一个结构往往可以列出几种等价的极限状态方程,不同的极限状态函数在运用中心点法计算时,其结果可能不一致。将非线性功能函数在随机变量的平均值处展开不合理,由于平均值不在极限状态曲面上,展开后的线性极限状态面可能会较大程度地偏离原来的极限状态曲面。 基本变量不服从正态分布和对数正态分布时,计算出的结构可靠指标与结构的实际情况出入较大。,中心点法的线性示意图,该法选用的线性化点(即平均值点)不在失效边界上,4.2验算点法(改进的一次二阶矩法),4.2.1验算点法基本思想运
7、用泰勒级数进行线性化处理时,略去高阶项的误差随线性化点到失效边界的距离增加而增大(中心点法运用泰勒级数在中心点处展开进行线性化处理)。选择在失效边界Z0的某点处将极限状态函数展开为泰勒级数,以避免这种形式的误差。验算点法,4.2.1验算点法基本思想,关于设计验算点设点P(x1,x2,xn)为极限状态方程Z=0所对应的曲面上的点,d(P,M)为点P到中心点M( ) 的距离,则能使mind(P,M)的点P*称为设计验算点,简称为验算点。记为P*(x*1,x*2,x*n),显然验算点的坐标满足,验算点法示意,M,验算点法,4.2.2设计验算点法求可靠指标,理论推导确定验算点位置进而求可靠指标当线性化
8、点选在设计验算点xi*(i=1,2,n)上时 Z的均值为(求Z的数学期望得)由于设计验算点在失效边界上,故有 则有,假设各随机变量独立,则可求解Z的方差: 引入分离函数式,将上面的根式线性化,得 表示第个随机变量对整个标准差的相对影响,因此称为灵敏系数。,4.2.2设计验算点法求可靠指标,4.2.2设计验算点法求可靠指标,根据可靠指标的定义,有由于 , 必有 (对于所有i),将 乘到分母上整理后得,4.2.2设计验算点法求可靠指标,式(13)代表n个方程,再加上(14)共有n+1个方程,未知数有 和 ,也是n+1个。可联立求解。解出设计验算点P*( , , )及相应的常用迭代法求解,4.2.2
9、设计验算点法求可靠指标,计算步骤(1)选取设计验算点坐标的初值,一般取 (2)由式(11) 计算 的值,其中包括 (3)由式(13)得到 和 的关系 (4)由式(14)解出 值 (5)将该值代入式(13),求出 新值 以该 新重复进行(2)-(5)计算,直到 值与上次相等或误差不超过允许值,此时 即为所求的可靠指标, 即为所确定的设计验算点坐标。,4.2.3验算点法举例,设极限状态功能函数为,均为随机变量,并具有如下统计信息:,试用验算点法计算结构的可靠指标,1.取均值作为验算点初值,2.计算,0.2639,0.7037,-0.6597,3. 与 的关系,从而有新的,4.由(14)求 及新的,
10、4.2.3验算点法举例,5.以新 进行第二次迭代,重复4直到 值与上次相等或误差不超过允许值,此时即为所求的可靠指标 , 即为所确定的设计验算点坐标,本例第三次迭代求得,第四次迭代求得,即为最终结果,设计验算点法求可靠指标评述,验算点法对极限状态方程中服从正态分布的随机变量计算结果尚可,而非正态分布误差较大。,JC法求可靠指标,问题的提出:验算点法对极限状态方程中服从正态分布的随机变量计算结果尚可,而非正态分布误差较大。在工程结构可靠度分析中,永久荷载一般服从正态分布,风压、雪压、楼面活载服从其它类型(如极值型等),截面抗力R服从对数正态分布。因此在极限状态方程中,常包含非正态变量分布的基本变
11、量,对于这种极限状态方程的可靠度分析,一般要把非正态变量进行当量正态化。JC法的基本概念就是在应用前面所述方法(验算点法)时,将非正态的随机变量先行“当量正态化”。JC法是由Rackwitz-Fiessler、Hasofer-Lind等人先后提出来的,因为国际安全度联合委员会(JCSS)推荐采用这个方法而得名。,JC法求可靠指标计算步骤,设结构的极限状态方程为Z=g(x1,x2,xn)=0,基本变量xi的统计参数均值为 标准差为 ,(i=1,n),则JC法的具体实施步骤加下: (l)假定基本随机变量的设计验算点P*的坐标值 (一般取为 ) (2)对基本变量xi进行当量正态化处理,计算其当量正态
12、分布yi的统计参数 , ,并用来代替 和 并记为 , ( 1,n);(3)计算 的值,其中包括,JC法求可靠指标计算步骤,(4) 确定 和 的关系 (5)求(6)求出 新值 (7)以该新 重复进行(2)(5)计算 ,直到 值与上次相等或误差不超过允许值,此时 即为所求的可靠指标, 即为所确定的设计验算点坐标。,JC法计算实例,某轴心受压短柱,承受永久荷载产生的压力SG,汽车、人群可变荷载产生的压力SQ,截面抗力为R。各变量的统计信息如下表,试用JC法求可靠指标及对应的失效概率。,JC法计算实例,1.建立极限状态方程 Z=g(R , SG , SQ)=R-SG-SQ=02.假定初始的设计验算点,由已知条件得: x1*=R*=365.8;x2*=SG*=63;x3*=SQ*=84.03.将非正态变量SQ和R在均值处进行当量正态化处理抗力R对数正态分布,则,JC法计算实例,可变荷载SQ极值I型分布,则,JC法计算实例,4计算,Z=g(R , SG , SQ)=R-SG-SQ=0,JC法计算实例,5. 和 的关系,JC法计算实例,6.将求得的 代回5重新R*、,7. 以第6步得值进行第二次迭代,此时,R*108.16,R的当量正态化处理:,JC法计算实例,SQ的当量正态化处理:,