波动方程有限元解有关理论.docVIP

1.1控制方程 在经典的线弹性理论当中，线弹性均匀介质的运动方程为 () (1-1) 其中，应力张量和应变张量满足Hooke定律，即满足本构方程: (1-2) 应变张量可通过位移矢量确定，从而给出几何方程: (k,l=1,2,3) (1-3) 各向同性线弹性材料的弹性特征可以用下列材料常数进行描述:弹性模量E,剪切 模量、泊松比、体积模量K和Lame常数等。这些常数中任意两个独立，而其它 几个常数能够用它们间接进行表示。在各向异性的情形下，则需要采用21个相互独立 的弹性常数，进行表示。 介质表面S上的边界条件为 (1-4) (1-5) 由方程(2-1)-(2-3)这三个线弹性动力学的控制方程所描述的线弹性动力问题的初始 条件为 (1-6) (1-7) 其中，，，和是已知量。 由弹性动力学问题解的唯一性定理可知:若弹性体(体积为Y，表面为S)的解能 满足方程(1-1)-式(1-7)，则其位移场、应力场和应变场的解答是唯一的。 将三个控制方程合并，可以得到用位移表示的运动方程: () (1-8) 方程(1-8)所示的就是著名的Navie-Cauchy运动方程。而在各向同性线弹性的情形 下，则有 (1-9) 因此，各向同性材料的本构关系也为表示为 (1-10) 1.2有限元方程 现在我们考虑一个任意形状的封闭区域，并将其用一定形式的网格划分为相应的有 限单元，各单元的单元矩阵可由如下所示的虚功方程求得。 (1-11) 其中，为作用在边界上的外力矩阵; 和为惯性力矩阵，其形式由式(1-12)和(1-13) (1-12) (1-13) 给出。 在单元内任一点的位移，可以用单元形函数N和单元节点位移u近似表示为 (1-14) 由此，可相应的得到 (1-15) (1-16) (1-17) (1-18) (1-19) 其中，D为弹性矩阵，L为微分算子。 将以上各个经过插值后的量(式(1-12)-式(1-19))，代入式(1-11)中，可得 (1-20) 由于虚位移是任意的，因此整理可得 (1-21) 方程(1-21)也可写为有限元方程的标准形式: (1-22) 其中，K为刚度矩阵，C为一致性阻尼矩阵(对于无阻尼问题，该矩阵为零阵)，M为 一致性质量矩阵，F为外力矩阵，分别可表示为 (1-23) (1-24) (1-25) (1-26) 这就表明，对于该封闭区域，通过单元网格划分可以确定其总体刚度矩阵和总体质 量矩阵。因此，只要给定问题的边界条件和初始条件，就能对所列出的线性方程组(1-22) 进行求解，求解得到的结果即为该问题的解。由于材料的粘滞系数难以确定(一般只能 用实验方法测定)，在工程计算中通常采用Rayleigh阻尼法确定阻尼矩阵，即假定阻尼 矩阵为刚度矩阵和质量矩阵的线性组合 (1-27) 特别要注意的是，对于某些无阻尼问题，有限元方程的形式仅为 (1-28) 1