有限元理论基础
有限元的理论基础
有限元的理论基础有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。
采用不同的权函数和插值函数形式，便构成不同的有限元方法。
1.加权余量法：是指采用使余量的加权函数为零求得微分方程近似解的方法称为加权余量法。
（Weigh ted residual method WRM ）是一种直接从所需求解的微分方程及边界条件出发，寻求边值问题近似解的数学方法。
加权余量法是求解微分方程近似解的一种有效的方法。
设问题的控制微分方程为:在V 域内 在S 边界上式中 :L 、B ——分别为微分方程和边界条件中的微分算子；f 、g ——为与未知函数u 无关的已知函数域值；u ——为问题待求的未知函数 ()0B u g -=(5.1.2)()0L u f -=(5.1.1)混合法对于试函数的选取最方便，但在相同精度条件下，工作量最大。
对内部法和边界法必须使基函数事先满足一定条件，这对复杂结构分析往往有一定困难，但试函数一经建立，其工作量较小。
无论采用何种方法，在建立试函数时均应注意以下几点：(1)试函数应由完备函数集的子集构成。
已被采用过的试函数有幂级数、三角级数、样条函数、贝赛尔函数、切比雪夫和勒让德多项式等等。
(2)试函数应具有直到比消除余量的加权积分表达式中最高阶导数低一阶的导数连续性。
(3)试函数应与问题的解析解或问题的特解相关联。
若计算问题具有对称性，应充分利用它。
显然，任何独立的完全函数集都可以作为权函数。
按照对权函数的不同选择得到不同的加权余量计算方法，主要有：配点法、子域法、最小二乘法、力矩法和伽辽金法。
其中伽辽金法的精度最高。
2、虚功原理——平衡方程和几何方程的等效积分“弱”形式虚功原理包含虚位移原理和虚应力原理，是虚位移原理和虚应力原理的总称。
有限元理论基础及应用
有限元理论基础及应用有限元理论是应用于工程计算领域的一种数值分析方法，它是通过将连续的结构或物体分割成有限数量的离散单元，然后在每个单元上进行近似计算，最终得到整个结构或物体的近似解。
有限元理论广泛应用于结构分析、流体力学、电磁场分析等领域，是工程计算的重要工具。
有限元理论的基础是有限元方法，它将连续的结构或物体以网格的形式划分成一系列有限的单元，通过在每个单元内进行节点位移或其他物理量的近似表示，建立起离散的数学模型。
在有限元方法中，常用的单元形状包括线元、三角形单元、四边形单元等。
每个单元的节点之间通过连接的方式形成整个结构的网格。
有限元理论的基本原理是将连续的物理问题转化为离散的代数问题，通过求解代数方程组得到数值结果。
其基本步骤包括：1.离散化：将连续的结构或物体划分为离散的单元，并在每个单元上建立近似解。
2.建立单元方程：根据结构或物体的本构关系、边界条件等，建立每个单元的方程。
3.组装：根据单元之间的连接方式，将每个单元的方程组装成整个结构或物体的方程。
4.边界条件处理：考虑边界条件对结构或物体的约束作用，修改方程。
5.求解代数方程组：将边界条件处理后的方程组进行求解，得到数值解。
有限元理论的应用非常广泛，主要包括：1.结构分析：有限元方法在结构力学领域的应用非常广泛，可以用于预测结构的应力、变形、疲劳寿命等。
例如，在建筑工程中，可以使用有限元方法对建筑结构进行静力分析，以确保结构的稳定性和安全性。
2.流体力学：有限元方法在流体力学领域的应用包括流体流动、传热、空气动力学等方面。
通过将流体分割成离散的单元，并建立流体的动量方程、能量方程等，可以模拟和预测流体的各种特性。
3.电磁场分析：有限元方法可以用于模拟和分析电磁场的分布、辐射、散射等现象。
在电子器件设计中，有限元方法可以用于预测电磁场的影响和优化设计。
此外，有限元方法还应用于声学、热力学、生物力学等领域。
它的优势包括模拟结果的准确性、适用于复杂几何形状和边界条件、计算速度较快等。
02-01有限元分析基础-理论基础
Kq=f——————（1） 其中：K是整体刚度矩阵；
q是节点位移矩阵； f是载荷矩点位移 解有限元方程Kq=f可得到位移。在根据方
程组的特点来选择合适的计算方法。
通过上述分析了解到，有限元分析的基本 思路是“先离散在组装”，离散为了进行单 元分析，组装为了对整体结构进行分析。
σ=Eε—————（2-4） 将式（2-2）、式（2-3）代入到式（2-4） 后简化得到：
F=(AE/l)Δl—————（2-5) 式（2-5）与弹簧方程F=kx很相似。因此， 受轴向力作用的等截面杆看做一个弹簧，则：
keq=AE/l——————（2-6）
一、有限元分析理论基础
根据上述分析，杆件的截面面积都是在 一个方向上变化的。可以将杆件近似地看做 是由4个弹簧串联起来的模型。
(2)假定一个近似描述单元特性解 为研究典型单元的力学特性，不妨先考虑
横截面积为A、长度为l的杆件在外力F作用下 构件的变形。
杆件的平均应力由下式给出： σ=F/A————（2-2） 杆件的平均正应变ε为
ε=Δl/l————（2-3)
一、有限元分析理论基础
在弹性区域内，应力和应变服从胡克定 律，即：
1.2 定义单元特性 (2)定义单元的力学关系
根据单元的材料、形状、尺寸、节点数目、 位置等参数，找出单元节点力和节点位移的 关系式。 (3)计算等效节点力
物理模型离散化后，假定力是通过节点在 单元间进行传递的，但对于实际连续体，力 是通过单元的公共界面在单元间进行传递。
一、有限元分析理论基础
1.3 组装单元 利用结构中力的平衡条件和边界条件将各
利用以上模型，假定力施加在各节点上。 可根据有图中节点1~节点5的受力情况, 得到各节点上力的静平衡： 节点1：R1-k1（u2-u1）=0 节点2：k1（u2-u1）-k2（u3-u2）=0 节点3：k2（u3-u2）-k3（u4-u3）=0 节点2：k3（u4-u3）-k4（u5-u4）=0 节点2：k4（u5-u4）-P=0
有限元基础及应用
（3）若采用位移作为首先求解的基本变量，则 可以使问题的求解变得更规范一些，下面就基 于 A、B、C 三个点的位移 来进行以上问题的 求解。
方法二：节点位移求解及平衡关系
要求分别针对每个连接节点，基于节点的位移来构建 相应的平衡关系，然后再进行求解。
课程介绍
一、课程内容： 1、有限元法理论基础； 2、应用ANSYS有限元软件对汽车/机械结构进
行分析。 二、学习方法：
理论与实践相结合，即通过应用有限元分析实 际问题来掌握有限元理论。 三、学时数：54学时（36学时理论+18学时实 验） 四、考核方式：平时成绩+上机考试+笔试成绩
第一章 绪论
1.1 有限元法概述 有限元法诞生于20世纪中叶（1943
PA
C
A
Rc
a
b
A'
A
C
A
图（a）所示一平衡的杠杆，对C点写力矩
PB
平衡方程：
(a)
B
PB b PA a
图（b）表示杠杆绕支点C转动时的刚体位
移图：
B b A a
综合可得：
(b)
B
即：
B
PA b B PB a A PAA PBB 0
上式是以功的形式表述的。表明：图a的
B'
平衡力系在图b的位移上作功时，功的总
目前应用较多的通用有限元软件如下表：
软件名称 MSC/Nastran MSC/Dytran MSC/Marc ANSYS ADINA ABAQUS
简介 著名结构分析程序，最初由NASA研制 动力学分析程序 非线性分析软件 通用结构分析软件 非线性分析软件 非线性分析软件
CAE课有限元分析理论基础
类型。
精度要求
03
根据问题对精度的要求，选择足够高阶的有限元以保证求解精
度。
常用有限元的介绍
四面体有限元
适用于解决三维问题，具有较高的计算效率 和适应性。
壳体有限元
适用于解决薄壁结构问题，能够模拟结构的 弯曲和变形。
六面体有限元
适用于解决二维和三维问题，精度较高但计 算效率较低。
梁有限元
适用于解决细长结构问题，能够模拟结构的 轴向拉伸和弯曲。
CAE课有限元分析理论基础
目 录
• 引言 • 有限元分析的基本原理 • 有限元的分类和选择 • 有限元分析的实现过程 • 有限元分析的应用实例 • 结论与展望
01 引言
目的和背景
目的
有限元分析（FEA）是一种数值分析方法，用于解决复杂的工程问题，如结构 分析、热传导、流体动力学等。本课程旨在使学生掌握有限元分析的基本原理 和应用。
弯曲有限元
适用于解决大变形问题，如结 构动力学、流体动力学等。
非线性有限元
适用于解决非线性问题，如塑 性力学、断裂力学等。
耦合有限元
适用于解决多物理场耦合问题 ，如流体-结构耦合、电磁-热
耦合等。
有限元的选择
问题特性
01
根据问题的物理特性、边界条件和求解精度要求选择合适的有
限元类型。
计算资源
02
考虑计算资源的限制，选择计算效率高、内存占用小的有限元
04 有限元分析的实现过程
建立模型
确定分析对象和边界条件
首先需要明确分析的对象和所受的边界条件， 这是建立有限元模型的基础。
几何建模
根据分析对象的特点，利用CAD软件建立几何 模型。
模型简化
有限元法的基础理论
一、里兹法与迦辽金法（摘自电磁场有限元方法 金建铭） 1. 里兹法里兹法是一种变分方法，其中边值问题用变分表达式（也称泛函）表示，泛函的极小值对应于给定边界条件下的控制微分方程。
通过求泛函相对于其变量的极小值可得到近似解。
2. 伽辽金法伽辽金法属于残数加权方法类型，它通过对微分方程的残数求加权的方法得到方程的解。
若u是方程的近似解，将u 代入方程可得到非零的残数： r Luf =- u的最佳近似应能使残数r 在Ω内所有点上有最小值。
残数加权方法要求： 0i i R rd ωΩ=Ω=⎰这里i R 表示残数的加权积分，i ω是所选的加权函数。
在伽辽金法中，加权函数与近似解展开中所用的函数相同。
通常，这样可得到最精确的解。
二、有限元方法里兹法和伽辽金法中，在整个解域内找出能表示或至少近似表示问题真实解的试探函数是非常重要的。
然而对于许多问题，这个步骤是十分困难的，对二维和三维问题尤其如此。
为此，我们可将整个区域划分成小子域，并应用定义在每个子域上的试探函数。
因为子域是小区域，因而在每一子域内函数的变化不大，所以定义在子域上的试探函数通常比较简单。
这正是有限元法的基本思想。
应用里兹法的过程通常称为里兹有限元法或变分有限元法，而应用伽辽金方法的过程通常称为伽辽金有限元方法。
有限元法与经典里兹法和伽辽金法的不同之处是在试探函数的公式上。
在经典里兹法和伽辽金法中，试探函数由定义在全域上的一组基函数组成。
这种组合必须能够（至少近似）表示真实解，也必须满足适当的边界条件。
在有限元法中，试探函数是由定义在组成全域的子域上的一组基函数构成。
因为子域很小，所以定义在子域上的基函数能够十分简单。
三、关于形函数（摘自有限元法在电磁计算中的应用 张榴晨）对于一个待求的微分方程，用一组线性独立的尝试函数i ψ和待定系数i C 来表示方程的近似解，并用加权余数法（迦辽金法）来求解这些待定系数。
求解待定系数的代数方程组为：1[]1,2,,ni j i j i d C q d j n ψψψΩΩ=∇∇Ω=Ω=∑⎰⎰这里j ψ为所选择的加权函数，应用迦辽金法时，所选取的加权函数即为尝试函数。
有限元法的理论基础
有限元法的理论基础有限元法是一种离散化的数值计算方法，对于结构分析而言，它的理论基础是能量原理。
能量原理表明，在外力作用下，弹性体的变形、应力和外力之间的关系受能量原理的支配，能量原理与微分方程和定解条件是等价的。
下面介绍有限元法中经常使用的虚位移原理和最小势能原理。
1.虚位移原理虚位移原理又称虚功原理，可以叙述如下：如果物体在发生虚位移之前所受的力系是平衡的(物体内部满足平衡微分方程，物体边界上满足力学边界条件），那么在发生虚位移时，外力在虚位移上所做的虚功等于虚应变能（物体内部应力在虚应变上所做的虚功）。
反之，如果物体所受的力系在虚位移（及虚应变）上所做的虚功相等，则它们一定是平衡的。
可以看出，虚位移原理等价于平衡微分方程与力学边界条件。
所以虚位移原理表述了力系平衡的必要而充分的条件。
虚位移原理不仅可以应用于弹性性力学问题，还可以应用于非线性弹性以及弹塑性等非线性问题。
2.最小势能原理最小势能原理可以叙述为：弹性体受到外力作用时，在所有满足位移边界条件和变形协调条件的可以位移中，真实位移使系统的总势能取驻值，且为最小值。
根据最小势能原理，要求弹性体在外力作用下的位移，可以满足几何方程和位移边界条件且使物体总势能取最小值的条件去寻求答案。
最小势能原理仅适用于弹性力学问题。
2.2有限元法求解问题的基本步骤弹性力学中的有限元法是一种数值计算方法，对于不同物理性质和数学模型的问题，有限元法的基本步骤是相同的，只是具体方式推导和运算求解不同，有限元求解问题的基本步骤如下。
2.2.1问题的分类求解问题的第一步就是对它进行识别分析，它包含的更深层次的物理问题是什么？比如是静力学还是动力学，是否包含非线性，是否需要迭代求解，要从分析中得等到什么结果等。
对这些问题的回答会加深对问题的认识与理解，直接影响到以后的建模与求解方法的选取等。
2.2.2建模在进行有限元离散化和数值求解之值，我们为分析问题设计计算模型，这一步包括决定哪种特征是所要讨论的重点问题，以便忽略不必要的细节，并决定采用哪种理论或数学公式描述结果的行为。
有限元分析理论基础
16
第17页/共39页
三、物理方程（本构关系）
1、有限元本构关系的矩阵形式为：
s De
对于三维情况有：
1 0
0
0
1
0
0
0
1 0
0
0
De
E
0
(1 )(1 2 )
0
0 1 2 0 2 1 2
0
0 0 0 0
0
2
0
0
0
0
0
1 2 2
2024/7/27
17
第18页/共39页
15
第16页/共39页
二维问题的应变-位移关系可简化为：
u
xx yy xy
x v
y
u x
v y
x
0
y
0
y
uv
u
x
一维问题的应变-位移关系可进一步简化为：
xx
u x
x
u
u
则应变-位移关系可以简记为统一的矩阵形式：
u
2024/7/27
由于有限元采用的多项式位移插值函数全部满 足相容条件，只要求了解这一概念，具体形式不作 要求。
2024/7/27
19
第20页/共39页
虚功原理及虚功方程
PA
C
A
Rc
a
b
A'
A
C
A
图 1-8
PB (a)
B
图1-8a示一平衡的杠杆，对C点
写力矩平衡方程：
PA b
PB
a
图1-8b表示杠杆绕支点C转动时
还要注意，当位移是在某个约束条件下发生时，则在该约束力方向
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
有限元理论基础有限元理论基础2.1 数值模拟技术2.1.1数值模拟技术简介在工程技术领域中许多力学问题和场问题，实质上就是在一定的边界条件下求解一些微分方程。
对于少数简单问题，人们可以通过建立它们的微分方程与边界约束求出该问题的解析解。
但是对于比较复杂的数学方程问题以及不规则的边界条件通过激吻戏法往往难以求解，而需要借助各种数值模拟方法活的相应的工程数值解，这就是所谓的数值模拟技术。
在实际工程领域中，用数值模拟技术可以对复杂的工程结构进行受力和响应分析，这样可以在设计或者加工前预知实体结构工作状态下的大概情况。
目前在工程实际应用中，常用的数值求解方法有：有限单元法、有限差分法、边界元等但从实用性和使用范围来说，有限单元法则是随着计算机技术的发展而被广泛应用的一种行之有效的数值计算方法。
2.2.2 有限元法有限元法是一种基于能量原理的数值计算方法，是解决工程实际问题的一种有效的数值计算工具。
它是里茨法的另一种表示形式，它可应用里茨法分析的所有弹性理论。
限元法是处理连续的结构体离散或有限个单元集合，也就是将连续的求解域离散为一定数量的单元集合体。
且每个单元都具有一定的节点，相邻单元通过节点相互连续，同时使用等效节点力代替作用于单元上的力和选定场函数的节点值作为基本未知量。
并在每一单元中假设一个近似插值函数以表示单元中场函数的分布规律：进而利用力学中的某些变分原理去建立用以求解节点未知量的有限元法方程，从而将一个连续域中的无限自由度问题化为离散域中的有限自由度问题。
求解后，可利用解出的节点值和设定的插值函数确定整个单元集体上的场函数。
有限元求解问题中的单元分析：t tt akF=式中：:t F单元节点作用力。
tK:单元刚度矩阵。
ta：单元节点位移。
通过单元分析确定单元刚度矩阵，建立单元节点作用力和单元为伊关系。
有限元求解问题时建立的结构整体平衡方程：PKU=式中：P—结构整体等效点力载荷K—结构总体刚度矩阵U—结构节点位移阵列单元内力的计算：tDBa=σ式中：D—弹性矩阵P—应变矩阵整个结构的有限元分析就是一句上述方程而进行的具体的有限元求解过程如图结构离散生引入约束输出节计算并输结构总体结构节点求解线性2.2有限元法的基础理论2.2.1 有限元法理论在有限元法中，单元的应变—位移关系可表示为：But=式中：t—应变向量u—位移向量B—应变—位移变换矩阵单元的应力—应变关系表示为：εσD=式中：σ—应力向量D—材料相关系数在线弹性材料条件下，D矩阵是一个常量：在非线性弹性材料中，D矩阵上市应变t的函数。
有限元刚度方程为：PKu=式中：P—结构总体刚度矩阵K—单元刚度矩阵其中单元刚度矩阵K为：⎰=vt DBdV BK式中：V—积分域对于非线性弹性材料而言，D矩阵和单元刚度矩阵K均是应变t和位移u的函数。
在小变形问题中，矩阵B与位移u没有相关关系。
而在大变形问题时，矩阵B和单元刚刚度矩阵K则均是位移u的函数。
在轮胎分析中，轮胎由于充气和垂直载荷等作用，轮胎结构会产生较大的变形，轮胎集合结构的这种变形属于集合非线性问题。
轮胎结构本身又是多种材料构成的复合体，其材料属性既有各向同性又有各项异性，这属于材料非线性问题。
在轮胎的静态解触、自由滚动和动态接触状态下，轮胎与路面之间的接触，又涉及到接触非线性问题提。
2.2.2应力—应变理论有限元是里茨法的另一种表示形式，它可应用里茨法分析的所有弹性理论，而应力—应变理论则是里茨分析法的弹性理论的基础。
因此在有限元分析中，一般使用弹性理论俺就载荷作用下物体中的内力状态和变形规律。
1.应力物体收到外力的作用时发生变形，这种变形改变了物体内各分子间的间距，在物体内形成了一个内立场。
当内力和外力相互平衡时，变形不再继续，物体达到稳定平衡状态。
这种由于物体受外力的作用因其物体的变形，而导致内部各部分之间因相对位置改变而因其的相互作用，这种相互作用成为内力。
所谓应力，就是指分布内力系在物体内某一点处的强弱的相互作用。
为了研究物体内某一点C 处的内力，假设用以经过点C 的截面mm 将物体分开，在这选取包含点C 的一个部分进行研究。
如图2.2所示。
围绕点C 取微笑面积A ∆，A ∆上存在着分布内力系的合力F ∆，如图2.2b 所示。
F ∆的大小和方向与点C 的位置和A ∆的大小密切相关。
A F ∆∆与的比值成为平均应力。
A F P m∆∆= m P 为一矢量，表示在A ∆范围内，单位面积上内力的平均集度，称作平均应力。
随着A ∆的逐渐变小，m P 的大小和方向都将逐渐变化。
当A ∆趋近于零时，m P 的大小和方向都将趋近于一定极限P ，即C 点应力P 为A Fp p A m A ∆∆==→∆→∆00lim lim应力P 是分布内力系F ∆在点C 处的集度，反应了分布内力系F ∆在点C 处的强弱成都。
对于应力P ，威力表征其与无提议的形变或者材料的相关性，通常将应力P 分解成垂直于截面的分量正应力σ和切与截面的分量切应力τ。
2. 应变应变表示物体收到外力的作用时发生变形的强弱程度。
在图2.2a 中，物体中的M 点因变形位移到'M 点，'MM 为无体统变形时M 点的位移。
这里假设物体受到约束，没有刚性位移，M 点的位移全是由变形引起的，假想在M 点附近取平面与坐标平面平行的正六面体（当正六面体的边长趋于无限小时为单元体），设该六面体的棱边边长分别为zy x l l l ,,。
变形后其边长和棱边夹角都发生微小变化将单元体投影到xy 平面，如图2.3b 所示。
变形前单元体平行于x 轴的边长长度为x l ，变形后，点M 和点N 分别位移到'M 点和'N 点。
''N M 的长度为x lx ∆+，且MN N M x -=∆''比值：x m l x MNMN N M ∆=-=''ε 表示单元体边长x l 的平均长度变化，成为平均应变。
当MN 长度趋近于0时，则mε的极限为：x N MN m l x MN MN N M ∆=-=→→00lim ''lim ε2.2.3 大变形的有限元数学描述在轮胎结构分析中，轮胎变形问题属于集合非线性范畴，此时小变形情况下的几何方程不再适合轮胎的结构分析，为此必须重新定义新的平衡方程。
这主要是方程和平衡方程不能消除刚性运动的影响，从而无法度量大变形物体的形态。
在度量物体的变形时，需要选取一个特定的结构为基准在变形问题中，一般有两种参考构型。
设初始时刻t=0,时质点坐标为()3,2,1=i X i ，iX 表示任意时刻t 质点的位置。
当选择初始构型为参考构型时则iX 为自变量。
u x X -= 应变表达式为：⎪⎪⎭⎫ ⎝⎛∂∂∂∂+∂∂+∂∂=⎪⎪⎭⎫ ⎝⎛-∂∂∂∂=j k i k i j j i v j k i k v X u X u X u X u X x X x E 2121δ式中：vE 为格林应变 j kk i v x X X x ∂∂∂∂=δ为克朗内记号。
若选择现时构型，即i x 作为自变量，其应变表达式为：⎪⎪⎭⎫ ⎝⎛∂∂∂∂-∂∂+∂∂=⎪⎪⎭⎫ ⎝⎛∂∂∂∂-=j k i k i j j i j ki k v v x u x u x u x u x x X x e 2121δ 式中：v e 为阿尔曼西应变。
2.2.4 大变形状态下的平衡方程描述在变形情况下，应力是和应变相关的，在外部载荷的作用下，物体发生形变，并在变形后重新达到平衡。
因此，用变形后状态即现时构型构成的Cauchy 应力平衡方程更是合适。
设物体在现时构型中占据的区域为V ，它的边界为A ，其中A 由ut A A ,两部分构成，同事设现时构型的单位体积的体积载荷和表面载荷分别为ii q p ,，如图所示。
这时在区域V 的Cauchy 应力平衡方程为：0=+∂∂j j u p x σ应力边界条件：j j u q n=σ 式中：j n 为现时构型边界i A 的外法线方向余弦。
结合Cauchy 应力张量和Kirchohoff 应力张量可以将大变形问题的平衡方程表述为：01001001dV u p dA u q E S j V j A A v v ⎰⎰⎰+=δδδ 2.2.5 大变形增量问题轮胎的充气和与路面接触分析都与其变形过程有着密切关系，对于这些非线性问题要想获得一个合适解，就必须采用增量法。
所谓增量法，即将时间历程离散或时间序列，即132,1,,,0+=n n t t t t t t 在t t t i i ∆+→的增量阶段，需要选用一个参考构型，一般使用两种方式：一种是完全拉格朗日方程，即选取某一时刻t 的构型为参考构型。
通常在轮胎这类大变形问题中，一般选用拉格朗日方程求解比较合适。
2.2.6接触问题在轮胎的静态接触、自由滚动和动态接触状态下，轮胎与路面之间发生接触行为。
在ANSYS 中，接触属于一种高度的状态非线性行为，它是状态非线性类型中一个特殊而又比较重要的子集部分。
接触问题主要体现在（1）接触区域的大小、相互接触物体的相对位置以及接触的具体状态事先都是都是未知的；（2）接触无提议之间不可相互侵入，接触界面间的法向作用只存在压力，切向可以滑移，这些条件共同构成高度非线性的单边性不等式约束。
在状态非线性分析中，刚度矩阵[K]相关于位移矩阵[x]。
()[]{}[]F x x K =式中：力与位移的关系是非线性的。
其关系如图所示接触一般分为刚—柔接触和柔柔接触，接触方式一般分为三种形式：面面接触、点面接触和点点接触。
轮胎路面接触属于典型的刚—柔面面接触。
在刚柔面面接触中，需要把一个个面设置为“目标”面，而把另一个面设置为“接触面”，通常选择刚性路面为“目标面”，而接触面选择柔性的胎面花纹表面。
这两个面共同组成“接触对”，使用相同的实常数号，有限元程序则是通过这组相同的实常数号来识别该接触对。
接触分析中，接触区域需要同时满足运动和力学两方面约束条件。
运动约束条件是指接触的两个物体在接触边界上法向无相对运动，即无提议之间不存在相互渗透。
力学约束是指接触面法向方向上只能存在接触压力，而在接触面的切向方向上会产生大小相等方向相反的摩擦力。
两个接触物体在接触过程中可能会出现三种情形：粘式、滑移和分离。
粘式指接触的两个无提议在接触点粘在一处，不存在相对运动。
滑移是指接触的两个物体在接触边界上发生切向运动。
分离是指两个接触的物体在变形前存在接触，但在变形后发生分离。
23本章小结：。