下载文档原格式
接触问题是一种高度非线性行为，需要较大的计算资源，为了进行实为有效的计算，理解问题的特性和建立合理的模型是很重要的。
接触问题存在两个较大的难点：其一，在你求解问题之前，你不知道接触区域，表面之间是接触或分开是未知的，突然变化的，这随载荷、材料、边界条件和其它因素而定；其二，大多的接触问题需要计算摩擦，有几种摩擦和模型供你挑选，它们都是非线性的，摩擦使问题的收敛性变得困难。
一般的接触分类接触问题分为两种基本类型：刚体─柔体的接触，半柔体─柔体的接触，在刚体─柔体的接触问题中，接触面的一个或多个被当作刚体，（与它接触的变形体相比，有大得多的刚度），一般情况下，一种软材料和一种硬材料接触时，问题可以被假定为刚体─柔体的接触，许多金属成形问题归为此类接触，另一类，柔体─柔体的接触，是一种更普遍的类型，在这种情况下，两个接触体都是变形体（有近似的刚度）。
ANSYS接触能力ANSYS支持三种接触方式：点─点，点─面，平面─面，每种接触方式使用的接触单元适用于某类问题。
为了给接触问题建模，首先必须认识到模型中的哪些部分可能会相互接触，如果相互作用的其中之一是一点，模型的对立应组元是一个结点。
如果相互作用的其中之一是一个面，模型的对应组元是单元，例如梁单元，壳单元或实体单元，有限元模型通过指定的接触单元来识别可能的接触匹对，接触单元是覆盖在分析模型接触面之上的一层单元，至于ANSTS使用的接触单元和使用它们的过程，下面分类详述。
点─点接触单元点─点接触单元主要用于模拟点─点的接触行为，为了使用点─点的接触单元，你需要预先知道接触位置，这类接触问题只能适用于接触面之间有较小相对滑动的情况（即使在几何非线性情况下）如果两个面上的结点一一对应，相对滑动又以忽略不计，两个面挠度（转动）保持小量，那么可以用点─点的接触单元来求解面─面的接触问题，过盈装配问题是一个用点─点的接触单元来模拟面─与的接触问题的典型例子。
点─面接触单元点─面接触单元主要用于给点─面的接触行为建模，例如两根梁的相互接触。
基于ANSYS软件的接触问题分析及在工程中的应用基于ANSYS软件的接触问题分析及在工程中的应用一、引言接触问题是工程领域中常见的一个重要问题，它在很多实际应用中都具有关键作用。
接触分析能够帮助工程师设计和改进各种产品和结构，从而提高其性能和寿命，减少故障和事故的发生。
ANSYS作为一款强大的工程仿真软件，提供了多种接触分析方法和工具，为工程师们解决接触问题提供了便利。
本文将重点介绍基于ANSYS软件的接触问题分析方法和其在工程中的应用。
二、接触问题的分析方法接触问题的分析方法主要包括两种：解析方法和数值模拟方法。
解析方法基于一系列假设和理论分析，能够给出理论解析解，但局限于简单的几何形状和边界条件。
数值模拟方法通过建立几何模型和边界条件，利用数值计算的方法求解接触过程的力学行为和变形情况，可以适用于复杂的几何形状和边界条件。
ANSYS软件采用的是数值模拟方法，它基于有限元法和多体动力学原理，可以使用接触元素来建立模型，模拟接触过程中的相互作用，得到接触点的应力、应变以及变形信息，从而分析接触的性能和行为。
接下来将介绍ANSYS软件中的接触分析方法和其在工程中的应用。
三、接触分析方法1. 接触元素：ANSYS软件提供了多种接触元素供用户选择，包括面接触元素、体接触元素和线接触元素。
用户可以根据具体的接触问题选择合适的接触元素，建立几何模型来模拟接触行为。
2. 接触定义：在ANSYS软件中，用户可以通过定义接触性质、接触参数和接触约束来描述接触问题。
接触性质包括摩擦系数、接触行为模型等；接触参数包括接触初始状态、接触刚度等；接触约束包括接触面间的约束条件等。
3. 接触分析：通过在ANSYS软件中建立模型，定义接触参数和加载条件，进行接触分析，得到接触点的应力、应变和变形信息。
可以通过分析结果来评估接触性能，发现可能存在的问题，并进行改进和优化。
四、ANSYS软件在工程中的应用1. 机械工程领域：在机械工程中，接触问题广泛存在于各种设备和结构中，如轴承、齿轮、支撑结构等。
接触问题(参考ANSYS的中文帮助文件)当两个分离的表面互相碰触并共切时,就称它们牌接触状态。
在一般的物理意义中,牌接触状态的表面有下列特点:1、不互相渗透；2、能够互相传递法向压力和切向摩擦力；3、通常不传递法向拉力。
接触分类:刚性体－柔性体、柔性体－柔性体实际接触体相互不穿透,因此,程序必须在这两个面间建立一种关系,防止它们在有限元分析中相互穿过。
――罚函数法。
接触刚度――lagrange乘子法,增加一个附加自由度(接触压力),来满足不穿透条件――将罚函数法和lagrange乘子法结合起来,称之为增广lagrange法。
三种接触单元:节点对节点、节点对面、面对面。
接触单元的实常数和单元选项设置:FKN:法向接触刚度。
这个值应该足够大,使接触穿透量小；同时也应该足够小,使问题没有病态矩阵。
FKN值通常在0.1~10之间,对于体积变形问题,用值1.0(默认),对弯曲问题,用值0.1。
FTOLN:最大穿透容差。
穿透超过此值将尝试新的迭代。
这是一个与接触单元下面的实体单元深度(h)相乘的比例系数XX省为0.1。
此值太小,会引起收敛困难。
ICONT:初始接触调整带。
它能用于围绕目标面给出一个“调整带”,调整带内任何接触点都被移到目标面上；如果不给出ICONT值,ANSYS根据模型的大小提供一个较小的默认值(＜0.03＝PINB:指定近区域接触范围(球形区)。
当目标单元进入pinball区时,认为它处于近区域接触,pinball区是围绕接触单元接触检测点的圆(二维)或球(三维)。
可以用实常数PINB调整球形区(此方法用于初始穿透大的问题是必要的) PMIN和PMAX:初始容许穿透容差。
这两个参数指定初始穿透范围,ANSYS 把整个目标面(连同变形体)移到到由PMIN和PMAX指定的穿透范围内,而使其成为闭合接触的初始状态。
初始调整是一个迭代过程,ANSYS最多使用20个迭代步把目标面调整到PMIN和PMAX范围内,如果无法完成,给出警告,可能需要修改几何模型。
《基于ANSYS软件的接触问题分析及在工程中的应用》篇一一、引言随着现代工程技术的快速发展，接触问题在各种工程领域中扮演着越来越重要的角色。
ANSYS软件作为一种强大的工程仿真工具，被广泛应用于解决各种复杂的工程问题，包括接触问题。
本文将详细介绍基于ANSYS软件的接触问题分析，并探讨其在工程中的应用。
二、ANSYS软件接触问题分析1. 接触问题基本理论接触问题是一种高度非线性问题，涉及到两个或多个物体在力、热、电等作用下的相互作用。
在ANSYS软件中，接触问题主要通过定义接触对、设置接触面属性、设定接触压力等参数进行模拟。
2. ANSYS软件中接触问题的分析步骤（1）建立模型：根据实际问题，建立相应的几何模型和有限元模型。
（2）定义接触对：在ANSYS软件中，需要定义主从面以及相应的接触类型（如面-面接触、点-面接触等）。
（3）设置接触面属性：根据实际情况，设置接触面的摩擦系数、粘性等属性。
（4）设定载荷和约束：根据实际情况，设定载荷和约束条件。
（5）求解分析：进行求解分析，得到接触问题的解。
3. 接触问题分析的难点与挑战接触问题分析的难点主要在于高度的非线性和不确定性。
此外，还需要考虑多种因素，如接触面的摩擦、粘性、温度等。
这些因素使得接触问题分析变得复杂且具有挑战性。
三、ANSYS软件在工程中的应用1. 机械工程中的应用在机械工程中，ANSYS软件被广泛应用于解决各种接触问题。
例如，在齿轮传动、轴承、连接件等部件的设计和优化中，ANSYS软件可以模拟出部件之间的接触力和应力分布，为设计和优化提供有力支持。
2. 土木工程中的应用在土木工程中，ANSYS软件可以用于模拟土与结构之间的接触问题。
例如，在桥梁、大坝、建筑等结构的分析和设计中，ANSYS软件可以模拟出结构与土之间的相互作用力，为结构的设计和稳定性分析提供依据。
3. 汽车工程中的应用在汽车工程中，ANSYS软件被广泛应用于模拟汽车零部件之间的接触问题。
ANSYS⾼级接触问题ANSYS⾼级接触问题第⼀章接触问题概述在⼯程中会遇到⼤量的接触问题，如齿轮的啮合、法兰联接、机电轴承接触、卡头与卡座、密封、板成形、冲击等等。
接触是典型的状态⾮线性问题，它是⼀种⾼度⾮线性⾏为。
接触例⼦如图1：分析中常常需要确定两个或多个相互接触物体的位移、接触区域的⼤⼩和接触⾯上的应⼒分布。
接触分析存在两⼤难点：在求解之前，你不知道接触区域、表⾯之间是接触或分开是未知的，表⾯之间突然接触或突然不接触会导致系统刚度的突然变化。
⼤多数接触问题需要计算摩擦。
摩擦是与路径有关的现象，摩擦响应还可能是杂乱的，使问题求解难以收敛。
1.1 接触分类1.1.1 刚－柔⼀个表⾯是完全刚性的—除刚体运动外⽆应变、应⼒和变形，另⼀表⾯为软材料构成是可变形的。
只在⼀个表⾯特别刚硬并且不关⼼刚硬物体的应⼒时有效。
1.1.2 柔－柔两个接触体都可以变形。
1.2 接触单元ANSYS采⽤接触单元来模拟接触问题：跟踪接触位置；保证接触协调性（防⽌接触表⾯相互穿透）；在接触表⾯之间传递接触应⼒（正压⼒和摩擦）。
接触单元就是覆盖在分析模型接触⾯上的⼀层单元。
在ANSYS中可以采⽤三种不同的单元来模拟接触：⾯—⾯接触单元；点—⾯接触单元；点—点接触单元。
不同的单元类型具有完全不同的单元特性和分析过程。
1.2.1 ⾯—⾯接触单元⽤于任意形状的两个表⾯接触，不必事先知道接触的准确位置；两个⾯可以具有不同的⽹格；⽀持⼤的相对滑动；⽀持⼤应变和⼤转动。
例如：⾯⼀⾯接触可以模拟⾦属成型，如轧制过程。
1.2.2 点—⾯接触单元⽤于某⼀点和任意形状的⾯的接触，可使⽤多个点－⾯接触单元模拟棱边和⾯的接触；不必事先知道接触的准确位置；两个⾯可以具有不同的⽹格；⽀持⼤的相对滑动；⽀持⼤应变和⼤转动。
例：点⾯接触可以模拟棱边和⾯之间的接触。
1.2.3 点－点接触单元⽤于模拟单点和另⼀个确定点之间的接触。
建⽴模型时必须事先知道确切的接触位置；多个点－点接触单元可以模拟两个具有多个单元表⾯间的接触；每个表⾯的⽹格必须是相同的；相对滑动必须很⼩；只对⼩的转动响应有效。
接触问题的关键在于接触体间的相互关系（废话，），此关系又可分为在接触前后的法向关系与切向关系。
??法向关系：??在法向，必须实现两点：1）接触力的传递。
2）两接触面间没有穿透。
??ANSYS通过两种算法来实现此法向接触关系：罚函数法和拉格朗日乘子法。
??1．罚函数法??是通过接触刚度在接触力与接触面间的穿透值（接触位移）间建立力与位移的线性关系：?? 接触刚度*接触位移=法向接触力??对面面接触单元17*，接触刚度由实常数FKN来定义。
穿透值在程序中通过分离的接触体上节点间的距离来计算。
接触刚度越大，则穿透就越小，理论上在接触刚度为无穷大时，可以实现完全的接触状态，使穿透值等于零。
但是显而易见，在程序计算中，接触刚度不可能为无穷大（否则病态），穿透也就不可能真实达到零，而只能是个接近于零的有限值。
??以上力与位移的接触关系可以很容易地合并入整个结构的平衡方程组K*X=F中去。
并不改变总刚K的大小。
这种罚函数法有以下几个问题必须解决：??1）接触刚度FKN应该取多大??2）接触刚度FKN取大些可以减少虚假穿透，但是会使刚度矩阵成为病态。
??3）既然与实际情况不符合的虚假穿透既然是不可避免的，那么可以允许有多大为合适??因此，在ANSYS程序里，通常输入FKN实常数不是直接定义接触刚度的数值，而是接触体下单元刚度的一个因子，这使得用户可以方便地定义接触刚度了，一般FKN取0.1到1中间的值。
当然，在需要时，也可以把接触刚度直接定义，FKN输入为负数，则程序将其值理解为直接输入的接触刚度值。
??对于接近病态的刚度阵，不要使用迭代求解器，例如PCG等。
它们会需要更多的迭代次数，并有可能不收敛。
可以使用直接法求解器，例如稀疏求解器等。
这些求解器可以有效求解病态问题。
??穿透的大小影响结果的精度。
用户可以用PLESOL,CONT,PENE来在后处理中查看穿透的数值大小。
接触问题的关键在于接触体间的相互关系（废话，），此关系又可分为在接触前后的法向关系与切向关系。
法向关系：在法向，必须实现两点：1）接触力的传递。
2）两接触面间没有穿透。
ANSYS通过两种算法来实现此法向接触关系：罚函数法和拉格朗日乘子法。
1．罚函数法是通过接触刚度在接触力与接触面间的穿透值（接触位移）间建立力与位移的线性关系： 接触刚度*接触位移=法向接触力对面面接触单元17*，接触刚度由实常数FKN来定义。
穿透值在程序中通过分离的接触体上节点间的距离来计算。
接触刚度越大，则穿透就越小，理论上在接触刚度为无穷大时，可以实现完全的接触状态，使穿透值等于零。
但是显而易见，在程序计算中，接触刚度不可能为无穷大（否则病态），穿透也就不可能真实达到零，而只能是个接近于零的有限值。
以上力与位移的接触关系可以很容易地合并入整个结构的平衡方程组K*X=F中去。
并不改变总刚K的大小。
这种罚函数法有以下几个问题必须解决：1）接触刚度FKN应该取多大2）接触刚度FKN取大些可以减少虚假穿透，但是会使刚度矩阵成为病态。
3）既然与实际情况不符合的虚假穿透既然是不可避免的，那么可以允许有多大为合适因此，在ANSYS程序里，通常输入FKN实常数不是直接定义接触刚度的数值，而是接触体下单元刚度的一个因子，这使得用户可以方便地定义接触刚度了，一般FKN取到1中间的值。
当然，在需要时，也可以把接触刚度直接定义，FKN输入为负数，则程序将其值理解为直接输入的接触刚度值。
对于接近病态的刚度阵，不要使用迭代求解器，例如PCG等。
它们会需要更多的迭代次数，并有可能不收敛。
可以使用直接法求解器，例如稀疏求解器等。
这些求解器可以有效求解病态问题。
穿透的大小影响结果的精度。
用户可以用PLESOL,CONT,PENE来在后处理中查看穿透的数值大小。
如果使用的是罚函数法求解接触问题，用户一般需要试用多个FKN值进行计算，可以先用一个较小的FKN值开始计算，例如。
接触问题（参考ANSYS的中文帮助文件）当两个分离的表面互相碰触并共切时，就称它们牌接触状态。
在一般的物理意义中，牌接触状态的表面有下列特点：1、不互相渗透；2、能够互相传递法向压力和切向摩擦力；3、通常不传递法向拉力。
接触分类：刚性体－柔性体、柔性体－柔性体实际接触体相互不穿透，因此，程序必须在这两个面间建立一种关系，防止它们在有限元分析中相互穿过。
――罚函数法。
接触刚度――lagrange乘子法，增加一个附加自由度（接触压力），来满足不穿透条件――将罚函数法和lagrange乘子法结合起来，称之为增广lagrange法。
三种接触单元：节点对节点、节点对面、面对面。
接触单元的实常数和单元选项设置：FKN：法向接触刚度。
这个值应该足够大，使接触穿透量小；同时也应该足够小，使问题没有病态矩阵。
FKN值通常在0.1~10之间，对于体积变形问题，用值1.0（默认），对弯曲问题，用值0.1。
FTOLN：最大穿透容差。
穿透超过此值将尝试新的迭代。
这是一个与接触单元下面的实体单元深度（h）相乘的比例系数，缺省为0.1。
此值太小，会引起收敛困难。
ICONT：初始接触调整带。
它能用于围绕目标面给出一个“调整带”，调整带内任何接触点都被移到目标面上；如果不给出ICONT值，ANSYS根据模型的大小提供一个较小的默认值（＜0.03＝PINB：指定近区域接触范围（球形区）。
当目标单元进入pinball区时，认为它处于近区域接触，pinball区是围绕接触单元接触检测点的圆（二维）或球（三维）。
可以用实常数PINB调整球形区（此方法用于初始穿透大的问题是必要的）PMIN和PMAX：初始容许穿透容差。
这两个参数指定初始穿透范围，ANSYS 把整个目标面（连同变形体）移到到由PMIN和PMAX指定的穿透范围内，而使其成为闭合接触的初始状态。
初始调整是一个迭代过程，ANSYS最多使用20个迭代步把目标面调整到PMIN和PMAX范围内，如果无法完成，给出警告，可能需要修改几何模型。
a n s y s接触问题牛人的
经验之谈
集团企业公司编码：（LL3698-KKI1269-TM2483-LUI12689-ITT289-
接触问题的关键在于接触体间的相互关系（废话，），此关系又可分为在接触前后的法向关系与切向关系。
法向关系：在法向，必须实现两点：1）接触力的传递。
2）两接触面间没有穿透。
ANSYS通过两种算法来实现此法向接触关系：罚函数法和拉格朗日乘子法。
1．罚函数法是通过接触刚度在接触力与接触面间的穿透值（接触位移）间建立力与位移的线性关系：&n b s p;&n b s p;&n b s p;&n b s p;接触刚度*接触位移=法向接触力对面面接触单元17*，接触刚度由实常数FKN来定义。
穿透值在程序中通过分离的接触体上节点间的距离来计算。
接触刚度越大，则穿透就越小，理论上在接触刚度为无穷大时，可以实现完全的接触状态，使穿透值等于零。
但是显而易见，在程序计算中，接触刚度不可能为无穷大（否则病态），穿透也就不可能真实达到零，而只能是个接近于零的有限值。
以上力与位移的接触关系可以很容易地合并入整个结构的平衡方程组K*X=F中去。
并不改变总刚K的大小。
这种罚函数法有以下几个问题必须解决：1）接触刚度F K N应该取多大2）接触刚度FKN取大些可以减少虚假穿透，但是会使刚度矩阵成为病态。
3）既然与实际情况不符合的虚假穿透既然是不可避免的，那么可以允许有多大为合适因此，在ANSYS程序里，通常输入FKN实常数不是直接定义接触刚度的数值，而是接触体下单元刚度的一个因子，这使得用户可以方便地定义接触刚度了，一般FKN取0.1到1中间的值。
当然，在需要时，也可以把接触刚度直接定义，FKN输入为负数，则程序将其值理解为直接输入的接触刚度值。
对于接近病态的刚度阵，不要使用迭代求解器，例如PCG等。
它们会需要更多的迭代次数，并有可能不收敛。
可以使用直接法求解器，例如稀疏求解器等。
这些求解器可以有效求解病态问题。
穿透的大小影响结果的精度。
用户可以用PLESOL,CONT,PENE来在后处理中查看穿透的数值大小。
如果使用的是罚函数法求解接触问题，用户一般需要试用多个FKN值进行计算，可以先用一个较小的FKN值开始计
算，例如0.1。
因为较小的FKN有助于收敛，然后再逐步增加FKN值进行一系列计算，最后得到一个满意的穿透值。
FKN的收敛性要求和穿透太大产生的计算误差总会是一对矛盾。
解决此矛盾的办法是在接触算法中采用扩展拉格朗日乘子法。
此方法在接触问题的求解控制中可以有更多更灵活的控制。
可以更快的实现一个需要的穿透极限。
2．拉格朗日乘子法与扩展拉格朗日乘子法拉格朗日乘子法与罚函数法不同，不是采用力与位移的关系来求接触力，而是把接触力作为一个独立自由度。
因此这里不需要进行迭代，而是在方程里直接求出接触力（接触压力）来。
K x=F+F c o n t a c t 从而，拉格朗日乘子法不需要定义人为的接触刚度去满足接触面间不可穿透的条件，可以直接实现穿透为零的真实接触条件，这是罚函数法所不可能实现的。
使用拉格朗日乘子法有下列注意事项：1）刚度矩阵中将有零对角元，使有些求解器不克使用。
只能使用直接法求解器，例如波前法或系数求解器。
而PCG之类迭代求解器是不能用于有零主元问题的。
2）由于增加了额外的自由度，刚度阵变大了。
3）一个可能发生的严重问题，就是在接触状态发生变化时，例如从接触到分离，从分离到接触，此时接触力有个突变，产生chattering（接触状态的振动式交替改变）。
如何控制这种chattering，是纯粹拉格朗日法所难以解决的。
因此，为控制chattering，ANSYS采用的是罚函数法与拉格朗日法混合的扩展拉格朗日乘子法。
在扩展拉格朗日法中，可以采用实常数TOLN来控制最大允许穿透值。
还有最大允许拉力FTOL。
这两个参数只对扩展拉格朗日乘子法有效。
在扩展拉格朗日乘子法里，程序按照罚函数法开始，与纯粹拉格朗日法类似，用TOLN来控制最大允许穿透值。
如果迭代中发现穿透大于允许的TOLN值，（对178单元是TOLN，而对面面接触单元171-174则是FTOLN）则将各个接触单元的接触刚度加上接触力乘以拉格朗日乘子的数值。
因此，这种扩展拉格朗日法是不停更新接触刚度的罚函数法，这种更新不断重复，直到计算的穿透值小于允许值为止。
尽管与拉格朗日法相比，扩展拉格朗日法的穿透并不是零，与罚函数法相比，可能迭带次数会更多。
扩展拉格朗日法有下列优点：1）较少病态，个接触单元的接触刚度取值可能更合理。
2）与罚函数法相比较少病态，与单纯的拉格朗日法相比，没有刚度阵零对角元。
因此在选择求解器上没有限制，PCG等迭代求解器都可以应用。
3）用户可以自由控制允许的穿透值TOLN。
（如果输入了TOLN，而使用
罚函数法，则程序忽略它）切向关系：摩擦的处理与法向接触力类似。
由于摩擦是非对称的，使问题变的更为复杂。
ANSYS缺省是做对称求解，即使用对称求解器作近似求解。
但是可以改变几个选项强迫做非对称求解。
非对称求解更精确，但是计算量大许多。
