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

pengweicai

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

原来一直有种映像，认为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

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

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


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

sophia_shi

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

pengweicai

原帖由 sophia_shi 于 2007-5-29 21:07 发表

                               
登录/注册后可看大图

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

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

shitty

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

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

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

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

zhangjuhua

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

thatsky

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

yaotao

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

zhangjiae

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

syq4223

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

sgive

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

jxyfiting

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

sunyang1010

这是个好贴，顶

lx24816

请教楼主  你的后处理照片怎么这么清晰哦