压电分析建模,理应该分析结果位移曲线是抛物线,中点位移最大,可是自己画出的图形却是直线,而且中间的位移也很奇怪。
现将code列出,求大家指点
/BATCH
/title,pztplate
!/filnam,pztplate
/prep7
et,1,5,3
et,2,45
!material properties of the PZT actuator
!mp,dens,1,7800
MP,DENS,1,7800 !DENSITY OF PIEZ
MP,PERX,1,3.186e-8 !PERMITTIVITY (X Y AND Z DIRECTION)
MP,PERY,1,3.186e-8
MP,PERZ,1,3.363e-8
TB,PIEZ,1,,,1 !DEFINE PIEZ. TABLE
TBDATA,16,741e-12 !E61 PIEZOELECTRIC CONSTANT
TBDATA,14,741e-12 !E52 PIEZOELECTRIC CONSTANT
TBDATA,3,-310e-12 !E13 PIEZOELECTRIC CONSTANT
TBDATA,6,-310e-12 !E23 PIEZOELECTRIC CONSTANT
TBDATA,9,593e-12 !E33 PIEZOELECTRIC CONSTANT
TB,ANEL,1,,,1 !DEFINE STRUCTURAL TABLE
TBDATA,1,16.5e-12,-4.78e-12,-8.45e-12 !INPUT [c] MATRIX
TBDATA,7,16.5e-12,-8.45e-12
TBDATA,12,20.7e-12
TBDATA,16,43.5e-12
TBDATA,19,43.5e-12
TBDATA,21,42.6e-12
!material properties of the aluminium
mp,dens, 2,2.71e3
mp,ex, 2,70e9
mp,nuxy, 2,0.33
length1=.019
width1=.019
length2=.038
width2=.038
thick1=.000127
thick2=.000508
volts=-15
k,1, 0,0,0
k,2, length2,width2,0
k,3, length2,0,0
k,4, 0,width2,0
k,5, 0,0,-thick2
k,6, length2,width2,-thick2
k,7, length2,0,-thick2
k,8, 0, width2,-thick2
k,9, 0.0095,0.0095,0
k,10,0.0095+length1,0.0095+width1,0
k,11,0.0095,0.0095+width1,0
k,12,0.0095+length1,0.0095,0
k,13, 0.0095,0.0095,thick1
k,14,0.0095+length1,0.0095+width1,thick1
k,15,0.0095,0.0095+width1,thick1
k,16,0.0095+length1,0.0095,thick1
elsize=.001
v,9,11,10,12,13,15,14,16
type,1
mat,1
esize,elsize
mshape,0,3d
mshkey,1
vmesh,all
V,1,3,2,4,5,7,6,8
type,1
mat,2
esize,elsize
mshape,0,3d
mshkey,1
vmesh,all
nummrg,node
nsel,s,loc,x,0
nsel,r,loc,z,-thick2,thick1
d,all,ux,0
d,all,uy,0
d,all,uz,0
!d,all,rotx,0
!d,all,rotz,0
alls
nsel,s,loc,x,length2
nsel,r,loc,z,-thick2,thick1
d,all,ux,0
d,all,uy,0
d,all,uz,0
!d,all,rotx,0
!d,all,rotz,0
alls
nsel,s,loc,y,0
nsel,r,loc,z,-thick2,thick1
d,all,ux,0
d,all,uy,0
d,all,uz,0
!d,all,roty,0
!d,all,rotz,0
alls
nsel,s,loc,y,length2
nsel,r,loc,z,-thick2,thick1
d,all,ux,0
d,all,uy,0
d,all,uz,0
!d,all,roty,0
!d,all,rotz,0
alls
inner1=0*length2
outer1=1*length2
inner2=0*length2
outer2=0*length2
nsel,s,loc,y,0.0095,0.0285
nsel,r,loc,x,0.0095,0.0285
nsel,r,loc,z,0
d,all,volt,volts
alls
nsel,s,loc,y,0,0.038
nsel,r,loc,x,0,0.038
nsel,r,loc,z,0
d,all,volt,volts
alls
nsel,s,loc,y,0.0095,0.0285
nsel,r,loc,x,0.0095,0.0285
nsel,r,loc,z,thick1
d,all,volt,0
alls
/SOLU
antyp,static
!pstres,on
solve
!nsel,s,loc,z,thick1
!nsel,s,loc,y,0.0095,0.0285
!nsel,s,loc,x,length2/2
/post1
path,disp,2,,10,
ppath,1,,0,length2/2,thick1
ppath,2,,0.038,length2/2,thick1
pdef,zdisp,u,z
/output,output.txt
plpath,zdisp
pasave
fini |