时间相关的断裂理论及应用()

时间相关的断裂理论及应用

Time-dependentfracturemechanicsandapplication华東理工大学EastChinaUniversityofScience&Technology轩福贞机械与动力工程学院2014年3月17日主要内容什么是蠕变（时间相关的断裂）工程中的蠕变现象蠕变变形及断裂理论高温结构的剩余寿命分析目前的问题及挑战作业22023/1/31蠕变的定义蠕变定义在服役温度（35-70%熔点）下，应力远低于屈服应力时结构/材料中长度（变形）随时间增长而增加，卸载后变形无回复；破坏时无明显颈缩蠕变过程及失效机理42023/1/31ExplanatorypictureoftensiletestbyGalileoTestingbarWeight蠕变现象蠕变过程不仅表现为外形尺寸的改变，而且内部衍生孔洞，聚合长大导致部件/材料的断裂Pyzalla,Camin.SCIENCE2005,308高能束同步辐射(HESR)观察下的铜合金蠕变损伤过程

5蠕变定义与短时拉伸曲线的区别：低应力，时间长CreepstrainStressrupturecreepruptureTransientcreepTimeSteady-statecreepacceleratingcreep62023/1/31蠕变的特征两面性：变形与强度变形随时间温度增加而增大72023/1/31蠕变的特征两面性：变形与强度Creeptemperatureregion;TemperatureregionatwhichCreepmustbeconsideredNon-creeptemperatureregionCreeprupturestrengthShort-termtensilestrengthTemperature破断强度Creepdesignofstructuralcomponent强度随时间温度增加而降低82023/1/31蠕变的特征2023/1/319时间增加强度下降2023/1/3110Rupturestress应力准则TimeFailurepointσ0σdestdes蠕变的特征2023/1/3111发生蠕变的温度蠕变三要素温度、应力/载荷、时间蠕变现象三特征failuretertiarycreepd/dt

primarycreepd/dtAndradelawtimestrainsecondarycreepd/dt=constanttf()stresstimetimed/dtinelasticelastic,reversible应力恒定应变三阶段：增加、稳定、增加应变速率三阶段：减速、稳定、增速蠕变的特征蠕变研究的历史2023/1/31121830年，法国工程师Louis-JosephVicat最早报道了桥梁绳索蠕变超过弹性变形导致失效的现象1910年，英国伦敦大学EdwardAndrade教授最早建议了描述金属材料的粘性方程1922年，英格兰工程师Dickenson发现了高温蠕变强度低于短时拉伸强度的现象Riedel教授于1977年提出了蠕变条件下裂纹张开与扩展Dugdale模型Harper同年提出了C*断裂参量用于蠕变研究美国西屋公司最早编制汽轮机蠕变裂纹扩展寿命评价标准2023/1/3113蠕变研究的历史航空工业的发展出尽了高温断裂问题的研究工程中的蠕变现象2023/1/3114完美结构产品三要素：需体现成本效率、安全和经济寿命的理想匹配产品三要素服役结构的两类非正常现象ProcessoffractureLowpressurerotorMiddlepressurerotor美国TVAGallatin电站转子服役早期发生断裂美国Sabine热电厂2号机组(1979)主蒸汽管道弯头5万小时蠕变破裂1.结构早于设计寿命发生破坏服役结构的两类非正常现象英国1962年建造Parsons汽轮机组超过20万小时设计寿命10余年，目前仍然安全服役大庆石化总厂炼油厂F－401加氢反应器1966年投入运行至今，超过了设计寿命30余年2.结构安全服役超过数倍的设计寿命原因：材料？制造工艺？对蠕变问题认识不足！新技术使得温度压力日益提升涂善东，高温结构完整性原理，科学出版社，2003TuST,HighTemperatureStructuralIntegrity,SciencePublisher,2003蠕变失效的宏微观表现蠕变失效宏观形式□鱼骨状含Nb碳化物骨架状含Cr碳化物新材料使用16万小时航空发动机叶片蠕变断裂2023/1/3120焊接结构的蠕变10万小时的焊缝蠕变强度削弱系数不同材料，不同温度，对应不同寿命下，接头的蠕变强度减弱焊接结构的蠕变断裂特点晶界为蠕变敏感区，低应力破坏源于热影响区的细晶粒区

SchubertJ.,A.KlenkandK.Maile.IntConfCreepFractureinHighTemperatureComponents,ECCC),2005传统的低合金铁素体钢WSRF=0.8蠕变强化钢9-12%CrWSRF=0.6蠕变断裂的复杂性2023/1/3123蠕变变形及损伤理论2023/1/3124蠕变研究的任务蠕变变形的本构方程蠕变断裂参量及裂纹扩展预测2023/1/3125蠕变变形的描述方程26Norton，1929McVetty，1943Dorn，1955Evansetal，1985投影法Kachanov-Rabotnov，蠕变损伤，1968蠕变过程中的内部变化27GrainelasticdeformationSubgrainformationGrainboundaryslipGrainplasticdeformation三种主要变形过程：

晶粒变形、晶界滑移、亚晶粒形成蠕变变形的机理28位错蠕变（Dislocationcreep）或Powerlawcreep蠕变速率表达式Deff：EffectivediffusioncoefficientStressindex;n≥3DislocationClimb析出粒子,precipitate平移平移上升上升蠕变变形的机理29扩散蠕变（diffusioncreep）Latticediffusion(Nabarro-Herring)creep:VacanciesmoveinthelatticeGrainboundarydiffusion(Coble)creep:VacanciesmovealongthegrainboundaryStressindex;n=1,

Graindiameterindex;p=2Graindiameterindex;p=3CreepdeformationmechanismmapofNickel,Grainsize;1mm电厂多部件受到蠕变-疲劳交互作用的影响31超超临界汽轮机关键部件高温反应堆反应堆压力容器、管路系统、热交换器等火电设备锅炉、管路系统、汽轮机转子、汽缸、螺栓等联合循环设备燃气轮机的热部件何谓交互作用？二者叠加不等于线性之和32Creep-fatigueinteractionFatigueaftercreepPurefatigue载荷谱的简化33e0ExampleforloadspectrumTCTCCreepfatigue蠕变疲劳分析的基本要素34蠕变疲劳交互作用的影响因素保载时间：会导致循环次数下降加载速率效应、平均应力控制模式：应力？应变？热机疲劳效应（TMF）多轴应力的影响环境的影响（氧化、腐蚀等）35蠕变-疲劳交互作用导致独特的破坏机制断口裂纹扩展路径晶界316L600度疲劳蠕变-疲劳37蠕变-疲劳交互作用存在敏感窗口FatiguedominatedFatigue-creepinteractionCreepdominated蠕变-疲劳交互作用研究的历程基于Coffin-Manson方程的蠕变疲劳模型发展历程现有蠕变疲劳主要预测方法DuctilityexhaustionapproachStrainsensitivematerialsLineardamagesummationDefinitionof“creepductility”TruerupturestrainReductionofareaRuptureelongation39Timefractionmethodassumption:twodamageprocessesareindependent: -CreepDamageisaninternalprocessdominatedbycavitationsatgrainboundaries -Fatigueprocessisasurfacephenomenon标准中的简化处理规则40多轴应力的影响2023/1/3141应力状态对变形的影响NuclearEngineeringandDesign,2006,237:1969产生同样变形所需应力更大Ds多轴应力的影响2023/1/3142marginεdesDesignlifemargin基于分析（FE）的蠕变强度设计结构几何与载荷参数（温度、压力）材料行为如NortonBalley方程FE蠕变分析蠕变强度校核危险点应力≤蠕变破断强度危险点应变≤蠕变破断韧度核心问题本构关系和材料参数核心问题失效判据高温结构的蠕变断裂理论44蠕变裂纹扩展45

高温下裂纹扩展额特征裂尖蠕变区域及应力场蠕变区蠕变区蠕变区（b）瞬态蠕变c）稳态蠕变（a）小范围蠕变蠕变裂尖应力场平面应变

平面应力

Goldman和Hutchinson指出，在按幂指数硬化的对速率不敏感材料中，最后得到的方程与控制渐进特性的方程具有相同的形式，即应变与应力的奇异性是HRR型的，相应的渐进场可写为裂尖强度因子随时间的变化2axyoC*的定义式中

为应变能速率的密度

对服从Norton律的材料：

稳态蠕变下C*与路径无关蠕变裂纹行为Goldman和Hutchinson指出，在按幂指数硬化的对速率不敏感材料中，最后得到的方程与控制渐进特性的方程具有相同的形式，即应变与应力的奇异性是HRR型的，相应的渐进场可写为C(t)的定义注意：随着时间的增加，裂尖的小范围蠕变也逐渐演变成大范围蠕变，C(t)即等同于C*并与积分路径无关。小范围蠕变时裂尖的应力应变场定义裂尖蠕变区域的尺寸：(s)HRR=(s)elastic

蠕变区域的尺寸为

C(t)的另一表示在平面应变条件下

从小范围蠕变向大范围蠕变的过渡时间可以表示为

小范围蠕变下的Ct载荷蠕变位移率试样厚度试样半宽度

F是q和n的函数，具体变化如图（RiedelandRice,1980）大范围蠕变下的Ct与试样尺寸（a/w）以及蠕变指数n有关。具体数值可参考A.Saxena等人的的著作。C*的获取当构件中裂纹前端的韧带区都进入稳态蠕变阶段（蠕变第二阶段）时，弹性的应变率将消失，材料显示出非线性粘性流的特性。在这种情况下，C*值与积分路径无关，并且可用于裂尖应力应变场的表征。因此如何在理论和实验中获取该值，对蠕变断裂力学具有重要意义。有限元求C*Shih等人在虚裂纹闭合研究的基础上，发展了所谓的域积分技术，并且这一积分格式已在商用有限元软件（如ABAQUS）中得到了应用通过实验求C*半经验法求C*其中P为加载点载荷，为加载点位移速率，B、W分别为试样的厚度与宽度，h为a/W,n,n1的函数。

CT试样的η值可用下式表示参考应力法求C*

Ainsworth提出PL为塑性破坏载荷

为参考应力作用下的蠕变应变率，R为长度参量，可通过有限元得到。高温结构安全评价方法问题蠕变损伤模型试验断裂力学参量经验公式的准确性无法保证缺乏高温下的预测模型传统模型参数繁多，难以标定多轴应力状态对蠕变本构和损伤方程的影响亟待修正耗时费力观测困难全寿命预测模型损伤断裂全寿命预测模型高温结构剩余寿命分析2023/1/3163实际结构常见的载荷形式2023/1/316465单纯疲劳载荷下的裂纹扩展数学模型大量的金属疲劳裂纹扩展实验表明，裂纹扩展数率与裂纹尖应力应变场的控制参量（应力强度因子波动幅度K）有关。描述疲劳裂纹扩展速率的表达式中，最简单、应用最广泛的是Paris于1961年提出的表达式：C0和n0分别被称为裂纹扩展系数和裂纹扩展指数，是与材料有关的常数。△KⅠ表示应力强度因子幅。66单纯疲劳载荷下的裂纹扩展数学模型应力强度因子幅的表达式为：

Q—为缺陷形状参数，包含了塑性引起的修正

Y—为几何结构修正系数，对于内部椭圆形裂纹（埋藏缺陷）Y=1，对于表面椭圆裂纹（表面缺陷）Y=1.12。67单纯疲劳载荷下的裂纹扩展数学模型裂纹的临界裂纹尺寸的计算公式为：—为转子材料的断裂韧性

—为裂纹（或缺陷）尖端的公称应力

所以转子的疲劳寿命的表达式为：68在蠕变、疲劳交互作用下裂纹扩展模型为：

上式右边第一项是由交变载荷引起的裂纹扩展率，第二项为由蠕变及蠕变、疲劳交互作用引起的裂纹扩展率。考虑到汽轮机的启动时间比带负荷时间要小的多，可以不计蠕变对疲劳裂纹扩展的影响，但交变载荷对蠕变裂纹扩展的影响需要考虑。也就是说，由于不计蠕变对疲劳裂纹扩展的影响，第一项可以利用材料的低周疲劳结果进行计算。蠕变－疲劳交互作用下的裂纹扩展数学模型69蠕变－疲劳交互作用下的裂纹扩展数学模型(Ct)ave参量可利用下式来表示：工程中仅考虑稳态蠕变阶段（第二阶段蠕变）下的计算公式70蠕变－疲劳交互作用下的裂纹扩展数学模型B和q为材料常数，对于CrMoV钢有：B=0.071,q=0.75

对上式进行积分就可得出循环次数为：

71多裂纹情况的干涉与合并（表面裂纹）当-0.1<S/D<0.4且h/a<0.8时，可以将这两个裂纹合并为一个裂纹；否则就不合并。S=2b-4a,s/2b=1-2a/b=s/D.其中的2a表示的表面裂纹的长度，也就是简化为半椭圆后的2c,在这里取两个裂纹中较大的那条裂纹长度.D表示这两个裂纹总长度.寿命评定流程2023/1/3172202