4、常用有限元软件操作入门：以ANSYS Workbench或Abaqus为例，界面、前处理、求解、后处理基本流程

光机结构有限元分析的核心目标，是评估结构在机械载荷、温度载荷及光学元件自重作用下的变形、应力与模态特性，进而判断光学面形变化是否超出公差范围。ANSYS Workbench与Abaqus是光机行业最主流的两种通用有限元平台，前者以流程化、参数化及多物理场耦合便捷著称，后者以非线性求解能力和单元算法丰富见长。本章围绕这两款软件，从界面认知、前处理建模、求解设置到后处理提取，梳理一套适用于光机结构分析的标准化操作流程，帮助工程师在最短时间内建立“模型—网格—边界—求解—判读”的完整闭环。

上图给出了光机结构有限元分析的整体操作链路。实际工程中，前处理阶段往往占据整个分析周期60%以上的时间，而其中网格划分与接触定义又是最容易引入人为误差的环节。下面按照软件操作的自然顺序，分别对ANSYS Workbench与Abaqus的界面布局、前处理操作、求解控制及后处理技巧进行展开说明。

4.1 软件界面与工程环境认知

核心原理：ANSYS Workbench采用“项目级”管理思路，将几何建模（DesignModeler/SpaceClaim）、网格划分（Meshing）、求解设置（Mechanical）及后处理整合在同一工作台内，通过“组件系统”拖拽方式搭建分析流程。Abaqus/CAE则采用“模块化”操作逻辑，依次进入Part、Property、Assembly、Step、Interaction、Load、Mesh、Job、Visualization等模块，每个模块对应一个独立的功能界面。

操作流程对比：

实战场景：对于光机结构中的镜筒、镜座、支撑板等零部件，建议在ANSYS Workbench中使用SpaceClaim进行几何清理——去除倒角、小孔、螺纹等不影响全局刚度的特征，然后直接进入Mechanical进行网格划分。若分析对象涉及镜片与镜座的粘接层、柔性支撑的预紧接触等非线性行为，则Abaqus的接触算法和收敛控制更为稳健。

4.2 前处理：几何准备与材料属性

核心原理：前处理的核心任务是将物理对象转化为可供求解器计算的数学模型。对于光机结构，几何模型需要准确表达光学元件的面形（球面、非球面）、镜座的支撑方式（周边支撑、背部支撑、柔性支撑）以及连接关系（粘接、压圈、螺钉）。材料属性则需涵盖弹性模量、泊松比、热膨胀系数、密度等参数，对于温度载荷分析还需输入不同温度下的材料曲线。

操作流程：

1. 几何导入与清理：在ANSYS Workbench中，通过Geometry组件导入STEP文件，进入SpaceClaim进行“修复”和“简化”。在Abaqus中，通过File→Import→Part导入几何，使用Part模块中的Repair工具修复破损面。
2. 材料参数输入：光机结构常用材料包括铝合金（6061-T6）、钛合金（TC4）、因瓦合金（4J32）、微晶玻璃（Zerodur）、碳化硅（SiC）等。以Zerodur为例，其弹性模量约为90.3GPa，泊松比0.24，热膨胀系数接近0±0.02×10⁻⁶/K，这些参数必须精确输入。
3. 接触与连接定义：镜片与镜座之间若采用粘接，可定义为绑定接触（Bonded）；若采用压圈固定，则需定义摩擦接触并设置预紧力。在Abaqus中，接触对的定义需要指定主面与从面，并设置接触刚度。

案例拆解：某反射镜组件由微晶玻璃镜坯、因瓦合金柔性支撑和铝合金镜座组成。在ANSYS Workbench中，将镜坯与柔性支撑的粘接面设为Bonded接触，柔性支撑与镜座之间设为No Separation接触（允许微小滑移但不分离）。材料参数分别输入三种材料的弹性模量与热膨胀系数。网格划分时，对镜坯光学面设置0.5mm的单元尺寸，镜座主体设置2mm单元尺寸，过渡区域采用自动过渡。

4.3 网格划分策略与质量检查

核心原理：网格划分的本质是用离散单元逼近连续体。对于光机结构，六面体单元（Hex20）在弯曲变形和应力梯度较大的区域具有更高的精度，而四面体单元（Tet10）在复杂几何中更容易生成。网格密度需要根据应力梯度或变形梯度进行自适应调整——光学镜面附近需要较密网格以准确提取面形变化，而远离光学面的结构件可适当稀疏。

判断方法：

• 单元质量指标：偏斜度（Skewness）应小于0.85，正交质量（Orthogonal Quality）应大于0.15，翘曲率（Warping Factor）应小于30°。
• 网格无关性验证：取两倍加密网格计算同一工况，若镜面最大变形变化小于2%，则认为当前网格密度满足要求。
• 单元类型选择：光机结构首选二阶六面体单元（C3D20R或Solid186），对于薄壁结构可选用壳单元（S8R或Shell181），对于粘接层等薄层结构需保证厚度方向至少2层单元。

操作流程（以ANSYS Workbench为例）：

1. 在Mechanical中选中几何体，右键插入Method，选择Hex Dominant或Multizone。
2. 对镜面、粘接面等关键面插入Sizing，设置单元尺寸为0.5~1mm。
3. 对薄壁区域插入Thin Sweep，设置层数为3~5层。
4. 生成网格后，查看Mesh Metric中的Skewness和Orthogonal Quality。

实战场景：某卡塞格林望远镜的主镜口径300mm，厚度40mm，采用背部三点支撑。在Abaqus中，将主镜划分为六面体网格，径向单元数30，轴向单元数8，背部支撑区域局部加密至0.3mm。镜座采用四面体网格，与主镜接触区域设置面网格过渡。计算模态时，一阶固有频率为385Hz，与理论估算值378Hz偏差1.8%，验证了网格精度。

4.4 求解设置与收敛控制

核心原理：求解设置包括分析类型选择、边界条件施加、载荷步定义和求解器参数控制。光机结构分析常见类型有：静力分析（自重+温度）、模态分析（固有频率与振型）、热-结构耦合分析（温度场→变形）、屈曲分析（薄壁镜筒稳定性）。边界条件必须真实模拟结构的安装方式——过约束会低估变形，欠约束则导致刚体位移。

操作流程：

1. 分析类型：在ANSYS Workbench中，从Toolbox拖入Static Structural、Modal或Thermal-Structural耦合分析。在Abaqus中，在Step模块中创建Static, General或Frequency分析步。
2. 边界条件：固定镜座底面或安装法兰的螺栓孔位置，施加标准地球重力加速度9806.65mm/s²，温度载荷根据实际工况输入均匀温度场或梯度温度场。
3. 求解器参数：ANSYS中设置Large Deflection为On（若变形超过厚度的1/10），Abaqus中设置Nlgeom为ON。对于接触问题，设置Contact Stabilization以辅助初始收敛。
4. 求解监控：在ANSYS中查看Solution Information中的力收敛曲线，在Abaqus中监控Job Monitor中的迭代残差。

案例拆解：某空间相机主镜组件需承受1g自重和20℃均匀温升。在ANSYS Workbench中，设置两个载荷步：第一步施加自重，第二步叠加温度载荷。求解时发现接触不收敛，原因是镜座与柔性支撑之间的摩擦接触初始穿透过大。解决方案：在Contact Tool中调整Pinball Region为0.1mm，并开启Automatic Contact Stabilization，阻尼系数设为0.01。重新求解后，经过12次迭代收敛，镜面最大变形为0.023μm，满足设计指标0.05μm。

4.5 后处理：变形提取与面形评价

核心原理：后处理的目标是从海量节点数据中提取光学设计所需的关键指标。对于光机结构，最重要的输出是光学镜面的节点位移场，需要将其转化为面形误差（PV值、RMS值）和刚体位移（平移、倾斜、离焦）。此外，还需检查结构应力是否超过材料许用值，模态频率是否避开卫星或地面设备的激励频率。

操作流程：

1. 变形提取：在ANSYS中，插入Deformation→Total，选择镜面节点集，导出节点坐标与位移。在Abaqus中，进入Visualization模块，创建Contour Plot，通过Tools→Query→Probe Values提取节点位移。
2. 面形拟合：将镜面节点位移导出为文本文件，使用MATLAB或Python脚本进行Zernike多项式拟合。去除平移项（Z1）、倾斜项（Z2、Z3）和离焦项（Z4）后，剩余高阶项即为面形误差。
3. 应力校核：查看Equivalent Stress（Von Mises）云图，对比材料屈服强度。对于脆性材料（微晶玻璃、SiC），应使用最大主应力准则，安全系数取3~5。
4. 模态评价：提取前6阶固有频率，对比基频是否大于100Hz（空间相机通常要求基频＞100Hz以避免与卫星平台共振）。

实战场景：某离轴三反光学系统的主镜在1g自重下的镜面变形云图显示，最大变形出现在镜面边缘，数值为0.18μm。通过Zernike拟合，去除刚体位移后，面形误差PV=0.032μm，RMS=0.006μm，满足光学设计公差（PV＜λ/10，λ=632.8nm）。同时，镜座最大应力为45MPa，远低于铝合金屈服强度276MPa，结构安全裕度充足。

常见误区与复盘：

• 误区一：网格越密越好。实际上，当网格密度达到一定程度后，计算精度提升有限，但计算时间成倍增加。建议先做粗网格试算，再逐步加密关键区域。
• 误区二：忽略接触刚度的影响。在Abaqus中，接触刚度过低会导致穿透过大，过高则引起收敛振荡。建议初始接触刚度设为相邻单元刚度的1~10倍。
• 误区三：直接使用默认求解器设置。对于光机结构的大变形或接触问题，应开启大变形开关并调整求解器参数，否则结果可能严重偏离实际。
• 误区四：后处理只关注最大应力。光机结构更应关注光学面的变形形态，而非局部应力峰值。应力集中可通过局部圆角优化解决，而面形误差则需要从支撑方案上重新设计。

通过本章的学习，读者应能够独立完成一个典型光机结构从几何导入到面形评价的全流程分析。ANSYS Workbench与Abaqus虽然操作界面不同，但核心分析逻辑一致——前处理决定模型准确性，求解设置决定收敛性，后处理决定结果可用性。建议初学者先以简单镜座模型为对象，完整走通上述流程，再逐步增加接触、温度、模态等复杂工况。