Antony21 发表于 2007-6-23 17:32

三项耦合的分析--Help都不知道该找哪部分了

提供我所建的一个MODEL，其中最底部是PZT压电材料附在铜片上的压电蜂鸣片，
结合不锈钢环形结构，里面再注满水，量测蜂鸣片的共振频率，依据实验结果，
不同体积的水会有不同的频率，但我做的结果频率都一样，
参考HELP的MODEL ANALYSIS和FLUID ANALYSIS中的Acou stics部分，似乎还是不够，
不知道有人可以解决吗？

/PREP7
ET,1,SOLID5   ! DOF:UX UY UZ
KEYOPT,1,1,2
KEYOPT,1,3,0
KEYOPT,1,5,0
ET,2,SOLID5   ! DOF:UX UY UZ VOLT
KEYOPT,2,1,3
KEYOPT,2,3,0
KEYOPT,2,5,0
ET,3,FLUID30    ! Acoustic elementinterface   
KEYOPT,3,2,0
ET,4,FLUID30 ! Acoustic elementno-interface   
KEYOPT,4,2,1
MP,EX,1,100E9      ! Brass MODULUS OF ELASTICITY
MP,NUXY,1,0.34   ! Brass POISSON RATIO
MP,DENS,1,8500   ! Brass DENSITY
MP,DENS,2,7650   ! BM500=PZT-5A DENSITY   
MP,EX,3,193E9      ! Stainless steel MODULUS OF ELASTICITY
MP,NUXY,3,0.3      ! Stainless steel POISSON RATIO
MP,DENS,3,8000   ! Stainless steel DENSITY
MP,VISC,4,0.001    ! Water
MP,DENS,4,998   
MP,SONC,4,1495
MP,MU,4,0.04
MP,EX,4,2.2E9   
TB,PIEZ,2   
TBDATA,3,-5.4      ! DEFINE PIEZO. TABLE FOR 500=PZT-5A
TBDATA,6,-5.4      ! PIEZO MATRIX CONSTANTS
TBDATA,9,15.8   
TBDATA,11,12.3
TBDATA,13,12.3
MP,PERX,2,916      ! PERMITTIVITY (X direction)
MP,PERY,2,916      ! PERMITTIVITY (Y direction)
MP,PERZ,2,830      ! PERMITTIVITY (Z direction)
TB,ANEL,2       ! DEFINE STRUCTURAL TABLE FOR BM500=PZT-5A
TBDATA,1,121E9,75.4E9,75.2E9! INPUT MATRIX FOR BM500=PZT-5A   
TBDATA,7,121E9,75.2E9   
TBDATA,12,111E9
TBDATA,16,21.1E9
TBDATA,19,21.1E9
TBDATA,21,22.6E9
K,1,0,0,0       !Build model
K,2,0.0095,0,0
K,3,0.0095,0,0.0002
K,4,0,0,0.0002
K,5,0,0,0.0004
K,6,0.0095,0,0.0004
K,7,0.011,0,0.0004
K,8,0.011,0,0.0002
K,9,0.012,0,0.0002
K,10,0.012,0,0.0004
K,11,0.0135,0,0.0004
K,12,0.0135,0,0.0002
K,13,0.0135,0,0.0014
K,14,0.012,0,0.0014
K,15,0.011,0,0.0014
K,16,0.0095,0,0.0014
K,17,0,0,0.0014
K,18,0,0,0.0074
K,19,0.0095,0,0.0074
K,20,0.011,0,0.0074
K,21,0.012,0,0.0074
K,22,0.0135,0,0.0074
K,23,0,0,0.0084
K,24,0.0095,0,0.0084
K,25,0.011,0,0.0084
K,26,0.012,0,0.0084
K,27,0.0135,0,0.0084
A,1,2,3,4   
A,3,4,5,6   
A,3,6,7,8   
A,7,8,9,10
A,9,10,11,12
A,10,11,13,14   
A,7,10,14,15
A,6,7,15,16
A,5,6,16,17
A,16,17,18,19   
A,15,16,19,20   
A,14,15,20,21   
A,13,14,21,22   
A,18,19,24,23   
A,19,20,25,24   
A,20,21,26,25   
A,21,22,27,26   
ASEL,ALL
LSEL,A,LINE,,ALL
K,900,0,0,0
K,901,0,0,0.01
VROTAT,ALL,,,,,,900,901,360,4   
TYPE,   4   
MAT,       4
REAL,   
ESYS,       0   
SECNUM,
!*
ESIZE,0.001,0,
! VPLOT
/UI,MESH,OFF
FLST,5,8,6,ORDE,8   
FITEM,5,10
FITEM,5,-11
FITEM,5,27
FITEM,5,-28
FITEM,5,44
FITEM,5,-45
FITEM,5,61
FITEM,5,-62
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y
!*
VMESH,_Y1   
!*
CMDELE,_Y   
CMDELE,_Y1
CMDELE,_Y2
TYPE,   3   
MAT,       4
REAL,   
ESYS,       0   
SECNUM,
!*
ESIZE,0.001,0,
FLST,5,28,6,ORDE,20
FITEM,5,7   
FITEM,5,-9
FITEM,5,12
FITEM,5,14
FITEM,5,-16
FITEM,5,24
FITEM,5,-26
FITEM,5,29
FITEM,5,31
FITEM,5,-33
FITEM,5,41
FITEM,5,-43
FITEM,5,46
FITEM,5,48
FITEM,5,-50
FITEM,5,58
FITEM,5,-60
FITEM,5,63
FITEM,5,65
FITEM,5,-67
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y
!*
VMESH,_Y1   
!*
CMDELE,_Y   
CMDELE,_Y1
CMDELE,_Y2
TYPE,   1   
MAT,       3
REAL,   
ESYS,       0   
SECNUM,
!*
ESIZE,0.001,0,
FLST,5,12,6,ORDE,12
FITEM,5,6   
FITEM,5,13
FITEM,5,17
FITEM,5,23
FITEM,5,30
FITEM,5,34
FITEM,5,40
FITEM,5,47
FITEM,5,51
FITEM,5,57
FITEM,5,64
FITEM,5,68
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y
!*
VMESH,_Y1   
!*
CMDELE,_Y   
CMDELE,_Y1
CMDELE,_Y2
TYPE,   1   
MAT,       1
REAL,   
ESYS,       0   
SECNUM,
!*
ESIZE,0.001,0,
FLST,5,16,6,ORDE,8
FITEM,5,2   
FITEM,5,-5
FITEM,5,19
FITEM,5,-22
FITEM,5,36
FITEM,5,-39
FITEM,5,53
FITEM,5,-56
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y
!*
VMESH,_Y1   
!*
CMDELE,_Y   
CMDELE,_Y1
CMDELE,_Y2
TYPE,   2   
MAT,       2
REAL,   
ESYS,       0   
SECNUM,
!*
ESIZE,0.001,0,
MSHAPE,0,3D
MSHKEY,1
FLST,5,4,6,ORDE,4   
FITEM,5,1   
FITEM,5,18
FITEM,5,35
FITEM,5,52
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y
!*
VMESH,_Y1   
!*
CMDELE,_Y   
CMDELE,_Y1
CMDELE,_Y2
FLST,5,4,5,ORDE,4   
FITEM,5,18
FITEM,5,73
FITEM,5,128
FITEM,5,183
ASEL,S, , ,P51X
NSEL,ALL
NSLA,R,1
! NPLOT   
FLST,2,1009,1,ORDE,2
FITEM,2,13691   
FITEM,2,-14699
/GO
!*
D,P51X,VOLT,10
ALLSEL,ALL
FLST,5,4,6,ORDE,4   
FITEM,5,1   
FITEM,5,18
FITEM,5,35
FITEM,5,52
VSEL,S, , ,P51X
ASEL,ALL

ASEL,ALL
ASLV,R
FLST,5,4,5,ORDE,4   
FITEM,5,20
FITEM,5,75
FITEM,5,130
FITEM,5,185
ASEL,R, , ,P51X
NSEL,ALL
NSLA,R,1
! NPLOT   
FLST,2,1009,1,ORDE,8
FITEM,2,12322   
FITEM,2,-12592
FITEM,2,12688   
FITEM,2,-12939
FITEM,2,13030   
FITEM,2,-13281
FITEM,2,13372   
FITEM,2,-13605
/GO
!*
D,P51X,VOLT,0   
ALLSEL,ALL
FLST,5,16,6,ORDE,16
FITEM,5,5   
FITEM,5,-6
FITEM,5,13
FITEM,5,17
FITEM,5,22
FITEM,5,-23
FITEM,5,30
FITEM,5,34
FITEM,5,39
FITEM,5,-40
FITEM,5,47
FITEM,5,51
FITEM,5,56
FITEM,5,-57
FITEM,5,64
FITEM,5,68
VSEL,S, , ,P51X
NSEL,ALL
NSLV,R,1
! NPLOT   
FLST,2,2160,1,ORDE,33   
FITEM,2,8072
FITEM,2,8092
FITEM,2,-8128   
FITEM,2,8419
FITEM,2,-8532   
FITEM,2,8842
FITEM,2,-8860   
FITEM,2,8879
FITEM,2,-8914   
FITEM,2,9185
FITEM,2,-9292   
FITEM,2,9581
FITEM,2,-9598   
FITEM,2,9617
FITEM,2,-9652   
FITEM,2,9923
FITEM,2,-10030
FITEM,2,10319   
FITEM,2,-10336
FITEM,2,10354   
FITEM,2,-10387
FITEM,2,10639   
FITEM,2,-10740
FITEM,2,11009   
FITEM,2,-12321
FITEM,2,12631   
FITEM,2,-12687
FITEM,2,12976   
FITEM,2,-13029
FITEM,2,13318   
FITEM,2,-13371
FITEM,2,13640   
FITEM,2,-13690
!*
/GO
D,P51X, ,0, , , ,UX,UY,UZ, , ,
ALLSEL,ALL
! EPLOT   
ACEL,0,0,-9.81,
FLST,5,28,6,ORDE,20
FITEM,5,7   
FITEM,5,-9
FITEM,5,12
FITEM,5,14
FITEM,5,-16
FITEM,5,24
FITEM,5,-26
FITEM,5,29
FITEM,5,31
FITEM,5,-33
FITEM,5,41
FITEM,5,-43
FITEM,5,46
FITEM,5,48
FITEM,5,-50
FITEM,5,58
FITEM,5,-60
FITEM,5,63
FITEM,5,65
FITEM,5,-67
VSEL,S, , ,P51X
ASEL,ALL
ASLV,R
FLST,5,36,5,ORDE,35
FITEM,5,22
FITEM,5,25
FITEM,5,31
FITEM,5,39
FITEM,5,56
FITEM,5,62
FITEM,5,65
FITEM,5,67
FITEM,5,-68
FITEM,5,77
FITEM,5,80
FITEM,5,86
FITEM,5,94
FITEM,5,111
FITEM,5,117
FITEM,5,120
FITEM,5,122
FITEM,5,-123
FITEM,5,132
FITEM,5,135
FITEM,5,141
FITEM,5,149
FITEM,5,166
FITEM,5,172
FITEM,5,175
FITEM,5,177
FITEM,5,-178
FITEM,5,186
FITEM,5,188
FITEM,5,193
FITEM,5,199
FITEM,5,210
FITEM,5,214
FITEM,5,216
FITEM,5,-218
ASEL,R, , ,P51X
FLST,2,36,5,ORDE,35
FITEM,2,22
FITEM,2,25
FITEM,2,31
FITEM,2,39
FITEM,2,56
FITEM,2,62
FITEM,2,65
FITEM,2,67
FITEM,2,-68
FITEM,2,77
FITEM,2,80
FITEM,2,86
FITEM,2,94
FITEM,2,111
FITEM,2,117
FITEM,2,120
FITEM,2,122
FITEM,2,-123
FITEM,2,132
FITEM,2,135
FITEM,2,141
FITEM,2,149
FITEM,2,166
FITEM,2,172
FITEM,2,175
FITEM,2,177
FITEM,2,-178
FITEM,2,186
FITEM,2,188
FITEM,2,193
FITEM,2,199
FITEM,2,210
FITEM,2,214
FITEM,2,216
FITEM,2,-218
/GO
SFA,P51X, ,FSIN,1,
ALLSEL,ALL
FINISH
/SOL
!*
ANTYPE,2      ! Modal analysis
!*
MSAVE,0
!*
MODOPT,UNSYM,6
EQSLV,FRONT
MXPAND,6, , ,0
LUMPM,0
PSTRES,0
!*
MODOPT,UNSYM,6,1000,15000, ,OFF

Antony21 发表于 2007-6-25 14:56

看来问题好像有点严重，没有人做这方面的吗？

Antony21 发表于 2007-6-27 14:30

要怎么把水对于蜂鸣片频率改变的影响作出来呀？

Antony21 发表于 2007-6-29 10:58

发现一个典型的耦合分析事例，分享给大家 不过有几个问题想跟大家讨论一下


1.在液面上将压力设为0的意义为何？

2.假设我将水箱底部的一小部分挖去，改成压电材料， 我只想看压电材料这个部分的模态及频率，分析的方式是否应该有所改变？

3.我的模型尺寸比例比较小，必须使用SOLID5的元素，分析方式一样吗？
----------------------------------------------------------------------------------------------------------------------------
如何进行盛水水箱的模态分析
这是一个流固耦合模态分析的典型事例，采用ANSYS/MECHANICAL可以完成。处理过程中需要注意以下几个方面的问题：
1、单元的选择；
2、流体材料模式；
3、流固耦合关系的定义；
4、模态提取方法。

一实例，水箱采用SHELL63单元，水箱中的水采用FLUID30单元，以下即为整个流固耦合模态计算的命令流文件：
length=1
width=0.6
height=0.8
/prep7
et,1,63
et,2,30 !选用FLUID30单元，用于流固耦合问题
r,1,0.01
mp,ex,1,2e11
mp,nuxy,1,0.3
mp,dens,1,7800
mp,dens,2,1000 !定义Acoustics材料来描述流体材料-水
mp,sonc,2,1400
mp,mu,2,
!
block,,length,,width,,height
esize,0.1
mshkey,1
!
type,1
mat,1
real,1
asel,u,loc,y,width
amesh,all
alls
!
type,2
mat,2
vmesh,all
fini

/solu
antype,2
modopt,unsym,10 !非对称模态提取方法处理流固耦合问题
eqslv,front
mxpand,10,,,1
nsel,s,loc,x,
nsel,a,loc,x,length
nsel,r,loc,y
d,all,,,,,,ux,uy,uz,
nsel,s,loc,y,width,
d,all,pres,0
alls
asel,u,loc,y,width,
sfa,all,,fsi !定义流固耦合界面
alls
solv
fini

/post1
set,first
plnsol,u,sum,2,1
fini
页: [1]
查看完整版本: 三项耦合的分析--Help都不知道该找哪部分了