/clear
/units,si
/prep7
et,1,plane55
et,2,surf151
keyopt,2,8,1
mp,dens,1,7850
mp,kxx,1,71
mp,c,1,460.24
mp,alpx,1,11.3e-6
mp,ex,1,215e9
mp,nuxy,1,0.26
mp,kxx,2,87.86
mp,c,2,510
mp,alpx,2,13.36e-6
mp,ex,2,204e9
mp,nuxy,2,0.28
mptemp,2,25,355,460,565,854,1024
mpdata,enth,2,1,8.48e8,1.45e9,1.64e9,1.81e9,2.43e9,4.28e9
p=1600
pi=3.1415926
r=0.003
kk=0.478
Qmax=p*kk/(pi*r**2)
lx=0.06
ly=0.02
v=0.004
lsize=0.0005
tm=lx/v
tinc=lsize/v
max_time=1+tm/tinc
max_x=1+lx/lsize
max_y=1+ly/lsize
rectng,0,lx,0,0.001
rectng,0,lx,0.001,ly
aglue,all
numcmp,all
esize,lsize
mshape,0,2D
mshkey,1
type,1
mat,2
amesh,1
esize,lsize*4
type,1
mat,1
amesh,2
esize,lsize
type,2
lsel,s,line,,1,5
lmesh,all
allsel,all
eplot
finish
*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)*lsize
*enddo
*do,j,1,max_y,1
flux2(0,j,k)=(j-1)*lsize
*enddo
*do,i,1,max_x,1
*do,j,1,max_y,1
xcenter=v*(k-1)*tinc
ycenter=0
distance=sqrt(((i-1)*lsize-xcenter)**2+((j-1)*lsize-ycenter)**2)
*if,distance,le,r,then
flux2(i,j,k)=Qmax*exp(-distance**2/(r**2))
*else
flux2(i,j,k)=0
*endif
*enddo
*enddo
flux2(0,0,k)=(k-1)*tinc
*enddo
/solu
antype,trans
tunif,100
kbc,1
timint,on
autots,off
esel,s,type,,,2
nsel,s
sf,all,conv,-1,20
allsel
*do,t,1e-6,tm,tinc
time,t
autots,on
deltim,,,,1
kbc,0
nsel,s,loc,x,0,0.06
nsel,r,loc,y,0,0.02
sf,all,hflux,%flux2%
outres,all,all
solve
*enddo
allsel,all
finish
为什么怎么加上表面效应载荷后仍没有??
敬请各位大侠指点 |