首页/文章/ 详情

ANSYS中单元类型、材料表字段变量、非线性自适应网格定义命令详解!

7月前浏览602


在ANSYS workbench中,有时不能完全通过GUI形式对分析进行设置和处理,即GUI不能完全图形化处理APDL中的功能。为此,有时我们在使用workbench时需要插入一些命令流。


本公众 号在《ANSYS中关键字、实常数、材料属性和材料模型定义常用命令流!》一文中介绍了部分ANSYS中的命令。今天我们再讲讲几个常用命令的定义,包括:


① ET命令定义单元类型;

② TBFIELD定义材料表字段变量值;

③ NLADAPTIVE定义非线性自适应网格划分准则。



—01—

定义单元类型


ET命令,即Element Type,是在ANSYS中用于定义单元类型的命令。


ET命令的语法格式:

ET,ITYPE,Ename,KOP1,KOP2,KOP3,KOP4,KOP5,KOP6,INORP


其中:

 ITYPE—给选取的单元类型编号,默认为当前最大

              编号加1;

Ename—ANSYS单元类型库中的单元类型名称或

              编号,如beam188或单元编号,如188

              Ename=0,则定义一个空单元。

KOP1~KOP6用于控制单元的特定行为或特性的

             关键字,比如:激活或禁用某些功能、选

             择不同的计算模式等等。具体取值需参考

             单元库中对应单元的文档说明,如下各

             图。若需设置KEYOPT(7)及更高编号的

             选项,需使用KEYOPT命令单独定义。

INORP如果等于1,则抑制该单元类型所有的单

              元解输出,即不显示该类型单元的计算结

              果;默认为0(输出结果)。

 

在ANSYS的单元库帮助文档中,列出了各单元检索列表,如下图。可以通过检索列表查看各单元的具体内容。如KEYOPT关键字选项、自由度、单元特性等等。

举例,下图为SOLID185单元的KEYOPT(2)选项,用于设置该单元的计算模式:

输入ET命令举例:

ET,,BEAM188  !定义单元类型为beam188,取的

                            单元类型采用默认的编号

ET,1,MASS21  !定义单元类型为mass21,取的

                            单元类型采用编号为1




—02—

定义材料表字段变量值


TBFIELD命令是ANSYS中用于定义材料数据表(Table)的字段变量值的命令。


TBFIELD的主要作用是为材料后续的属性值定义提供字段变量类型和值(如温度、时间等物理变量和其值),以便动态调整材料在该字段变量类型和值下的材料数据(通过TBDATA命令定义材料数据)。


TBFIELD命令的语法格式:

TBFIELD, field_type, VALUE


其中:

field_type—为字段变量类型。常用的包括:
                   TEMP:温度,用于定义与温度相关

                   的材料属性(如热膨胀系数、弹性模

                   量随温度变化)。

                   TIME:时间,用于定义随时间变化的 

                   材料行为(如蠕变、磨损等)。
                   其他类型可能包括应力、应变等,

                   体可查询软件帮助文档

    VALUE—为指定字段变量的具体取值


输入命令举例:

①建立温度相关的材料属性:

TBFIELD, TEMP, 50    ! 定义温度为50℃
TBDATA, 1, 200          ! 在50℃下定义某属性值为200
TBFIELD, TEMP, 100  ! 定义温度为100℃

TBDATA, 1, 300         ! 在100℃下定义某属性值为300


②建立时间相关的材料属性(以接触磨损为例):

!材料参数磨损系数K、材料硬度H、压力指数m、速度指数n定义

K1=0

K2=1E-8

H=1000

m=1

n=1

! 定义磨损模型,其中MATID是与接触单元关联的材料ID,workbench中插入命令流时,特定接触面可用CID代替

TB,WEAR,MATID,,,ARCD     !激活接触磨损模型

TBFIELD,TIME,0             !字段变量TIME为0 时

TBDATA,1,K1,H,m,n       !定义材料磨损相关参数

TBFIELD,TIME,1            !字段变量TIME为1时

TBDATA,1,K1,H,m,n      !定义材料磨损相关参数

TBFIELD,TIME,1.01     !字段变量TIME为1.01时

TBDATA,1,K2,H,m,n     !定义材料磨损相关参数

TBFIELD,TIME,4           !字段变量TIME为4时

TBDATA,1,K2,H,m,n     !定义材料磨损相关参数


上述定义中,字段变量为TIME,当其值为0~1时为第一个载荷步,此时定义接触磨损系数为K1=0,当TIME值为1.01~4时,为第二个载荷步,此时定义接触磨损系数为K2=1E-8.


ANSYS帮助文档中关于TBFIELD命令内容更全面,具体应用建议查阅有关帮助文档,如下图。



—03—

定义非线性自适应网格划分准则


NLADAPTIVE命令是ANSYS中用于定义非线性自适应网格重新划分准则的命令。


该命令常用于大变形、接触、磨损等等复杂非线性分析问题中。通过该命令可自动调整网格质量,提高分析计算的计算精度和收敛性。


NLADAPTIVE命令语法格式:

NLADAPTIVEComponentActionCriterionOptionVAL1VAL2VAL3VAL4VAL5VAL6VAL7


其中:

component—指定参与自适应网格重划分的单元的

                      组件(element component),可通

                      过ESEL(element select)定义。
action—定义命令行为。包括:
              ADD:添加新的自适应准则。
                ON:启用自适应功能,并可指定检查时

                        间范围(val2为起始时间,val3为

                        结束时间)。
             LIST:列出当前定义的自适应准则。
criterion—定义自适应的判断标准:
                ENERGY:基于单元能量变化,当单元

                                能量超过阈值时网格重划分。
                      BOX:基于位置区域,当单元节点落

                                入指定的矩形或六面体区域内

                                时触发网格划分。
            CONTACT:基于接触单元数量或磨损

                                量,适用于接触或磨损模拟。

      MESH, SKEW:基于网格扭曲度,当单元扭曲

                                度超过阈值时触发。
option—根据不同的criterion选择相应的选项:对

             于ENERGY,option可省略,val1定义能

             量阈值。对于BOX,option为XYZRANGE,

             val1为区域尺寸,x1,x2,y1,y2,z1,z2定义区

             域边界。对于CONTACT,option为

             NUMELEM或WEAR,val1定义接触单元数

             量阈值或磨损量与单元高度的比率。对于

             MESH, SKEW,val1定义扭曲度阈值。
val1~val7—为具体数值根据criterion和option的组合而定,用于设置阈值、时间范围等。


输入命令举例:

注意以下举例中CP1~CP5为定义的组件的名称,具体分析项目中按自己定义的组件确定。


!基于能量的自适应准则
NLADAPTIVE, CP1, ADD, ENERGY, , 1.0e-3
!基于位置区域的自适应准则
NLADAPTIVE, CP2, ADD, BOX, XYZRANGE, 0.1, 0, 1, 0, 1, 0, 1
!基于接触单元数量的自适应准则
NLADAPTIVE, CP3, ADD, CONTACT, NUMELEM, 5

!基于接触磨损量的自适应准则

NLADAPTIVE,CP4,ADD, CONTACT,WEAR,0.50

!启用自适应功能并设置检查时间范围
NLADAPTIVE, CP5, ON, , , 1.0, 10.0


ANSYS帮助文档中关于NLADAPTIVE命令内容更全面,具体应用建议查阅有关帮助文档,如下图。




来源:薛定谔的Cube
ACTLS-DYNAWorkbenchAbaqus非线性UM材料控制
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-01-13
最近编辑:7月前
巴郡撸猫人
硕士 签名征集中
获赞 15粉丝 25文章 73课程 0
点赞
收藏
作者推荐

ANSYS摩擦接触理论知识详解!!

最近有点忙,许久未更新。今天给大家分享一篇ANSYS中的摩擦接触理论知识(本文源于ANSYS官方讲解)~--01--摩擦的理论模型两个接触体的切向滑动可能是无摩擦的,也有可能需要考虑摩擦。无摩擦时,两个物体相对滑动而没有任何阻力,摩擦存在时两个物体上会产生剪应力。我们知道,摩擦作用下,系统会沿着摩擦路径做功,消耗能量,因此摩擦是摩擦是一种与路径相关、具有能量耗散的物理现象。为准确地对摩擦问题进行分析,与摩擦相关的载荷的加载方式应该和实际情况相同,且子步要充分。原因就在于子步越充分,加载路径越精确,计算也就越准确。如下图所示:图1摩擦路径与时间步关系但需要注意,与塑性分析不同,摩擦分析中自动时间步长不考虑摩擦响应的增加。摩擦是一种机理比较复杂的物理现象,即使一个简单的摩擦力试验中,施加恒定的压力,其摩擦力-位移曲线也非常不规则,如下图2。通常摩擦与接触材料硬度、表面粗糙度、温度、湿度、相对速度等有关。因此摩擦分析中,我们仅考虑其主要数学物理原理,通常只能在数值上取得近似解。通常摩擦主要与以下物理参数相关:①接触材料(包括润滑材料);②表面粗糙度;③温度;④摩擦物体之间的相对速度。图2摩擦力—位移关系在ANSYS中,摩擦分析采用coulomb模型(库伦摩擦)施加摩擦力,并考虑了剪切摩擦力和粘合效应。库仑摩擦模型,分为弹性库仑摩擦(允许黏接和滑动)和刚性库仑摩擦(只允许滑动)。一、弹性库伦摩擦模型在ANSYS中,弹性库伦摩擦模型允许粘合(有摩擦力但还未发生滑移)和滑移,摩擦接触的两物体都存在反向的剪应力(TAU)。库伦摩擦准则要满足如下关系式:其中:μ—为摩擦系数;Fₜ—接触面切向力或剪应力;Fn—接触面的法向力。上述公式的Coulomb模型中,可以看到随着法向压力的增大,传递的最大剪应力TAUMAX(摩擦力)也随之增大,如下图所示。如果Fₜ超过上述准则,就会产生滑动。即:①粘合条件:Fₜ<μFₙ②滑移条件:Fₜ>μFₙ或TAU>TAUMAX(剪应力最大值)。图3库伦摩擦模型在ANSYS中,实常数TAUMAX默认值很大,为1E20,但一些单元可以通过命令定义TAUMAX的值。一般TAUMAX的上限值为σᵧ/(3)½,σᵧ为接触材料的米塞斯屈服应力。原因在于接触面之间的剪切屈服限制了其剪应力的大小。如下图,橡胶的剪切屈服使出现打滑,并且这也限制了轮胎和道路之间可以传逆的摩擦力。图4摩擦力超过轮胎材料最大剪应力工程上,TAUMAX一般使用经验值,以满足实际情况。下图为弹性库伦摩擦的压力—剪应力图,可以看到在TAUMAX之前,摩擦力或剪应力与压力P成正比。如果我们在计算时不修改最大的TAUMAX为实际值,由于程序的默认值很大,程序实际上的摩擦力就等于摩擦系数乘以法向力(虚线)。图5压力与剪应力关系从上可知,在ANSYS中粘合区域被视作弹性处理,其切向刚度为KT。切向刚度和法向刚度类似:刚度越大精度越高,刚度越小则越容易收敛。软件可以指定KT的值,或使用程序指定的值(KT=1%KN)。图6摩擦力与位移的关系二、带cohesive内聚力的库伦摩擦模型在某些情况下,接触表面粘合在一起,如胶粘结构,即使没有法向压力的作用也能提供一定大小的滑移阻力。ANSYS中的一些单元能够模拟这种现象(使用内聚力COHE)。图7带内聚力的压力与剪应力关系带cohesive内聚力的库伦摩擦模型满足如下关系式:其中:①P—为接触面法向压力;②COHE—为接触部位内聚力。三、刚性库伦摩擦模型刚性库伦摩擦没有考虑静摩擦,因此该模型的接触无“粘合”效应。该模型仅适合模型在某固定方向连续滑动时。如使用研磨轮对部件进行成型加工时(图9),u=0处的不连续等效于无穷大的刚度值。如果滑动停止或方向改变,则会出现收敛困难。图8刚性库伦模型的摩擦力与位移关系图9研磨加工--02--摩擦系数和与其相关的实常数一、摩擦系数对于所有的ANSYS接触单元,摩擦系数μ通过属性MU来指定,(缺省时为μ=0)。滑动时的摩擦系数比静止时的小,我们称作动摩擦系数,静止时,称作静摩擦系数。ANSYS中,面一面单元(CONTA171-CONTA174)和点一面单元(CONTA175)可以指定一个和表面滑动速度相关的动摩擦系数。静摩擦到动摩擦的过程为避免收敛问题,需要考虑其整体数学描述时的下降过程,ANSYS中定义摩擦系数方程如下:其中:MU—为动摩擦系数,用户定义;MSU—为静摩擦系数;FACT—为静摩擦系数与动摩擦系数之比,即MSU/MU,缺省为1;DC—为衰减系数用户定义,缺省为0;Vrel—为表面间滑动的相对速度。如果我们知道静、动摩擦系数和至少一个数据点(μ1,Vrel),则可以确定摩擦衰减系数为:从上述内容可以看到,动摩擦系数和表面速度相关,这使得动静态之间可以平滑过度,如下图。图10摩擦系数与相对滑动速度的关系需要注意,在workbench中,接触选项中(下图)我们输入的摩擦系数通常为动摩擦系数MU。二、摩擦接触的一些实常数ANSYS软件中为摩擦接触提供了一些实常数,以便我们对数值模型有关参数进行调整。关于库伦摩擦接触模型的一些实常数如下:在workbench中,可插入Command,使用RMODIF设置接触实常数,命令如下:RMODIF,NSET,STLOC,VALUENSET—表示接触对的实常数号,特定面可以用系统默认的CID或TID;STLOC—表示实常数序号,即上表中实常数表格的位置No.,对应了具体实常数Name;VALUE—表示实常数定义的值。对于一个接触其上述实常数可按如下方式定义:RMODIF,CID,9,200E6!设置最大摩擦应力TAUMAXRMODIF,CID,13,100E5!设置最大粘结应力COHERMODIF,CID,21,1.2!设置静动态摩擦系数之比FACTRMODIF,CID,21,1.2!设置摩擦衰减系数DC--03--摩擦生热摩擦现象可以将机械能转换为热能,为了模拟摩擦耗散能量的热生成,在ANSYS中应当执行瞬态热—结构耦合分析。在分析时,可以应用TIMINT,STRUC,OFF命令关闭结构自由度上的瞬态效应,但必须保留热自由度上的瞬态效应。摩擦生热需要两个实常数,即:①FHTG—摩擦耗散能转换成热量的比例系数,即机械能转化为热能的效率(默认为1.0);②FWTG—接触面和目标面之间热量分布的权重系数(默认为0.5)。在模拟热—结构耦合的问题中,ANSYS通过下式定义摩擦生热总耗散量:其中:τ—等效摩擦应力;υ—相对滑动速度;FHGT—摩擦耗散能转换成热量的比例系数,默认为1.0,并且可以作为实常数输入。对FHGT=0的情况,软件中应该输入一个很小的数值,如果输入0软件将采用默认值。在接触面和目标面上摩擦耗散的能量按计算公式分配,如下:其中:qᴄ—为接触面上的摩擦耗散能;qᴛ—为目标面上的摩擦耗散能。摩擦生热中,实常数FWGT和FHGT在ANSYS中可按如下形式定义:RMODIF,CID,18,0.5!设置摩擦热分布的权重系数RMODIF,CID,20,0.6!设置摩擦生热能量转化系数--04--关于摩擦问题的求解在力学求解中,无摩擦时,刚度矩阵法向力与切向力相互独立,形成的刚度矩阵都是对称矩阵,但若考虑摩擦接触问题,切向力与法向力通过摩擦系数μ耦合,会导致不对称的刚度矩阵,如下图。这时在求解时,我们就需要调用非对称求解器,但使用非对称求解器比对称求解器更费时间。因此,ANSYS使用了对称化的算法,这样可使得包括摩擦的接触问题也使用对称矩阵进行求解。来源:薛定谔的Cube

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈