hbutren 发表于 2011-5-23 20:46:00

请教各位大大及版主:ANSYS后处理换能器导纳的计算

请教各位:
      最近做了一个换能器的ANSYS模拟,谐响应分析后,计算导纳的过程时,根据莫喜平老师说的方法,提取耦合部序号最低的节点(正极耦合部的)的反作用力AMPS_2后,出现了问题:得到的AMPS_2是一个不变的值,没有虚部。这样G/B就无法继续下去。
请各位大大及版主看看哪里出了问题。

FINISH
/CONFIG,NRES,20000
!define const
R=5E-3
R_PZT=1E-3
R1=3E-3
T1=2E-3
T_Al=6E-3
T_Steel=6E-3
T_Pzt=1.5E-3
T_Epoxy=1.5E-3
R_Pzt=1.5E-3
W_Outerspace=0.2E-3

/PREP7
ET,1,PLANE13,7
KEYOPT,1,3,1
ET,2,PLANE42
KEYOPT,2,3,1
!*****************************************************************
!Material properties for PZT-5,   
TB,ANEL,1,,,0                           
TBDATA,1,1.2035e11,7.7519e10,7.509e10   
TBDATA,7,1.2035e11,7.509e10
TBDATA,12,1.1087e11
TBDATA,16,2.2584e10
TBDATA,19,2.1053e10
TBDATA,21,2.1053e10
MP,DENS,1,7750
MP,PERX,1,919   
MP,PERY,1,919   
MP,PERZ,1,826   
TB,PIEZ,1   
TBDATA,3,-5.3512
TBDATA,6,-5.3512
TBDATA,9,15.784
TBDATA,14,12.295
TBDATA,16,12.295
!*****************************************************************   
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,4800                   ! STEEL DENSITY

!*****************************************************************
!Material properties for epoxy resin
MP,EX,4,2.5E9
MP,DENS,4,1430
MP,PRXY,4,0.3

!********CREAT back layer********
RECTNG,,2*R,,T_Steel
!********CREAT piezoelectric layer********
RECTNG,,2*R,T_Steel,T_Steel+T_Epoxy
RECTNG,,2*R,T_Steel+T_Epoxy,T_Steel+T_Epoxy+T_Pzt
RECTNG,,2*R,T_Steel+T_Epoxy+T_Pzt,T_Steel+T_Epoxy+2*T_Pzt
RECTNG,,2*R,T_Steel+T_Epoxy+2*T_Pzt,T_Steel+2*T_Epoxy+2*T_Pzt
!********CREAT front layer********
RECTNG,,2*R,T_Steel+2*T_Epoxy+2*T_Pzt,T_Steel+2*T_Epoxy+2*T_Pzt+T_Al
RECTNG,R-R1,R+R1,T_Steel+2*T_Epoxy+2*T_Pzt+T_Al,T_Steel+2*T_Epoxy+2*T_Pzt+T_Al+T1
!********Glue all layers***********************************
ALLSEL,ALL
AGLUE,ALL
!********CREAT piezoelectric composite layer********
K,29,W_Outerspace,T_Steel
K,30,W_Outerspace+2*R_Pzt,T_Steel
K,31,W_Outerspace+2*R_Pzt,T_Steel+T_Epoxy
K,32,W_Outerspace,T_Steel+T_Epoxy
K,33,R-R_Pzt,T_Steel
K,34,R+R_Pzt,T_Steel
K,35,R+R_Pzt,T_Steel+T_Epoxy
K,36,R-R_Pzt,T_Steel+T_Epoxy
K,37,2*R-W_Outerspace-2*R_Pzt,T_Steel
K,38,2*R-W_Outerspace,T_Steel
K,39,2*R-W_Outerspace,T_Steel+T_Epoxy
K,40,2*R-W_Outerspace-2*R_Pzt,T_Steel+T_Epoxy
K,41,W_Outerspace,T_Steel+T_Epoxy+2*T_Pzt
K,42,W_Outerspace+2*R_Pzt,T_Steel+T_Epoxy+2*T_Pzt
K,43,W_Outerspace+2*R_Pzt,T_Steel+2*T_Epoxy+2*T_Pzt
K,44,W_Outerspace,T_Steel+2*T_Epoxy+2*T_Pzt
K,45,R-R_Pzt,T_Steel+T_Epoxy+2*T_Pzt
K,46,R+R_Pzt,T_Steel+T_Epoxy+2*T_Pzt
K,47,R+R_Pzt,T_Steel+2*T_Epoxy+2*T_Pzt
K,48,R-R_Pzt,T_Steel+2*T_Epoxy+2*T_Pzt
K,49,2*R-W_Outerspace-2*R_Pzt,T_Steel+T_Epoxy+2*T_Pzt
K,50,2*R-W_Outerspace,T_Steel+T_Epoxy+2*T_Pzt
K,51,2*R-W_Outerspace,T_Steel+2*T_Epoxy+2*T_Pzt
K,52,2*R-W_Outerspace-2*R_Pzt,T_Steel+2*T_Epoxy+2*T_Pzt
L,29,32
L,30,31
L,33,36
L,34,35
L,37,40
L,38,39
L,41,44
L,42,43
L,45,48
L,46,47
L,49,52
L,50,51
FLST,3,6,4,ORDE,5   
FITEM,3,5   
FITEM,3,-6
FITEM,3,8   
FITEM,3,-10
FITEM,3,12
ASBL,       8,P51X
FLST,3,6,4,ORDE,5   
FITEM,3,13
FITEM,3,-14
FITEM,3,16
FITEM,3,-18
FITEM,3,20
ASBL,      11,P51X

!*************************划分网格*********************************
ALLSEL,ALL
LESIZE,ALL,0.1E-3

TYPE,1
MAT,1
AMESH,2
AMESH,5
AMESH,13
AMESH,8
AMESH,17
AMESH,19
FLST,2,7,4,ORDE,7   
FITEM,2,22
FITEM,2,24
FITEM,2,42
FITEM,2,44
FITEM,2,46
FITEM,2,48
FITEM,2,50
LCCAT,P51X
AMESH,9
FLST,2,7,4,ORDE,7   
FITEM,2,3   
FITEM,2,51
FITEM,2,53
FITEM,2,55
FITEM,2,57
FITEM,2,59
FITEM,2,61
LCCAT,P51X
AMESH,10
allsel,all
lsel,r,lcca
ldele,all

TYPE,2
MAT,4
AMESH,3
AMESH,4
AMESH,6
AMESH,14
AMESH,15
AMESH,16
AMESH,18
AMESH,20

TYPE,2
MAT,3
ALLSEL,ALL
FLST,2,7,4,ORDE,7   
FITEM,2,21
FITEM,2,23
FITEM,2,41
FITEM,2,43
FITEM,2,45
FITEM,2,47
FITEM,2,49
LCCAT,P51X
AMESH,1
TYPE,2
MAT,2
ALLSEL,ALL
FLST,2,7,4,ORDE,7   
FITEM,2,7   
FITEM,2,52
FITEM,2,54
FITEM,2,56
FITEM,2,58
FITEM,2,60
FITEM,2,62
LCCAT,P51X
FLST,2,3,4,ORDE,3   
FITEM,2,25
FITEM,2,39
FITEM,2,-40
LCCAT,P51X
AMESH,12
allsel,all
lsel,r,lcca
ldele,all
AMESH,7

!**********************模态分析谐响应分析*************************
ALLSEL,ALL
/SOL
ANTYPE, 2
MODOPT, LANB, 10
EQSLV, SPAR
MXPAND, 10, , , 0
LUMPM, 0
PSTRES, 0
MODOPT, LANB, 10, 0, , , OFF
!耦合负电极为0
NSEL,S,LOC,Y,T_Steel
D,ALL,VOLT,0
NSEL,S,LOC,Y,T_Steel+T_Epoxy+T_Pzt
D,ALL,VOLT,0
NSEL,S,LOC,Y,T_Steel+2*T_Epoxy+2*T_Pzt
D,ALL,VOLT,0
!耦合正电极
NSEL,S,LOC,Y,T_Steel+T_Epoxy
CP,1,VOLT,ALL
D,all,VOLT,100
NSEL,s,LOC,Y,T_Steel+T_Epoxy+2*T_Pzt
CP,2,VOLT,ALL
*get,n1,node,,num,min
D,n1,VOLT,100
ALLSEL,ALL
SOLVE
FINISH
ALLSEL,ALL
/SOL
ANTYPE,HARM
HROPT,FULL
HROUT,OFF   
LUMPM,0
EQSLV,FRONT,0
PSTRES,0
HARFRQ,50000,150000,
NSUBST,250,
KBC,1
DMPRAT,0.006
OUTRES,ALL,ALL
ALLSEL,ALL
SOLVE
FINISH

hbutren 发表于 2011-5-24 10:24:20

呃发现问题了 但是得到的导纳图却怪怪的,不是很理想:L

tomance 发表于 2012-3-14 20:44:20

给个图看看先。

dream1376 发表于 2014-8-31 23:04:35

l楼主发现的问题能说一下吗,我也在做这个结果差不多

cnpclym 发表于 2014-9-9 19:02:46

同问,求指导

江南烟雨天272 发表于 2014-12-9 08:57:15

楼主我也是这个问题,请问你是怎么解决的?
页: [1]
查看完整版本: 请教各位大大及版主:ANSYS后处理换能器导纳的计算