牛顿运动三定律 — 从F=ma到有限元运动方程[M]ü+[C]u̇+[K]u=F
目录
1. 为什么从牛顿定律学CAE？
教授，我想学FEM有限元分析，但感觉不知道从哪里入手，教材一翻就是矩阵和偏微分方程，完全看不懂……
别慌，其实FEM的核心就是牛顿第二定律。你把一辆车的车身拆成成千上万个小单元，对每个节点写一遍F=ma，然后组装成矩阵——这就是FEM的运动方程。高中物理学扎实了，FEM的物理本质其实很好理解。
有限元分析（FEM）表面上看起来是复杂的线性代数，但它的物理根基非常简单：对每一个离散化的质点，牛顿第二定律必须成立。理解了牛顿三定律，你就掌握了CAE仿真的物理语言。
现代工程仿真中，无论是汽车碰撞安全性评估（LS-DYNA）、飞机机翼颤振分析（MSC Nastran），还是土木工程的地震响应计算，其控制方程的物理根源都是牛顿力学。这一章我们从最基础的概念出发，逐步推导到FEM的运动方程。
2. 第一定律：惯性定律
牛顿第一定律（惯性定律）表述为：
这看似简单，但在工程中有深刻含义：
- 静止结构：桥梁、建筑物处于合力为零的平衡状态
- 惯性载荷：飞机急转弯时乘客"被甩向外侧"，是因为身体的惯性倾向于保持直线运动
- FEM静力分析：线性静力分析的基础假设就是加速度为零，即 $[K]\{u\} = \{F\}$
那为什么卫星在太空中能一直绕地球转，不会飞走也不会坠落？这和第一定律有关系吗？
完全正确！卫星受到重力（向心力），这个力不断改变速度方向，使卫星做圆周运动而不是直线飞出去。如果突然重力消失，按第一定律，卫星会沿切线方向飞走。在卫星结构的CAE分析里，这个离心"惯性力"必须作为载荷输入——轨道机动时的结构强度设计就是这么算的。
在CAE中，惯性力（pseudo-static）的处理是结构分析的基本操作。对于旋转机械，惯性载荷通过旋转角速度 $\omega$ 施加：
3. 第二定律：F=ma的深层含义
牛顿第二定律是整个经典力学和FEM动力学的核心：
这个公式包含了三个关键物理量之间的关系：
| 物理量 | 符号 | 单位（SI） | 在FEM中的对应 |
|---|---|---|---|
| 力 | $F$ | N (kg·m/s²) | 载荷向量 $\{F\}$ |
| 质量 | $m$ | kg | 质量矩阵 $[M]$ |
| 加速度 | $a$ | m/s² | 节点加速度 $\{\ddot{u}\}$ |
需要注意的是，在FEM中，"质量"不是集中在一个点上，而是分布在整个单元体积内。这就需要引入质量矩阵的概念：
其中 $\rho$ 是材料密度，$[N]$ 是形函数矩阵，描述单元内部位移如何由节点位移插值得到。
4. 第三定律：作用力与反作用力
牛顿第三定律："每一个作用力都有一个大小相等、方向相反的反作用力。"
在CAE中，第三定律体现在：
- 接触分析：两个零件接触时，接触力满足第三定律。汽车碰撞中，A车对B车施加的力，等于B车对A车的反作用力
- 单元间力传递：FEM中相邻单元共享节点，节点处力的传递自动满足第三定律
- 约束反力：地基对建筑物的支持力，等于建筑物对地基的压力
等等，既然作用力等于反作用力，那撞击时两辆车受到的力相同，为什么小车比大车更容易被撞坏？
问得很好！力确实相同，但加速度不同：a = F/m。小车质量小，同样的力产生更大的减速度（冲击更猛），因此变形和损伤更严重。这也是为什么NCAP碰撞安全测试里，小型车需要特别设计吸能结构——用更长的变形行程来降低乘员舱的峰值减速度。
5. 力的分解与叠加
在实际工程问题中，载荷往往是多方向、多类型的。力的叠加原理（线性系统中）允许将复杂载荷分解为简单载荷的叠加：
以飞机机翼为例，作用在翼面上的载荷包括：
- 气动升力（法向）：$L = \frac{1}{2}\rho v^2 S C_L$
- 气动阻力（切向）：$D = \frac{1}{2}\rho v^2 S C_D$
- 结构自重（重力）
- 发动机推力（轴向）
- 热载荷（高温区域的温度应力）
在FEM分析中，这些载荷可以分别施加，最终的响应是各载荷响应的叠加（仅限线性分析）。
| 载荷类型 | FEM施加方式 | 典型量级 |
|---|---|---|
| 集中力 | 节点力 | 发动机推力：10-500 kN |
| 分布压力 | 面载荷 | 气动压力：0.01-0.1 MPa |
| 体力（重力） | 体积力 | 自重：g = 9.81 m/s² |
| 热载荷 | 温度场输入 | 发动机叶片：600-1200°C |
| 惯性载荷 | 加速度场 | 起降：2-3g |
6. 从F=ma到FEM运动方程
这是本章最关键的推导。将牛顿第二定律应用于每个有限单元，考虑弹性恢复力、阻尼力和外载荷，得到有限元动力学方程：
各项物理含义如下：
| 项 | 矩阵/向量 | 物理含义 | 类比 |
|---|---|---|---|
| 惯性力项 | $[M]\{\ddot{u}\}$ | 质量×加速度 | F=ma中的"ma" |
| 阻尼力项 | $[C]\{\dot{u}\}$ | 阻尼×速度 | 空气阻力、结构内阻尼 |
| 弹性恢复力项 | $[K]\{u\}$ | 刚度×位移 | 弹簧力F=kx |
| 外载荷 | $\{F(t)\}$ | 时变外力向量 | F=ma中的"F" |
各矩阵的构成：
其中 $[B]$ 是应变-位移矩阵，$[D]$ 是材料本构矩阵（弹性张量）。
这个运动方程很大，一般怎么求解呢？是直接联立方程组吗？
分情况。如果只关心结构在静载荷下的变形，加速度和速度都是零，方程简化为 [K]{u}={F}，直接求解线性方程组就好。但如果是碰撞分析、地震响应这类动态问题，就要做时间积分——把时间轴切成很多小步，一步一步往前推进。常用的有Newmark法、中心差分法等，LS-DYNA就主要用显式中心差分。
7. 工程实例：汽车碰撞分析
汽车正面碰撞是牛顿定律在CAE中最典型的应用场景之一。以下是一个典型的碰撞工况参数：
| 参数 | 典型值 | 说明 |
|---|---|---|
| 碰撞初速度 | 56 km/h (15.6 m/s) | NCAP正面碰撞标准工况 |
| 车身质量 | 1500 kg | 含乘员 |
| 碰撞持续时间 | ~100 ms | 从接触到最大变形 |
| 峰值减速度 | 30-50 g | 取决于吸能结构设计 |
| 车头压溃量 | 400-600 mm | 前舱吸能区变形量 |
对乘员舱施加第二定律：
这就是为什么安全带预紧器和安全气囊需要在碰撞开始后的极短时间内（约20-30 ms）就触发——必须在乘员相对于车身继续前冲之前完成保护动作。
在LS-DYNA中，这类显式动力学分析的时间步长由CFL（Courant-Friedrichs-Lewy）条件限制：
其中 $L_e$ 是最小单元尺寸，$c_s$ 是材料中的弹性波速。对于钢（$E = 210$ GPa，$\rho = 7850$ kg/m³），$c_s \approx 5170$ m/s，若最小单元为5 mm，则 $\Delta t \approx 10^{-6}$ s。100 ms的碰撞需要约10万个时间步。
8. 航空航天中的动力学
飞机在飞行中所受的载荷是多轴、时变的，必须用FEM运动方程全面评估。以民用客机机翼颤振分析为例：
颤振（Flutter）是气动弹性耦合不稳定现象。当气动力频率与结构固有频率耦合时，振幅会指数增长，可在数秒内导致结构破坏。控制方程：
其中下标 $s$ 表示结构项，$a$ 表示气动项（与频率相关）。颤振速度 $V_F$ 是使系统失稳的临界速度，FAA要求颤振速度须高于设计限速的1.15倍（即15%颤振裕度）。
| 飞机类型 | 颤振分析软件 | 典型自由度数 | 关键模态 |
|---|---|---|---|
| 民用客机（B737级） | MSC Nastran / ZAERO | 10⁶ ~ 10⁷ | 弯扭耦合 |
| 战斗机 | STARS / NASTRAN | 10⁵ ~ 10⁶ | 机翼/外挂颤振 |
| 直升机旋翼 | CAMRAD II | 10⁴ ~ 10⁵ | 地面共振/空中共振 |
9. CAE实践要点
实际做FEM分析的时候，牛顿定律在哪些地方会"出问题"或者需要特别注意？
有几个坑要注意。第一，牛顿力学只适用于低速（远低于光速）和宏观物体，量子效应和相对论效应不在其中。第二，FEM运动方程里的质量矩阵有"一致质量矩阵"和"集中质量矩阵"两种，后者精度低但计算快，选哪个会影响动力分析结果。第三，阻尼矩阵[C]最难准确建模，工程上常用瑞利阻尼C=αM+βK，但α、β要靠实验数据标定。
在实际CAE工作中，基于牛顿定律的FEM分析要注意以下几点：
- 单位一致性：最常见错误之一。推荐使用SI单位或mm-N-s单位制（后者在汽车行业广泛使用），混用单位会导致质量矩阵数量级错误
- 边界条件的物理含义：约束点实际上是在施加反作用力。过约束（hyperstatic）会导致人为应力集中
- 质量矩阵类型：一致质量矩阵（consistent）精度高，集中质量矩阵（lumped）计算效率高。显式动力学（碰撞）通常用集中质量矩阵
- 阻尼模型：瑞利（Rayleigh）阻尼是最常用的近似。$\alpha$ 影响低频段，$\beta$ 影响高频段
- 非线性检验：大变形（位移/特征尺寸 > 10%）或材料进入塑性时，线性假设失效，需要启用几何或材料非线性
其中阻尼系数通过以下关系由两个已知模态的阻尼比 $\xi_i, \xi_j$ 确定：
总结
牛顿三定律是整个CAE力学分析的物理基础：
- 第一定律 → 静力平衡，FEM静力分析 $[K]\{u\}=\{F\}$
- 第二定律 → 动力学，FEM运动方程 $[M]\{\ddot{u}\}+[C]\{\dot{u}\}+[K]\{u\}=\{F(t)\}$
- 第三定律 → 接触力、单元间力传递、约束反力
掌握这些基础后，FEM的数学框架将变得有血有肉，而不再是抽象的矩阵符号游戏。