单元的生和死

1、单元的生和死何为单元的生和死？如果模型中加入（或删除）材料，模型中相应的单元就存在“（或消亡）。单元生死选项就用于在这种情况下杀死或重新激活选择的单元。（可用的单元类型在表6-1中列出。）本选项主要用于钻孔（如开矿和挖通道等），建筑物施工过程（如桥的建筑过程），顺序组装（如分层的计算机芯片组装）和另外一些用户可以根据单元位置来方便的激活和不激活它们的一些应用中。单元生死功能只适用于ANSYS/Multiphysics,ANSYS/Mechanical和ANSYS/Structural产品。Table6-1ElementswithbirthanddeathcapabilityLINK1SURF1

2、9SHELL41SOLID64PLANE831SHELL143PLANE2PIPE20PLANE421SOLID65SOLID87:SURF151BEAM3MASS21SHELL43PLANE67SOLID90SURF152BEAM4SURF22BEAM44LINK68SOLID92SURF153SOLID5BEAM23SOLID451SOLID69SHELL93SURF154LINK8BEAM24PLANE53SOLID70SOLID95SHELL157LINK10PLANE25BEAM541MASS71SOLID96丁ARGE169LINK11MATRIX27PLANE55SOLID72S

3、OLID97TARGE170PLANE13LINK31SHELL571SOLID73SOLID98CONTA171COMBIN14LINK32PIPE59PLANE75SHELL99CONTA172PIPE16LINK33PIPE60PLANE77PLANE121CONTA173PIPE17LINK34SOLID62PLANE78SOLID122CONTA174PIPE18PLANE35SHELL631PLANE82SOLID123在一些情况下，单元的生死状态可以根据ANSYS勺计算数值决定，如温度,应力，应变等。可以用ETABLE命令（MainMenu>GeneralPostproc&

4、gt;ElementTable>DefineTable）和ESEL命令（UtilityMenu>Select>Entities）来确定选择的单元的相关数据，也可以改变单元的状态（溶和，固结，俘获等）。本过程对于由相变引起的模型效应（如焊接过程中原不生效的熔融材料变为生效的模型体的一部分），失效扩展和另外一些分析过程中的单元变化是有效的。单元生死是如何工作的？要激活单元死”的效果，ANSY邻序并不是将杀死”的单元从模型中删除，而是将其刚度（或传导，或其他分析特性）矩阵乘以一个很小的因子ESTIFo因子缺省值为1.0E-6,可以赋为其他数值（详见施加载荷并求解”一章）。死单元的单

5、元载荷将为0,从而不对载荷向量生效（但仍然在单元载荷的列表中出现）。同样，死单元的质量，阻尼，比热和其他类似效果也设为0值。死单元的质量和能量将不包括在模型求解结果中。单元的应变在杀死”的同时也将设为00与上面的过程相似，如果单元出生”，并不是将其加到模型中，而是重新激活它们。用户必须在PREP7中生成所有单元，包括后面要被激活的单元。在求解器中不能生成新的单元。要加入”一个单元，先杀死它，然后在合适的载荷步中重新激活它。当一个单元被重新激活时，其刚度，质量，单元载荷等将恢复其原始的数值。重新激活的单元没有应变记录（也无热量存储等）。但是，初应变以实参形式输入（如LINK1单元）的不为单元生死

6、选项所影响。而且，除非是打开了大6-1变形选项NLGEOMQNL些单元类型将以它们以前的几何特性恢复（大变形效果有时用来得到合理的结果）。单元在被激活后第一个求解过程中同样可以有热应变（等于a*（T-TREF）,如果其承受热量体载荷。如何使用单元生死特性可以在大多数静态和非线性瞬态分析中使用单元生死，其基本过程与相应的分析过程是一致的。对于其他分析来说，这一过程主要包括以下三步：建模施加载荷并求解查看结果修改基本分析步骤如下以包括单元生死特征：建模：在PREP汁，生成所有单元，包括那些只有在以后载荷步中才激活的单元。在PREP为卜不能生成新的单元。施加载荷并求解：在SOLUTIO版完成以下操作

7、：定义第一个载荷步：在第一个载荷步中，用户必须选择分析类型和所有的分析选项。用下列方法指定分析类型：Command:ANTYPEGUI:MainMenu>Solution>-AnalysisType-NewAnalysis在结构分析中，大变形效果应打开。用下列命令设置该选项：Command:NLGEOM,ONGUI:MainMenu>Solution>AnalysisOptions对于所有单元生死应用，在第一个载荷步中应设置牛顿-拉夫森选项，因为程序不能预知EKILL命令出现在后面的载荷步中。用下列命令完成该操作：Command:NROPTGUI:MainMenu>

8、;Solution>AnalysisOptions杀死EKILL所有要加入到后续载荷步中的单元，用下列命令：Command:EKILLGUI:MainMenu>Solution>-LoadStepOpts-Other>KillElements单元在载荷步的第一个子步被杀死（或激活），然后在整个载荷步中保持该状态。要注意保证使用缺省的矩阵缩减因子不会引起一些问题。有些情况下要考虑用严格的缩减因子。用下列方法指定缩减因子数值：Command:ESTIFGUI:MainMenu>Solution>Other>StiffnessMult不与任何激活的单元相连的

9、结点将漂移”，或具有浮动的自由度数值。在一些情况下，用户可能想约束不被激活的自由度D,CP等以减少要求解的方程的数目，并防止出现位置错误。约束非激活自由度，在重新激活的单元要有特定的（或温度等）时很有影响，因为在重新激活单元时要删除这些人为的约束。同时要删除非激活自由度的结点载荷（也就是不与任意激活的单元相连的结点）。同样，用户必须在重新激活在自由度上施加新的结点载荷。6-2下面是第一个载荷步中命令输入示例：!第一个载荷步TIME,.!设定时间值（静力分析选项）NLGEOM,ON!打开大位移效果NROPT,FULL!设定牛顿拉夫森选项ESTIF,.!设定非缺省缩减因子（可选）ESEL,.!选择

10、在本载荷步中将不激活的单元EKILL,.!不激活选择的单元ESEL,S,LIVE!选择所有活动单元NSLE,S!选择所有活动结点NSEL,INVE!选择所有非活动结点（不与活动单元相连的结点）D,ALL,ALL,0!约束所有不活动的结点自由度（可选）NSEL,ALL!选择所有结点ESEL,ALL!选择所有单元D,.!施加合适的约束F,.!施加合适的活动结点自由度载荷SF,.!施加合适的单元载荷BF,.!施加合适的体载荷SAVESOLVE请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D,F,SFBF命令得到更详细的解释。后继载荷步在后继载荷步中，

11、用户可以随意杀死或重新激活单元。象上面提到的，要正确的施加和删除约束和结点载荷。用下列命令杀死单元：Command:EKILLGUI:MainMenu>Solution>-LoadStepOpts-Other>KillElements用下列命令重新激活单元：Command:EALIVEGUI:MainMenu>Solution>-LoadStepOpts-Other>ActivateElem!第二个（或后继）载荷步：TIME,.ESEL,.EKILL,.!杀死选择的单元ESEL,.EALIVE,.!重新激活选择的单元.FDELE,.!删除不活动自由度的结点载

12、荷D,.!约束不活动自由度.F,.!在活动自由度上施加合适的结点载荷6-3DDELE,.!删除重新激活的自由度上的约束SAVESOLVE请参阅TIME,ESEL,EKILL,EALIVE,FDELE,D,F和DDEL瑜令得至U更详细的解释。查看结果对于大多数部分来说，用户在对包含不激活或重新激活的单元操作时应按照标准的过程来做。但是必须清楚的是，杀死”的单元仍在模型中，尽管对刚度（传导）矩阵的贡献可以忽略。因此，它们将包括在单元显示，输出列表等操作中。例如，不激活的单元在结点结果平均（PLNSOL#令或MainMenu>GeneralPostproc>PlotResults>

13、NodalSolu）时将污染”结果。整个不激活单元的输出应当被忽略，因为很多项带来的效果都很小。建议在单元显示和其他后处理操作前用选择功能将不激活的单元选出选择集。使用ANSYS吉果控制单元生死在许多时候，用户并不清楚的知道杀死和重新激活单元的确切位置。例如，用户要在热分析中杀死”熔融的单元（在模型中移去溶化的材料），事先不会知道这些单元的位置；用户必须根据ANSYS+算出的温度确定这些单元。当决定杀死或重新激活单元依靠ANSYS+算结果时（如温度，应力，应变等），用户可以使用命令识别并选择关键单元。用下列方法识别关键单元：Command:ETABLEGUI:MainMenu>Gener

14、alPostproc>ElementTable>DefineTable用下列方法选择关键单元：Command:ESELGUI:UtilityMenu>Select>Entities然后用户可以杀死或重新激活选择的单元。（也可以用ANSYSAPD语言编写宏以完成这些操作。）用下列方法杀死选择的单元：Command:EKILL,ALLGUI:MainMenu>Solution>-LoadStepOpts-Other>KillElements用下列方法重新激活选择的单元：Command:EALIVE,ALLGUI:MainMenu>Solution&g

15、t;-LoadStepOpts-Other>ActivateElem下面的例子是杀死总应变超过许用值的单元：/SOLU!进入求解器.!标准的求解过程SOLVEFINISH!/POST1!进入POST1SET,.6-4ETABLE,STRAIN,EPTO,EQV冬总应变存入ETABLEESEL,S,ETAB,STRAIN,0.20!选择所有总应变大于或等于0.20的单元FINISH!/SOLU!重新进入求解器ANTYPE,RESTEKILL,ALL!杀死选择（超过允许值）的单元ESEL,ALL!读入所有单元.!继续求解请参阅ETABLE,ESEL,ANTYPEEKILL命令得至IJ更详细的

16、解释。进一步的说明不活动的自由度上不能施加约束方程CE,CEINTF等。（不活动的自由度当结点不与活动的单元相连时出现）。可以通过先杀死然后重新激活单元的方法做应力松弛（如退火）操作。在非线性分析中，注意不要因为杀死或重新激活单元引起奇异性（如结构分析中的尖角）或刚度突变。这将使得收敛困难。在有单元生死的分析中打开FULL牛顿-拉夫森方法的适应下降选项将得到好的结果。用下列方法：Command:NROPT,FULL,ONGUI:MainMenu>Solution>AnalysisOptions可以通过一个参数值来指示单元生死状态*GET,Par,ELEM,n,ATTR,LIVE（U

17、tilityMenu>Parameters>GetScalarData）。该参数可以用于APDLL辑分支（*IF等），或其他要控制单元生死的应用场合中。用户可能想通过改变材料特性来杀死或重新激活单元MPCHG（MainMenu>Preprocessor>MaterialProps>ChangeMatNum）但是，在这个过程中要特别小心。软件保护系统和限制使得杀死”的单元在求解器中改变材料特性时将不生效。（单元集中载荷不能自动删除；应变，质量，比热等也不能删除。）不当的使用MPCHG令将带来许多问题。例如，如果将单元的刚度缩减到近于0,而保留其质量，在有加速度和惯性

18、载荷的问题中将产生奇异性。一个MPCHGJ应用是在建立模型时涉及出生”单元的应变历程的情况下。使用MPCHGT以得到单元在变形的结点构造中的初始应变。在单元生死中不能用多载荷步求解LSWRITE,因为不激活或重新激活的单元状态将不写入载荷步文件中。有多个载荷步的生死单元分析应该用一系列的SOLV命令（MainMenu>Solution>CurrentLS）来做。单元生死应用实例（命令行格式）问题描述等截面杆两端固定，承受均匀的温度载荷时将其中间1/3段移去。过程是将其应变自由化并移去均匀温度。分析其热应力和应变情况。材料特性和几何模型参数见下图。6-5MaterialpropertiesE=30x1gpsja=0.0D005in/in-TGeometricPropertiesL=10in