依旧是双梁的问题。我分别左右电极接地以后发现输出的最大电压有差别,谁知道这是什么原因造成的啊。。那位帮我看看。
命令流如下。。short-circuit哪里分别先左接地求解,再右
/title,PZT 2 CANTILEVER BEAM
/UNIT,SI ! Unit setting (International Standard)
/prep7
mp,dens,1,7500
mp,perx,1,1700 ! dielectric constant
mp,pery,1,1700
mp,perz,1,1470
tb,ANEL,1 ! stiffness C, N/m^2
tbdata,1,12.6E10,7.95E10,8.41E10 ! c11,c12,c13
tbdata,7,12.6E10,8.41E10 ! c11,c33
tbdata,12,11.7E10 ! c33
tbdata,16,2.30E10
tbdata,19,2.30E10
tbdata,21,2.35E10
tb,PIEZ,1 ! PIEZ constant, C/m^2
tbdata,3,-6.5 ! e31
tbdata,6,-6.5 ! e31
tbdata,9,23.3 ! e33
tbdata,11,17.0 ! e15
tbdata,13,17.0 ! e15
!/define main beam material properties
mp,dens,2,7800
mp,EX,2,209e9
mp,nuxy,2,0.3
!/define dimension of PZT
L=2E-4
W=2E-5
H =2e-6
!/define dimension of main beam
L1=3E-4
H1=2E-6
block,,L,0,W,0,H !left PZT
block,0,L1,0,W,H,H+H1 !MAIN BEAM
wpoffs,L1,0,H+H1
block,0,2e-5,,W,,-1E-6
wpoffs,2e-5,,,
block,0,L1,,W,,-H1
wpoffs,L1,,,
block,,-L,,W,,H
vglue,all
et,1,solid5,3 !define the element of PZT
et,2,solid45 !define the element of beam
wpoffs,-L,,-H
wprota,,,90
VSBW,ALL
wpoffs,-H1
wpoffs,,,L-L1
wpoffs,,,-2e-5
wpoffs,,,L-L1
VSBW,ALL
WPOFFS,1E-6
WPOFFS,,,L1-L
wprota,,,-90
vsbw,all
!MESHING
LESIZE,ALL,1e-5, , , ,1, , ,1
mat,1 !left PZT
vmesh,1
mat,1 !right PZT
vmesh,7
mat,2
vmesh,8,15 !main beam
mat,2
vmesh,6
!define the electrode
nsel,s,loc,z,0 !define the left electrode
cp,1,volt,all
*get,n_left,node,0,num,min
nsel,s,loc,z,H+H1+H !define the right electrode
cp,2,volt,all
*get,n_right,node,0,num,min
nsel,s,loc,z,H !coupling the connected areas between the left PZT and MAIN BEAM
nsel,r,loc,x,0,L
cm,CM_3,node
cp,3,volt,all
nsel,s,loc,z,H+H1 !coupling the connected areas between the right PZT and the MIAN Beam
nsel,r,loc,x,2*L1-L+2e-5,2*L1+2e-5
cm,CM_4,node
cp,4,volt,all
fini
/solu
antype,static
nsel,s,loc,x,0
d,all,ux,0,,,,uy,uz
nsel,s,loc,x,2*L1+2e-5
d,all,ux,0,,,,uy,uz
nsel,all
d,n_left,volt,0
d,n_right,volt,0
F,690,FZ,1E-6
solve |