如题,命令流如下
/title,PZT 2 CANTILEVER BEAM
/prep7
mp,dens,7.5e-15
mp,perx,1,1.7e9 ! permittivity
mp,pery,1,1.7e9
mp,perz,1,1.47e9
tb,ANEL,1 ! stiffness C, UN/um^2
tbdata,1,12.6E4,7.95E4,8.41E4 ! c11,c12,c13
tbdata,7,12.6E4,8.41E4 ! c11,c33
tbdata,12,11.7E4 ! c33
tbdata,16,2.30E4
tbdata,19,2.30E4
tbdata,21,2.35E4
tb,PIEZ,1 ! PIEZ constant, pC/um^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,7.8E-15 !Kg/
mp,EX,2,209e3 !KPa
mp,nuxy,2,0.3
!/define dimension of PZT
L=200
W=20
H =2
!/define dimension of main beam
L1=300
H1=2
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,20,,W,,-1
wpoffs,20,,,
block,0,L1,,W,,-H1
wpoffs,L1,,,
block,,-L,,W,,H
vglue,all
wpoffs,-L,,-H
wprota,,,90
VSBW,ALL
wpoffs,-H1
wpoffs,,,L-L1
wpoffs,,,-20
wpoffs,,,L-L1
VSBW,ALL
WPOFFS,1
WPOFFS,,,L1-L
wprota,,,-90
vsbw,all
block,-90,110,,w,1,3
wpoffs,-90
wprota,,,90
vsbw,all
wpoffs,,,200
vsbw,all
wpoffs,,,-90
vsbw,all
wpoffs,,,-20
vsbw,all
vglue,all
et,1,solid5,3 !define the element of PZT
et,2,solid45 !define the element of beam
LESIZE,ALL,10, , , ,1, , ,1
!meshing
mat,1 !left PZT
vmesh,1
mat,1 !right PZT
vmesh,7
mat,1 !middle PZT
vmesh,2
mat,1
vmesh,21
mat,1
vmesh,25
mat,2
vmesh,3,5,2
mat,2
vmesh,15,17,2
mat,2
vmesh,8,10
mat,2
vmesh,12,14,2
mat,2
vmesh,20
mat,2
vmesh,22,24
!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
nsel,r,loc,x,2*L1+20-L,2*L1+20
cm,CM_2,node
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+20,2*L1+20
cm,CM_4,node
cp,4,volt,all
nsel,s,loc,z,H+H1 !coupling the connected areas between the middle PZT and the MIAN Beam
nsel,r,loc,x,L1+10-L/2,L1+10+L/2
cm,CM_5,node
cp,5,volt,all
fini
/solu
antype,static
nsel,s,loc,x,0
d,all,ux,0,,,,uy,uz
nsel,s,loc,x,2*L1+20
d,all,ux,0,,,,uy,uz
nsel,all
d,n_left,volt,0
d,n_right,volt,0
nsel,all
F,753,FZ,-1
solve
fini
/post1
输出的最大电压见附件。。由于压电陶瓷对形变敏感,所以我估计电压是要大点,谁帮我看看,是不是大到几V的电压。。。命令流里面用的是微米的单位。。加载的力是1微牛
另外谁帮推荐个做极限载荷的例子给我看看,我想对中间那块地方做极限载荷分析 |