!!!Single Trans126 element
/prep7
mwidth=40 ! square membrane
melectrode=40
height=0.85 !
Al_height=0.53
bracelength=10
bracewidth=2
k_sio=3.9
d_sio=1.7
hgap=0.53
vlt=42!Bias voltage
deta=3.9968e-14
et,1,solid45
r,1
et,2,trans126,,2
c0=8.854e-6*mwidth*mwidth ! C0/x constant for Capacitance equation
r,2,0,0,(hgap+d_sio/k_sio),0.05,.1
rmore,c0 ! Real constant C0
mp,ex,1,0.69e5 !SiO2 Membrane
mp,nuxy,1,0.17
mp,dens,1,2.2e-15
mp,ex,2,0.7e5 !Al electrode
mp,nuxy,2,0.36
mp,dens,2,2.7e-15
k,1,mwidth/2-bracewidth,0,mwidth/2
k,2,mwidth/2,0,mwidth/2
k,3,mwidth/2+bracelength,0,mwidth/2+bracelength
k,4,mwidth/2-bracewidth+bracelength,0,mwidth/2+bracelength
k,5,mwidth/2-bracewidth,-(height+Al_height),mwidth/2
k,6,mwidth/2,-(height+Al_height),mwidth/2
k,7,mwidth/2+bracelength,-(height+Al_height),mwidth/2+bracelength
k,8,mwidth/2-bracewidth+bracelength,-(height+Al_height),mwidth/2+bracelength
v,1,2,3,4,5,6,7,8
k,9,mwidth/2,0,mwidth/2-bracewidth
k,10,mwidth/2+bracelength,0,mwidth/2+bracelength-bracewidth
k,11,mwidth/2,-(height+Al_height),mwidth/2-bracewidth
k,12,mwidth/2+bracelength,-(height+Al_height),mwidth/2+bracelength-bracewidth
v,2,9,10,3,6,11,12,7
vsymm,x,all,,1,0,0
vsel,all
vsymm,z,all,,4,0,0
block,-mwidth/2,mwidth/2,0,-(height+Al_height),-mwidth/2,mwidth/2
vsel,all
wprota,,-90
wpoffs,,,-Al_height
rect,-60,60,-60,60
vsba,all,51!51 is the rectng's number
wpoffs,,,Al_height
wprota,,90
vglue,all
nummrg,all
numcmp,all
vsel,all
vsel,s,loc,y,-height-Al_height,-Al_height
vatt,1,1,1
vsel,s,loc,y,-Al_height,0
vatt,2,1,1
vsel,all
vsel,u,volu,,1
vsel,u,volu,,2
cm,BRACE_VOLU,volu
cmsel,s,BRACE_VOLU
cmplots
lsel,s,loc,y,-height/2-Al_height
lesize,all,,,height
lsel,s,loc,y,-Al_height/2
lesize,all,,,Al_height
asel,s,loc,y,0
asel,a,loc,y,-height-Al_height
asel,a,loc,y,-Al_height
lsla,s
lesize,all,bracewidth/2
mshape,0,3D
mshkey,1
vmesh,BRACE_VOLU
!Must ccat area firstly,then ccat line
asel,s,loc,x,-mwidth/2,-mwidth/2
asel,r,loc,y,-height-Al_height,-Al_height
accat,all
asel,s,loc,x,-mwidth/2,-mwidth/2
asel,r,loc,y,-Al_height,0
accat,all
asel,s,loc,z,-mwidth/2,-mwidth/2
asel,r,loc,y,-height-Al_height,-Al_height
accat,all
asel,s,loc,z,-mwidth/2,-mwidth/2
asel,r,loc,y,-Al_height,0
accat,all
asel,s,loc,x,mwidth/2,mwidth/2
asel,r,loc,y,-height-Al_height,-Al_height
accat,all
asel,s,loc,x,mwidth/2,mwidth/2
asel,r,loc,y,-Al_height,0
accat,all
asel,s,loc,z,mwidth/2,mwidth/2
asel,r,loc,y,-height-Al_height,-Al_height
accat,all
asel,s,loc,z,mwidth/2,mwidth/2
asel,r,loc,y,-Al_height,0
accat,all
lsel,s,loc,y,0,0
lsel,r,loc,x,-mwidth/2,-mwidth/2
lccat,all
lsel,s,loc,y,0,0
lsel,r,loc,z,-mwidth/2,-mwidth/2
lccat,all
lsel,s,loc,y,0,0
lsel,r,loc,x,mwidth/2,mwidth/2
lccat,all
lsel,s,loc,y,0,0
lsel,r,loc,z,mwidth/2,mwidth/2
lccat,all
lsel,s,loc,y,-Al_height,-Al_height
lsel,r,loc,x,-mwidth/2,-mwidth/2
lccat,all
lsel,s,loc,y,-Al_height,-Al_height
lsel,r,loc,z,-mwidth/2,-mwidth/2
lccat,all
lsel,s,loc,y,-Al_height,-Al_height
lsel,r,loc,x,mwidth/2,mwidth/2
lccat,all
lsel,s,loc,y,-Al_height,-Al_height
lsel,r,loc,z,mwidth/2,mwidth/2
lccat,all
lsel,s,loc,y,-height-Al_height,-height-Al_height
lsel,r,loc,x,-mwidth/2,-mwidth/2
lccat,all
lsel,s,loc,y,-height-Al_height,-height-Al_height
lsel,r,loc,z,-mwidth/2,-mwidth/2
lccat,all
lsel,s,loc,y,-height-Al_height,-height-Al_height
lsel,r,loc,x,mwidth/2,mwidth/2
lccat,all
lsel,s,loc,y,-height-Al_height,-height-Al_height
lsel,r,loc,z,mwidth/2,mwidth/2
lccat,all
vsel,s,volu,,1
vsel,a,volu,,2
cm,MEM_VOLU,volu
cmsel,s,MEM_VOLU
cmplot
lsel,s,loc,y,-height/2-Al_height
lesize,all,,,height
lsel,s,loc,y,-Al_height/2
lesize,all,,,Al_height
asel,s,loc,y,0
asel,a,loc,y,-height-Al_height
asel,a,loc,y,-Al_height
lsla,s
lesize,all,bracewidth/2
mshape,0,3D
mshkey,1
vmesh,MEM_VOLU
nummrg,all
numcmp,all
asel,s,loc,x,-(mwidth/2+bracelength)
asel,a,loc,x,(mwidth/2+bracelength)
asel,a,loc,z,-(mwidth/2+bracelength)
asel,a,loc,z,(mwidth/2+bracelength)
lsla,s,1
nsll,s
d,all,ux
d,all,uy
d,all,uz !constraintion
n,50000,-deta,-(height+Al_height)-(hgap+d_sio/k_sio),-deta
type,2
real,2
e,50000,5183
nsel,s,NODE,,50000
d,50000,ux
d,50000,uy
d,50000,uz
d,50000,volt,0
nsel,s,NODE,,5183
d,5183,volt,vlt
allsel,all
/solu
antype,static ! Static analysis
pstres,on ! turn on prestress effects
solve
save
finish
/post1
!===get the average uy of the Membrane nodes
nsel,s,loc,y,-(height+Al_height)
nsel,r,loc,x,-mwidth/2,mwidth/2
nsel,r,loc,z,-mwidth/2,mwidth/2
*GET, nnodemin, NODE, 0, num, min
*GET, nnodecount, NODE, 0, count
*DIM,nuy,,nnodecount
*vget,nuy,node,nnodemin,u,y
ynodesum=0
*do,i,1,nnodecount
ynodesum=ynodesum+nuy(i)
!k_y(i)=nsy(i)/nuy(i)
*enddo
ynodeavg=ynodesum/nnodecount
ynodeavgperv=ynodeavg/vlt
yobj=1+ynodeavgperv
!get the max von mise stress
nsel,all
NSORT, S, EQV
*GET, SMAX, SORT,, MAX
PLNSOL,UY
finish |