为什么只能加载一面热流密度，另一面加载不上？？？用的高斯热源，高手支招！！！
命令流如下/clear
/units,si
/config,nres,10000
/filname,welding,1
/prep7
et,1,solid90
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,250
MPTEMP,3,500
MPTEMP,4,750
MPTEMP,5,1000
MPTEMP,6,1500
MPTEMP,7,1700
MPTEMP,8,2750
MPDATA,DENS,1,,7800
MPDATA,DENS,1,,7700
MPDATA,DENS,1,,7610
MPDATA,DENS,1,,7550
MPDATA,DENS,1,,7490
MPDATA,DENS,1,,7350
MPDATA,DENS,1,,7300
MPDATA,DENS,1,,7090
MPTEMP,,,,,,,,,
MPTEMP,1,20
MPTEMP,2,250
MPTEMP,3,500
MPTEMP,4,750
MPTEMP,5,1000
MPTEMP,6,1500
MPTEMP,7,1600
MPTEMP,8,1700
MPTEMP,9,2500
MPDATA,KXX,1,,50
MPDATA,KXX,1,,47
MPDATA,KXX,1,,40
MPDATA,KXX,1,,27
MPDATA,KXX,1,,30
MPDATA,KXX,1,,35
MPDATA,KXX,1,,200
MPDATA,KXX,1,,500
MPDATA,KXX,1,,1500
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,250
MPTEMP,3,500
MPTEMP,4,750
MPTEMP,5,1000
MPTEMP,6,1500
MPTEMP,7,1700
MPTEMP,8,2750
MPDATA,C,1,,460
MPDATA,C,1,,480
MPDATA,C,1,,530
MPDATA,C,1,,675
MPDATA,C,1,,670
MPDATA,C,1,,660
MPDATA,C,1,,2000
MPDATA,C,1,,2500
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,600
MPTEMP,8,1200
MPTEMP,9,2500
MPDATA,EX,1,,208000
MPDATA,EX,1,,204000
MPDATA,EX,1,,199000
MPDATA,EX,1,,191000
MPDATA,EX,1,,183000
MPDATA,EX,1,,174000
MPDATA,EX,1,,164000
MPDATA,EX,1,,80000
MPDATA,EX,1,,80000
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,600
MPTEMP,8,1500
MPTEMP,9,2500
UIMP,1,REFT,,,25
MPDATA,ALPX,1,,1.37e-5
MPDATA,ALPX,1,,1.37e-5
MPDATA,ALPX,1,,1.38e-5
MPDATA,ALPX,1,,1.41e-5
MPDATA,ALPX,1,,1.44e-5
MPDATA,ALPX,1,,1.48e-5
MPDATA,ALPX,1,,1.5e-5
MPDATA,ALPX,1,,1.5e-5
MPDATA,ALPX,1,,1.5e-5
TB,BISO,1,6,2,
TBTEMP,20
TBDATA,,370,732,,,,
TBTEMP,250
TBDATA,,289,782,,,,
TBTEMP,500
TBDATA,,230,315,,,,
TBTEMP,650
TBDATA,,30,50,,,,
TBTEMP,1200
TBDATA,,3,10,,,,
TBTEMP,2500
TBDATA,,3,10,,,,
TBDE,BISO,1,,,
TB,BISO,1,6,2,
TBTEMP,20
TBDATA,,370,732,,,,
TBTEMP,250
TBDATA,,289,782,,,,
TBTEMP,500
TBDATA,,230,315,,,,
TBTEMP,650
TBDATA,,30,50,,,,
TBTEMP,1200
TBDATA,,3,10,,,,
TBTEMP,2500
TBDATA,,3,10,,,,
k,1,0,0,0
k,2,0,0,0.12
k,3,0,0.006,0.12
k,4,0,0.006,0
k,5,0,0,0.14
k,6,0,0.006,0.14
k,7,0,0.006,0.146
k,8,0,0,0.146
k,9,0,0,0.154
k,10,0,0.006,0.154
k,11,0,0.006,0.16
k,12,0,0,0.16
k,13,0,0,0.18
k,14,0,0.006,0.18
k,15,0,0.006,0.3
k,16,0,0,0.3
k,17,0,0.012,0.154
k,18,0,0.012,0.146
k,19,0,0.032,0.146
k,20,0,0.032,0.154
k,21,0,0.156,0.154
k,22,0,0.156,0.146
a,1,2,3,4
a,2,5,6,3
a,5,6,7,8
a,7,8,9,10
a,9,10,11,12
a,11,12,13,14
a,13,14,15,16
a,6,7,18
a,7,10,17,18
a,10,11,17
a,17,18,19,20
a,19,20,21,22
aglue,all
/view,1,1,1,1
/replot
vext,all,,,0.3,0,0
vglue,all
lesize,70,,,50,,,,,1
lesize,67,,,20,,,,,1
mat,1
type,1
vmesh,7
lesize,40,,,50,,,,,1
lesize,36,,,20,,,,,1
mat,1
type,1
vmesh,1
lesize,64,,,5,,,,,1
mat,1
type,1
vmesh,6
lesize,44,,,5,,,,,1
mat,1
type,1
vmesh,2
lesize,57,,,2,,,,,1
mat,1
type,1
vmesh,5
lesize,54,,,2,,,,,1
mat,1
type,1
vmesh,4
lesize,47,,,2,,,,,1
mat,1
type,1
vmesh,3
vmesh,9
lesize,81,,,5,,,,,1
mat,1
type,1
vmesh,11
lesize,84,,,31,,,,,1
mat,1
type,1
vmesh,12
lesize,73,,,2,,,,,1
vmesh,8
lesize,78,,,2,,,,,1
vmesh,10
nummrg,all
numcmp,all
/image,save,vmesh,jpeg
save
rb=0.006
zb=rb/6
length=0.3
width=0.005
/solu
antype,4
tunif,25
TRNOPT,FULL
OUTRES,NSOL,all
NROPT,FULL,,ON
AUTOTS,1
KBC,1
lnsrch,1
TSRES,ERASE
TINTP,,,,1,,,,,
TIMINT,0,struct
TIMINT,1,THERM
TIMINT,0,MAG
allsel,all
nsel,s,ext
sf,all,conv,15,25
v=0.006
tinc=rb/v
tm=length/v
U=15
I=100
ef=0.7
p=U*I*ef
r=0.005
qmax=3*p/(3.14159*r**2)
max_time=length/rb+1
max_x=length/rb+1
max_y=width/zb+1
local,11,0,0,0.009,0.143,,45,
*dim,flux2,table,max_x,max_y,max_time,x,y,time
*do,k,1,max_time,1
*do,i,1,max_x,1
flux2(i,0,k)=(i-1)*rb
*enddo
*do,j,1,max_y,1
flux2(0,j,k)=(j-1)*zb
*enddo
*do,i,1,max_x,1
*do,j,1,max_y,1
xcenter=v*(k-1)*tinc
ycenter=0
distance=sqrt(((i-1)*rb-xcenter)**2+((j-1)*zb-ycenter)**2)
*if,distance,le,r,then
flux2(i,j,k)=qmax/exp(3*distance**2/(r**2))
*else
flux2(i,j,k)=0
*endif
*enddo
*enddo
flux2(0,0,k)=(k-1)*tinc
*enddo
t=0
*do,i,1,5,1
/solu
t=t+1
time,t
NSUBST,2,,,1
allsel,all
asel,s,area,,44
nsla,s,1
sfdele,all,conv
sf,all,hflux,%flux2%
allsel,all
solve
asel,s,area,,44
nsla,s,1
sfdele,all,hflux
sf,all,conv,15,25
allsel,all
*enddo
*DO,I,1,20
/solu
allsel,all
t=t+1
TIME,t
NSUBST,2,,,1
SOLVE
*ENDDO
*DO,I,1,20
/solu
allsel,all
t=t+1*10
TIME,t
NSUBST,2,,,1
SOLVE
*ENDDO
save
/solu
antype,4
tunif,25
TRNOPT,FULL
OUTRES,NSOL,all
NROPT,FULL,,ON
AUTOTS,1
KBC,1
lnsrch,1
TSRES,ERASE
TINTP,,,,1,,,,,
TIMINT,0,struct
TIMINT,1,THERM
TIMINT,0,MAG
rb=0.006
zb=rb/6
length=0.3
width=0.005
v=0.006
tinc=rb/v
tm=length/v
U=15
I=100
ef=0.7
p=U*I*ef
r=0.005
qmax=3*p/(3.14159*r**2)
max_time=length/rb+1
max_x=length/rb+1
max_y=width/zb+1
csdele,11,,1
local,11,0,0,0.009,0.157,,135,
*dim,flux3,table,max_x,max_y,max_time,x,y,time
*do,k,1,max_time,1
*do,i,1,max_x,1
flux3(i,0,k)=(i-1)*rb
*enddo
*do,j,1,max_y,1
flux3(0,j,k)=(j-1)*zb
*enddo
*do,i,1,max_x,-1
*do,j,1,max_y,-1
xcenter=v*(k-1)*tinc
ycenter=0
distance=sqrt(((i-1)*rb-xcenter)**2+((j-1)*zb-ycenter)**2)
*if,distance,le,r,then
flux3(i,j,k)=qmax/exp(3*distance**2/(r**2))
*else
flux3(i,j,k)=0
*endif
*enddo
*enddo
flux3(0,0,k)=(k-1)*tinc
*enddo
*do,i,1,5,1
/solu
t=t+1
time,t
NSUBST,2,,,1
allsel,all
asel,s,area,,49
nsla,s,1
sfdele,all,conv
sf,all,hflux,%flux2%
allsel,all
solve
asel,s,area,,49
nsla,s,1
sfdele,all,hflux
sf,all,conv,15,25
allsel,all
*enddo
其中，44面可以加载上，我为了调程序就只取到5秒，另一面总是加载在末端的一个点，摸索了很久，不知道为什么？请高手指点
我没能放图片在上面
页:
[1]