单元生死技术.docVIP

单元生死是如何工作的？ 要激活“单元死”的效果，ANSYS程序并不是将“杀死”的单元从模型中删除，而是将其刚度（或传导，或其他分析特性）矩阵乘以一个很小的因子[ESTIF]。因子缺省值为1.0E-6，可以赋为其他数值。死单元的单元载荷将为0，从而不对载荷向量生效（但仍然在单元载荷的列表中出现）。同样，死单元的质量，阻尼，比热和其他类似效果也设为0值。 死单元的质量和能量将不包括在模型求解结果中。单元的应变在“杀死”的同时也将设为0。 与上面的过程相似，如果单元“出生”，并不是将其加到模型中，而是重新激活它们。用户必须在PREP7 中生成所有单元，包括后面要被激活的单元。在求解器中不能生成新的单元。要“加入”一个单元，先杀死它，然后在合适的载荷步中重新激活它。 当一个单元被重新激活时，其刚度，质量，单元载荷等将恢复其原始的数值。重新激活的单元没有应变记录（也无热量存储等）。但是，初应变以实参形式输入（如LINK1 单元）的不为单元生死选项所影响。而且，除非是打开了大变形选项[NLGEOM,ON],一些单元类型将以它们以前的几何特性恢复（大变形效果有时用来得到合理的结果）。单元在被激活后第一个求解过程中同样可以有热应变（等于a*(T-TREF)）,如果其承受热量体载荷。 实例1 实例2 /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 ! ! 假设各材料都是双线性随动硬化弹塑性本构关系 ! 下面是各材料在各温度点下的屈服应力和屈服后的弹性模量 tb,bkin,1,5 tbtemp,20,1 tbdata,1,1200e6,0.193e11 tbtemp,500,2 tbdata,1, 933e6,0.150e11 tbtemp,1000,3 tbdata,1, 435e6,0.070e11 tbtemp,1500,4 tbdata,1, 70e6,0.010e11 tbtemp,2000,5 tbdata,1, 7e6,0.001e11 ! tb,bkin,2,5 tbtemp,20,1 tbdata,1,800e6,0.102e11 tbtemp,500,2 tbdata,1,400e6,0.050e11 tbtemp,1000,3 tbdata,1, 70e6,0.008e11 tbtemp,1500,4 tbdata,1, 1e6,0.0001e11 tbtemp,2000,5 tbdata,1,0.1e6,0.00001e11 ! tb,bkin,3,5 tbtemp,20,1 tbdata,1,900e6,0.117e11 tbtemp,500,2 tbdata,1,700e6,0.090e11 tbtemp,1000,3 tbdata,1,230e6,0.030e11 tbtemp,1500,4 tbdata,1, 40e6,0.005e11 tbtemp,2000,5 tbdata,1, 4e6,0.0005e11 ! ! 材料密度(假设为常值) mp,dens,1,8030 mp,dens,2,4850 mp,dens,3,8900 ! 热膨胀系数(假设为常值) mp,alpx,1,1.78e-5 mp,alpx,2,9.36e-6 mp,alpx,3,1.66e-5 ! 泊松比(假设为常值) mp,nuxy,1,0.29 mp,nuxy,2,0.30 mp,nuxy,3,0.30 ! 热传导系数(假设为常值) mp,kxx,1,16.3 mp,kxx,2,7.44 mp,kxx,3,393