应力与应变入门 — σ=F/A·ε=ΔL/L 到FEM结构分析
1. 应力、应变与FEM — Von Mises应力的真相
在设计结构时，最终要问的问题是"这个部件会不会破裂？"为了获得答案，需要理解应力（Stress）和应变（Strain）这两个概念。FEM分析软件是计算这两个量的工具，而正确解读结果需要对概念的深刻理解。
2. 法向应力与应变的定义
将力 $F$ [N] 除以截面积 $A$ [m²] 就是法向应力（Normal Stress） $\sigma$。
$$\sigma = \frac{F}{A} \quad \text{[Pa = N/m}^2\text{]}$$一般约定拉伸为正（$\sigma > 0$），压缩为负（$\sigma < 0$）。
将零件的伸长量 $\Delta L$ 除以原始长度 $L$ 就是法向应变（Normal Strain） $\varepsilon$。
$$\varepsilon = \frac{\Delta L}{L} \quad \text{[无量纲，或mm/mm]}$$胡克定律（弹性范围）
在弹性范围内，应力和应变成正比。比例常数为杨氏模量（Young's Modulus） $E$。
$$\sigma = E\varepsilon$$| 材料 | 杨氏模量 E [GPa] | 说明 |
|---|---|---|
| 结构用碳钢（SS400） | 200〜210 | 最常用的结构材料 |
| 铝合金（A6061-T6） | 68〜70 | 航空、汽车轻量化材料 |
| 钛合金（Ti-6Al-4V） | 110〜115 | 航空、医疗植入物 |
| 碳纤维增强树脂（CFRP） | 70〜300（纤维方向） | 高强度、各向异性 |
| 混凝土 | 20〜40 | 按压缩强度设计 |
| 橡胶（天然橡胶） | 0.001〜0.01 | 超弹性、大变形 |
3. 屈服应力与安全系数
将材料拉伸时，当应力达到某个水平时会产生永久变形（塑性变形）。这个应力点称为屈服应力（Yield Stress） $\sigma_Y$。继续拉伸，材料会在硬化的同时最终断裂，此时的应力为抗拉强度（Ultimate Tensile Strength） $\sigma_B$（或 $\sigma_{UTS}$）。
设计中要保证许用应力 $\sigma_{allow}$ 达到的安全系数 $S$ 满足要求。
$$S = \frac{\sigma_Y}{\sigma_{allow}} \geq S_{min}$$JIS、ASME等规范要求的最小安全系数（$S_{min}$）因用途而异，但一般机械设计的静荷载情况下通常为1.5〜3.0。
| 材料 | 屈服应力 σY [MPa] | 抗拉强度 σB [MPa] |
|---|---|---|
| SS400（一般结构钢） | 245 | 400〜510 |
| S45C（机械结构用碳钢） | 490 | 690 |
| A6061-T6（铝合金） | 276 | 310 |
| Ti-6Al-4V（钛合金） | 880〜1000 | 950〜1100 |
| CFRP（单向材、纤维方向） | 600〜1500 | 700〜1800 |
| 混凝土（压缩） | —（无屈服） | -20〜-40（压缩强度） |
4. 泊松比与横向应变
拉伸杆时，在轴向（拉伸方向）伸长的同时，垂直方向（横向）会收缩。横向应变与轴向应变之比（绝对值）称为泊松比（Poisson's Ratio） $\nu$。
$$\varepsilon_{lateral} = -\nu \varepsilon_{axial}$$理论上泊松比的范围为 $-1 \leq \nu \leq 0.5$，但一般各向同性材料通常满足 $0 < \nu < 0.5$。
- 金属（钢、铝）：$\nu \approx 0.27〜0.34$
- 橡胶（非压缩性近似）：$\nu \approx 0.499$
- 混凝土：$\nu \approx 0.15〜0.20$
- CFRP（面内方向）：$\nu \approx 0.3$（纤维方向具有各向异性）
体积变化
一般3D应力状态下的体积应变为：
$$\frac{\Delta V}{V} = \varepsilon_x + \varepsilon_y + \varepsilon_z = \frac{1-2\nu}{E}(\sigma_x + \sigma_y + \sigma_z)$$当 $\nu = 0.5$ 时 $\Delta V = 0$（体积不变）——这解释了为什么橡胶和水几乎呈不可压缩状态。在FEM分析中使用接近0.5的泊松比时，刚度矩阵在数值上会变得不稳定（锁定现象）。需要采用特殊的单元定式（B-bar法、混合型等）。
5. 剪切应力与剪切应变
平行于表面的力（剪切力 $V$ [N]）除以截面积 $A_s$ 就是剪切应力（Shear Stress） $\tau$。
$$\tau = \frac{V}{A_s} \quad \text{[Pa]}$$剪切应力对应的变形是剪切应变 $\gamma$（剪切角 [rad]），在弹性范围内满足：
$$\tau = G\gamma$$$G$ 是剪切弹性系数（横弹性系数）。对于各向同性材料，$G$、$E$、$\nu$ 之间的关系为：
$$G = \frac{E}{2(1+\nu)}$$对于钢铁（$E = 206\,\text{GPa}$，$\nu = 0.3$）：$G = 206/(2 \times 1.3) \approx 79\,\text{GPa}$
6. 应力张量与多轴应力状态
三维空间中任意点的应力状态可用3×3对称张量 $\sigma_{ij}$ 表示。
$$\boldsymbol{\sigma} = \begin{pmatrix} \sigma_{xx} & \tau_{xy} & \tau_{xz} \\ \tau_{yx} & \sigma_{yy} & \tau_{yz} \\ \tau_{zx} & \tau_{zy} & \sigma_{zz} \end{pmatrix}$$由于对称性（$\tau_{xy} = \tau_{yx}$ 等），独立分量共6个。FEM计算这6个分量的分布。
主应力与主应力方向
应力张量的特征值为主应力（$\sigma_1 \geq \sigma_2 \geq \sigma_3$），对应的特征向量为主应力方向。在主应力面上不存在剪切应力。
最大剪切应力：$\tau_{max} = (\sigma_1 - \sigma_3)/2$
二维的莫尔应力圆
在平面应力状态下（$\sigma_{zz} = \tau_{xz} = \tau_{yz} = 0$），可以用莫尔应力圆几何方法求主应力和最大剪切应力。圆的半径等于 $\tau_{max}$，圆心在 $(\sigma_x + \sigma_y)/2$。虽然CAE求解器会自动数值计算，但手工理解这个方法能帮助快速判断。
7. Von Mises应力与屈服准则
Von Mises应力（等效应力）由主应力 $\sigma_1, \sigma_2, \sigma_3$ 定义为：
$$\sigma_{vM} = \sqrt{\frac{1}{2}\left[(\sigma_1-\sigma_2)^2 + (\sigma_2-\sigma_3)^2 + (\sigma_3-\sigma_1)^2\right]}$$或用应力分量表示为：
$$\sigma_{vM} = \sqrt{\sigma_{xx}^2 + \sigma_{yy}^2 + \sigma_{zz}^2 - \sigma_{xx}\sigma_{yy} - \sigma_{yy}\sigma_{zz} - \sigma_{zz}\sigma_{xx} + 3(\tau_{xy}^2 + \tau_{yz}^2 + \tau_{zx}^2)}$$屈服条件： 当 $\sigma_{vM} \geq \sigma_Y$ 时开始屈服（塑性变形开始）
物理含义（偏应力能）
Von Mises准则基于这样的假设："形状变化所消耗的应变能超过一定值时就会屈服"。纯静水压（$\sigma_1 = \sigma_2 = \sigma_3$）无论多大也不会引起屈服——这与金属材料的实验观察相符。
与Tresca准则的对比
| 准则 | 屈服条件 | 特点 |
|---|---|---|
| Von Mises | $\sigma_{vM} = \sigma_Y$ | 与实验值吻合好。金属材料标准。 |
| Tresca（最大剪切应力） | $\tau_{max} = \sigma_Y/2$ | 保守（安全）。用于延性材料。 |
| Mohr-Coulomb | （含摩擦角的式子） | 适用于土、岩石等内摩擦重要的材料。 |
8. FEM节点应力的注意事项
FEM计算得到的应力有几个重要注意点。不理解这些会导致分析结果的误读。
积分点到节点的外推
FEM单元内的应变和应力在积分点（Gauss点）被精确计算。之后把积分点的值外推到节点，再对相邻单元的节点值进行平均化（平均处理），最后绘制云图。这是标准步骤。
在网格粗糙的区域（特别是应力集中区），外推误差很大，可能导致实际的最大应力被低估。
通过应力跳跃检查网格质量
相邻单元在公共边界处应力值差很大（跳跃）说明此处网格不足。如果看未平均化（Unaveraged）的应力云图，单元边界处的跳跃就会显现出来。跳跃越小说明解的收敛性越好。
自由表面的应力边界条件
与外界接触的自由表面（比如暴露在空气中的表面），其法向应力应为零（若无外压）。FEM结果中如果在自由表面看到很大的法向应力，那说明边界条件设置有问题或网格有问题。
9. 拉伸试验的理解
拉伸试验（JIS Z 2241）是最直接地测量材料机械性质的基本试验。用恒定速度拉伸试样，记录荷载与位移。
工程应力-应变曲线的读法
- 比例限（Proportional Limit）：应力-应变呈线性（胡克定律成立）的最大应力。
- 弹性限：卸载后无永久变形的最大应力。
- 上屈服点、下屈服点：在软钢中明显。转位开始活动产生应力尖峰（上），之后平稳的屈服应力（下）。
- 缩颈开始点：最大荷载点（对应 $\sigma_B$）。从这里开始试样局部变细（缩颈）。
- 断裂点：试样断成两截时的应力和应变。
工程应力 vs 真应力、真应变
工程应力用原始截面积计算（$\sigma = F/A_0$），但缩颈后截面积大幅减小，与实际应力偏离。
$$\sigma_{true} = \sigma_{eng}(1 + \varepsilon_{eng}), \quad \varepsilon_{true} = \ln(1 + \varepsilon_{eng})$$在有限变形FEM分析（大变形）中要用真应力和真应变。向Abaqus等软件输入材料数据时，需要提供真应力-塑性应变的曲线。
10. 实践：汽车悬架臂的强度评估
汽车悬架臂在行驶中要承受路面传来的上下、前后、左右等多个方向的复合荷载。以最严苛荷载（最大颠簸：垂直加速度3G）为条件进行FEM强度评估，这是典型的工程流程。
分析流程
- 荷载设置：轮胎接地点施加垂直9G、横向1.5G、前后1G的复合荷载。
- 边界条件：车身连接点（衬套）全部约束。
- 材料：铝合金A6061-T6（$E = 69\,\text{GPa}$，$\sigma_Y = 276\,\text{MPa}$）
- 网格：在应力集中部位（支架、圆角）采用细密网格（单元大小1mm）。
- 结果检查：Von Mises最大应力目标为屈服应力的70%（安全系数1.43）以下。
轻量化与刚度的权衡
臂厚度降低可以减轻重量，但会增加应力集中。当代的做法是使用拓扑优化（Abaqus/Tosca或Altair OptiStruct）："在保证应力约束的前提下去除荷载路径中不必要的材料"。量产车通过CAE驱动的设计，从初期样品起就能达到95%以上的强度目标。