梳齿谐振器机电耦合分析请教

albert-guo

请哪位大牛帮忙修改一下我的程序，万分感谢/batch,list
/filname,h10v100,1
/prep7

!==================parameter===========
dcv=0               !zhi liu dian ya
acv=100            !jiao liu dian ya

cbh=60   !===changable                     
clh=50   !===changable     
tbh=100   !===changable
tl=30    !===changable                    !centre block width and length

bb=6
bg=16     !===changable
bl=150
bn=21      !===changable only singular   
                         !beams dimension and number

jbl=30  !===changable
jbbh=(bn-1)*(bb+bg)/2+bb/2  !movable comb joint block

guo=500    !===changable
sl=4      !===changable
jian=10    !===changable
jslh=clh+sl
lai=20   
sg=(clh-sl)*2/3
lailh=jslh-sl-sg                !sustained beams

sjbl=jbl
sjbbh=(bn+1)*(bb+bg)/2-bg/2  !static comb joint block

cg=50   !===changable                    !comb gap
h=10    !===changable         !thickness
bbg=2   !===changable         !basic board gap
sh=h+bbg                      !static comb thinkness

!==================modeling==========

blc4,0,0,cbh,clh              !create 1/4 model, centre area
blc4,0,clh,tbh,tl             !step area
blc4,0,clh+tl,jbbh,jbl        !joint area

blc4,0,clh+tl+jbl+cg+bl,sjbbh,sjbl
blc4,0,clh+tl+jbl,bb/2,bl

*do,i,1,(bn-1)/2
blc4,bg*i+bb*(i-1)+bb/2,clh+tl+jbl,bb,bl
*enddo                        !movable beams area

*do,j,1,(bn+1)/2
blc4,sjbbh-(bg+bb)*(j-1),clh+tl+jbl+cg,-bb,bl
*enddo                       !static beams area

blc4,tbh,clh,guo,sl
blc4,tbh+guo,clh+sl,jian,-jslh
!blc4,tbh,0,-lai,lailh        !static block
blc4,tbh,lailh,guo,-sl        !sustained beam area

asel,s,loc,y,clh+tl+jbl+cg+bl/2,clh+tl+jbl+cg+bl+sjbl
aadd,all
!asel,s,area,,1,20
asel,s,loc,y,0,clh+tl+jbl+bl/2
!asel,u,loc,y,0,lai/2
asel,u,loc,x,tbh+guo/2,tbh+guo+jian
asel,a,loc,y,clh+tl+jbl,clh+tl+jbl+bl/2
aadd,all

asel,s,loc,y,0,lailh
asel,a,loc,x,tbh+guo/2,tbh+guo+jian
!asel,s,area,,21,23
aadd,all

asel,s,,,all
arsym,x,all                   !reflect
arsym,y,all

asel,s,loc,y,clh+tl+jbl+cg+bl/2,clh+tl+jbl+cg+bl+sjbl
aadd,all
asel,s,loc,y,-(clh+tl+jbl+cg+bl/2),-(clh+tl+jbl+cg+bl+sjbl)
aadd,all

asel,s,loc,x,-sjbbh/2,sjbbh/2
asel,u,loc,y,clh+tl+jbl+cg+bl/2,clh+tl+jbl+cg+bl+sjbl
asel,u,loc,y,-(clh+tl+jbl+cg+bl/2),-(clh+tl+jbl+cg+bl+sjbl)
aadd,all

asel,s,loc,x,tbh+guo/2,tbh+guo+jian      !sustain beam
aadd,all
asel,s,loc,x,-(tbh+guo/2),-(tbh+guo+jian)
aadd,all

allsel,all
vext,1,,,,,h
vext,3,,,,,h
vext,5,,,,,sh
vext,10,,,,,h
vext,12,,,,,sh

vglue,1,2,4

!============air model===============
block,-sjbbh-bb,sjbbh+bb,-(clh+tl+jbl+cg+bl+sjbl),clh+tl+jbl+cg+bl+sjbl,-10,sh

vsel,s,,,3,6
vovlap,all
allsel,all
vglue,3,5,6,7
finish


!==================meshing
/prep7

esize1=20
esize2=30
esize3=50

et,1,122
et,2,98
et,3,45

emunit,epzro,8.854e-6
mp,ex,1,170e3
mp,nuxy,1,0.34
mp,perx,2,1
MP,DENS,1,2.329e-15

VSEL,S, , ,7
vatt,2,,1
!smrtsize,10
esize,esize3
!aesize,all,10
mshape,0

VSEL,S, , ,3
vatt,1,,2
esize,esize1
mshape,1

VSEL,S, , ,5
vatt,1,,2
esize,esize1
mshape,1

VSEL,S, , ,6
vatt,1,,2
esize,esize3
mshape,1

VSEL,S, , ,1,2
vatt,1,,3
esize,esize2
!mshape,0

ALLSEL,ALL
!mshape,1
!smrtsize,10,5,2     !mesh size
vmesh,all

!====================load and physics
/prep7

VSEL,S, , ,7
cm,air,volu

!VSEL,s,,,3
!VSEL,A, , ,6
!VSEL,A, , ,5
!VSEL,A, , ,1
!VSEL,A, , ,2
ASLV,S
cm,struc,area   !define nonmorphing area

VSEL,S, , ,3
ASLV,S
da,all,volt,acv       ! initial volt
  
VSEL,S, , ,5
ASLV,S
da,all,volt,-acv
   
VSEL,S, , ,6
ASLV,S  
da,all,volt,0

et,2,0  
et,3,0  
ALLSEL,ALL
physics,write,ELECTROS  !name changed

physics,clear
et,2,98
et,3,45
et,1,0

mp,ex,1,170e3
mp,nuxy,1,0.34
MP,DENS,1,2.329e-15
!mp,kxx,1,1.5e8
!mp,rsvx,1,2.3e-3
!mp,perx,1,11.8

da,14,ux,0
da,14,uy,0
da,14,uz,0
da,19,ux,0
da,19,uy,0
da,19,uz,0
da,38,ux,0
da,38,uy,0
da,38,uz,0
da,43,uy,0
da,43,ux,0
da,43,uz,0

vsel,s,,,3,5
asel,r,loc,z,sh
da,all,uz,0
da,all,ux,0
da,all,uy,0

allsel,all
finish
physics,write,STRUCTURE !name changed
physics,clear

!==========the following is solve command
allsel,all
ESSOLV,'ELECTROS','STRUCTURE',3,2,'air','struc',,,10,2