非线性有限元
线性和非线性有限元

目
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 进行数值模拟计算运用有限元法的基本原理,将每个单元的机械性质进行计算,根据力学分析的情况,可以得到结构节点的位移、应变和应力等参数。
- 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 小结
材料非线性问题有限元分析的基本问题有 两个方面,即材料本构关系的建立和非线 性方程组的解法。 材料本构关系通常分为两类,即全量型本 构关系和增量型本构关系。 前者适用于应力和应变之间存在不依赖于 变形历史和路径的一一对应关系的情况。 由于塑性变形和蠕变变形等是不可恢复的 非弹性变形,应力状态通常不能由变形的