有限元分析中常用单元类型与单位制

1、ansys结构分析之单元速查 link1可承受单轴拉压的单元，不能承受弯矩作用plane22 维6节点三角形实体结构单元，可用作平面单元 (平面应力或平面应变)，也可以用作轴对称单元beam3可承受拉、压、弯作用的单轴单元，每个节点有三个自由度，即沿x,y方向的线位移及绕z轴的角位移beam4 承受拉、压、弯、扭的单轴受力单元，每个节点上有六个自由度：x、y、z三个方向的线位移和绕x,y,z三个轴的角位移solid5三维耦合场体单元，8个节点，每个节点最多有6个自由度link8三维杆（或桁架）单元，用来模拟：桁架、缆索、连杆、弹簧等等，是杆轴方向的拉压单元，每个节点具有三个自由度：沿节点坐标系

2、x、y、z方向的平动plane132 维耦合场实体单元，有 4 个节点，每个节点最多有 4 个自由度plane254 节点轴对称谐波结构单元，用于承受非轴对称载荷的 2 维轴对称结构的建模link32二维热传导杆单元，应用在二维（板或轴对称）稳态或瞬态热分析plane352 维 6 节点三角形热实体单元，用作平面单元或轴对称单元plane422 维实体结构单元，作平面单元 (平面应力或平面应变)，也可以用作轴对称单元。本单元有 4 个节点，每个节点有 2 个自由度，分别为 x 和 y 方向的平移shell43 4 节点塑性大应变单元，适合模拟线性、弯曲及适当厚度的壳体结构。单元中每个节点具有六

3、个自由度：沿xy和z 方向的平动自由度以及绕x、y和z 轴的转动自由度 plane532 维 8 节点磁实体单元，用于 2 维 (平面和轴对称) 磁场问题的建模plane552 维 4 节点热实体单元，作为平面单元或轴对称环单元，用于 2 维热传导分析。本单元有 4 个节点，每个节点只有一个自由度 温度shell63弹性壳单元，具有弯曲能力和又具有膜力，可以承受平面内荷载和法向荷载。本单元每个节点具有6个自由度：沿节点坐标系x、y、z方向的平动和沿节点坐标系x、y、z轴的转动solid643-d 各向异性结构实体单元，用于各向异性实体结构的3d建模。单元有8个结点，每个结点3个自由度，即沿x、

4、y、z的平动自由度solid65用于含钢筋或不含钢筋的三维实体模型。该实体模型可具有拉裂与压碎的性能plane672 维热-电耦合实体单元，有 4 各节点，每个节点两个自由度：温度和电压plane754 节点轴对称谐波热单元，作轴对称环单元，具有 3 维热传导能力。本单元有 4 个节点，每个节点只有一个自由度 温度 tempplane772 维 8 节点热实体单元，2 维 4 节点热单元 (plane55) 的高阶版本。每个节点只有一个自由度 温度plane788 节点轴对称-谐波热单元，轴对称环单元，具有 3 维热传导能力。本单元有 8 个节点，每个节点只有一个自由度 温度 tempplan

5、e822 维 8 节点结构实体单元,是 2 维 4 节点单元 (plane42) 的高阶版本。对于四边形和三角形混合网格，它有较高的结果精度；可以适应不规则形状而较少损失精度。本 8 节点单元具有一致位移形状函数，能很好地适应曲线边界plane838 节点轴对称谐波结构实体单元，非轴对称载荷的 2 维轴对称结构的建模，有 8 个节点，每个节点有三个自由度 节点在 x, y 和 z 方向的平移。对于未转动的节点坐标，其方向分别对应径向、轴向和切线方向 (圆周方向)plane1212 维 8节点静电单元。本单元每个节点只有一个自由度：电压，用于 2 维静电场分析plane1452 维四边形结构实体

6、 p 单元，支持最多 8 阶多项式，用作平面单元 (平面应力或平面应变) 或作为轴对称单元plane1462 维三角形结构实体 p 单元，支持最多 8 阶多项式，作平面单元 (平面应力或平面应变) 或作为轴对称单元plane162显式动力 2 维结构实体，用于平面问题，也可用于轴对称问题。本单元有4 个节点，每个节点 6 个自由度：节点在 x 和 y 方向的平移、速度和加速度plane182单元限制 2 维 4 节点结构单元，可用作平面单元 (平面应力、平面应变或广义平面应变)，也可作为轴对称单元，具有塑性、超弹性、应力刚度、大变形和大应变能力，并具有力-位移混合公式的能力，可以模拟接近不可压

7、缩的弹塑性材料的变形plane1832 维 8 节点实体结构单元，用作平面单元 (平面应力、平面应变和广义平面应变)，也可用作轴对称单元。本单元具有塑性、蠕变、应力刚度、大变形及大应变的能力。并具有力-位移混合公式的能力，可以模拟接近不可压缩的弹塑性材料的变形beam188 3 维线性有限应变梁单元，适合于分析从细长到中等粗短的梁结构plane2232 维 8 节点耦合场单元， 2 维结构、电、压组和压电分solid453-d结构实体单元产品：mp me st pr pp edsolid45单元说明solid45单元用于构造三维实体结构.单元通过8个节点来定义,每个节点有3个沿着xyz方向平移

8、的自由度.单元具有塑性,蠕变,膨胀,应力强化,大变形和大应变能力。有用于沙漏控制的缩减积分选项。有关该单元的细节参看 ansys, 理论参考中的solid45部分。类似的单元有适用于各向异性材料的solid64单元。solid45单元的更高阶单元是solid95。图45.1 solid45几何描述solid45输入数据 该单元的几何形状、结点位置、坐标系如图45.1: solid45 几何描述所示。该单元可定义8个结点和正交各向异性材料。正交各向异性材料方向对应于单元坐标方向。单元坐标系方向参见坐标系部分。 单元荷载参见结点和单元荷载部分。压力可以作为表面荷载施加在单元各个表面上，如图45.1

9、: solid45 几何描述所示。正压力指向单元内部。可以输入温度和流量作为单元节点处的体载荷。节点 i 处的温度 t(i) 默认为 tunif。如果不给出其它节点处的温度，则默认等于 t(i)。对于任何其它的输入方式，未给定的温度默认为 tunif。对于流量的输入与此类似，只是默认值用零代替了 tunif。keyopt(1)用于指定包括或不包括附加的位移形函数。keyopt(5)和keyopt(6)提供不同的单元输出选项（参见单元输出部分）。 当keyopt(2)=1时，该单元也支持用于沙漏控制的均匀缩减（1点）积分。均匀缩减积分在进行非线性分析时有如下好处： 相对于完全积分选项而言，单元刚

10、度集成和应力（应变）计算需要更少的cpu时间，而仍能获得足够精确的结果。 当单元数量相同时，单元历史存储记录（.esav 和 .osav）的长度约为完全积分（222）的1/7。 非线性分析的收敛性通常远比采用额外位移形状的完全积分要好；即，keyopt(1) = 0, keyopt(2) = 0。 分析结果不会受（由塑性或其它不可压缩材性引起的）体积锁死的影响。采用均匀缩减积分有以下缺点： 当采用相同网格进行弹性分析时，结果显然不如完全积分方法精确。 采用单层单元时不能很好的得到结构的弯曲特性（例如，一根悬臂梁，受横向集中力，采用单层单元）。建议采用4层单元。当采用均匀缩减积分选项时（keyo

11、pt(2) = 1 这和 solid185 用 keyopt(2) = 1是一样的），应对总能量 (etable 命令，sene 标识符)和沙漏造成的伪能量 (etable 命令，aene 标识符) 进行比较以检查结果的精度。如果沙漏能与总能量之比小于 5%，结果一般是可以接受的。如果该比值超过5%，则需细化网格。也可以在求解阶段用 outpr,veng 命令控制总能量和沙漏能。更多说明见 ansys 理论手册。 可用istress 或者 isfile 命令给单元施加初始应力状态。更多信息参见ansys基础分析纸指南中的施加初始应力部分。你也可以用keyopt(9) = 1来从用户子程序中读取

12、初应力ustress。关于用户子程序的细节，参见ansys upf指南。 在进行几何非线性分析时，可以使用 solcontrol,incp 命令来包含压力的影响。在线性特征值屈曲分析中自动包括压力载荷刚度效应。如果需要非对称的压力载荷刚度效应矩阵，使用 nropt,unsym 命令。该单元的输入概要参见 solid45 输入数据摘要. 单元输入数据的一般性描述参见单元输入部分。solid45单元输入数据摘要节点 i，j，k，l，m，n，o，p自由度 ux，uy，uz实常数hgstf沙漏控制因子，仅当keyopt(2) = 1时需要设置。注：有效值为任意正数，默认为1.0。建议值为1到10之间。

13、材料参数 ex, ey, ez, prxy, pryz, prxz (或 nuxy, nuyz, nuxz), alpx, alpy, alpz (或 ctex, ctey, ctez or thsx, thsy, thsz), dens, gxy, gyz, gxz, damp表面载荷 压力表面1（j-i-l-k），表面2（i-j-n-m），表面3（j-k-o-n），表面4（k-l-p-o）,表面5（l-i-m-p），表面6（m-n-o-p）体载荷 温度t（i），t（j），t（k），t（l），t（m），t（n），t（o），t（p）；流量fl（i），fl（j），fl（k），fl（l），fl（m

14、），fl（n），fl（o）fl（p）特殊功能 塑性蠕变膨胀应力刚化大变形大应变单元死活自适应下降初始应力输入keyopt (1) 0 包括附加的位移形函数 1 不包括附加的位移形函数keyopt (2) 积分选项 0 依据keyopt(1)带或者不带附加位移形函数，执行完全积分 1 带砂漏控制的均匀缩减积分,不带附加的位移形函数(keyopt(1)自动设置为1)keyopt(4) 单元坐标系0单元坐标系平行于整体坐标系 1单元坐标系基于单元ij边keyopt(5) 额外单元输出0 基本单元解 1 在所有积分点上重复基本解 2节点应力解 keyopt(6) 额外表面输出0基本单元解1附带表面i-

15、j-n-m的表面解2表面i-j-n-m和表面k-l-p-o的表面解(表面解只对线性材料可用) 3附带每个积分点的非线性解4非零压力表面的表面解keyopt(9) 初始应力子程序选项（仅适用于直接用keyopt命令输入时）0没有用户子程序提供初应力(默认) 1从用户子程序ustress中读入初始应力数据(有关用户子程序参见ansys upf指南)solid45输出数据与单元结果相联系的结果输出主要有两种方式：l 节点位移和所有节点结果。l 附加的单元输出，如表45.1: solid45 单元输出定义所示。 图 45.2: solid45应力输出显示了几项结果。单元应力方向平行于单元坐标系。表面应

16、力输出在表面坐标系上，各面上的结果都可得到（keyopt(6)）。面ijnm和klpo的坐标系如图45.1: solid45几何描述所示。其他的表面坐标系遵从类似的定位，即由受压面结点关系确定。表面应力输出仅当满足单元输出的条件时可用。一般性的描述参见结果输出。关于如何查看结果，参见ansys 基础分析指南。图45.2: solid45应力输出(这里显示的应力方向相应于 keyopt(4) = 0)当keyopt(2)=1（单元采用均匀缩减积分），单元积分点上所有的输出和完全积分方法的输出形式相同。为保证输出的一致性，在相同单元类型中完全积分方法的点号。单元输出定义表使用如下标记:在名称列中的

17、冒号 (：) 表示该项可以用分量名方法 etable, esol 处理；o 列表示该项可用于 jobname.out 文件；r 列表示该项可用于结果文件。无论 o列或 r 列，y 表示该项总是可用的，一个数字表示表的一个注解,其中说明了使用该项的条件；而减号 - 表示该项不可用。表45.1 solid45单元的结果输出定义名称定义orei单元号yynodes节点-i,j,k,l,m,n,o,pyymat材料编号yyvolu体积yyxc,yc,zc结果输出位置y3pres节点j,i,l,k的p2;jkon的p3;klpo的p4;limp的p5;mnop的p6yytemp温度t(i) ,t(j),

18、t(k),t(l),t(m),t(n),t(o),t(p)yyfluen流量fl(i) fl(j) fl(k) fl(l) fl(m) fl(n) fl(o) fl(p)yys:x,y,z,xy,yz,xz应力yys;1,2,3主应力yys:int应力强度yys:eqv等效mises应力yyepel:x,y,z,xy,yz,xz弹性应变yyepel:1,2,3主弹性应变y-epel:eqv等效弹性应变yyepth:x,y,z,xy,yz,xz平均热应变11epth:eqv等效热应变11eppl:x,y,z,xy,yz,xz平均塑性应变11eppl:eqv等效塑性应变11epcr:x,y,z,x

19、y,yz,xz平均蠕变应变11epcr:eqv等效蠕变应变11epsw:平均膨胀应变11nl:epeq平均等效塑性应变11nl:srat屈服表面上的迹应力和应力之比11nl:sepl从应力-应变曲线平均等效mises应力11nl:hpres静水压力1face表面lable22area表面面积22temp表面平均温度22epel表面弹性应变(x,y,xy)22press表面压力22s(x,y,xy)表面应力(x轴平行于定义该表面的前面两个结点连接)22s(1,2,3)表面主应力22sint表面应力强度22seqv表面等效mises应力22loci:x,y,z积分点位置-y1 仅当单元具有非线性材

20、料时，有非线性解输出。2 表面输出(如果keypoint(6)是1、2或者4) 。3 用*get可以获得质心位置的结果.。 4等效应变用一个等效泊松比来计算：对于弹性和热问题，这个值由(mp,prxy)输入；对于塑性和蠕变问题，这个值取0.5表45.2solid45综合单元输出描述输出项名or非线性积分点结果eppl, epeq, srat, sepl, hpres, epcr, epsw1-积分点应力结果temp, s(x, y, z, xy, yz, xz), sint, seqv, epel2-结点应力结果 temp, s(x, y, z, xy, yz, xz), sint, seqv

21、, epel3-1. 若单元有非线性材料且keyopt(6)=3，在8个积分点处输出结果。2. 若keyopt(5)=1，在每个积分点输出结果。3. 若keyopt(5)=2，在每个结点输出结果。表 45.3 solid45 输出项和序列号 列出了可以通过 etable 命令，用序列号方法输出的内容列表。更多信息见 ansys 基本分析指南 中 一般后处理 (post1) 部分和本手册中 输出项和序列号表 部分。在 表 45.3 solid45 输出项和序列号 中使用如下标识符：name 与表 45.1: solid45 单元输出定义 中相同定义的输出量；item 用于 etable 命令的预

22、先定义的输出项；e对于单值或常数型单元数据的序列号；i,j,p 节点i,j,.,p 处数据的序列号；表45.3solid45 输出项和序列号输出量名etable和 esol 命令输入itemijklmnopp1smisc2143-p2smisc56-87-p3smisc-910-1211-p4smisc-1314-1615p5smisc18-1719-20p6smisc-21222324s:1nmisc16111621263136s:2nmisc27121722273237s:3nmisc38131823283338s:intnmisc49141924293439s:eqvnmisc51015

23、2025303540fluennmisc4142434445464748 有关适用etable命令的表面结果输出参见本手册的表面输出部分。solid45假定和限制 体积等于0的单元是不允许的。 单元结点编号可参照图45.1: solid45几何描述，面ijkl和 mnop也可互换。 单元不能扭曲，这样单元就会有两个独立的体。这通常发生在当单元结点编号不当时。 所有单元都必须有8个结点。 可以通过定义重合的k和l、o和p来形成棱柱形单元（参见三角形、棱柱形和四面体单元）。 四面体形状也是允许的。对四面体单元，额外形状被自动删除。solid45产品限制对于以下产品，将在上述一般假设和限制的基础上再

24、增加一定的限制：。ansys professional. 不允许有damp材性。 不能施加流量体荷载。 唯一允许的特殊性能是应力钢化。 keyopt(6)=3不可用。mpc184多点约束单元：刚性杆，刚性梁，滑块，球铰，销轴和万向联轴器mp me st pr pp edmpc184单元描述mpc184包括使用拉格朗日乘子法实现运动约束的一类常用的多点约束单元。这些单元可以简单地分为“约束单元”或“连接单元”。用户可以在一些需要施加运动约束的场合中使用这些单元。这些约束可以简单到铰链上的具有相同的位移值，也可以复杂到包括模型的刚性部分，或者在柔性体之间以某一特定方式传递运动的运动约束。例如，结构

25、中可能包含一些刚性部件或者通过转动或滑块约束连接在一起的运动部件。结构的刚性部分可以使用mpc184的刚性杆或刚性梁单元来模拟，运动部分可以使用mpc184的滑块，球铰，销轴和万向联轴器单元模拟。因为这些单元使用拉格朗日乘子法实现，ansys能够输出约束反力和力矩。使用的约束类型取决于实际的应用场合。约束单元如果没有其它说明，使用这些单元时，三维单元选项(keyopt(2) = 0)为默认值。刚性杆/梁图184.1mpc184几何mpc184可以用来模拟两个变形体之间的刚性约束，或者在工程模型中模拟传递力和力矩的刚性部件。这个单元也可以应用在由其它刚性结构热膨胀引起的变形分析中。如果keyop

26、t(1) = 0 (默认值)，单元为刚性杆，有二个节点，每个节点有三个自由度（仅有平移）。如果keyopt(1) = 1，单元为刚性梁，有二个节点，每个节点有六个自由度（x，y，z三个方向的平移和转动）。该单元适合线性，大转动，大应变非线性情况。 滑块设置keyopt(1) = 3来定义三节点的滑块单元。3维滑块单元(keyopt(2) = 0)每个节点有三个自由度(x，y和z方向平移)。2维滑块单元(keyopt(2) = 1)每个节点有二个自由度(x，y方向平移)。滑块单元的运动约束施加方式：节点k被定义为从节点，该节点的运动轨迹始终在两个主节点（i和j）的连线上。k节点只允许在i和j节点

27、连线上滑动。球铰设置keyopt(1) = 5来定义二节点的球铰。两个节点必须重合。3维球铰每个节点有三个自由度(x，y和z方向平移)。2维球铰单元(keyopt(2) = 1)每个节点有二个自由度(x，y方向平移)。球铰单元的运动约束施加方式：组成单元的两个节点平移位移保持一致；不约束转动自由度（如果存在）。说明也可以使用ce或cp命令指定二个节点具有相同的位移值，在那种情况下，约束自由度被删除。然后，对球铰单元而言，约束通过拉格朗日乘子法施加，可以得到约束反力。使用ce或cp命令施加位移约束通常更有效率，如果可能，尽量使用ce或cp命令代替mpc184球铰单元。链接单元数值模拟中经常涉及两

28、个零件连接的问题。这些链接所承受的运动约束可以简单到两个零件连接处的具有相同的位移值，也可以复杂到模型中两个变形体之间运动的传递。复杂的链接也可以包括一些控制机构，如限制器或挡块，两个物体相对运动分量的锁定器。在许多情况下，这些链接还可以包括在两个物体相对运动的未约束分量之间的刚度，阻尼或摩擦行为特性。销轴链接和万向联轴器广泛应用于汽车，机器人，生物工程和其它行业。这些链接单元有二个节点。根据链接的不同定义，在二个节点上施加不同的运动约束。约束通过拉格朗日乘子法施加。链接单元在每个节点上有六个自由度，定义六个相对运动分量：三个相对平移和三个相对转动。这六个相对运动分量主要用于模拟链接单元运动行

29、为。对于链接单元，一些分量可能被运动约束限制住了，而另一些分量是“自由的”或者“无约束的”。万向联轴器和销轴链接单元中，二个节点假设为重合，相对位移为零。销轴链接单元只有一个相对旋转运动分量绕着销轴旋转是没有约束住的，而万向联轴器有两个旋转分量是自由的。这些单元包括一些控制特性，如挡块，锁定器。一些激励载荷或边界条件也可以施加到单元的二个节点的相对运动分量上去。例如，在销轴链接中，可以在销轴旋转方向上设置挡块，限制旋转角度在一个特定范围内。也可以在二个节点旋转分量上施加位移或力边界条件模拟链接的驱动特性。驱动力或位移由实际机构电动或液压系统提供。链接单元的未约束相对运动分量上可以施加线性刚度和

30、阻尼特性。如果需要，可以施加与温度相关的刚度和阻尼特性。除了ansys中已经存在的输出选项，链接单元的相对运动分量也可以输出。该单元适合线性，大转动，大应变非线性情况。 销轴链接设置keyopt(1) = 6定义二节点销轴链接。销轴单元的二个节点必须有相同的空间坐标。mpc184销轴链接单元只有一个基本自由度绕着轴或销相对旋转。单元能够包括控制特性，如未约束自由度上的挡块，锁定器。旋转边界条件也可以施加到相对运动分量上。另外，相对旋转方向可以施加线性刚度和阻尼特性。ansys单元库中另一种旋转链接单元是combin7。combin7单元允许用户控制连接柔度，摩擦，阻尼和其它控制特性。局部坐标系

31、固定在链接上，并与之一起运动，所以，该单元可以应用于大变形分析中，(更多详细内容参见combin7 - 销轴链接单元)。mpc184销轴链接单元的节点上施加了运动约束确保节点有相同的位移，另外，只有一个相对旋转轴的转动自由度未被约束，另外二个转动被固定。万向联轴器设置keyopt(1) = 7来定义二节点方向联轴器单元。单元的二个节点必须有相同的空间坐标。mpc184万向联轴器单元有两个自由的相对转动自由度。该单元包括控制特性：相对运动旋转分量上的挡块和锁定器。也可以在单元两个相对旋转分量上施加旋转边界条件。另外，可以在单元相对旋转分量上指定线性刚度和阻尼特性。mpc184输入数据输入数据因为

32、mpc184约束和链接的不同而不同。约束输入数据刚性杆/刚性梁图184.1mpc184几何显示了该单元的几何，节点位置和单元坐标系。单元需定义二个节点。单元的x轴从节点i指向节点j。单元截面面积假设为一个单位。ansys自动选择截面坐标系；选择方法参见beam4中的描述。当单元作为刚性梁使用时，截面坐标系仅仅与弯矩输出有关。因为单元模拟一个刚性约束或刚性部件，材料刚度特性不必输入。如果要考虑热膨胀效应，热膨胀系数必须指定。如果分析中要考虑刚性单元的质量，密度必须指定。如果密度指定，ansys在单元使用集中质量矩阵。该单元支持单元生死特性，使用ealive和ekill实现该功能。节点和单元载荷描

33、述单元载荷。可以在节点上输入温度作为单元体载荷。节点i温度t(i)，默认值为tunif。节点j温度默认值为t(i)。mpc184输入概述 总结了单元的输入数据。mpc184输出数据 提供了单元输出的一般描述。滑块图184.2: mpc184滑块约束几何 显示单元的几何形状和节点位置。单元由三个节点(i，j和k)定义。假设节点i在节点j和k的连线上。图184.2: mpc184滑块约束几何该单元不需要定义材料刚度特性。目前不支持单元生死选项。 mpc184滑块单元输入概要总结了单元输入参数。 mpc184滑块单元输出数据提供单元输出的常用描述。球铰图184.2: mpc184球铰约束几何显示单元

34、的几何形状和节点位置。单元由二个节点(i，j和k)定义。假设节点二个节点(i和j)具有相同的空间坐标。图184.2: mpc184球铰约束几何该单元不必输入材料刚度特性，目前不支持单元生死。mpc184球铰单元输入概要总结了单元输入参数。 mpc184球铰单元输出数据提供单元输出的常用描述。链接输入数据所有mpc184链接单元有一些通用的输入要求。每个链接单元的特殊要求将在相应章节着重显示。截面定义每一个链接单元必须有相应的截面定义。使用sectype定义截面类型和子类型。局部坐标系说明经常在每个节点的局部坐标系定义链接单元的运动约束。使用命令完成。局部坐标系和方向对每种链接单元来说可能会不一样。相应的每个链接单元的输入数据中会注明要求。典型的局部坐标系总是定义在链接单元的第一个节点上。第二个节点的局部坐标系是可选的，如果没有定义它，假定和第一节点的坐标系一致。链接单元二个节点相对运动的旋转分量用bryant角(或cardan角)表示，且在局部坐标系下表