铸钢件凝固过程热应力场数值模拟研究.docVIP

1998年第19卷第4期华　北　工　学　院　学　报V o l.19　N o.4　1998 (总第64期JOURNAL OF NORTH CH INA INSTITUTE OF TECHNOLOG Y(Sum N o.64 铸钢件凝固过程热应力场数值模拟研究 ΞΞΞ 赵宇宏　　程　军　　徐　宏　　侯　华　　赵建平 　(华北工学院,太原030051 摘　要　用数值模拟方法研究了铸钢件凝固过程的热应力场,并计算了标准应力框的应力和位移,应力框 变形模拟结果与实际情况基本一致,为进一步研究热裂提供了依据. 关键词　热应力场;数值模拟;应力框 中图分类号　T G213 铸件凝固过程中热胀冷缩的差异,铸型或型芯对收缩的阻力以及在低温区域出现的相变,造成了铸件形成过程中的热应力、机械阻碍应力及相变应力[1].当应力值超过屈服极限或断裂极限时,将引起工件变形以致开裂,其中热应力多为残余应力.所以建立热应力场数学模型进行有限元分析,得到不同时刻热应力分布情况及变形情况,对于进一步制定热裂判据研究热裂位置及分布具有重要意义. 合金的物理性能、力学性能都是随温度变化的,如弹性模量E,线膨胀系数Α,以及屈服应力Ρs,强度极限Ρb等都是温度的函数,从而决定了铸造热应力计算的复杂性[2].目前铸造应力的计算模型有纯弹性模型、弹性 蠕变模型和弹塑性 蠕变模型等[3].本文采用精度较高的弹塑性模型,在准固态区间考虑粘性因素. 数值模拟方法主要有三种:有限差分法、有限元法和边界元法.有限元法单元划分形状限制小,能更好地逼近铸件形状,并且有限元法适合应力场数值模拟计算,所以选用有限元法开发此软件. 1　热应力场计算数学模型及模拟方法 111　概述 (1建立作用于单元上的结点力和结点位移之间的关系式,即单元的平衡方程,用矩阵表示 k e?e=F e(1式中　k e为单元刚度矩阵;?e为单元结点位移向量;F e为单元热负荷向量;e为单元. (2集合所有单元的平衡方程,建立整个结构的平衡方程得到 k?=F 　　考虑边界条件后,解出所有的未知结点位移?. (3几何方程[4]　Ε=B?e　　Ε为应变列矩阵;B为单元几何矩阵. (4本构方程[4]Ρ=DΕ　　Ρ为应力列矩阵;D为与单元材料有关的弹性矩阵. 112　单元刚度矩阵分析 11211　准固态区间的单元刚度矩阵分析 根据准固态区材料的流变模型,可得三维问题应力应变本构关系 D0Ε=Ρ D0=(D1-1+D2-1-1 Ξ 收稿日期:1998207209　第一作者　女　1974年出生　硕士生 ΞΞ 式中　D 0为准固态区间与粘性有关的弹性矩阵;D 1为粘弹性单元的切线模量矩阵;D 2为粘塑性单元的切线模量矩阵.有了D 0,即可通过常规则方法求得单元刚度矩阵k e ,为了推导D 0,可从D 1和D 2入手.1.2.1.1　求粘弹性单元切线模量矩阵D 1 本单元采用三参数标准粘弹性模型,推导可得微分型本构方程[5] Ρ+(1ΒΡ =g 0Ε+(g 0+g 1Β Ε g 0= E 1E 2E 1+E 2　　g 1=E 1+g 0　　Β= E 1+E 2Γ 式中　Ρ? 为应力率;Ε? 为应变率.Ρ为应力;Ε为应变;E 1,E 2为杨氏弹性模量;Γ为粘性系数. 假设粘弹性过程对体积没有影响以及时间步增量?t 很小,可得应力偏量增量?Ρ与应变偏量增量?Ε的关系　　　　　　　　　　?Ρ=2g ?Ε 式中　g =g 0+g 1(1-e -Β?t Β?t 则粘弹性单元切线模量矩阵的显式表达式为 D 1= Η2Η1Η1 000Η1 Η2Η1000Η1Η1Η20 00000g 00000g g 式中　Η1=K -(2 3g ;Η2=K +(4 3g ;K 为体弹性模量.1.2.1.2　求粘塑性单元的切线模量矩阵D 2 采用玻任纳本构方程,由虎克定律和粘塑性流动法则得本构关系为 Ε =D -1 Ρ +Χ[5(f f 0]( 5f Ρ 　　f =2J ′2+R 　　f 0=R J ′ 2= 16 〔(Ρx -Ρy 2+(Ρy -Ρz 2+(Ρz -Ρx 2〕+Σ2x y +Σ2y z +Σ2 z x 式中　Χ为材料常数;J ′2为应力偏量第二不变量;R 为屈服半径 .得出本构关系后,通过分析Ε和Ρ之间的关系来求得D 2.11212　固相线以下的弹塑性问题 一般情况下,无法建立起最终的应力状态和应变状态之间的全量关系,而只能建立起反映加载路径的应力应变之间的增量关系,增量关系包括加载和卸载过程. 应力应变的增量形式本构关系为 d Ρ=D 3d Ε　　D 3=D 4-D 5式中　D 3为固相线以下后弹塑性矩阵;　D 4为固相线以下后弹性矩阵;　D 5为固相线以下后塑性矩阵.推导得D 3=K ′ (Ρ1x 2Ρ1x Ρ1y Ρ1x Ρ1z Ρ1x Σ