非线性有限元
线性和非线性有限元
目
CONTENCT
录
• 线性有限元方法 • 非线性有限元方法 • 线性与非线性有限元的比较 • 线性与非线性有限元的实例分析 • 未来研究方向与展望
01
线性有限元方法
定义与原理
定义
线性有限元方法是一种数值分析方法，用于求解偏微分方程的近 似解。它将复杂的求解区域离散化为有限个小的、简单的子区域 ，即有限元，然后对每个有限元进行求解，最终得到原偏微分方 程的近似解。
THANK YOU
感谢聆听
在实际应用中，应根据问题的特性和需求选择合适 的有限元方法。对于复杂的问题，可能需要结合多 种有限元方法进行求解。
05
未来研究方向与展望
线性有限元方法的改进与优化
80%
高效求解算法
研究更快速、稳定的线性有限元 求解算法，提高计算效率。
100%
自适应网格生成
发展更智能、自动的网格生成技 术，以适应复杂几何形状和边界 条件。
线性有限元
由于线性有限元基于线性方程组进行求解，因此计算复杂度 相对较低，适用于求解一些较简单的问题，如弹性力学问题 。
非线性有限元
非线性有限元需要求解非线性方程组，计算复杂度较高，但 能够处理更复杂的问题，如塑性力学、流体力学等领域的问 题。
精度比较
线性有限元
对于一些简单的问题，线性有限元可以给出较为精确的结果。然而，对于一些 复杂的问题，线性有限元可能无法准确描述非线性行为。
80%
多物理场耦合
研究线性有限元在多物理场耦合 问题中的应用，如流体-结构、电 磁-热等。
非线性有限元方法的改进与优化
高阶非线性有限元
发展高阶非线性有限元方法， 以更精确地描述复杂非线性行 为。
非线性有限元介绍1
非线性有限元介绍1.为什么使用FEA解决有限元问题（1）理解设计的意图。
有限元分析（FEA）是研究不同力学设计的有力工具。
（2）降低产品成本和开发周期。
1) FEA通过以下方式降低产品成本和开发周期；2) 在模具制造前识别成型问题；3) 使模具制造返工成本最低；4) 提前识别设计中缺陷减小样机的成本；5) 使用最少的材料；（3）获得结果的唯一办法。
1) FEA可以用来预测产品在极端工况下的性能，这些在实验中无法复现；2) 能在设计阶段提前考虑这些工况。
（4）很多工况在设计阶段无法预料。
2.收敛定义（1）在有限元中收敛有多重意义：1) 网格收敛；2) 时间积分精度；3) 非线性程序收敛；4) 求解精度；（2）网格收敛1) 增加模型单元数量会使仿真解趋于解析解。
网格收敛对线性和非线性问题都适用；在Abaqus中使用H网格自适应技术；2)进一步加密网格时，结果变化很小或不变时，认为网格达到收敛。
3)网格收敛规则的例外：网格奇异解；材料损伤累计在模型特定区域的局部问题；在Abaqus中使用H自适应网格技术；Abaqus提供特殊技术来减小网格依赖性，解决材料软化局部影响。
4)Abaqus提供评估网格收敛工具。
在打印输出文件（.dat）和结果文件（.fil）输出的节点（SJP）处应变跳变。
5)Abaqus后处理云图设置。
应力云图；不连续云图。
6)自适应网格误差设计。
见Abaqus/Standard 自适应网格课程介绍。
（3）瞬态问题时间积分精度。
对于具有物理时间尺度的瞬态问题，Abaqus提供用户定义参数，以控制对相关方程的积分精度：半增量容差。
1)评估当前增量中途点的最大不平衡力；2)一个增量中允许的最高温度变化；3)增量中根据开始和结束速率条件计算出的蠕变应变增量的最大容差。
（4）非线性求解收敛：见后续。
（5）求解精度。
获得精确解需要满足以下条件：准确解需要工程经验来建立合适有限元模型：材料、载荷、边界和求解程序。
非线性结构有限元分析概论
一、线性问题的基本方程
由复杂结构受力平衡问题的虚功方程有:
v T dv vuT qvdv suT qsds u0T R0
vmu
T
••
u dv
v
Du
T
•
u
dv
（10-1）
上式左端为内力的虚功，右端为外力的功。
由于： u N u Bu C
式中 u 为单元体内的位移； u为节点位移； N 形函数阵；
t t t
T
S t t t
dvt
W t t
（10-18）
返回
其中：
W tt o
tv
u
T
q tt tv
中推荐采用BFGS法。
程序对几何非线性的考虑可采用完全的拉格朗
日公式或改进的拉格朗日公式。在非线性动态分析
中采用隐式时间积分（Newmarli法和Wilson- 法） 或显式时间积分（中心差分法）的方法。隐式时间
积分通常用来分析结构的振动问题，显式时间积分
主要用来分析波传布现象。
返回
第一节 有限元基本方程
解此方程也用隐式时间积分，显式时间积分或振形迭加
法求解。
返回
二、非线性问题的基本方程
对于非线性问题通常不能用一步直接求解方案，必须分成
若干步加载,按各个阶段不同的非线性性质逐步求解，即增量求
解方案。
1.增量形式的平衡方程：
已知设：0，△t，2△t‥‥的位移和应力（各载荷步的）
要求出：t+△t步时的位移和应力。
ov oe T o
o e dv
ov
o
T
t o
SdvtW t o来自ovoe Tt o
S
dv
非线性有限元分析1
非线性问题的类型和求解特点1 非线性问题的类型1. 1 线性分析的含义在有限元分析中的线性假设包含下列含义：即结点位移为无限小量，材料为线弹性，加载时边界条件的性质保持不变。
于是，静力平衡方程可以表示为：[]{}{}R U K = (2.1)其中，[]K 为刚度矩阵，{}R 为荷载矢量。
由于[]K 和{}R 的元素为常数，故位移响应{}U 是荷载矢量{}R 的线性函数。
也就是说，如果{}R 变为{}R α，则{}U 变为{}U α，其中，α为常数。
这就是所谓的线性有限元分析。
如果上述假设中的任何一条不能得到满足，那么就属于非线性有限元分析。
1. 2 非线性分析的必要性结构力学问题，从本质上讲都是非线性的，线性假设只是实际工程问题的一种简化。
当然，任何实际工程问题的求解都避免不了适当地简化，简化是否合理主要应根据求解效果和实际经验来判断。
对于目前工程实际中的很多问题，如地震作用下结构的弹塑性动力响应，高层建筑抗风，大跨度网壳结构动力稳定性，索膜结构找形荷载与裁减分析，大型桥梁风致振动等问题的研究，仅仅假设为线性问题是很不够的，常常需要进一步考虑为非线性问题。
因此，对各种工程结构的非线性分析就是必不可少且日趋重要了。
对于结构力学的非线性问题来说，有限单元法是最为有效的数值分析方法。
1. 3 非线性问题的类型通常，把非线性问题分为两大类，即分为几何非线性和材料非线性。
但从建立基本方程和程序设计的方便出发，又可分为三种类型：1．材料非线性：非线性效应仅由应力应变关系的非线性引起，位移分量仍假设为无限小量，故仍可采用工程应力和工程应变来描述，即仅材料为非线性。
非线性的应力应变关系是结构非线性的常见原因，许多因素都可以影响材料的应力应变性质，包括加载历史（如在弹塑性响应状况下），环境状况（如温度），加载的时间总量（如在蠕变响应状况下）等。
2．几何非线性：如果结构经受大变形，则变化了的几何形状可能会引起结构的非线性响应，这又可以分为两种情形：第一种情形，大位移小应变。
《有限元非线性》课件
本课件介绍《有限元非线性》课程的重要概念和应用领域，帮助学习者深入 了解非线性有限元分析的基本原理和解决方案。
有限元分析基础概念
介绍有限元分析的基本原理，包括离散化方法、单元类型和刚度矩阵的计算。
进一步学习非线性有限元方法
深入讨论非线性有限元方ቤተ መጻሕፍቲ ባይዱ的应用和优缺点，以及适用场景。
常见的非线性问题类型
弹性-塑性耦合模型
讨论弹性和塑性耦合的模型，以及其在结构分析和变形分析中的应用。
本构方程的求解方法
详细介绍求解非线性本构方程的数值方法和迭代策略，包括线性化方法和增量迭代法。
探讨材料非线性、几何非线性和边界条件非线性等常见问题类型，并提供解决方案。
经典弹塑性模型
介绍经典弹塑性模型及其在非线性有限元分析中的应用，包括塑性流动准则和硬化规律。
渐进式塑性模型
探讨渐进式塑性模型的特点及其在复杂材料行为建模中的应用。
黏塑性模型
介绍黏塑性模型及其在某些材料和地质工程分析中的应用，如粘土和岩石。
非线性有限元法综述
非线性有限元法综述摘要：本文针对非线性有限元法进行综述，分别从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法则忽略了其影响，只考虑了弹性刚度矩阵和初应力矩阵的影响。
非线性有限元方法
非线性有限元方法非线性有限元方法是大量应用于工程领域的计算方法，它主要用于求解复杂结构的力学问题，例如材料的变形、破坏和变形控制等。
与线性有限元方法不同，非线性有限元方法考虑因为载荷和边界条件的非线性导致问题的非线性本质，以及材料的非线性行为。
在这篇文章中，我们将讨论非线性有限元方法，包括其应用、工作原理以及其在工程领域中的重要性等内容。
首先，我们来研究一下非线性有限元方法的应用。
非线性有限元方法在许多方面都有应用。
其中最重要的领域是结构力学，包括建筑、航空航天、汽车等领域。
由于这些结构需要承受复杂的载荷，因此非线性有限元方法可以很好地模拟这些结构的行为，预测它们的性能和寿命。
此外，非线性有限元方法还可以应用于材料力学研究中，例如破碎、断裂和塑性变形等方面。
其次，我们来了解一下非线性有限元方法的工作原理。
与线性有限元方法类似，非线性有限元方法通过将结构分成小块进行离散，然后在每个小块中进行力学分析，最后将分析结果合并为整个结构的行为。
但是，与线性有限元方法不同的是，非线性有限元方法考虑到材料的非线性行为，采用迭代的方法计算结构的响应。
通常，在每一次迭代中，我们都将结构的当前状态作为一个初始猜测，然后求解出该状态下的切应力和位移场。
然后我们将这个位移场的结果代入底部，从而更新结构的状态。
如果解决方案收敛，则完成计算，否则就将新的状态再次代入求解。
这种方法的本质是将非线性问题转化为一系列线性问题的求解，通过迭代求解来逼近非线性问题的解。
最后，我们来讨论一下非线性有限元方法在工程领域中的重要性。
非线性有限元方法已成为现代工程设计和分析的不可或缺的工具。
它允许工程师们模拟和预测各种工程机构的行为，以及设计和优化各种结构。
例如，它可以帮助我们了解在不同载荷下建筑和桥梁行为的变化，预测材料的破坏和失效，以及优化汽车和飞机的结构以提高其性能。
总之，非线性有限元方法是一种复杂但十分有用的计算方法，它可以模拟各种结构的行为并预测其性能和寿命。
如何利用非线性有限元法进行力学分析
如何利用非线性有限元法进行力学分析非线性有限元法是一种用于数值分析问题的计算方法，其主要应用于力学分析领域。
这种方法在于其对于复杂结构的建模能力和高精度数值计算能力而备受推崇。
在本文中，将介绍如何对力学问题进行分析，以及如何应用非线性有限元法对力学分析进行模拟。
1. 引言力学分析整体上分为两种类型：静力学分析和动力学分析。
静力学分析研究对于物体的力和静止条件进行研究，其中力一般会造成物体的运动。
而动力学分析则研究运动物体的变化，特别是再一定条件下物体的振动问题等。
因为力学分析问题具有很高的复杂性，很多时候需要使用非线性有限元法来得到更准确的结果。
下面我们将详细介绍使用非线性有限元法进行力学分析的方法和流程。
2. 有限元法简介有限元法是一种现代数值计算方法，它将大工程结构分割为小的有限元。
在每个有限元内，结构的物理性质可以被认为是常量。
(具体内容可以自己百度)3. 如何利用非线性有限元法进行力学分析使用非线性有限元法进行力学分析的核心是将宏观问题转变为微观问题来进行模拟计算。
其中需要注意下面几点：3.1 确定力学分析的类型根据要进行分析的结构本身的性质和应用场景，可能涉及到静力学分析或者动力学分析。
其中静力学分析的计算主要涉及到结构在平衡状态下的情况，而动力学分析主要涉及到结构在某种条件下的运动和振动情况。
因此，在进行力学分析之前需要确定其类型，以便进行后续的计算。
3.2 建立结构模型根据具体情况，需要对结构进行建模。
建模可以通过一定的工具软件实现，或者手工建立结构模型。
模型的建立需要考虑到其复杂性和具体的应用场景。
构建好结构模型之后，需要对其进行精细化剖分得到单元网格，并进行编号。
3.3 确定边界条件在进行力学分析时，还需要考虑结构的边界条件。
边界条件可以通过指定某些点的坐标或者某些角度的变化来确定。
因此，在进行计算时需要根据具体情况设定边界条件，以便进行后续的计算。
3.4 进行数值模拟计算运用有限元法的基本原理，将每个单元的机械性质进行计算，根据力学分析的情况，可以得到结构节点的位移、应变和应力等参数。
非线性有限元之非线性求解方法
非线性有限元之非线性求解方法平衡回顾✧静态平衡是内力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处不平衡力。
03非线性分析要点
第三部分非线性分析第一章非线性有限元概述1.1非线性行为1、 非线性结构的基本特征是结构刚度随载荷的改变而变化。
如果绘制一个非线 性结构的载荷一位移曲线，则 力与位移的关系是非线性函数。
2、 引起结构非线性的原因：a 几何非线性：大应变，大位移，大旋转 （例如钓鱼竿的变形）b 材料非线性：塑性，超弹性，粘弹性，蠕变c 状态改变非线性：接触，单元死活3、 非线性行为一一分析方法特点A 不能使用叠加原理！B 结构响应与路径有关，也就是说加载的顺序可能是重要的。
C 结构响应与施加的载荷可能不成比例。
1.2非线性分析的应用1、 一些典型的非线性分析的应用包括： 非线性屈曲失稳分析金属成形研究碰撞与冲击分析制造过程分析（装配、部件接触等）材料非线性分析 （塑性材料、聚合物）2、 橡胶底密封：一个包含几何非线性（大应变与大变形），材料非线性（橡胶）， 及状态非线性（接触）的例子。
2.1非线性方程组的解法1、求解一个结构的平衡问题通常等于求解结构的总位能的驻值 问题。
结构总位能n ： 口 "3弋门心 2、 增量法：就是将荷载分成一系列的荷载增量，即 ANSYS 中的荷载步或荷载子 步。
A 要点：在每一个荷载增量求解完成后，继续进行下一个荷载增量之前， 刚度矩阵以反映结构刚度的变化。
B 增量法的优点：可以追踪结构变形历程，这对于材料或几何非线性（特别是 极限值屈曲分析）十分有用。
C 增量法的缺点：随着荷载步增量的增加而产生积累误差，导致荷载-位移曲 线飘移。
D 对飘移进行平衡修正，可以大大提高增量法的精度。
应用最广的就是在每一 级载荷增量上用Newton-Raphsor 或其变形的迭代法。
3、 迭代法：割线刚度法：收敛性差，因此很少应用切线刚度法Newto n-Ra phsor 迭代法：切向刚度法中 2.2 Newto n-Ra phsor 迭代法 1、 优点：对于一致的切向刚度矩阵有 二次收敛速度。
非线性有限元解法
(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
第13讲 非线性有限元问题的分类与一般解法-11_35620112
10.1 引言
c）应力-应变关系为{σ}=[D]{ε} ，式中[D]为常数矩阵。 ⇒ 基于线弹性假设。 d）边界约束条件在加载过程中保持不变。 如果在加载中位移边界条件发生改变，（如某自由的自 由度在一定载荷水平下成为被约束自由度），则系统成为非 线性。这种情况在接触分析中出现。 * 如果不满足上述①、②、③假设，结构系统的力学行为将 出现非线性。
{ΔPi } ( i = 1, 2,", m ) （增量可不等）
{P} = ∑ {ΔPi }
i =1 m
P
P
ΔPi
0
δ
δLeabharlann ②逐步施加载荷增量，逐步求解。每一步计算中，将 刚度矩阵[K(δ)]处理为常数（线性化），在不同载荷 步中，刚度矩阵具有不同值。即 由
{ΔPi } 线性化 JJJJJJJ J M {Δδ i } → {Δε i } → {Δσ i }
大位移（大转 动），大应变
接触非线性
与其他组合出现
汽车工程系
结构分析与CAE研究室
10.2 非线性问题的分类
用以下例子说明各类非线性.
Δ
P/2
σ ,ε
P/2 L 求 P∼Δ 之关系 A
汽车工程系
结构分析与CAE研究室
10.2 非线性问题的分类
① 线弹性（小位移） 应力： σ = P A 位移： 应变： ε
有限元法应用
(第十三讲)
清华大学汽车工程系 结构分析与CAE研究室
第10章 非线性有限元问题的分类与 一般解法
10.1 引言 10.2 非线性问题的分类 10.3 非线性问题的一般解法
10.1 引言
在线性有限元分析中作了以下假设： ①小位移(小变形)； ②线弹性； ③在加载过程中边界条件不变化。 由此，得静力学有限元方程：[ K ]{δ } = {P} 其中[K]为常数矩阵. 该方程为线性方程 即,
非线性有限元
（三）混合法 如对同一非线性方程组混合使用增量
法和迭代法，则称为混合法或逐步迭代法。 一般在总体上采用Euler增量法，而在
同一级荷载增量内，采用迭代法。
Ki-1
刚度的取值可根据给定的应力-应变曲 线导出。若每级计算都采用上一级增量计算 终了时的刚度值，则称为始点刚度法。
Ki-1
始点刚度法类似于解微分方程初值问题 的欧拉(Euler)折线法，计算方法简单但计算 精度较低，容易“漂移”。
若采用中点刚度法则可以提高精度。该 法类似于解常微分方程初值问题的龙格-库塔 （Runge-Kutta）法，包括中点切线刚度法 和中点平均刚度法。
(1) 直接迭代法 对非线性方程组
设其初始的近似解为 ，由此确定近似的
矩阵
可得出改进的近似解
重复这一过程，以第i次近似解求出第i＋1 次近似解的迭代公式为直接迭代法
对非线性方程组
直到 变得充分小，即近似解收敛时，终止迭代。
在迭代过程中,得到的近似解一般不会满足 作为对平衡偏离的一种度量，称为失衡力。
q-Newton—Raphson迭代法的计算过程
（2）初应力法 如果在弹性材料内确实存在初应力 ，则材料的应力应变关系为
由上式及虚功原理可导出单元的结点力为
集合单元得出以下的有限元方程 式中， 为由初应力 引起的等效结点荷载
初应力法就是将初应力看作是变化的， 以此来反映应力和应变之间的非线性关系。 通过不断地调整初应力，使线弹性解逼近非 线性解。
接触非线性 由于接触体的变形和接触边界的摩擦作用，
使得部分边界条件随加载过程而变化，且不 可恢复。这种由边界条件的可变性和不可逆 性产生的非线性问题，称为接触非线性。
材科非线性有限元法 材料非线性是由本构关系的非线性引
非线性有限元1_非线性方程求解及收敛控制
2014-12-3
18
2.1 非线性方程组的解法－－增量法
增量法：就是将荷载分成一系列的荷载增量，即ANSYS中的 荷载步或荷载子步。 要点：在每一个荷载增量求解完成后，继续进行下一个荷载 增量之前，调整刚度矩阵以反映结构刚度的变化。
公式（ 4） ui 1 ui ui 1
2014-12-3 6
1.1 非线性行为——材料非线性
非线性的应力－应变关系是产生结构非线性的一个 普遍原因。
应力
应力
应变
应变
钢
橡胶
2014-12-3 7
1.1 非线性行为——状态改变非线性
许多非线性问题是与状态相关的。例如一段 绳索可以是松驰的或拉紧的。一个装配件的两部分可 能接触或脱离接触。
在这个接触例题中 ，接触面积未知， 它取决与施加载荷 的大小。
2014-12-3
8
1.1 非线性行为——分析方法特点
不能使用叠加原理! 结构响应与路径有关，也就是说加载的顺序可能是重要的。 结构响应与施加的载荷可能不成比例。
2014-12-3
9
第一章 非线性有限元概述
1.1 非线性行为 1.2 非线性分析的应用
2014-12-3 25
2.2 Newton-Raphson迭代法－－力平衡
Newton-Raphson 法需要一个收敛的度量以决定何 时结束迭代。 给定外部载荷（Fa），内部载荷（ Fnr ,由单元应力 产生并作用于节点），在一个体中，外部载荷必须与内 力相平衡。
Fa - Fnr = 0
收是平衡的度量。
载荷
F
收敛半径 如果 ustart 在收敛半径内将收 敛，否则将发散。
ustart ?
非线性有限元分析
非线性有限元分析1 概述在科学技术领域内，对于许多力学问题和物理问题，人们已经得到了它们所应遵循的基本方程（常微分方程或偏微分方程）和相应的定解条件（边界条件）。
但能够用解析方法求出精确解的只是少数方程性质比较简单，并且几何形状相当规则的问题。
对于大多数工程实际问题，由于方程的某些特征的非线性性质，或由于求解区域的几何形状比较复杂，则不能得到解析的答案。
这类问题的解决通常有两种途径。
一是引入简化假设，将方程和几何边界简化为能够处理的情况，从而得到问题在简化状态下的解答。
但是这种方法只是在有限的情况下是可行的，因为过多的简化可能导致误差很大甚至是错误的解答。
因此人们多年来一直在致力于寻找和发展另一种求解途径和方法——数值解法。
特别是五十多年来，随着电子计算机的飞速发展和广泛应用，数值分析方法已成为求解科学技术问题的主要工具。
已经发展的数值分析方法可以分为两大类。
一类以有限差分法为代表，主要特点是直接求解基本方程和相应定解条件的近似解。
其具体解法是将求解区域划分为网格，然后在网格的结点上用差分方程来近似微分方程，当采用较多结点时，近似解的精度可以得到改善。
但是当用于求解几何形状复杂的问题时，有限差分法的精度将降低，甚至发生困难。
另一类数值分析方法是首先建立和原问题基本方程及相应定解条件相等效的积分提法，然后再建立近似解法并求解。
如果原问题的方程具有某些特定的性质，则它的等效积分提法可以归结为某个泛函的变分，相应的近似解法实际上就是求解泛函的驻值问题。
诸如里兹法，配点法，最小二乘法，伽辽金法，力矩法等都属于这一类方法。
但此类方法也只能局限于几何形状规则的问题，原因在于它们都是在整个求解区域上假设近似函数，因此，对于几何形状复杂的问题，不可能建立合乎要求的近似函数。
1960年，发表了有限单元法的第一篇文献“The Finite Element Method in Plane Stress Analysis”，这同时也标志着有限单元法（FEM）的问世。
第6章 非线性有限元法(几何非线性)分析
FkiFkj ij dxidxi 2eijdxidxi
由于大变形问题有
２、限A元lm方an程sh主i应要变采用张量
T.L列式法或U.L列式 Alm法an建sh立i应，变因张此量应采在用初Eular运动 描述始方状法态，下即定按义当应前变状张态下的构 形定量义，应即变采张用量G。reen应
变ds张2 量d。s2 dxidxi dxidxi
dxidxi dxi Fki1Fkj1dx j
ij Fki1Fkj1 dxidxi 2Eij dxidxi
eij
1 2
FkiFkj ij
式中，eij称为Green应变张量或 Green-Lagrangian应变张量。
Eij
第六章 非线性有限元法（几何非线性）
1、变几形何非体线性的的有运限动元方描程一述 般采用T.L或U.L列式法建立！
变形体上的质点的运动状态 可以随不同的坐标选取以下几 种描述方法：
1、全拉格朗日列式法（T.L列式 法—Total Lagrangian Formulation)：
选取t0=0时刻未变形物体的构 形A0作为参照构形进行分析。
uk xj
ij
ij
式中：
ij
1
ui
2 xj
u j xi
为小变形应变张量；
ij
1 2
uk xi
uk xj
为非线性二次项
2、Green变形张量也可写为：
eij
1 2
Cij
ij
式中，Cij是Cauchy变形张量
Cij FkiFkj
由于Cauchy变形张量是正定对称 阵，因此该张量有三个实特征值； 这些特征值的平方根记为材料的 主轴拉伸。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
初始屈服条件
此条件规定材料开始塑性变形时的应力状态。 对于初始各向同性材料，在一般应力状态下开始 进入塑性变形的条件是：
F F ( , ) 0
O O i, j 0
对于金属材料，通常采用的屈服条件有： （1) V.Mises 条件 （2）Tresca 条件 返回
流动法则
流动法则用来规定材料进入塑性应变后的 塑性应变增量在各个方向上的分量和应力 分量以及应力增量之间的关系。
16.4 有限元求解方程及解法
1. 2. 16.3.4 方程解法 静力分析 动力分析 16.4.4 平衡路径的追踪和载荷步长的选择 16.4.5 依赖于变形的载荷
16.5 大变形条件下的本构关系
大变形问题分为： 1.大位移、大转动和小应变问题—变形很大， 但应变很小，甚至还保持在弹性应变范围 内。 2.大位移、大转动和大应变问题—应变很大， 从材料角度考虑，分为弹性问题和塑性问 题。
16.1 引言
这种由于大位移和大转动引起的非线性问称 为几何非线性问题。 在几何非线性问题的有限单元法中，通常采用增量分析方 法。 增量分析方法一般采用两种表达格式 1. 完全的Lagrange格式：静力学和运动学变量总是参考初 始位形，即整个分析过程中参考位形保持不变。 2. 更新的Lagrange格式：静力学和运动学变量参考 于每一载荷或时间步长开始时的位形，即在分析过 程中参考位形不断在更新。
16.3 几何非线性问题的表达格式
16.3.4 平衡方程的线性优化
1. 物理方程的线性化。 2. 求解格式的进一步线性化。
16.4 有限元求解方程及解法
1. 2. 16.4.1 有限元求解方程 静力分析问题 动力分析问题 16.4.2 用于几何分线性的单元及担忧矩阵 和向量举例 1. 实体元 2. 板壳元
15.8 小结
当前状态确定，而必须的路径和历史确定。 为了适用这种一般情况，所以在有限元分析 中通常采用的是增量型本构关系。
有限单元法
第16章 几何非线性问题
16.1 引言
在以前各章所讨论的问题中都是基于小变 形的假设，即假定问题所发生的位移远小 于物体自身的几何尺度，同时材料的应变 远小于1。 实际上，我们会遇到河道不符合小变形假 设的问题，例如板和壳等薄壁结构在一定 载荷作用下，尽管应变很小，甚至未超出 弹性极限，但是位移较大，材料线元素会 有较大的位移和转动。
图 应力的度量
16.3 几何非线性问题的表达格式
在涉及几何非线性问题的有限元方法中， 通常采用增量分析的方法，这不仅是因为 问题可能涉及依赖于变形历史的材料的非 弹性，而且因为即使问题不涉及材料非弹 性，但为了得到加载过程中应力和变形的 演变历史，已经保证求解的精度和稳定， 通常也需要采用增量方法求解。
有限单元法
第15章 材料非线性问题
15.1 引言
① ② ③
我们在本章之前所讨论的内容均属于线性 问题。下面我们将进入非线性问题的研究。 首先对比线性问题我们来了解一下什么是 非线性问题： 线弹性力学的基本特点是： 平衡方程是不依赖变形状态的线性方程 几何方程的应变和位移的关系是线性的 物理方程的应力和应变关系是线性的
16.3 几何非线性问题的表达格式
16.3.1 虚位移原理 16.3.2 完全Lagrange格式 这种格式中所有变量以时间0的位形为参考 位形。参考于初始位形的并与平衡方程等 效的虚位移原理。 16.3.3 更新的Lagrange格式 这种格式中所有变量以时间t的位形为参考 位形。在求解过程中参考位形是不断改变 的。
I. 直接迭代法 II. N-R方法 III. 增量法 1.欧拉方法 2.N-R方法
15.2 非线性方程组的解法
加速收敛的方法
15.3 材料弹塑性本构关系
15.3.1 材料弹塑性行为的描述
{
1.单调加载 2.反向加载 3.循环加载
15.3 材料弹塑性本构关系
1. 2. 3. 4. 15.3.2 塑性力学的基本法则 初始屈服条件 流动法则 硬化法则 加载、卸载准则
15.4 弹塑性增量有限元分析
15.4.1 弹塑性问题的增量方程 15.4.2 增量的有限元格式
15.5 弹塑性增量分析数值方法中的几 个问题
15.5.1 非线性增量方程组的求解方案 15.5.2 载荷增量步长的自动选择 15.5.3 弹塑性状态的决定和本构关系的积 分 15.5.4 单元刚度矩阵的数值积分 15.5.5 线性方程组的求解
16.6 结构稳定性和屈曲问题
分析的目的是求解结构从稳定平衡过渡到 不稳定平衡的临界载荷和失稳后的屈曲形 态。 1. 线性稳定分析 2. 非线性稳定分析
16.7 小结
本章讨论了建立于非线性连续介质力学基 础上的非线性有限元分析的基本理论和方 法。 关于本构关系，在大变形总的前提下，首 先应区分是大位移、大转动、小应变情况 还是大应变情况；其次应区分是弹性变形 还是非弹性变形。
15.3 材料弹塑性本构关系
15.3.3 弹塑性增量的应力应变关系 1.建立弹塑性应力应变关系需遵循的原则 2.各向同性硬化材料的应力应变关系 3.用于不同问题的具体表达形式 4.其它硬化材料的应力应变关系
15.3 材料弹塑性本构关系
15.3.4 弹塑性全量的应力应变关系 1.材料的 曲线 2.塑性应变表达式 3.弹塑性全量应力应变表达式
16.2 大变形条件下的应变和应力的度量
16.2.1 应变的度量 固定在笛卡尔坐标系内的物体，在外力作用下连续改变位 形。
图 笛卡尔坐标系内物体的运动和变形
16.2 大变形条件下的应变和应力的度量
16.2.2 应力的度量 在大变形问题中，是从变形后的物体中截取出微元体建立 平衡方程和与之相等效的虚功原理，所以应从变形后的物 体内截取单元体定义应力张量－－欧拉应力张量
终于结束啦！ 谢谢大家。
材料非线性问题可分两类
一类是：不依赖于时间的弹塑性问题，其 特点是，当荷载作用以后，材料变形立即 发生并且不再随时间而变化。 另一类是依赖于时间的的粘（弹、塑）性 问题，其特点是，当荷载作用以后，材料 不仅立即发生相应的弹（塑）性变形，而 且变形随时间而继续变化。 返回
15.1 引言
④ 力的边界上的外力和位移边界上的位移是 独立或线性依赖于变形状态的。 非线性问题的特点：不符合任何一个 上述线性力学特点的方程或边界条件的， 则此问题就是非线性的。 非线性问题可分为3类，即 ① 材料非线性问题 ② 几何非线性问题
15.1 引言
③ 边界非线性问题
15.2 非线性方程组的以判别从一塑性状态出发是继续 塑性加载还是弹性卸载，这是计算过程中 判定是否继续塑性变形以及决定是采用弹 塑性本构关系还是弹性本构关系所必须的。
返回
硬化法则
硬化法则是用来规定材料进入塑性变形后 的后继屈服函数（又称加载函数或加载曲 面）在应力空间中变化的规则。 对于硬化材料，通常采用的硬化法则有： （1）各向同性硬化法则 （2）运动硬化法则 （3）混合硬化法则 返回
15.6 弹塑性全量有限元分析
15.6.1 弹塑性全量分析的有限元方程 15.6.2 非线性有限元方程的求解方法
15.7 热弹塑性-蠕变有限元分析
15.7.1 材料的蠕变行为 15.7.2 稳态蠕变分析 15.7.3 热弹塑性-蠕变增量分析
15.8 小结
材料非线性问题有限元分析的基本问题有 两个方面，即材料本构关系的建立和非线 性方程组的解法。 材料本构关系通常分为两类，即全量型本 构关系和增量型本构关系。 前者适用于应力和应变之间存在不依赖于 变形历史和路径的一一对应关系的情况。 由于塑性变形和蠕变变形等是不可恢复的 非弹性变形，应力状态通常不能由变形的