- 1、本文档共25页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
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
- 河北网通可视化数据挖掘系统的设计与实现:技术、应用与优化.docx
- 昆山农村商业银行信用风险管理:问题剖析与优化路径.docx
- 水溶性芯材微胶囊化:制备、性能与应用探索.docx
- 金车公司人力资源战略:困境、突破与展望.docx
- 创新型中小企业融资结构对投资行为的影响:基于多维度视角的剖析.docx
- 电力国企应收账款质量:评价体系构建与管控策略探究.docx
- 产业技术分析中产业政策分析方法的多维探究与实践.docx
- 预埋止水带位置偏差控制技术专题.pptx
- 冰片对小儿难治性癫痫血清和脑脊液丙戊酸钠浓度影响的深度剖析.docx
- 数字时代下国际展览产业联盟(阿里巴巴)运营模式的创新与发展研究.docx
文档评论(0)