有限元求解非线性问题

2）对数应变和真实应力 对数应变/自然应变/真实应变是度量大应变的方法，计算公式如下：
它是非线性应变的度量，因此是关于最终长度的非线性函数。与线性应变相比，对数应变（或真实应变）是可加
的。考虑一个初始长度为1m的杆经过下面3步的变形： 第1步： 从1m 变形至1.2m 第2步：从1.2m 变形至1.5m 第3步：从1.5m变形至2m 在下表中我们比较了工程应变和真实应变。可以清楚地看到，只有真实应变是可加的，因此在非线性分析中应该
大位移和大转角（小应变；线性或非线性材料）
大位移、大转角和大应变（线性或非线性材料）
K.J. Bathe, Finite Elemente Methoden 在线性FEA中，应变，如x方向应变可写为εx = ∂u/∂x，也就是说在表达式εx = ∂u/∂x + ...[(∂u/∂x)z + (∂v/∂x)z + (∂w/∂x)z]中只考虑了一次项的影响。在大位移（非线性）中，表达式的二次项也要考虑。另外，材料的应力-应变关 系也不一定是线性的。 2）材料非线性
材料非线性的特点
非线性材料（小位移）
K.J. Bathe, Finite Elemente Methoden 所有的工程材料本质上都是非线性的，因为无法找到单一的本构关系满足不同的条件比如加载、温度和应变率。 可以对材料特性进行简化，只考虑对分析来说重要的相关因素。线弹性材料（胡克定律）假设是最简单的一种。如果 变形可恢复，则材料为线弹性，如果变形不可恢复，则为塑性。如果温度效应对材料属性影响较大，则应该通过热弹性或热-塑性关系考虑结构和热之间的耦合效应。如果应变率对材料有明显影响，则应使用粘-弹性或粘-塑性理论。 上图是一个材料非线性的示例。 材料非线性的简单分类： 1. 非线性弹性 2. 超弹性 3. 理想弹-塑性 4. 弹性-时间无关塑性 5. 时间相关塑性（蠕变） 6. 应变率相关弹-塑性 7. 温度相关的弹性和塑性 如果考察上图中的应力-应变曲线，则材料非线性可以分为以下几类： 1. 线弹性-理想塑性 2. 线弹性-塑性。应力-应变曲线的塑性段与时间无关，还可细分为两种：

在程序中，对增量方程求解的平衡迭代采用修正 的牛顿迭代法或BFGS法。 1. 修正的牛顿迭代法。它与完全的牛顿法的不同在 于迭代过程中系数矩阵保持不变，因此不需要重新形 成和分解刚度阵，从而大大减少了计算量。但是这样 又带来了收敛速度慢和发散问题，对此程序中加入了 加速收敛和发散处理的措施。这些措施并不明显地增 加求解的时间，但却会对修正的牛顿迭代法的性能有 所改进。 2. BFGS法。又称矩阵修正迭代，是拟牛顿法的一 种。它实际上是完全的牛顿法与修正的牛顿法之间的 一种折中方法。因为它在迭代过程中，并不重新形成
0 t t t k xi N k0 xik， xi N kt xik， xi N kt t x（ i 10-28） k 1 k 1 k 1 n n n
0 k t k t t k 其中： xi ， xi ， xi 为节点k，i方向上在0，t， t+△t时刻的
返回
取位移插值函数为： n
t
写成矩阵形式：
t i
ui N u
k 1
t k k i
；
ui N k uik
k 1
n
（10-26） （10-27）
u [N ] u
t k i
；
ui [ N ]uik
其中：Nk为插值函数，[N]为形函数矩阵； t k ui ，uik 为k点i方向上t时刻的位移和位移增量； n为单元节点数。 取坐标变换为：
v
v s
{R} [ N ]T qv dv [ N ]T qs ds {R0}
{u}

外载荷阵 （10-6） 为节点位移对时间的二 次导数；
为节点位移对时间的一 次导数。
{u}

有限元非线性
本课件介绍《有限元非线性》课程的重要概念和应用领域，帮助学习者深入 了解非线性有限元分析的基本原理和解决方案。
有限元分析基础概念
介绍有限元分析的基本原理，包括离散化方法、单元类型和刚度矩阵的计算。
进一步学习非线性有限元方法
深入讨论非线性有限元方ቤተ መጻሕፍቲ ባይዱ的应用和优缺点，以及适用场景。
常见的非线性问题类型
弹性-塑性耦合模型
讨论弹性和塑性耦合的模型，以及其在结构分析和变形分析中的应用。
本构方程的求解方法
详细介绍求解非线性本构方程的数值方法和迭代策略，包括线性化方法和增量迭代法。
探讨材料非线性、几何非线性和边界条件非线性等常见问题类型，并提供解决方案。
经典弹塑性模型
介绍经典弹塑性模型及其在非线性有限元分析中的应用，包括塑性流动准则和硬化规律。
渐进式塑性模型
探讨渐进式塑性模型的特点及其在复杂材料行为建模中的应用。
黏塑性模型
介绍黏塑性模型及其在某些材料和地质工程分析中的应用，如粘土和岩石。

非线性有限元法综述摘要：本文针对非线性有限元法进行综述，分别从UL列式及TL列式、CR列式、几何精确梁、壳理论三个方面介绍其分析思路和发展动态，旨在为相关学者提供一些思路参考。

关键词：几何非线性；UL列式；TL列式；CR列式；几何精确梁、壳理论1引言几何非线性是由于位置改变引起了结构非线性响应。

进行结构几何非线性分析，实质上就是要得到结构真实的变形与受力情况。

有限元方法是进行结构几何非线性分析的最成熟的方法，也是应用最广泛的分析方法.2非线性有限元法研究思路非线性有限元法主要指UL列式法、TL列式法、CR列式法和几何精确梁、壳理论等，它们有着基本相同的思路，即利用虚功原理建立平衡方程。

方程中充分考虑了非线性因素对结构应变和应力的影响，也就是将线性应变和非线性应变都代入到表达式中，然后确定单元的本构关系并选取合适的形函数，导出单元对应的弹性刚度矩阵和几何刚度矩阵，再选取合适的增量-迭代算法进行求解，由此就完成了结构的整个几何非线性分析求解过程。

非线性有限元法将结构的变形过程划分为三个主要阶段：C0状态、C1状态和C2状态，如图1所示。

图1 单元的变形C0状态是单元的初始状态，C1状态是单元受力变形后上一次处于平衡的状态；C2状态是单元的当前状态，也就是所求的状态。

2.1UL法和TL法研究思路UL法和TL法为几何非线性问题提供了新的分析思路。

这两种方法本质上没有很大区别，但是方程建立的参考状态有所不同。

完全拉格朗日法(TL法)是以结构变形前C0状态为参考建立平衡方程的，考虑结构从C0状态到C2状态之间的变形；而更新的拉格朗日法(UL法)以结构变形后C1状态为参考建立平衡方程的[2]，考虑结构从C1状态到C2状态之间的变形。

两种拉格朗日法的主要形式如下：（1）TL列式（2）UL列式从上面两式可以看出：TL法和UL法的另一个不同是TL法的增量平衡方程中考虑了初位移矩阵的影响，而UL法则忽略了其影响，只考虑了弹性刚度矩阵和初应力矩阵的影响。

第7章 非线性有限元
——材料非线性问题
1. 材料非线性问题的求解方法 2. 塑性应力应变关系 3. 弹塑性矩阵的表达式 4. 弹塑性问题的求解方法 5. 弹塑性问题的实例计算
7.1 材料非线性问题的求解方法
前面各章中，我们所讨论的问题都是线弹性力学问题。在线弹 性力学中，位移与应变的关系（几何方程）是线性的，应变与 应力的关系（本构方程）也是线性的。
(x
xn
)＝0
它的解是
Δx
Y
( xn
)
/(
dY dx
)n；
xn1 xn＋Δxn1
这就是牛顿－拉斐逊方法的迭代公式。
牛顿－拉斐逊方法的迭代过程如图7.1(a)所示，它要求在每 次迭代时计算 Y '(x) dY (x) / dx ，因此计算工作量巨大。
修正的牛顿－拉斐逊方法 迭代公式是
Δx
Y
( xn
dY ＝dF
d d
KT
式中，KT是曲线F＝Kδ的斜率，代表切 线刚度。第二步，从B1点作曲线F＝Kδ 的切线交直线F＝R于A2点，取A2点的横 坐标是δ2。从图中看出
A1B1 /(2 1)＝(KT )1
由于 A1B1＝R F(1)；2 1 (Δ )2
得
A1B1＝(KT )（1 Δ）2＝R F(1)
f ({σ},{ε}) 0
(7.2)
必须注意，由于小变形的关系，应力形式的平衡方程仍然是线 性的，但是以结点位移列阵{δ} 表示的平衡方程则不再是线性的 了。因为应力{σ}和应变{ε} 之间是非线性的，从而应力{σ}与位 移{δ}之间也是非线性的；于是(7.1)式可以写成
K({δ}){δ}{R} 0
(7.3)
1 牛顿-拉斐逊（Newton-Raphson）法

0
t t t k xi N k0 xik， xi N kt xik， xi N kt t x（ i 10-28） k 1 k 1 k 1
n
n
n
0 k t k t t k 其中： xi ， xi ， xi 为节点k，i方向上在0，t， t+△t时刻的 节点坐标值。
（10-25）

T T t T t t t e C e dv dv W e t tv t t t tv t tv t dv
此为改进的拉格朗日（ U·L ）公式。 三、非线性问题有限元基本方程 有了方程（10-19），（10-25）式，就可以按通常的方 法进行有限元离散，从而得到非线性问题的有限元基本方程。
第一节
有限元基本方程
一、线性问题的基本方程 由复杂结构受力平衡问题的虚功方程有:
T T T v v v s s
dv u q dv u q ds u R
T 0 0
mu u dv Du u dv
[M ]
t t
{u} [ D]

t t
{u} [ K ]t t {u} t t {R} （10-8）

解此方程也用隐式时间积分，显式时间积分或振形迭加 法求解。
二、非线性问题的基本方程 对于非线性问题通常不能用一步直接求解方案，必须分成 若干步加载,按各个阶段不同的非线性性质逐步求解，即增量求 解方案。 1.增量形式的平衡方程： 已知设：0，△t，2△t‥‥的位移和应力（各载荷步的） 要求出：t+△t步时的位移和应力。 ①全拉格朗日（T·L）公式 以t=0时刻状态为度量基准，求t+△t时刻的值。 由虚功方程： 其中：

工程力学中的非线性分析方法有哪些？在工程力学领域，非线性问题的研究至关重要。

与线性问题相比，非线性问题更加复杂，需要采用专门的分析方法来准确描述和解决。

下面我们就来探讨一下工程力学中常见的非线性分析方法。

首先要提到的是有限元法。

这是一种非常强大且广泛应用的数值分析方法。

在处理非线性问题时，它能够有效地模拟材料的非线性行为，比如塑性、蠕变等。

通过将复杂的结构离散为有限个单元，并对每个单元进行分析，最终得到整个结构的响应。

对于几何非线性问题，如大变形、大转动等，有限元法能够通过更新坐标和刚度矩阵来准确捕捉结构的变化。

而对于材料非线性，如弹塑性问题，通过定义合适的本构关系，可以精确地模拟材料在不同应力状态下的行为。

再来看看边界元法。

它是另一种有效的数值方法，特别适用于处理无限域或半无限域问题。

在非线性分析中，边界元法可以结合迭代算法来求解非线性边界条件或非线性材料特性。

与有限元法相比，边界元法通常只需要对边界进行离散，从而降低了问题的维数，减少了计算量。

但在处理复杂的非线性问题时，其数学推导和编程实现可能会相对复杂。

还有一种方法是摄动法。

这是一种基于微扰理论的分析方法。

对于弱非线性问题，通过将非线性项视为对线性问题的小扰动，将问题的解表示为一个级数形式。

通过求解这个级数的各项，可以逐步逼近非线性问题的精确解。

摄动法在处理一些简单的非线性问题时非常有效，但对于强非线性问题，其精度可能会受到限制。

接下来是增量法。

在处理非线性问题时，将加载过程或变形过程分成一系列的小增量。

在每个增量步内，将问题近似为线性问题进行求解，然后逐步累加得到最终的结果。

这种方法适用于各种非线性问题，尤其是在考虑加载历史和路径相关性的情况下。

非线性有限差分法也是常用的手段之一。

它直接对控制方程进行离散，通过差分近似来表示导数项。

在处理非线性问题时，可以采用迭代的方式求解离散后的方程组。

这种方法对于简单的几何形状和边界条件的问题较为适用，但对于复杂的结构可能会面临网格划分和精度控制的挑战。

非线性有限元之非线性求解方法平衡回顾✧静态平衡是内力I和外载P力量平衡；✧在非线性问题中，模型的内力I可以是以下量的非线性函数；✧在非线性问题中，模型的外力P也可以是某些量的非线性函数，如位移u和时间t。

非线性求解方法1.已知一个分析，知道结构总载荷和初始刚度，目的是找到最后的位移。

线性分析中，一次计算就能求解出最终位移；非线性问题中不可能，因为结构刚度随着结构变形而改变。

2.求解这类非线性问题需要的是一种增量\迭代技术，获得的解是非线性问题准确的近似。

这些方程通常没有精确解。

3.Abaqus使用迭代求解该方程：使用牛顿拉普森方法求解近似解，使误差最小。

4.Abaqus用法：1）载荷历史被拆解为一系列的分析步；每个分析步拆解为一系列增量步；用户为初始时间增量猜测一个值；Abaqus使用自动增量算法确定其他的增量步。

在每个增量步结束时，Abaqus根据载荷与时间关系计算当前负载大小2）使用牛顿拉普森程序迭代求解每个增量结束时的解；根据收敛容差判断牛顿拉普森程序的收敛；如果迭代不收敛，减少增量步的大小；然后使用小增量步重新进行计算。

5.分析步、增量步、迭代步1）分析步仿真载荷历程含有一个或多个分析步。

2）增量步是分析步的一部分；在静态问题中，总载荷被分成很小的增量步。

以便可以沿着非线性路径求解。

3）迭代步迭代步是增量步中寻找平衡解得一次计算尝试。

5.牛顿拉普森方法Abaqus/Standard 基于牛顿拉普森方法的增量迭代求解技术，该方法是无条件稳定（任何大小的增量步都可以）。

增量步大小影响动态分析精度，每个增量步通常要求多次迭代才能满足收敛要求，每个分析步通常有多个增量步，牛顿拉普森定义了一个残差为0位移曲线。

6.牛顿拉普森方法基础。

平衡是u的非线性方程，牛顿拉普森迭代求解在Cu 处的线性方程，Cu是位移u的修正量。

7.残差定义为了得到线性方程组，重写一下平衡方程，R（u）是u的残差。

这个残差表示的是位移u处不平衡力。

于是方程的解为 u ( KT ( un ,n ))1( R ( un ,n )) ( KT ( un ,n ))1( n1R P( un ))
(9)
(10 )
•在增量方法中通常引入载荷因子λ,用 R R表示载荷， 于是非线性有限
元方程可写成： ( u, ) P( u ) R 0
(1)
用载荷因子λ系列： 0 0 1 2 M 1
(2)
相应于不同的载荷。
若相应于载荷因子 n 的解已经求得，记为 u un ，则 ( un ,n ) P( un ) n R 0
KT n
KT ( un
)
un
(8)
un1 un un
其收敛判据与直接迭代法的收敛判据雷同。
非线性有限元方程组的解法（增量法）
•求解非线性方程组的另一类方法是增量方法。使用增量方法的一个优点是 可以得到整个载荷变化过程的一些中间的数值结果。当问题的性质与加载的 历史有关时，例如弹塑性问题，则必须采用增量方法。
u1 ( K1 )1 R
据此容易写出直接迭代法的迭代公式：
Kn K( un )
un1 ( K n )|1 R
(2)
按照这种迭代公式可以得到一个解数列 { un } ，当这个数列收敛时停止计
算，其数列收敛值就是方程（1）的解。
非线性有限元方程组的解法（直接迭代法）
关于数列收敛的判据，可以采用各种各样的范数定义和收敛判据。若设第 n
( un ) K( un )un R 0
(7)
该值可作为对偏离平衡的一种度量（称为失衡力），收敛判据可相应地取为：
( un ) R
(8)
（失衡力收敛判据）
非线性有限元方程组的解法（牛顿法）
把非线性有限元方程记为： ( u ) P( u ) R 0 (1)

非线性结构有限元分析的步骤与流程
• 设定边界条件和载荷，如固定约束、压力 或力矩等。
非线性结构有限元分析的步骤与流程
01 步骤三：求解
02
选择合适的求解器，如Newton-Raphson迭代法或 直接积分法。
03 进行迭代计算，求解非线性结构的内力和变形。
非线性结构有限元分析的步骤与流程
01
步骤四：后处理
非线性有限元分析的基本概念
总结词
非线性有限元分析是一种数值分析方法，通过将复杂的结构或系统离散化为有限个小的单元，并建立 每个单元的数学模型，来模拟和分析结构的非线性行为。
详细描述
非线性有限元分析是一种基于离散化的数值分析方法，通过将复杂的结构或系统划分为有限个小的单 元（或称为有限元），并建立每个单元的数学模型，来模拟和分析结构的非线性行为。这种方法能够 考虑各种复杂的边界条件和材料特性，提供更精确的数值结果。
非线性有限元分析的常用方法
总结词
非线性有限元分析的常用方法包括迭代法、增量法、 降维法等。这些方法可以根据不同的非线性问题选择 使用，以达到更好的分析效果。
详细描述
在非线性有限元分析中，常用的方法包括迭代法、增量 法、降维法等。迭代法是通过不断迭代更新有限元的位 移和应力，逐步逼近真实解的方法；增量法是将总载荷 分成若干个小的增量，对每个增量进行迭代计算，最终 得到结构的总响应；降维法则是通过引入一些简化的假 设或模型，将高维的非线性问题降维处理，以简化计算 和提高计算效率。这些方法各有优缺点，应根据具体的 非线性问题选择使用。
03
02
弹性后效
材料在卸载后发生的变形延迟现象。
材料强化
材料在受力过程中发生的强度增加 现象。
04

Ki-1
（三）混合法 如对同一非线性方程组混合使用增量
法和迭代法，则称为混合法或逐步迭代法。 一般在总体上采用Euler增量法，而在
同一级荷载增量内，采用迭代法。
Ki-1
刚度的取值可根据给定的应力-应变曲 线导出。若每级计算都采用上一级增量计算 终了时的刚度值，则称为始点刚度法。
Ki-1
始点刚度法类似于解微分方程初值问题 的欧拉(Euler)折线法，计算方法简单但计算 精度较低，容易“漂移”。
若采用中点刚度法则可以提高精度。该 法类似于解常微分方程初值问题的龙格-库塔 （Runge-Kutta）法，包括中点切线刚度法 和中点平均刚度法。
(1) 直接迭代法 对非线性方程组
设其初始的近似解为 ，由此确定近似的
矩阵
可得出改进的近似解
重复这一过程，以第i次近似解求出第i＋1 次近似解的迭代公式为直接迭代法
对非线性方程组
直到 变得充分小，即近似解收敛时，终止迭代。
在迭代过程中,得到的近似解一般不会满足 作为对平衡偏离的一种度量，称为失衡力。
q-Newton—Raphson迭代法的计算过程
（2）初应力法 如果在弹性材料内确实存在初应力 ，则材料的应力应变关系为
由上式及虚功原理可导出单元的结点力为
集合单元得出以下的有限元方程 式中， 为由初应力 引起的等效结点荷载
初应力法就是将初应力看作是变化的， 以此来反映应力和应变之间的非线性关系。 通过不断地调整初应力，使线弹性解逼近非 线性解。
接触非线性 由于接触体的变形和接触边界的摩擦作用，
使得部分边界条件随加载过程而变化，且不 可恢复。这种由边界条件的可变性和不可逆 性产生的非线性问题，称为接触非线性。
材科非线性有限元法 材料非线性是由本构关系的非线性引

【问题描述】如图I所示的模型，纵向尺寸均为100mm，水平尺寸均为30mm，圆角半径均为10mm，模型厚度为4mm。

图I 本例中所使用的模型【要求】在ANSYS Workbench软件平台上，通过改变材料属性，分别对该模型进行线性材料静力分析以及非线性材料的静力分析，并加以对比。

1.分析系统选择（1）运行ANSYS Workbench，进入工作界面，首先设置模型单位。

在菜单栏中找到Units下拉菜单，依次选择Units>Metric（kg,mm,s,℃，mA，N,mV）命令。

（2）在左侧工具箱【Toolbox】下方“分析系统”【Analysis Systems】中双击“静力结构分析”【Static Structural】系统，此时在右侧的“项目流程”【Project Schematic】中会出现该分析系统共7个单元格。

相关界面如图1所示。

图1 Workbench中设置静力分析系统2.输入材料属性（1）在右侧窗口的分析系统A中双击工程材料【Engineering Data】单元格，进入工程数据窗口。

（2）我们首先进行的是线性材料问题，选用系统默认的结构钢作为材料即可。

（3）可以看见，系统本身默认结构钢【Structural Steel】已在备选材料窗口中，在此不必再另行选择，直接单击【Project】选项卡回到项目流程界面即可。

3.导入几何模型（1）双击分析系统A中的“几何”【Geometry】单元格。

（2）找到菜单栏中的文件【File】选项，依次选择【File】>【Import External Geometry File】，在弹出的对话框中找到模型文件“non-linear.igs”并打开。

（3）单击工具栏中的【Generate】选项，即选项，确认生成导入的模型。

导入完成后的模型如图2所示。

（4）至此，模型导入步骤完成。

图2 导入的模型3.网格划分（1）双击Workbench界面中系统A的第四个单元格，模型【Model】单元格，进入【Static Structural】的静力分析模块。

材料非线性有限元分析材料非线性有限元分析是一种重要的计算力学方法，用于研究在载荷作用下，材料会发生非线性行为的情况。

这种分析方法已经被广泛应用于工程领域，例如建筑结构、航空航天以及汽车工业等。

本文将详细介绍材料非线性有限元分析的原理、方法和应用。

首先，我们来介绍一下材料非线性。

在工程领域，材料的非线性行为主要包括弹塑性、损伤、断裂、破坏等。

这些非线性行为往往在高载荷作用下会显著增加结构的应力和应变，从而导致结构的失效。

因此，准确地预测和分析这些非线性行为对于工程设计和结构优化具有重要意义。

材料非线性有限元分析是一种基于有限元方法的计算机模拟技术，用于模拟和分析复杂结构在非线性载荷下的力学行为。

它通过将结构离散为许多小的有限元单元，并以数学模型描述每个单元的材料行为，从而建立了结构的有限元模型。

然后，结构的力学行为可以通过求解相应的离散形式的力学方程得到。

在材料非线性有限元分析中，有两个关键问题需要解决。

首先是材料本构模型的建立。

材料本构模型是描述材料应力和应变关系的数学模型，常用的包括弹性模型、塑性模型、损伤模型等。

选择合适的材料本构模型对准确预测和分析结构的非线性行为至关重要。

其次是数值方法的选择。

对于材料非线性问题，通常需要使用迭代算法，如牛顿-拉夫森法，来求解非线性方程。

此外，还需要选择适当的数值积分方法，以解决离散形式的力学方程。

材料非线性有限元分析在许多领域都有广泛的应用。

在结构工程领域，它可以用于分析钢筋混凝土结构、大跨度桥梁以及高层建筑等的受力性能。

在航空航天领域，材料非线性有限元分析可用于研究飞机机翼、航天器的结构强度和振动特性。

在汽车工业中，它可以用于分析车辆的碰撞、耐久性和振动特性。

总结起来，材料非线性有限元分析是一种重要的计算力学方法，能够准确地模拟和分析结构在非线性载荷下的力学行为。

它在工程领域有着广泛的应用，能够为工程设计和结构优化提供科学依据。

未来随着计算机硬件和数值方法的不断发展，材料非线性有限元分析将在更多领域得到应用，并为解决工程实际问题提供更准确和高效的方法。

增量有限元法
增量有限元法是一种用于求解非线性问题的数值方法，它基于有限元法的思想，通过将求解域划分为有限个单元，并在每个单元上建立近似的函数表达式，来逼近问题的真实解。

增量有限元法的基本思路是将整个时间区间划分为一系列的时间增量，在每个时间增量内，将非线性问题线性化，并通过求解线性方程组来获得该时间增量内的近似解。

然后，将每个时间增量内的近似解进行叠加，得到整个时间区间内的近似解。

增量有限元法具有许多优点，例如可以有效地处理非线性问题、适用于大变形问题、可以提高计算效率等。

它已经被广泛应用于工程领域，如结构力学、流体力学、热力学等。

在实际应用中，增量有限元法需要考虑许多因素，如时间增量的大小、线性化的方法、数值稳定性等。

同时，还需要选择合适的有限元网格和单元类型，以保证计算精度和效率。

总之，增量有限元法是一种非常有效的数值方法，它可以帮助工程师和科学家解决许多复杂的非线性问题，为工程设计和科学研究提供了重要的工具。

• 在加工、密封、撞击等问题中，接触和摩擦的作 用不可忽视，接触边界属于高度非线性边界。
• 平时遇到的一些接触问题，如齿轮传动、冲压成 型、轧制成型、橡胶减振器、紧配合装配等，当 一个结构与另一个结构或外部边界相接触时通常 要考虑非线性边界条件。 • 实际的非线性可能同时出现上述两种或三种非线 性问题。
由于从理论上还丌能提供能普遍接受的据有时非线性材料特性可用数学模型迚行模拟尽管这些模型总有他们的局限性
有限元求解非线性问题
元计算科技发展有限公司
1）材料非线性问题
• 材料的应力和应变是非线性的，但应力与应变却 很微小，此时应变与位移呈线性关系，这类问题 属于材料的非线性问题。由于从理论上还不能提 供能普遍接受的本构关系，所以，一般材料的应 力与应变之间的非线性关系要基于试验数据，有 时非线性材料特性可用数学模型进行模拟，尽管 这些模型总有他们的局限性。在工程实际中较为 重要的材料非线性问题有：非线性弹性（包括分 段线弹性）、弹塑性、粘塑性及蠕变等。
2）几何非线性问题，几何非线性问题是由于 位移之间存在非线性关系引起的
• 当物体的位移较大时，应变与位移的关系是非线 性关系。研究这类问题一般都是假定材料的应力 和应变呈线性关系。它包括大位移大应变及大位 移小应变问题。如结构的弹性屈曲问题属于大位 移小应变问题，橡胶部件形成过程为大应变问题 。
3）非线性边界问题
Thanks！