固体中声波的传播

pipinia

大家好，
在论坛上看到了一篇帖子：固体中各种声波（弹性波）传播 http://forum.simwe.com/viewthread.php?tid=155094&extra=page%3D1%26amp%3Bfilter%3Dtype%26amp%3Btypeid%3D142&page=1（很老的贴了 ）。对照该帖中给出的命令流，我自己运行了一遍，但是仿真结果与原帖结果存在差异——
原帖中，某一时间下声波传播结果为：

而我运行的结果为：

与原帖的结果相比，下图中的波动怎么跑到圆柱外面去了（像长了两个小角）？？似乎在圆柱边缘产生了很大的形变。其他时刻的图形结果存在这个问题。而原帖贴的其他的几张图都没有这种畸形的东西出现。请问是什么原因呢？？
以下是原帖给出的命令流：
*SET,_FNCNAME,'myload'                !定义变量名
*DIM,%_FNCNAME%,TABLE,6,10,3        !定义变量维度
/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

sfwei01

这只是显示比例设置的问题，你把显示比例设为1，结果就是你想看到的

pipinia

2# sfwei01
非常感谢！！！的确是显示比例设置的问题。
Plotctrls/style/displacement scaling 里面把设置改过来就好了。再次谢谢

xhd8215

哈哈，终于找到答案了，我也试试。多谢了。

zhsqzzu

各位兄弟姐妹，能否给解释一下这个例子中波的加载问题呀，特别是函数是怎么通过TABLE定义的：
! 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})   
对这个定义不明确，恳请高手指点一下！

amanty

同来向高手学习

bookhere

我想问一个问题，波的传播规律是如此的，能否将波的波动曲线显示出来

yangzhen614

结果怎么查看啊？各位高手们

tzhdacy

后处理1中，plot result 选择节点结果，和自己需要的结果既可

bingbing1478

1# pipinia
不错

yangzhen614

我做了一个立体的
/prep7
et,1,solid65
mp,prxy,1,0.2
mp,ex,1,3e10
mp,1,dens,2600

建立一个立方体
划分网格

对底面施加位移约束0
对侧面上一节点127施加位移约束wo=0.05*sin(2*{PI}/(5*10^-7)*{TIME})
求解
/solu
antype,4
trnopt,full
time,5e-7
d,1,all,0,,4225,1
lswrite,1
ddele,1,all,4225,1
*do,i,2,9
time,5e-7+(i-1)/(8e6)
nsubst,100
d,127,uz,%wo%,,127,1
lswrite,i
*enddo
d,127,uz,0,,127,1
*do,i,10,15
time,25e-7+(i-9)*8/(1e5)
nsubst,100
lswrite,i
*enddo
lssolve,1,15,1
但在127相对的侧面上任意一节点上为什么测不到规则的波动

yangzhen614

接上面
在127上施加位移约束
http://home.simwe.com/space-4989 ... 4156-goto-down.html
127对面上一点测得的波动图
http://home.simwe.com/space-4989 ... 4342-goto-down.html

hitclaire

好像这个table是通过函数加载后得到的，如果一开始选择分段函数加载形式，就会得到这样的table了

yangyangbjiu

各位，谁知道这里面的这个公式 4e-5*sin(2*{PI}*5e6*{TIME})，里面的各个函数代表的是什么呢？
给解释解释吧。。

wmb110119

像高手学习了，自己也弄来试试

馨甜

学习了，谢谢各位，我也都试试

lx24816

终于找到原因了  十分感谢啊