ANSYS非线形分析指南 接触分析 Rl和R2定义目标单元几何形状 FKN定义法向接触刚度因子 FTOLN定义最大的渗透范围 ICONT定义初始靠近因子 PINB定义“ Pinball"区域 PMN和PMAⅹ定义初始渗透的容许范围 TAUMAR指定最大的接触摩擦 命令:R GUI: main menu> preprocessor>real constant 对实常数FKN, FTOLN, ICONT,PINB,PMAX,和PMN,你既可以定 义一个正值也可以定义一个负值,程序将正值作为比例因子,将负值作为真实值, 程序将下面覆盖原单元的厚度作为ICON, FTOLN,PINB,PMAX和PMN的 参考值,例如对ICON,0.1表明初始间隙因子是0.1*下面覆盖层单元的厚度 然而,0.1表明真实缝隙是0.1,如果下面覆盖层单元是超单元,则将接触单元 的最小长度作为厚度。 单元关键字 每种接触单元都饭知好几个关键字,对大多的接触问题缺省的关键字是合适 的,而在某些情况下,可能需要改变缺省值,来控制接触行为 接触算法(罚函数+拉格郎日或罚函数)( KEYOPT(2)) 出现超单元时的应力状态( DEYOPT(3)) 接触方位点的位置 (KEYOPI (4)) 刚度矩阵的选择 (KEYOPT (6)) 时间步长控制 (KEYOPT (7)) 初始渗透影响 (KEYOPT(9)) 接触表面情况 (KEYOPT(12)) 命令: KEYOPT ET GUI: main menu>preprocessor>Elemant Type>Add/Edit/Delete 选择接触算法 对面一面的接触单元,程序可以使用扩增的拉格朗日算法或罚函数方法,通 第11页
ANSYS 非线形分析指南 接触分析 第11页 R1 和 R2 定义目标单元几何形状 FKN 定义法向接触刚度因子 FTOLN 定义最大的渗透范围 ICONT 定义初始靠近因子 PINB 定义“Pinball"区域 PMIN 和 PMAX 定义初始渗透的容许范围 TAUMAR 指定最大的接触摩擦 命令：R GUI：main menu> preprocessor>real constant 对实常数 FKN，FTOLN，ICONT，PINB，PMAX，和 PMIN，你既可以定 义一个正值也可以定义一个负值，程序将正值作为比例因子，将负值作为真实值， 程序将下面覆盖原单元的厚度作为 ICON，FTOLN，PINB，PMAX 和 PMIN 的 参考值，例如对 ICON，0.1 表明初始间隙因子是 0.1*下面覆盖层单元的厚度。 然而，-0.1 表明真实缝隙是 0.1，如果下面覆盖层单元是超单元，则将接触单元 的最小长度作为厚度。 单元关键字 每种接触单元都饭知好几个关键字，对大多的接触问题缺省的关键字是合适 的，而在某些情况下，可能需要改变缺省值，来控制接触行为。 接触算法（罚函数+拉格郎日或罚函数）（KEYOPT（2）） 出现超单元时的应力状态（DEYOPT（3）） 接触方位点的位置 （KEYOPI（4）） 刚度矩阵的选择 （KEYOPT（6）） 时间步长控制 （KEYOPT（7）） 初始渗透影响 （KEYOPT（9）） 接触表面情况 （KEYOPT（12）） 命令：KEYOPT ET GUI：main menu>preprocessor>Elemant Type>Add/Edit/Delete 选择接触算法： 对面─面的接触单元，程序可以使用扩增的拉格朗日算法或罚函数方法，通
ANSYS非线形分析指南 接触分析 过使用单元关键字 KETOPT(2)来指定。 扩张的拉格朗日算法是为了找到精确的拉格朗日乘子而对罚函数修正项进 行反复迭代,与罚函数的方法相比,拉格朗日方法不易引起病态条件,对接触刚 度的灵敏度较小,然而,在有些分析中,扩增的拉格朗日方法可能需要更多的迭 代,特别是在变形后网格变得太扭曲时。 使用拉格朗日算法的同时应使用实常数 FTOLN FTOLN为搠格朗日算法指定容许的最大渗艉,如果程序发现渗透大于此值 时,即使不平衡力和位移增量已经满足了收敛准则,总的求解仍被当作不收敛处 理, FTLON的缺省值为0.1,你可以改变这个值,但要注意如果此值太小可能会 造成太多的迭代次数或者不收敛 决定接触刚度 所有的接触问题都需要定义接触刚度,两个表面之间渗量的大小取决了接 触刚度,过大的接触刚度可能会引起总刚矩阵的病态,而造成收敛困难,一般来 谘,应该选取足够大的接触刚度以保证接触渗透小到可以接受,但同时又应该让 接触刚度足够小以使不会引起总刚矩阵的病态问题而保证收敛性。 程序会根据变形体单元的材料特性来估计一个缺省的接触刚度值,你能够用 实常数FKN来为接触刚度指定一个比例因子或指定一个真正的值,比例因子 般在001和10之间,当避免过多的迭代次数时,应该尽量使渗透到达极小值 为了取得一个较好的接触刚度值,又可需要一些经验,你可以按下面的步骤 1、开始时取一个较低的值,低估些值要比高估些值好因为由一个较低的接 触刚度导致的渗透问题要比过高的接触刚度导致的收敛性困难,要容 易解决。 2、对前几个子步进行计算 3、检查渗透量和每一子步中的平衡迭代次数,如果总体收敛困难是由过大 的渗透引起的(而不是由不平衡力和位移增量引起的),那么可能低估 了FKN的值或者是将 FTOLN的值取得大小,如果总体的收敛困难是由 于不平衡力和位移增量达到收敛值需要过多的迭代次数,而不是由于过 大的渗透量,那么FKN的值可能被高估 4、按需要调查FKN或 FTOLN的值,重新分析。 第12页
ANSYS 非线形分析指南 接触分析 第12页 过使用单元关键字 KETOPT（2）来指定。 扩张的拉格朗日算法是为了找到精确的拉格朗日乘子而对罚函数修正项进 行反复迭代，与罚函数的方法相比，拉格朗日方法不易引起病态条件，对接触刚 度的灵敏度较小，然而，在有些分析中，扩增的拉格朗日方法可能需要更多的迭 代，特别是在变形后网格变得太扭曲时。 使用拉格朗日算法的同时应使用实常数 FTOLN FTOLN 为搠格朗日算法指定容许的最大渗艉，如果程序发现渗透大于此值 时，即使不平衡力和位移增量已经满足了收敛准则，总的求解仍被当作不收敛处 理，FTLON 的缺省值为 0.1，你可以改变这个值，但要注意如果此值太小可能会 造成太多的迭代次数或者不收敛。 决定接触刚度 所有的接触问题都需要定义接触刚度，两个表面之间渗 量的大小取决了接 触刚度，过大的接触刚度可能会引起总刚矩阵的病态，而造成收敛困难，一般来 谘，应该选取足够大的接触刚度以保证接触渗透小到可以接受，但同时又应该让 接触刚度足够小以使不会引起总刚矩阵的病态问题而保证收敛性。 程序会根据变形体单元的材料特性来估计一个缺省的接触刚度值，你能够用 实常数 FKN 来为接触刚度指定一个比例因子或指定一个真正的值，比例因子一 般在 0.01 和 10 之间，当避免过多的迭代次数时，应该尽量使渗透到达极小值。 为了取得一个较好的接触刚度值，又可需要一些经验，你可以按下面的步骤 过行。 1、 开始时取一个较低的值，低估些值要比高估些值好因为由一个较低的接 触刚度导致的渗透问题要比 过高的接触刚度导致的收敛性困难，要容 易解决。 2、 对前几个子步进行计算 3、 检查渗透量和每一子步中的平衡迭代次数，如果总体收敛困难是由过大 的渗透引起的（而不是由不平衡力和位移增量引起的），那么可能低估 了 FKN 的值或者是将 FTOLN 的值取得大小，如果总体的收敛困难是由 于不平衡力和位移增量达到收敛值需要过多的迭代次数，而不是由于过 大的渗透量，那么 FKN 的值可能被高估。 4、 按需要调查 FKN 或 FTOLN 的值，重新分析
ANSYS非线形分析指南 接触分析 选择摩擦类型。 在基本的库仑摩擦模型中,两个接触面在开始相互滑动之前,在它们的界面 上会有达到某一大小的剪应力产生,这种状态则作粘合状态( stick)库仑摩擦模 型定义了一个等效剪应力。)一旦剪应力超过此值后,两个表面之间将开始相互 滑动,这种状态,叫作滑动状态( Sliding)粘合滑动计算决定什么时候一个点 从粘合状态到滑动状态或从滑动状态变到粘合状态,摩擦系数可以是任一非负 值。程序缺省值为表面之间无摩擦,对 rough或 bonded接触( KEYOPT(2)=1 (或3),程序将不管给定的MV值而认为摩擦阻力无限大。 程序提供了一个不管接触压力的故而人为指定最大等效剪应力的选项,如果 等效剪应力达到此值时,滑动发生。看图,4-1,为了指定接触界面上最大许可 剪应力,设置常数 TAUMAX(缺省为10E20),这种限制剪应力的情况一般用 于接触压力非常大的时候,以至于用库仑理论计算出的界面剪应力超过了材料的 屈服极限。一对 TAUMAX的一个合理高估为(一是材料的 mises屈服应 力) Sliding Sticker 图4-1摩擦模式 对无摩擦 rough和 bonded接触,接触单元刚度矩阵是对称的,而涉及到摩 擦的接触问题产生一个不对称的刚度,而在每次迭代使用不对称的求解器比对称 的求解器需要更多的计算时间,因此 ANSYS程序采用对称化算法。通过采用这 种算法大多的摩擦接触问题能够使用对称系统的求解器来求解。如果摩擦应力在 整个位移范围内有相当大的影响,并且摩擦应力的大小高度依赖于求解过程。对 刚度阵的任何对称近似都可能导致收敛性的降低,在这种情况下,选择不对称求 第13页
ANSYS 非线形分析指南 接触分析 第13页 选择摩擦类型。 在基本的库仑摩擦模型中，两个接触面在开始相互滑动之前，在它们的界面 上会有达到某一大小的剪应力产生，这种状态则作粘合状态（stick)库仑摩擦模 型定义了一个等效剪应力。）一旦剪应力超过此值后，两个表面之间将开始相互 滑动，这种状态，叫作滑动状态（Sliding）粘合\滑动计算决定什么时候一个点 从粘合状态到滑动状态或从滑动状态变到粘合状态，摩擦系数可以是任一非负 值。程序缺省值为表面之间无摩擦，对 rough 或 bonded 接触（KEYOPT（2）=1 （或 3），程序将不管给定的 MV 值而认为摩擦阻力无限大。 程序提供了一个不管接触压力的故而人为指定最大等效剪应力的选项，如果 等效剪应力达到此值时，滑动发生。看图，4─1，为了指定接触界面上最大许可 剪应力，设置常数 TAUMAX（缺省为 1.0E20），这种限制剪应力的情况一般用 于接触压力非常大的时候，以至于用库仑理论计算出的界面剪应力超过了材料的 屈服极限。一对 TAUMAX 的一个合理高估为 y 3 （ y 3 是材料的 mises 屈服应 力）。 图 4—1 摩擦模式 对无摩擦 rough 和 bonded 接触，接触单元刚度矩阵是对称的，而涉及到摩 擦的接触问题产生一个不对称的刚度，而在每次迭代使用不对称的求解器比对称 的求解器需要更多的计算时间，因此 ANSYS 程序采用对称化算法。通过采用这 种算法大多的摩擦接触问题能够使用对称系统的求解器来求解。如果摩擦应力在 整个位移范围内有相当大的影响，并且摩擦应力的大小高度依赖于求解过程。对 刚度阵的任何对称近似都可能导致收敛性的降低，在这种情况下，选择不对称求
ANSYS非线形分析指南 接触分析 解选项( KEYOPT(6)=1)来改善收敛性。 选择检查接触与否的位置 接触检查点位于接触单元的积分点上,在积分点上,接触单元不渗透进入目 标面,然而,目标面能渗透进入接触面,看图4-2。 Def ble solid Gauss integration point Contact segment Target s 图4—2接触检查点位于高斯积分点上 ANSYS面一接触单元使用 GAUSS积分点作为缺省值, GAUSS积分点通常 会比 Newton- Cotes/robatto结点积分项产生更精确的结果, Newton- cotes/, lobatto 使用结点本身作为积分点,通过 KEYOPT(4)来选择,你想使用的方法,然而 使用结点本身作为积分点仅应该用于角接触问题(看图4-3)。 Deformed bod Gauss integration point 图4-3接触检查点位于高斯结点上 然而,使用结点作为接触发现点,可能会导致其它的收敛性问题,例如“滑 脱”(结点滑下目标面的边界)看图4-4,对大多的点一面的接触问题,我们 第14页
ANSYS 非线形分析指南 接触分析 第14页 解选项（KEYOPT（6）=1）来改善收敛性。 选择检查接触与否的位置 接触检查点位于接触单元的积分点上，在积分点上，接触单元不渗透进入目 标面，然而，目标面能渗透进入接触面，看图 4—2。 图 4—2 接触检查点位于高斯积分点上 ANSYS 面─接触单元使用 GAUSS 积分点作为缺省值，GAUSS 积分点通常 会比 Newton-Cotes/robatto 结点积分项产生更精确的结果，Newton-cotes/lobatto 使用结点本身作为积分点，通过 KEYOPT（4）来选择，你想使用的方法，然而， 使用结点本身作为积分点仅应该用于角接触问题（看图 4─3）。 图 4—3 接触检查点位于高斯结点上 然而，使用结点作为接触发现点，可能会导致其它的收敛性问题，例如“滑 脱”（结点滑下目标面的边界）看图 4─4，对大多的点─面的接触问题，我们
ANSYS非线形分析指南 接触分析 推荐使用其它的点一面的接触单元,例如 CONTA26、 CONTA48和 CONTA49。 图4-4结点滑脱 调整初始接触条件 在动态分析中,刚体运动一般不会引起问题,然而在静力分析中,当物体没 有足够的约束时会产生刚体运动,有可能引起错误而终止计算 在仅仅通过接触的出现来约束刚体运动时,必须保证在初始几何体中,接触 对是接触的,换句话说,你要建立模型以便接触对是“刚好接触”的,然而这样 作可能会遇到以下问题 刚体外形常常是复杂的,很难决定第一个接触点发生在哪儿 既使实体模型是在初始接触状态,在网格划分后余于数值舍入误差;两个 面的单元网格之间也可能会产生小的缝隙。 接触单元的积分点和目标单元之间可能有小的缝隙 同理,在目标面和接触面之间可能发生过大的初始渗透,在这种情况下,接 触单元可能会高估接触力,导致不收敛或民接触面之间脱离开接触关系。定义 初始接触也许是建立接触分析模型时最重要的方面,因此,程序提供了几种方法 来调整接触对的初始接触条件 注意:下面的技巧可以在开始分析时独立执行成几个联合起来执行,它们是 为了消除由于生成网格造成的数值舍入误差而引起的小缝隙或渗透,而不是为了 改正网格或几何数据的错误。 1、使用实常数 ICONT来指定一个好的初始接触环,初始接触环是指沿着 目标面的“调整环”的深度,如果没有人为指定 ICONT的值,程序会 根据几何尺寸来给 ICONT提供一个小值,同时输出一个表时什么值被 第15页
ANSYS 非线形分析指南 接触分析 第15页 推荐使用其它的点─面的接触单元，例如 CONTA26、CONTA48 和 CONTA49。 图 4—4 结点滑脱 调整初始接触条件 在动态分析中，刚体运动一般不会引起问题，然而在静力分析中，当物体没 有足够的约束时会产生刚体运动，有可能引起错误而终止计算。 在仅仅通过接触的出现来约束刚体运动时，必须保证在初始几何体中，接触 对是接触的，换句话说，你要建立模型以便接触对是“刚好接触”的，然而这样 作可能会遇到以下问题： ·刚体外形常常是复杂的，很难决定第一个接触点发生在哪儿 ·既使实体模型是在初始接触状态，在网格划分后余于数值舍入误差；两个 面的单元网格之间也可能会产生小的缝隙。 ·接触单元的积分点和目标单元之间可能有小的缝隙。 同理，在目标面和接触面之间可能发生过大的初始渗透，在这种情况下，接 触单元可能会高估接触力，导致不收敛或 民接触面之间脱离开接触关系。定义 初始接触也许是建立接触分析模型时最重要的方面，因此，程序提供了几种方法 来调整接触对的初始接触条件。 注意：下面的技巧可以在开始分析时独立执行成几个联合起来执行，它们是 为了消除由于生成网格造成的数值舍入误差而引起的小缝隙或渗透，而不是为了 改正网格或几何数据的错误。 1、 使用实常数 ICONT 来指定一个好的初始接触环，初始接触环是指沿着 目标面的“调整环”的深度，如果没有人为指定 ICONT 的值，程序会 根据几何尺寸来给 ICONT 提供一个小值，同时输出一个表时什么值被