找回密码
 注册
Simdroid-非首页
查看: 172|回复: 5

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

[复制链接]
发表于 2011-5-23 20:46:00 | 显示全部楼层 |阅读模式 来自 湖南湘潭
请教各位:
      最近做了一个换能器的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
 楼主| 发表于 2011-5-24 10:24:20 | 显示全部楼层 来自 湖南湘潭
Simdroid开发平台
呃  发现问题了 但是得到的导纳图却怪怪的,不是很理想
回复 不支持

使用道具 举报

发表于 2012-3-14 20:44:20 | 显示全部楼层 来自 黑龙江哈尔滨
给个图看看先。
回复 不支持

使用道具 举报

发表于 2014-8-31 23:04:35 | 显示全部楼层 来自 江苏扬州
l楼主发现的问题能说一下吗,我也在做这个结果差不多
回复 不支持

使用道具 举报

发表于 2014-9-9 19:02:46 | 显示全部楼层 来自 陕西西安
同问,求指导
回复 不支持

使用道具 举报

发表于 2014-12-9 08:57:15 | 显示全部楼层 来自 中国
楼主我也是这个问题,请问你是怎么解决的?
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-3-29 07:30 , Processed in 0.038167 second(s), 13 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表