弹塑性问题有限元分析

合集下载

弹塑性有限元法

弹塑性有限元法
第六章 弹塑性有限元法
当变形体同时存在大的弹性和塑性变形时,必 须采用弹塑性力学进行分析,相应的有弹塑性 有限元法,其较一般弹性有限元复杂得多。
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时刻的变形为

弹塑性力学土木工程应用有限元ABAQUS分析课件

弹塑性力学土木工程应用有限元ABAQUS分析课件

A
A0
l0 l
l 0 未变形的长度 A 0 未变形的平面面积
FF l
A A0 l0
nom(ll0)
nominal
n o m 名义应力
真实应力
弹塑性力学土木工程应用 有限元ABAQUS分析
名义、真实应力(变) 名义应变,每单位未变形长度的伸长。
noml0l
ll0 l0
l l0
1
l l0
1 nom
塑性性能的材料实验数据,提供的应变包括塑性应变和弹性应 变,是材料的总体应变。所以总体应变分解为弹性和塑性应变两 项。
弹性应变等于真实应力与弹性模量的比值。
t pl el
el / E
p lte lt/E
p l 真实塑性应变
t 总体真实应变
弹塑性力学土木工程应用 有限元ABAQUS分析
l0d lllnll0
lnl lnl0l
l0
l0
nom
l l0
lnl0 l0lln1nom
弹塑性力学土木工程应用 有限元ABAQUS分析
名义、真实应力(变) 真实应力与名义应力的关系
nom(1nom)
真实应变与名义应变的关系
ln1nom
弹塑性力学土木工程应用 有限元ABAQUS分析
名义、真实应力(变)
弹塑性力学的发展
早期 精确算法 线性问题
如今 数字分析法 非线性问题
实际的需要,软件应用计算 ANSYS、ABAQUS
弹塑性力学土木工程应用 有限元ABAQUS分析
PART.02
名义应力(变)与真实应力(变)
弹塑性力学土木工程应用 有限元ABAQUS分析
名义、真实应力(变)
在ABAQUS中必须 用真实应力和真实应 变定义塑性。

基于有限元的弹塑性裂纹数值分析

基于有限元的弹塑性裂纹数值分析

就 属 于 大 范 围屈 服断 裂 问 题 。 此 时 线弹 般 脆 性 金属 材 料 , 如 铸 铁 等 在裂 纹 扩 展 前 , 其 端部 都 将 出现 甚 至 超 过 裂 纹 本 身 的尺寸, 个 塑性 区 。 当此 塑性 区 尺寸 很 小 , 即远 小 于 裂 纹 时 , 线 弹性 断 裂 性 断 裂力 学理 论 已 不再 适用 , 用 应 力强 度 因子衡 量 裂 尖 应 力场 强度 这 种 塑性 变形 占较 大 比 重 的断 裂问 题 就要 用弹 塑性 断 力 学 仍 有 足够 的 精度 。 此 类 断 裂 称 为小 范 围屈 服 断 裂 , 可 以 采 用 对 将 失 去意 义 。 目前广为应 用的 是 C OD 原理. J 积分 理论 等方 法… 。 线 弹 性 力学 导 出的 应 力 强 度 因子 进行 修正 的 方 法 来 处 理 。 然 而 对 裂 理论 来 解决 ,

Q : §
工 程 技 术
Sc i e nce an d Tec h no l o gy I n no vat i on Her al d
基于 有 限 元的弹 塑 性 裂纹 数 值分 析
邢文金 ( 吉林大学机械工程 学院 吉林长春
1 3 0 0 0 0 )
摘 要 : 在 线弹性断裂力学和D - M 模型的基础上 , 推导出了 受单向拉仲含 中心穿透裂纹的理 想弹塑性材料J 积分的解析式- 通过A N S Y S 对弹塑 性J 积分进行数值计算, 与推导 出的解析解此较 , 表明 了 用 有限元方法计算弹塑性J 积分具有相当高的精度; 分析 了 J 积分与裂纹 初始 长度及 外荷 或的关系, 对理想弹塑性 材料塑性区大小进行 了 探讨, 结果表明, 塑性 区尺寸随外荷戴增大而增大, 并且外荷裁接近屈服应 力时, 裂纹塑性区尺

塑性成形过程中的有限元法

塑性成形过程中的有限元法

塑性成形过程中的有限元法金属塑性成形技术是现代化制造业中金属加工的重要方法之一。

它是金属材料在模具和锻压设备作用下发生变形,获得所需要求的形状、尺寸和性能的制件的加工过程。

金属成形件在汽车、飞机仪表、机械设备等产品的零部件中占有相当大的比例。

由于其具有生产效率高,生产费用低的特点,适合于大批量生产,是现代高速发展的制造业的重要成形工艺。

据统计,在发达国家中,金属塑性成形件的产值在国民经济中的比重居行业之首,在我国也占有相当大的比例。

随着现代制造业的快速发展,对塑性成形工艺分析和模具设计提出了更高的要求。

如果工艺分析不完善、模具设计不合理或选材不当,产品将不符合质量要求,导致大量不良品和废品,增加模具的设计制造时间和成本。

为了防止缺陷,提高产品质量,降低产品成本,国内外许多大公司、企业、高校和研究机构对塑料成型件的性能进行了大量的理论分析、实验研究和数值计算,通过对成形过程中应力应变分布及变化规律的研究,试图找出各零件在产品成形过程中遵循的共同规律和机械失效所反映的共同特征。

由于影响塑性成形过程的因素很多,一些因素,如摩擦和润滑、变形过程中材料的本构关系等,还没有被人们充分理解和掌握。

因此,到目前为止,还无法对各种材料和形状零件的成形过程做出准确的定量判断。

由于大变形机理非常复杂,塑性成形研究领域一直是一个充满挑战和机遇的领域。

一般来说,产品研究与开发的目标之一就是确定生产高质量产品的优化准则,而不同的产品要求不同的优化准则,建立适当的优化准则需要对产品制造过程的全面了解。

如果不掌握诸如摩擦条件、材料性能、工件几何形状、成形力等工艺参数对成形过程的影响,就不可能正确地设计模具和选择加工设备,更无法预测和防止缺陷的生成。

在传统工艺分析和模具设计中,主要还是依靠工程类比和设计经验,经过反复试模修模,调整工艺参数以期望消除成形过程中的产品缺陷如失稳起皱、充填不满、局部破裂等。

仅仅依靠类比和传统的经验工艺分析和模具设计方法已无法满足高速发展的现代金属加工工业的要求。

弹塑性力学与有限元-应变分析

弹塑性力学与有限元-应变分析
位移—由于外部因素如载荷或温度 变化,物体内部各点空间位置发生的 变化 ;
如果各点的位移完全相同,物体发 生刚体平移;如果各点的位移不同, 但各点间的相对距离保持不变,物 体发生刚体转动等刚体移动;
《弹塑性力学与有限元》
应变分析
应变—位移关系
连续体内如果各点(或部分点)间的相对距离发生变化, 则物体发生了变形,这时的位移是变形体位移。此物体 被称为有变形或有应变。
《弹塑性力学与有限元》
应变分析
主应变和主剪应变
I1 x y z
I 2
x y
y z
z x
2 xyБайду номын сангаас
2 yz
2 zx
x
y
y
z
z
x
1 4
(
2 xy
2 yz
2 zx
)
I
3
x
y z
2 xy yz zx
(
x
2 yz
y
2 zx
z
2 xy
)
x
y z
1 4
xy
yz
zx
1 4
(
x
2 yz
个 Mohr圆一起沿 轴平移一个距离
,该距离等于所叠加的静水应力,
O P3 O M P2 s3
P1
并不改变Mohr圆的大小。
➢ τ轴的位置与屈服及塑性变形无关 ,决定屈服与塑性变形的只是Mohr 圆本身的大小。
m
s2
s1
图 3-4
《弹塑性力学与有限元》
应力分析
应力的Mohr圆
若将τ轴平移到O' ,并使
应变分析
应变—位移关系(几何方程)
同理可得另外两个剪应变 xy, yz ,即有剪应变的表达式:

弹塑性有限元法基本理论与模拟方法

弹塑性有限元法基本理论与模拟方法
流体动力学
用于模拟流体流动和传热问题 ,如流体机械、航空航天和化 工等领域。
电磁场
用于分析电磁场问题和电气设 备性能,如电机、变压器和天 线等。
声学
用于模拟声音传播和噪声控制 问题,如声学器件和声学环境
等。
04 弹塑性有限元法的基本原 理
弹塑性有限元法的离散化方法
有限元离散化
将连续的物理场或结构体离散为有限个小的单元体, 每个单元体之间通过节点相互连接。
结构强度分析的模拟
结构强度评估
通过弹塑性有限元法模拟,可以对结构的强度进行评估,预测结构在不同载荷下的响应, 确保结构的安全性和稳定性。
疲劳寿命预测
利用弹塑性有限元法,可以模拟结构的疲劳载荷历程,预测结构的疲劳寿命,为结构的维 护和更换提供依据。
结构优化设计
通过模拟结构的应力分布和变形,可以优化结构设计,降低结构重量,提高结构效率。
边界条件和初始条件
在平衡方程中考虑边界条件和初始条件,以确保模拟的准确性和收 敛性。
弹塑性有限元法的边界条件和初始条件
边界条件的处理
01
根据实际情况,将边界条件转化为节点约束或单元载荷的形式。
初始条件的设置
02
在非稳态问题中,需要考虑初始条件的设置,以模拟问题的初
始状态。
边界条件和初始条件的实施
03
随着计算机技术的不断发展,弹塑性 有限元法在各个工程领域中得到了广 泛应用,如机械、航空航械设计中,弹塑性有限元法可用于分析各种复杂结构 的应力分布、变形和疲劳寿命等,提高产品的可靠性和安 全性。
航空航天
在航空航天领域,弹塑性有限元法可用于分析飞行器结构 在各种载荷下的响应,优化结构设计,提高飞行器的性能 和安全性。

第四章__弹塑性有限元法基本理论与模拟方法讲解

第四章__弹塑性有限元法基本理论与模拟方法讲解

Pi( k ) P(k 1) Pi(k )
q
(k ) 1
k) (k ) qi( k ) qi( q 1 i
q(k )
q( k 1)
第四章 弹塑性有限元法基本理论与模拟方法
(3) 所有载荷段循环,并将结果进行累加
第四章 弹塑性有限元法基本理论与模拟方法
4.2 材料非线性问题及分类
为了与初始屈服应力相区别,我们称之为后继屈服应力。 与初始屈服应力不同,它不是一个材料常数,而是依赖 于塑性变形的大小和历史。 后继屈服应力是在简单拉伸下,材料在经历一定塑性变形 后再次加载时,变形是按弹性还是塑性规律变化的界限。
第四章 弹塑性有限元法基本理论与模拟方法 第四章 弹塑性有限元法基本理论与模拟方法
第四章 弹塑性有限元法基本理论与模拟方法
F
ห้องสมุดไป่ตู้ nom
s0
F nom A0 L nom L0
L
nom

s0
nom (1 nom ) p e ln(1 nom ) E
x1 x x 2 , xn
F(x)=0
f1 (x) f ( x) F ( x) 2 , f n ( x)
0 0 0 0
第四章 弹塑性有限元法基本理论与模拟方法 第四章 弹塑性有限元法基本理论与模拟方法
和简单应力状态相似,材料在复杂应力状态下同样 存在初始屈服和后继屈服的问题。
材料在复杂应力状态下,在经历初始屈服和发生塑性 变形后,此时卸载,将再次进入弹性状态(称为后继弹 性状态)。
第四章 弹塑性有限元法基本理论与模拟方法 第四章 弹塑性有限元法基本理论与模拟方法

塑性材料的有限元分析

塑性材料的有限元分析

针对复杂材料和结构,需要深入研究材料的非线 性行为和多场耦合效应,建立更加完善的物理模 型和数值算法。
此外,应加强与实验研究的结合,通过实验验证 和修正有限元模型,提高模拟结果的可靠性。同 时,实验研究也能够为有限元分析提供更加真实 和全面的材料性能数据。
THANK YOU
03
有限元分析方法
有限元分析的基本原理
离散化
将连续的物理系统离散为有限个小的单元,每个 单元称为有限元。
近似解法
通过数学方法求解每个有限元的近似解,再通过 组合所有有限元的解得到整个系统的近似解。
平衡方程
建立每个有限元的平衡方程,通过求解平衡方程 得到每个节点的位移和应力。
有限元分析的实现过程
然而,塑性材料的有限元分析仍存在 一些挑战和限制,如模型的简化、边 界条件的确定、材料参数的获取等, 需要进一步研究和改进。
研究展望
未来研究应致力于发展更加精确和高效的有限元 分析方法,提高模拟结果的可靠性和精度。
在实际工程应用中,应加强有限元分析与其他数 值方法(如边界元、有限体积等)的结合,实现 优势互补,提高计算效率。
塑性变形的微观机

塑性变形是通过位错的滑移和攀 移等微观机制实现的,这些机制 在宏观上表现为塑性变形。
塑性变形的温度效

温度对塑性变形的影响较大,温 度升高会使材料的屈服强度降低, 塑性变形能力增强。
塑性变形的加工硬

在塑性变形过程中,材料的屈服 强度会随着变形程度的增加而逐 渐提高,这种现象称为加工硬化。
背景
随着计算机技术的不断发展,有限元分析已成为工程领域中解决复杂问题的常 用方法。通过有限元分析,可以模拟材料的变形、应力分布、应变等,为实际 工程提供重要的理论依据。
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
弹塑性问题的有限元分析
专硕-
1
材料的弹塑性行为实验
2
材料塑性行为的屈服准则
3
材料塑性行为的流动法则
4
材料塑性行为的强化准则
5
材料塑性行为的模型
研究弹塑性问题的关键在于物理方程的处理。下面主要讨论小 变形情形下的弹塑性问题。
1、材料的弹塑性行为实验
典型的材料性能实验曲线是通过标准试样的单向拉伸与压缩获 得的,如下图所示
但不发生新的塑性流动
4、塑性强化准则 该准则用来描述屈服面是如何改变的,以确定后续屈服面的新 状态,一般可以有几种模型: 等向强化模型 随动强化模型 混合强化模型 5、材料塑性行为的模型 基于以上准则,在根据各种材料的应力应变曲线、经过归纳和 分类给出以下几种典型的描述材料弹塑性行为的模型 (1)、双线性Bauschinger随动强化 (2)、多线性Bauschinger随动强化 (3)、双线性等向强化 (4)、多线性等向强化 (5)、非等向强化 (6)、Drucker-Prager模型 所谓Bauschinger效应为反向屈服点到卸载点的数值为 2 yd 。
I1 1 2 3
I2 1 2 2 3 31(2)
I3 1 2 3
基于主应力空间,由等倾面组成的八面体的平面上的正应力和剪应力具有
一些特殊的性质。
设某一点的应力状态为 ij ,其中三个主应力为 1、 2、 3 ,并且1> 2> 3
如果坐标轴与主方向重合,则应力不变量如式(2)
其中 yd 为临界屈服剪应力,将由实验来确定,一般通过单拉实
验获得,由于单拉实验获得的是临界屈服拉应力 yd ,所以通过
以下关系来换算:
如果定义等效应力为
eq
3 2
y
d
yd
2 3
yd
1 2
(
xx
yy )2
(yyzz)2(zz
xx )2
6(
2 xy
2 yz
2 xz )
则初始屈服条件可以写成
eq yd 将等效应力写成更一般的形式,有
eq
f(

ij
则屈服函数为
F( ij) f ( ij ) yd 0 各种屈服面函数详情可见P319
3、塑性流动法则 该法则用来确定塑性应变分量在塑性变化时的大小和方向,相应的分 量沿着一个势函数的法线方向增长即:
d ij pl
Q ij
其中 dij pl为塑性应变增量, 为塑性增长乘子,Q为塑性势函数,若
弹性 极限
应 力
加 载
卸 载
塑性应变 弹性应变
断裂 应变
在实际结构中,真实的情况是材料处于复杂 的受力状态,ij 即中 的各个分量都存在,如何基 于材料的单拉应力-应变实验曲线,来描述复杂 应力状态下材料的真实弹塑性行为,就必须涉及 屈服准则、塑性流动法则、塑性强化法则这三个 方面的描述,有了这三个方面的描述就可以完全 确定出复杂应力状态下材料的真实弹塑性行为
PPT模板下载:www.1p pt.co m/ moban / 节日PPT模板:www.1p pt.co m/ jieri/ PPT背景图片:www.1p pt.co m/ beijing / 优秀PPT下载:www.1p pt.co m/ xiazai/ Word教程: /word/ 资料下载:www. 1ppt.co m/zilia o/ 范文下载:www. 1ppt.co m/fan wen/ 教案下载:www. 1ppt.co m/jiao an/
n( 2 3)
由于
n p1nx p2ny p3nz 1n2 x 2n2 y 3n2( z 4)
由(3)、(4)可得 :
n
2 1
n
2
x
22n2 y
32n2z
( 1n2 x
2n2 y
3n2z )2
由于静水压力是三个主应力的平均值,也就是说在该斜面上,其 正应力为
px xxnx yxny zxnz
同理可得
py yyny yxnx yznz pz zz nz zxnx zyny
px nx n , py ny n , pz nz n
代入上式可得:
nx ( xx n ) ny xy nx yz ny ( yy n )
为关联塑性流动,则Q就是屈服函数,即
Q( ij ) F ( ij )
当材料从一个塑性状态出发是继续塑性加载还是弹性卸载,通过以下
关系来判断:
如果F=0,并且 Q ij
d ij>0,则继续塑性加载
如果F=0,并且 Q
ij
d ij<0
,则继续弹性卸载
如果F=0,并且Qij d ij=0 ,则对于理想弹塑性材料 是塑性加载,对于硬化材料,此情况为中性变载,即继续为塑性状态,
nz nz
xz yz
0 0
nx zx
ny zy
nz ( zz
n)
0
这是关于nx , ny , nz的齐次线性方程组,其非零解的条件为行列式
等于零
展开可得:
n3
I1
2 n
I 2
n
I3
0(1)
其中
I1 xx yy zz
I2
xx
yy
xx zz
zz
yy
xy2
2 yz
2 zx
设该点有一斜面的应力矢量为p,它与 ij 保持平衡,该斜面的法线n的方
向为p余1 弦 为1nnxx、, pn2y、nz ,2n由y , 合p3 力 平3衡nz 可,以于得是到该p面在上坐的标与轴p方等向价的的三正个应投力影分n 和别剪
应力 n 的关系为:
2 n
p2
n2
2 1
nx
22ny
32nz
n 1/ 3(1 2 3() 5)
又为等倾面八面体
若要得到(5), 根据(4)可得:
nx ny nz 1/ 3
该八面体上的切应力为
8
1 3
( xx
yy )2
(
yy
zz )2
(
zz
xx )2
6(
2 xy
2 yz
2 xz )
它是决定材料是否屈服的力学参量,因此初始屈服条件为
8 yd
xx xy xz
I3 xy
yy
yz
xz yz zz
式(1)是关于 n 的三次方程,它的三个根,即为3个主应力。
在给定的应力状态下,由于物体内任一点的主应力不会随着坐标 系的改变而改变,所以 I1、I2、I3 也不会随着坐标系而改变,称 I1、I2、I3 分别为第一,第二,第三应力张量不变量。当坐标轴与主应力方 向重合,则应力不变量可以写成下式:
谢谢观赏
2、材料塑性行为的屈服准则 该准则用来确定材料产生的屈服时的临界应力状态。大量实
验表明:多数材料的塑性屈服与静水压力无关,对于复杂应力状 态,由等倾面组成的八面体上的正应力恰好就是静水压力。 这里说明下何为静水压力:
对于空间中的一个斜面的力进行分解
如左图,微小体元的一个斜面ABC, 其外法线上的为n,设该斜面上的剪 应力为零,则该斜面上只有正应力, 其方向沿着n 法线n,沿着坐标轴三个 方向的分解为
px nx n , py ny n , pz nz n
其中 nx , ny , nz 为斜面外法线n的方向 余弦
△ABC △S △BOC nx△S △COA ny△S △AOB nz△S
由 Fx 0
px△S xxnx△S yxny △S zxnz △S Fx△V 0
当OABC P :
行业PPT模板:www.1p pt.co m/ hang ye / PPT素材下载:/sucai/ PPT图表下载:www.1p pt.co m/ tubiao/ PPT教程: /powerpoint/ Excel教程:www.1ppt.c om/excel/ PPT课件下载:www.1p pt.co m/ kejian/ 试卷下载:www.1ppt.c om/shiti /
相关文档
最新文档