ABAQUS基础第56讲：热传导、热应力与膜接触一块讲

打开sta文件全是U，Too many attempts made to have numerically severe contact，步才推进到1.2e-3就被杀了，这种开局热成形类模型十次有八次。

先把分析步类型定死，全耦合走Coupled temp-displacement，CAE里Step→Create→Coupled temp-displacement，纯温度场走Heat transfer，热应力想省事就顺序耦合，先跑热传导，再在同一个model里建Static步，Predefined Field里Create→Temperature→Read，从上一个job的odb把节点温度场读进来，材料里膨胀系数和参考温度必须有，参考温度在Initial conditions→Temperature给，默认零的话升温过程整个件从0度往上胀，应变直接爆。

单元选择，温度位移耦合C3D8RT，R是减缩积分，纯热传导C3D8T够用，薄膜用M3D4R或SFM3D4，厚度小于特征尺寸十分之一的时候拿S4R壳顶替膜单元更稳，壳有弯曲刚度，膜一点没有，受压就皱成一团。

接触面的热参数全塞在Interaction Property里，三样东西，Thermal conductance给间隙导热表，clearance对conductance两三个点线性外推，Gap radiation给发射率，view factor默认1，对流不在这个property里，用Interaction→Surface Film Condition直接加载到暴露面上，film coefficient加sink temperature，系数做成温度相关的表格。Marc里这三项一个对话框全给齐，Abaqus散在三个入口，第一次配的人基本都要翻半天文档。

去年11月跑真空成形覆膜那个算例，那几天办公室暖气没来，穿着羽绒服挂机，旁边同事在嗦螺蛳粉，声音特别响，模型是半球冲头压0.05mm的PET膜，第一次跑膜整个起皱，job推到30%报Excessive element distortion，把网格从2mm加密到1mm再跑还是崩，才反应过来是膜单元法向翻了，接触搜索找不到主面，Mesh模块Edit→Flip Normals翻过来，再在接触属性Normal behavior里勾Automatic stabilization，耗散能分数2e-4，褶皱才压得住，整个case跑了四十多分钟，盯着进度条的心情跟看倒计时一样。

膜和模具的主从，刚性模具当master，膜当slave，穿透检查在slave侧，离散方式选surface to surface配Finite sliding，比node to face稳一个量级，pressure-overclosure用Hard，法向接触刚度scale factor默认10会把增量步压得很碎，降到1步能顺不少，切向Penalty，摩擦系数0.1到0.3看材料对。

增量步和输出，耦合分析温度自由度收敛判据是温度增量幅值的0.5%，Step→Other里能看到，场输出请求记得勾NT11、HFL、TEMPC，默认输出只有位移和应力，后处理画不了温度云图，顺序耦合那步也读不到温度场，Read直接报No temperature results found。

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