Abaqus笔帽接触分析用轴对称简化，计算快还准

3D全模型网格划到80万节点，跑一次过盈配合插拔四个多小时，换轴对称之后节点数掉到1.8万，单核二十分钟出结果，精度对比过接触压力曲线几乎重合。

笔帽这种回转体天生就是轴对称的理想对象，CAE里Part创建时把Solid、Axisymmetric勾上，草图画在XZ平面，X必须是正的，轴线是Y轴不是X轴，这一点记错过两次，负坐标区域直接报The nodes of axisymmetric elements must lie in positive X，检查半天才发现草图画反了。

单元选CAX4R，减缩积分四节点，接触区局部加密到0.1mm量级，远场可以放到0.5mm，沙漏控制用默认enhanced就够，笔帽壁厚1.2mm的话厚度方向铺4层单元起步，3层以下弯曲应力明显偏低，CAX8R二阶单元精度更高但接触收敛慢，过盈配合这种大变形接触还是R单元省心。

接触定义和3D一样，帽内壁当主面笔杆外壁当从面，法向Hard，切向Penalty摩擦系数0.25，过盈量0.08mm用*CONTACT INTERFERENCE在Interaction里给，配合Surface to Surface，轴对称模型的接触搜索是二维的，收敛比三维稳得多，增量步很少被切。

边界条件只能约束U1和U2，U3在轴对称单元里根本不存在，自由度编号是1、2、6，6是绕轴旋转，笔帽拔出力算出来是单位弧度上的值，乘以2π才是整圈的力，这一步忘了换算数值差六倍多，之前对着结果纳闷了半天量级不对。

位移加载走Amplitude线性，分析步Static General开几何非线性Nlgeom，初始增量0.01，最小1e-6，拔出行程按实际插拔深度给，比如3mm，历史输出勾CAREA和CFS接触面积和接触力，反力从笔杆约束点的RF2读。

材料是PP，弹性模量1.2GPa泊松比0.4，塑性段给真实应力应变，工程应力应变换算记得乘(1+e)和ln(1+e)，不换算的话大变形段刚度偏低，接触压力峰值差15%以上，HyperWorks做前处理导轴对称模型路径麻烦，还是Abaqus原生CAE画二维草图来得直接。

结果后处理View→ODB Display Options里把镜像打开，Deform Mirror可以绕轴展开看整圈效果，云图是二维的一条带，配合镜像至少能看出周向分布，应力集中位置在帽口倒角根部，von Mises峰值18MPa，低于PP屈服26MPa，插拔寿命按磨损算不按强度算。

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删