连续介质力学是物理学与工程学的核心分支，研究物质在宏观尺度下的运动与变形行为。它将物质视为连续填充空间的介质，通过场方程描述其力学行为，广泛应用于流体力学、固体力学、地球物理学和生物力学等领域。
连续介质力学（Continuum Mechanics）是研究物质在连续介质假设下的力学行为的学科。它不关注物质的微观粒子结构（如分子、原子），而是将物质理想化为连续填充整个空间的介质，用连续函数描述其物理量（密度、速度、应力、温度等）的分布与变化。
这一学科的核心思想是：当观察尺度远大于分子平均自由程时，物质的离散分子结构可以被忽略，物质可视为连续的整体。连续介质力学因此成为连接微观物理与宏观工程的桥梁。
| 时期 |
代表人物 |
贡献 |
| 17世纪 |
牛顿 (Newton) |
提出粘性流体模型（牛顿流体） |
| 18世纪 |
欧拉 (Euler) |
建立理想流体运动方程 |
| 19世纪 |
纳维 (Navier)、斯托克斯 (Stokes) |
推导粘性流体运动方程（Navier-Stokes方程） |
| 19世纪 |
柯西 (Cauchy) |
建立应力与应变理论，提出柯西应力张量 |
| 19世纪 |
圣维南 (Saint-Venant) |
发展弹性力学理论 |
| 20世纪 |
普朗特 (Prandtl) |
提出边界层理论，连接粘性流体与理想流体 |
| 20世纪 |
里夫林 (Rivlin)、特鲁斯德尔 (Truesdell) |
建立理性力学框架 |
连续介质假设是连续介质力学的基石。该假设认为物质在空间中是连续分布的，每个空间点对应一个物质点，物质点的物理性质（密度、速度、应力等）是空间位置和时间的连续函数。
假设的有效性条件：当特征尺度 L 远大于分子平均自由程 λ 时成立，即：
Kn=Lλ≪1
其中 Kn 为克努森数（Knudsen number）。
具体例子：在标准大气压下，空气的分子平均自由程 λ≈68 nm。对于一个直径 10 cm 的管道流动：
Kn=0.168×10−9=6.8×10−7≪1
连续介质假设完全适用。但在地球大气层边缘（海拔 100 km 以上），分子平均自由程可达数米，此时连续介质假设失效，需要改用分子动力学方法（如玻尔兹曼方程）。
| 场景 |
特征尺度 L |
平均自由程 λ |
Kn |
连续介质假设？ |
| 水杯中水流 |
0.1 m |
∼0.1 nm |
10−9 |
✅ 完全成立 |
| 管道空气流动 |
0.1 m |
68 nm |
6.8×10−7 |
✅ 成立 |
| 微流控芯片 |
10 µm |
68 nm |
6.8×10−3 |
⚠️ 边界滑移需修正 |
| 高真空系统 |
1 cm |
1 m |
100 |
❌ 失效 |
┌──────────────────────────────────────────────────┐
│ 连续介质力学 │
├───────────────────────┬──────────────────────────┤
│ 运动学 (Kinematics) │ 动力学 (Dynamics) / 平衡律 │
│ • 变形梯度 │ • 质量守恒 │
│ • 应变张量 │ • 动量守恒 (柯西运动方程) │
│ • 速度梯度 │ • 能量守恒 │
│ • 旋度与散度 │ • 熵不等式 (热力学第二定律) │
├───────────────────────┴──────────────────────────┤
│ 本构关系 (Constitutive Relations) │
│ • 线弹性 (胡克定律) • 粘性 (牛顿粘性定律) │
│ • 塑性 (屈服准则) • 粘弹性 │
├──────────────────────────────────────────────────┤
│ 具体分支与应用领域 │
│ ┌─────────┐ ┌──────────┐ ┌───────────────────┐ │
│ │ 固体力学 │ │ 流体力学 │ │ 交叉学科 │ │
│ │ 弹性力学 │ │ 理想流体 │ │ 生物力学 │ │
│ │ 塑性力学 │ │ 粘性流体 │ │ 地球物理学 │ │
│ │ 断裂力学 │ │ 可压缩流 │ │ 材料科学 │ │
│ └─────────┘ └──────────┘ └───────────────────┘ │
└──────────────────────────────────────────────────┘
连续介质力学的数学语言是张量分析。由于物理量（应力、应变、应变率等）不依赖于坐标系的选择，张量提供了描述这些物理量的自然框架。
| 阶数 |
名称 |
例子 |
分量数量 (3D) |
| 0 |
标量 (Scalar) |
温度 T、密度 ρ、压力 p |
1 |
| 1 |
向量 (Vector) |
位移 u、速度 v、力 f |
3 |
| 2 |
二阶张量 (2nd-order tensor) |
应力 σ、应变 ε |
9 |
| 4 |
四阶张量 (4th-order tensor) |
弹性刚度张量 C |
81 |
在笛卡尔坐标系中，向量 v 的分量记为 vi（i=1,2,3）。二阶张量 T 的分量记为 Tij。
求和约定：重复指标表示求和：
a⋅b=aibi=i=1∑3aibi
张量运算：
- 点积：(A⋅B)ij=AikBkj
- 双点积：A:B=AijBij
- 散度：(∇⋅σ)i=∂xj∂σij
张量的本质特征是其在坐标变换下的行为。当从坐标系 {xi} 变换到 {xi′} 时，变换矩阵为 Qij=∂xj∂xi′：
- 一阶张量：vi′=Qijvj
- 二阶张量：Tij′=QikQjlTkl
具体数值例子：设二维应力张量：
σ=[10050500] MPa
绕 z 轴旋转 30∘ 后，Q=[cos30∘−sin30∘sin30∘cos30∘]=[0.866−0.50.50.866]
σ11′=Q1kQ1lσkl
=0.8662×100+0.866×0.5×50+0.5×0.866×50+0.52×0
=75.0+21.65+21.65+0=118.3 MPa
旋转后应力分布发生变化，这说明应力是张量而非标量——它的值依赖于观察方向。
连续介质力学有两种描述运动的方式：
| 描述方式 |
拉格朗日描述 (Lagrangian) |
欧拉描述 (Eulerian) |
| 关注点 |
追踪物质点的运动 |
观察空间点的物理量变化 |
| 自变量 |
初始坐标 X，时间 t |
空间坐标 x，时间 t |
| 典型应用 |
固体力学 |
流体力学 |
| 加速度 |
a=∂t2∂2u(X,t) |
a=∂t∂v+v⋅∇v |
物质导数（Material Derivative）是连接两种描述的关键概念：
DtD(⋅)=∂t∂(⋅)+v⋅∇(⋅)
具体数值例子：考虑管道中的水流，速度场 v=(vx,0,0)，其中 vx=2+0.1x (m/s)。在某点 x=5 m 处，温度场 T(x,t)=20+2t (°C)，则温度的物质导数为：
DtDT=∂t∂T+vx∂x∂T=2+(2+0.5)×0=2 °C/s
即随流体运动的观察者感受到 2 °C/s 的温升。如果温度场同时有空间梯度，比如 T(x,t)=20+0.5x，则：
DtDT=0+(2+0.1×5)×0.5=2.5×0.5=1.25 °C/s
变形梯度张量 F 描述了物质点从初始构型到当前构型的映射：
Fij=∂Xj∂xi
极分解：F=RU=VR，其中 R 是旋转张量，U 和 V 分别是右和左伸长张量。
| 应变张量 |
定义 |
适用场景 |
| 无穷小应变 |
εij=21(∂xj∂ui+∂xi∂uj) |
小变形（≪1%应变） |
| Green-Lagrange应变 |
Eij=21(FkiFkj−δij) |
有限变形（橡胶、生物组织） |
| 变形率张量 |
Dij=21(∂xj∂vi+∂xi∂vj) |
流体力学 |
具体数值例子——立方体单向拉伸：一个边长为 1 m 的立方体，沿 x 方向拉伸至 1.1 m，y 和 z 方向收缩至 0.95 m。
变形梯度：
F=1.10000.950000.95
无穷小应变（适用于小变形）：
ε=0.1000−0.05000−0.05
Green-Lagrange应变（适用于大变形）：
E11=21(F112−1)=21(1.21−1)=0.105
E22=21(F222−1)=21(0.9025−1)=−0.04875
两者的差异在应变较大时开始显著。
速度梯度张量 L=∇v，可分解为对称部分和反对称部分：
L=D+W
- 变形率张量 D=21(L+LT)：描述形状变化速率
- 旋度张量 W=21(L−LT)：描述刚性旋转
质量既不能创生也不能消灭。在欧拉描述下：
∂t∂ρ+∇⋅(ρv)=0
或等价地用物质导数表示：
DtDρ+ρ∇⋅v=0
对于不可压缩流体（密度恒定），简化为：
∇⋅v=0
具体数值例子：考虑一个圆管收缩段，入口直径 D1=0.2 m，出口直径 D2=0.1 m，入口速度 v1=1 m/s。假设不可压缩流动，由质量守恒：
A1v1=A2v2
π(0.1)2×1=π(0.05)2×v2
v2=0.0025π0.01π=4 m/s
出口速度是入口的4倍，这与日常经验（捏住水管口水流加速）一致。
牛顿第二定律应用于连续介质：
ρDtDv=∇⋅σ+ρb
其中 σ 为柯西应力张量，b 为单位质量体力（如重力）。
在笛卡尔坐标系中展开（i方向）：
ρDtDvi=∂xj∂σij+ρbi
角动量守恒要求应力张量是对称的：σij=σji。这意味着9个应力分量中只有6个是独立的。物理意义上，这排除了连续介质中存在"体力矩"的可能性。
ρDtDe=σ:D+ρr−∇⋅q
其中 e 是单位质量内能，r 是热源项，q 是热通量向量，σ:D=σijDij 是应力功率（单位体积的功率输入）。
ρDtDs+∇⋅(Tq)−Tρr≥0
其中 s 是单位质量熵，T 是温度。这约束了本构关系的可能形式。
本构关系（Constitutive Relations）是描述材料响应（应力或热流）与变形/温度之间关系的数学方程。不同材料具有截然不同的本构行为。
各向同性线弹性材料的本构关系：
σ=λ(trε)I+2με
或用分量形式：
σij=λεkkδij+2μεij
其中 λ 和 μ（也称为 G）是拉梅常数（Lamé constants）。
完整的弹性常数对照表：
| 常数名称 |
符号 |
与 E, ν 的关系 |
典型值（钢材） |
| 杨氏模量 |
E |
E |
200 GPa |
| 泊松比 |
ν |
ν |
0.3 |
| 剪切模量 |
μ (或 G) |
2(1+ν)E |
76.9 GPa |
| 体积模量 |
K |
3(1−2ν)E |
166.7 GPa |
| 拉梅第一常数 |
λ |
(1+ν)(1−2ν)Eν |
115.4 GPa |
具体数值例子——钢棒拉伸：一根钢棒（E=200 GPa，ν=0.3）受轴向应力 σ11=100 MPa。
轴向应变为：
ε11=Eσ11=200×109100×106=0.0005=0.05%
横向应变为：
ε22=ε33=−νε11=−0.3×0.0005=−0.00015=−0.015%
一根 1 m 长的钢棒在 100 MPa 应力下伸长 0.5 mm，直径减少 0.015%。
对于牛顿流体（如空气、水、甘油），应力与变形率张量成正比：
σ=−pI+2μD+λv(∇⋅v)I
其中 p 是静水压力，μ 是动力粘度，λv 是体积粘度。
几种常见流体的粘度对比：
| 流体 |
动力粘度 μ (Pa·s) |
运动粘度 ν (m²/s) |
| 空气 (20°C) |
1.8×10−5 |
1.5×10−5 |
| 水 (20°C) |
1.0×10−3 |
1.0×10−6 |
| 发动机油 (SAE 30, 20°C) |
0.29 |
3.2×10−4 |
| 甘油 (20°C) |
1.41 |
1.12×10−3 |
| 蜂蜜 (20°C) |
∼10 |
∼7×10−3 |
牛顿内摩擦定律的数值例子：两块平行平板，间距 h=1 mm，下板固定，上板以 U=0.5 m/s 运动。板间充满水（μ=10−3 Pa⋅s）。速度梯度近似为：
dydv≈hU=0.0010.5=500 s−1
剪切应力：
τ=μdydv=10−3×500=0.5 Pa
如果换用机油（μ=0.29 Pa⋅s），剪切应力变为：
τ=0.29×500=145 Pa
是水中的290倍——这就是为什么机油比水"粘稠"得多。
当应力超过弹性极限时，材料发生不可逆的塑性变形。常用的屈服准则：
冯·米塞斯屈服准则 (von Mises)：
23s:s=σy
其中 s=σ−31(trσ)I 是偏应力张量，σy 是屈服应力。
具体数值例子：一根纯铝棒受到轴向拉伸应力 σ11=150 MPa，且 σ22=σ33=0。铝的屈服应力 σy=100 MPa。
偏应力张量为：
s11=150−3150=150−50=100 MPa
s22=s33=0−50=−50 MPa
冯·米塞斯等效应力：
σvm=23(s112+s222+s332)=23(1002+502+502)=23×15000=22500=150 MPa
由于 σvm=150 MPa>σy=100 MPa，铝棒将发生不可逆的塑性变形（屈服）。
| 材料类型 |
本构关系 |
特征参数 |
典型现象 |
| 线弹性固体 |
σ=C:ε |
E, ν, λ, μ |
弹簧变形（可恢复） |
| 牛顿流体 |
σ=−pI+2μD |
μ, λv |
水的流动（剪切阻抗） |
| 弹塑性固体 |
屈服准则+流动法则 |
σy, Et |
金属塑性成形（不可恢复） |
| 粘弹性材料 |
/ |
松弛模量 G(t) |
聚合物蠕变（时间相关） |
| 非牛顿流体 |
τ=μ(γ˙)γ˙ |
幂律指数 n |
牙膏、血液（剪切变稀） |
将牛顿流体的本构关系代入柯西运动方程，得到著名的Navier-Stokes方程：
ρ(∂t∂v+v⋅∇v)=−∇p+μ∇2v+(λv+3μ)∇(∇⋅v)+ρb
对于不可压缩牛顿流体（∇⋅v=0），方程简化为：
ρ(∂t∂v+v⋅∇v)=−∇p+μ∇2v+ρb
各项的物理意义：
局部加速度ρ∂t∂v+对流加速度ρv⋅∇v=压力梯度−∇p+粘性耗散μ∇2v+体力ρb
对Navier-Stokes方程进行无量纲化，引入特征速度 U、特征长度 L：
∂t∗∂v∗+v∗⋅∇∗v∗=−∇∗p∗+Re1∇∗2v∗
其中雷诺数 (Reynolds number) 是流体力学中最重要的无量纲参数：
Re=μρUL=νUL
雷诺数表示惯性力与粘性力的比值。
不同流动场景的雷诺数：
| 流动场景 |
特征速度 U |
特征长度 L |
运动粘度 ν |
Re |
流动特征 |
| 细菌游动 |
30 µm/s |
1 µm |
10−6 |
3×10−5 |
层流，粘性主导 |
| 毛细血管血流 |
1 mm/s |
10 µm |
3.5×10−6 |
3×10−3 |
层流 |
| 汽车行驶 (50 km/h) |
14 m/s |
2 m (车长) |
1.5×10−5 |
1.9×106 |
湍流 |
| 飞机巡航 (900 km/h) |
250 m/s |
50 m (翼展) |
3.5×10−5 |
3.6×108 |
湍流，边界层效应 |
| 高尔夫球飞行 |
80 m/s |
0.043 m |
1.5×10−5 |
2.3×105 |
临界区（表面凹坑诱导转捩） |
圆管层流（泊肃叶流动）：当 Re < 2000 时，圆管流动保持层流。速度分布为抛物线形：
vz(r)=4μR2dzdp(1−R2r2)
其中 R 是管径，dzdp 是压力梯度。
数值计算示例：一根直径 D=2 mm 的毛细管，长 L=10 cm，两端压差 Δp=500 Pa，流过水（μ=10−3 Pa⋅s）。
压力梯度：dzdp=−0.1500=−5000 Pa/m
中心线速度（r=0）：
vmax=4×10−3(0.001)2×5000=4×10−310−6×5000=2.5×10−4×5000=1.25 m/s
流量：
Q=8μπR4dzdp=8×10−3π×(0.001)4×5000=8×10−3π×10−12×5000=1.96×10−6 m3/s≈2 mL/s
将线弹性本构关系（胡克定律）代入平衡方程，得到以位移表示的纳维方程（Navier's equation）：
(λ+μ)∇(∇⋅u)+μ∇2u+ρb=ρ∂t2∂2u
静力平衡的特殊情况（忽略惯性项）：
(λ+μ)∇(∇⋅u)+μ∇2u+ρb=0
从纳维方程可以推导出弹性体中波的传播。将位移场分解为 u=∇ϕ+∇×ψ（标量势 + 向量势）：
- P波（纵波）：∂t2∂2ϕ=cp2∇2ϕ，波速 cp=ρλ+2μ
- S波（横波）：∂t2∂2ψ=cs2∇2ψ，波速 cs=ρμ
具体数值例子——钢材中的波速：钢材 ρ=7800 kg/m3，E=200 GPa，ν=0.3。
μ=2(1+ν)E=2×1.3200=76.9 GPa
λ=(1+ν)(1−2ν)Eν=1.3×0.4200×0.3=115.4 GPa
P波波速：
cp=7800115.4×109+2×76.9×109=7800269.2×109=34.5×106≈5870 m/s
S波波速：
cs=780076.9×109=9.86×106≈3140 m/s
这解释了为什么地震时P波总是先于S波到达——P波在钢（和地壳岩石）中的传播速度约为S波的1.87倍。
常见材料的弹性波速：
| 材料 |
ρ (kg/m³) |
E (GPa) |
ν |
cp (m/s) |
cs (m/s) |
| 钢 |
7800 |
200 |
0.3 |
5870 |
3140 |
| 铝 |
2700 |
70 |
0.33 |
5320 |
3120 |
| 混凝土 |
2400 |
30 |
0.2 |
3800 |
2280 |
| 橡胶 |
1100 |
0.01 |
0.49 |
~140 |
~2 |
| 水 |
1000 |
— |
— |
1480 (压缩波) |
0 (无剪切) |
1904年，普朗特（Ludwig Prandtl）在海德堡国际数学大会上提出了边界层概念，这是流体力学史上最重要的突破之一。他认识到：对于高雷诺数流动（Re≫1），粘性效应只集中在固体壁面附近的一个薄层内——这就是边界层。
边界层厚度 δ 随位置 x 和雷诺数变化：
层流边界层（Blasius解）：
δ(x)≈Rex5.0x
湍流边界层：
δ(x)≈Rex1/50.37x
其中 Rex=νUx 是基于距前缘距离 x 的雷诺数。
数值例子：空气（ν=1.5×10−5 m2/s）以 U=10 m/s 流过一块平板。
在 x=0.1 m 处：Rex=1.5×10−510×0.1=6.7×104
层流边界层厚度：
δ(0.1 m)=6.7×1045.0×0.1=2590.5≈0.0019 m=1.9 mm
在 x=1.0 m 处：Rex=1.5×10−510×1.0=6.7×105
此时 Rex>5×105，边界层已转捩为湍流：
湍流边界层厚度：
δ(1.0 m)=(6.7×105)1/50.37×1.0=14.60.37≈0.025 m=25 mm
如果仍用层流公式计算：δ=6.7×1055.0×1.0=8195.0=0.0061 m=6.1 mm
湍流边界层厚度约为层流的4倍——湍流混合更剧烈，将动量更深入地输运到主流中。
当逆压梯度存在时，边界层可能从壁面分离，形成尾流区。分离条件是壁面处速度梯度为零：
∂y∂vxy=0=0
分离导致压差阻力急剧增加，这在钝体绕流中特别重要——高尔夫球的凹坑正是为了延迟分离、减小尾流区。
连续介质力学的有限元方法（FEM）是工程结构分析的基石。从桥梁、建筑到飞机机翼、汽车车身，结构工程师使用连续介质力学方程计算应力分布和变形。
桥梁荷载分析：一座简支钢梁桥，跨度 L=30 m，承受均布荷载 w=50 kN/m（包括自重和活载）。I型截面钢梁的惯性矩 I=0.02 m4。
最大弯矩（跨中）：
Mmax=8wL2=850×103×302=5.625×106 N⋅m
跨中截面顶部最大正应力：
σmax=IMmaxy=0.025.625×106×0.5=140.6 MPa
对于Q345钢（屈服强度345 MPa），安全系数约为2.45，设计合理。
地震波传播（P波和S波）的分析完全基于连续介质力学。地壳结构探测利用地震波的反射和折射数据反演地球内部材料性质。
地壳分层速度模型（典型值）：
| 层位 |
深度 (km) |
cp (km/s) |
cs (km/s) |
密度 (g/cm³) |
| 沉积层 |
0-2 |
2.0-4.0 |
1.0-2.3 |
2.1-2.5 |
| 上地壳花岗岩 |
2-15 |
5.5-6.0 |
3.2-3.5 |
2.6-2.7 |
| 下地壳玄武岩 |
15-35 |
6.5-7.0 |
3.6-3.9 |
2.8-3.0 |
| 地幔顶部 |
35-70 |
8.0-8.2 |
4.5-4.7 |
3.3-3.4 |
地震定位（震源确定）就是通过各台站P波和S波到达时间差，反演波速结构中的震源位置。
人体组织（血管壁、皮肤、肌肉、软骨）的力学行为需要复杂的本构模型来描述：
- 血管壁：非线性粘弹性材料，在脉压作用下周期性扩张
- 软骨：双相介质（固体基质+孔隙水），应力由固液两相共同承担
- 红细胞：在毛细血管中变形通过，需要超弹性膜模型
动脉血压与血管膨胀：主动脉直径 D=2.5 cm，壁厚 h=2 mm，弹性模量 E=0.5 MPa。收缩压（120 mmHg = 16 kPa）和舒张压（80 mmHg = 10.7 kPa）之间的压力变化 Δp=5.3 kPa。
环向应力变化（薄壁简近似）：
Δσh=2hΔp⋅D=2×0.0025330×0.025=33.3 kPa
环向应变变化：
Δεh=EΔσh=5×10533300=0.0666=6.66%
血管直径变化：
ΔD=D⋅Δεh=0.025×0.0666=1.67 mm
即每个心跳周期，主动脉直径从2.5 cm膨胀到约2.67 cm——这就是我们感受到的脉搏！
大气和海洋的运动由连续介质力学方程描述，在大尺度上需要考虑科里奥利力（地球自转效应）。
地转平衡：在大尺度大气运动中，压力梯度力与科里奥利力平衡：
−ρ1∂x∂p+fv=0
−ρ1∂y∂p−fu=0
其中 f=2Ωsinϕ 是科里奥利参数，Ω 是地球自转角速度（7.27×10−5 rad/s），ϕ 是纬度。
在北纬 30∘ 处：
f=2×7.27×10−5×sin30∘=7.27×10−5 s−1
连续介质力学方程（偏微分方程组）在大多数实际问题中无法解析求解，因此数值方法至关重要。
| 方法 |
原理 |
适用场景 |
优点 |
缺点 |
| 有限差分法 (FDM) |
用差分代替微分 |
简单几何形状 |
实现简单，理论成熟 |
复杂边界处理困难 |
| 有限元法 (FEM) |
变分原理+单元插值 |
固体力学，复杂结构 |
几何适应性强 |
对流占优问题不适用 |
| 有限体积法 (FVM) |
控制体积分 |
流体力学 |
天然守恒 |
高精度格式复杂 |
| 边界元法 (BEM) |
将偏微分方程转化为边界积分方程 |
无限域、断裂力学 |
降维，计算效率高 |
仅适用线性问题 |
| 光滑粒子法 (SPH) |
粒子无网格离散 |
大变形、自由表面流 |
无网格，自然处理大变形 |
精度和稳定性较差 |
| 格子玻尔兹曼法 (LBM) |
玻尔兹曼方程离散 |
多相流、渗流 |
并行性好 |
可压缩效应限制 |
以静力弹性问题为例，有限元法将连续域 Ω 离散为 N 个单元，在每个单元内用形函数 Ni 插值位移：
u(e)=i=1∑nNiui(e)
单元刚度矩阵：
K(e)=∫Ω(e)BTCBdΩ
其中 B 是应变-位移矩阵，C 是弹性刚度矩阵。整体刚度矩阵由单元刚度矩阵装配而成：
KU=F
对于有100,000个节点的3D模型，每个节点3个自由度，K 是 300,000×300,000 的稀疏矩阵。现代FEM求解器利用稀疏矩阵技术可在几分钟内求解。
| 核心概念 |
描述 |
| 连续介质假设 |
将物质视为连续体，忽略微观结构 |
| 变形梯度 |
描述物质点从初始到当前构型的映射 |
| 应力张量 |
描述内力分布的二阶对称张量 |
| 守恒律 |
质量、动量、能量的守恒方程 |
| 本构关系 |
材料响应的数学模型 |
| 雷诺数 |
惯性力与粘性力的比值 |
| 边界层 |
壁面附近粘性效应显著的薄层 |
| 有限元法 |
偏微分方程的数值求解方法 |
-
多尺度连续介质力学：将微观结构信息（晶粒尺寸、位错密度）耦合到宏观连续模型中，桥接分子动力学与连续介质力学
-
广义连续介质理论：引入偶应力、微极理论等非经典连续介质模型，解释尺度效应（如微梁的尺寸依赖性刚度）
-
数据驱动计算力学：使用机器学习直接从实验数据中学习本构关系，替代传统的解析本构模型
-
软物质力学：研究凝胶、泡沫、液晶等软材料的非线性大变形行为，需要超弹性-粘弹性耦合本构模型
-
湍流的高精度模拟：直接数值模拟（DNS）和大涡模拟（LES）的算法改进，推动湍流基础理论发展
| 阶段 |
学习内容 |
推荐资源 |
| 基础 |
向量分析、线性代数、偏微分方程 |
本科数学教材 |
| 入门 |
弹性力学、流体力学基础 |
Timoshenko《弹性理论》、White《流体力学》 |
| 进阶 |
连续介质力学系统理论 |
Lai & Krempl《连续介质力学》 |
| 高级 |
非线性连续介质力学 |
Truesdell & Noll《理性力学》 |
| 计算 |
有限元方法、CFD |
Hughes《有限元方法》、Ferziger《CFD》 |
连续介质力学作为连接宏观与微观的力学桥梁，是材料科学、结构工程、流体力学和多个交叉学科的理论基石。掌握其核心概念和数学工具，对于深入理解自然现象和解决工程问题具有不可替代的价值。