1. 课程导入与工程问题界定
1.1 从一次断裂谈起:为什么需要损伤本构
某石化装置的一段主蒸汽管道在服役约8万小时后发生爆裂。事后检查发现,管壁并未出现明显减薄,外表面也无可见裂纹,但断口呈脆性、无明显塑性变形。金相分析显示,材料晶界处已聚集大量蠕变微孔洞,局部孔洞体积分数接近3%。这一案例揭示了一个关键事实:材料的失效往往不是从宏观裂纹开始的,而是从内部微缺陷的萌生、长大与聚合逐步演化而来。传统的强度校核——只比较工作应力与屈服强度或许用应力——无法解释这种"无预兆"的破坏。
损伤力学正是为描述这一过程而发展起来的。它引入一个连续场变量D(损伤变量),把微孔洞、微裂纹的分布效应"抹平"到连续介质框架中。当D=0时材料完好,D=1时材料完全失去承载能力。损伤本构方程则把D的演化与应力、应变、温度、时间等可测量联系起来,使"失效演化"从定性描述变为可计算的工程问题。
1.2 三类典型失效案例与损伤机制
工程中的损伤失效可归纳为三种典型模式,它们对应不同的微观机制和本构描述方式。
案例一:蠕变损伤(高温承压部件)。在550℃、内压16MPa的蒸汽管道中,晶界孔洞在持续应力下形核、长大。损伤演化常用Kachanov-Rabotnov形式描述:dD/dt = A·(σ/(1-D))^n,其中A、n为材料常数。当D从0.1增长到0.3时,有效承载面积下降,蠕变速率显著加快,最终进入第三阶段并断裂。这类问题的关键是时间—温度—应力三参数耦合。
案例二:疲劳损伤(循环载荷部件)。涡轮盘榫槽、桥梁拉索等承受交变载荷。损伤以微裂纹形式在应力集中区萌生,D随循环数N累积。常用Chaboche非线性疲劳损伤模型:dD/dN = [σ_a/(M(1-D))]^β,其中σ_a为应力幅。与线性Miner法则不同,该模型能反映载荷顺序效应——高—低加载比低—高加载更危险。
案例三:韧性损伤(冲压、碰撞部件)。汽车B柱在碰撞中经历大塑性变形,微孔洞在夹杂物处形核、长大、聚合。Gurson-Tvergaard-Needleman(GTN)模型用孔洞体积分数f作为损伤变量,当f达到临界值f_c时孔洞聚合,承载能力骤降。这类问题的关键是应力三轴度对孔洞长大的强烈影响。
1.3 课程知识地图:从变量到决策
本课程共30章,按"基础—模型—实现—应用"四层递进组织。第一层(第1~6章)建立损伤变量、有效应力、热力学框架等基本概念,回答"损伤如何定义、如何测量";第二层(第7~16章)分别展开蠕变、疲劳、韧性、蠕变—疲劳交互等损伤本构模型,回答"不同机制下D如何演化";第三层(第17~24章)进入数值实现,包括隐式积分算法、切线刚度推导、UMAT/VUMAT编写与验证,回答"如何把本构方程变成可运行的代码";第四层(第25~30章)面向工程应用,涵盖参数标定、加速试验设计、寿命评估与可靠性分析,回答"如何用模型做决策"。
贯穿全课程的一条主线是参数标定与验证。损伤本构方程中的材料常数(如A、n、M、β、f_c)不能凭经验随意选取,必须通过标准试验(蠕变持久试验、低周疲劳试验、缺口拉伸试验)结合反演优化确定。例如,蠕变损伤参数n通常由不同应力水平下的持久寿命数据拟合得到:在双对数坐标中,log t_r与log σ呈线性关系,斜率即为−n。这一过程将在第25章详细展开。
1.4 工程问题界定:边界与假设
在进入具体模型之前,必须明确本课程的工程问题边界。我们讨论的对象是金属材料在连续介质尺度(特征尺寸远大于微缺陷间距)下的损伤演化。这意味着:不涉及原子尺度位错动力学,不涉及离散裂纹的扩展路径追踪(那是断裂力学的范畴),也不涉及复合材料层间损伤等特殊问题。损伤变量D或f是唯象的连续场变量,其演化方程基于热力学耗散势构建,而非从微观第一性原理推导。
另一个关键界定是准静态与小变形假设。本课程大部分模型建立在小变形框架下,应变与位移呈线性关系。对于大变形问题(如碰撞、冲压),需采用有限变形框架下的损伤本构,这部分将在第16章简要介绍。此外,温度场假设为均匀或缓变,不涉及热冲击引起的瞬态损伤。
明确这些边界不是限制,而是为了让学习者在正确的范围内使用工具。当你面对一个实际工程问题时,第一步永远是判断:它属于蠕变、疲劳还是韧性损伤主导?是否有多机制交互?变形是否超出小变形范围?只有界定清楚,才能选对模型、标对参数、给出可信的寿命预测。