Abaqus 残余应力仿真全流程 :冷扩衬套如何把耳片疲劳寿命提升约 5 倍
冷扩衬套如何把耳片疲劳寿命提升约 5 倍：一套 Abaqus 残余应力仿真全流程
起落架主接头的耳片，几乎每一次地面慢速转弯都在经受弯曲叠加拉压的循环载荷，疲劳是它最头疼的设计瓶颈。空客军用飞机（AIM）想回答一个问题：把衬套改成冷扩（Cold Expansion）安装，到底能让寿命提升多少？能不能在做试验之前就算出来？
先看结论
对象：7010-T74 铝合金锻造耳片 + 青铜衬套，对比间隙配合、低温收缩配合、冷扩（Force-Mate）三种安装方式。
变幅谱（代表实际使用）下，冷扩相对其他方案的寿命提升系数（LIF）约为 5。
残余应力：2D 平面应变模型可作初步估算，3D 模型能给出沿厚度的分布。
寿命预测：应力-寿命法整体优于应变-寿命法；冷扩构型的预测可用于初步评估，间隙配合预测偏不保守。
一、问题从哪来
起落架结构按安全寿命（safe-life）设计，依据 CS 25.571 通过“分析 + 试验”取证。试验成本高，还会拖慢研制进度，所以工程上很希望在试验前就有一个靠谱的分析预判。
耳片的危险位置在内孔与外表面的交界处，方向垂直于载荷线。团队评估过几种改进路线：
加大耳片尺寸：空间受限，还要满足收放等系统要求，收益很有限。
换材料：钛合金成本明显上升；换钢则整体增重太多。
高强钢空心销：销壁可以做薄，但变形更大，对降低应力集中帮助小，增重不划算。
改衬套安装方式：采用冷扩高过盈衬套，被证明最有效，但业内没有公认的疲劳收益分析方法。
最终设计是局部改型 + 销加强 + 冷扩衬套的组合，本文研究的是其中最有效的一项：FTI 公司的 Force-Mate 冷扩衬套。
二、研究思路与试验
方法分三步：① 评估残余应力 → ② 评估应力集中系数 → ③ 计算疲劳寿命（应力-寿命、应变-寿命两种路线），并用试验验证。
试验件与加载
试验件是实际几何的简化版，材料为 7050-T7451 板材（原始厚度 75–80 mm，加工至 65 mm），内孔半径 105 mm，外半径 202 mm，表面做铬酸阳极化加底漆。在 500 吨（5000 kN）万能试验机上只加拉力以避免屈曲，并让载荷相对试件中心线偏置，用来模拟实际零件的弯曲应力。三种衬套构型均进行了常幅与变幅谱试验。
冷扩衬套怎么装
冷扩衬套的外径略小于耳片内孔，先装入孔内，再在常温下用一根外径大于衬套内径的芯棒拉过去，把衬套撑开，与耳片形成高过盈配合，孔周产生塑性变形。
三、残余应力：先用解析解给 FEM 把关
塑性应变：厚壁筒内压扩张
作者先推导了解析解：轴对称、平面应变、理想弹塑性（无硬化），塑性应变由 Prandtl-Reuss 方程和 von Mises 屈服准则得到。然后用 1/4 圆 FEM 做对比，在精度与耗时之间选定网格：二阶四边形单元，每 2.5° 一个单元（整圈 144 个）。
接触：用过盈配合验证
同样推导了过盈配合（轴对称、平面应变、弹性）的解析解，与 FEM 圆环模型对比，验证了接触设置。
解析解为什么不够用
完整安装过程的解析解已经做出来，但无法直接用于后续疲劳计算：一是塑性硬化会让应力分布由尖峰变成圆滑曲线（图 6）；二是真实形状不是轴对称的；三是必须能表达后续加载过程。所以必须上 FEM。
四、2D 有限元模型
虽然问题本质是三维的，但 2D 模型能很快做参数敏感性分析（边界条件、材料行为、加载方式），再用更耗时的 3D 去校验。
收缩配合：把衬套外径设得略大于耳片内孔，利用 Abaqus 的接触面调整功能处理过盈。
冷扩分两步：扩张——在衬套内径施加均匀径向位移，等效芯棒穿过厚度；松弛——去掉该边界条件，模拟芯棒移除后的卸载。
接触：耳片-衬套、衬套-销均为面对面、无摩擦、小滑移，衬套面作从面；初始过闭合通过调整消除。
网格：对称性只建 1/4，平面应变二阶缩减积分单元 CPE8R，半圆孔周 72 个单元。
材料方面，铝（耳片）和青铜（衬套）按弹塑性处理，应力-应变曲线采用 Ramberg-Osgood 公式，钢销（4340）只取弹性。硬化模型虽然比较过随动硬化，最终选用了各向同性硬化。
五、3D 有限元模型：看清沿厚度的分布
3D 模型把芯棒沿轴向真正“拉过去”，得到残余应力沿厚度的分布，同时能模拟销加载，评估弯曲效应。代价是内存和 CPU 时间都很高，因此对试件几何做了简化，并只在孔周附近细化网格。
芯棒为解析刚体，无需划分网格或定义材料；材料定义沿用 2D 模型。
新增分析步：步 1、2 先后建立衬套法兰面与圆柱面和耳片的接触；步 3、4 为芯棒拉过衬套孔完成扩张，期间耳片下表面支撑在拉拔机上，以边界条件表示（专门验证过，引入拉拔机刚度并不会改善结果）；之后释放附加边界条件。
接触拆成圆柱体与法兰两部分，并引入 0.1 的摩擦系数帮助收敛；销与芯棒只与衬套内表面接触。
利用对称性建半个耳片，单元为 8 节点缩减积分砖单元 C3D8R。
六、仿真能信吗：孔钻法实测对比
团队评估了中子衍射和钻孔应变释放法，前者技术上更优，但受项目限制最终选了后者：在表面贴多个应变片，在附近钻小孔，根据应力释放量反算残余应力（按 ASTM 技术说明 TN-503-6 计算）。
预算和时间有限，没能在试验前测量，只能测已经做过试验的试件（选取两个耳片中未失效的那个）。这意味着试验载荷本身可能已经改变了残余应力场，对比自然更困难。
从图 13 可以读到：2D 预测与实测点吻合较好，孔边为较大的压应力；3D 不同路径（Path-2、Path-7）给出的孔边压应力差异明显（读图约在 −270 至 −440 MPa 之间），说明残余应力沿厚度并不均匀。作者也坦率指出这种吻合“可能是巧合”，并给出了三点保留：
七、从残余应力到疲劳寿命
疲劳评估沿用 AIM 的常规流程：应力集中系数 Kt + 合适的 S-N 数据 + Miner 线性累积损伤。由于冷扩改变了危险点的局部应力，这里直接使用局部应力，而不是参考应力。Kt 定义为：
损伤累积采用 Miner 准则（通用形式）：D = Σ ni / Ni。
三种构型的局部应力有何不同
加载过程仿真里的一个小技巧
实际装配中，冷扩松弛后衬套内径仍小于销径，最后再铰孔得到与销的间隙配合。FEM 里不模拟铰孔，而是让销外径等于衬套内径。对冷扩构型，松弛后衬套内径反而比销略大，销与衬套之间出现缝隙，因此要先加一个很小的过渡载荷步让销贴合衬套，并启用 Abaqus 的稳定化（stabilization）：相当于在销与衬套间加一层粘性材料，随增量加载其刚度递减，到最后一个增量帧刚度为零。之后再对销的剪切面施加正式载荷。
两条计算路线
应力-寿命法：先把载荷谱（吨）转成危险点缺口应力（初始残余应力来自 FEM，曲线斜率体现 Kt），再用多 R 比的 S-N 曲线配合 Miner 准则累积损伤。
应变-寿命法：间隙与收缩配合同上；冷扩则通过“先施加一次过载”来模拟，使其产生的环向残余应力等于 FEM 结果，再用 Neuber 规则计算缺口应力应变，并用 ε-N 曲线与 Miner 准则累积损伤，材料数据来自 MMPDS-04。
八、预测 vs 试验
实际耳片按安全寿命设计，关注的是裂纹萌生寿命；多处断口分析显示，萌生寿命约占总寿命的三分之二。冷扩相对间隙或收缩配合的寿命提升系数与载荷水平有关；考虑各试件实测弯曲比（载荷偏心造成）后，变幅谱下 LIF 约为 5。
应力-寿命法整体更接近试验，变幅谱的相关性好于常幅谱。
间隙配合差距更大，预测偏不保守，推测原因是分析没有计入微动磨损；冷扩因过盈高，微动影响小得多。
冷扩构型的应力-寿命预测合理，可用于初步评估。
九、结论与展望
试验证明大尺寸（内径 105 mm）铝合金耳片采用冷扩衬套有明显疲劳收益，变幅谱下 LIF 约 5。
2D 模型适合做初步估算；3D 模型才能反映厚度方向的差异，但与实测相关性仍待深入。
基于该研究，AIM 正在建立冷扩类零件的疲劳预测方法；后续还需要更明确的裂纹萌生定义，以及更多残余应力测量与试验相关性工作。
