ANSYS经典版残余应力计算(含APDL命令流)

ANSYS经典版残余应力计算（含APDL命令流）
目录
简介： (2)
附件1 钢轨建模 (2)
附件2 构造残余应力文件 (6)
附件3 读入残余应力文件并云图显示残余应力 (8)
简介：
ANSYS 经典版残余应力计算，首先建立钢轨模型，然后构造残余应力文件，之后导入残余应力文件，并云图显示残余应力。

图1 残余应力云图
附件1 钢轨建模
! ***************环境设置************************
finish
/clear
/filn,file
/title,
MN MX
X
Y
Z
15.416730.779946.143261.506476.869792.2329107.596122.959
/PREP7
ET,1,185
mp,ex,1,2e5 ! 定义材料1的弹性模量mp,nuxy,1,0.3 ! 定义材料1的泊松比
TB,BISO,1,1,2
TBTEMP,0
TBDATA,,275 !材料屈服强度
*SET,H1,30.5
*SET,H2,48.5
*SET,H3,48.5
*SET,B1,45.7
*SET,B2,29.25
*SET,B3,16.5
*SET,B4,36.5
*SET,B5,10
*SET,XD1,1/3
*SET,XD2,1/9
*SET,XD3,1/3
*SET,XD4,1/20
*SET,R0,400
*SET,R1,2
*SET,R2,4
*SET,R3,40
*SET,R4,20
*SET,R5,25
*SET,R6,8
*SET,R7,13
*SET,R8,80
*SET,R9,300
*SET,R10,2
*SET,FAI,43 ! R6 圆心角
*SET,H,H1+2*H2+H3
K,1,
K,2,B1+B2
K,3,B1+B2,H1-B1*XD1-B2*XD2
K,4,B1,H1-B1*XD1
K,5,0,H1
*DO,I,1,4
L,I,I+1
*ENDDO
*SET,CTA,ATAN(XD3)
*SET,BTA,(90-FAI)*ACOS(-1)/180-CTA
*SET,A,1/TAN(CTA)
*SET,C1,B3/2+R0-R6*SIN(CTA)-(R5-R6)*COS(BTA)
*SET,C2,R5*SIN(BTA)+R6*COS(CTA)-R6*SIN(BTA)-H2
*SET,C,(C1-C2/TAN(CTA))/(R0-R5)
*SET,C3,A*C-SQRT(A*A-C*C+1)
*SET,C3,C3/(A*A+1)
*SET,REFA,ASIN(C3)
*SET,C1,R6*SIN(CTA)+(R5-R6)*COS(BTA)+(R0-R5)*COS(REFA) *SET,ROUX,B3/2+R0-C1
*SET,ROUY,ROUX*TAN(CTA)
*SET,ANGI,ASIN((H1+H2)/R0)
K,6,B3/2+R0-R0*COS(ANGI)
K,7,B3/2+R0-R0*COS(REFA),H1+H2+R0*SIN(REFA)
LARC,6,7,2,R0
*SET,DX1,R6*(COS(BTA)-SIN(CTA))
*SET,DY1,R6*(COS(CTA)-SIN(BTA))
K,8,ROUX-DX1,H1+2*H2+ROUY-DY1 K,9,ROUX,H1+2*H2+ROUY
LARC,7,8,2,R5
LARC,8,9,2,R6
K,10,B4,H1+2*H2+B4*XD3
K,11,B4-(H3-B4*XD3)*XD4,H
L,9,10
L,10,11
*SET,ANGI,ASIN(B5/R9)
K,12,0,H
K,13,B5,H-(R9-R9*COS(ANGI)) LARC,12,13,1,R9
*SET,DX1,B4-B5+R8*SIN(ANGI)
*SET,ANGI1,ASIN(DX1/R8)
*SET,Y1,R8*COS(ANGI)-R8*COS(ANGI1) K,14,B4,KY(13)-Y1
LARC,13,14,1,R8
LFILLT,1,2,R1
LFILLT,2,3,R2
LFILLT,3,4,R3
LPTN,4,5
LFILLT,17,18,R4
LFILLT,8,9,R10
LPTN,9,11
LFILLT,20,22,R7
LDELE,15,16,1,1
LDELE,19,21,2,1
LSYMM,X,ALL
NUMMRG,ALL
NUMCMP,ALL
AL,ALL
vext,1,,,,,200
et,2,63
aesize,1,5
amesh,1
lesize,77,5
vsweep,1
aclear,all
etdele,2
NUMMRG,ALL NUMCMP,ALL
allsel,all
finish
SAVE,jm,db
附件2 构造残余应力文件
/prep7
RESUME,jm,db
!约束
/solu
asel,s,,,3,36,33
da,all,all
allsel,all
!计算设置
!设定输出选项
asel,s,,,18,21,1
sfa,all,,pres,100
allsel,all
OUTRES,ALL,ALL ! 输出所有结果TIME,1 ! 设定时间步
NSUB,100 ! 设定子步数
KBC,0 !线性加载
NROPT,FULL ! 将N-L选项设置为FULL PRED,ON ! 打开时间步长预测LNSRCH,ON ! 打开线性搜索AUTOTS,ON ! 激活自动时间步长
!NLGEOM,ON ! 打开大变形
ALLSEL ! 全选择
lswrite,1
/solu
asel,s,,,18,21,1
sfa,all,,pres,0
allsel,all
!设定输出选项
OUTRES,ALL,ALL ! 输出所有结果
TIME,2 ! 设定时间步
NSUB,100 ! 设定子步数
KBC,0 !线性加载
NROPT,FULL ! 将N-L选项设置为FULL PRED,ON ! 打开时间步长预测LNSRCH,ON ! 打开线性搜索
AUTOTS,ON ! 激活自动时间步长
!NLGEOM,ON ! 打开大变形
ALLSEL ! 全选择
lswrite,2
INISTATE, WRITE, 1, , , , 0, S !形成ist文件lsSOLVE,1,2,1
FINI
附件3 读入残余应力文件并云图显示残余应力
/prep7
RESUME,jm,db
!约束
/solu
asel,s,,,3,36,33
da,all,all
allsel,all
INISTATE, read, file,ist !导入ist文件
solve
/post1
plnsol,s,eqv。