求助：钢管自应力钢管砼应力分析

shiwei_zy

问题描述：
作钢管自应力钢管砼分析，对钢管由常温(30度)，加到60度，浇筑砼，钢管自然降温到30度，分析砼产生的应力第1、2、3主应力
Ansys建模：
选取1/4模型，砼用solid65单元，钢管用Solid45单元，加温度荷载用两种方式，一种用Bfv命令，直接让钢管降温30度分析结果,一种是采用两个载荷步，先将砼的Solid65单元杀死，给钢管的solid45加温度载荷，然后重新将Solid65单元激活，再计算，然后分析砼的应力，但发现两种计
算结果差异较大，显然第一种方法符合实际情况，不知问题出在哪里，请各位指教；还有本人刚学习Ansys,分析方法不知是得当，还请各位指教。
/prep7                  
/title,Anlyse of CFST   
et,1,solid65            
et,2,solid45            
et,3,shell63            
mp,ex,1,3.45e10         
mp,prxy,1,0.2           
mp,alpx,1,1.0e-5        
mp,dens,1,2400         
mp,ex,2,2.06e11         
mp,prxy,2,0.3           
mp,alpx,2,1.2e-5        
mp,dens,2,7850         
!建立物理模型           
k , 1  ,0     ,1        
k , 2  ,0.5   ,0.96     
k , 3  ,1     ,0.84     
k , 4  ,1.5   ,0.64     
k , 5  ,2     ,0.36     
k , 6  ,2.5   ,0        
bspline,1,2,3,4,5,6     
kdele,2,5,1            
numcmp,all              
wpoff,,1,               
wprota,,,-90            
cswpla,11,0            
cyl4,,,0.016,90,,180   
cyl4,,,0.096,90,0.016,180
cyl4,,,0.10,90,0.096,180
nummrg,all
numcmp,all
lesize,6,,,5,0.4
lesize,7,,,5,0.4
lesize,2,,,6
lesize,5,,,6
type,3
amesh,2
lesize,3,,,6
lesize,4,,,6
amesh,1
lesize,9,,,1
lesize,10,,,1
lesize,8,,,6
amesh,3
arsym,y,1,3,1,,0,0
nummrg,all
numcmp,all
wpcsys,-1,1
csys,0
type,1
mat,1
extopt,esize,90,0,
extopt,aclear,1
vdrag,1,2,4,5,,,1
type,2
mat,2
vdrag,3,6,,,,,1
nummrg,all
numcmp,all
finish
/solu
da,10,all
da,14,all
da,24,all
da,17,all
da,20,all
da,27,all
da,25,symm
da,18,symm
da,15,symm
da,9,symm
da,13,symm
da,23,symm
da,1,symm
da,2,symm
da,3,symm
da,4,symm
da,5,symm
da,6,symm
nlgeom,on                     !打开大变形开关      
nropt,full                    !完全牛顿-拉普森迭代
time,30
nsubst,100
outres,all
!加温度载荷方法1
vsel,s,,,5,6
bfv,all,temp,-30
vsel,all
esel,all
nsel,all
solve
!加温度载荷方法2
esel,s,type,,1
ekill,all
allsel,all
tref,60
tunif,30
esel,all
nsel,all
solve
esel,s,type,,1
ealive,all
esel,all
nsel,all
solve