求助 壓電薄膜PVDF的聲壓分析
本帖最后由 chenian925 于 2014-4-13 17:45 编辑最近再做一個寬2mm，深10mm的PVDF薄膜的聲壓計算
但是聲壓輻射圖，看不出來有輻射狀況，
另外，我如果想要算一段頻率，在30cm的聲壓曲線，請問該怎麼做？
跪求各路高手前輩指點迷津~
下面是我寫的ADPL
FINISH
/CLEAR,START
/NERR,,99999
!*==============================================================================
!* Input Variable
!*==============================================================================
!* Mechanical Constat
!* Young's Modulus
YEX=(2.0E9)*1 !PVDF
!* Poisson's Ratio
VXY=0.29
!* Shear Modulus G=E/2(1+v)
SGXY=YEX/(2*(1+VXY))
!* Density (Kg/m^3)
DENSITY=1780
!* Geometry
!* WIDTH
W=10E-3
!* RADIUS
RAD=2E-3
!* Thickness
TPVDF=28E-6
!*==============================================================================
!* Dielectric Constant, Permittivity
!*==============================================================================
EP33=12
!*==============================================================================
!* Piezoelectric Matrix: :stress matrix, :strain matrix
!*===============================================================================
ED31=23E-12
!ED32=3E-12
ED33=-33E-12
!*==============================================================================
/PREP7
!*
ET,1,SOLID5,3 !PVDF
ET,3,FLUID30
KEYOPT,3,2,0
ET,4,FLUID130
R,4,300E-3,0,0
ET,5,FLUID30
KEYOPT,5,2,1
!*********************
!* Material Constant *
!*********************
MP,EX,1,YEX
MP,NUXY,1,VXY
MP,GXY,1,SGXY
!********************
!* Density (Kg/m^3) *
!********************
MP,DENS,1,DENSITY
!********************
!* Speed of Sound *
!* Density of Air *
!********************
MP,DENS,3,1.2
MP,SONC,3,344
!************************************
!* Dielectric Matrix (Permittivity) *
!* PVDF: EP11=EP22=EP33 *
!************************************
MP,PERX,1,EP33
MP,PERY,1,EP33
MP,PERZ,1,EP33
!**********************************
!* Piezoelectric Strain Matrix *
!**********************************
TB,PIEZ,1
TBDATA,3,ED31
TBDATA,9,ED33
!************
!* Modeling *
!************
!*
LOCAL,11,0,0,0,0,0, 0,0,0,0
LOCAL,12,0,0,0,0,0,-90,0, ,
CSYS,11
!* PVDF
BLC4,-RAD/2,0,RAD/2,TPVDF,-W
BLC4,0,0,RAD/2,TPVDF,-W
!*
WPROTA,0,-90,0
SPHERE,,300E-3,0,180
SPHERE,,100E-3,0,180
WPROTA,0,90,0
!*
VPTN,ALL
VGLUE,All
VPLOT,ALL
!*
!********
!* Mesh *
!********
!*
!* PVDF
CSYS,12
TYPE,1
MAT,1
SMRTSIZE,1
MSHKEY,0
MSHAPE,1,3D
ESYS,12
VMESH,1
VMESH,5
!*
ALLSEL,ALL
!*
CSYS,11
!*
!****************
!* MESH OF AIR
!****************
CSYS,11
VATT,3,,3,11
!*
!DESIZE,,5,15
SMRTSIZE,1
MSHKEY,0
MSHAPE,1,3D
ESYS,11
VMESH,6,7
ALLSEL,ALL
!*
CSYS,11
ESEL,S,TYPE,,1
NSLE,S,ALL
ESLN,S,0
NSLE,S
ESEL,INVE
NSLE,S,1
EMODIF,ALL,TYPE,5
ALLSEL,ALL
!*
CSYS,11
ASEL,S,AREA,,13
ASEL,A,AREA,,14
NSLA,S,1
ESYS,11
TYPE,4
REAL,4
MAT,3
ESURF
ALLSEL,ALL
!*
CSYS,11
ASEL,S,LOC,Y,0 !RAD
ASEL,R,LOC,X,-RAD/2,RAD/2 !RAD
ASEL,A,LOC,Y,TPVDF !(RAD+TPVDF)
NSLA,S,1
ESEL,S,TYPE,,3
SF,ALL,FSI
ALLSEL,ALL
CSYS,11
!*
!**********************
!* Boubdary Condition *
!**********************
!* DOF
CSYS,11
LSEL,S,LOC,Y,0,
LSEL,U,LOC,X,0,
LSEL,U,LOC,Z,0,
LSEL,U,LOC,Z,-W,
NSLL,S,1
D,ALL,ALL,0
ALLSEL,ALL
!*
!* TOP SURFACE ELECTRODE
CSYS,11
ASEL,S,LOC,Y,TPVDF !(RAD+TPVDF)
NSLA,S,1
CP,1,VOLT,ALL
D,ALL,VOLT,1
ALLSEL,ALL
!*
!* GROUND BOTTOM ELECTRODE
CSYS,11
ASEL,S,LOC,Y,0 !RAD
ASEL,R,LOC,X,-RAD/2,RAD/2 !RAD
NSLA,S,1
CP,2,VOLT,ALL
D,ALL,VOLT,0
ALLSEL,ALL
!*
CSYS,11
!*
!*SYMM
CSYS,11
NSEL,S,LOC,Z,0
DSYM,SYMM,Z
ALLSEL,ALL
!*
CSYS,11
!*
/SOLU
EQSLV,SPARSE
!EQSLV,JCG
ANTYPE,HARMIC
HROPT,FULL
NSUBST,1
KBC,1
HARF,6856,6856
SOLVE
!*
/POST1
!CSYS,12
SET,FIRST
NSEL,ALL
PLNSOL,PRES
是不是我的問題描述不清楚？
都沒有人回應～～@@ 是不是在空气汇总传播？你设置了一个空气球包围在其中？
页:
[1]