柱坐标三维有限差分.docx

参考文献 柱坐标三维有限差分 速度应力方程组 圆柱坐标中一阶速度应力方程组， ( STYLEREF 1 \s 11) 其中，为速度，为应力。是体积外力分量，是应力变化速度的体积源。具体在下面有加载方法。 吸收边界条件 完美匹配层PML 首先利用波动方程的平面波解来阐述PML层的作用。以振动速度为例，在二维情况下，其波动方程平面波解的形式一般为: 以x轴为例，在吸收层内引入复拉伸因子： 其中，的脚标表示x方向，为角频率，为衰减函数： 其中，。 等式两侧对x积分，得到： 在吸收层内用代替x，则式变为： 由此可见，在PML层内入射波将在x方向呈指数衰减。 接下来考虑PML方法在波动方程上的实现，以式中的关系式为例，其频域表达形式如下： 其中，分别为的频域形式。在PML层内引入了复拉伸因子，因此需要将上式中的空间偏微分算子分别用来代替，则式变为： 分裂式完美匹配层SPML 使用PML将该频域表达式变换回时域中，必然会涉及到卷积运算。为了避免卷积运算，SPML方法将一阶波动方程式中的每个速度分量和应力分量进一步分裂成与坐标轴方向平行和垂直的三部分。例如Vx可以表示为：Vx = Vxx +Vxz，则方程可以展开成以下形式： 即： 可以十分容易的将（4-9）式变换到时域而不需卷积运算： 不同位置的PML吸收层的参数设置 在FDTD的实现时，除了在几个PML的角