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

学习目标
--------

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

核心概念
--------

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

本章导读
--------

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

5.1有限元分析技术
-------------------
5.1.1有限元分析方法概述
~~~~~~~~~~~~~~~~~~~~~~~~~~
有限元分析方法的思路和作法可归纳如下：

1. 物体离散化

2. 单元特性分析
 
 1）选择位移模式
 
 2）分析单元的力学性质
 
 3）计算等效结点力

3. 单元组集

4. 求解未知结点位移

5.1.2有限元分析方法中单元特性的导出方法
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
建立刚度矩阵的方法可采用：直接方法；虚功原理法；能量变分原理法。

5.1.3有限元法的解体步骤
~~~~~~~~~~~~~~~~~~~~~~~~~
1. 单元剖分和插值函数的确定

2. 单元特性分析

3. 单元组集

4. 解有限元方程

5. 计算应力

有限元基本方程
================

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

.. math::

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

其中：

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

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

5.2机械优化设计方法
--------------------
5.2.1机械优化设计问题的数学模型
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
1. 设计变量
2. 约束条件
3. 目标函数
4. 优化问题的数学模型

优化设计数学模型
------------------

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

.. math::

    \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

其中：

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

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

5.2.2机械优化设计问题的基本解法
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
1. 解析解法与数值解法

2. 优化准则法与数学规划法

3. 迭代终止条件

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 等虚拟样机软件中建立刚柔耦合多体动力学模型，预测机械振动、磨损和疲劳寿命。

复习问题
--------

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

延伸学习建议
------------

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