1. 应力与应变到底在说什么——基本概念解构1.1 应力的定义单位面积上的内力先聊一个最基础的问题一根钢丝绳吊着200公斤的重物钢丝绳内部到底发生了什么变化直觉告诉我们它被“拉紧了”。但“拉紧”这个词太模糊工程师需要的是精确描述拉紧到什么程度才算安全多粗的钢丝绳才够用答案是——应力。应力的物理含义是单位面积上的内力公式非常简单[ \sigma \frac{F}{A} ]其中F是截面上的内力A是受力面积。但这里有一个初学者最容易忽略的点应力不是直接测量出来的它是我们“人为定义”的一个物理量。真实世界里你能测到的只有力、位移、变形而应力是通过“截面法”在脑海中把一个物体切开用平衡条件反推出来的内力分布结果。截面上任意一点的应力又可以分解成两个分量垂直于截面的叫正应力用σ表示通常拉应力为正、压应力为负平行于截面的叫切应力用τ表示。为什么要这样分解因为材料的破坏模式不一样正应力会把材料拉断或压碎切应力会让材料沿截面“错开滑移”。一个螺栓的断裂可能是拉断的也可能是剪切破坏的不区分这两个量你根本没法判断失效原因。再说一个工程习惯上的小细节应力单位在国际单位制中是帕斯卡Pa但实际工程中大家几乎都用兆帕MPa或牛每平方毫米N/mm²。这两个单位是等价的1 MPa 1 N/mm²只是叫法不同。我看过不少新人把MPa换算成N/m²的时候少了一个数量级结果整个强度校核全部出错这个坑后面专门说。1.2 应变的定义变形程度的度量应力描述的是“内部受力状态”但受力之后物体会怎样会变形。应变就是用来度量变形程度的物理量。线应变的定义是长度变化量与原始长度的比值[ \varepsilon \frac{\Delta L}{L} ]这个量是无量纲的但在工程图纸和规范里常用微应变με来表示1 με 1×10⁻⁶。比如一根100毫米长的钢杆拉伸后变为100.05毫米应变量就是0.05/100 5×10⁻⁴也就是500微应变。还有一个概念叫切应变用γ表示它描述的是原本互相垂直的两条微小线段在变形后夹角的变化量单位是弧度。注意这里有个易混点切应力τ对应的是切应变γ正应力σ对应的是线应变ε千万不要把对应关系搞反。我在检查报告时经常看到有人把剪切问题用正应力公式去算结果偏差几十个百分点。实际工程中应变可以直接用电阻应变片测出来。应变片的原理很巧妙金属丝拉伸时电阻变化通过惠斯通电桥把电阻变化转成电压信号标定后换算成应变值。结构试验、桥梁监测、飞机机翼加载测试用的都是这一套。明白“应变是直接可测的应力是间接算出来的”这个事实对理解后面的本构关系非常有帮助——因为应力-应变曲线实验里纵坐标的应力其实也是通过力传感器测量后除以面积得到的横坐标的应变才是真正的“原始数据”。1.3 六个应力分量与六个应变分量空间应力状态的完整描述到了三维空间里情况比单轴拉伸复杂得多。物体内部任意一点实际上是无限小的一个微元体的应力状态需要六个独立的应力分量才能完整描述三个正应力分量σx、σy、σz以及三个切应力分量τxy、τyz、τzx。为什么六个就够了因为切应力互等定理τxy τyxτyz τzyτzx τxz原本看起来有九个分量实际上对称后只剩六个独立量。这个定理不是数学游戏它是微元体的力矩平衡条件推出来的——如果微元体两侧的切应力不相等微元体就会发生旋转。对应的也有六个应变分量εx、εy、εz、γxy、γyz、γzx。这个“分量体系”经常让初学者头大为什么不能像中学物理那样只用一个力、一个位移来描述打个比方你在一个房间里推一个柜子柜子同时受到来自前、后、左、右、上、下各个方向的力每个方向的受力情况可能都不一样。要完整描述柜子的受力状态只报告“它受力了”显然没用必须知道每个面的受力大小和方向。应力分量就是干了这件事。刚接触有限元软件的人最容易忽视这个分量体系——在ANSYS或Abaqus里查看结果时默认显示的往往是最直观的“冯·米塞斯等效应力”但这个值是六个分量综合折算出来的标量它丢失了方向信息。真正做结构校核时正应力和切应力要分开看分别在对应的强度条件里去比对只盯一个米塞斯应力是不够的。2. 从单轴到三轴——应力状态的全貌2.1 单轴拉伸与杆件的基本假设材料力学课本第一课基本都从单轴拉伸试验讲起一根圆棒试件夹在万能试验机上缓慢加载记录力-位移曲线。在这个最简单的受力状态下试件横截面上的应力可以认为是均匀分布的这也是材料力学最基本的平截面假设——变形前为平面的截面变形后仍然保持为平面。这个假设极其重要它是工程梁理论、杆件应力公式的基石。但要注意它的适用边界只有在“远离端部”的截面区域才成立。在夹头附近、在截面突然变化的区域应力分布根本不均匀这也就是后面要讲的圣维南原理和应力集中问题。单轴拉伸试验给出了一条最重要的曲线工程应力-应变曲线。横坐标是应变ε纵坐标是应力σF/A₀A₀为初始截面积。这条曲线可以分为几个典型阶段弹性阶段、屈服阶段、强化阶段、颈缩阶段。钢材、铝合金、铜合金都符合这个基本规律而铸铁、陶瓷等脆性材料则没有明显的屈服阶段直接断裂。这里必须强调一个很多人栽过的坑工程应力用的是初始面积A₀而不是瞬时面积Ai。在弹性阶段两者差别不大但进入塑性阶段特别是颈缩之后试件实际截面积急剧减小真实应力F/Ai实际上还在上升而工程应力反而在下降。所以曲线末段的“下降”是假象材料本身的抗力并没有下降。做材料成型仿真时一定要用真实应力-应变数据否则回弹和破裂预测会偏得离谱。2.2 切应力互等定理为什么切应力总是成对出现很多教材把切应力互等定理当作一个简单的推导结论提一下但它在工程中的意义远比表面看起来大。想象一个微元体只在上表面作用一个向右的切应力τ。如果下表面没有反向的切应力这个微元体就会受到一个力矩而开始旋转。为了保证微元体平衡不转动必须同时存在左侧面向下的切应力和右侧面向上的切应力形成反向力偶来抵消。这就是切应力互等定理的本质——切应力总是成对出现在互相垂直的平面上大小相等、方向指向或背离交线。这个定理的直接推论是如果你在一个受扭的圆轴上测得某点的切应力方向是某个角度那么与该点相邻的另一个垂直截面上必然存在对应的切应力。很多做机械设计的同行在校核键槽或花键连接时只算了主剪切面上的应力忘了检查与之垂直的截面结果疲劳裂纹从另一个方向萌生这就是对切应力互等定理理解不深。2.3 主应力与最大切应力斜截面上的应力变化单轴拉伸时如果你切一个垂直于轴线的截面截面上只有正应力σ没有切应力。但如果你切一个与轴线成45°的斜截面呢斜截面上的正应力和切应力都会发生变化。经过推导可以得到单轴拉伸状态下与轴线夹角为α的斜截面上的应力为[ \sigma_\alpha \sigma \cos^2\alpha ] [ \tau_\alpha \frac{\sigma}{2} \sin 2\alpha ]当α45°时切应力达到最大值τmax σ/2。这个结论解释了铸铁拉伸试验中一个经典的宏观现象——断口并不是垂直于轴线的平齐断面而是呈约45°的斜茬。为什么因为铸铁的抗拉强度远低于抗剪强度但最大主应力先达到极限时发生的是正断而韧性材料扭转时最大切应力首先达到剪切屈服强度所以表现为沿横截面的剪切破坏。每次看到有人拿一个扭转断口来问“为什么是螺旋状断裂的”我都会让他回去翻一下这个公式——45°方向上的最大拉应力或最大切应力决定了一切。由此可以引申出更一般的问题一个复杂应力状态下过同一点的不同截面上的应力是不同的。那到底哪个截面上应力最大哪个方向的应力最危险这就引出了主应力的概念——三个互相垂直的截面上切应力为零这三个截面上的正应力就是主应力σ1、σ2、σ3按代数值大小排列σ1 ≥ σ2 ≥ σ3。求主应力的工具叫摩尔圆。在σ-τ坐标系里画一个圆圆上每个点对应某一个截面的应力状态圆心的位置是平均正应力(σxσy)/2半径是最大切应力。摩尔圆用一个圆把平面应力状态的全部信息都囊括了你根本不用背主应力公式画个圆量一下就行。做混凝土结构抗裂验算和岩土工程边坡稳定分析时摩尔圆是每天的必备工具。2.4 平面应力与平面应变两类简化模型的工程取舍实际结构大多数是三维的但直接做三维分析代价太大——计算量大、建模复杂、结果解读困难。工程上常用两类简化模型平面应力问题和平面应变问题。平面应力适用于薄板类结构板厚远小于其他两个方向尺寸比如薄壁容器壳体、压力容器封头。因为板面自由厚度方向的正应力σz 0但厚度方向的应变εz一般不为零——板会变薄或变厚。平面应变适用于横截面不变的细长结构或厚壁结构比如水坝、隧道衬砌、长管道。因为沿长度方向被约束长度方向的应变εz 0但长度方向的正应力σz不为零。选择错误的模型会导致严重误判。举个例子有人把一段厚壁高压管道按平面应力去算算出来壁厚方向的应力为零觉得特别安全实际上厚壁管道内压作用下的环向应力分布是内壁大、外壁小而且存在径向应力梯度平面应力完全无法描述。选模型时先问自己一个问题这个结构的厚度方向是自由变形还是被约束的自由的用平面应力被约束的用平面应变这个是原则性判断。3. 应力与应变的桥梁——胡克定律3.1 胡克定律线弹性阶段的应力-应变关系材料在弹性阶段内应力和应变成线性关系这就是胡克定律[ \sigma E \cdot \varepsilon ]E是弹性模量杨氏模量单位与应力相同MPa或GPa它表征材料抵抗弹性变形的能力。某种意义上E就是材料“刚性”的量度E越大同样的应力产生的应变越小。钢材的E约为206 GPa铝合金约为70 GPa混凝土约为30 GPa。每次做结构选材时我都会先看这个数——很多刚入行的设计人员只顾着比强度忘了比刚度结果做出的结构强度富余很多但变形超标返工时被折腾得够呛。剪切工况下有对应的剪切胡克定律[ \tau G \cdot \gamma ]G是剪切模量它与弹性模量E和泊松比ν之间存在一个固定的关系[ G \frac{E}{2(1\nu)} ]这个公式不是随便凑的它是由各向同性弹性体的张量性质推导出来的。对各向同性材料三个弹性常数E、G、ν中只有两个是独立的。也就是说你只要从拉伸试验测出E和ν剪切模量G可以算出来不用单独做扭转试验。做复合材料的同行要特别注意各向异性材料的独立弹性常数数量完全不同正交各向异性有9个独立常数上面的关系不成立别套错公式。3.2 泊松比拉伸时为什么会变细拉伸一根橡皮筋它不只是变长还会明显变细。这个横向收缩效应用泊松比ν来描述[ \nu -\frac{\varepsilon_{\text{横向}}}{\varepsilon_{\text{纵向}}} ]绝大多数金属材料的ν在0.25到0.35之间钢材约0.3橡胶接近0.5不可压缩软木接近0。泊松比的意义容易被低估。它直接影响三类问题第一配合公差计算——过盈配合的轴和孔压装后内孔的径向变形估算需要用到ν第二密封设计——法兰螺栓预紧后垫片的横向膨胀量与ν直接相关第三热应力分析——温度变化引起的约束应变会产生内应力而ν决定了这个内应力的大小。这里有个很经典的单轴拉伸扩展公式单向拉伸时轴向应变为εx σ/E而横向应变εy εz -ν·σ/E。很多初学者算薄板单轴拉伸的宽度变化时忘了乘负号或者忘了ν结果方向都对、数量级错了。3.3 广义胡克定律多向应力状态下的本构关系实际结构很少是纯粹的单向受力。一个压力容器筒体内压作用下环向应力大约是轴向应力的两倍这就是一个双向应力状态。这种情况下x方向的应变不仅由σx产生还受σy、σz的影响具体的广义胡克定律为[ \varepsilon_x \frac{1}{E}[\sigma_x - \nu(\sigma_y \sigma_z)] ] [ \varepsilon_y \frac{1}{E}[\sigma_y - \nu(\sigma_x \sigma_z)] ] [ \varepsilon_z \frac{1}{E}[\sigma_z - \nu(\sigma_x \sigma_y)] ]表面上看只是叠加了个ν的项但物理意义深刻。一个方向的应变是三个方向正应力共同作用的结果——这就解释了为什么压力容器爆破试验时筒体的周向应变和轴向应变的组合关系可以反推出材料真实的弹塑性特性也解释了为什么做一个只受x方向拉伸的试件如果约束了y方向的变形y方向会产生拉应力因为材料想收缩但被挡住了。回到工程应用上电阻应变片测出的是应变而强度校核需要的是应力。在复杂应力状态下把应变直接乘E换算成应力是错误的必须根据具体的边界条件判断是单轴、平面应力还是平面应变状态再选择对应的广义胡克定律公式。我在做桥梁健康监测时看到过不少报告用单向胡克定律处理双向应变数据算出的主应力方向都是错的这个教训值得记下来。3.4 真实应力-应变曲线屈服、强化与颈缩前面讲的胡克定律只适用于线弹性阶段。真实材料在进入塑性阶段后应力-应变关系变成曲线此时再用σ E·ε去求应力就完全不适用了。典型金属材料的拉伸应力-应变曲线有几个关键点需要记牢比例极限应力-应变保持线性关系的最大应力值弹性极限卸载后不留下残余变形塑性变形的最大应力值屈服极限材料开始明显产生塑性变形的应力值。低碳钢有明显的屈服平台上屈服点和下屈服点工程上常取下屈服点作为设计依据抗拉强度材料承受的最大工程应力值设计规范中塑性材料用屈服极限作基准取安全系数脆性材料用抗拉强度作基准。这个概念不清楚的话后面的强度校核全盘皆错——有人拿抗拉强度去校核塑性材料的许用应力结果安全系数取小了结构在屈服阶段就倒了。关于“工程应力”和“真实应力”还有个常见认知盲区。前面说过工程应力用初始面积计算真实应力用瞬时面积计算。在弹性阶段两者几乎重合但颈缩开始后分道扬镳。做冲压成型、轧制、锻造这类大变形工艺仿真时材料卡片里必须输入真实应力-应变数据否则仿真结果会严重失真。这个问题的本质是有限元软件里的本构模型需要的物理变量是柯西应力真实应力和真实应变对数应变不是工程应力-工程应变。只要看到材料库里用的是“塑性应变-真实应力表”就说明你需要的是真实应力曲线。4. 工程中那些最容易踩的坑4.1 应力集中为什么孔边最先开裂开一个小圆孔的无限大平板在远处施加均匀拉伸应力σ₀孔边最大应力不是σ₀而是3σ₀。这个“3”就是圆孔的理论应力集中系数Kt。应力集中的物理本质是力线在几何突变处“拥挤”流线类比最容易理解河流在狭窄处流速加快应力集中处则应力升高。弹性力学理论解表明圆孔边缘切向应力分布不均匀孔边最大切向应力出现在垂直于加载方向的孔壁顶点。但要注意Kt是针对理想弹性材料而言的。如果材料是塑性的孔边局部应力超过屈服点后会发生应力重分布峰值应力被“削平”实际应力集中效应会减弱。这就是为什么塑性材料比脆性材料对孔洞、缺口“宽容”的原因。反之铸铁、玻璃、陶瓷等脆性材料对应力集中极其敏感灰铸铁拉伸试验断口往往从铸造缺陷或截面突变处起始就是这个原因。工程实际中处理应力集中的常用手段包括圆角过渡代替尖角、开卸荷孔、调整截面形状使力线平滑流动。做轴类零件设计时轴肩处一定要加过渡圆角半径尽量大这是起码的机械设计素养。焊接结构中焊缝余高和母材的过渡区域天然就是应力集中源所以焊缝外观质量对疲劳寿命的影响极大。4.2 名义应力与真实应力的换算静不定问题也要算名义应力就是用“理论公式/截面法”算出来的应力比如拉伸应力σ F/A、弯曲应力σ M·y/I。但是真实结构中因为存在应力集中、残余应力、装配应力等因素实际应力往往高于名义应力。经典例子螺栓连接的强度校核。名义上螺栓只承受纯拉伸载荷但实际螺牙根部存在严重的应力集中加上螺纹制造过程的表面残余应力实测疲劳强度只有光滑试件疲劳强度的几分之一。设计规范里给出的螺栓许用应力值都综合了这些折减因素直接套用理论公式反而危险。还有一个经常被忽略的情况是热应力。一根两端完全固定的钢杆温度升高50°C如果材料膨胀被完全约束内部会产生多大的应力算一下ε α·ΔT 12×10⁻⁶ × 50 6×10⁻⁴σ E·ε 206 GPa × 6×10⁻⁴ ≈ 123.6 MPa。这个应力已经相当可观了。管道工程中如果不设置膨胀节、U形弯等补偿器管道热胀冷缩产生的应力足以把支架推倒、把法兰接头拉裂。4.3 应力单位换算的惨痛教训这个坑值得单独拿出来讲因为我在实际评审中真的遇到过因为单位换算出错的严重事故报告。常见错误链条是这样的计算书里用了kN和cm为单位套用了公式σ F/A算出来的应力单位是kN/cm²可规范里允许应力用的是MPaN/mm²中间差了10倍。换算关系是1 kN/cm² 10 N/mm² 10 MPa这个10倍只差一步但出错的人不在少数。更隐蔽的问题是压力单位1 MPa 1 N/mm² 10 kgf/cm²粗略换算。有人把液压系统的25 MPa当成25 kgf/cm²去算油缸缸筒壁厚结果多出近一个数量级的误差。规避这个问题的办法只有一个所有计算书统一使用一套单位制。我个人的习惯是长度一律用mm、力一律用N这样应力直接就是MPa不需要做任何中间换算。如果是用美标或欧标的图纸先量纲分析一遍再做数值计算。4.4 强度校核时该选哪个“应力”——失效准则入门一个零件受到复杂的多向应力怎样才能判断它是否安全不能只拿最大正应力去和许用应力比因为不同的应力组合会产生不同的失效模式。工程上常用四个经典强度理论第一强度理论最大拉应力理论认为只要最大主应力σ1达到材料极限就失效。适用于脆性材料拉伸断裂比如铸铁。第二强度理论最大伸长线应变理论认为最大伸长线应变达到极限时失效。适用于脆性材料的压缩-拉伸组合工况但工程中已较少使用。第三强度理论最大切应力理论特雷斯卡准则认为最大切应力τmax (σ1-σ3)/2是失效判据。适用于塑性材料的屈服。计算偏保守是传统机械设计的主力。第四强度理论畸变能密度理论冯·米塞斯准则认为材料屈服与畸变能密度有关得到等效应力σe √[((σ1-σ2)²(σ2-σ3)²(σ3-σ1)²)/2]。适用于塑性材料与实验结果吻合度最高是现代有限元分析和压力容器设计的主流判据。为什么有限元软件里默认显示冯·米塞斯应力因为大部分结构材料是金属塑性材料而米塞斯准则能较好预测金属屈服。但如果你在分析脆性材料玻璃、陶瓷、铸铁、或者分析带裂纹构件的断裂行为米塞斯应力就不是合适的判据了。每次做分析前先确认材料类型再选择合适的失效准则这是基本流程。4.5 有限元软件里的应力解读与“奇异点”陷阱用有限元软件做分析时初学者最喜欢做的一件事就是放大应力云图、找到最大应力点然后拿它去跟许用应力比。这个做法有个致命陷阱——应力奇异点。在尖锐缺口、直角拐角、点载荷和约束位置理论上应力趋于无穷大有限元网格只要不断细化应力值就会不断飙升。这个位置的最大应力值是数值假象根本不代表真实物理应力。如果你拿这个数去校核强度永远都是超标的。识别奇异点的经验法则有两条第一看网格加密后该点应力是否持续增长而不收敛第二看最大应力是否出现在几何突变处尖角、刚性约束点。如果是这个“最大应力”应该放弃改用热点应力法或缺口应力法来评估。我个人的处理习惯是先查看远离应力集中区域的名义应力分布确认整体应力水平是否合理再对重点区域如圆角单独细化网格分析。对于焊缝、螺栓孔等位置参考相应设计规范比如欧洲规范里的FAT疲劳等级做疲劳评估而不是直接比有限元峰值应力。另外一个常见的软件操作误区是单位设置。Abaqus的经典默认单位系统是mm-N-s应力单位是MPa但CAD里导出的模型可能是m制的不换算的话几何尺寸直接差1000倍应力结果完全是天文数字。每次导入网格后第一步永远是对比尺寸量纲确认无误后再设置材料和荷载。检查的方法很简单量一下模型里一根已知长度的杆件看看数字是否符合预期单位。结尾回到第001篇的主题。应力与应变这两件事表面上是最基础的教材知识但真正吃透它们的人在结构设计、失效分析、仿真计算中的判断力是完全不同的档次。我见过太多工作三五年的人一遇到断裂事故就只会猜测“是不是材料不过关”却连最基本的断口方向与主应力方向的关系都搞不清楚——这恰恰是当年没把45°最大切应力方向记在心里的结果。最后分享一个我在实际工作中的小习惯做任何强度校核或仿真分析之前先用截面法手算一遍名义应力量级对了再上软件。不是因为手算更准而是这个“粗算→精算”的核对过程能挡住大部分单位错误、方向错误和模型错误。等哪一天你拿到一份有限元报告第一反应不是“这软件真厉害”而是“这个应力分布的规律对不对、这个峰值位置是不是奇异点”那你对基础理论的理解才算真正到位了。