第5章计算机辅助工程分析¶

学习目标¶

理解有限元分析（FEA）的基本思想、单元特性推导方法与求解步骤
掌握机械优化设计的数学模型构建方法及常用数值解法
了解虚拟样机技术在产品设计中的角色和实现方式

核心概念¶

有限元法（FEM/FEA）：将连续体离散为有限个单元，通过单元刚度矩阵组集和方程求解得到近似场量的数值方法
刚度矩阵 [K]：反映单元节点力与节点位移之间关系的矩阵，整体刚度矩阵由单元刚度按节点编号装配
设计变量 / 目标函数 / 约束：机械优化设计数学模型的三要素——决策量、性能指标、限制条件
优化准则法：基于工程准则（如满应力准则）迭代求优，收敛快但适用面较窄
虚拟样机（Virtual Prototyping）：用软件建立机械系统三维力学模型，在物理样机制造前完成性能仿真与改进
CAE：计算机辅助工程分析的统称，FEA、优化、动力学仿真均属其范畴

本章导读¶

本章聚焦 CAE（计算机辅助工程分析）的三大支柱：有限元分析求解结构场量、机械优化设计寻找最佳参数、虚拟样机技术在物理样机制造前完成性能评估。三者结合使工程设计从"经验+试验"演进为"仿真驱动设计"，是现代产品研发不可或缺的环节。

5.1有限元分析技术¶

5.1.1有限元分析方法概述¶

有限元分析方法的思路和作法可归纳如下：

物体离散化
单元特性分析

1）选择位移模式

2）分析单元的力学性质

3）计算等效结点力

单元组集
求解未知结点位移

5.1.2有限元分析方法中单元特性的导出方法¶

建立刚度矩阵的方法可采用：直接方法；虚功原理法；能量变分原理法。

5.1.3有限元法的解体步骤¶

单元剖分和插值函数的确定
单元特性分析
单元组集
解有限元方程
计算应力

有限元基本方程¶

有限元分析的基本方程可以表示为：

\[[K]\{u\} = \{F\}\]

其中：

\([K]\) — 整体刚度矩阵（Global Stiffness Matrix）
\(\{u\}\) — 节点位移向量（Nodal Displacement Vector）
\(\{F\}\) — 节点载荷向量（Nodal Force Vector）

单元刚度矩阵 \([k^e]\) 通过虚功原理或能量变分法导出，整体刚度矩阵由各单元刚度矩阵按节点编号组装而成。

5.2机械优化设计方法¶

5.2.1机械优化设计问题的数学模型¶

设计变量
约束条件
目标函数
优化问题的数学模型

优化设计数学模型¶

机械优化设计问题的标准数学模型可表述为：

\[ \begin{align}\begin{aligned}\min f(\mathbf{x})\\\text{s.t.} \quad g_i(\mathbf{x}) \leq 0, \quad i = 1, 2, \ldots, m\\\quad \quad h_j(\mathbf{x}) = 0, \quad j = 1, 2, \ldots, l\\\quad \quad x_k^{\min} \leq x_k \leq x_k^{\max}, \quad k = 1, 2, \ldots, n\end{aligned}\end{align} \]

其中：

\(\mathbf{x} = [x_1, x_2, \ldots, x_n]^T\) — 设计变量向量
\(f(\mathbf{x})\) — 目标函数（如重量最小、成本最低、刚度最大等）
\(g_i(\mathbf{x}) \leq 0\) — 不等式约束（如应力约束、变形约束）
\(h_j(\mathbf{x}) = 0\) — 等式约束（如几何约束、平衡条件）
\(x_k^{\min}, x_k^{\max}\) — 设计变量的上下界

常用的求解方法包括：数学规划法（如梯度法、牛顿法）、优化准则法、遗传算法等。

5.2.2机械优化设计问题的基本解法¶

解析解法与数值解法
优化准则法与数学规划法
迭代终止条件

5.3虚拟样机技术¶

利用软件建立机械系统的三维实体模型和力学模型，分析和评估系统的性能，从而为物理样机的设计和制造提供依据。

5.3.1虚拟样机技术的基本概念¶

5.3.2虚拟样机技术的形成和发展¶

5.3.3虚拟样机技术的相关技术¶

关键术语¶

有限元法（Finite Element Method, FEM）
单元刚度矩阵（Element Stiffness Matrix）
整体刚度矩阵（Global Stiffness Matrix）
节点位移向量（Nodal Displacement Vector）
虚功原理（Principle of Virtual Work）
能量变分原理（Energy Variational Principle）
设计变量（Design Variable）
目标函数（Objective Function）
约束条件（Constraints）
数学规划法（Mathematical Programming）
优化准则法（Optimality Criteria Method）
虚拟样机（Virtual Prototype）
刚度/强度分析（Stiffness/Strength Analysis）
动力学仿真（Dynamic Simulation）

工程应用场景¶

结构强度仿真: 对汽车保险杠、飞机机翼翼盒等大型薄壁件进行网格划分与边界条件设置，通过 FEA 计算应力/应变分布，验证是否满足强度要求。
拓扑优化: 在给定设计空间与载荷条件下，借助数学规划法优化材料分布，得到轻量化的最优结构（广泛应用于机械臂、卫星支架）。
机械系统动力学仿真: 在 ADAMS、RecurDyn 等虚拟样机软件中建立刚柔耦合多体动力学模型，预测机械振动、磨损和疲劳寿命。

复习问题¶

有限元法的基本求解步骤是什么？整体刚度矩阵是如何由单元刚度矩阵组集而成的？
机械优化设计数学模型由哪三个要素构成？分别举例说明。
什么是虚拟样机？它与传统物理样机相比有哪些优势？

延伸学习建议¶

在 ANSYS Workbench 或开源的 CalculiX 中完成一个简单的悬臂梁静力学分析，亲历建模-网格-后处理全流程
学习 Python 中 scipy.optimize 模块，掌握连续/离散优化的基本解法（梯度法、遗传算法等）
了解机器学习代理模型（Kriging、神经网络）如何加速优化仿真，提升 CAE 效率

第5章 计算机辅助工程分析¶