模拟单个试件在一变温度载荷条件下的温度场分布。
温度加载点为-1<X<1的所有Node
计算命令流如下:
/PREP7
ET,2,SOLID227,11
MP,EX,1,2e14 !young
MP,PRXY,1,0.29 !possion ratio
MP,DENS,1,7740 !density
MP,KXX,1,26.7692 !thermal conductivity
MP,C,1,477.1424 !Specify Cp
MP,ALPX,1,2e-6
MP,REFT,1,288
FINISH
/SOLU
LSCLEAR, ALL
ANTYPE,4
TRNOPT,FULL
LUMPM,0
TUNIF,288.15 ! set uniform temp for every node
TREF,288.15
SFA,1,1,CONV,0.1,288
SFA,2,1,CONV,0.1,288
SFA,3,1,CONV,0.1,288
SFA,4,1,CONV,0.1,288
SFA,5,1,CONV,0.1,288
SFA,6,1,CONV,0.1,288
Temp_1=291.983
tinc_1=2
tim_1=2
inum=0
*DO,inum,1,6
NSEL,ALL
NSEL,S,LOC,X,-1,1
D,ALL,TEMP,Temp_1
TIME,tim_1
AUTOTS,-1
DELTIM, , , ,1
KBC,0
ALLSEL
solve
tim_1=tim_1+tinc_1
*ENDDO
tinc_2=1.5
tim_2=13.5
inum=0
*DO,inum,1,55
Temp=0.000034281279*tim_2**3-0.008534289983*tim_2**2+0.836829421614*tim_2+290.068064269667
NSEL,ALL
NSEL,S,LOC,X,-1,1
D,ALL,TEMP,Temp
TIME,tim_2
AUTOTS,-1
DELTIM, , , ,1
KBC,0
ALLSEL
solve
tim_2=tim_2+tinc_2
*ENDDO
tinc_3=0.8
tim_3=95.3
inum=0
*DO,inum,1,7
Temp=0.14021351707*tim_3**2-29.12643848650*tim_3+1824.65007273490
NSEL,ALL
NSEL,S,LOC,X,-1,1
D,ALL,TEMP,Temp
TIME,tim_3
AUTOTS,-1
DELTIM, , , ,1
KBC,0
ALLSEL
solve
tim_3=tim_3+tinc_3
*ENDDO
FINISH
计算截图如下:
问题1: 为什么加载温度的节点处产生了膨胀?
问题2:温度场分布明显过于均匀,是因为时间步长过小吗?
PS:本试件材料为钢材,总长度10cm
希望高手帮忙看下命令流,多谢指教。 |