- 1、本文档共25页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
13-2静力隐式-显式和动力显式有限元列式
弹塑性大变形有限元方法
胡平
§ 大变形弹塑性本构方程
最简单的材料大变形弹塑性本构方程是不考虑加载历史的形变理论的全量本构方程。这类本构方程基于加载历史是比例加载的基本假定。对于许多金属冲压成形问题,比例加载的基本假定是近似正确的,因此,基于这类理论的所谓one step inverse algorithm,由于其高效高速的计算效率,近年来被广泛应用于汽车概念设计阶段的冲压工艺性快速校核。关于全量理论的One-step inverse algorithm的理论知识将另行介绍(参见理论手册)。
然而,更精细和准确的成形问题数值分析应该是能够考虑加载历史的增量本构理论。为了反映与加载历史的相关性,大变形弹塑性问题需要采用速率型(增量型)的本构方程。
首先,需要在大变形条件下,按照符合客观性的要求,建立增量型本构方程。由(7.26)式可知,变形率张量[d]是客观张量,但由(7.36)式可知Cauchy应力张量的物质导数不是客观张量,所以若用变形率张量和Cauchy应力张量来建立本构方程,则将不满足本构方程的客观性条件(7.42)式。为此,需要另外定义Cauchy应力张量的导数。
1.Cauchy应力张量的Jaumann导数(率)
而Cauchy应力的Jaumann导数定义为
(1)
推导后得到
(2)
其矩阵形式为
2. 第一Piola应力(名义应力)张量的本构导数
(3)
用应力导数表示的第一应力的本构导数为
(4)
其矩阵形式为
(5)
3. 第二Piola应力张量的本构导数
仿照第一应力本构导数的推导,有
(6)
把(7.179)式代入(7.182)式,便得用应力导数表示的第二应力的本构导数为
(7)
其矩阵形式为
(8)
4. 大变形弹塑形本构方程
三维大变形弹塑形问题属于大变形问题。对于大应变问题, 流动理论关于应变增量等于弹性与塑性应变增量之和的假定一般是不成立的。但对于金属材料的大应变问题,弹性应变比塑性应变小得多,所以仍然可以近似的用增量流动理论。
大变形弹塑性加——卸载准则与小变形弹塑性加——卸载准则在形式上相同,只要把其中的应力认为是应力即可。大变形弹塑形速率型(或增量型)本构方程也与小变形弹塑性增量本构方程在形式上类似,只要把其中的应力和应变分别改用Cauchy应力的Jaumann导数和变形率,应力偏量和等效应力分别改用Cauchy应力偏量和等效的Cauchy应力便可。 本构方程的张量分量形式为
(9)
若令
(10)
(11)
则对各向同性材料有矩阵形式,即
(12)
考虑到弹性变形微小,塑性变形时体积不变,(由(4)式便得
(13)
(14)
大变形弹塑性静力隐式和静力显式有限元法
对于弹塑性大变形有限元而言,也象非线性弹性大变形问题一样,可以有TL和UL两种算法。这取决于本构方程用何种应力应变去描述。一般而言,能够用第二Poila应力和Green应变表示的本构方程,可以用TL或UL算法建立有限元方程。
1. 第一种作法(静力隐式UL列式理论基础)
(1) 增量形式本构方程
把(8)式和(12)式联立可得
(15)
由于其中和与时间无关,所以在上式两端乘以,并写成有限增量形式为
(16)
若把上式展开,并考虑到、的对称性和材料的各向同性,则可写成矩阵形式,即
(17)
其中
(18)
(2) 单元增量刚度方程
由于时刻的单元应变和单元应力为
(19)
(20)
所以把上述各式式虚功方程,有
(21)
于是仿照几何非线性三维大变形有限元法的U.L法,便得单元增量刚度方程,即
(22)
其中小位移弹塑性刚度矩阵
(23)
初应力刚度矩阵
(24)
大位移弹塑性刚度矩阵
(25)
初应力节点力
(26)
2. 第二种作法(静力显式有限元理论基础)
(1)持续平衡方程
把虚功率方程(5.233)式用于时刻和时刻两
您可能关注的文档
- (三)大气环流形势.doc
- (乙酸异戊酯)天那水MSDS.doc
- (人教版选修3-1)36带电粒子在匀强磁场中的运动.doc
- (北邮)运筹学课本答案.doc
- (投标用)施工组织设计.doc
- (新)第2章(2节)模拟题.doc
- (扩区)集中胶带大巷东揭煤措施3.doc
- (文献)青少年肥胖及代谢异常与INSIG2基因多态性.doc
- (铝)生产设备操作规程.doc
- (计算器).doc
- 2025年押题宝典演出经纪人之《演出经纪实务》题库附完整答案详解【精选题】.docx
- 2025年押题宝典演出经纪人之《演出市场政策与法律法规》题库及完整答案详解(考点梳理).docx
- 2025年湖南省初中学业水平考试仿真密卷-道德与法治(B卷).pptx
- 第3章 一次方程(组) 综合素质评价(含答案)2025-2026学年湘教版数学七年级上册.docx
- 2025年湖南省初中学业水平考试仿真密卷-英 语(B卷).pptx
- DB11_T1293.4-2015_卫生应急最小工作单元装备技术要求第4部分:核与辐射事故处置类_北京市.pdf
- DB11_T1322.68-2019_安全生产等级评定技术规范第68部分:设施蔬菜生产企业及专业合作社_北京市.pdf
- DB11_T1322.71-2018_安全生产等级评定技术规范第71部分:社会福利机构_北京市.pdf
- DB11_T1313-2015_薄抹灰外墙外保温用聚合物水泥砂浆应用技术规程_北京市.pdf
- DB11_T1320-2023_危险场所电气防爆安全检测技术规范_北京市.pdf
文档评论(0)