关于弹性波的FDTD算法的问题

liujunhit

file:///C:/Documents%20and%20Settings/Administrator/桌面
这是我划分的网格的简要模型
格数：5×5
各个分量的坐标范围：
X方向速度：vx(1:6,1:6)
Y方向速度：vy(1:5,1:5)
应力：txx(1:5,1:6)
tyy(1:5,1:6)
txy(1:6,1:5)
问题：
1、不知道这样合不合理？
2、我是用pml边界条件，但不知道pml最外层界面处的vx,txx,txy,tyy该怎样赋值？

1）vx在四个pml最外层界面处的值都为零吗？还是只在垂直于x轴的两个界面上设为零？但不知在垂直于y轴的两个界面上该如何设置？
2）txx,txy,tyy,在四个pml最外层界面处该怎么设置？
3、关于加载平面波的问题：可不可以直接将激励源值赋给vx的一列呢？：例如vx(2,1:6)=pulse; 是否需要用电磁场中加平面波的方法呢？
希望能有高手指点迷津，小弟在这里先谢过了