本帖最后由 steven710 于 2010-9-20 13:43 编辑
我现在在仿真压电换能器在空气中的声压辐射图,可是,提取到的声压却全是零,请各位高手帮帮忙呀。。。。附代码如下:
/PREP7
/TITLE
S11=16.1E-12
S12=-5.74E-12
S13=-7.22E-12
S33=20.4E-12
S44=52.6E-12
S66=44.3E-12
D15=7.58E-10
D31=-2.1E-10
D33=3.74E-10
EP11=2270
EP33=2130
RHO=7650
TB,ANEL,1,,,1
TBDATA,1,S11,S13,S12
TBDATA,7,S33,S13
TBDATA,12,S11
TBDATA,16,S44
TB,PIEZ,1,,,1
TBDATA,2,D31
TBDATA,5,D33
TBDATA,8,D31
TBDATA,10,D15
TB,DPER,1,,,1
TBDATA,1,EP11,EP33
MP,DENS,1,RHO
MP,DMPR,1,0.005
ET,1,PLANE223,1001,,1
ET,2,FLUID29
ET,3,129
R,3,0.05,0,0
ET,4,FLUID29,,1,0
MP,DENS,2,1.21
MP,SONC,2,344
R=7.5E-3
T=15.5E-3
RL=0.05
CYL4,0,0,0,90,RL,-90
RECT,0,R,-T/2,T/2
AOVLAP,1,2
ASEL,S,AREA,,3
AATT,1,1,1,0
LESIZE,6,,,2,1
LESIZE,4,,,2,1
LESIZE,5,,,4,1
LESIZE,8,,,2,1
LESIZE,9,,,2,1
MSHKEY,0
MSHAPE,0,2D
AMESH,3
ASEL,S,AREA,,4
AATT,2,1,2,0
LESIZE,10,,,10,1
LESIZE,11,,,10,1
LESIZE,1,,,30,1
MSHKEY,0
MSHAPE,0,2D
AMESH,4
ESEL,S,TYPE,,1
NSLE,S
ESLN,S,0
NSLE,S
ESEL,INVE
NSLE,S
EMODIF,ALL,TYPE,4
ESEL,ALL
NSEL,ALL
CSYS,1
NSEL,S,LOC,X,RL
TYPE,3
REAL,3
MAT,2
ESURF
ESEL,ALL
NSEL,ALL
LSEL,S,,,6
LSEL,A,,,4
LSEL,A,,,5
NSLL,S,1
ESEL,S,TYPE,,2
SF,ALL,FSI,1
NSEL,ALL
ESEL,ALL
LSEL,S,,,6
NSLL,S,1
CP,1,VOLT,ALL
*GET,N_TOP,NODE,,NUM,MIN
D,ALL,VOLT,1
LSEL,S,,,4
NSLL,S,1
CP,2,VOLT,ALL
*GET,N_BOT,NODE,,NUM,MIN
D,ALL,VOLT,0
NSEL,ALL
ESEL,ALL
FINISH
/SOLU
EQSLV,SPARSE
ANTYPE,HARMIC
HROPT,FULL
HARFRQ,10000,500000
NSUBST,500
KBC,1
SOLVE
/POST1
SET,FIRST
ESEL,S,MAT,,2
PLESOL,NMISC,1,0,1 |