做了一个压电的声学例子 在求解的时候总是出现两个错误
1. Finer mesh required near center line for element 143924. 基本都是这个错误很多,好多条,这个是因为网格划得不够细吗?
2. Real constant 1 referenced by at least element types 3 and 4. 这个实常数只定义了单元类型,不晓得怎么老提示被两个单元引用了?
下面是程序语句 求各位大大及版主看看 着急啊
FINISH
/CONFIG,NRES,20000
!define const
R=5E-3
R1=8E-3
T1=1E-3
T_Al=5E-3
T_Steel=6E-3
T_Pzt=5E-3
R_Pzt=0.8E-3
W_Outerspace=1.6E-3
/PREP7
ET,1,PLANE13,7
KEYOPT,1,3,1
ET,2,PLANE42
KEYOPT,2,3,1
ET,3,FLUID29
KEYOPT,3,2,0
KEYOPT,3,3,1
ET,4,FLUID29
KEYOPT,4,2,1
KEYOPT,4,3,1
ET,5,FLUID129
KEYOPT,5,3,1
! Material properties for PZT-5, Y
mp,DENS,1,7750 ! Density, kg/m**3
MP,PERX,1,919 ! PERMITTIVITY (X AND Y DIRECTION)
MP,PERY,1,826
MP,PERZ,1,919 ! PERMITTIVITY (Z DIRECTION)
TB,PIEZ,1 ! DEFINE PIEZ. TABLE
TBDATA,2,-5.3512
TBDATA,5,15.784
TBDATA,8,-5.3512
TBDATA,10,12.295
TBDATA,15,12.295
TB,ANEL,1 ! DEFINE STRUCTURAL TABLE
TBDATA,1,1.2035E11,7.509E10,7.519E10 ! INPUT [C] MATRIX
TBDATA,7,1.1087E11,7.509E10
TBDATA,12,1.2035E11
TBDATA,16,2.1053E10
TBDATA,19,2.1053E10
TBDATA,21,2.2584E10
!*****************************************************************
MP,EX,2,7.1E10 ! ALUMINUM MODULUS OF ELASTICITY
MP,NUXY,2,.33 ! ALUMINUM POISSON RATIO
MP,DENS,2,2700 ! ALUMINUM DENSITY
!*****************************************************************
MP,EX,3,2.1E11 ! STEEL MODULUS OF ELASTICITY
MP,NUXY,3,.3 ! STEEL POISSON RATIO
MP,DENS,3,7850 ! STEEL DENSITY
!*****************************************************************
!Material properties for epoxy resin
MP,EX,4,2.5E9
MP,DENS,4,1430
MP,PRXY,4,0.3
!*****************************************************************
!Material properties for water0
MP,DENS,5,1000
MP,SONC,5,1500
MP,MU,5,0
!*****************************************************************
!Material properties for water1
MP,DENS,6,1000
MP,SONC,6,1500
MP,MU,6,1
!Real constant
R,1,0.05,0,0,
!********CREAT back layer********
RECTNG,,R,-(T_Steel+T_Pzt+T_Al+T1)/2,T_Steel-(T_Steel+T_Pzt+T_Al+T1)/2
!********CREAT piezoelectric layer********
RECTNG,,R,T_Steel-(T_Steel+T_Pzt+T_Al+T1)/2,T_Steel+T_Pzt-(T_Steel+T_Pzt+T_Al+T1)/2
!********CREAT front layer********
RECTNG,,R,T_Steel+T_Pzt-(T_Steel+T_Pzt+T_Al+T1)/2,T_Steel+T_Pzt+T_Al-(T_Steel+T_Pzt+T_Al+T1)/2
RECTNG,,R1,T_Steel+T_Pzt+T_Al-(T_Steel+T_Pzt+T_Al+T1)/2,T_Steel+T_Pzt+T_Al+T1-(T_Steel+T_Pzt+T_Al+T1)/2
!********Glue all layers********
ALLSEL,ALL
AGLUE,ALL
K,17,R_Pzt,T_Steel-(T_Steel+T_Pzt+T_Al+T1)/2
K,18,R-W_Outerspace-2*R_Pzt,T_Steel-(T_Steel+T_Pzt+T_Al+T1)/2
K,19,R-W_Outerspace,T_Steel-(T_Steel+T_Pzt+T_Al+T1)/2
K,20,R_Pzt,T_Steel+T_Pzt-(T_Steel+T_Pzt+T_Al+T1)/2
K,21,R-W_Outerspace-2*R_Pzt,T_Steel+T_Pzt-(T_Steel+T_Pzt+T_Al+T1)/2
K,22,R-W_Outerspace,T_Steel+T_Pzt-(T_Steel+T_Pzt+T_Al+T1)/2
k,23,R1,-(T_Steel+T_Pzt+T_Al+T1)/2
L,17,20
L,18,21
L,19,22
A,2,11,14,23
FLST,3,3,4,ORDE,3
FITEM,3,5
FITEM,3,-6
FITEM,3,8
ASBL, 5,P51X
cyl4,,,0,-90,0.05,90
!********Glue all areas********
ALLSEL,ALL
AOVLAP,ALL
ALLSEL,ALL
AGLUE,ALL
K,24,0,-(T_Steel+T_Pzt+T_Al+T1)
K,25,0,T_Steel+T_Pzt+T_Al+T1
K,26,2*R1,0
K,27,0.05,0
LSTR,23,24
LSTR,23,26
LSTR,15,25
LSTR,15,26
LSTR,26,27
FLST,3,5,4,ORDE,5
FITEM,3,7
FITEM,3,9
FITEM,3,18
FITEM,3,29
FITEM,3,34
ASBL, 12,P51X
ALLSEL,ALL
AGLUE,ALL
NUMCMP,ALL
NUMMRG,ALL,,,,LOW
SAVE
!********Attribute of areas********
ALLSEL,ALL
ASEL,S,,,8
ASEL,A,,,10
AATT,1,0,1,0
ASEL,S,,,1
AATT,3,0,2,0
ASEL,S,,,3
ASEL,A,,,9
AATT,4,0,2,0
ASEL,S,,,6
ASEL,A,,,7
AATT,2,0,2,0
ASEL,S,,,2
ASEL,A,,,4
ASEL,A,,,5
ASEL,A,,,11
ASEL,A,,,12
ASEL,A,,,13
AATT,5,0,4,0
!*************************划分网格*********************************
ALLSEL,ALL
LESIZE,1,0.1E-3
LESIZE,2,0.1E-3
LESIZE,3,0.1E-3
LESIZE,4,0.1E-3
LESIZE,5,0.1E-3
LESIZE,7,0.1E-3
LESIZE,9,0.1E-3
LESIZE,10,0.1E-3
LESIZE,11,0.1E-3
LESIZE,12,0.1E-3
LESIZE,13,0.1E-3
LESIZE,14,0.1E-3
LESIZE,15,0.1E-3
LESIZE,16,0.1E-3
LESIZE,18,0.1E-3
LESIZE,19,0.1E-3
LESIZE,20,0.1E-3
LESIZE,21,0.1E-3
LESIZE,22,0.1E-3
LESIZE,23,0.1E-3
LESIZE,24,0.1E-3
LESIZE,25,0.1E-3
LESIZE,26,0.1E-3
LESIZE,27,0.1E-3
LESIZE,29,0.1E-3
LESIZE,30,0.1E-3
LESIZE,32,0.1E-3
LESIZE,33,0.1E-3
LESIZE,6,,,100
LESIZE,17,,,100
LESIZE,8,,,100
LESIZE,28,,,100
LESIZE,31,,,120
LESIZE,34,,,120
LESIZE,36,,,120
LESIZE,35,,,200
LESIZE,37,,,200
ASEL,S,,,3
ASEL,A,,,8
ASEL,A,,,9
MSHKEY,1
MSHAPE,0,2D
AMESH,3
AMESH,8
AMESH,9
ALLS
ASEL,S,,,10
FLST,2,2,4,ORDE,2
FITEM,2,29
FITEM,2,-30
LCCAT,P51X
MSHKEY,1
MSHAPE,0,2D
AMESH,10
ASEL,S,,,1
FLST,2,4,4,ORDE,4
FITEM,2,12
FITEM,2,22
FITEM,2,24
FITEM,2,26
LCCAT,P51X
MSHKEY,1
MSHAPE,0,2D
AMESH,1
ASEL,S,,,6
FLST,2,4,4,ORDE,4
FITEM,2,15
FITEM,2,23
FITEM,2,25
FITEM,2,27
LCCAT,P51X
MSHKEY,1
MSHAPE,0,2D
AMESH,6
ALLSEL,ALL
ASEL,S,,,7
FLST,2,2,4,ORDE,2
FITEM,2,10
FITEM,2,20
LCCAT,P51X
MSHKEY,1
MSHAPE,0,2D
AMESH,7
ALLSEL,ALL
ASEL,S,,,11
FLST,2,3,4,ORDE,3
FITEM,2,2
FITEM,2,16
FITEM,2,18
LCCAT,P51X
MSHKEY,1
MSHAPE,0,2D
AMESH,11
ALLSEL,ALL
lsel,r,lcca
ldele,all
!###############手动划分面2、4、5、12、14 ########################
!ASEL,S,,,2
!ASEL,A,,,4
!ASEL,A,,,5
!ASEL,A,,,12
!ASEL,A,,,13
!MSHKEY,1
!MSHAPE,0,2D
!AMESH,2
!AMESH,4
!AMESH,5
!AMESH,12
!AMESH,13
!############改变流固接触面单元属性####################
LSEL,S,,,1
lsel,A,,,2
LSEL,A,,,13
lsel,A,,,14
LSEL,A,,,16
LSEL,A,,,18
LSEL,A,,,20
NSLL,S,1
ESLN,S
ESEL,u,TYPE,,1
NSLE,S
EMODIF,ALL,TYPE,3
alls
!################ Creat new element for fluid129 #########################
CSYS,1
LSEL,S,LOC,X,0.05
nsll,s,1
TYPE,5
MAT,5
REAL,1
ESURF
CSYS,0
alls
!###########################Boundary conditions################
ALLSEL,ALL
NSEL,S,LOC,Y,T_Pzt/2
NSEL,R,LOC,X,R
CP,1,VOLT,ALL
*GET,N1,NODE,0,NUM,MIN
D,ALL,VOLT,10
NSEL,R,LOC,X,R
NSEL,S,LOC,Y,-T_Pzt/2
CP,2,VOLT,ALL
D,ALL,VOLT,0
!######################水中谐振频率分析###########################
/solu
antype,harmonic
hropt,full
harfrq,30000,150000
nsubst,120
kbc,1
DMPRAT,0.042
OUTRES,ALL,ALL
ALLSEL,ALL
SOLVE
|