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

freund

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

freund

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
  
下面的是后处理了，就不用多说了吧...
祝大家春节快乐^_^

freund

2

freund

3

freund

4

freund

5

freund

6

simwe

谢谢分享

zdhgood

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

你声学采用的是什么单元？圆柱体用的是什么材料？边界条件怎么加的？参考论文是谁的？恳请您的帮助！谢谢你给我带来的希望！急切盼望您的答复！

freund

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

yymaverick

请尽快整理命令流！

yaxing66

已经看到了，可惜没有命令流。希望能将命令流贴上，以供参考。谢谢。

xiriyu

这个我急需，我怎样才能提高我的积分，我想看到它，谢谢！

ShapeMyself

怎么才能自己的积分大于1啊？我想分享这个帖子的内容

zzb001_0

麻烦   楼主  帮我发一份吧    我的分数不够啊     谢谢      zzb001_0@163.com

foxmail3771

你好，freund ！
我想学习你介绍的用ANSYS模拟圆柱体中各种声波（弹性波）传播，可惜我的积分不够，你能否给我发送一份？以求共同进步。谢谢！
发送到：foxmail3771@163.com

xiejin

奇怪,命令流怎么不能运行啊,怎么有如下错误: