2、有限元基础回顾：弹性力学基本方程、虚功原理、单元类型与形函数、刚度矩阵组装

光机结构有限元分析的核心，是将连续弹性体离散为有限个单元，通过节点自由度建立代数方程组求解。这一过程的理论根基，正是弹性力学基本方程与虚功原理。本章作为后续动力学、热变形、优化设计的共同基础，重点梳理从弹性力学控制方程到有限元离散格式的完整逻辑链，并针对光机结构常用单元类型与刚度矩阵组装流程做工程化解读。理解本章内容，是正确设置边界条件、判断网格质量、解读应力结果的前提。

2.1 弹性力学基本方程：光机结构分析的物理基础

核心原理

弹性力学基本方程描述连续弹性体在外部载荷作用下的力学响应，包含三组控制方程：平衡方程（应力与体力的关系）、几何方程（应变与位移的关系）、本构方程（应力与应变的关系）。对于光机结构，材料通常处于线弹性小变形范围，各向同性材料的本构方程由弹性模量 E 和泊松比 ν 两个独立常数确定。在笛卡尔坐标系下，平衡方程为 σ_ij,j + f_i = 0，几何方程为 ε_ij = (u_i,j + u_j,i)/2，本构方程为 σ_ij = λε_kkδ_ij + 2με_ij，其中 λ 和 μ 为拉梅常数。

光机结构分析的特殊性在于：光学镜面通常要求亚微米甚至纳米级的面形精度，而结构尺寸从几毫米到数米不等。因此，弹性力学方程在光机领域的应用需要特别关注边界条件的准确性——支撑方式（周边支撑、背部支撑、中心支撑）直接决定镜面变形模式。例如，一个直径 300mm 的反射镜，采用三点支撑与采用周边支撑，在相同重力载荷下的 PV 值可能相差 5~10 倍。

判断方法与工程应用

在实际光机结构分析中，判断弹性力学方程适用性的关键指标包括：

• 材料线性度：应力水平是否远低于材料屈服强度（通常要求安全系数 ≥ 2.5）
• 几何线性度：最大位移是否小于结构特征尺寸的 1/20（如镜体厚度的 5%）
• 应变范围：最大应变是否小于 0.1%（玻璃、微晶玻璃等脆性材料更严格）

2.2 虚功原理：从微分方程到积分弱形式

核心原理

虚功原理（Principle of Virtual Work）是有限元方法的理论基石。它表述为：对于处于平衡状态的弹性体，任意满足位移边界条件的虚位移 δu 下，内力虚功等于外力虚功。数学表达式为：

∫V δεᵀσ dV = ∫V δuᵀf dV + ∫S δuᵀT dS

其中，δε 为虚应变，σ 为真实应力，f 为体力，T 为面力。虚功原理将弹性力学中要求每一点都满足的强形式微分方程，转化为对整个求解域满足积分形式的弱方程。这一转化使得我们可以用有限个离散自由度近似描述连续位移场，而不必逐点求解偏微分方程。

从变分原理角度看，虚功原理等价于最小势能原理：在所有满足位移边界条件的可能位移场中，真实位移场使系统总势能取最小值。对于光机结构，这意味着有限元解总是偏刚的——即计算位移小于真实位移。这一特性在镜面变形分析中需要特别注意：如果有限元计算的镜面 PV 值接近光学设计公差，实际结构可能已经超差。

操作流程与工程映射

在有限元软件中，虚功原理的应用流程为：

1. 建立几何模型并划分网格，定义单元类型和材料属性
2. 施加边界条件（约束和载荷），软件自动构造虚位移场
3. 计算单元应变矩阵 B 和本构矩阵 D，形成单元刚度矩阵
4. 组装整体刚度矩阵，求解节点位移
5. 后处理计算应变、应力，并提取光学面变形数据

2.3 单元类型与形函数：光机结构网格划分的核心选择

核心原理

形函数 N_i 是有限元方法中连接节点位移与单元内任意点位移的桥梁。对于等参单元，几何坐标和位移场采用相同的形函数插值：x = ΣN_ix_i，u = ΣN_iu_i。形函数必须满足：在节点 i 处取值为 1，在其他节点处取值为 0；且单元内任意点的形函数之和恒等于 1（保证刚体位移的精确表示）。

光机结构分析中常用的单元类型包括：

• 六面体单元（C3D8/C3D20）：适用于镜体、镜筒等实体结构。C3D20（20节点二次单元）精度高，但计算成本约为 C3D8 的 4~6 倍
• 四面体单元（C3D4/C3D10）：适用于复杂几何的自动网格划分。C3D4 线性四面体过刚且存在体积锁定问题，不建议用于光机结构；C3D10 二次四面体精度可接受
• 壳单元（S4/S8）：适用于薄壁镜筒、加强筋等板壳结构。要求厚度方向尺寸远小于面内尺寸（通常 < 1/10）
• 梁单元（B31/B32）：适用于支撑杆、桁架结构。光机结构中常用于简化支撑腿分析

选择判断方法

案例拆解：某卡塞格林望远镜副镜支撑结构

副镜直径 120mm，采用三翼梁支撑结构。初始方案使用 C3D4 四面体单元自动网格划分，分析得到一阶固有频率 215Hz。改用 C3D10 二次四面体后，一阶频率降至 187Hz，差异达 13%。进一步采用 C3D20R 六面体单元精细网格，结果为 182Hz。最终确认 C3D4 因过刚导致频率高估约 18%。该案例说明：单元类型选择对动力学分析结果的影响远超预期，光机结构应避免使用线性四面体单元。

2.4 刚度矩阵组装：从单元到整体的集成逻辑

核心原理

整体刚度矩阵 K 由所有单元刚度矩阵 k^e 按照节点自由度编号规则叠加而成。组装过程遵循直接刚度法：将每个单元的刚度系数 k_ij^e 添加到整体矩阵的对应位置 K_I(i),J(j) 中，其中 I(i) 和 J(j) 是单元局部节点 i、j 对应的全局自由度编号。整体刚度矩阵具有对称性、带状分布和奇异性（需施加边界条件消除刚体位移）。

对于光机结构，刚度矩阵组装的关键在于多部件连接关系的正确处理。镜体与支撑结构之间通常采用胶接、螺栓连接或柔性环节，这些连接在有限元模型中表现为：

• 共节点：网格划分时保证连接处节点重合，适用于胶接或焊接
• Tie 约束：不同网格之间建立绑定关系，适用于面-面连接
• 弹簧/阻尼单元：模拟柔性环节或螺栓连接的等效刚度
• 接触对：考虑间隙和摩擦，适用于非永久连接

操作流程与工程要点

// 组装流程伪代码（以 3 节点三角形单元为例）
for each element e:
    compute k_e (6×6 矩阵)
    for i = 1 to 3:
        for j = 1 to 3:
            I = global_dof(node_i, dof_index)
            J = global_dof(node_j, dof_index)
            K[I][J] += k_e[i][j]

在光机结构分析中，刚度矩阵组装的常见工程问题包括：

1. 自由度不匹配：实体单元（3 个平动自由度）与壳单元（6 个自由度）连接时，需使用 MPC（多点约束）或 Solid-Shell 耦合
2. 刚性连接过强：使用 RBE2 刚性单元模拟螺栓连接时，会人为增加连接刚度，导致局部应力失真
3. 网格不连续：不同部件网格密度差异过大时，Tie 约束可能引入虚假应力集中

常见误区与总结复盘