6.1 塑性理论及有限元
6.1 塑性理论及有限元
![6.1 塑性理论及有限元](https://img.taocdn.com/s3/m/635f568ab9d528ea81c77926.png)
b B s p A’ A
C
O
E
p
f
F
e
由于加载、卸载规 律 不 同 , 导 致 -ε 关 系不唯一。只有知道 变形历史,才能得到 一 一 对 应 的 -ε 关 系 ,即塑性变形与变形 历史或路径有关。这 是第3个重要特征。
b B s p A’ A
C
O
E
p
f
F
e
> s 以后的点都可
铝合金的标称应力-应变曲线
真实应力-应变曲线(应使用真应力和真应变来讨论塑性变形状态)
真实应力:真实应力为轴力除以真实的横截面面积,对于均匀变形:
P A
式中
P A0
加载瞬间载荷; 同一瞬间试样横断面积。
真实应变:真实应变为长度改变量除以当时长度:
n ln dl l1 l0 l2 l1 l3 l2 ln ln1 li ln ln l0 l l0 l1 l2 ln1 l0 i 1 li
塑性理论及有限元
1 2 3 4 单向弹塑性应力应变关系 塑性条件方程 塑形本构关系 全量和增量理论
塑性力学的主要任务是研究固体发生塑性变形时的应力分布和应变分布的规 律。由于大量的工程设计问题(如锻造、压延、冲压、拉丝等)及许多自然 现象(如结构的弹塑性断裂与损伤、地壳的构造运动等)都与塑性变形有关 ,最近几十年来,塑性力学得到迅速的发展。 从物理上看,塑性变形过程属于不可逆过程,伴随着机械能的耗散,相比弹 性变形过程更复杂。从数学上看,塑性力学本质上是非线性的,所以比弹性 力学问题困难得多。塑性力学问题的研究不是探究材料塑性变形的内在机理 ,而是从材料的宏观塑性行为中抽象出力学模型,并建立相应的数学方程来 予以描述。 应力平衡方程和应变位移间的几何关系是与材料性质无关的,因此对弹性力 学与塑性力学都一样;两者的主要差别主要变现在应力与应变的物理关系, 即本构关系上。这也带来了处理方法上的不同,塑性力学首先要解决的问题 是在实验资料的基础上确立塑性本构关系,进而与平衡和几何关系一起去建 立弹塑性边值问题;然后根据不同的具体情况寻求数学计算方法求解弹塑性 边值问题。 塑性力学的主要研究对象是金属材料,人们对此进行了比较充分的实验研究 ,并总结出若干塑性理论。而对于非金属材料相对来说研究得不够充分,因 此将塑性理论应用到非金属材料和其它材料上时,需持非常审慎的态度。塑 性力学今后的研究重点是用宏观与微观相结合的方法研究材料的塑性变形机 理,吸收材料学的近代成果,丰富和发展原有的理论,推动研究向前发展。
弹塑性有限元法基本理论与模拟方法
![弹塑性有限元法基本理论与模拟方法](https://img.taocdn.com/s3/m/9b846b64e3bd960590c69ec3d5bbfd0a7956d5e3.png)
用于模拟流体流动和传热问题 ,如流体机械、航空航天和化 工等领域。
电磁场
用于分析电磁场问题和电气设 备性能,如电机、变压器和天 线等。
声学
用于模拟声音传播和噪声控制 问题,如声学器件和声学环境
等。
04 弹塑性有限元法的基本原 理
弹塑性有限元法的离散化方法
有限元离散化
将连续的物理场或结构体离散为有限个小的单元体, 每个单元体之间通过节点相互连接。
结构强度分析的模拟
结构强度评估
通过弹塑性有限元法模拟,可以对结构的强度进行评估,预测结构在不同载荷下的响应, 确保结构的安全性和稳定性。
疲劳寿命预测
利用弹塑性有限元法,可以模拟结构的疲劳载荷历程,预测结构的疲劳寿命,为结构的维 护和更换提供依据。
结构优化设计
通过模拟结构的应力分布和变形,可以优化结构设计,降低结构重量,提高结构效率。
边界条件和初始条件
在平衡方程中考虑边界条件和初始条件,以确保模拟的准确性和收 敛性。
弹塑性有限元法的边界条件和初始条件
边界条件的处理
01
根据实际情况,将边界条件转化为节点约束或单元载荷的形式。
初始条件的设置
02
在非稳态问题中,需要考虑初始条件的设置,以模拟问题的初
始状态。
边界条件和初始条件的实施
03
随着计算机技术的不断发展,弹塑性 有限元法在各个工程领域中得到了广 泛应用,如机械、航空航械设计中,弹塑性有限元法可用于分析各种复杂结构 的应力分布、变形和疲劳寿命等,提高产品的可靠性和安 全性。
航空航天
在航空航天领域,弹塑性有限元法可用于分析飞行器结构 在各种载荷下的响应,优化结构设计,提高飞行器的性能 和安全性。
NFEM06-粘塑性有限元
![NFEM06-粘塑性有限元](https://img.taocdn.com/s3/m/450932fd04a1b0717fd5ddec.png)
Chap 6 弹-粘塑有限元§6.1 粘性介质的本构关系与机械元件模型弹性变形和塑性变形都与时间无关,认为是瞬时发生的。
但实际上变形的发生与时间有关。
只是一般情况下,时间的影响可忽略不计,然而对某些情况,例如高温下的金属变形,外力作用下的岩体变形,可以随时间累积至比较大的量值。
因而对它们进行变形和应力分析,必须考虑时间因素。
介质这种变形和应力随时间变化的特性称为粘性。
(6-8)σ作用,并假定粘塑性具(6-15)vp=0,基本方程式成为ε(6-17)采用五元件粘弹性模型来对围压作用下片岩的粘弹性流变曲线进行辨识.工作应力水平,岩石均表现为典型的粘弹性特性。
由于对出现加速流变阶段的流变曲线进行辨识。
式中ε为岩石总的应变量;t 为经过的时间;E 1为瞬时弹性模量;E 2和E 3均为粘弹性模量;η1和η2均为粘滞系数,表示流变阶段趋向稳定的快慢程度。
而σ1和σ3分别为轴向应力和围压值。
五元件粘弹性模型参数求算的方法采用直接迭代法,需要求取5个流变参数:E 1,E 2,E 3,η1和η2。
根据试验得到n 组数据(ε,t ),首先给定一组初始近似值(E 10,E 20,E 30,η10,η20),以此为基础进行迭代,这样反复计算,直到满足所需要的精度。
(6-18)四. 粘弹性三维本构方程弹性介质的应力和应变服从广义虎克定律,令[]Tzxyz xy z y x e e e e e e 222=e 可有eD s 2G =式中⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=110122022D (6-19)(6-20)(6-21)材料的本构方程为εD I 1K m =′σ其中,⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=0011101111111D 牛顿粘性介质的应力和应变应服从下式eD s &2η=对于Kelvin模型的粘弹性介质,由于介质应力是弹性应力与粘性应力之和eD e D s &22η+=G (6-22)(6-23)(6-24)(6-25)ANSYS粘塑性ANSYS中可用的三种其它粘塑性本构模型。
弹塑性力学与有限元
![弹塑性力学与有限元](https://img.taocdn.com/s3/m/8056b700b52acfc789ebc994.png)
§1.3 发展与研究方法7
钱学森
钱伟长 胡海昌
§1.3 发展与研究方法8
杨桂通
徐芝伦
§1.3 发展与研究方法9
•弹性力学——促进数学和自然科学基本理论 的建立和发展; •广泛工程应用——造船、建筑、航空和机械 制造等。 •发展——形成了一些专门的分学科; •现代科学技术和工程技术——仍然提出新的 理论和工程问题。 •对于现代工程技术和科研工作者的培养—— 对于专业基础,思维方法以及独立工作能力都 有不可替代的作用。
§1.3
弹性力学的发展 和研究方法
弹性力学是一门有悠久历史的学科,早期 研究可以追溯到1678年,胡克 (R.Hooke)发现胡克定律。 这一时期的研究工作主要是通过实验方法 探索物体的受力与变形之间的关系。
§1.3 发展与研究方法2
•近代弹性力学的研究是 从19世纪开始的。 •柯西1828年提出应力、 应变概念,建立了平衡微 分方程,几何方程和广义 胡克定律。 •柯西的工作是近代弹性 力学的一个起点,使得弹 性力学成为一门独立的固 体力学分支学科。
•——物理关系或者物理方程 •线性弹性体和非线性弹性体
§1.1 弹性力学任务11
研究方法的差别造成弹性力学与材料力 学问题的最大不同。
•常微分方程,数学求解没有困难。 •偏微分方程边值问题,在数学上求解困难重 重,除了少数特殊问题,一般弹性体问题很难 得到解析解。 •这里并不是说弹性力学分析不再需要假设,事 实上对于任何学科,如果不对研究对象作必要 的抽象和简化,研究工作都是寸步难行的。
三维数学问题,综合分析的结果是偏微分 方程边值问题。
§1.1 弹性力学任务4
建筑工程
§1.1 弹性力学任务5
建筑工程
§1.1 弹性力学任务6
塑性材料的有限元分析
![塑性材料的有限元分析](https://img.taocdn.com/s3/m/c5ccb3f01b37f111f18583d049649b6648d709bb.png)
针对复杂材料和结构,需要深入研究材料的非线 性行为和多场耦合效应,建立更加完善的物理模 型和数值算法。
此外,应加强与实验研究的结合,通过实验验证 和修正有限元模型,提高模拟结果的可靠性。同 时,实验研究也能够为有限元分析提供更加真实 和全面的材料性能数据。
THANK YOU
03
有限元分析方法
有限元分析的基本原理
离散化
将连续的物理系统离散为有限个小的单元,每个 单元称为有限元。
近似解法
通过数学方法求解每个有限元的近似解,再通过 组合所有有限元的解得到整个系统的近似解。
平衡方程
建立每个有限元的平衡方程,通过求解平衡方程 得到每个节点的位移和应力。
有限元分析的实现过程
然而,塑性材料的有限元分析仍存在 一些挑战和限制,如模型的简化、边 界条件的确定、材料参数的获取等, 需要进一步研究和改进。
研究展望
未来研究应致力于发展更加精确和高效的有限元 分析方法,提高模拟结果的可靠性和精度。
在实际工程应用中,应加强有限元分析与其他数 值方法(如边界元、有限体积等)的结合,实现 优势互补,提高计算效率。
塑性变形的微观机
制
塑性变形是通过位错的滑移和攀 移等微观机制实现的,这些机制 在宏观上表现为塑性变形。
塑性变形的温度效
应
温度对塑性变形的影响较大,温 度升高会使材料的屈服强度降低, 塑性变形能力增强。
塑性变形的加工硬
化
在塑性变形过程中,材料的屈服 强度会随着变形程度的增加而逐 渐提高,这种现象称为加工硬化。
背景
随着计算机技术的不断发展,有限元分析已成为工程领域中解决复杂问题的常 用方法。通过有限元分析,可以模拟材料的变形、应力分布、应变等,为实际 工程提供重要的理论依据。
弹塑性力学与有限元绪论
![弹塑性力学与有限元绪论](https://img.taocdn.com/s3/m/81b7ade4a26925c52cc5bfa7.png)
《弹塑性力学与有限元》
矢量和张量
矢量代数、标量积、矢量积、三重积、标量场和矢量场
如下图,三维空间直角坐标系Oxyz,
x3
x1, x2 , x3
x, y, z
P
rV
e3
o
x1
e1
e2
n i 1
ai
xi
• 另外 ( x y z )( x y z ) 应写成 ii jj ,不能写作 ii ii ,因为 后者的标号重复了4次。
《弹塑性力学与有限元》
矢量和张量
指标记法与求和约定
求和约定
定义:凡在同一项内不重复出现的指标,如 a ji xi bj , j 为自由指标
《弹塑性力学与有限元》
弹塑性力学与有限元
建立力学模型
模型建立以后,对其采用适当的求解方法(解析 解和有限元)进行计算并对计算的结果进行分析整理 ,返回实际问题进行验证。
一般通过实验验证:直接实验验证:直接实验比 较简单时可以直接进行,但有时十分困难。相似模型 实验:相似实验的模型一般应与实际问题的边界条件 和形态是几何相似的。
《弹塑性力学与有限元》
弹塑性力学与有限元
教学内容:包括绪论和如下三部分内容:
(一)弹性力学理论
(三)有限单元法
1) 矢量与张量 2) 应力分析 3) 应变分析 4) 弹性应力-应变关系
(二)塑性力学理论
1) 有限元法的理论基础——加权余量法和变分原理 2) 弹性力学问题有限元方法的一般原理和表达格式 3) 单元和插值函数的构造 4) 等参元和数值积分 5) 有限元法应用中的若干实际考虑 6) 线性代数方程组的解法
a a21 a22 a23 a11a22a33 a12a23a31 a13a21a32 a31 a32 a33 a13a22a31 a12a21a33 a11a23a32
弹塑性力学与有限元-弹塑性应力-应变关系
![弹塑性力学与有限元-弹塑性应力-应变关系](https://img.taocdn.com/s3/m/0074a5eca417866fb84a8ef6.png)
《弹塑性力学与有限元》
弹塑性应力-应变关系
强化法则
Mises初始屈服条件
J2
2 s
3
0
3J 2 s 0
加载(后继屈服)条件
3J 2 0
3
2
sij
sij
0
( d p )0
函数可通过单轴拉伸下实验曲线确定.
《弹塑性力学与有限元》
弹塑性应力-应变关系
强化法则
➢ 随动强化
• 几何特点(在应力空间):形状和大小、方向保持不变,只 是中心位置发生改变,加载面作刚体移动。
量,硬化参量记为 .
《弹塑性力学与有限元》
弹塑性应力-应变关系
强化法则
目前常用的硬化参量有如下几种:
1.塑性功 w p, w p
ij
d
p ij
是目前岩土弹塑性理论中用得较多的。
2.有效塑性应变
p ij
3.等效塑性剪应变 p S
2 3
d
ijpd
p ij
4.塑性体应变
p v
p x
p y
p z
弹塑性应力-应变关系
强化材料的增量应力-应变关系
其中, 是弹塑性柔度张量,表示为 • 以应变增量表示的应力增量
考虑到式(7.29)和式(7.105)有
把dλ和应变增量联系起来,则有 其中 从式(7.4),式(7.179)和式(7.112)可得
《弹塑性力学与有限元》
弹塑性应力-应变关系
强化材料的增量应力-应变关系
• 物理意义:材料在强化后为各向异性。
•
数学表示:f (ij,ij)=f 0(ij-ij) k = 0
f 0(ij-ij) = k
ij 是一个表征加载面中心移动的应力值,称为反(背)应力
弹塑性力学与有限元:2 力学位移和应变分析T
![弹塑性力学与有限元:2 力学位移和应变分析T](https://img.taocdn.com/s3/m/1fe16c339e31433238689353.png)
O
x u u dx
x
u
u x
dx
u
dx
u x
PB的正应变:
u
P
dx
v P A
dy
x v v dx
x
y
v
v y
dy
dy
v
v y
P点的剪应变:
y
v v dy y
B
A
B
u u dy
P点两直角线v 段夹v角d的x 变 v化
tan
x dx u dx
xy
v x
u x
u
u
x dy
u
tan
y dy v dy
符号规定:u,v,w与坐标轴正方向一致为正,相反为负。
考虑外力作用下的两种状态: 平衡状态:M点只随位置变化,不随时间变化;位移分量(u,v,w)只随位置变化, 不随时间变化。 运动状态: M点不仅随位置变化,而且随时间变化;位移分量(u,v,w)随位置和 时间变化而变化。
本章仅考虑平衡状态。
根据连续性假设,物体上任一点M,当物体变形后, 都一一对应于相应的点M’;
考察P点邻域内线段的变形:
PA dx dy
y
u
P
dx
x u u dx x
v P A
dy
v v dx x
B
A
变形前 P
A
变形后
P
u v
v v dy y
B
u u dy y
u u dx x
A v v dx B
x
u u dy
B
y v v dy
y
注:这里略去了二阶以上高阶无穷小量。
PA的正应变:
基于广义塑性理论上界法的有限元法及其应用
![基于广义塑性理论上界法的有限元法及其应用](https://img.taocdn.com/s3/m/62dccc3d10661ed9ad51f3fb.png)
以原点为圆心的圆域.见图 2 。由于三角形单元中
采用线 性速 度场 形式 ,所 以,此 塑 性势 函数也 需要 进 行 线 性化 .为此 ,可 参 照文 [,5的方 法 ,用 ~ 2 ] 个外切 正多边 形来逼近 上述 圆域 。设正 多边形 的边 数为 m.则第 k 塑性 势 函数 的线 性 表达 式 为 边 Q =Aa +Bdv t”一 c +C f 3 t 7
式 中 : A :C S k B =一 O , C :2 i c O ̄ , k CS k n q, s
=
2nm,k 1 2….m k/ =,.
£T
图 l 上界法的三 角形单元
F g 1 T in u a lme t f p e o n t o i ra g l e e n u p r u dme h d r o b
的与机动容许的塑性变形速度场相对应 的荷载中, 极 限荷载为晟小。采用正交流动法 则时,有
f d+ “≤ , f V ( l L I: + , 2 d ed ) AL
采 用 非正交 流动 法 则 时 .有 I
I vV L棚 ≤f;d+ .; + Ed  ̄ V 4
2 基于广义塑性理论上界法 的原理
时不会产生塑性体变 ; , 分别为作用在物体 上的面 力和体力 ; . 分别 为速度场 中 的速 度 和 应变率;o 为静力场中的应力 ;e - , j , 为土体抗剪 强度 指标 ;△ 为速度 间断线两侧 切 向速度 变 化最 ; r , 分别为速度间断线 £上的剪应 力和正应力。
1 引
言
剪胀变 形 ,当采用 非正交 流动法 则时 , = ,此
2
2 纪5 0世 0年代 , ree和 Pae 把应力场 和 D ukr rgr 速度场 结合起 来 并提 出极值 理论 ,建立 了极限分析 理论 。文 【】 于关 联流 动法 则 提 出 了岩 土 的上 界 1基 法 ,此 后 。Son 又给 出了相应 该上 界法 的有限元 la
弹塑性力学基础与有限元分析-接触分析实例
![弹塑性力学基础与有限元分析-接触分析实例](https://img.taocdn.com/s3/m/ef8ffec1a1116c175f0e7cd184254b35eefd1af0.png)
06
结论与展望
结论
1
本文通过理论分析和有限元模拟,深入研究了弹 塑性力学基础与有限元分析在接触分析中的应用。
2
研究结果表明,弹塑性力学基础与有限元分析在 接触分析中具有较高的精度和可靠性,能够有效 地模拟复杂接触问题。
3
本文所采用的有限元分析方法在处理接触问题时 具有较好的通用性和扩展性,为进一步研究复杂 接触问题提供了有力支持。
弹塑性本构模型
弹塑性本构模型的定义
弹塑性本构模型是描述弹塑性材料力学行为的数学模型,它通过应力应变关系来描述材料的弹塑性行 为。
常见的弹塑性本构模型
常见的弹塑性本构模型包括Mohr-Coulomb模型、Drucker-Prager模型、Cam-Clay模型等。这些模 型在描述材料的弹塑性行为方面各有特点,适用于不同的材料和工程问题。
接触面完全贴合,无相对运动。
滑动状态
接触面部分贴合,存在相对运动。
混合状态
接触面同时存在分离、粘结和滑动。
接触检测与跟踪
初始接触检测
确定初始状态下接触面的位置和状态。
接触状态跟踪
实时监测接触面的运动状态和相互作用。
接触面更新
根据接触状态调整接触面的几何形状和参数。
接触刚度与阻尼
1 2
接触刚度
描述接触面间的相互作用力与相对位移的关系。
求解阶段主要进行有限元 方程的求解,得到各节点 的位移和应力等结果。
ABCD
前处理阶段主要完成有限元 模型的建立和网格划分,为 求解阶段提供输入数据。
后处理阶段主要对求解结果进 行可视化、分析和评估,为工 程设计和优化提供依据。
04
接触分析原理
接触状态描述
分离状态
塑性成形过程中的有限元法
![塑性成形过程中的有限元法](https://img.taocdn.com/s3/m/6bdda02eb4daa58da1114a05.png)
塑性成形过程中的有限元法金属塑性成形技术是现代化制造业中金属加工的重要方法之一。
它是金属材料在模具和锻压设备作用下发生变形,获得所需要求的形状、尺寸和性能的制件的加工过程。
金属成形件在汽车、飞机仪表、机械设备等产品的零部件中占有相当大的比例。
由于其具有生产效率高,生产费用低的特点,适合于大批量生产,是现代高速发展的制造业的重要成形工艺。
据统计,在发达国家中,金属塑性成形件的产值在国民经济中的比重居行业之首,在我国也占有相当大的比例。
随着现代制造业的高速发展,对塑性成形工艺分析和模具设计方面提出了更高的要求。
若工艺分析不完善、模具设计不合理或材料选择不当,则会造成产品达不到质量要求,造成大量的次品和废品,增加了模具的设计制造时间和费用。
为了防止缺陷的产生,以提高产品质量,降低产品成本,国内外许多大公司企业及大专院校和研究机构对塑性成形件的性能、成形过程中的应力应变分布及变化规律进行了大量的理论分析、实验研究与数值计算,力图发现各种制件、产品成形工艺所遵循的共同规律以及力学失效所反映的共同特征。
由于塑性成形工艺影响因素甚多,有些因素如摩擦与润滑、变形过程中材料的本构关系等机理尚未被人们完全认识和掌握,因而到目前为止还未能对各种材料各种形状的制件成形过程作出准确的定量判定。
正因为大变形机理非常复杂,使得塑性成形研究领域一直成为一个充满挑战和机遇的领域。
一般来说,产品研究与开发的目标之一就是确定生产高质量产品的优化准则,而不同的产品要求不同的优化准则,建立适当的优化准则需要对产品制造过程的全面了解。
如果不掌握诸如摩擦条件、材料性能、工件几何形状、成形力等工艺参数对成形过程的影响,就不可能正确地设计模具和选择加工设备,更无法预测和防止缺陷的生成。
在传统工艺分析和模具设计中,主要还是依靠工程类比和设计经验,经过反复试模修模,调整工艺参数以期望消除成形过程中的产品缺陷如失稳起皱、充填不满、局部破裂等。
仅仅依靠类比和传统的经验工艺分析和模具设计方法已无法满足高速发展的现代金属加工工业的要求。
第六章 塑性力学基本概念
![第六章 塑性力学基本概念](https://img.taocdn.com/s3/m/929fb3355a8102d276a22ffa.png)
第六章 塑形力学的基本概念
6.1绪论 绪论 什么是塑性力学? 什么是塑性力学? •塑性力学是相对于弹性力学而言。 塑性力学是相对于弹性力学而言。 塑性力学是相对于弹性力学而言 •在弹性力学中,物质微元的应力和应变之间具有单一 在弹性力学中, 在弹性力学中 的对应关系。然而, 的对应关系。然而,材料在一定的外界环境和加载条 件下,其变形往往会具有非弹性性质, 件下,其变形往往会具有非弹性性质,即应力和应变 之间不具有单一的对应关系。 之间不具有单一的对应关系。 非弹性变形包括塑性变形和粘性变形: 非弹性变形包括塑性变形和粘性变形:
塑性变形- 塑性变形-指物体在除去外力后所残留下来的永 久变形,习惯上按破坏时的变形大小分为塑性和脆性, 久变形,习惯上按破坏时的变形大小分为塑性和脆性, 如果材料的延性好,进入延性仍能承受荷载。 如果材料的延性好,进入延性仍能承受荷载。 塑性力学来研究这类问题。 塑性力学来研究这类问题。
粘性变形随时间而改变,例如蠕变、应力松弛等, 粘性变形随时间而改变,例如蠕变、应力松弛等, 这里不研究。 这里不研究。
σ ≤σs,
ε=
σ
E
σ > σs,
ε=
σs
E
+
σ −σ s
Et
O
A
E
刚性线性强化模型 问题:卸载 线怎样描述 线怎样描述? 问题:卸载BE线怎样描述?
εp
ε
εe
3、幂指数硬化模型: 、幂指数硬化模型:
将硬化阶段的曲线简化为一条幂指数曲线,
σ ≤ σs,
σ > σs,
ε=
σ
E
σ σb B σs σp A’ A C
6.2 材料实验结果
一、单轴拉伸实验 • 材料塑形变形性质通过试验研究获得。 材料塑形变形性质通过试验研究获得。 • 最简单实验是室温单轴拉压实验: 最简单实验是室温单轴拉压实验: •材料:金属多晶体材料 材料: 材料 •试件如图 试件如图 •名义应力和名义应变定义为 名义应力和名义应变定义为
有限元_塑性力学
![有限元_塑性力学](https://img.taocdn.com/s3/m/239c76e95022aaea998f0f64.png)
rg
连续体
三结点三角形单元
六结点三角形单元
矩形单元 任意四边形单元 8结点曲边 四边形单元
2. 单元分析
建立: {F}e=[k]{d}e
[k]:单元刚度矩阵
y vj j
vi , (Vi) i ui , (Ui)
uj
vm
m um
x
单元结点力向量{F}e=[Ui Vi Uj Vj Um Vm]T
利用单元的结点位移将整个单元的位移分量表
示为坐标的函数。
y
i
2. 三结点三角形单元的位移模式 j P
设:u=a1+a2 x+a3 y v=a4+a5 x+a6 y
m x
系数a1~a6由结点位移 ui , vi , uj , vj , um , vm确定。
• 将位移模式写成结点位移的显式:
i
u= Niui+ Njuj +Nmum
单元结点位移向量{d}e=[ui vi uj vj um vm]T
•
体力、面力
静力等效 ——
等效结点荷载
3. 整体分析
建立: {F}=[K]{d}, [K]:整体刚度矩阵
由各结点平衡{F}={R},得有限元方程:
[K]{d}= {R}
0.3 位移模式与解答的收敛性
1. 什么是位移模式(位移函数)
下标1,2,3表示坐标 x1, x2, x3(即x, y, z)方向.
2. 一点斜面上的应力
3
SNi i1l1 i2l2 i3l3 ijl j j1 (i 1, 2, 3)
x3 N
SN3
SN
SN1 O
弹塑性力学与有限元-屈服总则和弹塑性应力-应变关系
![弹塑性力学与有限元-屈服总则和弹塑性应力-应变关系](https://img.taocdn.com/s3/m/85e6925e5fbfc77da269b1a7.png)
《弹塑性力学与有限元》
屈服总则和弹塑性应力-应变关系
与静水压力无关材料的屈服总则
s
2
➢ Von Mises屈服条件
2 s1 3s1
2 s 2 3s 2
2 s3 3s3
0
s1
s3
x s1 s2 sin 30o s3 sin 30o
x
1 6
2s
1
s2
s3
y s2 cos 30o s3 cos 30o
屈服总则和弹塑性应力-应变关系
与静水压力无关材料的屈服总则
对于金属类材料,由于材料的屈服对静水压力不敏感,因而剪切应 力控制其屈服。
➢ Tresca屈服条件 1864年,Tresca作了一系列的挤压实验来研究屈服条件,金属材料在屈 服时,可以看到接近于最大剪应力方向的细痕纹(滑移线),因此塑性 变形可以是由于剪切应力所引起的晶体网格的滑移而引起的。
认为最大剪应力达到极限值时开始屈服: (s1 s 2 s3 )
max (s1 s 3 ) / 2 k
(6.18)
(材料力学的第三强度理论)
《弹塑性力学与有限元》
屈服总则和弹塑性应力-应变关系
如不规定s1 s 2 s3
xs
1 2
(s 1
s
3
)
ys
1 6
(2s
2
s1
s
3
)
s1 s2 s3
《弹塑性力学与有限元》
屈服总则和弹塑性应力-应变关系
屈服总则定义
➢ 应力空间和主应力空间 L直线:主应力空间中过原点并与坐标轴成等角的直线。
其方程为 s1 s2 s3 显然,
L直线上的点代表物体中承受静 水应力的点的状态,这样的应力 状态将不产生塑性变形。
弹塑性问题有限元分析概述
![弹塑性问题有限元分析概述](https://img.taocdn.com/s3/m/58f023f5172ded630b1cb653.png)
行业PPT模板:/hangye/ PPT素材下载:/sucai/ PPT图表下载:/tubiao/ PPT教程: www.1ppt.c om/powerpoint/ Excel 教程:www.1ppt.c om/excel/ PPT课件下载:/kejian/ 试卷下载:www.1ppt.c om/shiti /
2 2 n p 2 n 2 12 n x 2 2 n y 3 2 n z n (3 )
由于
n p1nx p2 n y p3 nz 1n 2 x 2 n 2 y 3 n 2 ( ) z 4
由( 3 )、( 4)可得 :
PPT模板下载:/moban/ 节日PPT模板:/jieri/ PPT背景图片:/beijing/ 优秀PPT下载:/xiazai/ Word教程: /word/ 资料下载:/ziliao/ 范文下载:/fanwen/ 教案下载:/jiaoan/
若要得到( 5) , 根据(4)可得: n x ny n z 1 / 3
该八面体上的切应力为 1 8 ( xx yy ) 2 ( yy zz ) 2 ( zz xx ) 2 6( 2 xy 2 yz 2 xz ) 3 它是决定材料是否屈服的力学参量,因此初始屈服条件为
p x nx n , p y n y n , p z nz n 代入上式可得: nx ( xx n ) n y xy nz xz 0 nx yz n y ( yy n ) nz yz 0 nx zx n y zy nz ( zz n ) 0 这是关于nx , n y , nz的齐次线性方程组,其 非零解的条件为行列式 等于零 展开可得:
弹塑性力学与有限元-等参元数值分析
![弹塑性力学与有限元-等参元数值分析](https://img.taocdn.com/s3/m/d72c728e866fb84ae45c8df6.png)
弹塑性力学与有限元 —等参元和数值积分
等参元和数值积分
等参元和数值积分
等参元
➢ 等参单元的基本概念和单元矩阵的变换 ➢ 等参变换的条件和等参元的收敛性 ➢ 等参元用于分析弹性力学问题的一般格式
数值积分
➢ 数值积分方法 ➢ 等参元计算中数值积分阶次的选择
《弹塑性力学与有限元》
x a 2 b c
y d 2 e f
(4-13)
可见在整体坐标系中,单元的边是一条抛物线或退化为一条直线。
《弹塑性力学与有限元》
等参元和数值积分
等参单元的矩阵变换
图4-6 ANSYS提供的Plane82单元
《弹塑性力学与有限元》
等参元和数值积分
等参单元的矩阵变换
如图4-6所示,ANSYS提供的PLANE82单元是一个四边形八结 点等参单元,局部坐标定义为s和t,如图4-7所示。PLANE82单元可 以退化为三角形六结点单元。
《弹塑性力学与有限元》
等参元和数值积分
等参单元的基本概念和单元矩阵的变换
形函数为,
1x y
Ni
(1 4
)(1 a
) b
Nj
1 (1 4
x )(1 y ) ab
Nm
1 (1 4
x )(1 a
y) b
Np
1 (1 4
x )(1 a
y) b
上述单元位移模式满足位移模式选择的基本要求:反映了单元的刚体位
因此给出任意四边形单元的结点位移就能得到整个单元上的位移,上 式给出的位移模式就是所要找的正确的位移模式。把局部坐标与整体坐 标的变换式也取为:
x N1x1 N2 x2 N3x3 N4 x4 (4-8) y N1 y1 N2 y2 N3 y3 N4 y4
弹塑性力学及有限元法_
![弹塑性力学及有限元法_](https://img.taocdn.com/s3/m/cde2b165ddccda38376bafc2.png)
写成矩阵形式
R11 cos 2 θ x 1 Ry1 EA cos θ sin θ 1 = Rx 2 l1 − cos 2 θ R1 2 − cos θ sin θ y cos θ sin θ sin 2 θ − cos θ sin θ − sin 2 θ − cos 2 θ − cos θ sin θ cos 2 θ cos θ sin θ
单元刚度矩阵的子矩阵 K ij 表示:当单元 e 中节点 j 取单 位位移,且其它节点位移为零时,对应于 i 节点的节点力。
第五章 有限元法简介
单元1的节点力和节点位移的关系可写成
R1 K11 = R2 K 21
1
K12 K 22
1
δ1 δ 2
1 θFx1(u1) 3 Fx3 (u3) Fy1(v1 ) Fy3 (v3) y 2 o x
1
Fy2 (v2) Fx2(u2)
2
图5-1 简例结构图
第五章
分析步骤:
有限元法简介
2
1
1 1 Ry2(v2) 1 1 Rx2(u2)
1. 离散结构物为有限个单元 分为2个单元,第一个单元的节点编号 为1和2,第二个单元的节点编号为2和3。 对于第一单元,在第1、2节点处的节点力 为 R 11 , R 11 , R 1 2 , R 1 2 ,表示节点施加在单元1上 x y x y
1 − cos θ sin θ u1 1 2 − sin θ v1 cos θ sin θ u1 2 1 si成
R11 k x 1 11 Ry1 k21 1 = Rx 2 k31 R1 k41 y2 k12 k22 k32 k42 k13 k23 k33 k43
弹塑性力学与有限元
![弹塑性力学与有限元](https://img.taocdn.com/s3/m/135b0968767f5acfa1c7cd76.png)
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)。
失稳破裂 弹性
金属材料单轴加载时的应力与应变特征:
b B s p A’ C
A
E
p
f F
e
(1)加载开始后,当 应力小于A点的应力 值时,应力与应变呈 线性关系。材料处于 线弹性变形阶段。A 点的应力称为比列极 限。在此阶段卸载, 变形沿OA线返回。
b B s p A’ A
b
B
C
s A’ p A
O
E
p
f
F
e
静水压力试验
Bridgman等人在不同的静水压力容器中做单向拉伸试验。结 果表明: 1. 关于体积变形。 Bridgman的研究表明,直到1500MPa,体积 变形仍然是弹性的。体应变 与静水压力 p 之间的关系为
ap bp2
金属材料的a、b值可参阅相关书籍。所以在塑性理论中常认为体 积变形是弹性的。试验还表明,这种弹性体积变化是很小的,对 于金属材料,当发生较大塑性变形时,可忽略弹性体积的变化, 即认为在塑性变形阶段,材料是不可压缩的。 2.材料的塑性变形与静水压力无关。从对比试验可以发现,静 水压力对初始屈服应力的影响很小,因而对钢材等金属材料, 可以认为塑性变形不受静水压力的影响。 注意:对铸铁、岩石等脆性材料,土壤多孔材料,静水压 力对屈服应力和塑形变形有明显影响。
b B s p A’ A
C
O
E
p
f
F
e
在简单压缩下,忽略摩擦影响, 得到的压缩 -ε与拉伸-ε基本相同。 但是若将拉伸屈服后的试样经卸载 并反向加载至屈服,反向屈服一般 低于初始屈服。同理,先压后拉也 有类似现象。这种正向变形强化导 致后继反向变形软化的现象称作 Bauschinger 效应。这是金属微观组 织变化所致。一般塑性理论分析不 考虑Bauschinger效应。 由于塑性应变不可恢复 ,所以 外力所做的塑性功具有不可逆性, 或称为耗散性,在一个加载卸载循 环中外力做功恒大于零,这一部分 能量被材料的塑性变形耗散掉了。
、ln 是每步载荷增量后的构件长度,因为上式是 其中, l1、l2、 用对数表示的,真实应变又称为对数应变,也称为自然应变。
利用对数应变可以对应变使用加法,如 1 ln(l1 是两个相继发生的应变,则总应变为
l0 ) 和2 ln(l2 l1 )
l2 l2 l1 ln ln ln 1 2 l0 l1 l0
b B s p A’ A
C
O
E
p
f
F
e
由于加载、卸载规 律 不 同 , 导 致 -ε 关 系不唯一。只有知道 变形历史,才能得到 一 一 对 应 的 -ε 关 系 ,即塑性变形与变形 历史或路径有关。这 是第3个重要特征。
b B f
F
e
> s 以后的点都可
e
p
对于 卸载
d 0 d 0
d Et d d Ed
从B点卸载到E点后,再重新加拉应 力(称为正向加载),这时应力应 变按卸载曲线BE变化。 当应力达到卸载前的B点应力,材 料才最新进入屈服。与B点对应的 应力,叫后续屈服应力。由于强化 效应,后续屈服应力通常比初始屈 服应力高,其提高程度与塑性变形 的历史有关,或者与塑性阶段的加 载历史有关。因此,后续屈服应力 可以写为
K
PK AK
lK ЄK ln l0
真实应力-对数应变曲线的确定
均匀变形
存在颈缩
o
真实应力-对数应变曲线
Є
真实与标称应力-应变曲线的比较
均匀变形
存在颈缩
o
—— 真实应力-对数应变曲线 —— 标称应力-对数应变曲线
Є
弹塑性变形基本特征(6个)
均匀塑性变形
如图 所示,金 属变形分 为弹性、 均匀塑性 变形、失 稳破裂三 个阶段。
1 单向弹塑性应力应变关系
基本概念 弹性(elasticity) 卸载后变形可以恢复特性,可逆性; 塑性(plasticity) 物体产生永久变形的能力,不可逆性; 屈服(yielding) 开始产生塑性变形的临界状态; 损伤(damage) 材料内部缺陷产生及发展的过程; 断裂(fracture) 宏观裂纹导致结构破断的现象。
C
O
E
p
f
F
e
应力在A~A’之间, 应力与应变关系不 再为线性关系。变 形仍为弹性。 若卸载,变形仍按 照原来的应力-应变 关系曲线返回初始 状态。
b B
C
s p
A’ A
O
E
由于一般材料的比例极 限、弹性极限和屈服应 力相差不大,不加区别, 统称屈服应力表示为。
p
f F
s BЄ
n
o
Є
线性强化刚塑性应力应变关系模型
s E1
o
Є
理想刚塑性应力应变关系模型
s
o
Є
• 应力状态与应变状态的进一步研究
• 前面我们已经阐明了有关应力与应变的 基本知识,为了今后论证问题的方便, 需要进一步补充相关知识 • 正八面体上的应力
弹塑性变形基本特征 •弹塑性共存 •加载卸载过程不同的-ε关系 •塑性变形与变形历史或路径有关 •硬化或强化 •Bauschinger效应 •静水压力
应力—应变关系的简化模型
金属材料的塑性行为相当复杂,为了使求解塑性问题成为可能,需要对材料 的塑性行为进行一些假设,同时还需要对应力-应变曲线加以理想化,也就是 根据不同的问题对实际材料特性作不同的简化。
常用应力应变关系模型
e
幂指数硬化应力应变关系模型
BЄ
n
o
Є
式中,B为常数,n可取0—0.1之间的任意数,一般由实际 的应力—应变曲线拟合而定
线弹性幂指数硬化应力应变关系模型
EЄ
s
n
s BЄ
s
o
Є
刚塑性幂指数硬化应力应变关系模型
塑性力学基本假设
材料为均匀连续,且各向同性;且拉伸和压缩的真应力-真应变曲线一致, 即不考虑Bauschinger效应; 材料具有无限的韧性,即认为材料可以进行无限的变形而不出现断裂; 材料的塑性行为与时间、温度无关,认为变形为准静态; 卸载时材料服从弹性规律,重新加载后的屈服应力等于卸载前的应力; e p 任何状态下的总应变可以分解为弹性与塑性两部分,即 , 且材料的弹性性质不因塑性变形而改变。 塑性变形是在体积不变(不可压缩)的条件下进行的,静水压力只产生 体积的弹性变化,不产生塑性变形。
在塑性理论中研究物体产生的塑性变 形条件时,除了用到最大 切应力外,还用到正八面体上的切应 力。正八面体的面就是通 过空间一点而和三个主平面夹角相等 的平面。取主平面为坐标 面,满足上述条件的八个面构成如图 所示的正八面体。 • 对于第一象限内的正八面体,它的方向余弦为 1 l0 m0 n0 3
从B点卸载到E点后,如加压应力(称 为反向加载),应力应变沿EF曲线变 化,材料在F点屈服。通常F点对应的 屈服应力明显低于比B点对应的应力值 (称为包辛格效应)。表明材料的后 续屈服性质不仅与它所经历的塑性变 形有密切关系,还受到它所经历的塑 性变形的方向有关。
E
C B
O
p
f F
e
若在均匀塑性变形阶段 出现卸载现象,一部分变 形得以恢复,另一部分则 成为永久变形。卸载阶段 -ε 呈 线 性 关 系 。 这 说 明 了塑性变形时,弹性变形 依然存在。 弹塑性共存、加载卸载 过程不同的-ε关系是塑 性变形的两个基本特征。
e
应力超过A’后,材料从 弹性状态进入塑形状态 。
b B A’ A
C
随应力的增加,变形不断 s 增加。应变硬化。
在硬化阶段,切线斜率 (硬化率)不断减小, 直至峰值应力 b 。
p
O
E
p
f
F
e
在应变进入硬化阶段后,如减小应力 (如在B点),应力与应变将不会沿 路径BAO返回到O点,而是沿BE路径 回到零应力。 b e p 弹性变形 恢复,塑形变形 保留
s A’ 在加载和卸载过程中应力和应变服从不同的 规律,这是材料在塑性阶段的一个重要特点。 p A 简单应力状态下的加载准则可以写成
C
B
加载 卸载
d 0 d 0
O f
E
p
F
e
此式也适用于 0 的压缩情况。有了这一
准则,我们可以把简单的拉伸试件在塑性阶 段的应力—应变关系归纳为
铝合金的标称应力-应变曲线
真实应力-应变曲线(应使用真应力和真应变来讨论塑性变形状态)
真实应力:真实应力为轴力除以真实的横截面面积,对于均匀变形:
P A
式中
P A0
加载瞬间载荷; 同一瞬间试样横断面积。
真实应变:真实应变为长度改变量除以当时长度:
n ln dl l1 l0 l2 l1 l3 l2 ln ln1 li ln ln l0 l l0 l1 l2 ln1 l0 i 1 li
金属受力变形过程
弹性变形-塑性变形-断裂
Smith W F. Foundations of Materials Science and Engineering. McGRAW.HILL.3/E
P
e
P e P
e
材料的单轴拉伸实验曲线有如图所示两种形态。
conditional yield limit 条件屈服极限