lpc8584 发表于 2013-1-8 16:20:37

请hy1966老师帮忙看一下这个仿真问题

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

不知道是哪里的问题，这样仿真出来的结果看不到应力的分布情况。我感觉是约束和载荷有问题，望老师指教

lpc8584 发表于 2013-1-9 11:38:51

hy1966老师或者其他大侠们，帮我看看吧。。。。困扰好几天了，毕业设计要用的，混凝土内的声场是不是都不能仿啊？我看好多帖子都是做的水或者空气，没有做混凝土的。。。。另外代码有问题么？我感觉问题应该不小，，，，，麻烦各位高人指点啊

lpc8584 发表于 2013-1-9 17:52:23

各位高手，混凝土内部的声压和应力是不是不能仿啊？谁帮我看看程序的问题啊？感激不尽啊
页: [1]
查看完整版本: 请hy1966老师帮忙看一下这个仿真问题