做一个瞬态热分析和静态结构分析耦合的例子,就是先按瞬态热分析计算出温度场分布,再转为结构分析,把温度作为载荷加上。分析的时候就会报错,如图所示。请大侠指点。帮忙看看前面的分析过程有问题没,万分感激
命令流如下:
!!!!!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
|