帮忙看一下这个程序中
aatt,1,,1和aatt,2,,2这两句是不是位置反了???
程序如下:
/batch,list
/prep7, Silicon beam deflection from an applied voltage
et,1,121 ! PLANE121 element for air region
et,2,121 ! Temporary element for beam region
emunit,epzro,8.854e-6 ! Free-space permittivity, [$micro]MKSV units
mp,perx,2,1 ! Relative permittivity for air
bl=150 ! Beam length ([$micro]m)
bh=2.0 ! Beam height
glc=bl/2 ! Center location of ground
gl=90 ! Ground conductor length
gh=1.5 ! Ground conductor thickness
gap=4.5 ! Air gap
vltg=120 ! Applied voltage
rectng,0,bl,gap,bh+gap ! Create model
rectng,glc-gl/2,glc+gl/2,-gh,0
rectng,-10,170,-20,30
aovlap,all
asel,s,area,,1 ! Area for beam
aatt,1,,1
asel,s,area,,4 ! Area for air elements
cm,air,area ! Group air area into component
aatt,2,,2
allsel,all
smrtsiz,2
amesh,1 ! Mesh beam
mshape,1
amesh,4 ! Mesh air with triangles
asel,s,area,,1
lsla,s
dl,all,,volt,vltg ! Apply voltage to beam
asel,s,area,,2
lsla,s
dl,all,,volt,0 ! Ground conductor (not meshed)
allsel,all
et,1,0 ! Set structure to null element type
physics,write,ELECTROS ! Write electrostatic physics file
physics,clear ! Clear Physics
et,1,82,,,2 ! Define beam element y type
et,2,0 ! Set air to null element type
mp,ex,1,170e3 ! Set Modulus [$micro]N/([$micro]m)**2
mp,nuxy,1,0.34
dl,4,,ux,0 ! Apply beam constraints
dl,4,,uy,0
allsel,all
finish
physics,write,STRUCTURE ! Write structural physics file
ESSOLV,'ELECTROS','STRUCTURE',2,0,'air',,,,10 ! Solve coupled-field problem
finish
physics,read,ELECTROS ! Read electrostatic physics file
/post1
set,last ! Retrieve results
esel,s,mat,,2 ! Select air elements
etable,fx,fmag,x ! Retrieve electrostatic forces
etable,fy,fmag,y
ssum ! Sum forces
finish |