[发明专利]一种基于三层预处理子的三维力学模态仿真模拟方法

权利要求书

1.一种基于三层预处理子的三维力学模态仿真模拟方法，包括以下步骤：

A.将目标飞行器结构进行有限元建模，引入位移边界条件或者应力边界条件建立对应的几何结构模型；

B.采用四面体网格剖分求解域；

C.选择二阶标量叠层基函数，并运用标准的有限元本征分析方法得到目标飞行器结构的有限元本征方程；

D.采用ARPACK软件包中的隐式重启Arnoldi迭代方法求解C步骤所获得的有限元本征方程，获得特征值λ和对应的特征向量α即振幅向量，Arnoldi迭代方法每一步迭代所产生的线性方程采用基于三层预处理子的预处理共轭梯度法求解；

所述三层预处理子的具体构造方法如下：

在隐式重启Arnoldi迭代的过程中需要反复求解以下线性方程

Px＝My_j (1)

其中矩阵P为待预处理的矩阵，M为系统质量矩阵，x为待求中间向量，y_j为迭代向量，采用预处理共轭梯度方法求解(1)式，在预处理共轭梯度法中需要求解如下预处理系统

Nz＝r (2)

其中N为预处理子，z待求预处理向量，r为残差向量；首先把矩阵P按照x,y,z三个方向分块排序，获得如下分块矩阵形式

其中P_xx为x方向的矩阵，P_yy为y方向的矩阵，P_zz为z方向的矩阵，P_xy和P_yx为x方向和y方向的耦合矩阵，P_xz和P_zx为x方向和z方向的耦合矩阵，P_yz和P_zy为y方向和z方向的耦合矩阵，取第一层预处理子的形式为

对(4)式中的块对角矩阵P_xx，P_yy，P_zz采用基于叠层基函数的性质来构造第二层预处理子；首先构造P_xx的第二层预处理子，按照一阶基函数、二阶基函数的排列顺序把P_xx写成如下分块形式

其中P₁₁为一阶基函数形成的矩阵，P₂₂为二阶基函数形成的矩阵，P₁₂和P₂₁为一阶基函数和二阶基函数之间的耦合矩阵，取第二层预处理子为如下形式

矩阵P_yy，P_zz所对应的第二层预处理子的构造方法和矩阵P_xx的第二层预处理子的构造方法一样，分别给他们所对应的第二层预处理子取名为N_y，N_z；

对于(6)式中的块对角矩阵P₁₁和P₂₂采用多波前不完全cholesky分解来构造第三层预处理子，矩阵P₁₁和P₂₂所对应的第三层预处理子分别为N_l和N_h，首先对矩阵P₁₁进行多波前不完全cholesky分解，获得如下第三层预处理子的形式

其中D₀为对角矩阵，P₀为重排序矩阵，L为矩阵P₁₁的cholesky分解下三角矩阵，N_h的构造方法和N_l的构造方法一样，完成三层预处理子的构造；

E.对D步骤获得的特征值和对应特征向量进行后处理获得振动模态频率和对应振动振型。