论文部分内容阅读
有限元方法是目前结构分析领域应用最为广泛的一种数值方法。自振频率计算作为结构分析中的一个重要内容,可归结为典型的特征值问题,例如:一维波动方程可描述杆结构振动,其特征值为杆结构自振频率;二维波动方程可描述膜结构振动,相应的特征值为膜结构自振频率;三维波动方程可描述声压在空间中的传播,其特征值为声压的振动频率。采用有限元法求解波动方程特征值问题时,质量矩阵的构造方法会直接影响振动频率的计算精度。因而,如何构造可使自振频率有更高收敛阶次及精度的高阶质量矩阵具有重要的研究价值。然而,目前的高阶质量矩阵多限于一维和二维低阶单元,并且二维高阶质量矩阵所得频率计算精度依赖于波的传播方向,不能同时提高各阶频率的收敛阶次及精度。此外,目前对于任意阶单元并没有统一的高阶质量矩阵构造方法。 本文针对波动方程特征值问题有限元分析,提出了一种适用于任意阶单元和维数的新型超收敛有限元分析方法。首先,在Lobatto单元的基础上提出了任意阶一维杆单元高阶质量矩阵的统一构造方法。该方法采用了一种新的质量矩阵构造模式—α矩阵,通过优化参数α可直接构造任意阶Lobatto杆单元的高阶质量矩阵。相较于传统的一致质量矩阵和集中质量矩阵,采用高阶质量矩阵的频率计算精度可以提高2阶。随后,对于二维平面膜结构振动问题和三维波动问题,提出了一种基于积分点的新型超收敛有限元分析方法,该方法中采用的新型积分点是通过对比分析一维高阶质量矩阵确定的。本文提出的新型超收敛有限元分析方法完全消除了超收敛计算中频率的波动方向依赖性,可使任意阶频率的收敛阶次同时提升2阶,同时新型积分点可以通过张量积的形式简便地推广到多维问题,数值实现简便直接。文中通过理论分析和数值计算详细验证了所提波动方程特征值问题新型超收敛有限元分析方法的有效性和计算精度。