求助：热应力计算出现的问题（附命令流）

qingshuihuayu

我在运行以下程序时，总是出现以下问题算不出结果，请问大家我该如何处理？十分感谢！！no results(RTH) available,the leread command is ignored.
命令流如下：
/title,pile
/units,si
/prep7
et,1,plane55         
mp,kxx,1, 0.488*24*3600
mp,kxx,2,0.972*24*3600      

!建立模型
K,1,-10,0
K,2,10,0
K,3,10,-20
K,4,-10,-20
K,5,-0.5,0
K,6,0.5,0
K,7,-0.5,-8.6
K,8,0.5,-8.6
K,9,-10,-8.6
K,10,10,-8.6
A,1,5,7,9
A,5,6,8,7
A,6,2,10,8
A,9,10,3,4
AADD,1,3,4
allsel,all
numcmp,all  

!赋属性并划分网格
ASEL,S,area,,1
AATT,2,1,1
lsla,s
lesize,all,0.2
mshape,0,2d
mshkey,1
amesh,all
ASEL,S,area,,2
AATT,1,1,1
lsla,s
lesize,all,0.2
mshape,0,2d
mshkey,2
amesh,all

!合并侧面土体上的重复单元及节点
allsel,all
aplot
esla,s
numcmp,all

/solu
allsel,all
NSLa,S,1
ic,ALL,TEMP,0

!施加边界条件
lSEL,S,LOC,y,-20
NSLl,S,1
F,all,heat,-12009.6

lSEL,S,LOC,X,10
lSEL,A,LOC,x,-10
NSLl,S,1
SF,ALL,HFLUX,0
lSEL,S,LOC,y,0
NSLl,S,1
D,ALL,TEMP,-1
ALLSEL,ALL
FINISH

!求解
/SOLVE
save
FINISH
/POST1
PLNSOL,TEMP
FINISH
/PREP7
ETCHG,TTS
mp,ex,1,2.5e8
mp,prxy,1,0.4
uimp,1,aply,,,82.48e-4
mp,ex,2,2.5e10
mp,prxy,2,0.2
uimp,2,aply,,,1.0e-5
FINISH
/SOLU
lSEL,S,LOC,y,-20
NSLl,S,1
D,all,ALL
lSEL,S,LOC,X,10
lSEL,A,LOC,x,-10
NSLl,S,1
D,ALL,UX,0
ALLSEL,ALL
LDREAD,temp,,,,,file.rst
TREF,0
SOLVE
FINISH

shaozhantao

!求解
/solu
solve
FINISH
/POST1
PLNSOL,TEMP
FINISH
/PREP7
ETCHG,TTS
mp,ex,1,2.5e8
mp,prxy,1,0.4
uimp,1,aply,,,82.48e-4
mp,ex,2,2.5e10
mp,prxy,2,0.2
uimp,2,aply,,,1.0e-5
FINISH
/SOLU
lSEL,S,LOC,y,-20
NSLl,S,1
D,all,ALL
lSEL,S,LOC,X,10
lSEL,A,LOC,x,-10
NSLl,S,1
D,ALL,UX,0
ALLSEL,ALL
LDREAD,temp,,,,,file.rst
TREF,0
SOLVE
FINISH

zhangyp

/title,pile
/units,si
/prep7
et,1,plane55         
mp,kxx,1, 0.488*24*3600
mp,kxx,2,0.972*24*3600      
!建立模型
K,1,-10,0
K,2,10,0
K,3,10,-20
K,4,-10,-20
K,5,-0.5,0
K,6,0.5,0
K,7,-0.5,-8.6
K,8,0.5,-8.6
K,9,-10,-8.6
K,10,10,-8.6
A,1,5,7,9
A,5,6,8,7
A,6,2,10,8
A,9,10,3,4
AADD,1,3,4
allsel,all
numcmp,all  
!赋属性并划分网格
ASEL,S,area,,1
AATT,2,1,1
lsla,s
lesize,all,0.2
mshape,0,2d
mshkey,1
amesh,all
ASEL,S,area,,2
AATT,1,1,1
lsla,s
lesize,all,0.2
mshape,0,2d
mshkey,2
amesh,all
!合并侧面土体上的重复单元及节点
allsel,all
aplot
esla,s
numcmp,all
/solu
allsel,all
NSLa,S,1
ic,ALL,TEMP,0
!施加边界条件
lSEL,S,LOC,y,-20
NSLl,S,1
F,all,heat,-12009.6
lSEL,S,LOC,X,10
lSEL,A,LOC,x,-10
NSLl,S,1
SF,ALL,HFLUX,0
lSEL,S,LOC,y,0
NSLl,S,1
D,ALL,TEMP,-1
ALLSEL,ALL
FINISH
!求解
/SOLU
solve
FINISH
/POST1
PLNSOL,TEMP
FINISH
/PREP7
ETCHG,TTS
mp,ex,1,2.5e8
mp,prxy,1,0.4
uimp,1,aply,,,82.48e-4
mp,ex,2,2.5e10
mp,prxy,2,0.2
uimp,2,aply,,,1.0e-5
FINISH
/SOLU
lSEL,S,LOC,y,-20
NSLl,S,1
D,all,ALL
lSEL,S,LOC,X,10
lSEL,A,LOC,x,-10
NSLl,S,1
D,ALL,UX,0
ALLSEL,ALL
ldread,temp,,,,,,rth
TREF,0
SOLVE
FINISH

qingshuihuayu

谢谢大家的热心，我先谢了。我再试试看

qingshuihuayu

为什么我的命令流计算完了以后，出不来材料的热胀冷缩特性呢？在结果分析力，温度分布有了，可是没有变形和应力？

zhangyp

没定义热膨胀系数，如：
mp,alpx,1,11e-6

qingshuihuayu

为什么我的命令流计算完了以后，出不来材料的热胀冷缩特性呢？在结果分析力，温度分布有了，可是没有变形和应力？
没定义热膨胀系数，如：
mp,alpx,1,11e-6
定义完了后
/title,pile
/units,si
/prep7
et,1,plane55         
mp,kxx,1, 0.488*24*3600
mp,kxx,2,0.972*24*3600      
!建立模型
K,1,-10,0
K,2,10,0
K,3,10,-20
K,4,-10,-20
K,5,-0.5,0
K,6,0.5,0
K,7,-0.5,-8.6
K,8,0.5,-8.6
K,9,-10,-8.6
K,10,10,-8.6
A,1,5,7,9
A,5,6,8,7
A,6,2,10,8
A,9,10,3,4
AADD,1,3,4
allsel,all
numcmp,all  
!赋属性并划分网格
ASEL,S,area,,1
AATT,2,1,1
lsla,s
lesize,all,0.2
mshape,0,2d
mshkey,1
amesh,all
ASEL,S,area,,2
AATT,1,1,1
lsla,s
lesize,all,0.2
mshape,0,2d
mshkey,2
amesh,all
!合并侧面土体上的重复单元及节点
allsel,all
aplot
esla,s
numcmp,all
/solu
allsel,all
NSLa,S,1
ic,ALL,TEMP,0
!施加边界条件
lSEL,S,LOC,y,-20
NSLl,S,1
F,all,heat,-12009.6
lSEL,S,LOC,X,10
lSEL,A,LOC,x,-10
NSLl,S,1
SF,ALL,HFLUX,0
lSEL,S,LOC,y,0
NSLl,S,1
D,ALL,TEMP,-1
ALLSEL,ALL
FINISH
!求解
/SOLU
solve
FINISH
/POST1
PLNSOL,TEMP
FINISH
/PREP7
ETCHG,TTS
mp,ex,1,2.5e8
mp,prxy,1,0.4
mp,1,aply,,,82.48e-4
mp,ex,2,2.5e10
mp,prxy,2,0.2
mp,2,aply,,,1.0e-5
FINISH
/SOLU
lSEL,S,LOC,y,-20
NSLl,S,1
D,all,ALL
lSEL,S,LOC,X,10
lSEL,A,LOC,x,-10
NSLl,S,1
D,ALL,UX,0
ALLSEL,ALL
ldread,temp,,,,,,rth
TREF,0
SOLVE
FINISH


出现
*** ERROR ***                           CP=     44.609   TIME= 09:02:30
Unknown material property label= 1                                      
  For nonlinear materials try using the TB command.                     
  The MP command is ignored.
*** ERROR ***                           CP=     44.672   TIME= 09:02:37
Unknown material property label= 2                                      
  For nonlinear materials try using the TB command.                     
  The MP command is ignored.  
请问如何修改，十分感谢！

高山亭子

mp,1,aply,,,82.48e-4 mp,2,aply,,,1.0e-5 FINISH

mp,aplx,1,82.48e-4
mp,aplx,2,1.0e-5

qingshuihuayu

还是一样的！不可以！