ANSYS工程应用教程.docVIP

《ANSYS工程应用教程——机械篇》 /cgi-bin/ut/topic_show.cgi?id=8687h=1#45806 讲了单元生死应用。 ??单元生死应用实例（命令行格式） 问题描述 等截面杆两端固定，承受均匀的温度载荷时将其中间1/3段移去。过程是将其应变自由化并移去均匀温度。分析其热应力和应变情况。 ????/PREP7? ???/TITLE, ELEMENT BIRTH/DEATH IN A FIXED BAR WITH THERMAL LOADING ??ET,1,LINK1? ?? ?? ? !二维杆单元 ??MP,EX,1,30E6? ?? ???!材料特性 ??MP,ALPX,1,.00005 ??MP,EX,2,30E6 ??MP,ALPX,2,.00005? ? !重新‘出生’单元的特性 ??MP,REFT,2,100? ?? ? !单元出生的参考温度 ??R,1,1.0 ??N,1????N,4,10? ???FILL ??E,1,2? ?????EGEN,3,1,-1? ?? ?? ?!生成三个单元 ??FINISH? ??? ??/SOLU? ? ??ANTYPE,STATIC ??D,1,ALL,,,4,3? ?? ? !固定杆的两端 ??TREF,0? ?? ?? ?? ???!参考温度0度 ??TUNIF,100? ?? ?? ???!均匀温度载荷 ??NROPT,FULL ??OUTPR,BASIC,ALL ??SOLVE? ? ??EKILL,2? ?? ?? ?? ? !‘杀死’中间的单元 ??SOLVE ??EALIVE,2? ?? ?? ?? ?!重新激活中间单元 ??MPCHG,2,2? ?? ?? ???!将材料特性改为2以生成自由应变单元 ??SOLVE ??TUNIF,0? ?? ?? ?? ? !删除温度载荷 ??SOLVE材料特性和几何模型参数单元生死是如何工作的？ 要激活“单元死”的效果，ANSYS程序并不是将“杀死”的单元从模型中删除，而是将其刚度（或传导，或其他分析特性）矩阵乘以一个很小的因子[ESTIF]。因子缺省值为1.0E-6，可以赋为其他数值。死单元的单元载荷将为0，从而不对载荷向量生效（但仍然在单元载荷的列表中出现）。同样，死单元的质量，阻尼，比热和其他类似效果也设为0值。 死单元的质量和能量将不包括在模型求解结果中。单元的应变在“杀死”的同时也将设为0。 与上面的过程相似，如果单元“出生”，并不是将其加到模型中，而是重新激活它们。用户必须在PREP7 中生成所有单元，包括后面要被激活的单元。在求解器中不能生成新的单元。要“加入”一个单元，先杀死它，然后在合适的载荷步中重新激活它。 当一个单元被重新激活时，其刚度，质量，单元载荷等将恢复其原始的数值。重新激活的单元没有应变记录（也无热量存储等）。但是，初应变以实参形式输入（如LINK1 单元）的不为单元生死选项所影响。而且，除非是打开了大变形选项[NLGEOM,ON],一些单元类型将以它们以前的几何特性恢复（大变形效果有时用来得到合理的结果）。单元在被激活后第一个求解过程中同样可以有热应变（等于a*(T-TREF)）,如果其承受热量体载荷。上面的例子就是ANSYS帮助里提供的例子吧？另外有人做焊接过程利用“单元生死”的吗？我也分享一个焊接用生死单元的例子吧： !下面的命令流进行的是一个简单的二维焊接分析, 利用ANSYS单元生死和热-结构耦合分析功能进 !行焊接过程仿真, 计算焊接过程中的温度分布和应力分布以及冷却后的焊缝残余应力。????/title,Weld Analysis by Element Birth and Death??/prep7??/unit,si ! 采用国际单位制??!??et,1,13,4 ! 13号二维耦合单元, 同时具有温度和位移自由度??et,2,13,4??!??! 1号材料是钢??! 2号材料是铝??! 3号材料是铜??! 铝是本次分析中的焊料, 它将钢结构部分和铜结构部分焊接起来??! 下面是在几个温度点下, 各材料的弹性模量??mptemp,1,20,500,1000,1500,2000??mpdata,ex,1,1,1.93e11,1.50e11,0.70e11,0.10e11,0.01e11??mpdata,ex,2,1,1.02e11,0.50e11,0.08e11,0.001e11,0.0001e11??mpdata,ex,3,1,1.17e11,0.90e11,0.30e11,0.05e11,0.005e11??!??! 假设各材料