大应变热弹塑性蠕变分析研究.pdfVIP

大应变热弹塑性蠕变分析 ‘ 陈火红 刘长乱 席源山 (航空工业总公司第六二八研究所) [摘要1本文结合工程实际，对大应变热弹塑性蠕变有限元分析进行一些探讨。采用更新拉 格朗日列式作为参照系统来描述大应变问题几何关系。在本构关系数值积分方面采用参数 积分法，分别采用牛顿一拉斐逊法进行应力各分量迭代和等效应力函数法对等效应力进行迭 代。编制了有限元程序并集成人非线性分析程序AROLANS之中，最后给出一个工程例题验 证算法及程序的可靠性 也证实了等效应力函数法可以大幅度提高计算效率。 关健词:大应变 几何非线性 蠕令 M;E件 I 引言 在工程实际中，常会遇到大应变热弹塑性蠕变问题，如长期在高温环境下工作的承力结 构、金属超塑性成形、压力加工过程等。对这些结构进行应力、应变分析及加工成形过程进行 数值模拟时，需要采用几何非线性来描述位移与应变之间的关系及合适的参照系统，并且需要 效率高、稳定性好的本构方程积分方法。 本文先讨论大应变下的有限元控制方程，接着给出了由文献仁21、[3]提出的二种不同的热 弹塑性蠕变本构方程的积分处理方法，编制了有限元程序，最后给出一个金属成形过程数值模 拟的工程实例，证实了采用等效应力函数法可以大幅度提高计算效率。 2 几何非线性有限元方程 大应变问题在力学上属于高度几何非线性问题，用传统的工程应变来描述超塑成形的变 形情况已经不适宜了。现在，力学上普遍采用了两种参照系统来描述大变形、大转动和大应变 的问题。一种是空间描述，即欧拉描述;再一种是物质描述，即拉格朗日描述，包括全拉格朗日 描述和更新拉格朗日描述。本文中，采用更新的拉格朗日描述(UpdatedLagangianFoanula- tion)t`i，它是以平衡时刻 的‘物体构形作为物体的参照构形，即认为在t时刻的物体是未变 形”物体，去求解下一时刻 +‘,,At时物体的各个未知变量。 因此 ，用更新拉格朗日法所定义的格林应变应为: △丙=2lr1aa}`xu;;a·a;`x,·aIAu`x,iaa·p`uzt; 。航空科学基金、人事部留学回国人员科技活动 D类资助项 目 二A城 +Q}E4 (1) 式中下标t表示参照构形选在。时刻。 用更新拉格朗日法所定义的KimJihoff应力在 十‘乙:时刻写成 }11}Se=:。，+4,}凡 (2) 式中，:，为:时刻的真实应力，也就是柯西应力，lli.Sv为参照:时刻物体的Kirchhoff应力。 有限元的最后控制方程为: (:仁K)o+,[K],)IAgI二JKIr)L9) 了 ‘ q 、 ! 二}F)+}T}一{P) 、 ︸ 厂 了 、 奋 连 I :[K]a=Ef,。,[BI[C],[B]dv 、 J 了 ‘ 气 、 !fK一“艺丁:，}x]丁[al}ax]·‘ ︺ “山:二’二+乙，山‘、