我先用ansys进行热分析,施加温度曲线,这时可以得到不同时刻的温度分布。我想直接将温度结果用LDREAD读入结构分析中,采用了循环,可是不同时刻的应力结果还是一样,请指教。下面附上命令流
fini
/CLEAR
!===================第一步进行稳态热分析======================
/PREP7
ET,1,SOLID70
MP,EX,1,3.45E10
MP,PRXY,1,0.2
MP,DENS,1,2500
MP,KXX,1,2.3
ZH=1.0
BLOCK,,1,,1,,ZH, !厚度暂取1.0
ALLSEL,ALL
LESIZE,ALL,0.2, , , , , , ,1
VMESH,ALL
/SOLU
ANTYPE,trans
trnopt,full
lumpm,0
timint,on
time,12
deltim,0.2
AUTOTS,1
*DIM,wendukx,TABLE,8,1,1
*SET,wendukx(1,0,1) , 0 !小时
*SET,wendukx(2,0,1) , 1
*SET,wendukx(3,0,1) , 2
*SET,wendukx(4,0,1) , 4
*SET,wendukx(5,0,1) , 6
*SET,wendukx(6,0,1) , 8
*SET,wendukx(7,0,1) , 10
*SET,wendukx(8,0,1) , 12
*SET,wendukx(1,1,1) , 8 !混凝土外侧的传热系数。
*SET,wendukx(2,1,1) , 8
*SET,wendukx(3,1,1) , 8
*SET,wendukx(4,1,1) , 8
*SET,wendukx(5,1,1) , 8
*SET,wendukx(6,1,1) , 8
*SET,wendukx(7,1,1) , 8
*SET,wendukx(8,1,1) , 8
*DIM,wendu,TABLE,8,1,1
*SET,wendu(1,0,1) , 0 !小时
*SET,wendu(2,0,1) , 1
*SET,wendu(3,0,1) , 2
*SET,wendu(4,0,1) , 4
*SET,wendu(5,0,1) , 6
*SET,wendu(6,0,1) , 8
*SET,wendu(7,0,1) , 10
*SET,wendu(8,0,1) , 12
*SET,wendu(1,1,1) , 170 !温度
*SET,wendu(2,1,1) , 145
*SET,wendu(3,1,1) , 120
*SET,wendu(4,1,1) , 120
*SET,wendu(5,1,1) , 100
*SET,wendu(6,1,1) , 80
*SET,wendu(7,1,1) , 60
*SET,wendu(8,1,1) , 50
ASEL,S,LOC,Z,ZH
SFA,ALL,1,CONV,%wendukx%,%wendu% !此处将%wendukx%换成上面数组的数值8,计算就会给错误,看了命令帮助,不知道原因
ASEL,S,LOC,Z,0
SFA,ALL,1,CONV,10e6,40
ALLSEL,ALL
outres,all,all
SOLVE
!===================第二步进行结构分析======================
/PREP7
ALLSEL,ALL
ETCHG,TTS
ET,1,SOLID45
MP,EX,1,3.45E10
MP,PRXY,1,0.2
MP,DENS,1,2500
MP,ALPX,1,1e-5
/SOL
ANTYPE,trans
trnopt,full
lumpm,0
timint,on
AUTOTS,1
ALLSEL,ALL
SFADELE,ALL
tref,40
NSEL,S,LOC,z,0
D,ALL,ALL
ALLSEL,ALL
*do,I,0.2,12,0.2
time,I
LDREAD,TEMP,,, ,I ,'001','rth',' ' !读入不同时刻温度计算
allsel,all
outres,all,all
solve
*enddo
fini
/POST1
!by picks 读入不同的结果,应力确都相同
!节点应力结果在每个子步都相同?
|