pengweicai 发表于 2007-5-27 15:08

【转贴】圆柱体中各种声波（弹性波）传播

【分享】圆柱体中各种声波（弹性波）传播

原来一直有种映像，认为ANSYS的声波（弹性波）方面功能很弱，甚至认为不能进行反射波的分析计算。今天按照一篇博士论文（采用FEM, 用MATLAB编程），做了类似的分析，结果发现分析结果和原作者计算的很相近，观察的现象完全一致。让我对ANSYS的功能有了更深的理解，在这我也特别要感谢原作者。
    下面贴几张图形，图形中显示的有纵波，横波，瑞立波，PP波，PS波等，与大家共享，希望能对大家有所启示。

回复： 【分享】圆柱体中各种声波（弹性波）传播

Sorry，太忙了，好几个朋友发邮件询问我命令流的事。现在手头的事，刚有了眉目，今天也要回家了，抽空稍微整理了一下。

/units,si
/prep7
et,1,plane82,,,2      
mp,dens,1,2700
mp,ex,1,70e9
mp,nuxy,1,0.34
pcirc,0,4e-3,0,180
esize,4e-5
mshape,1,2d
mshkeyk,0
amesh,1
finish

/solu
antype,trans
nsel,s,loc,y,0
dsym,symm,y
allsel
nsel,s,node,,2

! 下面为分段载荷函数，如果
! 采用菜单命令，会更简单
! 加载函数名称为 ：myload，具体为
!4e-5*sin(2*{PI}*5e6*{TIME})， 当 t<2e-7
!                     0                     ， 当 t>2e-7
! Begin of equation: {TIME}
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, 99, 0, 1, 1, 0, 0   
*SET,%_FNCNAME%(0,2,1),   0
*SET,%_FNCNAME%(0,3,1),   0
*SET,%_FNCNAME%(0,4,1),   0
*SET,%_FNCNAME%(0,5,1),   0
*SET,%_FNCNAME%(0,6,1),   0
*SET,%_FNCNAME%(0,7,1),   0
*SET,%_FNCNAME%(0,8,1),   0
*SET,%_FNCNAME%(0,9,1),   0
*SET,%_FNCNAME%(0,10,1),   0
! End of equation: {TIME}   
! Begin of equation: 4e-5*sin(2*{PI}*5e6*{TIME})
*SET,%_FNCNAME%(0,0,2), 2.0e-7, -999
*SET,%_FNCNAME%(2,0,2), 0.0
*SET,%_FNCNAME%(3,0,2), 0.0
*SET,%_FNCNAME%(4,0,2), 0.0
*SET,%_FNCNAME%(5,0,2), 0.0
*SET,%_FNCNAME%(6,0,2), 0.0
*SET,%_FNCNAME%(0,1,2), 1.0, -1, 0, 2, 0, 0, 0   
*SET,%_FNCNAME%(0,2,2), 0.0, -2, 0, 3.14159265358979310, 0, 0, -1   
*SET,%_FNCNAME%(0,3,2),   0, -3, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,4,2), 0.0, -1, 0, 5e6, 0, 0, -3   
*SET,%_FNCNAME%(0,5,2), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,6,2), 0.0, -1, 0, 1, -2, 3, 1
*SET,%_FNCNAME%(0,7,2), 0.0, -1, 9, 1, -1, 0, 0
*SET,%_FNCNAME%(0,8,2), 0.0, -2, 0, 4e-5, 0, 0, -1   
*SET,%_FNCNAME%(0,9,2), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,10,2), 0.0, 99, 0, 1, -3, 0, 0
! End of equation: 4e-5*sin(2*{PI}*5e6*{TIME})   
! Begin of equation: 0   
*SET,%_FNCNAME%(0,0,3), 6e-6, -999   
*SET,%_FNCNAME%(2,0,3), 0.0
*SET,%_FNCNAME%(3,0,3), 0.0
*SET,%_FNCNAME%(4,0,3), 0.0
*SET,%_FNCNAME%(5,0,3), 0.0
*SET,%_FNCNAME%(6,0,3), 0.0
*SET,%_FNCNAME%(0,1,3), 1.0, 99, 0, 0, 0, 0, 0   
*SET,%_FNCNAME%(0,2,3),   0
*SET,%_FNCNAME%(0,3,3),   0
*SET,%_FNCNAME%(0,4,3),   0
*SET,%_FNCNAME%(0,5,3),   0
*SET,%_FNCNAME%(0,6,3),   0
*SET,%_FNCNAME%(0,7,3),   0
*SET,%_FNCNAME%(0,8,3),   0
*SET,%_FNCNAME%(0,9,3),   0
*SET,%_FNCNAME%(0,10,3),   0
! End of equation: 0

d,all,ux,%myload%
allsel
time,2.0e-8
tintp,,,,1.0
auto,on
deltim,2e-9,2e-9,5e-9
outres,nsol,1
solve

下面的是后处理了，就不用多说了吧...

pengweicai 发表于 2007-5-27 15:15

回复： 【分享】圆柱体中各种声波（弹性波）传播

以上分析实际很简单，只不过对网格划分以及时间精度要求很高，具体可参考吴先梅博士论文（同济大学，2000）。这段时间很忙，我是用GUI和命令流结合做的，以后如果有必要我会将命令流整理，贴出的。再次对作者表示忠心的感谢，也祝大家共同进步。


原贴位置
http://www.simwe.com/forum/viewthread.php?tid=155094&extra=page%3D1%26amp%3Bfilter%3Dtype%26amp%3Btypeid%3D142&page=1

sophia_shi 发表于 2007-5-29 21:07

呀，吴先梅老师是俺们实验室的，嘿嘿。

pengweicai 发表于 2007-5-30 09:20

原帖由 sophia_shi 于 2007-5-29 21:07 发表 http://www.chinavib.com/forum/images/common/back.gif
呀，吴先梅老师是俺们实验室的，嘿嘿。

   不知道你是做那方面的，有什么好东西欢迎分享啊。

shitty 发表于 2007-6-8 20:22

我直接用这个命令流做了一下

但是结果是一片兰色没有任何变化

请问这个命令流 需要改什么设置么？

要用什么版本的ansys啊我8 和10 都试了一样的结果

zhangjuhua 发表于 2007-7-4 18:14

这个流的函数加载错误

这个流的函数加载错误,具体看仿真论坛上搜索有改进过的,有结果,不过结果没上面的好看,我运行过

thatsky 发表于 2008-10-31 14:42

这个流函数应该怎么加呢？
加载函数与结果的联系是什么样的？
哪位高手知道。。。

yaotao 发表于 2008-12-11 20:13

运行后的波怎么跑到圆板外边？

zhangjiae 发表于 2009-4-12 15:21

我也在研究，不知道是不是其他软件会比较好

syq4223 发表于 2009-5-11 19:56

回复 6楼 zhangjuhua 的帖子

能不能把你具体的命令流发上来？
大家共享一下！

sgive 发表于 2009-7-16 08:38

我现在在做超声波模具分析，这里所讨论的函数加载我不懂，希望各位给与指点！谢谢

jxyfiting 发表于 2011-5-18 16:22

楼主很是厉害啊，我是刚入门，学习中

sunyang1010 发表于 2011-9-4 08:53

这是个好贴，顶

lx24816 发表于 2013-9-21 22:15

{:{39}:}请教楼主你的后处理照片怎么这么清晰哦
页: [1]
查看完整版本: 【转贴】圆柱体中各种声波（弹性波）传播