1、第四讲 损伤理论,荷载谱处理,应变疲劳上节概述p-S-N曲线,失效率,存活率疲劳数据处理,正态分布,标准正态分布检验数据是否服从正态分布,正态分布坐标纸正态分布存在的问题:不能反映构件疲劳寿命有一个大于等于零的下限。威布尔分布,与正态分布的比较威布尔分布坐标纸回归方程,最小二乘法相关系数,起码值 回归分析的基本方法损伤理论疲劳积累损伤理论是构件在变幅疲劳荷载作用下疲劳损伤的积累规则和疲劳破坏的准则。疲劳积累损伤理论回答下述三个问题1一个荷载循环对材料或结构造成的损伤是多少?2多个荷载循环时,损伤是如何积累的?3失效时的临界损伤是多少?1)线性疲劳积累损伤理论线性疲劳积累损伤理论假设在循环荷载作
2、用下,疲劳损伤是可以线性累加,各应力之间相互独立,当累加损伤到某一数值构件发生疲劳破坏。PalmgrenMiner(PM)线性疲劳积累损伤准则若构件在某恒幅应力水平S作用下的疲劳寿命为N,则经受n次循环时的损伤为:N1N2S1S3SmaxNS2N3 构件在应力水平Si作用下经受ni次循环的损伤为Di = ni/Ni,在k个应力水平作用下的总损伤为: 破坏准则:Miner疲劳积累损伤理论没有考虑荷载作用的先后次序。材料“锻练效应”(training effect)一般,高低加载顺序临界积累损伤值DCR小于1,低高加载顺序临界积累损伤值DCR大于1。变幅拉压疲劳寿命疲劳加载形式序号n1n2两级高低
3、1641165200.7290.511264122600.2723641103190.531低高1624232221.2051.9942624228211.08036242112083.697对于随机荷载,按Miner理论计算破坏时的临界损伤值DCR接近于1。Manson双线性模型Manson在修正Miner准则时提出的两级疲劳加载时的疲劳寿命预测公式 :与疲劳加载顺序相关的参数当疲劳加载从高到低时0 1。即该模型实质上是将第一级应力水平下的疲劳损伤等效为第二级应力水平下引起的损伤。且在等幅加载或三级以上疲劳加载下并不适用。2)非线性疲劳积累损伤理论CartenDolan非线性疲劳积累损伤准则
4、CartenDolan从疲劳破坏过程的损伤微观物理模型出发,给出材料经受n次循环时的损伤为 m:材料损伤核数目r:损伤发展速率,正比于应力水平c、d:材料常数则构件在k个应力水平作用下的总损伤为 临界疲劳损伤:N1:作用的荷载系列中最大一级荷载所对应的疲劳寿命破坏准则:疲劳损伤核在后续加载过程中不会消失,因此mi = m1 损伤发展速率r正比于应力水平S,ri Si CartenDolan模型实质上是将任意应力水平Si下循环ni次所引起的材料损伤等效为最大一级荷载S1作用nie次所引起的损伤 CartenDolan基于疲劳实验数据建议 4.8 高强度钢 d = 5.8 其它疲劳积累损伤理论的应
5、用例一,某构件可用的S-N曲线为S2N = 2.51010,设计寿命期间的荷载谱如下表,试估计构件不发生疲劳破坏可承受的最大应力水平设计荷载Pi循环数ni (106)P0.050.8P0.10.6P0.50.4P5.0解:设选定构件后的最大应力为S1,各级应力分别为S2 = 0.8S1、S3 = 0.6S1、S4 = 0.4S1,相应破坏寿命 , 按Miner准则估算 解出:S = 151.17MPa 所选构件的最大应力应不超过151.17MPa。例二,某构件的S-N曲线为S2N = 2.51010,如实测一年内所承受的典型应力谱如下表,试估计其寿命。Si (MPa)ni (106)Ni (1
6、06)ni/Ni1500.011.1110.0091200.051.7360.029900.103.0860.033600.356.9440.050解:将典型应力谱作为一个循环块,则各年构件所承受的循环荷载是该典型应力谱的重复,设构件寿命为年,则总损伤 按Miner准则估算 解出: = 8.27年疲劳积累损伤理论解决的两类问题1)已知设计寿命期间的荷载谱,确定应力水平2)已知典型周期内的应力谱,估算使用寿命2疲劳荷载谱处理作用在结构上的荷载按结构的反应可分为静力荷载和动力荷载。动力荷载包括偶然荷载和大部分可变荷载,如风载、波浪荷载、吊车荷载等。疲劳荷载按幅值和频率可分为等幅、变幅和随机荷载。问
7、题:如何将随机荷载谱等效转换为变幅或恒幅荷载谱,以便利用以前的方法处理问题。循环计数法:将不规则的、随机的荷载时间历程转化为一系列循环的方法。雨流计数法适用于以典型荷载谱段为基础的重复历程。雨流计数法方法如下1)由随机荷载谱中选取适合计数的、最大峰或谷处起止的典型段,如图11段(最大峰起止)或22段(最大谷起止)。St典型段11222)将谱曲线旋转900放置并将荷载历程看作多层屋顶,假想有雨滴沿最大峰或谷处开始往下流。如无屋顶阻挡则雨滴反向。42-2-4tS42-2-4tS42-2-4tS3)记下雨滴流过的最大峰、谷值,作为一个循环。图示第一次雨流,循环荷载变程S 5-(-4) = 9,平均荷
8、载 Sm = 5+(-4)/2 = 0.54)从荷载历程中删除流过的部分,对剩余历程段重复以上雨流记数,直至无剩余历程为止。上述雨流法结果如下表循环12345S94732Sm0.510.5-0.5-1荷载谱如是应力,则雨流记数法得到应力变程S和平均应力Sm,因此雨流记数法是二参数记数。与其它记数法相比,雨流记数法的记数结果均为全循环。4不同荷载间的转换记数后的多级荷载可按需要进一步简化为有限的荷载级。荷载间的转换应遵守损伤等效原则。设构件在S1下循环n1次所造成的损伤与构件在S2下循环n2次所造成的损伤相等,由PalmgrenMiner线性疲劳积累损伤准则有 荷载间的转换将造成与真实情况的差别
9、,因此荷载转换次数越少越好。应变疲劳1应力应变关系RembergOsgood模型 稳态循环曲线的数学描述实验结果表明稳态循环曲线中循环应力与塑性应变可用幂函数近似描述 a:循环应力幅值;ap:循环塑性应变幅值;K:循环强度系数;n:循环应变硬化指数。由此得到的稳态循环-曲线的近似表达式 , 2滞后环曲线Massing假设假设滞后环曲线与稳态循环应力应变曲线几何相似,即在a-a坐标系中的a、a分别为D-D坐标系中的D/2和D/2,滞后环曲线为 或: , 满足Massing假设的材料称为Massing材料3变幅循环下的-响应一般,对随机荷载可用典型谱描述,如下例01段,首次加载,应力应变路径为稳定
10、循环应力应变曲线 由此可解出应力112段,1处荷载反向变程,应力应变路径为滞后环曲线 2处的应力应变 , 23段,2处荷载反向变程,D23按滞后环曲线求解,3处的应力应变 , 34段,232形成封闭环,其后的响应按124计算 , 其余类推变幅循环下的响应计算1)第一次加载以循环应力应变曲线描述 2)后续反向以滞后环曲线描述 且: 加载变程用“+”,卸载变程用“”。3应变寿命(-N)曲线按标准实验方法,在R = -1对称循环下得到的应变寿命曲线。荷载以应变幅a表示,寿命以反向次数2N表示。1)将总应变幅分为弹性应变幅ae和塑性应变幅ap , 则lgaelg(2N),lgaplg(2N)呈线性关系
11、,即a2Nae2Nap2Nlg2Nlga2Nt f:疲劳强度系数(应力量纲)b:疲劳强度指数f:疲劳延性系数c:疲劳延性指数一般金属材料,b = -0.06-0.14,估计时取0.1。c = -0.5-0.7,常取-0.6作为典型值。-N曲线(Manson-Coffin公式) 长寿命段,以弹性应变幅ae为主,ap 可忽略 , 或: (幂函数式)短寿命段,以塑性应变幅ap为主,ae可忽略 , 或:2)当ae = ap时有 2Nt:转变寿命寿命大于2Nt,弹性应变为主,属于应力疲劳。寿命小于2Nt,塑性应变为主,属于低周应变疲劳。3)循环性能和疲劳性能各参数的关系循环性能应力应变关系式 , 疲劳性
12、能应变寿命关系式 , 有: 对比有:, 各系数由实验结果拟合得出,故数值上往往不能严格满足。4-N曲线的近似估计应变控制下一般金属材料的-N曲线有下图所示特征,即在应变 = 0.01处有大致相同的寿命。Manson经验公式(仅对恒幅对称应变循环) f:材料断裂真应变 RA:断裂时面积收缩率 5平均应力修正Morrow弹性应力线性修正公式 m:平均应力Gerber弹性应力曲率修正公式 Sachs弹性修正公式 例:已知构件的E = 210103MPa,n = 0.2,K = 1220MPa,f = 930MPa,b = -0.095,c = -0.47,f = 0.26。试估计构件在下图所示三种应
13、变历程下的寿命。0.020.005-0.005-0.02t0.020.005-0.005-0.02t0.020.005-0.005-0.02t解:1)恒幅应变对称循环 求出:N = 5858次 2)计算-响应 01段, 1 = 542MPa 12段,D12 = 972MPa , 12, 43 23段, 3 = 0.005,3 = 342MPa 34以后形成封闭环 估算寿命 ,N = 6170次 拉伸后引入了残余压应力,提高了寿命12, 43 3)类似地,34以后形成封闭环 估算寿命 N = 5565次 压缩后引入了残余拉应力,降低了寿命构件应变疲劳的寿命与加载历史有关。内容总结(1)第四讲 损伤理论,荷载谱处理,应变疲劳上节概述p-S-N曲线,失效率,存活率疲劳数据处理,正态分布,标准正态分布检验数据是否服从正态分布,正态分布坐标纸正态分布存在的问题:不能反映构件疲劳寿命有一个大于等于零的下限(2)与其它记数法相比,雨流记数法的记数结果均为全循环16