3. 非线性接触仿真基础

在光机结构分析中，很多工程师习惯把整机简化为线性静力学问题：给定载荷，求解刚度矩阵 K·u = F，得到位移与应力。这个做法在螺栓预紧、镜座过盈、柔性铰链大变形、胶层塑性等场景下会迅速失效。原因不是求解器不够强，而是问题本身同时包含接触非线性、几何非线性与材料非线性，三者往往不是独立叠加，而是相互耦合、彼此放大。本节先厘清三类非线性的物理来源，再说明它们在光机结构中的耦合方式，最后给出判断"能否用线性静力学替代"的工程准则。

一、三类非线性各自的物理来源

接触非线性来自边界条件随状态变化。两个零件在未接触时互不影响，一旦接触，法向产生压力、切向可能产生摩擦或黏滞，接触面积、接触刚度、接触位置都随载荷变化。典型对象包括：镜座与镜框的贴合面、螺钉头与法兰的承压面、柔性支撑的球头与锥窝、光学元件与压圈的边缘接触。接触问题的核心特征是"状态依赖"——同一结构在不同载荷步下，边界条件集合不同。

几何非线性来自变形量相对于结构尺寸不可忽略。当位移或转角较大时，应变-位移关系不再线性，刚度矩阵随构形变化，载荷方向也可能随变形转动。光机结构中常见于：大口径镜面的重力变形、长悬臂支撑、薄壁镜筒的失稳、柔性铰链的大转角（可达几度到十几度）、以及预紧力作用下的细长杆件。此时必须区分"变形前构形"与"变形后构形"，采用总拉格朗日或更新拉格朗日格式。

材料非线性来自应力-应变关系不再线性。金属在超过屈服后进入塑性，胶粘剂、橡胶垫、聚氨酯等呈现超弹性或黏弹性，某些光学胶还存在固化收缩与温度依赖。光机结构中，胶层连接、减振垫、密封圈、以及局部应力集中区的微塑性，都属于材料非线性范畴。

二、耦合关系：为什么不能分开处理

三类非线性在光机结构中很少单独出现。以典型的反射镜组件为例：镜体通过三点柔性支撑与镜座连接，螺钉施加预紧力。预紧过程中，螺钉头与法兰面从分离到接触，接触面积随预紧力增大而扩展，这是接触非线性；柔性支撑在预紧力下产生较大转角，几何构形改变，这是几何非线性；胶层在局部高压区进入塑性或黏弹性，这是材料非线性。三者同时发生，且相互影响：接触面积变化改变局部刚度，局部刚度改变变形分布，变形分布又反过来改变接触压力分布。

这种耦合意味着不能"先算接触、再算几何、最后算材料"。正确的做法是在同一个增量-迭代框架内同时更新接触状态、构形和材料本构。Newton-Raphson 方法在每个载荷增量步内反复组装切线刚度矩阵，直到位移残差和力残差同时满足收敛容差。接触状态用罚函数或拉格朗日乘子法处理，几何非线性用总拉格朗日格式，材料非线性用径向返回算法更新应力。

三、光机结构为何不能只用线性静力学替代

线性静力学替代非线性仿真，在以下场景会产生显著误差甚至错误结论：

• 预紧力传递失真：线性模型无法模拟螺钉从分离到接触的过程，预紧力被简化为节点力，导致镜座贴合面压力分布与实际不符，进而影响面形预测。
• 柔性铰链转角被低估：线性小变形假设下，柔性铰链的等效刚度偏高，转角偏小，支撑点位移计算错误，最终影响光学元件的位置精度。
• 胶层应力被高估：线弹性模型下胶层边缘应力集中可达实际值的数倍，若据此判断胶层失效，会得出过于保守的结论；反之，若忽略塑性耗能，又可能低估长期蠕变。
• 接触分离无法捕捉：当载荷方向改变或温度变化导致零件分离时，线性模型仍保持连接，刚度矩阵不变，计算结果完全失真。
• 失稳现象被遗漏：薄壁镜筒或长悬臂在压力载荷下可能发生屈曲，线性静力学无法预测分叉点，而非线性屈曲分析可以给出临界载荷。

四、工程判断与建模建议

实际建模时，不必对所有零件都开启非线性。合理的做法是分级处理：对远离接触区、变形小、应力低的零件保持线弹性小变形假设；对接触区、柔性环节、胶层等关键区域开启非线性。这样既保证精度，又控制计算成本。接触对的定义应优先使用面-面接触，法向采用硬接触或罚函数，切向根据实际是否允许滑移选择摩擦模型或无摩擦。对于预紧力，建议用螺栓截面预紧单元或过盈配合模拟，而不是直接施加节点力。

收敛性是非线性接触仿真的核心难点。常见策略包括：将载荷分为 10～20 个增量步，开启自动时间步，设置合理的残差容差（力残差约 0.5%，位移残差约 5%），对接触刚度进行自适应调整。若出现不收敛，先检查接触对是否初始穿透、刚体位移是否约束充分、材料参数是否合理，再考虑调整求解器设置。

总结来说，光机结构的非线性接触仿真不是"更复杂的线性分析"，而是对物理过程的更真实描述。接触非线性处理边界状态变化，几何非线性处理构形变化，材料非线性处理本构变化，三者在增量-迭代框架内耦合求解。线性静力学只在特定简化条件下成立，一旦涉及预紧、过盈、大变形、胶层或分离，就必须切换到非线性分析。理解这一点，是建立可信光机仿真模型的第一步。