hy1966老师,我想用ANSYS做压电陶瓷振动时在混凝土内部的声场分布或者应力的分布。我的压电陶瓷用的SOLID5,混凝土用的SOLID45,并把压电陶瓷和混凝土面接触的地方做了粘接,现在已经建模和划分完网格了,但是不知道需要用哪种分析方法,怎么加载载荷和约束,希望老师指点。代码如下:
!压电陶瓷圆片直径12.5mm,厚2mm
!混凝土模块长100mm,宽100mm,高50mm
fini
/clear
/prep7
k,1
k,2,,,50
k,3,,,50+2
l,1,2
l,2,3
k,
Blc5,0,0,100,100
circle,1,6.25
csys,1
lsel,s,loc,z,0
lsel,r,loc,x,0,20
lgen,2,all,,,,45,,,,1
allsel
lplot
csys,0
LSTR, 8, 10
LSTR, 10, 1
LSTR, 1, 12
LSTR, 12, 6
LSTR, 5, 11
LSTR, 11, 1
LSTR, 9, 1
LSTR, 9, 7
allsel
lsel,u,line,,1
lsel,u,line,,2
asbl,all,all
CM,A1,area
Vdrag,A1,,,,,,1
csys,1
asel,s,loc,z,50
asel,r,loc,x,0,20
Vdrag,All,,,,,,2
csys,0
allsel
vglue,all
Vsel,s,loc,z,0,50
CM,V2,volume
Vsel,s,loc,z,50,52
CM,V1,volume
allsel
!赋属性
ET,1,solid5 !压电陶瓷
KEYOPT,1,1,3
KEYOPT,1,3,0
KEYOPT,1,5,2
ET,2,solid45 !混凝土
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,7500
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,PERX,1,,370
MPDATA,PERY,1,,370
MPDATA,PERZ,1,,635
TB,ANEL,1,1,21,0
TBTEMP,0
TBDATA,,1.39e+011,7.78e+010,7.43e+010,0,0,0
TBDATA,,1.39e+011,7.43e+010,0,0,0,1.15e+011
TBDATA,,0,0,0,3.06e+010,0,0
TBDATA,,2.56e+010,0,2.56e+010,,,
TB,PIEZ,1,,,0
TBMODIF,1,1,0
TBMODIF,1,2,0
TBMODIF,1,3,-5.2
TBMODIF,2,1,0
TBMODIF,2,2,0
TBMODIF,2,3,-5.2
TBMODIF,3,1,0
TBMODIF,3,2,0
TBMODIF,3,3,15.1
TBMODIF,4,1,0
TBMODIF,4,2,0
TBMODIF,4,3,0
TBMODIF,5,1,0
TBMODIF,5,2,12.7
TBMODIF,5,3,0
TBMODIF,6,1,12.7
TBMODIF,6,2,0
TBMODIF,6,3,0
TBDE,PIEZ,1
TB,PIEZ,1,,,0
TBMODIF,1,1,0
TBMODIF,1,2,0
TBMODIF,1,3,-5.2
TBMODIF,2,1,0
TBMODIF,2,2,0
TBMODIF,2,3,-5.2
TBMODIF,3,1,0
TBMODIF,3,2,0
TBMODIF,3,3,15.1
TBMODIF,4,1,0
TBMODIF,4,2,0
TBMODIF,4,3,0
TBMODIF,5,1,0
TBMODIF,5,2,12.7
TBMODIF,5,3,0
TBMODIF,6,1,12.7
TBMODIF,6,2,0
TBMODIF,6,3,0
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,24000
MPDATA,PRXY,2,,0.2
TBDE,ANEL,1,,,
TB,ANEL,1,1,21,0
TBTEMP,0
TBDATA,,1.39E+011,7.78E+010,7.43E+010,0,0,0
TBDATA,,1.39E+011,7.43E+010,0,0,0,1.15E+011
TBDATA,,0,0,0,3.06E+010,0,0
TBDATA,,2.56E+010,0,2.56E+010,,,
TB,CONC,2,1,9,
TBTEMP,0
TBDATA,,0.35,1,3.1125,-1,0,0
TBDATA,,0,0,0,,,
Vsel,s,Volume,,V1
VATT,1,,1,0
! the rubber matrix
Vsel,s,volume,,V2
VATT,2,,2,0
allsel
!网格划分
MSHAPE,0,3D
MSHKEY,1
csys,0
Lsel,all
Lsel,u,Line,,1
Lsel,u,Line,,2
Lsel,r,loc,z,1,49
Lesize,all,,,10
Lsel,all
Lsel,u,Line,,1
Lsel,u,Line,,2
Lsel,r,loc,z,50.1,51.9
Lesize,all,,,2
Lsel,s,Loc,z,0
Lsel,a,Loc,z,50
Lsel,a,Loc,z,52
Lesize,all,,,10
!two.划分网格
allsel
Vsweep,all
!************************************几何模型 网格划分完成************************
!加载
asel,s,loc,z,0
da,all,all
csys,1
nsel,s,loc,x,0
nsel,r,loc,z,50
d,all,Ux,0
d,all,Uy,0
nsel,s,loc,x,0
nsel,r,loc,z,52
d,all,Ux,0
d,all,Uy,0
csys,1
asel,s,loc,z,52
DA,all,VOLT,10
asel,s,loc,z,50
asel,r,loc,x,0,10
DA,all,VOLT,0
allsel
finish
!求解
/solu
solve
不知道是哪里的问题,这样仿真出来的结果看不到应力的分布情况。我感觉是约束和载荷有问题,望老师指教
|