弹性波仿真分离纵、横波的方法

hadrons

固体中的弹性波含有纵波、横波，也有称为膨胀波、剪切波等等，不同领域叫法不同。
特别是弹性波遇到裂缝、空洞、界面发生散射时，纵波和横波往往会同时产生，如果单单从位移场或应力场是无法区分的。

根据弹性波理论，
位移场的散度 就是膨胀波
位移场的旋度就是剪切波

所以对位移场计算一下散度和旋度就可以分别获得膨胀波和剪切波

我只算2D的情况，这里给一个例子（matlab）
dx=0.1e-3;
dy=0.1e-3;
[x,y]=meshgrid(0:dx:1e-1,-16e-3:dy:-6e-3);
%时间
t=24e-6;
%总位移
vv=postinterp(fem,'disp_acpn',[x(:)';y(:)'],'T',t);
vv=reshape(vv,size(x));
subplot(3,1,1);surf(x,y,vv);axis normal;shading interp;view(0,90);title('所有波')
% return;

%纵波--膨胀波only
ux=postinterp(fem,'ux',[x(:)';y(:)'],'T',t);
vy=postinterp(fem,'vy',[x(:)';y(:)'],'T',t);
ux=reshape(ux,size(x));
vy=reshape(vy,size(x));
subplot(3,1,2);surf(x,y,ux+vy);shading interp;view(0,90);title('膨胀波')

%横波--旋变only
uy=postinterp(fem,'uy',[x(:)';y(:)'],'T',t);
vx=postinterp(fem,'vx',[x(:)';y(:)'],'T',t);
uy=reshape(uy,size(x));
vx=reshape(vx,size(x));
subplot(3,1,3);surf(x,y,uy-vx);shading interp;view(0,90);title('切变波')

附件中容易看到散射波是剪切波（横波）

xhfzjlx

很好的思路，很受启发。
此处，fem是从comsol到过来的数据吗

hadrons

2# xhfzjlx
是的。
fem导入到matlab中

sunnywell

这个相当不错，受益匪浅

bmezhou

未考虑耦合波吧。

burger

学到东西了谢谢lz

qianping

能不能不用在matlab中二次开发，直接在后处理里面搞定。

shenchenghit

不知道能不能将弯曲波和横波，纵波分离开？比如Euler梁单元，怎么得到纯弯曲模式的固有频率？

duanqiao86

学习了，受启发，谢谢。