粘塑性材料结构的有限元分析方法.pdfVIP

第 13 卷　第 4 期 航空动力学报 V o l13 N o 4 1998 年 10 月 Journa l of Aerospace Power O ct. 　1998 粘塑性材料结构的有限元分析方法 北京航空航天大学　　杨晓光 　耿　瑞 【摘要】　本文针对W alker 粘塑性本构方程, 通过发展一自适应变步长积分算法, 用户子程序的方 法把该本构方程结合进了国际著名的非线性有限元软件ABAQU S, 并进行了算例验证。它包括了 四种单元类型, 可考虑材料变形的率相关, 加载的循环、非等温、变温等情况。适用于发动机热端部 件非线性结构分析。 　主题词: 　粘塑性　本构方程　有限元 　分类号: 　 231 9　 346 V O 1　引　言 [ 1 ] 　　很多的研究及应用表明 : 目前基于经典弹塑性和蠕变理论的非线性结构分析方法不足 以正确设计、改进和分析航空航天及核动力设备关键热端零件的强度及疲劳寿命。其主要原因 之一是结构承受各种复杂恶劣的循环加载和瞬态加载条件, 这些条件使材料经受了与应力、温 度和时间相关的并且是耦合的非线性行为。基于经典弹塑性及蠕变理论的本构方程, 通常是把 非弹性应变分为与时间相关的蠕变和与时间无关的塑性。但是大量的实验表明: 在高温下其本 质上是时间相关的, 并且两者有交互作用。这表明非弹性变形是由一单一的机理在控制, 宜用 统一的方法, 即把塑性及蠕变相联系起来的弹粘塑性本构方程来处理。 目前, 这类统一型的本构理论在近 20 年得到了很大的发展, 已有十多种方法。尤其是与有 限元结构分析方法后, 逐渐在实际的工程分析得到了应用。这些本构模型大都根据不可逆热力 学及位错动力学的研究成果, 恰当地选择单个或多个内变量及其演变方程用以描述非弹性变 形的诸方面, 如塑性、蠕变、松驰等。本文将根据W alker 粘塑性本构理论, 通过发展自适应显 式积分方法, 把该本构方程结合进通用的大型非线性有限元程序ABAQU S, 并进行了该本构 理论率相关, 循环, 热恢复及非等温等分析能力的验证。 2　 粘塑性本构理论 W a lker 　　一般应变率可分为弹性应变率和非弹性应变率, 并且弹性部分遵循Hooke 定律。 e in ( 1) = + ij ij ij e ( ) ( ) ( ) = [ 1 + E ] - E 2 ij ij kk ij 式中: 为杨氏模量, 为泊松比, ij 为 符号。非弹性部分根据对材料硬化和恢复的 E K ronecker 不同处理将有不同的形式。 1998 年 5 月收稿; 1998 年 6 月收到修改稿。 男　36 岁　教授　北京航空航天大学 405 教研室　100083 © 1995-2005 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 第　4　期 粘塑性材料结构的有限元分析方法