第三章 力学基础：弹性力学基本方程、热应力概念、热应变计算、材料力学属性
各位同学，欢迎来到热力耦合仿真的力学基础课。说实话，很多做热分析的朋友一碰到力学就头疼，觉得公式太多。但我的经验是——你不需要成为力学专家，但必须搞懂几个核心概念。否则仿真算出来结果，你都不知道是对是错。
这一章，我们就来啃下这几块硬骨头：弹性力学基本方程、热应力到底是个啥、热应变怎么算，以及材料力学属性那些事儿。我保证，尽量不堆公式，多讲点实际工程里的体会。
3.1 弹性力学基本方程：结构分析的“宪法”
弹性力学基本方程，说白了就是描述一个物体在受力后怎么变形的“游戏规则”。一共就三组方程，我一个个说。
3.1.1 平衡方程
这个最简单。物体内部任何一个微元体，受到的力必须平衡。用大白话说就是：你推它左边，它右边就得有东西顶着，不然它就飞出去了。在三维问题里，这个方程长这样（别怕，看看就行）：
∂σx/∂x + ∂τxy/∂y + ∂τxz/∂z + fx = 0
∂τyx/∂x + ∂σy/∂y + ∂τyz/∂z + fy = 0
∂τzx/∂x + ∂τzy/∂y + ∂σz/∂z + fz = 0其中σ是正应力，τ是剪应力，f是体积力（比如重力）。
我的经验： 做热力耦合时，温度变化会产生等效的“热载荷”，这个载荷必须满足平衡方程。我曾经见过一个案例，有人把热膨胀系数设错了，结果算出来的应力场根本不满足平衡，一看就是错的。
3.1.2 几何方程
这个描述的是“位移”和“应变”之间的关系。说白了，你挪动了多少距离，材料内部就产生了多大的变形（应变）。
εx = ∂u/∂x
εy = ∂v/∂y
εz = ∂w/∂z
γxy = ∂u/∂y + ∂v/∂x
γyz = ∂v/∂z + ∂w/∂y
γzx = ∂w/∂x + ∂u/∂z这里u、v、w是三个方向的位移，ε是正应变，γ是剪应变。
注意： 这个方程是纯几何关系，跟材料没关系。不管是钢还是塑料，只要位移一样，应变就一样。
3.1.3 物理方程（本构关系）
这才是材料的“性格”体现。对于线弹性材料，就是大家熟悉的胡克定律：
εx = (1/E)[σx - ν(σy + σz)] + αΔT
εy = (1/E)[σy - ν(σx + σz)] + αΔT
εz = (1/E)[σz - ν(σx + σy)] + αΔT
γxy = τxy/G
γyz = τyz/G
γzx = τzx/G看到最后那个 + αΔT 了吗？这就是热力耦合的关键——热应变直接加进来了。
核心要点： 弹性力学三组方程（平衡、几何、物理）加上边界条件，就构成了一个完整的边值问题。有限元软件本质上就是在求解这个方程组。你想想看，理解了这三组方程，你就看懂了仿真软件在后台干了什么。
3.2 热应力概念：温度变化带来的“内伤”
热应力，简单说就是：温度变了，但材料不能自由伸缩，于是内部产生了应力。
为什么会这样？我给你举个例子。一根钢棒，两端固定。你把它加热，它想膨胀，但两端被卡死了，动不了。于是它内部就产生了压应力。反过来，冷却就会产生拉应力。
我在项目中遇到过： 有一次做电子芯片的封装仿真，芯片和基板的热膨胀系数不一样。温度一变化，两者之间就产生了巨大的热应力，直接导致焊点开裂。这就是典型的热应力问题。
避坑指南： 我曾经以为只要温度均匀，就没有热应力。错！只要结构有约束，或者不同材料的热膨胀系数不匹配，即使温度均匀，照样会产生热应力。这一点在做多材料装配体分析时尤其要注意。
3.3 热应变计算：公式很简单，但陷阱不少
热应变的计算公式，其实就一行：
ε_thermal = α × ΔT其中α是线膨胀系数（单位：1/℃），ΔT是温度变化量（单位：℃）。
但这里有几个坑，我一个个说：
- α不是常数： 很多材料的α随温度变化。我习惯在仿真中设置成随温度变化的表格，而不是用一个固定值。否则高温区误差很大。
- 参考温度要选对： ΔT = T_current - T_reference。参考温度通常是“零应力温度”，比如材料的成型温度、装配温度。选错了，整个应力场都偏了。
- 各向异性： 有些材料（比如复合材料、单晶硅）在不同方向上的α不一样。这时候热应变就得用张量形式计算，不能简单用一个标量。
警告： 在热力耦合分析中，热应变是作为“初始应变”或“等效载荷”施加到结构上的。如果你在软件里设置了热膨胀系数，但忘了施加温度场，那热应力就是零。嗯，我刚开始学的时候犯过这个低级错误。
3.4 材料力学属性：弹性模量、泊松比、屈服强度
做热力耦合仿真，材料属性是输入的基础。我挑三个最重要的说。
3.4.1 弹性模量（E）
弹性模量，也叫杨氏模量，衡量材料抵抗弹性变形的能力。E越大，材料越“硬”。
关键点： E是温度的函数。温度升高，E通常会下降。比如钢材在500℃时，E可能只有室温下的70%。如果你用室温的E去算高温工况，结果会偏保守（应力算大了）。
3.4.2 泊松比（ν）
泊松比描述的是：你拉长一个材料，它在横向会收缩多少。大多数金属的ν在0.3左右，橡胶接近0.5（不可压缩）。
我的经验： 泊松比对温度不太敏感，一般取常数就行。但如果你做的是大变形或超弹性材料（比如橡胶密封件），泊松比的影响就很大了。
3.4.3 屈服强度（σy）
屈服强度是材料从弹性进入塑性的门槛。热力耦合分析中，如果热应力超过了屈服强度，材料就会发生塑性变形，甚至破坏。
注意： 屈服强度也随温度变化。高温下，屈服强度会显著降低。比如铝合金在200℃时，屈服强度可能只有室温的一半。所以做高温工况的强度评估时，一定要用高温下的屈服强度。
| 材料属性 | 符号 | 单位 | 温度依赖性 | 工程影响 |
|---|---|---|---|---|
| 弹性模量 | E | GPa | 强（随温度升高而降低） | 影响刚度、热应力大小 |
| 泊松比 | ν | 无量纲 | 弱（通常取常数） | 影响横向变形、体积变化 |
| 屈服强度 | σy | MPa | 强（随温度升高而降低） | 决定是否发生塑性变形 |
| 线膨胀系数 | α | 1/℃ | 中等（部分材料变化明显） | 决定热应变大小 |
总结一下： 做热力耦合仿真，材料属性一定要用随温度变化的曲线，不要偷懒用常数。我个人的习惯是，从材料手册或实验数据中至少取5个温度点的数据，然后插值使用。这样算出来的结果才靠谱。
3.5 本章知识体系图
下面我用一张图把本章的核心逻辑串起来，方便你理解：
这张图把本章的四个核心模块串在了一起。你从上往下看：弹性力学方程是基础框架，热应力和热应变是热力耦合的“桥梁”，材料属性是输入参数。最终落脚点是强度评估——也就是判断结构会不会坏。
好了，力学基础就讲到这里。记住，这些概念是热力耦合仿真的“内功心法”。招式（软件操作）可以学，但内功不扎实，迟早要出问题。