单晶滑移模型弹塑性蠕变计算有限元程序考核研究.pdfVIP

单晶滑移模型弹塑性蠕变计算有限元程序考核 成晓鸣 尹泽勇 杨治国 (株洲航空动力机械研究所 412002) 岳珠峰 (西北工业大学工程力学系 710072) 摘【要〕本文对单品滑移模型弹塑性有限元计算程序进行了考核，考核内容包括经典滑移模 型峨变计算、塑性计算、弹塑性峨变计算、11}环加教娇变计算及中间温度峨变插值计算· 经典滑移模型蠕变计算考核 首先给出均匀拉仲杆蠕变应变的理论解。 设材料主轴坐标系下的应力为13ii，则滑移系二p的剪应力:“，为 :”·=33。。naa!p (1) i=I j=t 其中Iti、1尹分别为滑移面法向和滑移方向单位矢呈。 根据陈化理论，滑移系ap的姗变剪应变为 子p yap二B(}ep卜”to (2) }州 式((2)中B.fl,p为滑移系的蠕变参数. 滑移系蠕变剪应变率为 (为简便起见，省略上标，9) P-1 p T (3) }:! 可以采用差分法求解方程(3)，但由于:=。时Y,=00，给差分求解造成一定困难。设: 随时间变化的规律为:一1“借母(一‘，,)=a,+st(一t,)娜tst2),直接求解(3)式得 Y}t()二}spPnB(h})IP、二p」。 (4) 在;二au (t,sc5t2)情况下有 :(!一‘卜}aX战一)·Itttc}V·hI (5) (4),(5)两式中c为待定系数。 将滑移系蠕变剪应变转换成材料主轴坐标系下蠕变应变为 伞。)c YP((沙ntl;ena) (6) 二矍+E=Ip,u,22I 式(6)中na,n3(根据不同的滑移系族取不同的值。 考题1:经典滑移模型十二面体滑移系开动蠕变计算 取轴线方向与11001取向重合的拉杆为计算模型，杆长 60mm。杆所受拉应力为 a=100MPa，取十二面体滑移系蠕变参数为:B=U237X10$,n=2.3427,p=0.579. 蠕变时问为2小时，分为 10个时间步，每个时间步步长均为0.2小时。加载函数为 r-o.o,a=0.0;t--0.2,a=100.0;t=2.0,a=100.0;各时刻轴向蠕变应变。。及轴向位移u}程序计算 值及理论解见表1. 表t 十二面体滑移系开动蠕变计算结果(eax10-,u,x10mm) 如) 0.2 0.4 0.6 0名 1.0 、2 1.4 1.6 l.8 2.0 ‘亡计算值 1.572 4.985 7.169 8.949 10.50 11.90 13.19 14.40 15.53 16.61 ￡。理论值 1.812 5.127 7之83 9.049 10.59 11.99 13.27 14.47 15.60 16.68