材料属性与温度的关系影响求解？热应力分析遇见的问题

alwayszi

做一个瞬态热分析和静态结构分析耦合的例子，就是先按瞬态热分析计算出温度场分布，再转为结构分析，把温度作为载荷加上。分析的时候就会报错，如图所示。请大侠指点。帮忙看看前面的分析过程有问题没，万分感激
命令流如下：
!!!!!1, fr,2,copper,3,plcc,4,solder!!!!!!!!
finish
/clear
/file, pcb thermalmechanical analysis
/title, pcb thermalmechanical analysis
/prep7
lpcb=200e-3
wpcb=180e-3
uh=0.2e-3
wplcc=30e-3
thplcc=1e-3
thsd=0.4e-3
l1=80e-3
l2=40e-3
w1=40e-3
hair=20
tair=25
hng=30
t1=230
t2=180
t3=260
t4=25
et,1,70
mp,dens,1,1870                                               
mp,kxx,1,0.25
mptemp,1,25,70,100,130,220
mpdata,c,1,1,1200,1380,1452,1512.5,1625
mp,dens,2,8933                                                  
mp,kxx,2,394
mp,c,2,390
mp,dens,3,1800                                                
mp,kxx,3,18
mptemp,1,25,80,160,170,220
mpdata,c,3,1,740,850,880,960,970
mp,dens,4,8350
mp,kxx,4,50.9
mp,c,4,150

   
block,0,lpcb,0,wpcb,0,0.5*uh                                       
block,0,lpcb,0,wpcb,0.5*uh,1.5*uh
block,0,lpcb,0,wpcb,1.5*uh,2.5*uh
block,0,lpcb,0,wpcb,2.5*uh,3.5*uh
block,0,lpcb,0,wpcb,3.5*uh,4.5*uh
block,0,lpcb,0,wpcb,4.5*uh,5.5*uh
block,0,lpcb,0,wpcb,5.5*uh,6.5*uh
block,0,lpcb,0,wpcb,6.5*uh,7*uh
/pnum,volu,1
vplot

block,l1,l1+wplcc,w1,w1+wplcc,0,7*uh+thsd
block,l2,l2+wplcc,wpcb-w1-wplcc,wpcb-w1,0,7*uh+thsd
block,lpcb-l2-wplcc,lpcb-l2,wpcb-w1-wplcc,wpcb-w1,0,7*uh+thsd

block,l1,l1+wplcc,w1,w1+wplcc,7*uh+thsd,7*uh+thsd+thplcc
block,l2,l2+wplcc,wpcb-w1-wplcc,wpcb-w1,7*uh+thsd,7*uh+thsd+thplcc
block,lpcb-l2-wplcc,lpcb-l2,wpcb-w1-wplcc,wpcb-w1,7*uh+thsd,7*uh+thsd+thplcc
vsel,all
vovlap,all
vglue,all
   
vsel,all
wprot,0,90
wpoffs,0,0,-w1
vsbw,all
wpave
wpoffs,0,0,-w1-wplcc
vsbw,all
wpave
wpoffs,0,0,-(wpcb-w1-wplcc)
vsbw,all
wpave
wpoffs,0,0,-(wpcb-w1)
vsbw,all
wpave
wpcsys,0

wprot,0,0,90
wpoffs,0,0,l2
vsbw,all
wpoffs,0,0,wplcc
vsbw,all
wpoffs,0,0,l1-l2-wplcc
vsbw,all
wpoffs,0,0,wplcc
vsbw,all
wpave
wpoffs,0,0,lpcb-l2
vsbw,all
wpoffs,0,0,-wplcc
vsbw,all
wpave
wpcsys
vsel,s,loc,z,0,0.5*uh
vsel,a,loc,z,1.5*uh,2.5*uh
vsel,a,loc,z,4.5*uh,5.5*uh
vsel,a,loc,z,6.5*uh,7*uh
vatt,2,,1
allsel
vsel,s,loc,z,0.5*uh,1.5*uh
vsel,a,loc,z,2.5*uh,4.5*uh
vsel,a,loc,z,5.5*uh,6.5*uh
vatt,1,,1
allsel
vsel,s,loc,z,7*uh,7*uh+thsd
vatt,4,,1
allsel
vsel,s,loc,z,7*uh+thsd,7*uh+thsd+thplcc
vatt,3,,1
allsel
shpp,modify,1,60
mshape,0,3d
mshkey,1
esize,5e-3
vmesh,all

allsel
/title,model mesh
save
/solu
antype,trans
solcontrol,on
tref,25
tunif,25
ic,all,temp,25
kbc,1
tim=0   
*do,i,1,10
tim=tim+6
time,tim
deltim,2               
nsel,s,loc,z,0
nsel,a,loc,z,7*uh
nsel,a,loc,z,7*uh+thsd+thplcc
nsel,r,loc,x,0,0.020*i
sf,all,conv,hng,t1
allsel
nsel,s,loc,z,0
nsel,a,loc,z,7*uh
nsel,a,loc,z,7*uh+thsd+thplcc
nsel,r,loc,x,0.020*i,0.200
sf,all,conv,hair,tair
allsel
solve
sfdele,all,conv                 
*enddo

/post1
set,first
set,last
plnsol,temp
finish
/prep7
et,1,185

mp,nuxy,1,0.28
mptemp,1,25,70,100,130,220
mpdata,ctex,1,1,13.79e-6,18.27e-6,24.95e-6,31.62e-6,37.6e-6
mpdata,ex,1,1,6.85e9,6.72e9,6.42e9,5.31e9,3.42e9
mp,nuxy,2,0.3
mp,ex,2,80e6
mptemp,1,25,77,127,177,220
mpdata,ctex,2,1,16.65e-6,16.70e-6,17.12e-6,17.51e-6,17.85e-6

mp,nuxy,3,0.3
mp,ex,3,6e9
mp,ctex,3,15e-6
mptemp,1,25,50,75,100,125
mpdata,ex,4,1,29.8e9,27.6e9,25.4e9,23.2e9,21.0e9
mpdata,nuxy,4,1,0.383,0.385,0.390,0.397,0.400
mpdata,ctex,4,1,25.6e-6,25.8e-6,26.9e-6,27.3e-6,28.7e-6
nsel,s,loc,z,0
nsel,r,loc,x,0
nsel,r,loc,y,0
nplot
d,all,all,0
allsel
nsel,s,loc,z,0
nsel,r,loc,x,0
nsel,r,loc,y,wpcb
d,all,uz,0
allsel
nsel,s,loc,z,0
nsel,r,loc,y,0
nsel,r,loc,x,lpcb
d,all,uz,0
allsel
nsel,s,loc,z,0
nsel,r,loc,x,lpcb
nsel,r,loc,y,wpcb
nplot
d,all,uz,0
allsel
eplot
finish
/solu
!!!!!!!!!!!!!!!!!!!!!!
nlgeom,on                        
antype,static        
solcontrol,on
kbc,0
time,10                       
deltim,2
!!!!!!!!!!!!!!!!!!
nlgeom,on                        
antype,static        
solcontrol,on
kbc,0
tref,25
LDREAD,TEMP,,,60,2,'pcbthermalmechanicalanalysis','rth',' '
allsel
solve

alwayszi

没人遇到过这种问题么？

yuxuerui2008

alwayszi 发表于 2012-2-24 10:19
没人遇到过这种问题么？

不用一个结构转换命令吗？  比如ETCHG,TTS

alwayszi

yuxuerui2008 发表于 2012-2-25 09:33
不用一个结构转换命令吗？  比如ETCHG,TTS

找到原因在哪里了，nlgeom命令应该写在/prep7 里面，谢谢你

buaacmd

alwayszi 发表于 2012-2-28 10:48
找到原因在哪里了，nlgeom命令应该写在/prep7 里面，谢谢你

我做热结构耦合分析时， nlgeom命令写在/solu 后面的，没有问题 。
大变形分析选项，我觉得也应该写在求解命令流里面呀

alwayszi

buaacmd 发表于 2012-4-23 09:43
我做热结构耦合分析时， nlgeom命令写在/solu 后面的，没有问题 。
大变形分析选项，我觉得也应该写在求 ...

help里说明过,写在/prep7里才有效