求助帮忙 急

ymz.124

在做一个钢管混凝土温度应力，从室温加载到1000，还要得到200,400,600,800不同温度点下的温度场温度曲线和温度应力变化曲线。
  没有热学基础，本想这几个温度分别加载，但是还是有很多问题，程序看的都晕了，时间很紧，希望哪位朋友加我QQ多多请教！非常感谢！拜托各位了！

ymz.124

fini
/clear

/PREP7  

ET,1,SOLID70

ET,2,SOLID70

/NOPR   
/PMETH,OFF,0
KEYW,PR_SET,1   
KEYW,PR_STRUC,1
KEYW,PR_THERM,1
KEYW,PR_FLUID,0
KEYW,PR_ELMAG,0
KEYW,MAGNOD,0   
KEYW,MAGEDG,0   
KEYW,MAGHFE,0   
KEYW,MAGELC,0   
KEYW,PR_MULTI,1
KEYW,PR_CFD,0   
/GO

MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,2.06E+005  
MPDATA,PRXY,1,,0.3  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,7850
MPTEMP,,,,,,,,  
MPTEMP,1,20
MPTEMP,2,200
MPTEMP,3,400
MPTEMP,4,600
MPTEMP,5,800
MPTEMP,6,900
MPTEMP,7,1000   
MPDATA,KXX,1,,47.56
MPDATA,KXX,1,,43.6  
MPDATA,KXX,1,,39.2  
MPDATA,KXX,1,,34.8  
MPDATA,KXX,1,,30.4  
MPDATA,KXX,1,,28.3  
MPDATA,KXX,1,,28.3  
MPTEMP,,,,,,,,  
MPTEMP,1,20
MPTEMP,2,200
MPTEMP,3,400
MPTEMP,4,600
MPTEMP,5,800
MPDATA,C,1,,481.82  
MPDATA,C,1,,513.49  
MPDATA,C,1,,609.45  
MPDATA,C,1,,769.39  
MPDATA,C,1,,993.31  
MPTEMP,,,,,,,,  
MPTEMP,1,20
MPTEMP,2,200
MPTEMP,3,400
MPTEMP,4,600
MPTEMP,5,800
MPTEMP,6,1000   
UIMP,1,REFT,,,  
MPDATA,ALPX,1,,1.21E-005
MPDATA,ALPX,1,,1.28E-005
MPDATA,ALPX,1,,1.36E-005
MPDATA,ALPX,1,,1.44E-005
MPDATA,ALPX,1,,1.52E-005
MPDATA,ALPX,1,,1.6E-005
TB,MELA,1,1,4,  
TBTEMP,0

TBPT,,0.002,412
TBPT,,0.003,412
TBPT,,0.004,412
TBPT,,0.005,412
MPTEMP,,,,,,,,  
MPTEMP,1,0
MPDATA,EX,2,,30000  
MPDATA,PRXY,2,,0.2      
TB,MELA,2,1,5,  
TBTEMP,0
TBPT,,0.0005,15
TBPT,,0.001,24
TBPT,,0.002,30
TBPT,,0.0025,29
TBPT,,0.0038,22
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,2,,2400
MPTEMP,,,,,,,,  
MPTEMP,1,20
MPTEMP,2,200
MPTEMP,3,293
MPTEMP,4,400
MPTEMP,5,600
MPTEMP,6,800
MPDATA,KXX,2,,1.355
MPDATA,KXX,2,,1.355
MPDATA,KXX,2,,1.335
MPDATA,KXX,2,,1.22  
MPDATA,KXX,2,,0.972
MPDATA,KXX,2,,0.724
MPTEMP,,,,,,,,  
MPTEMP,1, 20
MPTEMP,2,200
MPTEMP,3,400
MPTEMP,4,600
MPTEMP,5,800
MPDATA,C,2,,913.22  
MPDATA,C,2,,890.56  
MPDATA,C,2,,858.89  
MPDATA,C,2,,1200.00
MPDATA,C,2,,1388.90  
MPTEMP,,,,,,,,  
MPTEMP,1,20
MPTEMP,2,200
MPTEMP,3,400
MPTEMP,4,600
MPTEMP,5,800
MPTEMP,6,1000   
UIMP,2,REFT,,,  
MPDATA,ALPX,2,,6.2E-006
MPDATA,ALPX,2,,7.6E-006
MPDATA,ALPX,2,,9.2E-006
MPDATA,ALPX,2,,1.08E-005
MPDATA,ALPX,2,,1.24E-005
MPDATA,ALPX,2,,1.4E-005

!建模型
block,0, 0.305,0,1.05,0, 0.305
BLOCK,0.00635, 0.29865,0,1.05, 0.00635, 0.29865
vovlap,1,2

wpoff,0.00635
wprota,,,90
vsbw,all
wpoff,-0.00635
wprota,,,90
vsbw,all
wpoff,-0.2923
wprota,,,90
vsbw,all
wpoff,0.2923
wprota,,,90
vsbw,all
wpcsys,-1
vglue,all
vsel,s,loc,x,0.00635,0.2923
vsel,r,loc,z,0.00635,0.2923
vatt,2,,2
vsel,inve
vatt,1,,1
allsel
lesize,all,0.025
vmesh,all

!加载
/solu
Antype,TRANS
TIMINT,ON
AUTOTS,ON
OUTRES,,ALL
TIME,50
NSUBST,25
timint,on       
BFUNIF,TEMP,20
KBC,0
asel,s,,,1
asel,a,,,3,4,1
asel,a,,,18,19,1
asel,a,,,26
asel,a,,,31
asel,a,,,34,35
asel,a,,,38
asel,a,,,45,46
NSLA,,,1
D,ALL,TEMP,200
allsel
asel,s,,,1
asel,a,,,3,4,1
asel,a,,,18,19,1
asel,a,,,26
asel,a,,,31
asel,a,,,34,35
asel,a,,,38
asel,a,,,45,46
NSLA,,,1
SF,ALL,CONV,10
allsel

SOLVE
FINISH

/POST1
SET,LAST
/EDGE,,1
/plopts,info,on
/title,temperature of steel about CFST
Esel,s,mat,,1
Nsle,s
Plnsol,temp
Allsel
/title,temperature of concrete about CFST
Esel,s,mat,,2
Nsle,s
Plnsol,temp
Allsel
Plnsol,temp
Finish



/post26
/title, between temperature and time about steel of CFST
/title, between temperature and time about concrete of CFST
/plopts,info,on
/axlab,x,time,(min)
/axlab,y,temperature
/gthk,axis,3
/gthk,curve,3
/replot
/xrange,0,50
/yrange,20,200
!混凝土中间点
Nsol,2,5346,temp
plvar,2
!混凝土37mm
nsol,3,3075,temp
plvar,3
!钢管角点
nsol,4,9633,temp
plvar,4
!钢管中间角点
nsol,5,9655,temp
plvar,5
!钢管中间-中间
nsol,6,8762,temp
plvar,6

Nsol,3,,t
emp
0.,3
Finish
这是本人的做的热应力分析一块，转换为结构还没有做，这个里面很多问题，我做的这个不要求太高，只要得到大概曲线就可以。
麻烦各位帮忙看一下了！
我想用多荷载步求解
查了一晚上的资料
脑袋糊涂了，因为急用，拜托大家！

ymz.124

难道斑竹也不在吗？