ansys压电有限元分析求助
浏览：133865
我目前在做一个压电周期性结构的ansys有限元分析，但是结果总是和理论结果相差很大。程序如下，希望大神有看到的能帮我看一下程序 看看哪里有什么问题。
模型描述：整体结构是一个圆筒，分别是一段压电管一段铝管交替相接，压电管是上下两层，极化方向不一致。构建之后，我发现在管的最左端施加一个径向位移（见代码求解部分），然后考察最右端面上某一节点位移，计算其频率响应函数【此函数为20*log10（abs（Ub/Ua））】，就是这个函数曲线一直和理论解对不上，请大神前辈帮忙分析一下，我的命令流有哪些写的不合适 或者错误的地方
finish
/clear
/prep7
WPSTYLE,,,,,,,,1 !显示工作平面，1表示显示，0表示取消
r1=0.080
r0=0.0805 !定义内外圆半径
r2=0.081
L=0.030 !弹性管=压电管长度
num0=8 !总周期数
num1=2*num0-1
f1=500
f2=4000
ET,1,solid5,3 !小应变
ET,2,SOLID45
ET,3,CIRCU94,1 !电感
ET,4,CIRCU94,0
/COM--------------------------
!定义两个局部坐标来定义两个反向的极化方向
LOCAL,11,1
LOCAL,12,1,0,0,0,180
CYL4,0,0,r1,,r0,,L
CYL4,0,0,r0,,r2,,L
*DO,m,1,num1
WPOFFS,0,0,L !平移工作平面
CYL4,0,0,r1, ,r0, ,L
CYL4,0,0,r0, ,r2, ,L
*ENDDO
WPOFFS,0,0,-num1*L !平移工作平面至初始位置
VGLUE,ALL
!启动圆柱坐标系. (X,Y,Z----R,theta,Z)
!对压电材料进行参数设置
S11=52.0E-12
S12=-18.9E-12
S13=-31.1E-12
S33=67.7E-12
S44=14.0E-12
S66=15.2E-12
/COM, - PIEZOELECTRIC STRAIN COEFFICIENTS, C/N
D31=-921E-12
D33=1981E-12
D15=190E-12
TB,ANEL,1,,,1 ! 弹性柔顺系数 刚度矩阵为0(第一个1是指对材料1赋予)
TBDA,1,S33,S13,S13
TBDA,7,S11,S12
TBDA,12,S11
TBDA,16,S44
TBDA,19,S66
TBDA,21,S44
TB,PIEZ,1,,,1 ! PIEZOELECTRIC STRAIN MATRIX 1 压电系数
TBDA,1,D33
TBDA,4,D31
TBDA,7,D31
TBDA,11,D15
TBDA,18,D15
/COM, - RELATIVE PERMITTIVITY AT CONSTANT STRESS
!
EP11=3307
EP33=1242
MP,PERX,1,EP33
MP,PERY,1,EP11
MP,PERZ,1,EP11
!EMUNIT,EPZRO,8.854e-12 !真空介电常数
!TB,DPER,1,,,1 !常应力下介电常数输入法 仅仅对plane223
!TBDA,1,EP33 !solid 226 solid227可以
!TBDA,2,EP11
!TBDA,3,EP11
!TBLIST,ALL
/COM, - DENSITY, KG/M3
RHO=8040
MP,DENS,1,RHO
!弹性材料参数
MP,DENS,2,2700 !Density
MP,EX,2,7.76e+10 !Young's Modulus
MP,PRXY,2,0.3 !Poisson's ratio
*afun,deg
wprota,,90
vsbw,all
wprota,,,90
vsbw,all
!圆柱管1进行实体单元设置
*Do,J,1,num0
CSYS,1 !下层管
VSLA,,0
ASEL,S,LOC,X,r1,
ASEL,R,LOC,Z,(J-1)*2*L+L,(J-1)*2*L+2*L
VATT,1,,1,11
!上层管
CSYS,1
VSLA,,0 !内外层圆筒同心，vsel已经不能准确选择， !用面选择体的方式来解决 VSLA
ASEL,S,LOC,X,r2
ASEL,R,LOC,Z,(J-1)*2*L+L,(J-1)*2*L+2*L
VATT,1,,1,12
!弹性管
CSYS,1
VSEL,S,LOC,Z,(J-1)*2*L,(J-1)*2*L+L
VATT,2,,2,11
*ENDDO
ALLSEL,ALL
NUMCMP,ALL
nummrg,ALl
/COM---------------------
/COM-----------------------
!划分网格
esize,0.005
MESHAPE,0,3D
MSHKEY,1
VSWEEP,ALL
WPCSYS,-1,0 !将工作平面坐标系旋转与全局坐标系一致
!定义电极：
*DIM,n_bot,array,num0
*DIM,n_mid,array,num0
*DIM,n_top,array,num0
CSYS,1
*DO,g,1,num0
NSEL,S,LOC,X,r1
NSEL,R,LOC,Z,(g-1)*2*L+L,(g-1)*2*L+2*L !定义下层电极
cp,1,volt,all
*get,n_bot(g),node,0,num,min
NSEL,S,LOC,X,r0
NSEL,R,LOC,Z,(g-1)*2*L+L,(g-1)*2*L+2*L !定义中间层电极
cp,next,volt,all
*get,n_mid(g),node,0,num,min
NSEL,S,LOC,X,r2
NSEL,R,LOC,Z,(g-1)*2*L+L,(g-1)*2*L+2*L !定义上层电极
cp,next,volt,all
*get,n_top(g),node,0,num,min
R,3,0.005 !连接电感
type,3 $ real,3
e,n_top(g),n_mid(g)
R,4,1E-06
type,4 $ real,4
e,n_mid(g),n_bot(g)
*ENDDO
ALLSEL,ALL
/prep7
NROTAT,ALL
/SOLU
antype,harmic ! Harmonic analysis
HROPT,FULL ! FULL HARMONIC RESPONSE
*DO,i,f1,f2,10
omega=2*3.1415926*i
UXR=1/omega/omega
CSYS,1
NSEL,S,LOC,X,r0
NSEL,R,LOC,Z,0
D,ALL,UX,UXR
harfrq,i
dmprat,0
ALLSEL,ALL
SAVE
SOLVE ! Solve for resonance frequency
*ENDDO
模型描述：整体结构是一个圆筒，分别是一段压电管一段铝管交替相接，压电管是上下两层，极化方向不一致。构建之后，我发现在管的最左端施加一个径向位移（见代码求解部分），然后考察最右端面上某一节点位移，计算其频率响应函数【此函数为20*log10（abs（Ub/Ua））】，就是这个函数曲线一直和理论解对不上，请大神前辈帮忙分析一下，我的命令流有哪些写的不合适 或者错误的地方
finish
/clear
/prep7
WPSTYLE,,,,,,,,1 !显示工作平面，1表示显示，0表示取消
r1=0.080
r0=0.0805 !定义内外圆半径
r2=0.081
L=0.030 !弹性管=压电管长度
num0=8 !总周期数
num1=2*num0-1
f1=500
f2=4000
ET,1,solid5,3 !小应变
ET,2,SOLID45
ET,3,CIRCU94,1 !电感
ET,4,CIRCU94,0
/COM--------------------------
!定义两个局部坐标来定义两个反向的极化方向
LOCAL,11,1
LOCAL,12,1,0,0,0,180
CYL4,0,0,r1,,r0,,L
CYL4,0,0,r0,,r2,,L
*DO,m,1,num1
WPOFFS,0,0,L !平移工作平面
CYL4,0,0,r1, ,r0, ,L
CYL4,0,0,r0, ,r2, ,L
*ENDDO
WPOFFS,0,0,-num1*L !平移工作平面至初始位置
VGLUE,ALL
!启动圆柱坐标系. (X,Y,Z----R,theta,Z)
!对压电材料进行参数设置
S11=52.0E-12
S12=-18.9E-12
S13=-31.1E-12
S33=67.7E-12
S44=14.0E-12
S66=15.2E-12
/COM, - PIEZOELECTRIC STRAIN COEFFICIENTS, C/N
D31=-921E-12
D33=1981E-12
D15=190E-12
TB,ANEL,1,,,1 ! 弹性柔顺系数 刚度矩阵为0(第一个1是指对材料1赋予)
TBDA,1,S33,S13,S13
TBDA,7,S11,S12
TBDA,12,S11
TBDA,16,S44
TBDA,19,S66
TBDA,21,S44
TB,PIEZ,1,,,1 ! PIEZOELECTRIC STRAIN MATRIX 1 压电系数
TBDA,1,D33
TBDA,4,D31
TBDA,7,D31
TBDA,11,D15
TBDA,18,D15
/COM, - RELATIVE PERMITTIVITY AT CONSTANT STRESS
!
EP11=3307
EP33=1242
MP,PERX,1,EP33
MP,PERY,1,EP11
MP,PERZ,1,EP11
!EMUNIT,EPZRO,8.854e-12 !真空介电常数
!TB,DPER,1,,,1 !常应力下介电常数输入法 仅仅对plane223
!TBDA,1,EP33 !solid 226 solid227可以
!TBDA,2,EP11
!TBDA,3,EP11
!TBLIST,ALL
/COM, - DENSITY, KG/M3
RHO=8040
MP,DENS,1,RHO
!弹性材料参数
MP,DENS,2,2700 !Density
MP,EX,2,7.76e+10 !Young's Modulus
MP,PRXY,2,0.3 !Poisson's ratio
*afun,deg
wprota,,90
vsbw,all
wprota,,,90
vsbw,all
!圆柱管1进行实体单元设置
*Do,J,1,num0
CSYS,1 !下层管
VSLA,,0
ASEL,S,LOC,X,r1,
ASEL,R,LOC,Z,(J-1)*2*L+L,(J-1)*2*L+2*L
VATT,1,,1,11
!上层管
CSYS,1
VSLA,,0 !内外层圆筒同心，vsel已经不能准确选择， !用面选择体的方式来解决 VSLA
ASEL,S,LOC,X,r2
ASEL,R,LOC,Z,(J-1)*2*L+L,(J-1)*2*L+2*L
VATT,1,,1,12
!弹性管
CSYS,1
VSEL,S,LOC,Z,(J-1)*2*L,(J-1)*2*L+L
VATT,2,,2,11
*ENDDO
ALLSEL,ALL
NUMCMP,ALL
nummrg,ALl
/COM---------------------
/COM-----------------------
!划分网格
esize,0.005
MESHAPE,0,3D
MSHKEY,1
VSWEEP,ALL
WPCSYS,-1,0 !将工作平面坐标系旋转与全局坐标系一致
!定义电极：
*DIM,n_bot,array,num0
*DIM,n_mid,array,num0
*DIM,n_top,array,num0
CSYS,1
*DO,g,1,num0
NSEL,S,LOC,X,r1
NSEL,R,LOC,Z,(g-1)*2*L+L,(g-1)*2*L+2*L !定义下层电极
cp,1,volt,all
*get,n_bot(g),node,0,num,min
NSEL,S,LOC,X,r0
NSEL,R,LOC,Z,(g-1)*2*L+L,(g-1)*2*L+2*L !定义中间层电极
cp,next,volt,all
*get,n_mid(g),node,0,num,min
NSEL,S,LOC,X,r2
NSEL,R,LOC,Z,(g-1)*2*L+L,(g-1)*2*L+2*L !定义上层电极
cp,next,volt,all
*get,n_top(g),node,0,num,min
R,3,0.005 !连接电感
type,3 $ real,3
e,n_top(g),n_mid(g)
R,4,1E-06
type,4 $ real,4
e,n_mid(g),n_bot(g)
*ENDDO
ALLSEL,ALL
/prep7
NROTAT,ALL
/SOLU
antype,harmic ! Harmonic analysis
HROPT,FULL ! FULL HARMONIC RESPONSE
*DO,i,f1,f2,10
omega=2*3.1415926*i
UXR=1/omega/omega
CSYS,1
NSEL,S,LOC,X,r0
NSEL,R,LOC,Z,0
D,ALL,UX,UXR
harfrq,i
dmprat,0
ALLSEL,ALL
SAVE
SOLVE ! Solve for resonance frequency
*ENDDO
