solution是什么意思

传热分析中的生死单元法是什么含义?

众所周知,在固体力学里面,杀死单元的意思是将其单元刚度矩阵变为接近于零

其物理意义为单元失去承载能力.

那么在传热分析里面杀死单元又代表什么呢,肯定也将单元刚度矩阵变为接近于零,则其物理意义是什么呢.

热分析里杀死单元也是使其刚度（传热和其它分析特性）矩阵近似为零

从而单元的载荷将为零，使单元不对载荷向量生效

死单元的质量、比热、导热等类似效果也设为零值，死单元的质量和能量就不包括在模型求解结果中

我看可否这样理解：物理意义就是单元失去传热以及热容能力

失去传热能力，是不是绝缘了呢？

何为单元的生和死？

如果模型中加入（或删除）材料，模型中相应的单元就“存在”（或消亡）。单元生死选项就用于在这种情况下杀死或

重新激活选择的单元。（可用的单元类型在表6－1中列出。）本选项主要用于钻孔（如开矿和挖通道等），建筑物施

工过程（如桥的建筑过程），顺序组装（如分层的计算机芯片组装）和另外一些用户可以根据单元位置来方便的激活

和不激活它们的一些应用中。单元生死功能只适用于ANSYS/Multiphysics，ANSYS/Mechanical和

ANSYS/Structural产品。

Table6-1Elementswithbirthanddeathcapability

LINK1SURF19SHELL41SOLID64PLANE83SHELL143

PLANE2PIPE20PLANE42SOLID65SOLID87SURF151

BEAM3MASS21SHELL43PLANE67SOLID90SURF152

BEAM4SURF22BEAM44LINK68SOLID92SURF153

SOLID5BEAM23SOLID45SOLID69SHELL93SURF154

LINK8BEAM24PLANE53SOLID70SOLID95SHELL157

LINK10PLANE25BEAM54MASS71SOLID96TARGE169

LINK11MATRIX27PLANE55SOLID72SOLID97TARGE170

PLANE13LINK31SHELL57SOLID73SOLID98CONTA171

COMBIN14LINK32PIPE59PLANE75SHELL99CONTA172

PIPE16LINK33PIPE60PLANE77PLANE121CONTA173

PIPE17LINK34SOLID62PLANE78SOLID122CONTA174

PIPE18PLANE35SHELL63PLANE82SOLID123

在一些情况下，单元的生死状态可以根据ANSYS的计算数值决定，如温度，应力，应变等。可以用ETABLE命令（Main

Menu>GeneralPostproc>ElementTable>DefineTable）和ESEL命令（UtilityMenu>Select>Entities）来

确定选择的单元的相关数据，也可以改变单元的状态（溶和，固结，俘获等）。本过程对于由相变引起的模型效应（如

焊接过程中原不生效的熔融材料变为生效的模型体的一部分），失效扩展和另外一些分析过程中的单元变化是有效的。

单元生死是如何工作的？

要激活“单元死”的效果，ANSYS程序并不是将“杀死”的单元从模型中删除，而是将其刚度（或传导，或其他分析

特性）矩阵乘以一个很小的因子[ESTIF]。因子缺省值为1.0E-6，可以赋为其他数值（详见“施加载荷并求解”一

章）。死单元的单元载荷将为0，从而不对载荷向量生效（但仍然在单元载荷的列表中出现）。同样，死单元的质量，

阻尼，比热和其他类似效果也设为0值。死单元的质量和能量将不包括在模型求解结果中。单元的应变在“杀死”的

同时也将设为0。

与上面的过程相似，如果单元“出生”，并不是将其加到模型中，而是重新激活它们。用户必须在PREP7中生成所有

单元，包括后面要被激活的单元。在求解器中不能生成新的单元。要“加入”一个单元，先杀死它，然后在合适的载荷

步中重新激活它。

当一个单元被重新激活时，其刚度，质量，单元载荷等将恢复其原始的数值。重新激活的单元没有应变记录（也无热量

存储等）。但是，初应变以实参形式输入（如LINK1单元）的不为单元生死选项所影响。而且，除非是打开了大变形

选项[NLGEOM,ON],一些单元类型将以它们以前的几何特性恢复（大变形效果有时用来得到合理的结果）。单元在被

激活后第一个求解过程中同样可以有热应变（等于a*(T-TREF)),如果其承受热量体载荷。

如何使用单元生死特性

可以在大多数静态和非线性瞬态分析中使用单元生死，其基本过程与相应的分析过程是一致的。对于其他分析来说，这

一过程主要包括以下三步：

建模

施加载荷并求解

查看结果

修改基本分析步骤如下以包括单元生死特征：

建模：

在PREP7中，生成所有单元，包括那些只有在以后载荷步中才激活的单元。在PREP7外不能生成新的单元。

施加载荷并求解：

在SOLUTION中完成以下操作：

定义第一个载荷步：

在第一个载荷步中，用户必须选择分析类型和所有的分析选项。用下列方法指定分析类型：

Command:ANTYPE

GUI:MainMenu>Solution>-AnalysisType-NewAnalysis

在结构分析中，大变形效果应打开。用下列命令设置该选项：

Command:NLGEOM,ON

GUI:MainMenu>Solution>AnalysisOptions

对于所有单元生死应用，在第一个载荷步中应设置牛顿－拉夫森选项，因为程序不能预知EKILL命令出现在后面的载荷

步中。用下列命令完成该操作：

Command:NROPT

GUI:MainMenu>Solution>AnalysisOptions

杀死[EKILL]所有要加入到后续载荷步中的单元，用下列命令：

Command:EKILL

GUI:MainMenu>Solution>-LoadStepOpts-Other>KillElements

单元在载荷步的第一个子步被杀死（或激活），然后在整个载荷步中保持该状态。要注意保证使用缺省的矩阵缩减因子

不会引起一些问题。有些情况下要考虑用严格的缩减因子。用下列方法指定缩减因子数值：

Command:ESTIF

GUI:MainMenu>Solution>Other>StiffnessMult

不与任何激活的单元相连的结点将“漂移”，或具有浮动的自由度数值。在一些情况下，用户可能想约束不被激活的自

由度[D,CP等]以减少要求解的方程的数目，并防止出现位置错误。约束非激活自由度，在重新激活的单元要有特定的

（或温度等）时很有影响，因为在重新激活单元时要删除这些人为的约束。同时要删除非激活自由度的结点载荷（也

就是不与任意激活的单元相连的结点）。同样，用户必须在重新激活在自由度上施加新的结点载荷。

下面是第一个载荷步中命令输入示例：

!第一个载荷步

TIME,...!设定时间值（静力分析选项）

NLGEOM,ON!打开大位移效果

NROPT,FULL!设定牛顿－拉夫森选项

ESTIF,...!设定非缺省缩减因子（可选）

ESEL,...!选择在本载荷步中将不激活的单元

EKILL,...!不激活选择的单元

ESEL,S,LIVE!选择所有活动单元

NSLE,S!选择所有活动结点

NSEL,INVE!选择所有非活动结点（不与活动单

元相连的结点）

D,ALL,ALL,0!约束所有不活动的结点自由度（可

选）

NSEL,ALL!选择所有结点

ESEL,ALL!选择所有单元

D,...!施加合适的约束

F,...!施加合适的活动结点自由度载荷

SF,...!施加合适的单元载荷

BF,...!施加合适的体载荷

SAVE

SOLVE

请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D,F,SF和BF命令得到更详细的解释。

后继载荷步

在后继载荷步中，用户可以随意杀死或重新激活单元。象上面提到的，要正确的施加和删除约束和结点载荷。

用下列命令杀死单元：

Command:EKILL

GUI:MainMenu>Solution>-LoadStepOpts-Other>KillElements

用下列命令重新激活单元：

Command:EALIVE

GUI:MainMenu>Solution>-LoadStepOpts-Other>ActivateElem

!第二个（或后继）载荷步：

TIME,...

ESEL,...

EKILL,...!杀死选择的单元

ESEL,...

EALIVE,...!重新激活选择的单元

...

FDELE,...!删除不活动自由度的结点载荷

D,...!约束不活动自由度

...

F,...!在活动自由度上施加合适的结点载荷

DDELE,...!删除重新激活的自由度上的约束

SAVE

SOLVE

请参阅TIME,ESEL,EKILL,EALIVE,FDELE,D,F和DDELE命令得到更详细的解释。

查看结果

对于大多数部分来说，用户在对包含不激活或重新激活的单元操作时应按照标准的过程来做。但是必须清楚的是，“杀

死”的单元仍在模型中，尽管对刚度（传导）矩阵的贡献可以忽略。因此，它们将包括在单元显示，输出列表等操作

中。例如，不激活的单元在结点结果平均（PLNSOL命令或MainMenu>GeneralPostproc>PlotResults>Nodal

Solu）时将“污染”结果。整个不激活单元的输出应当被忽略，因为很多项带来的效果都很小。建议在单元显示和其他

后处理操作前用选择功能将不激活的单元选出选择集。

使用ANSYS结果控制单元生死

在许多时候，用户并不清楚的知道杀死和重新激活单元的确切位置。例如，用户要在热分析中“杀死”熔融的单元（在

模型中移去溶化的材料），事先不会知道这些单元的位置；用户必须根据ANSYS计算出的温度确定这些单元。当决

定杀死或重新激活单元依靠ANSYS计算结果时（如温度，应力，应变等），用户可以使用命令识别并选择关键单元。

用下列方法识别关键单元：

Command:ETABLE

GUI:MainMenu>GeneralPostproc>ElementTable>DefineTable

用下列方法选择关键单元：

Command:ESEL

GUI:UtilityMenu>Select>Entities

然后用户可以杀死或重新激活选择的单元。（也可以用ANSYSAPDL语言编写宏以完成这些操作。）

用下列方法杀死选择的单元：

Command:EKILL,ALL

GUI:MainMenu>Solution>-LoadStepOpts-Other>KillElements

用下列方法重新激活选择的单元：

Command:EALIVE,ALL

GUI:MainMenu>Solution>-LoadStepOpts-Other>ActivateElem

下面的例子是杀死总应变超过许用值的单元：

/SOLU!进入求解器

...

...!标准的求解过程

SOLVE

FINISH

!

/POST1!进入POST1

SET,...

ETABLE,STRAIN,EPTO,EQV!将总应变存入ETABLE

ESEL,S,ETAB,STRAIN,0.20!选择所有总应变大于或等于0.20的单元

FINISH

!

/SOLU!重新进入求解器

ANTYPE,,REST

EKILL,ALL!杀死选择（超过允许值）的单元

ESEL,ALL!读入所有单元

...

...!继续求解

请参阅ETABLE,ESEL,ANTYPE和EKILL命令得到更详细的解释。

进一步的说明

不活动的自由度上不能施加约束方程[CE,CEINTF等]。（不活动的自由度当结点不与活动的单元相连时出现）。

可以通过先杀死然后重新激活单元的方法做应力松弛（如退火）操作。

在非线性分析中，注意不要因为杀死或重新激活单元引起奇异性（如结构分析中的尖角）或刚度突变。这将使得收敛困

难。

在有单元生死的分析中打开FULL牛顿－拉夫森方法的适应下降选项将得到好的结果。用下列方法：

Command:NROPT,FULL,,ON

GUI:MainMenu>Solution>AnalysisOptions

可以通过一个参数值来指示单元生死状态[*GET,Par,ELEM,n,ATTR,LIVE]

(UtilityMenu>Parameters>GetScalarData)。该参数可以用于APDL逻辑分支（*IF等），或其他要控制单元生

死的应用场合中。

用户可能想通过改变材料特性来杀死或重新激活单元[MPCHG](MainMenu>Preprocessor>Material

Props>ChangeMatNum)。但是，在这个过程中要特别小心。软件保护系统和限制使得“杀死”的单元在求解器中

改变材料特性时将不生效。（单元集中载荷不能自动删除；应变，质量，比热等也不能删除。）不当的使用MPCHG

命令将带来许多问题。例如，如果将单元的刚度缩减到近于0，而保留其质量，在有加速度和惯性载荷的问题中将产生

奇异性。

一个MPCHG的应用是在建立模型时涉及“出生”单元的应变历程的情况下。使用MPCHG可以得到单元在变形的结点

构造中的初始应变。

在单元生死中不能用多载荷步求解[LSWRITE]，因为不激活或重新激活的单元状态将不写入载荷步文件中。有多个载

荷步的生死单元分析应该用一系列的SOLVE命令（MainMenu>Solution>CurrentLS）来做。

单元生死应用实例（命令行格式）

问题描述

等截面杆两端固定，承受均匀的温度载荷时将其中间1/3段移去。过程是将其应变自由化并移去均匀温度。分析其

热应力和应变情况。材料特性和几何模型参数见下图。

/PREP7

/TITLE,ELEMENTBIRTH/DEATHINAFIXEDBARWITHTHERMALLOADING

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