[发明专利]通过修正奇异位移分量消除褶皱计算时刚度矩阵奇异性的方法

权利要求书

1.通过修正奇异位移分量消除褶皱计算时刚度矩阵奇异性的方法，其特征是：所述方法的过程为：

步骤一：在皱曲特征矢量空间中，引入修正参数修正奇异刚度矩阵，获取非奇异刚度矩阵形式；

步骤二：结合非奇异刚度矩阵形式获取真实位移解；

步骤三：判断真实位移解的奇异性，并获取奇异位移分量；

步骤四：在皱曲临界点处引入皱曲模态等式、皱曲临界点处皱曲模态对应的载荷消除后的关系等式及奇异位移解存在条件等式，并利用所述的三个等式修正奇异位移分量，获得非奇异位移解，通过修正奇异位移解进而消除皱褶计算时刚度矩阵的奇异性。

2.根据权利要求1所述的通过修正奇异位移分量消除褶皱计算时刚度矩阵奇异性的方法，其特征在于：步骤一中所述的引入修正参数修正奇异刚度矩阵，获取非奇异刚度矩阵形式的过程具体为：

特征矢量空间中第一阶特征矢量的列向量形式为：

                φ＝(1，0，0，…，0)^T

所述的特征矢量空间中第一阶特征矢量φ为第1位为1，其余位为0的列向量，

引入一个修正参数β＞0，并结合特征矢量空间中第一阶特征矢量φ得到修正后的非奇异刚度矩阵形式为：

KT*=KT+Kφ=KT+βφφT]]>

其中，K_T为奇异刚度矩阵，K_φ为引入的修正矩阵，φ^T为特征矢量空间中第一阶特征矢量φ的转置。

3.根据权利要求1所述的通过修正奇异位移分量消除褶皱计算时刚度矩阵奇异性的方法，其特征在于：步骤二中所述的结合非奇异刚度矩阵形式获取真实位移解的过程具体为：

假定真实位移解为u，将其与修正后的非奇异刚度矩阵相乘，得到：

KT*u=KTu+β(φTu)φ=q+β(φTu)φ]]>

其中，q为对应切线刚度矩阵的载荷，

此时无奇异，得到真实位移解u：

u=uq+βφTuq1-βφTuφuφ]]>

其中，uφ=(KT*)-1φ,]]>uq=(KT*)-1q.]]>