[发明专利]基于Householder变换的无约束结构静力分析方法

摘要

本发明公开了基于Householder变换的无约束结构静力分析方法，其包括：步骤1、构建结构的刚体位移模式矩阵X；步骤2、根据刚体位移模式矩阵X构造6个相应的n×n阶Householder矩阵Pi；步骤3、利用步骤2中构造的Householder矩阵Pi对结构原始刚度矩阵做正交相似变换，得到去除刚体模态后的结构的刚度矩阵Kp；步骤4、采用修正共轭梯度法，求解用Householder矩阵去除结构总体刚体位移模态后的刚度方程KPUP=FP。本发明通过去除结构的刚体位移模态，并控制共轭梯度法的指定最小误差阈值，精确求解结构响应；计算实施步骤简洁，不需要修改目前通用的有限元计算框架；采用修正共轭梯度法使得结构总体刚度矩阵的稀疏特性能够得到很好的利用，分步求解使得整体求解过程占用空间少、计算效率高。

主权项

基于Householder变换的无约束结构静力分析方法，其特征在于，其包括：步骤1、构建结构的刚体位移模式矩阵H：H=[h1 h2 h3 h4 h5 h6]    （1）其中，h1～h3为结构的平动刚体位移模式，h4～h6为结构的转动刚体位移模式，h1～h6均为n维列向量，n为结构整体自由度；步骤2、根据h1～h6构造6个相应的n×n阶Householder矩阵Pi，i=1，2……6；设结构离散化后计算得到的未指定位移边界条件的有限元整体刚度阵为K，则此刚度阵K具有6重零特征值，且具有6个对应于零特征值的正交特征向量，这6个特征向量张成的空间为刚度阵K的零空间，这个零空间与6个结构刚体位移模式张成的空间一致，即满足K*hi=0    （2）所述步骤2包括以下步骤：步骤21、首先将平动刚体位移模式h1单位化，且与单位向量e1一起用来构造Householder矩阵P1，使得：P1*h1=e1    （3）所述Householder矩阵P1为对称正交矩阵；使用Householder矩阵P1对刚度阵K进行相似变换，将刚度阵K的第一行和第一列均化为0，即：K1=P1*K*P1    （4）其中，K1是刚度阵K的第一行和第一列均为0的相似刚度矩阵；步骤22、构造Householder矩阵P2，使用Householder矩阵P2对相似刚度矩阵K1进行相似变换，可使相似刚度矩阵K1的前两行和前两列均为0，即：K2=P2*K1*P2    （5）其中，K2是相似刚度矩阵K1和刚度阵K的前两行和前两列均为0的相似刚度矩阵；构造Householder矩阵P2的方法是：根据公式（2）以及正交对称矩阵特性可知：K*hi=0    （6）P1*P1=In    （7）其中，In为n×n阶单位矩阵；由公式（6）、（7）可得：K1*P1*h2=0    （8）令：h2′=P1*h2    （9）h2′为新构造的平动刚体位移模式，由于相似刚度矩阵K1的第一列元素均为零，且根据特征向量两两相互正交的特点，将h2′的第一个元素赋值为零，即：h2′(1)=0    （10）将新构造的平动刚体位移模式h2′单位化与单位向量e2一起构造Householder矩阵P2：P2*h2′=e2    （11）步骤23、根据步骤22中构造Householder矩阵P2的方法逐个对平动刚体位移模式h3以及转动刚体位移模式h4～h6进行操作，构造新的平动刚体位移模式h3′以及新的转动刚体位移模式h4′～h6′，进一步构造Householder矩阵P3～P6；步骤3、利用步骤2中构造的Householder矩阵Pi对结构原始刚度矩阵做正交相似变换，得到去除刚体模态后的结构的刚度矩阵Kp：Kp=P6P5P4P3P2P1KP1P2P3P4P5P6    （12）简化得：Kp=PTKP    （13）其中，P=P1P2P3P4P5P6；步骤4、通过有限元离散并用常规方法得到结构静力问题分析的线性方程：KU=F    （14）其中K为刚度阵，U为待求解的位移向量，F为外载右端项；通过Householder矩阵的正交相似变换去除结构总体刚体位移模态后上述方程（14）为：KPUP=FP    （15）其中，KP=PTKP，UP=PTU，FP=PTF，P=P1P2P3P4P5P6。