降温法模拟预应力

本人正在学习ansys，在做一个实体简支梁的分析，加载预应力时，按初应变方法没问题，但是同一个模型，我改用降温法缺不起作用，算也能算，但是结果都是0，应该是预应力没有加进去，找了好长时间，也没发现哪里的问题。现在把两种模型的命令流都贴出来，请各位大侠帮忙看看是什么原因
一、初应变方法
fini
/clear
/prep7
eg=1.95e11
ag=1.4e-4  
eh=3.45e10
!r0=9345
yyl=182000
et,1,link8
et,2,solid95
r,1,15*ag,yyl/eg/ag*1.036258
r,2
mp,ex,1,eg
mp,prxy,1,0.3
mp,ex,2,eh
mp,prxy,2,0.2
!------------定义实体
!产生关键点
k,1,-4.651,0
k,2,-5.251,0
k,3,-5.95,0
k,4,-7.95,0
k,5,-1.005,-0.22
k,6,-2.805,-0.22
k,7,-4.651,-0.22
k,8,-5.251,-0.42
k,9,-5.95,-0.4
k,10,-7.95,-0.2
k,11,-1.605,-0.42
k,12,-2.205,-0.42
k,13,-5.251,-0.6
k,14,-5.883,-0.6
k,15,-1.605,-0.95
k,16,-2.205,-0.95
k,17,-5.134,-0.95
k,18,-5.767,-0.95
k,19,-1.205,-1.15
k,20,-2.605,-1.15
k,21,-4.734,-1.15
k,22,-5.7,-1.15
k,23,-5.55,-1.6
k,24,-2.805,0       
k,25,-1.005,0       
k,26,-2.605,-1.6       
k,27,-1.205,-1.6       
k,28,-4.914,-1.6


k,3001,-4.914,0,6
k,3002,-5.514,0,6
k,3003,-5.95,0,6
k,3004,-7.95,0,6
k,3005,-1.13,-0.22,6
k,3006,-2.68,-0.22,6
k,3007,-4.914,-0.22,6
k,3008,-5.514,-0.42,6
k,3009,-5.95,-0.4,6
k,3010,-7.95,-0.2,6
k,3011,-1.73,-0.42,6
k,3012,-2.08,-0.42,6
k,3013,-5.514,-0.6,6
k,3014,-5.883,-0.6,6
k,3015,-1.73,-1.2,6
k,3016,-2.08,-1.2,6
k,3017,-5.314,-1.2,6
k,3018,-5.683,-1.2,6
k,3019,-1.33,-1.4,6
k,3020,-2.48,-1.4,6
k,3021,-4.914,-1.4,6
k,3022,-5.617,-1.4,6
k,3023,-5.55,-1.6,6
k,3024,-2.68,0,6
k,3025,-1.13,0,6
k,3026,-2.48,-1.6,6
k,3027,-1.33,-1.6,6
k,3028,-4.734,-1.6,6

!进行关键点的对称复制
ksym,x,1,28,1,100
ksym,x,3001,3028,1,100
kgen,3,1,28,1,0,0,0.5,1000
kgen,3,101,128,1,0,0,0.5,1000
kgen,2,1,28,1,0,0,29.96,5000
kgen,2,101,128,1,0,0,29.96,5000
kgen,3,5001,5028,1,0,0,0.5,1000
kgen,3,5101,5128,1,0,0,0.5,1000
kgen,2,3001,3028,1,0,0,18.96,1000
kgen,2,3101,3128,1,0,0,18.96,1000


              
*do,i,1000,6000,1000   !循环生成各个阶段的小体积
v,5+i-1000,i-1000+25,i-1000+125,i-1000+105,5+i,25+i,125+i,105+i
v,5+i-1000,i-1000+6,i-1000+24,i-1000+25,5+i,6+i,24+i,25+i
v,105+i-1000,i-1000+106,i-1000+124,i-1000+125,105+i,106+i,124+i,125+i
v,6+i-1000,i-1000+7,i-1000+1,i-1000+24,6+i,7+i,1+i,24+i
v,106+i-1000,i-1000+107,i-1000+101,i-1000+124,106+i,107+i,101+i,124+i
v,28+i-1000,i-1000+23,i-1000+22,i-1000+21,28+i,23+i,22+i,21+i
v,128+i-1000,i-1000+123,i-1000+122,i-1000+121,128+i,123+i,122+i,121+i
v,26+i-1000,i-1000+28,i-1000+21,i-1000+20,26+i,28+i,21+i,20+i
v,126+i-1000,i-1000+128,i-1000+121,i-1000+120,126+i,128+i,121+i,120+i
v,27+i-1000,i-1000+26,i-1000+20,i-1000+19,27+i,26+i,20+i,19+i
v,127+i-1000,i-1000+126,i-1000+120,i-1000+119,127+i,126+i,120+i,119+i

v,1+i-1000,i-1000+2,i-1000+8,i-1000+7,1+i,2+i,8+i,7+i
v,101+i-1000,i-1000+102,i-1000+108,i-1000+107,101+i,102+i,108+i,107+i
v,2+i-1000,i-1000+3,i-1000+9,i-1000+8,2+i,3+i,9+i,8+i
v,102+i-1000,i-1000+103,i-1000+109,i-1000+108,102+i,103+i,109+i,108+i
v,3+i-1000,i-1000+4,i-1000+10,i-1000+9,3+i,4+i,10+i,9+i
v,103+i-1000,i-1000+104,i-1000+110,i-1000+109,103+i,104+i,110+i,109+i
v,8+i-1000,i-1000+9,i-1000+14,i-1000+13,8+i,9+i,14+i,13+i
v,108+i-1000,i-1000+109,i-1000+114,i-1000+113,108+i,109+i,114+i,113+i
v,13+i-1000,i-1000+14,i-1000+18,i-1000+17,13+i,14+i,18+i,17+i
v,113+i-1000,i-1000+114,i-1000+118,i-1000+117,113+i,114+i,118+i,117+i
v,17+i-1000,i-1000+18,i-1000+22,i-1000+21,17+i,18+i,22+i,21+i
v,117+i-1000,i-1000+118,i-1000+122,i-1000+121,117+i,118+i,122+i,121+i
v,27+i-1000,i-1000+19,i-1000+119,i-1000+127,27+i,19+i,119+i,127+i
v,5+i-1000,i-1000+6,i-1000+12,i-1000+11,5+i,6+i,12+i,11+i
v,105+i-1000,i-1000+106,i-1000+112,i-1000+111,105+i,106+i,112+i,111+i
v,11+i-1000,i-1000+12,i-1000+16,i-1000+15,11+i,12+i,16+i,15+i
v,111+i-1000,i-1000+112,i-1000+116,i-1000+115,111+i,112+i,116+i,115+i
v,15+i-1000,i-1000+16,i-1000+20,i-1000+19,15+i,16+i,20+i,19+i
v,115+i-1000,i-1000+116,i-1000+120,i-1000+119,115+i,116+i,120+i,119+i
*enddo      
allsel
!vglue,all
/view,1,1,1,1
/ang,1
vplot

!------------定义力筋线
!!!!!!!!!!!!!!!!!!!!!腹板变弯束F1
numcmp,all                                          
*get,kpmax,kp,0,num,max
*get,lmax,line,0,num,max
lsel,none         
k,kpmax+1,1.905,-0.25,0.16
k,kpmax+2,1.905,-1.18,6.032
k,kpmax+3,1.905,-1.18,24.928
k,kpmax+4,1.905,-0.25,30.8
l,kpmax+1,kpmax+2
l,kpmax+2,kpmax+3              
L,kpmax+3,kpmax+4
lfillt,lmax+1,lmax+2,8
lfillt,lmax+2,lmax+3,8
nummrg,all
lcomb,all            
lsymm,x,all               
cm,F1,line                    
allsel
numcmp,all
!!!!!!!!!!!!!!!!!!!!!腹板变弯束F2
numcmp,all                                          
*get,kpmax,kp,0,num,max
*get,lmax,line,0,num,max
lsel,none         
k,kpmax+1,1.905,-0.6,0.16
k,kpmax+2,1.905,-1.33,5.354
k,kpmax+3,1.905,-1.33,25.606
k,kpmax+4,1.905,-0.6,30.8
l,kpmax+1,kpmax+2
l,kpmax+2,kpmax+3              
L,kpmax+3,kpmax+4
lfillt,lmax+1,lmax+2,8
lfillt,lmax+2,lmax+3,8
nummrg,all
lcomb,all            
lsymm,x,all               
cm,F2,line                    
allsel
numcmp,all

!-----------单元划分

cmsel,s,F1
cmsel,a,F2  
latt,1,1,1
esize,1
lmesh,all
vsel,all
vatt,2,2,2
esize,1
vmesh,all
allsel,all
!-------------定义约束
nsel,s,loc,y,-1.6
nsel,r,loc,x,-3.3,-2.7
nsel,r,loc,z,0.2,0.8
d,all,ux,,,,,uy

nsel,s,loc,y,-1.6
nsel,r,loc,x,2.7,3.3
nsel,r,loc,z,0.2,0.8
d,all,ux,,,,,uy

nsel,s,loc,y,-1.6
nsel,r,loc,x,-3.3,-2.7
nsel,r,loc,z,30.16,30.76
d,all,ux,,,,,uy,uz

nsel,s,loc,y,-1.6
nsel,r,loc,x,2.7,3.3
nsel,r,loc,z,30.16,30.76
d,all,ux,,,,,uy,uz
!耦合自由度
cmsel,s,F1
cmsel,a,F2
nsll,s,1
cm,cmljnod,node
*get,max1,node,0,count
*dim,ojd,,max1
*dim,jd,,max1
*get,nod1,node,0,num,min
ojd(1)=nod1
*do,i,2,max1
ojd(i)=ndnext(ojd(i-1))
*enddo
allsel,all
nsel,all
cmsel,u,cmljnod
*do,i,1,max1
nod1=ojd(i)
j=nnear(nod1)
jd(i)=j
*enddo
nsel,all
ji=1
*do,i,1,max1
cp,ji,ux,ojd(i),jd(i)
cp,ji+1,uy,ojd(i),jd(i)
cp,ji+2,uz,ojd(i),jd(i)
ji=ji+3
*enddo
allsel,all
finish
/solu
solve
finish

二、降温法模拟预应力
fini
/clear
/prep7
eg=1.95e11
ag=1.4e-4  
eh=3.45e10
!r0=9345
yyl=182000
et,1,link8
et,2,solid95
r,1,15*ag
r,2
mp,ex,1,eg
mp,prxy,1,0.3
mp,dens,1,7850
mp,ex,2,eh
mp,prxy,2,0.2
!------------定义实体
!产生关键点
k,1,-4.651,0
k,2,-5.251,0
k,3,-5.95,0
k,4,-7.95,0
k,5,-1.005,-0.22
k,6,-2.805,-0.22
k,7,-4.651,-0.22
k,8,-5.251,-0.42
k,9,-5.95,-0.4
k,10,-7.95,-0.2
k,11,-1.605,-0.42
k,12,-2.205,-0.42
k,13,-5.251,-0.6
k,14,-5.883,-0.6
k,15,-1.605,-0.95
k,16,-2.205,-0.95
k,17,-5.134,-0.95
k,18,-5.767,-0.95
k,19,-1.205,-1.15
k,20,-2.605,-1.15
k,21,-4.734,-1.15
k,22,-5.7,-1.15
k,23,-5.55,-1.6
k,24,-2.805,0       
k,25,-1.005,0       
k,26,-2.605,-1.6       
k,27,-1.205,-1.6       
k,28,-4.914,-1.6


k,3001,-4.914,0,6
k,3002,-5.514,0,6
k,3003,-5.95,0,6
k,3004,-7.95,0,6
k,3005,-1.13,-0.22,6
k,3006,-2.68,-0.22,6
k,3007,-4.914,-0.22,6
k,3008,-5.514,-0.42,6
k,3009,-5.95,-0.4,6
k,3010,-7.95,-0.2,6
k,3011,-1.73,-0.42,6
k,3012,-2.08,-0.42,6
k,3013,-5.514,-0.6,6
k,3014,-5.883,-0.6,6
k,3015,-1.73,-1.2,6
k,3016,-2.08,-1.2,6
k,3017,-5.314,-1.2,6
k,3018,-5.683,-1.2,6
k,3019,-1.33,-1.4,6
k,3020,-2.48,-1.4,6
k,3021,-4.914,-1.4,6
k,3022,-5.617,-1.4,6
k,3023,-5.55,-1.6,6
k,3024,-2.68,0,6
k,3025,-1.13,0,6
k,3026,-2.48,-1.6,6
k,3027,-1.33,-1.6,6
k,3028,-4.734,-1.6,6

!进行关键点的对称复制
ksym,x,1,28,1,100
ksym,x,3001,3028,1,100
kgen,3,1,28,1,0,0,0.5,1000
kgen,3,101,128,1,0,0,0.5,1000
kgen,2,1,28,1,0,0,29.96,5000
kgen,2,101,128,1,0,0,29.96,5000
kgen,3,5001,5028,1,0,0,0.5,1000
kgen,3,5101,5128,1,0,0,0.5,1000
kgen,2,3001,3028,1,0,0,18.96,1000
kgen,2,3101,3128,1,0,0,18.96,1000


              
*do,i,1000,6000,1000   !循环生成各个阶段的小体积
v,5+i-1000,i-1000+25,i-1000+125,i-1000+105,5+i,25+i,125+i,105+i
v,5+i-1000,i-1000+6,i-1000+24,i-1000+25,5+i,6+i,24+i,25+i
v,105+i-1000,i-1000+106,i-1000+124,i-1000+125,105+i,106+i,124+i,125+i
v,6+i-1000,i-1000+7,i-1000+1,i-1000+24,6+i,7+i,1+i,24+i
v,106+i-1000,i-1000+107,i-1000+101,i-1000+124,106+i,107+i,101+i,124+i
v,28+i-1000,i-1000+23,i-1000+22,i-1000+21,28+i,23+i,22+i,21+i
v,128+i-1000,i-1000+123,i-1000+122,i-1000+121,128+i,123+i,122+i,121+i
v,26+i-1000,i-1000+28,i-1000+21,i-1000+20,26+i,28+i,21+i,20+i
v,126+i-1000,i-1000+128,i-1000+121,i-1000+120,126+i,128+i,121+i,120+i
v,27+i-1000,i-1000+26,i-1000+20,i-1000+19,27+i,26+i,20+i,19+i
v,127+i-1000,i-1000+126,i-1000+120,i-1000+119,127+i,126+i,120+i,119+i

v,1+i-1000,i-1000+2,i-1000+8,i-1000+7,1+i,2+i,8+i,7+i
v,101+i-1000,i-1000+102,i-1000+108,i-1000+107,101+i,102+i,108+i,107+i
v,2+i-1000,i-1000+3,i-1000+9,i-1000+8,2+i,3+i,9+i,8+i
v,102+i-1000,i-1000+103,i-1000+109,i-1000+108,102+i,103+i,109+i,108+i
v,3+i-1000,i-1000+4,i-1000+10,i-1000+9,3+i,4+i,10+i,9+i
v,103+i-1000,i-1000+104,i-1000+110,i-1000+109,103+i,104+i,110+i,109+i
v,8+i-1000,i-1000+9,i-1000+14,i-1000+13,8+i,9+i,14+i,13+i
v,108+i-1000,i-1000+109,i-1000+114,i-1000+113,108+i,109+i,114+i,113+i
v,13+i-1000,i-1000+14,i-1000+18,i-1000+17,13+i,14+i,18+i,17+i
v,113+i-1000,i-1000+114,i-1000+118,i-1000+117,113+i,114+i,118+i,117+i
v,17+i-1000,i-1000+18,i-1000+22,i-1000+21,17+i,18+i,22+i,21+i
v,117+i-1000,i-1000+118,i-1000+122,i-1000+121,117+i,118+i,122+i,121+i
v,27+i-1000,i-1000+19,i-1000+119,i-1000+127,27+i,19+i,119+i,127+i
v,5+i-1000,i-1000+6,i-1000+12,i-1000+11,5+i,6+i,12+i,11+i
v,105+i-1000,i-1000+106,i-1000+112,i-1000+111,105+i,106+i,112+i,111+i
v,11+i-1000,i-1000+12,i-1000+16,i-1000+15,11+i,12+i,16+i,15+i
v,111+i-1000,i-1000+112,i-1000+116,i-1000+115,111+i,112+i,116+i,115+i
v,15+i-1000,i-1000+16,i-1000+20,i-1000+19,15+i,16+i,20+i,19+i
v,115+i-1000,i-1000+116,i-1000+120,i-1000+119,115+i,116+i,120+i,119+i
*enddo

         
allsel
!vglue,all


/view,1,1,1,1
/ang,1
vplot

!------------定义力筋线
!!!!!!!!!!!!!!!!!!!!!腹板变弯束F1
numcmp,all                                          
*get,kpmax,kp,0,num,max
*get,lmax,line,0,num,max
lsel,none         
k,kpmax+1,1.905,-0.25,0.16
k,kpmax+2,1.905,-1.18,6.032
k,kpmax+3,1.905,-1.18,24.928
k,kpmax+4,1.905,-0.25,30.8
l,kpmax+1,kpmax+2
l,kpmax+2,kpmax+3              
L,kpmax+3,kpmax+4
lfillt,lmax+1,lmax+2,8
lfillt,lmax+2,lmax+3,8
nummrg,all
lcomb,all            
lsymm,x,all               
cm,F1,line                    
allsel
numcmp,all
!!!!!!!!!!!!!!!!!!!!!腹板变弯束F2
numcmp,all                                          
*get,kpmax,kp,0,num,max
*get,lmax,line,0,num,max
lsel,none         
k,kpmax+1,1.905,-0.6,0.16
k,kpmax+2,1.905,-1.33,5.354
k,kpmax+3,1.905,-1.33,25.606
k,kpmax+4,1.905,-0.6,30.8
l,kpmax+1,kpmax+2
l,kpmax+2,kpmax+3              
L,kpmax+3,kpmax+4
lfillt,lmax+1,lmax+2,8
lfillt,lmax+2,lmax+3,8
nummrg,all
lcomb,all            
lsymm,x,all               
cm,F2,line                    
allsel
numcmp,all

!-----------单元划分

cmsel,s,F1
cmsel,a,F2  
latt,1,1,1
esize,1
lmesh,all
vsel,all
vatt,2,2,2
esize,1
vmesh,all
allsel,all
!-------------定义约束
nsel,s,loc,y,-1.6
nsel,r,loc,x,-3.3,-2.7
nsel,r,loc,z,0.2,0.8
d,all,ux,,,,,uy

nsel,s,loc,y,-1.6
nsel,r,loc,x,2.7,3.3
nsel,r,loc,z,0.2,0.8
d,all,ux,,,,,uy

nsel,s,loc,y,-1.6
nsel,r,loc,x,-3.3,-2.7
nsel,r,loc,z,30.16,30.76
d,all,ux,,,,,uy,uz

nsel,s,loc,y,-1.6
nsel,r,loc,x,2.7,3.3
nsel,r,loc,z,30.16,30.76
d,all,ux,,,,,uy,uz
!耦合自由度
cmsel,s,F1
cmsel,a,F2
nsll,s,1
cm,cmljnod,node
*get,max1,node,0,count
*dim,ojd,,max1
*dim,jd,,max1
*get,nod1,node,0,num,min
ojd(1)=nod1
*do,i,2,max1
ojd(i)=ndnext(ojd(i-1))
*enddo
allsel,all
nsel,all
cmsel,u,cmljnod
*do,i,1,max1
nod1=ojd(i)
j=nnear(nod1)
jd(i)=j
*enddo
nsel,all
ji=1
*do,i,1,max1
cp,ji,ux,ojd(i),jd(i)
cp,ji+1,uy,ojd(i),jd(i)
cp,ji+2,uz,ojd(i),jd(i)
ji=ji+3
*enddo
allsel,all
finish
/solu
E=1.95e5       !MPa
alpx=1.263e-5   !线膨胀系数
cmsel,s,F1
cmsel,a,F2
bfl,all,temp,-1300/(E*alpx)
allsel
solve
finish