CAE 计算机辅助工程 - 工业仿真的核心 / Computer-Aided Engineering as the Core of Industrial Simulation
日期：2026-04-29 标签：#CAE #有限元 #FEA #仿真 #多物理场 前置知识：CAD 计算机辅助设计
本章目标
- 理解 CAE 的定义和分类
- 掌握有限元分析(FEA)的基本原理
- 了解多物理场仿真的类型
- 理解常见仿真分析方法
- 认识 CAE 在产品开发中的价值
第1部分：CAE 概述
1.1 什么是 CAE？
CAE = Computer-Aided Engineering（计算机辅助工程）
使用计算机软件对产品设计进行性能分析和验证的技术。
┌─────────────────────────────────────────────────────────────────────────────┐
│ CAE 的核心价值 │
├─────────────────────────────────────────────────────────────────────────────┤
│ │
│ 传统方式： │
│ ┌─────────┐ ┌─────────┐ ┌─────────┐ ┌─────────┐ │
│ │ 设计 │ ──→ │ 制造样机 │ ──→ │ 物理测试 │ ──→ │ 发现问题 │ │
│ └─────────┘ └─────────┘ └─────────┘ └────┬────┘ │
│ │ │
│ ↑───────────────────────────────────────────────────┘ │
│ │
│ CAE方式： │
│ ┌─────────┐ ┌─────────┐ ┌─────────┐ │
│ │ 设计 │ ──→ │ 虚拟仿真 │ ──→ │ 发现问题 │ │
│ └─────────┘ └─────────┘ └────┬────┘ │
│ │ │
│ ↓ │
│ 修改设计 │
│ │ │
│ ┌─────────────────────────────────┘ │
│ ↓ │
│ ┌─────────┐ ┌─────────┐ ┌─────────┐ │
│ │ 验证通过 │ ←── │ 物理测试 │ ←── │ 样机制造 │ │
│ └─────────┘ └─────────┘ └─────────┘ │
│ │
│ CAE 让"试错"发生在虚拟世界，大幅节省时间和成本 │
│ │
└─────────────────────────────────────────────────────────────────────────────┘1.2 CAE 的分类
┌─────────────────────────────────────────────────────────────────────────────┐
│ CAE 仿真类型 │
├─────────────────────────────────────────────────────────────────────────────┤
│ │
│ 结构仿真 (Structural) │
│ ├── 静力学分析 (Static) │
│ ├── 动力学分析 (Dynamic) │
│ │ ├── 模态分析 (Modal) │
│ │ ├── 响应谱分析 (Response Spectrum) │
│ │ ├── 谐响应分析 (Harmonic Response) │
│ │ └── 瞬态动力学 (Transient) │
│ ├── 屈曲分析 (Buckling) │
│ ├── 疲劳分析 (Fatigue) │
│ └── 非线性分析 (Nonlinear) │
│ │
│ 热分析 (Thermal) │
│ ├── 稳态热分析 (Steady-State) │
│ └── 瞬态热分析 (Transient) │
│ │
│ 流体力学 (CFD) │
│ ├── 内流分析 (Internal Flow) │
│ ├── 外流分析 (External Flow) │
│ ├── 多相流 (Multi-Phase) │
│ └── 传热耦合 (Heat Transfer) │
│ │
│ 电磁仿真 (Electromagnetic) │
│ ├── 静磁场分析 (Magnetostatic) │
│ ├── 静电场分析 (Electrostatic) │
│ ├── 电磁场分析 (EM Field) │
│ └── 信号完整性 (Signal Integrity) │
│ │
│ 多物理场耦合 (Multi-Physics) │
│ ├── 流固耦合 (FSI) │
│ ├── 热结构耦合 (Thermo-Mechanical) │
│ └── 电磁热耦合 (Electro-Thermal) │
│ │
└─────────────────────────────────────────────────────────────────────────────┘第2部分：有限元分析基础
2.1 什么是有限元？
有限元 = Finite Element Method (FEM)
一种数值计算方法，用于求解复杂工程问题。
核心思想：
连续体问题 → 离散化 → 求解代数方程组 → 得到近似解
把一个复杂形状分割成很多小"积木"
每个积木假设是简单形状（三角形、四边形、四面体、六面体）
求解每个积木，然后把结果拼起来类比理解：
假设要测量一个大陆的面积：
方法1（解析法）：
找一本地理书，上面写着精确面积
（相当于解析解，不是总能找到）
方法2（有限元法）：
把地图画上格子
数有多少格子
格子面积 × 格子数 = 估计面积
格子越细，结果越准
（这就是有限元的思想！）2.2 有限元分析流程
┌─────────────────────────────────────────────────────────────────────────────┐
│ 有限元分析流程 │
├─────────────────────────────────────────────────────────────────────────────┤
│ │
│ ┌─────────────────────────────────────────────────────────────────────┐ │
│ │ 1. 前处理 (Pre-processing) │ │
│ │ ├── 导入几何模型 │ │
│ │ ├── 几何清理与简化 │ │
│ │ ├── 网格划分（离散化） │ │
│ │ ├── 定义材料属性 │ │
│ │ └── 施加边界条件和载荷 │ │
│ └─────────────────────────────────────────────────────────────────────┘ │
│ ↓ │
│ ┌─────────────────────────────────────────────────────────────────────┐ │
│ │ 2. 求解 (Solution) │ │
│ │ ├── 组装刚度矩阵 │ │
│ │ ├── 施加边界条件 │ │
│ │ ├── 求解线性/非线性方程组 │ │
│ │ └── 计算应力、应变、位移等 │ │
│ └─────────────────────────────────────────────────────────────────────┘ │
│ ↓ │
│ ┌─────────────────────────────────────────────────────────────────────┐ │
│ │ 3. 后处理 (Post-processing) │ │
│ │ ├── 结果可视化（云图、矢量图） │ │
│ │ ├── 结果动画演示 │ │
│ │ ├── 数据提取与报告 │ │
│ │ └── 结果验证与评估 │ │
│ └─────────────────────────────────────────────────────────────────────┘ │
│ │
└─────────────────────────────────────────────────────────────────────────────┘2.3 单元类型详解
一维单元：
梁单元 (Beam)
●━━━━━━━━━━━━━━━━━━●
用途：桁架、桥梁、建筑框架
变形：弯曲二维单元：
三角形单元 (Tria) 四边形单元 (Quad)
● ●───●
╱ ╲ ╱ ╲
╱ ╲ ●─────●
╱ ╲ ╲ ╱
●───────● ╲ ╱
●───●
3节点 4节点或8节点
用途：板壳、薄壁结构三维单元：
四面体单元 (Tet) 六面体单元 (Hex)
● ●───●
╱│╲ ╱│ ╲│
╱ │ ╲ ╱ │ ╲ │
╱ │ ╲ ╱ │ ╲ │
╱ │ ╲ ●───┼────┼─●
●───┼───● │╲ │ │╱╲
╲ │ ╱ │ ╲─┼────┼─╲
╲ │ ╱ │ ╲│ │ ╲
╲│╱ ●───●────┼──●
● │ ╱ │ ╱
│ ╱ │ ╱
●──────●
用途：实体结构
注意：六面体精度高但难生成2.4 网格质量
好的网格 vs 差的网格：
┌─────────────────────────────────────────────────────────────────────────────┐
│ 网格质量评估 │
├─────────────────────────────────────────────────────────────────────────────┤
│ │
│ 好的网格： │
│ ├── 单元形状规则（正三角形、正方形） │
│ ├── 过渡平缓，不突变 │
│ └── 关键区域加密适当 │
│ │
│ ✓ ✗ │
│ /\ /│\ │
│ / \ / │ \ │
│ / \ / │ \ │
│ ─────── / │ \ │
│ │
│ 差的网格： │
│ ├── 单元畸形（狭长、尖角） │
│ ├── 角度扭曲严重 │
│ └── 单元大小突变 │
│ │
│ ✗ ✗ │
│ /│\ /╲/╲ │
│ / │ \ /╲ ╱╲ │
│ / │ \ ╱╲ ╱╲╱ │
│ ─────── ╲╱╲╱ │
│ │
└─────────────────────────────────────────────────────────────────────────────┘网格质量指标：
| 指标 | 描述 | 优良范围 |
|---|---|---|
| 长宽比 (Aspect Ratio) | 最长边/最短边 | < 3 |
| 歪斜度 (Skewness) | 偏离正多边形程度 | < 0.5 |
| 雅可比 (Jacobian) | 单元畸变程度 | > 0.6 |
| 锥度 (Taper) | 四边形对角线比 | > 0.5 |
第3部分：结构仿真详解
3.1 静力学分析
定义：分析结构在静止载荷作用下的响应。
基本方程：
[K] {u} = {F}
其中：
[K] = 刚度矩阵
{u} = 位移向量
{F} = 载荷向量分析步骤：
1. 建立几何模型
2. 定义材料（弹性模量E、泊松比ν、密度ρ）
3. 划分网格
4. 施加边界条件（固定约束）
5. 施加载荷（力、压力、引力）
6. 求解
7. 查看结果（位移、应力、应变）实例：悬臂梁分析
问题：
一根钢梁，一端固定，另一端受向下力F=1000N
已知：
长度 L = 1m
截面 矩形 b×h = 0.05m × 0.1m
材料 钢 E = 200 GPa
求：
最大位移和最大应力
F = 1000N
↓
┌──────────────────────────┐
│ │
│ ══════════════════════ │ ← 梁
│ │
└──────────────────────────┘
▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓ ← 固定端
▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓▓
求解结果：
最大位移 = F×L³/(3×E×I) ≈ 2.67 mm
最大应力 = F×L/(S×Z) ≈ 60 MPa3.2 模态分析
定义：分析结构的固有振动特性（频率和振型）。
为什么重要？
避免共振！
每个物体都有自己的"音乐"（固有频率）
当外界振动频率 = 固有频率 → 共振 → 剧烈振动 → 损坏
著名的共振例子：
- 塔科马海峡大桥被风吹垮（1940）
- 航天飞机助推器振动问题固有频率计算：
[K] {φ} = λ [M] {φ}
其中：
[K] = 刚度矩阵
[M] = 质量矩阵
λ = 特征值（ω² = 固有频率的平方）
{φ} = 振型向量振型示例：
一阶振型 二阶振型 三阶振型
(第一音) (第二音) (第三音)
↓↓↓ ↓ ↓ ↓ ↓ ↓
███ █ █ █ █
███ █ █ █ █
███ █ █ █ █
═══════════ ═══════════ ═══════════
███ █ █ █ █
███ █ █ █ █
███ █ █ █ █
f₁ ≈ 12 Hz f₂ ≈ 35 Hz f₃ ≈ 65 Hz
(频率越高，振型越复杂)3.3 非线性分析
什么是非线性？
线性分析假设：
- 材料遵循胡克定律（应力∝应变）
- 小变形（位移不影响几何形状）
非线性分析处理：
- 材料非线性：塑性、粘弹性、超弹性
- 几何非线性：大变形、屈曲、接触
- 边界非线性：接触、摩擦、螺栓预紧
┌─────────────────────────────────────────────────────────────────────────────┐
│ 线性 vs 非线性 │
├─────────────────────────────────────────────────────────────────────────────┤
│ │
│ 线性问题： │
│ ┌──────────────┐ │
│ │ F │ │
│ │ ↑ │ │
│ │ │ ╱ │ P = K·u │
│ │ │╱ │ (直线关系) │
│ │ └──────────→ u │
│ │ 0 │ │
│ └──────────────┘ │
│ 特点：直接求解，一次搞定 │
│ │
│ 非线性问题： │
│ ┌──────────────┐ │
│ │ F │ │
│ │ ↑ │ │
│ │ │ ╱ │ │
│ │ │ ╱ │ 载荷需要分步施加 │
│ │ │ ╱ │ 用增量-迭代法求解 │
│ │ │ ╱ │ │
│ │ │╱ │ │
│ │ └──────────→ u │
│ │ 0 │ │
│ └──────────────┘ │
│ 特点：需要迭代，可能不收敛 │
│ │
└─────────────────────────────────────────────────────────────────────────────┘第4部分：其他仿真类型
4.1 计算流体动力学 (CFD)
CFD = Computational Fluid Dynamics
模拟流体（气体、液体）的流动和传热。
典型应用：
航空航天：飞机、汽车外流场分析
能源行业：换热器、管道流动
建筑行业：通风、空调、气流组织
水利工程：河流水库、阀门流动CFD vs 实验：
| 对比项 | CFD | 实验 |
|---|---|---|
| 成本 | 低 | 高 |
| 周期 | 短 | 长 |
| 数据 | 完整全场 | 有限测点 |
| 精度 | 依赖模型 | 真实可靠 |
| 可视化 | 丰富 | 有限 |
4.2 热分析
类型：
- 稳态热分析：温度场不随时间变化
- 瞬态热分析：温度场随时间变化
热传递方式：
┌─────────────────────────────────────────────────────────────────────────────┐
│ 热传递三种方式 │
├─────────────────────────────────────────────────────────────────────────────┤
│ │
│ 1. 导热 (Conduction) │
│ ┌─────────┐ 热 ┌─────────┐ │
│ │ 热端 │ ───────→ │ 冷端 │ │
│ └─────────┘ └─────────┘ │
│ 公式：q = -k·A·dT/dx (傅里叶定律) │
│ │
│ 2. 对流 (Convection) │
│ ┌─────────┐ ┌─────────┐ │
│ │ 壁面 │ ═══ 热 + 流动 ═══ │ 流体 │ │
│ └─────────┘ └─────────┘ │
│ 公式：q = h·A·(Tw - T∞) (牛顿冷却定律) │
│ │
│ 3. 辐射 (Radiation) │
│ ┌─────────┐ ┌─────────┐ │
│ │ ☀️ │ ════════════════ │ │ │
│ │ (热源) │ 电磁波 │ (受体) │ │
│ └─────────┘ └─────────┘ │
│ 公式：q = ε·σ·A·T⁴ (斯特藩-玻尔兹曼定律) │
│ │
└─────────────────────────────────────────────────────────────────────────────┘4.3 电磁仿真
应用领域：
电机设计：永磁电机、变压器
天线设计：手机天线、雷达
电子产品：PCB电磁兼容
医疗设备：MRI、CT分析类型：
- 静磁场：永磁体、电磁铁
- 静电场：电容、电场分布
- 低频电磁：感应加热、涡流
- 高频电磁：天线、微波
第5部分：CAE 软件体系
5.1 CAE 软件架构
┌─────────────────────────────────────────────────────────────────────────────┐
│ CAE 软件架构 │
├─────────────────────────────────────────────────────────────────────────────┤
│ │
│ ┌─────────────────────────────────────────────────────────────────────┐ │
│ │ 前处理 (Preprocessor) │ │
│ │ 几何处理 | 网格划分 | 边界条件 | 材料定义 │ │
│ └─────────────────────────────────────────────────────────────────────┘ │
│ ↓ │
│ ┌─────────────────────────────────────────────────────────────────────┐ │
│ │ 求解器 (Solver) │ │
│ │ ┌──────────┐ ┌──────────┐ ┌──────────┐ ┌──────────┐ │ │
│ │ │结构求解器│ │流体求解器│ │热求解器 │ │电磁求解器│ │ │
│ │ └──────────┘ └──────────┘ └──────────┘ └──────────┘ │ │
│ │ 线性求解 | 非线性求解 | 特征值求解 | 优化算法 │ │
│ └─────────────────────────────────────────────────────────────────────┘ │
│ ↓ │
│ ┌─────────────────────────────────────────────────────────────────────┐ │
│ │ 后处理 (Postprocessor) │ │
│ │ 结果可视化 | 云图动画 | 数据提取 | 报告生成 │ │
│ └─────────────────────────────────────────────────────────────────────┘ │
│ │
└─────────────────────────────────────────────────────────────────────────────┘5.2 主流 CAE 软件一览
| 软件 | 公司 | 专长 | 特点 |
|---|---|---|---|
| Ansys Mechanical | Ansys | 通用结构 | 行业标准 |
| Abaqus | DS SIMULIA | 非线性仿真 | 强大显式 |
| ANSYS Fluent | Ansys | CFD | 流体霸主 |
| ANSYS Maxwell | Ansys | 电磁 | 低频权威 |
| STAR-CCM+ | Siemens | CFD | 网格自动 |
| NX Nastran | Siemens | 结构 | NASA出身 |
| OptiStruct | Altair | 拓扑优化 | 轻量化 |
| HyperWorks | Altair | 多学科 | 开放平台 |
核心总结
总结1：CAE 分类
CAE = 结构 + 热 + 流 + 电磁 + 多物理场耦合
结构：静力学、动力学、非线性、疲劳
热：稳态、瞬态、辐射、耦合
CFD：内流、外流、多相、传热
电磁：静磁、静电、低频、高频总结2：有限元核心
FEM = 离散化 + 矩阵求解
连续体 → 分割成单元 → 组装方程组 → 求解 → 结果总结3：CAE 流程
前处理（80%时间）→ 求解（10%时间）→ 后处理（10%时间）章节测试
测试1：CAE 分类
CFD、CAE、FEA 三者是什么关系？
测试2：有限元思想
用一句话解释有限元的核心思想。
测试3：静力学方程
有限元静力学分析的基本方程 [K]{u}={F} 中各符号含义？
测试4：为什么需要网格
为什么有限元分析需要网格？
测试5：CAE vs 实验
CAE 仿真相比物理实验有哪些优势和劣势？
参考答案
测试1答案
- CAE：计算机辅助工程的统称
- FEA：有限元分析，是 CAE 的一个子集
- CFD：计算流体动力学，也是 CAE 的一个子集
测试2答案
把复杂连续体分割成有限个简单单元，在每个单元上用近似函数描述物理量，然后组装求解。
测试3答案
- [K]：整体刚度矩阵
- {u}：节点位移向量
- {F}：外载荷向量
测试4答案
因为实际的微分方程很难直接求解，有限元法通过离散化把连续问题变成离散问题，从而能用计算机数值求解。
测试5答案
优势：成本低、周期短、数据完整、可视化丰富、可做"不可能"实验 劣势：依赖数学模型、计算精度有限、复杂问题收敛困难
相关笔记
- [[02 - CAD 计算机辅助设计]] - 仿真对象
- [[04 - 前处理与后处理]] - 仿真细节
- [[05 - 求解器基础]] - 核心计算引擎
下一步学习
- [ ] 阅读 04 - 前处理与后处理
学习状态：完成