Ansys接触分析和设置

Ansys非线性接触分析和设置

设置实常数和单元关键选项

程序利用20个实常数和数个单元关键选项，来操纵面一面接触单元的接触。

参见《ANSYSElementsReference））中对接触单元的描述。

实常数

在20个实常数中，两个（R1和R2）用来概念目标面单元的几何形状。剩下的

用来操纵接触面单元。

R1和R2概念目标单元几何形状。

FKN概念法向接触刚度因子。

FTOLN是基于单元厚度的一个系数，用于计算许诺的穿透。

ICONT概念初始闭合因子。

PINB概念“Pinball”区域。

PMIN和PMAX概念初始穿透的允许范围。

TAUMAR指定最大的接触摩擦。

CNOF指定施加于接触面的正或负的偏移值。

FKOP指定在接触分开时施加的刚度系数。

FKT指定切向接触刚度。

COHE制定滑动抗力粘聚力。

TCC指定热接触传导系数。

FHTG指定摩擦耗散能量的热转换率。

SBCT指定Stefan-Boltzman常数。

RDVF指定辐射观看系数。

FWGT指定在接触面和目标面之间热散布的权重系数。

FACT静摩擦系数和动摩擦系数的比率。

DC静、动摩擦衰减系数。

命令：R

GUI：mainmenu>preprocessor>realconstant

对实常数FKN,l-'TOLN,ICONT,PINB,PMAX,PMIN,FKOP和FKT,用户既

能够概念一个正值，也能够概念一个负值。程序将正值作为比例因子，将负值作

为绝对值。程序将下伏单元的厚度作为ICON,FTOLN,PINB,PMAX和PMIN的

参考值。例如ICON二说明初始闭合因子是“*基层单元的厚度”。但是，ICON

二那么表示真实调整带是单位。若是下伏单元是超单元，那么将接触单元的

最小长度作为厚度。参见图5-8o

o1attKhed

•lement

图5-8基层单元的厚度

在模型中，若是单元尺寸转变专门大，而且在实常数如ICONT,FTOLN,PINB,

PMAX,PMIN中应用比例系数，那么可能会显现问题。因为从比例系数取得的实

际结果，取决于基层单元的厚度，这就可能引发大、小单元之间的重大转变。若

是显现这一问题，请用绝对值代替比例系数。

TCC,FHTG,SBCT,RDVF和FWGT仅用于热接触分析［KEYOPT(1)=1］。

单元关键选项

每种接触单元都包括数个关键选项。对大多的接触问题，缺省的关键选项是

适合的。而在某些情形下，可能需要改变缺省值。下面是能够操纵接触行为的一

些关键选项：

自由

度KEYOPT(l)

接触算法(罚函数+拉格朗日乘子或罚函数)KEYOPT(2)

存在超单元时的应力状态(仅2D)KEYOPT(3)

接触检测点的位置(仅低阶接触单元)KEYOPT(4)

CNOF自动调

整KEYOPT(5)

时刻步操纵KEYOPT(7)

伪接触预防KEYOPT(8)

初始穿透或间隙的阻碍KEYOPT(9)

法向和切向接触刚度修正方式操纵KEYOPT(10)

壳的厚度阻碍KEYOPT(11)

接触面行为(粗糙、绑定等)KEYOPT(12)

命令：KEYOPT

ET

GUT：mainmenu>preprocessor>ElemantType>Add/Edit/Delete

选择接触算法

对面一面接触单元，程序能够利用增进的拉格朗日方式或罚函数方式。通过

单元关键字KEY0PT（2；来指定。

增进的拉格朗日方式是为了找到精准的拉格朗日乘子（即接触力），而时罚

函数进行一系列修正迭代。与罚函数的方式相较，拉格朗日方式容易患到良态条

件，对接触刚度的灵敏性较小。但是，在有些分析中，增进的拉格朗日方式可能

需要更多的迭代，专门是在变形后网格变得太扭曲时。

利用拉格朗日方式的同时应利用实常数FTOLNoFT0LN为拉格朗日方式指定

允许的最大穿透。若是程序发觉穿透大于此值时，即便不平稳力和位移增量已

经知足了收敛准则，总的求解仍被看成不收敛处置。FTL0N的缺省值为。用户

能够改变那个值，但要注意，若是此值过小，可能会造成太多的迭代次数或不。

收敛

确信接触刚度

所有的接触问题都需要概念接触刚度，两个表面之间穿透量的大小取决于接

触刚度。过大的接触刚度可能会引发总刚矩阵的病态，从而造成收敛困难。一样

来讲，应该选取足够大的接触刚度以保证接触穿透小到能够同意，但同时又应该

让接触刚度足够小以不致引发总刚矩阵的病态而保证收敛性。

ANSYS程序依照下伏柔体单元的材料特件,夹估量一个缺省的接触刚度值.

用户可用实常数FKN来为接触刚度指定一个比例因子或指定一个绝对值。比例

因子一样在和10之间；关于大变形问题，选1是比较好的；而关于弯曲为主的

问题，一样为〜。用户应当老是查验以使穿透抵达极小值，而又幸免过量的迭代

次数。

注意一FTOLN和FKN从一个荷载步到另一人荷载步中，都能够修改。也能

够在重启动中修改。这时，必需概念KEY0PT（10）=l,2。

为了确信一个较好的接触刚度值，可能需要一些体会。用户能够按下面的步

骤来进行尝试：

一、开始时取一个较低的值。低估值要比高估值好，因为由一个较低的接

触刚度致使的穿透问题，比太高的接触刚度致使的收敛性困难，要容易

解决。

二、对前几个子步进行计算分析•，直到最终荷载的一个比例（恰好完全成立

接触）。

3、检查每一子步中的穿透量和平稳迭代次数。若是整体收敛困难是由过大

的穿透引发的（而不是由不平稳力和位移增量引发的），那么可能低估了

FKN的值，或是将FTOLN的值取得大小。若是整体的收敛困难是由于

不平稳力和位移增量达到收敛值时需要过量的迭代次数，而不是由于过

大的穿透量引发的，那么FKN的值可能被高估。

4、按需要调整FKN或FTOLN的值，从头进行完整的分析。

注意一若是穿透操纵变成整体平稳迭代中的主因（若是为使问题收敛到穿

透容差内，比收敛到不平稳力的容差内，需要更多的迭代），用户应该增大FTOLN

值，以许诺更多的穿透，或增大FKNo

选择摩擦类型

在大体的库仑摩擦模型中，两个接触面在开始彼此滑动之前，在它们的界面

上会有达到某一大小的剪应力产生。这种状态称为粘合状态（stick）。库仑摩擦

模型概念了一个等效剪应力T,在某一法向压应力p作用下剪应力达到此值时

表面开始滑动（T=UP+C0HE,其中u是摩擦系数一MU-作为材料特性概念，而

C0HE是粘聚力）。一旦剪应力超过此值后，两个表面之间将开始彼此滑动。这种

状态，叫作滑动状态（Sliding）。粘合/滑动计算决定何时一个点从粘合状态到滑

动状态，或从滑动状态变到粘合状态。摩擦系数能够是任一非负值，程序缺省值

为表面之间无摩擦。关于粗糙或绑定接触(KEY0PT(12)二一、3、五、6),程序将

不管给定的MU值而以为摩擦阻力无穷大。

程序提供了一个人为指定最大等效剪应力的选项，不管接触压力值的大小，

若是等效剪应力达到此值时，即发生滑动，见图5-9o为了指定接触界面上最

大允许剪应力，设置常数TAUMAX(缺省为。那个剪应力极限，通经常使用于在

接触压力超级大的时候(如在某些加工进程中)的一些情形，以至于用库仑理论计

算出的界面剪应力超过了材料的屈服极限。TAUMAX的一个合理上限估值为

。/石是表面周围材料的vonMises屈服应力)。体会数据有助于决定TAUMAX

的值。

I«I

TAUMAX-----/--------------------------

COME/

Zaror♦

图5-9摩擦模式

静、动摩擦系数

摩擦系数依托于接触面的相对滑动速度，通常静摩擦系数高于动摩擦系数。

ANSYS提供了如下表示的指数衰减摩擦模型：

U=MUX(1+(FACT-1)exp(-DCXVrel)

其中：

P为摩擦系数。

MU动摩擦系数，用MP命令输入。

FACT是静摩擦系数与动摩擦系数之比，缺省为最小值。

DC为衰减系数，缺省为，单位为time/length。因此，时刻在静态分

析中有一些意义。

V.1是ANS'S计算的滑动速度。

若是明白静、动摩擦系数和至少一个数据点（R-Vrol）,那么能够确信

摩擦衰减系数为：

-t（FACT-1）xMUj

若是不指定衰减系数，且FACT大于，当接触进入滑动状态时，摩擦

系数会从静摩擦系数突变到动摩擦系数，这种行为类似于C0NTAC46和C0NTAC49

单元所用的动摩擦模型，因为这会致使收敛困难，因此不建议采纳。

对称、不对称求解器

对无摩擦、粗糙和绑定接触，接触单元刚度矩阵是对称的。而涉及到摩擦的

接触问题产生一个不对称的刚度。在每次迭代利用不对称的求解器，比对称的求

解器需要更多的计算时刻。因此ANSYS程序采纳对称化算法。通过采纳这种算法

大多数摩擦接触问题，能够利用对称系统的求解器来求解。若是摩擦应力在整个

位移场内有相当大的阻碍，而且摩擦应力的大小高度依托于求解进程，那么对刚

度阵的任何对称近似都可能致使收敛性降低。在这种情形下，选择不对称求解选

项（NROPT,UNSYM）来改善收敛性。

选择接触检查的位置

接触检查点位于接触单元的积分点上。在积分点上，接触单元不穿透进入目

标面。但是，目标面能穿透进入接触而。见图5-10o

G*UMirvit9raiionpom

图5-10接触检查点位于高斯积分点上

图5T1接触检查点位于节点上

ANSYS面一面接触单元利用高斯积分点作为缺省值，高斯积分点通常会比

Newton-Cotes/Lobatto节点积分方案产生更精准的结果，

Newton-cotes/Lobatto用节点本身作为积分点。通过KEY0PT(4)来选择用户想

利用的方式。这一选项仅适用于低阶接触(C0NTAC171和CONTAC173)o但是，

利用节点本身作为积分点仅应该用于角接触问题;看图5T1)o

注意，利用节点作为接触检查点可能会致使其它收敛性问题，例如“滑

脱”(节点滑出目标面的边界)，见图5-12o对大多数的点一面的接触问题，

咱们推荐利用其它的点一面的接触单元，例如C0NTAC2六、C0NTAC48和

CONTAC49。见本书§o

图5-12节点滑脱

调整初始接触条件

在动态分析中，刚体运动一样可不能引发问题。但是在静力分析中，当物体

没有足够的约束时会产生刚体运动，有可能引发错误而终止计算。

在仅仅通过接触的显现来约束刚体运动时，必需保证在初始几何体中，接触

对是接触的。换句话说，用户要成立模型以便接触对是“恰好接触”的。但是如

此做，可能会碰到以下问题：

刚体外形常常是复杂的，很难决定第一个接触点发生在哪儿。

既使实体模型在初始时处于接触状态，在网格划分后由于数值舍入误

差，两个面的单元网格之间也可能会产生小裂缝。

接触单元的积分点和目标单元之间可能有小裂缝。

同理，在目标面和接触面之间可能发生过大的初始穿透。在这种情形下，接

触单元可能会高估接触力.致使不收敛或接触面之间脱开接触关系.

概念初始接触或许是成立接触分析模型时最亘要的方面。因此，程序提供了

几种方式来调整接触对的初始接触条件。

注意：下面的技术能够在开始分析时独立执行，或几个联合起来执行。它

们是为了排除由于生成网格造成的数值舍入误差而引发的小间隙或穿透，而不是

为了更正网格或儿何数据的错误。

一、应用实常数CNOF来指定一个接触面偏移。

指定正的值来使整个接触面偏向目标面。指定负的值来使接触面离开目标

面。

ANSYS能够自动提供CNOF值到恰好闭合间隙或减少初始穿透。如下设置

KEY0PT（5）：

=1：闭合间隙；

二2:减少初始穿透；

二3:闭合间隙或减少初始穿透。

若是设置了KEYOPT（5）＞0,那么TC0NT缺省值为0。

二、利用实常数TCONT来指定一个小的初始接触环，初始接触环是指沿着

目标面的“调整环”的深度。若是没有人为指定IC0NT的值，程序会依照几何

尺寸来给IC0NT提供一个小值（但成心义的值），同时输出一个表示什么值被指

定的警告信息。IC0M正值表示相关于基层单元厚度的比例因子；负值表示接

触环的绝对值。任何落在“调整环”域内的接触检查点被自动移到目标面上，（参

见图5T3a）。建议利用一个十分小的IC0NT值，不然可能会发生严峻不持续

（看图5-13b）

CortKt

(b)

图5T3用ICON进行接触面的调整。3）调整前；（b）调整后

CNOF与TCONT之间的不同，是前者把整个接触面移动CNOF的距离，而后

者把所有初始分开的（恰好位于调整环ICONT内的）接触点向目标面移动。

若是用户应用其他方式来平稳初始未约束的自由体（如FTOLN,PINB,PMAX

和PMIN）,大体上排除ICONT的阻碍（把它关联一个很小的值，如1E-20）,是一

个好方法。可是，设置IC0NTR,并非会关闭它。而是致使ANSYS用缺省值代

替。当与其他约束自由体的方式联历时，这反过来可能不能达到成效。

3、利用实常数PMIN和PMAX来指定初始允许的穿透范围。当指定PMAX或

PMIN后，在开始分析时，程序会将目标面移到初始接触状态，见图5T4。若

是初始穿透大于PMAX,程序会调整目标面来减少穿透。接触状态的初始调剂仅

仅通过平移来实现。

对给定载荷或给定位移的刚性目标面，将会执行初始接触状态的初始调剂。

对没有指定边界条件的目标面，也一样能够进行初始接触的调整。

当目标面上的所有节点,有给定的零位移值时,利用PMAX和PMIN的初始

调剂将可不能被执行。

图5T4接触面调整(PMIN,PMAX)

注意--ANSYS程序独立地处置目标面上节点的自由度。例如：若是用户指定

自由度UX值为“0”，那么沿着X方向就没有初始调查。但是，在Y和Z方

向仍然会激活PMAX和PMIN选项。

初始状态调整是一个迭代进程，程序最多进行20次迭代。若是目标面K能

进入可同意的穿透范围(即PMIN,PMAX范围)，程序将在原始几何实体上操作。

这时程序会给出一个警告信息，用户可能需要调整用户的初始几何模型。

图5-15给出了一个初始接触调整迭代失败的例子。目标面的UY被约束。

因此，初始接触唯一允许的调整是在X方向，但是，在那个问题中，刚性目标

面在X方向的任何运动都可不能引发初始接触。

关于柔体-柔体接触，这种方式不仅移动整个目标面，还同时移动与目标面

相连的整个柔体。请确保没有其他接触面或目标面与柔体相连。

图5-15一个初始调整失败的例子

4、设置KEY0PT⑼力来调整初始穿透或间隙，见图5-16。

图5T6忽略初始穿透,KEY0PT(9)=l

真正的初始穿透包括两部份：

几何模型产生的穿透或间隙；

用户概念的接触面偏移(CNOF)产生的穿透或间隙；

KEYOPT⑼提供以下功能：

包括由儿何和接触面偏移产生的初始穿透，设置KEYOPT(9)=0o这是

缺省。

忽略上面二者引发的初始穿透,设置KEYOPT(9)=1。在KEY0PT(12)=4

或5时，这一KEYOPT(9)=1,也将忽略间隙弹簧的初始力，如此，成立

了一个初始的“理想的”接触面一在接触截面上没有初始力的作用。

为了包括概念的接触面偏移(CNOF)，但忽略由于儿何模型引发的初始

穿透，设置KEYOPT(9)=3。在KEYOPT(12)=4或5时，这一KEYOPT⑼=3,

也将忽略打开间隙弹簧的初始力，如此，成立了一个初始的“理想的”

接触面一在接触截面上没有初始力的作用。

在某些情形下，例如过盈装配问题，期望有过度的穿透。如在第一个载荷步

施加阶跃初始穿透，可能造成收敛困难。为了减缓收敛性困难，在第一个载荷步

中设置渐变的初始穿透来克服。见图5-17o下面的KEY0PN9)设置用来泥供

渐变功能：

设置KEY0PT(9)=2,来渐变施加初始穿透(CNOF+由于几何模型造成的偏

移)。

设置KEY0PT(9)=4,来渐变施加接触面穿透，但忽略由于几何模型造成

的穿透。

关于上面两个KEY0PT(9)设置，用户还应该设置KBC,0,并在第一个载荷

步中不要给定任何其它外载荷。还要确保球形(Pinball)区域足够大以捕五到

初始过盈。

用户能够联合应用上面的技术。例如，用户可能希望设置十分精准的初始渗

透或间隙，但有限元节点的初始坐标可能无法提供足够的精度。这时，能够：

应用ICONT来移动初始张开的接触点恰好碰着目标面。

应用CNOF来指定穿透(正值)或间隙(负值)。

应用KEY0PT(9)=3来在第一个子步求解初始穿透，或应用

KEYOPT(9)=4来慢慢求解初始穿透。

Cortxt

StrUa

・/

EMM

图5-17渐变初始过盈

在开始分析时，程序会给出每一个目标面的初始接触状态的输出信息（输出

窗口或输出文件中），那个信息有助于决定每一个目标面的最大穿透或最小间隙。

关于给定的目标面，若是没有发觉接触，可能是目标面离接触面太远（超出

了球形区域），或是接触/目标单元已经被杀死。

决定接触状态和球形区域。

接触单元相关于其目标面的运动和位置，决定了接触单元的状态；程序检测

每一个接触单元，并给出一种状态：

STAT二0未闭合的远场接触

STAT=1未闭合的近场接触

STAT=2滑动接触

STAT=3粘合接触

当目标面进入球形区域后，接触单元就被看成未闭合的近场接触，球形区域

是以接触单元的积分点为中心的。利用实常数PINB来为球形区域指定一个比例

因子（正值），或其绝对值（负值）。缺省时，程序将球形区域概念为一个以“4*

基层单元厚度（关于刚体-柔体接触）”或“2*基层单元厚度（关于柔体-柔体接

触）”为半径的圆（对2-D问题）或球（对3-D问题

检查接触的计算时刻依托于球形区域的大小，远场接触单元计算简单，计算

时刻较少。近场接触计算（关于接近接触或实际接触的接触单元）较慢而且较复

杂。当单元己经接触时，计算最为复杂。

若是刚性面有好几个凸形区域，为了克服伪接触概念，设置一个适合的球形

区域是有效的。而对大多数问题，缺省值是适合的。

在自接触问题中幸免伪接触

在一些对称接触问题（包括自接触）中，ANSYS可能错误地假设在十分接近的

儿何位置上的接触面和目标面之间的接触。在自接触问题的角点会发生这一问

题。它可能是由单元的初始几何位置引发的，也可能在分析时通过变形而引发的。

在二个面位于球形区域内，而且它们之间夹角小于90度时，会产生这一问题。

在这种情形下,ANSYS程序以为发生了十分大的穿透。图5-18说明了这种情形。

rm

canueauvKhfrtorcintafoecvTiftnr

poknncrdjereur•

iwnerKSvOfe”LtucnsfUa

Rivers

图5-18检测伪接触

用户能够通过单元C0NTA171-C0NTA174的KE0PT(8)=l,来避免ANSYS程序

考虑这一问题。当应生那个KEYOPT时，ANSYS将忽略在如下情形下产生的“过

度穿透”接触：

初始检测到的穿透大于接触容差(FTOLN)的20机见图5-18a。

接触状态由图5-18a所示的远离接触突然转变到图5-18b所示的过

度穿透。

若是KEY0PT(8)激活，ANSYS各个荷步中第一次检测到伪接触时会发出一个

警告。如ANSYS在第1荷载步中发觉这种接触，用户可看到如下的信息：

aContactelementxhastoomuchpenetrationrelatedtotargetelement

y-

Weassumeit(maybemoreelements)isspuriouscontact.”

若是ANSYS检测到归类为伪接触的突变，用户看到如下信息：

“Contactelementxstatuschangedabruptlywithtargetelementy.

Weassumeit(maybemoreelements)isspuriouscontact.”

ANSYS在一个荷载步中仅发出一次如此的信息。在该荷载步中若是还存在其

他的伪接触，ANSYS再也不提示。

注意一在应用KEY0PT(8)=

l之前，请认真检查模型。ANSYS将忽略符合了过

度穿透准则的真正接触。这•技术应当仅用于在改变球形区域PINB不能阻止

伪接触的情形。

修正法向和切向接触刚度

在分析进程中，能够修正法向和切向接触刚度。能够自动修正(由于改变下

伏单元刚度的大应变效应产生)，也能够显式地修正(由用户从头指定FKN或

FKT值)。KEYOPT(IO)操纵法向和切向接触刚度如何修正：

KEYOPT(10)=0,禁止那些已经处于“闭合”状态的单元的接触刚度修

正。关于从“张开”转变到“闭合”状态的单元，将在每一个子步上修

正接触刚度。

KEYOPT(1O)=1,许诺已处于“闭合”状态的单元的接触刚度，在荷载

步之间或在重启动期间改变。关于从“张开”转变到“闭合”状态的单

元，将在每一个子步上修正接触刚度。

KEY0PT(10)=2与KEYOPT(10)=1相同,只是对所有单元(不论其状态)

将在每一个子步上，由程序决定自动修正。

选择表面作用模式

面一面接触单元支持法向单向接触模式及其他力学表面作用模式。通过设置

KEY0PT(12)来选择下面的某种作用模式：

KEYOPT(12)=0,法向单向接触，即在接触分开时，法向压力等于0。

KEYOPT(12)=1,理想粗糙接触，用来模拟无滑动的，表面完全粗糙的

摩擦接触问题，这种设置对应于摩擦系数无穷大，因此用户概念的摩擦

系数(MU)被忽略。

KEY0PT(12)=2,“不分离”接触,接触面和目标面一日接触,在其后

的分析中就连在一路(尽管许诺有相对滑动)。

KEY0PT(12)=3,绑定接触模式，目标面和接触面一旦接触，随后就在

所有方向上绑定。

KEY0PT（12）=4,不分离接触，其中的接触积分点，或初始在球形区域

内，或一旦接触，就老是与目标面沿接触面的法向连在一路，但许诺滑

动。调整FKOP（见下），可用“软弹簧”把这些区域联系在一路。

KEY0PT（12）=5,绑定接触模式，其中接触积分点，或初始在球形区域

内，或一旦接触，就老是与目标面沿接触面的法向和切向绑定在一路。

KEY0PT（12）=6,绑定接触模式，其中初始接触的接触积分点维持与目

标面接触，而初始处于打开状态的接触积分点，在整个分析期间维持打

开状态。那个选项与在初始接触的区域应用CEINTF类似。

关于模拟不分离或绑定接触，用户可能需要设置FKOP实常数。这在接触张

开时，提供一个刚度系数。若是FKOP为正值，那么真正的接触张开刚度等于

FKOP乘以接触闭合时施加的刚度。若是FKOP为负值，该值作为接触张开刚度

的绝对值。缺省的FKOP值为1。

不分离或绑定接触，在接触发生张开时，产生〃回拉〃力，那个力可能不足以

阻止分离。为了减小分离，概念一个较大的FKOP值。在有些时候，希望接触面

分离，但需要在接触面之间成立联系来阻止刚体运动，在这种情形下，能够指定

较小的FKOP值，以使接触面之间维持联系（这是〃软弹簧〃效应）。

用超单元成立接触模型

面一面的接触单元能模拟刚体（或一个线弹性体）和另一个有小位移的线弹

性体的接触。这些线弹性体可用超单元来建模，这大大降低了进行接触迭代的自

由度数。记住任何接触或目标节点都必需是超单元的主自由度。

由于超单元仅仅由一组保留的节点自由度组成，它没有效来概念接触面和目

标面的几何形状。因此，必需在形成超单元之前在原始单元表面上概念接触面和

目标面。来自超单元的信息，包括节点连结和组合刚度，可是没有材料特性和应

力状态（轴对称、平面应力或平面应变）。一个限制是接触单元的材料特性设置必

需与形成超单元之前的原始单元的材料特性相同。

利用KEY0PT(3)来提供2D接触分析的信息。C0NTA17—、C0NTA172单元的

选项如下：

不利用超单元(KET0PT(3)=0)o

轴对称(KEYOPT(3)=1)。

平面应变或单位厚度的平面应力(KEY0PT(3)=2)。

需要厚度输入的平面应力(KEY0PT(3)=3)。对这种情形，利用实常数R2

来指定厚度。

关于3D接触分析,C0NTA173,CONTA174单元的KEY0PT(3)选项忽略。ANSYS

将自动检查下伏单元是不是为超单元0

考虑厚度阻碍

程序能够用KEY0PX11)来考虑壳(2D和3D)、梁(2D)的厚度。关于刚体-柔体

接触，ANSYS将自动移动接触面到壳/梁的底面或顶面。关于柔体-柔体接触，

ANSYS将自动移动与壳/梁单元相连的接触面和目标面。缺省时，程序不考虑单

元厚度，用中面来表示梁和壳，而穿透距离从中面计算。当设置KFT0P1(1.)=1

时，那么考虑梁或壳的厚度。从指定的底面或顶面来计算穿透距离。

注意一仅在利用节点位于中面的壳或梁单元时(例如，KEYOPT(11)=O的

SHELL91单元)，用KEYOPT(11)=1来考虑厚度阻碍。

建模时如要考虑厚度，记住偏移可能来自接触面或目标面或二者。在和

KEYOPT(11)=1一路指定接触偏移(CNOF)时，CNOF从壳/梁的顶面和底面计算，而

不是中面。当和SHELL181一路利历时,还考虑变形进程中的厚度转变。

利历时刻步长操纵

时刻步长操纵是一个自动时刻步长特点，那个