弹塑性力学与有限元-材料非线性问题和几何非线性问题
弹塑性有限元法
当变形体同时存在大的弹性和塑性变形时，必 须采用弹塑性力学进行分析，相应的有弹塑性 有限元法，其较一般弹性有限元复杂得多。
1、塑性区中应力与应变之间为非线性关系，非线性问 题求解 — 增量法；
2、应力与应变关系不是一一对应的，加载与卸载关系 不同，必须判断是加载还是卸载状态；
3、多种材料硬化模型产生不同的有限元计算公式；
K u Q 非线性方程组
方程组
求解
与ij 有关
与ij 有关
u tt u t uu
和
三、弹塑性有限元处理的技术问题
1、加载增量步长的选定
计算精度与收敛性
加载的增量步长
tt P t P rmin P
增量步终止载荷
初始设定载荷增量
初始载荷 载荷约束因子
2、变形区弹塑性状态的判定
弹塑性变形过程中，变形体内部可能同时存在弹 性区、过渡区、塑性加载区和塑性卸载区等四种不同 状态的区域和单元，计算时必须分别进行处理。
x xy y xy z xy 2
xy
x yz y yz z yz xy yz 2
yz
x y
zx zx
z xy
zx zx
xy zx 2
zx
二、弹塑性有限元方程
由于 非线性的应力应变关系，只能按照增量法求解。
在小变形条件下，对t到t+Δt时刻的增量步进行 分析。设变形体为各向同性硬化材料、且服从Mises 屈服条件和Prandtl – Reuss方程的本构关系，并设t 时刻的变形条件为：单位体积的体积力为tpi；作用 在边界表面ST上的单位面积力为tTi；任一质点的位
移为tui，应变为tij，应力为tij。现以t时刻的变形为
材料非线性
( ) P( ) f 0 0
其中： 表示载荷变化的量。 dP d d f 0 KT f0 0 d d d d 1 K T ( ) f 0 d 切线矩阵
1 1 m1 m KT ( m ) f0m KT ( m )f m
一、材料弹塑性行为的描述
弹塑性材料进入塑性的特征是当载荷卸去后 存在不可恢复的永久变形，因而在涉及卸载的情 况下，应力和应变之间不再存在一一对应的关系， 这是区别于非线性弹性的基本属性。
11
单调加载 对于大多数材料存在屈服应力，应力低于屈服 应力时，材料为弹性，而当应力超出屈服应力时， 材料进入 弹塑性状态。 当应力达到屈服应力后，应力不再增加，而材料 变形可以继续增加—理想弹塑性材料。
第六章 材料非线性问题的有限元法
1
第一节
引言
线弹性力学基本方程的特点： 几何方程的位移和应变的关系是线性的； 物理方程的应力和应变的关系是线性的； 建立于变形前的平衡方程也是线性的。 几何非线性问题 结构的变形使体系的受力状态发生显著变化， 以致于不能用变形前的平衡方程分析，且位移和应 变的关系不是线性的。
K ( ) f 0
增量法 载荷分为若干步： f 0 , f1 , f 2 , f 3 位移分成若干步： 0 , 1 , 2 , 3 每两步之间增长量为增量。 增量解法的一般做法是： 假设第m步的载荷 f m 和位移 m ； 让载荷增加 f m1 ( f m f ) ，再求解 m1( m )。 如果每一步的增量 f 足够小，解的收敛性 可以得到保证
9
NmR-N方法求解非线性方程组时，收敛速度 较慢，特别是对于结构分析时载荷趋近极限载荷或突 然变软的情况下，收敛速度会很慢。为了加速收敛， 可以采用一些方法，比较常用和有效的是Aitken法。 该方法每隔一次迭代进行一次加速。
土木工程中的非线性问题概述
土木工程中的非线性问题概述张俊（空军工程大学机场建筑工程系，陕西西安710038）摘要：从四个方面对土木工程中的非线性问题进行了概述，分别是钢筋混凝土方面、地球科学方面、基础学科方面和其他相关方面。
钢筋混凝土中非线性包括材料非线性、几何非线性、双重非线性和结构非线性、钢筋混凝土破坏过程分析、结构内力分析、节点处理以及长期荷载分析等方面；地球科学中包括地震等自然灾害孕育过程、灾变临界条件及成灾规律等方面，尤其是地震科学中；基础学科中包括分形力学等方面。
关键词：土木工程；非线性；钢筋混凝土；地球科学中图分类号：U416文献标识码：A1 引言非线性科学是当代自然科学前沿理论，不仅自身发展迅猛，而且对其它学科也产生了极为深远的影响.国家科委已将《非线性科学》列为攀登计划项目。
在北大、南大、复旦等大学中已成立了非线性科学研究中心，非线性科学作为研究生课程己在部分院校讲授，中国科技大学还在本科生中专门开设了“非线性科学班”。
材料力学和结构力学是建立在小变形假设基础上，不考虑物体位置和形态的变化，用变形前的形态建立平衡条件，进行线性近似化可实际上在某些情沉下是不恰当的，不符合精度要求。
线性只是局部的，非线性才是自然界的最为木质的特征。
非线性科学是研究自然界和人类社会中非线性现象共性的基础学科，本文仅对其在土木工程及相关领域中的应用问题进行论述。
主要从四个方面对土木工程中的非线性问题进行了概述，分别是钢筋混凝土方面、地球科学方面、基础学科方面和其他相关方面。
2钢筋混凝土方面2.1 四种非线性包括材料非线性、几何非线性、双重非线性和结构非线性。
1 材料非线性。
材料的非线性是指材料的应为-应变不成线性比例。
材料非线性问题有非线性弹性问题和非线性弹塑性问题之分。
这两种主要区别是在加载历史上，后者是材料超过屈服极限后呈现出的非线性，常见于各种结构的弹塑性分析。
在简单加载过程中的非线性阶段两者并无木质区别，但卸载过程，前者是可逆过程;后者是不可逆的，将出现残余变形。
第四章 土木工程中的几何非线性问题
通过初时构型上的微元体定义的应力称为Kirchhoff应力，用 表示； 通过现时构型的微元体定义的应力称为现时（Updated）Kirchhoff 应力， 用 表示。
8/5/2024
非线性有限元
20
Euler应力张量:τij
在大变形问题中，是用从变形后的物体内截取的微元体来建立平衡方程及与之相 等效的虚功原理的。因此首先在变形后的物体内截取出的微元体上定义应力张量， 称为Euler应力张量； 此应力张量有明确的含义，即代表真实的应力张量。是现 时位形和变形相关的真实应力。
Case-1
同乘以时间增量
增量形式 …
Case-2
可以证明，这两个率都与转动无关
Jaumann 应力率
旋转率
现时Green应变的线性部分
可以证明，这两个率都与转动无关
8/5/2024
非线性有限元
27
三种本构关系间的关系
对于实际的大变形问题，上述三种本构关系并不等价。可以证明，弹性 材料是一种特殊的次弹性材料，超弹性材料是一种特殊的弹性材料。
19
应力是借助于微元体来定义的，但在大变形分析中，必须注意 微元体所在的构型。
与应变类似，连续介质力学理论具有严格的应力定义和多 种不同的应力概念。这里也只介绍后面将要用到的几种。
Euler应力：
从当前构型中取出微元体，在其上定义的应力称为Euler应力，用 表示。Euler应力代表物体的真实应力。然而，当前构型是待求的未知构型， 因而，有必要通过已知构型上的微元体再对应力进行描述。
Kirchhoff应力（增量）和Green应变（增量）。 优点：参考构型不发生变化，本构关系与虚功方程描述形式简单。
非线性有限元分析报告
非线性有限元分析1 概述在科学技术领域内，对于许多力学问题和物理问题，人们已经得到了它们所应遵循的基本方程（常微分方程或偏微分方程）和相应的定解条件（边界条件）。
但能够用解析方法求出精确解的只是少数方程性质比较简单，并且几何形状相当规则的问题。
对于大多数工程实际问题，由于方程的某些特征的非线性性质，或由于求解区域的几何形状比较复杂，则不能得到解析的答案。
这类问题的解决通常有两种途径。
一是引入简化假设，将方程和几何边界简化为能够处理的情况，从而得到问题在简化状态下的解答。
但是这种方法只是在有限的情况下是可行的，因为过多的简化可能导致误差很大甚至是错误的解答。
因此人们多年来一直在致力于寻找和发展另一种求解途径和方法——数值解法。
特别是五十多年来，随着电子计算机的飞速发展和广泛应用，数值分析方法已成为求解科学技术问题的主要工具。
已经发展的数值分析方法可以分为两大类。
一类以有限差分法为代表，主要特点是直接求解基本方程和相应定解条件的近似解。
其具体解法是将求解区域划分为网格，然后在网格的结点上用差分方程来近似微分方程，当采用较多结点时，近似解的精度可以得到改善。
但是当用于求解几何形状复杂的问题时，有限差分法的精度将降低，甚至发生困难。
另一类数值分析方法是首先建立和原问题基本方程及相应定解条件相等效的积分提法，然后再建立近似解法并求解。
如果原问题的方程具有某些特定的性质，则它的等效积分提法可以归结为某个泛函的变分，相应的近似解法实际上就是求解泛函的驻值问题。
诸如里兹法，配点法，最小二乘法，伽辽金法，力矩法等都属于这一类方法。
但此类方法也只能局限于几何形状规则的问题，原因在于它们都是在整个求解区域上假设近似函数，因此，对于几何形状复杂的问题，不可能建立合乎要求的近似函数。
1960年，R.W.CLOUGH发表了有限单元法的第一篇文献“The Finite Element Method in Plane Stress Analysis”，这同时也标志着有限单元法（FEM）的问世。
弧长法——弹塑性力学及有限元
Pm1
m 1
c ( )[( ) 2um ]
2 T 2 1
tg 1(K1 ) Tm
m
1 a um m
1
2 1
2 2
1 (K1 ) R Tm
1 1 (K1 ) Pm Tm
m m a m
am
a
5 增量弧长法
4)由 R 和 Pm 求
i 2 m
i 1
i m
i i m m
5 增量弧长法
1 1 ( ) 21
i i i 2 m i i i 21i u m m 2
i i i 2 m
i
2
i m
i
2
i
2
i
2
i i i i 2 u m 2m m (m ) 0
i i
2
(1 1 1 )( ) 2 (1
非线性代数方程组的数值解法
5 增量弧长法
用迭代法或增量法进行极限分析时，在极值点附近往往可能 不收敛。这时可用增量弧长法来解决。
5 增量弧长法
• 弧长法是一种用于得到不稳定(KT 0)或负刚度矩阵 (KT < 0)问题的数值稳定解的方法。
ri 弧长半径
F
ri
收敛子步
ri ri
ri 平衡路径 u
5 增量弧长法
i 2 i a(m ) 2b m c 0
式中系数为
T a 1 (1i )( 1i )
i T i i b m (1i )[( 2 ) um ]
c (2 )[(2 ) 2u ]
i T i i m
上述式子是从简单情况推出的，如果除 外 均理解为矩阵，即为一般情况的弧长法方程。
有限元方法中材料非线性计算综述
6
(non-associated flow)。关联流动中使用了屈服函数 作为流动势 。关 联流动用来描述由位错诱发的塑性流动， ABAQUS® 中除铸铁外的一般金属与 Cam-Clay 土力学模型采用了关联流动的格式。非关联流动在处理摩擦型塑性流 动方面比关联流动更好，Mohr-Coulomb 和 Drucker-Prager 等模型使用了非关联 流动。关联流动集成的刚度矩阵 K ep 是对称矩阵，在材料不出现软化现象的时候
, H 0
(3)
5
其中 H H1 , H 2 , 为后继屈服条件中的内变量，表征了材料的强化、粘性等各种 复杂性质，其表达形式也可以非常复杂。塑性力学中常用的 Tressca 屈服准则和 Mises 屈服准则是 的两种特殊形式。Mises 屈服准则的应用较多一些，其屈 服与金属拉伸试验结果吻合得更好，而且其函数形式比较光滑。 Mises 屈服准则为：
2
Newton-Raphson 方法要求在给定 u 的时候计算的切线刚度 K ep ， K ep f int u u 。
K ep 与材料的状态有关， K ep 的计算将在后文中提及，现在假定在给定 u 的情况下 K ep 已经算出。
(a) Newton-Raphson 法 图1
(b) Quasi-Newton 法
K ep 还是正定的，容易求解。而非关联流动集成的刚度矩阵是不对称的，容易导
致求解失败。 除了以上与率不相关的塑性流动外， 粘塑性计算中需要定义率相关的流动法 则，粘性流动率 通常是与应力相关的。常见的粘塑性流动法则有： Bingham 模型：
k , Mises 0 Mises Mises 0 0,
3
非线性理论与方法
工程中的非线性理论与方法Ref：1.冯康等编，数值计算方法，国防版。
2.何君毅，工程结构非线性问题的数值解法，国防工业版。
3.王德人编，非线性方程组解法与最优化方法，高教版4.李岳生编，数值逼近，人民教育出版社。
绪论一．非线性问题的广泛性工程中的非线性问题是普遍存在的。
严格地讲，工程中几乎绝大多数复杂问题都具有非线性本质或呈现出非线性现象，仅是在一定的条件之下，我们可将其理想化或简化为线性问题。
因此，曾有学者认为：在物质世界中，无论是宇观、宏观和微观，都是由一定层次结构和功能的非线性系统构成的，也即自然界和现实生活中几乎所有系统都是非线性的。
事实上，正是由于非线性的存在和作用，才孕育出大自然的五彩缤纷、万千气象和人类社会的风云变幻，人类思维的错综差异。
1.数学中的非线性问题：1）.代数插值2）.曲线曲面拟合3）.非线性回归4）.高次代数方程和超越方程5）.非线性方程组6）.非线性常微分方程（组）7）.非线性偏微分方程（初值、边值）8）.非线性规划（无约束、约束）。
2.机械与结构工程中的非线性问题：1）.柔性可变结构的计算（柔索计算）2）.材料非线性问题（弹塑性力学，塑性力学，蠕变力学）3）.复合材料力学4）.几何非线性问题（大变形问题，屈曲问题）5）.边界（接触）非线性6）.非线性动力学（定则振动，随机振动）柔性多体系统动力学7）.非线性系统控制问题8）.传热学中的非线性问题9）.流体力学和空气动力学中非线性问题等等二．非线性问题（系统）的特点尽管工程中的非线性问题涉及到许多学科，内容不尽相同。
但它们都具有如下非线性问题的共同特点：1.系统最终的控制方程均为非线性方程（代数、常微分、偏微分）2.线性迭加原理在整体上不成立，最多只在只局部近似成立例如：基于线性迭加原理的力法方程，杜哈美积分（卷积），振型迭加法等等，在整体上均不成立。
对于非线性问题应用线性问题中的这些求解方法将导致不真实甚至不合理的结果。
什么叫非线性
我们这么想象:如果发生了特征值屈曲,那么发生屈曲的这个荷载完全可以让结构发生非线性屈曲.那么我们就把线性屈曲分析失稳时的deflection缩小(乘以一个小于1的数),所为进行非线性屈曲分析时对结构初始缺陷的考虑.需要介绍的时,这个方法,是进行二阶计算的一个简化方法.另外一个二阶计算方法考虑的模型是刚塑性分析(把节点考虑为发生塑性变化,成为塑性铰,而结点以外梁柱其他地方仍然认为是刚性).
我们一般用非线性屈曲分析,和线性屈曲分析来进行判断求丫的分叉点,和类似正弦图象的最高点的值.
非线性屈曲分析是进行倒结构的限制荷载或最大荷载结束.分析中包含了塑性非线性的问题.非线性屈曲分析考虑了结构的初始缺陷问题,结构比特征值的屈
曲分析精确,是可以用在实际工程中的.
随动强化假定屈服面的大
小保持不变，而仅仅宰屈服的方向上移动，某方向的屈服应力升高，相反方向的屈服应力降低。
等向强化是屈服面以材料中所作塑性功的大小为基础宰尺寸上扩张。对于Von Mises屈服准则来说，屈服面宰所有方向上均匀扩张。
我们说屈曲分析是研究结构或构建的平衡状态是否稳定的问题.处于平衡位置的结构或构建在任意微小的外界扰动下,将偏离平衡位置,当扰动出去后,又恢复到平衡位置,这说明处置的平衡位置是稳定的,比若说小时候玩的不倒嗡,他最后还是会树立起来,可以相似的这么理解.如果不能回到初始的平衡位置,则说他是不稳定的,从初始平衡位置转变到另一个平衡位置,成为屈曲或者失稳.你可以这么想象,和人一样高的两个木桩放在水平地上,一个想手指头一样细,一个想沙发一样大的横截面,你说我对他们各踢一脚,谁会倒下去?但注意,这个时候他们都是完好的,我踢一脚,不能让他们损坏,但是可以让他失稳---倒下去.
非线性有限元9弹塑性本构关系ppt课件
对塑性变形基本规律的认识来自于实验： • 从实验中找出在应力超出弹性极限后材料的特性; • 将这些特性进行归纳并提出合理的假设和简化模型，
确定应力超过弹性极限后材料的本构关系; • 建立塑性力学的基本方程; 1) 求解这些方程，得到不同塑性状态下物体内的应力和
应变。
• 塑性阶段：继续加载，材料可承受 更大应力，称为材料强化，并伴随 出现塑性应变。至A点以前卸载， 路径接近直线，即处于弹性卸载状 态，其斜率等于加载斜率E。
1) 破坏点：继续加载至可承受的最大 极限应力，试件出现颈缩而破坏，
称为强度极限。
在日常生活中，随处都可以看到浪费 粮食的 现象。 也许你 并未意 识到自 己在浪 费，也 许你认 为浪费 这一点 点算不 了什么
1913年：泰勒（Taylor）的实验证明，LevyMises本构关系是真实情况的一阶近似。
1924年：提出塑性全量理论，伊柳辛（Ilyushin） 等苏联学者用来解决大量实际问题。
1930年：罗伊斯（Reuss）在普朗特（Prandtle） 的启示下，提出包括弹性应变部分的三维塑性应力 -应变关系。至此，塑性增量理论初步建立。
（屈服点），描写多维问题的屈服条件就需要应力或应变空间的一个临界曲面，该
曲面称为屈服面。
考虑到塑性变形与静
水压力无关的特点
f1,2,3C
FJ2,J3C
至今已出现许多屈服理论。俞茂宏教授在这方面做出了重要贡献。 屈服函数：
是描写屈服条件的函数。不同屈服条件，其屈服函数不尽相同。
在日常生活中，随处都可以看到浪费 粮食的 现象。 也许你 并未意 识到自 己在浪 费，也 许你认 为浪费 这一点 点算不 了什么
基本实验有两个： • 简单拉伸实验：实验表明，塑性力学研究的应力与应变
高层建筑结构抗震弹塑性分析方法及抗震性能评估的研究
高层建筑结构抗震弹塑性分析方法及抗震性能评估的研究一、本文概述本文旨在探讨高层建筑结构在地震作用下的弹塑性分析方法及其抗震性能评估。
地震是自然界中常见的灾害性事件，对人类社会和建筑结构产生深远影响。
高层建筑由于其特殊的结构特点和高度，使其在地震中更容易受到破坏。
因此，研究高层建筑结构的抗震性能，特别是在弹塑性阶段的分析和评估，对于提高建筑结构的抗震能力，减少地震灾害损失具有重要意义。
本文将首先介绍高层建筑结构抗震弹塑性分析的基本理论和方法，包括弹塑性力学基础、结构分析模型、地震动输入等。
在此基础上，探讨高层建筑结构在地震作用下的弹塑性响应特点，包括结构变形、内力分布、能量耗散等。
然后，本文将重点介绍高层建筑结构抗震性能评估的方法和技术，包括静力弹塑性分析、动力弹塑性分析、易损性分析等。
这些方法和技术可以用于评估高层建筑结构在地震中的安全性能和抗震能力。
本文还将对高层建筑结构抗震弹塑性分析方法和抗震性能评估的应用进行案例研究。
通过实际工程案例的分析，探讨不同分析方法和技术在实际工程中的应用效果，为高层建筑结构的抗震设计和评估提供参考和借鉴。
本文将对高层建筑结构抗震弹塑性分析方法和抗震性能评估的未来发展趋势进行展望，提出相关的研究建议和展望。
通过本文的研究，可以为高层建筑结构的抗震设计和评估提供更为科学、合理的方法和技术支持，有助于提高高层建筑结构的抗震能力，减少地震灾害损失。
二、高层建筑结构抗震弹塑性分析方法的研究高层建筑结构的抗震弹塑性分析是评估建筑在地震作用下的响应和性能的重要手段。
随着建筑高度的增加，结构的柔性和非线性特性愈发显著，因此，采用弹塑性分析方法可以更准确地模拟结构在地震中的实际行为。
材料本构关系的研究：高层建筑的抗震性能与其组成材料的力学特性密切相关。
研究材料在循环加载下的应力-应变关系、滞回特性以及损伤演化规律，是弹塑性分析的基础。
通过试验和数值模拟，可以建立更精确的材料本构模型，为结构分析提供数据支持。
XFEM在弹塑性断裂力学中的应用
XFEM在弹塑性断裂力学中的应用基础知识讲解弹塑性断裂力学简介13XFEM 在裂纹扩展中的应用扩展有限元（XFEM ）发展现状2Abaqus 中XFEM 功能的实现4•线弹性断裂力学是建立在线弹性力学基础上的，它没能考虑裂纹尖端附近塑性性区的影响，因而只适用于高强度（钢）脆性材料，对于工程中大量使用的中、低强度钢等具有较好塑性的材料是不适用的。
为了将应力强度因子推广到裂纹尖端有小范围塑性区的情况，人们推出了应力强度因子塑性区的修正方法，但适用性并不理想。
为了研究塑性材料的断裂问题，又产生了断裂力学的另一个分支——弹塑性断裂力学。
•弹塑性断裂力学研究范围：•（1）大范围屈服：端部的塑性区尺寸接近或超过裂纹尺寸，如：中低强度钢制成的构件。
•（2）全面屈服：材料处于全面屈服阶段，如：压力容器的接管部位。
•弹塑性断裂力学的任务：确定能定量描述裂纹尖端区域弹塑性应力，应变场强度的参量。
以便利用理论建立起这些参量与裂纹几何特性、外加载荷之间的关系，通过试验来测定它们，并最后建立便于工程应用的断裂准则。
•主要包括COD理论和J积分理论。
•1948年Irwin 和Orowan 将塑性应变能引入能量理论作为塑性材料裂纹的能量判据标志着弹塑性断裂力学研究的开始。
•1960年，Dugdale 建立了研究裂纹尖端塑性区的D-M 模型。
•1965年Wells 提出COD 准则：,其中为裂纹尖端张开位移，为开裂临界值，是由实验测得的材料常数，表征了材料的弹塑性断裂韧性。
但是裂纹开裂后，材料在到达失稳点并失效破坏前还可以继续承受更多的载荷，即裂纹在达到开裂临界状态后还有一定的承载能力，因此以为指标进行设计是偏于安全的。
•1968年，Rice 提出了J 积分理论建立了J 积分断裂准则。
•1968年，Hutchinson ，Rosengren 与Rice 提出了建立在塑性力学全量理论基础上的HRR 理论，其理论基础是J 积分，为弹塑性断裂力学奠定了理论基础。
弹塑性力学基础与有限元分析-接触分析实例
06
结论与展望
结论
1
本文通过理论分析和有限元模拟，深入研究了弹 塑性力学基础与有限元分析在接触分析中的应用。
2
研究结果表明，弹塑性力学基础与有限元分析在 接触分析中具有较高的精度和可靠性，能够有效 地模拟复杂接触问题。
3
本文所采用的有限元分析方法在处理接触问题时 具有较好的通用性和扩展性，为进一步研究复杂 接触问题提供了有力支持。
弹塑性本构模型
弹塑性本构模型的定义
弹塑性本构模型是描述弹塑性材料力学行为的数学模型，它通过应力应变关系来描述材料的弹塑性行 为。
常见的弹塑性本构模型
常见的弹塑性本构模型包括Mohr-Coulomb模型、Drucker-Prager模型、Cam-Clay模型等。这些模 型在描述材料的弹塑性行为方面各有特点，适用于不同的材料和工程问题。
接触面完全贴合，无相对运动。
滑动状态
接触面部分贴合，存在相对运动。
混合状态
接触面同时存在分离、粘结和滑动。
接触检测与跟踪
初始接触检测
确定初始状态下接触面的位置和状态。
接触状态跟踪
实时监测接触面的运动状态和相互作用。
接触面更新
根据接触状态调整接触面的几何形状和参数。
接触刚度与阻尼
1 2
接触刚度
描述接触面间的相互作用力与相对位移的关系。
求解阶段主要进行有限元 方程的求解，得到各节点 的位移和应力等结果。
ABCD
前处理阶段主要完成有限元 模型的建立和网格划分，为 求解阶段提供输入数据。
后处理阶段主要对求解结果进 行可视化、分析和评估，为工 程设计和优化提供依据。
04
接触分析原理
接触状态描述
分离状态
非线性弹性力学
1948年R.S.里夫林在任意形式的贮能函数下，得到不可压缩弹性体的几个简单而重要问题的精确解。将它们 应用于橡胶制品，即使橡胶的伸长为原长的两三倍，精度仍能达到百分之几。在这一成就的鼓舞下，学者们重新 开始探讨有限变形弹性理论，并导致了整个的蓬勃发展。此后，非线性弹性理论就成为理性力学的重要组成部分。 1952年起C.特鲁斯德尔、W.诺尔、B.D.科勒曼、J.L.埃里克森、M.E.格廷、A.C.爱林根以及美籍华人王钊诚在 非线性弹性力学方面作出较大贡献，中国的郭仲衡于1962～1963年连续发表了多篇论文。1972年奥登等人在用有 限元法进行数值解方面做了大量有成效的工作，从而使得非线性弹性力学在工程实际中得到较广泛的应用。但是 非线性弹性力学无论在理论方面、精确解方面还是数值近似解方面都比线性弹性力学难度大，所以至今远不如线 性弹性力学成熟，有许多问题尚需进一步探讨。非线性弹性力学的基本概念和方程比较复杂，在分析中大多采用 张量这一数学工具。
变形描述
变形描述在讨论非线性弹性力学问题时，取初始时刻物体在三维空间中所占的区域为参考构形（见）现时构形，在其上取笛卡儿坐标。
由方程 对于有单值逆变换的情形，存在 在时刻物质点的位置矢量为X，在运动过程中，该点在时刻的位置矢量为，则 在时刻物质点的位置矢量为X，在运动过程中，该点在时刻的位置矢量为，则 其中u是该物质点的位移矢量,它在和中的坐标分别记为和。 必须区分使用和坐标，这是非线性弹性力学区别于线性弹性力学的基本特征之一。 描述物体变形的量有变形梯度，在中，其定义为： 其中为克罗内克符号；为位移分量的偏导数，即变形梯度既包含纯变形又包含刚性转动，为把纯变形从其中 分解出来，须采用极分解定理，相应于左分解和右分解分别得到左柯西-格林应变(又称芬格应变)和右柯西-格林 应变(又称格林变)。而在中有逆应变(称为皮奥拉应变)和(称为柯西应变)。
三维弹塑性问题的比例边界有限元法
04
比例边界有限元法的实现 过程
网格划分与节点生成
网格划分
将三维空间离散化为有限个小的单元，每个单元由节点连接。
节点生成
根据几何形状和边界条件，在关键区域布置节点，确保计算的精确性。
比例边界条件的处理
边界条件转换
将比例边界条件转换为等效的节点力约束。
节点力平衡
确保所有节点力在平衡状态下，以实现真实比例边界条件的模拟。
材料属性
根据实际问题，设置材料属性，如弹性模量、泊松比、密度等 。
力学行为
考虑弹性和塑性行为，建立相应的本构关系和屈服条件。
边界条件与载荷施加
边界条件
根据实际问题，施加相应的边界条件，如固定边界的位移约束、滑动边界的 摩擦力约束等。
载荷施加
根据实际问题，施加相应的外部载荷，如重力、压力、扭矩等。同时考虑惯 性效应，如质量、阻尼等。
三维弹塑性问题的有限元 建模
有限元模型的建立
几何模型
根据实际物理模型，建立相应 的几何模型，包括三维实体、
表面等。
网格划分
根据模型复杂程度和计算精度要 求，选择合适的网格类型和密度 进行划分。
边界定义
根据实际问题，定义模型的边界条 件，如固定边界、滑动边界等。
单元选择与属性设置
单元类型
根据实际问题，选择合适的有限元单元类型，如四面体单元、 六面体单元等。
三维弹塑性问题的比例边界 有限元法
2023-11-06
目 录
• 引言 • 三维弹塑性理论基础 • 三维弹塑性问题的有限元建模 • 比例边界有限元法的实现过程 • 三维弹塑性问题的算例分析 • 结论与展望 • 参考文献
01
引言
研究背景与意义
材料非线性有限元分析
e p 1 d ij d ij d ij Dijkl d kl f , ij d
A f , p Dijkl f , kl M
ij
dij ( D
dij ( D
1 ijkl
H (l ) f , ij f , kl )d kl Dep1,ijkl d kl A
1 J 2 [( 1 2 ) 2 ( 2 3 ) 2 ( 3 1 ) 2 ] 6 2 2
随动强化的米塞斯屈服准则
这种材料的屈服面方程为
p ij 1 1 p p f ( sij , , k ) [ ( sij ij )( sij ij )] 2 0 0 2
kk pp
纯剪
单向拉伸
Gp是塑性剪切 模量
Ep是塑性拉伸 模量
A f , p Dijkl f , kl M
ij
f , kl Dijkl f , kl
s ij G 2J 2G s ij G 2 2 2
由此可得A=G+Gp（或A=G+Ep/3），又因 1 1G G p Dijkl Dijkl f , kl Dklij f , kl sij skl A A 1 G G G 2 s ij s ij s kl s 2 kl G Gp G Gp p p d ij Dijkl d kl 由此可得弹塑性矩阵为
J 2 sij sij / 2
，因此
由于偏张量第一不变量＝0
J1 sii 0
1 2 2 2 J 2 [( s11 s22 ) 2 ( s22 s33 ) 2 ( s33 s11 ) 2 ] s12 s23 s31 6
弹塑性力学与有限元
x yx zx X 0 x y z
《弹塑性力学与有限元》
应力分析
平衡微分方程
x xy xz 2u X 0( 2 ) x y z t yx y yz 2v Y 0( 2 ) x y z t zx zy z 2w Z 0( 2 ) x y z t
《弹塑性力学与有限元》
应变分析
应变—位移关系(几何方程)
根据泰勒级数展开式，可得：
f 1 ( x , y , z ) 1 2 f 1 ( x, y , z ) 2 u1 f 1 ( x , y , z ) dx dx 2 x 2! x
略去高阶项后得到：
u u1 u dx x
《弹塑性力学与有限元》
应变分析
应变—位移关系(几何方程)
同理可得另外两个剪应变 xy， yz ,即有剪应变的表达式:
xy
yz
zx
u v y x
v w z y u w z x
说明：剪应变的正负号
ij 0(i, j x, y , z )表示夹角变小 ij 0(i, j x, y , z )表示夹角变大
应变分析
应变—位移关系
位移—由于外部因素如载荷或温度
变化,物体内部各点空间位置发生的
变化 ；
如果各点的位移完全相同，物体发
生刚体平移；如果各点的位移不同， 但各点间的相对距离保持不变，物 体发生刚体转动等刚体移动；
《弹塑性力学与有限元》
应变分析
应变—位移关系
连续体内如果各点(或部分点）间的相对距离发生变化， 则物体发生了变形，这时的位移是变形体位移。此物体 被称为有变形或有应变。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题 材料非线性问题
➢ 塑性力学的基本法则 (i) Prager运动硬化法则 规定加载曲面中心的移动是在表征现时应力状态的应力点的法线方向。
Prager运动法则一般说只能应用于九维应力空间。
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题 材料非线性问题
(3)按单元内各个积分点计算D的预测值
1)计算屈服函数值
，然后区分三种情况
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题
材料非线性问题
➢ 弹塑性增量分析数值方法中的几个问题 弹塑性状态的决定和本构关系的积分 (i)
(ii) 若
，则该积分点为由弹性
进入塑性的过渡情况，计算比例因子m。
（iii）若
二． 应力的度量
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题
二． 应力的度量
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题
➢ 大变形情况下的本构关系
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题
➢ 大变形情况下的本构关系
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题
➢ 大变形条件下的应变和应力的度量 一． 应变的度量
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题
几何非线性问题
➢ 大变形条件下的应变和应力的度量 二． 应力的度量 在大变形问题中，是从变形后的物体中截取出微元体建立平衡方 程和与之相等效的虚功原理，所以应从变形后的物体内截取单元 体定义应力张量－－欧拉应力张量，tτij
➢ 大变形情况下的本构关系
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题
➢ 大变形情况下的本构关系
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题
➢ 几何非线性问题分析的有限元法表达格式
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题
➢ 几何非线性问题分析的有限元法表达格式
《弹塑性力学与有限元》
➢ 弹塑性增量的应力应变关系
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题 材料非线性问题
➢ 弹塑性问题的增量方程 将载荷分成若干个增量，然后对于每一载荷增量，将弹塑性方程 线性化。假设对于时刻t的解已经求得，要求解t+Δt时刻的解。
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题
材料非线性问题
如令M=1或M=0，混合硬化法则就分别蜕化为各向同性硬化法则和运 动硬化法则运动硬化法则。
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题 材料非线性问题
➢ 塑性力学的基本法则 4．加载、卸载准则
(i) 对于理想弹塑性材料，此情况是塑性加载。 (ii) 对于硬化材料，此情况是中性变载。
《弹塑性力学与有限元》
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题 材料非线性问题
➢ 材料弹塑性本构关系 材料弹塑性行为的描述 1．单调加载
2．反向加载
3．循环加载
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题 材料非线性问题
➢ 塑性力学的基本法则 1．初始屈服条件
(1) V．Mises条件
其中
有
《弹塑性力学与有限元》
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题
几何非线性问题
➢ 引言 ➢ 大变形条件下的应变和应力分量 ➢ 大变形条件下的本构关系 ➢ 几何非线性问题的表达格式 ➢ 有限元求解方程及解法
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题 几何非线性问题
➢ 引言
问题的类型： 大位移、小应变问题 大应变问题
令m=0
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题
作业：有一钢开孔板，其形状如下图所示，各尺寸a=1.6m ，b=0.8m ， r=0.05m，板厚 t=0.01m。设钢材的弹性模量为206Gpa ，屈服应力为235Mpa ，泊松比0.31 。 1、选取适当的单元利用ANSYS进行建模，并说明采用的单元、材料类型。 2、假设钢材为线性强化的弹塑性材料，其屈后刚度为弹性模量的1/10。对钢板进行位移 控制的往复加载。共有四级荷载，每一级先由ux=0 正向加载至位移为uxn ，再卸载并反向 加载至-uxn，最后回到ux=0, 完成一级的加载循环，四级荷载的位移值 uxn分别为0.5mm， 1.0mm，1.5mm，2.0mm。利用双线性随动强化和双线性等向强化两种模型进行分析并 绘制力-位移曲线。
《弹塑性力学与有限元》
弹塑性力学与有限元 —材料非线性问题和几何非线性问题
材料非线性问题和几何非线性问题 材料非线性问题和几何非线性问题 材料非线性问题 几何非线性问题
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题 材料非线性问题
➢非线性方程组的数值解法 ➢材料弹塑性本构关系 ➢弹塑性增量有限元分析 ➢弹塑性增量分析数值方法中的几个问题
如果应变是用变形前的坐标表示的 Green应变张量，则需要定义与之对 应的关于变形前位形的应力张量。
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题 几何非线性问题
二． 应力的度量
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题
二． 应力的度量
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题
➢ 弹塑性问题的增量方程 它们应满足的方程和边界条件是
应力应变关系应通过积分求得
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题 材料非线性问题
➢ 增量有限元格式 首先建立增量形式的虚位移原理如下：
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题 材料非线性问题
➢ 增量有限元格式 基于增量形式虚位移原理有限元表达格式
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题
几何非线性问题
➢ 引言 几何非线性问题： 板、壳等薄壁结构在一定载荷作用下，尽管应变很小，甚至未超过 弹性极限，但是位移较大。这时必须考虑变形对平衡的影响，即平 衡条件必须建立在变形后的位形上，同时应变表达式应包括位移的 二次项——平衡方程和几何条件都是非线性的； 金属成型材料在受载时都可能出现很大的应变，这时除了采用非线 性的平衡方程和几何关系外，还需要引入相应的应力应变关系。
材料非线性问题
➢ 弹塑性增量分析数值方法中的几个问题 非线性方程组的求解方案 1．欧拉法及其改进 2．变刚度迭代(N-R迭代) 迭代步骤： (1) 形成方程组 (2) 求解方程组，得到本次迭代的位移增量修正量
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题
材料非线性问题
➢弹塑性增量分析数值方法中的几个问题 (3) 计算各单元应变增量和应力增量修正量
➢ 塑性力学的基本法则 (ii) Zeigler修正运动硬化法则 规定加载曲面沿联结其中心和现时应力点的向量方向移动。
在九维应力空间，以及在包括三个正应力，或不包括任何正应力的 应力子空间，这两种法则是完全相同的。
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题 材料非线性问题
➢ 塑性力学的基本法则 (3) 混合硬化法则
材料非线性问题和几何非线性问题
➢ 几何非线性问题分析的有限元法表达格式
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题
➢ 几何非线性问题分析的有限元法表达格式 • 完全拉格朗日格式
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题
➢ 几何非线性问题分析的有限元法表达格式 • 更新拉格朗日格式
材料非线性问题和几何非线性问题 材料非线性问题
➢ 塑性力学的基本法则 (2)Tresca条件
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题
材料非线性问题
➢ 塑性力学的基本法则 2．流动法则 可从塑性势导出流动法则
对于关联塑性情况，流动法则表示为 3．硬化法则
(1) 各向同性硬化法则 (2) 运动硬化法则
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题
几何非线性问题
➢ 引言 在几何非线性问题的有限单元法中，通常采用增量分析方法。 增量分析方法一般采用两种表达格式： 完全的Lagrange格式：静力学和运动学变量总是参考初始位形， 即整个分析过程中参考位形保持不变。 更新的Lagrange格式：静力学和运动学变量参考于每一载荷或时 间步长开始时的位形，即在分析过程中参考位形不断在更新。
材料非线性问题和几何非线性问题
材料非线性问题
➢ 弹塑性增量的应力应变关系 1. 建立弹塑性增量的应力应变关系需遵循的原则 (1) 一致性条件 (2) 流动法则 (3) 弹性应力应变关系 2. 各向同性硬化材料的应力应变关系 以各向同性硬化材料为例
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题 材料非线性问题
(4) 根据收敛准则检验解是否满足收敛要求
3. 常刚度迭代(mN-R迭代)
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题
材料非线性问题
➢弹塑性增量分析数值方法中的几个问题 弹塑性状态的决定和本构关系的积分 决定弹塑性状态的一般算法步骤: (1)利用几何关系计算应变增量 (2)按弹性关系计算应力增量的预测值以及应力的预测值
《弹塑性力学与有限元》
材料非线性问题和几何非线性问题
材料非线性问题
➢ 弹塑性增量有限元分析 每一增量步包含下列三个算法步骤： 1．线性化弹塑性本构关系，并形成增量有限元方程。 2．求解有限元方程。 3．积分本构方程决定新的应力状态，检查平衡条件，并决定 是否进行新的迭代。
《弹塑性力学与有限元》