CN103902764B - 基于Householder变换的无约束结构静力分析方法 - Google Patents

基于Householder变换的无约束结构静力分析方法 Download PDF

Info

Publication number
CN103902764B
CN103902764B CN201410091321.7A CN201410091321A CN103902764B CN 103902764 B CN103902764 B CN 103902764B CN 201410091321 A CN201410091321 A CN 201410091321A CN 103902764 B CN103902764 B CN 103902764B
Authority
CN
China
Prior art keywords
matrix
householder
rigid body
stiffness matrix
body mode
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201410091321.7A
Other languages
English (en)
Other versions
CN103902764A (zh
Inventor
丁桦
何志年
张俊波
李庭婷
蔡智华
夏昌奇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Industry Technology Guangzhou of CAS
Original Assignee
Institute of Industry Technology Guangzhou of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Industry Technology Guangzhou of CAS filed Critical Institute of Industry Technology Guangzhou of CAS
Priority to CN201410091321.7A priority Critical patent/CN103902764B/zh
Publication of CN103902764A publication Critical patent/CN103902764A/zh
Application granted granted Critical
Publication of CN103902764B publication Critical patent/CN103902764B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Complex Calculations (AREA)

Abstract

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

Description

基于Householder变换的无约束结构静力分析方法
技术领域
本发明涉及有限元的结构分析,具体涉及一种无约束状态下的结构力学问题的有限元分析方法。
背景技术
通常做线性静力分析需要保证结构没有刚体位移,否则求解器没有办法计算。但是在空间中很多问题没有足够的约束,如飞机在飞行时,轮船在航行时,或者货品在起重机上吊装时,要想计算结构上的应力分布,需要采用特殊的无约束结构静力分析方法。现有的分析方法有以下几类。
1)惯性释放法,其基本思路是用结构的惯性(质量)力来平衡外力。尽管结构没有约束,分析时仍假设其处于一种“静态”的平衡状态。采用惯性释放功能进行静力分析时,需要对一个节点进行6个自由度的约束(虚支座)。针对该支座,程序首先计算在外力作用下每个节点在每个方向上的加速度,然后将加速度转化为惯性力反向施加到每个节点上,由此构造一个平衡的力系(支座反力等于零)。求解得到的位移描述所有节点相对于该支座的相对运动。
2)动力松弛法,动力松弛法是用显式动力学原理将结构的静力学问题转化为动力学问题近似求解线性或非线性系统平衡状态的数值方法,在计算时添加虚拟的质量和阻尼,构成质量-刚度-阻尼的平衡体系,添加的阻尼应足够大,以使得结构的动能在动力松弛分析完成后足够小。
随着有限元技术的推广,各种各样的无约束问题涌现出来,在用传统有限元法对其进行求解时,会遇到刚度阵奇异的问题。而现有的无约束问题求解法如惯性释放法需要引入虚支座的概念,由于虚支座不是真实的位移约束,与真实情况相比会产生较大的误差,且需要求解方程,效率有时偏低;动力松弛法由于采用显式分析方法,因此效率高但精度低。
对于有限元方程组的求解,通常结构有限元分析中形成的总体刚度阵都具有稀疏质,对于这类刚度阵对应的线性方程组非常适合采用共轭梯度法进行求解,共轭梯度法不仅是解决大型线性方程组最有用的方法之一,也是解大型非线性最优化最有效的优化方法之一。在各种优化方法中,共轭梯度法是非常重要的一种。其优点是所需存储量小,具有有限步收敛性,稳定性高,而且不需要任何外来参数。
发明内容
针对现有技术存在的问题,本发明的首要目的是提供一种基于Householder变换的无约束结构静力分析方法,该方法以空间无约束模型为基础,根据Householder变换理论,用结构的6个的刚体位移模式构造6个Householder矩阵,用此Householder矩阵去除结构的刚体模态,并采用修正后的共轭梯度法求解去刚体模态后的刚度方程。
为实现以上目的,本发明采取了以下的技术方案:
基于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(e1为第一个元素为1,其它元素都为0的单位向量)一起用来构造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
步骤1包括:假设将结构整体沿x方向平移一个单位位移,则在物体内不会产生任何应变,而此时结构的刚体位移模式为只在第6j+1个自由度上有位移值1其它都为0,具体可写为:
h1=[1 0 0 0 0 0 1 0 0 0 0 0 … 1 0 0 0 0 0]T (16)
同理,将结构整体分别沿y方向和z方向平移一个单位位移,构造平动刚体位移模式h2和h3
假设将结构上的节点m绕z轴旋转一个微小角度θ,此时m点x方向的位移为:
u = x m ′ - x m = - δ sin α = - r sin αθ = - y m θ - - - ( 17 )
同理,y方向的位移为:
v=xmθ (18)
z方向的位移为零:
w=0 (19)
此时m点的位移模式可写为:
h6=[-ym θ xmθ 0 0 0 1θ]T (20)
在上式中约去θ,并将其扩展到整个结构,即结构整体绕z轴旋转一个微小角度θ,而在物体内不产生任何应变,此时结构的刚体位移模式可写为:
h6=[-y1 x1 0 0 0 1 -y2 x2 0 0 0 1 … -yn/6 xn/6 0 0 0 1]T (21)
其中,x1~xn/6是节点1到节点n/6的初始x轴坐标,y1~yn/6是节点1到节点n/6的初始y轴坐标,1≤m≤(n/6);
同理,将结构整体分别绕x轴和y轴旋转一个微小角度θ,构造转动刚体位移模式h4和h5
Householder变换理论可以变换表述为:对于任意给定的两个n阶单位向量a和b,定义向量v满足:
c=a+b (22)
v=c/|c| (23)
其中|c|为向量c的长度。
Householder矩阵可以构造为
Pi=2vivi T-In (24)
其中vi为归一化后的n维Householder向量。矩阵P为对称正交矩阵,满足如下三个条件:
Pi=Pi T (25)
Pi=Pi -1 (26)
Pi*a=b (27)
在本方案中,Householder矩阵Pi均为对称正交矩阵,而由这些Householder矩阵Pi组成的Householder矩阵P则只满足正交性,不满足对称性。
步骤5、采用修正共轭梯度法求解线性方程组(15):
去除刚体模态后的结构的刚度矩阵Kp中前六行的元素均为零,因此相应的要把变换后的外载右端项的前六个元素设置为零,即:
FP(1)=FP(2)=...=FP(6)=0 (28)
为避免直接计算去除刚体模态后的结构的刚度矩阵Kp,而造成占用大量空间,并避开共轭梯度法求解的过程中计算矩阵向量乘法的时候需要花费大量的计算时间,对常规的共轭梯度法进行相应的修正,以便充分发挥共轭梯度法的优点。
对于占用空间大的问题,可以仅需存储具有稀疏特性的刚度阵K和Householder变换矩阵Pi以代替去除刚体模态后的结构的刚度矩阵Kp,由Householder变换的性质可知
Pi=2vivi T-In (29)
更进一步地也不需要存储6个Householder变换矩阵,仅需要存储Householder变换矩阵Pi对应的向量vi,避免存储满阵Kp占用大量空间(其中,刚度阵K具有稀疏特性,其为稀疏矩阵,而Kp为满阵)。这样的做法可以减少大量不必要的计算机存储空间。
另一方面,在共轭梯度法中计算矩阵向量乘法,可以采用分步乘积的方式,具体形式如下:
Kp*pk=PT*K*P*pk (30)
其中,pk为共轭梯度法中的修正方向。从上式可以看出,上式的矩阵向量乘法可以分解为三步进行:
sk=P*pk (31)
sk=K*sk (32)
sk=Kp*pk=PT*sk (33)
从上面的公式(32)可以看到,做矩阵向量乘法运算的时候,总体刚度矩阵K的稀疏特性可以大大提高计算效率。
更进一步,计算PT*sk时也可以均分解成6次乘法进行,每一次均对向量前乘一个Householder变换矩阵,由Householder变换矩阵的特性可知,对于任意Householder矩阵P和任一向量g,其乘积可以表示为
P*g=(2vvT-In)g=2v(vT g)-g (34)
由上式可以看出Householder矩阵P与向量做乘积运算时,实际上可以用一次点积、一次标量向量乘法和一次向量加减法运算代替,这无疑会大大减少计算量,特别是对于大规模、节点数庞大的结构分析问题,上述分析方法可以大大提高计算效率。
本发明与现有技术相比,具有如下优点:本发明可以在无需构造任何“虚”边界条件的情况下,通过去除结构的刚体位移模态,并控制共轭梯度法的指定最小误差阈值,精确求解结构响应;计算实施步骤简洁,不需要修改目前通用的有限元计算框架;采用修正共轭梯度法使得结构总体刚度矩阵的稀疏特性能够得到很好的利用,分步求解使得整体求解过程占用空间少、计算效率高。
附图说明
附图1为求旋转刚体位移模式的说明图;
附图2为实施例的壳单元示意图;
附图3为实施例的位移误差曲线。
具体实施方式
下面结合附图和具体实施方式对本发明的内容做进一步详细说明。
实施例:
本实施例为一个如图2所示的四边形壳单元构成的平面问题,一种材料。Abaqus单元类型:S4R。单元尺寸为1*1m,厚度为0.01m,材料参数为:弹性模量E=1.0,泊松比μ=0.3。
在Abaqus中,对上述单元的左端1、3节点施加固定边界条件,右端2、4节点施加向右的拉力。计算得出单元各节点的支反力和支反力矩,作为本方法求解的与Abaqus中有约束的静力等效的无约束问题的外载右端项F。
根据壳单元的位移公式及壳的平面应变和弯曲应变为0,写出上述单元的刚体位移,原理请参照图1所示。假设将结构整体沿x方向平移一个单位位移,则在物体内不会产生任何应变,而此时结构的刚体位移模式为只在第6j+1个自由度上有位移值1其它都为0,具体可写为:
h1=[1 0 0 0 0 0 1 0 0 0 0 0 … 1 0 0 0 0 0]T (35)
同理,将结构整体分别沿y方向和z方向平移一个单位位移,构造平动刚体位移模式h2和h3
假设将结构上的节点m绕z轴旋转一个微小角度θ,此时m点x方向的位移为:
u = x m ′ - x m = - δ sin α = - r sin αθ = - y m θ - - - ( 36 )
同理,y方向的位移为:
v=xmθ (37)
z方向的位移为零:
w=0 (38)
此时m点的位移模式可写为:
h6=[-ymθ xmθ 0 0 0 1θ]T (39)
在上式中约去θ,并将其扩展到整个结构,即结构整体绕z轴旋转一个微小角度θ,而在物体内不产生任何应变,此时结构的刚体位移模式可写为:
h6=[-y1 x1 0 0 0 1 -y2 x2 0 0 0 1 … -yn/6 xn/6 0 0 0 1]T (40)
其中,x1~xn/6是节点1到节点n/6的初始x轴坐标,y1~yn/6是节点1到节点n/6的初始y轴坐标,1≤m≤(n/6);
同理,将结构整体分别绕x轴和y轴旋转一个微小角度θ,构造转动刚体位移模式h4和h5,最后本实施例的刚体位移模式为:
H = h 1 h 2 h 3 h 4 h 5 h 6 = 1 0 0 0 0 - y 1 0 1 0 0 0 x 1 0 0 1 y 1 - x 1 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 · · · · · · · · · · · · · · · · · · 1 0 0 0 0 - y 4 0 1 0 0 0 x 4 0 0 1 y 4 - x 4 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 - - - ( 41 )
其中,x1~x4是节点1到节点4的初始x轴坐标,y1~y4是节点1到节点4的初始y轴坐标;
首先将平动刚体位移模式h1单位化,且与单位向量e1(e1为第一个元素为1,其它元素都为0的单位向量)一起用来构造Householder矩阵P1,使得:
P1*h1=e1 (42)
然后,构造Householder矩阵P2,值得注意的是,Householder矩阵P2的构造方式与Householder矩阵P1略有不同,构成新的特征向量h2′,并将h2′的第一个元素赋值为零,即
h2′=P1*h2 (43)
h2′(1)=0 (44)
将新构造的平动刚体位移模式h2′单位化与单位向量e2一起构造Householder矩阵P2
P2*h2′=e2 (45)
依次类推,根据构造Householder矩阵P2的方法逐个对平动刚体位移模式h3以及转动刚体位移模式h4~h6进行操作,构造新的平动刚体位移模式h3′以及新的转动刚体位移模式h4′~h6′,进一步构造Householder矩阵P3~P6
利用以上所构造的6个Householder矩阵去除结构的刚体模态,即用6个Householder矩阵对结构原始刚度矩阵做正交相似变换,得到去除刚体模态后的结构的刚度矩阵Kp
Kp=PTKP (46)
其中,K为原始的结构刚度矩阵;P=P1P2P3P4P5P6
根据上式生成的去除刚体模态后的结构的刚度矩阵Kp以及从Abaqus中生成的右端项F(为与Abaqus计算结果对比,采用Abaqus生成的外载F,包含支反力、支反力矩与施加的外力),做结构的静力分析:
KU=F (47)
其中,F为施加在结构上的力,即外载右端项;U为结构在力F下产生的位移。
将F与Kp代入上式:
KPUP=FP (48)
其中,UP=PTU;FP=PTF。
采用修正的共轭梯度法对上式进行求解,得到位移:
U=P*UP (49)
由于Kp的前6行6列都为0,并强制使FP的前6个值为0,从而不求解式(15)中的前6个方程,这样一来,UP的前6个值可以为任意值,进而求解出来的位移U就会包含刚体位移,最后只需减去其刚体位移,就可与真实值进行比较,即:
Un=U-kihi(i=1、2…6) (50)
其中,hi为第i个刚体位移模式,ki为系数。
将上式(50)求得的减去刚体位移后的位移与Abaqus求解的位移结果对比见表1所示。
表1位移求解结果与Abaqus求解的结果对比
把Abaqus求解的结果作为位移真值,去刚体位移的位移误差见图3所示。由此可见对于壳单元,本方法的结果是相当精确的,完全可以满足工程计算需要。
上列详细说明是针对本发明可行实施例的具体说明,该实施例并非用以限制本发明的专利范围,凡未脱离本发明所为的等效实施或变更,均应包含于本案的专利范围中。

Claims (7)

1.基于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
2.根据权利要求1所述的基于Householder变换的无约束结构静力分析方法,其特征在于,所述步骤1包括:
步骤11、假设将结构整体沿x方向平移一个单位位移,则在物体内不会产生任何应变,而此时结构的刚体位移模式为只在第6j+1个自由度上有位移值1,其它都为0,具体可写为:
h1=[1 0 0 0 0 0 1 0 0 0 0 0 … 1 0 0 0 0 0]T (16)
步骤12、根据步骤11的方法将结构整体分别沿y方向和z方向平移一个单位位移,构造平动刚体位移模式h2和h3
步骤13、假设将结构上的节点m绕z轴旋转一个角度θ,此时节点m的位移模式可写为:
h6=[-ymθ xmθ 0 0 0 1θ]T (17)
其中,xm是节点m的初始x轴坐标,ym是节点m的初始y轴坐标;
在上式(17)中约去θ,并将其扩展到整个结构,即结构整体绕z轴旋转一个角度θ,而在物体内不产生任何应变,此时结构的刚体位移模式可写为:
h6=[-y1 x1 0 0 0 1 -y2 x2 0 0 0 1 … -yn/6 xn/6 0 0 0 1]T (18)
其中,x1~xn/6是节点1到节点n/6的初始x轴坐标,y1~yn/6是节点1到节点n/6的初始y轴坐标,1≤m≤(n/6);
步骤14、根据步骤13的方法将结构整体绕x轴和y轴旋转一个角度θ,构造转动刚体位移模式h4和h5
3.根据权利要求1所述的基于Householder变换的无约束结构静力分析方法,其特征在于,步骤3中,去除刚体模态后的结构的刚度矩阵Kp中前六行的元素均为零,因此相应的要把变换后的外载右端项的前六个元素设置为零,即:
FP(1)=FP(2)=...=FP(6)=0 (19)。
4.根据权利要求3所述的基于Householder变换的无约束结构静力分析方法,其特征在于,对方程(15)求解的方法为修正的共轭梯度法。
5.根据权利要求4所述的基于Householder变换的无约束结构静力分析方法,其特征在于,已知Householder变换的性质为:
Pi=2vivi T-In (20)其中,vi为归一化后的n维Householder向量,
对方程(15)求解前,根据公式(12)和公式(20),仅需存储具有稀疏特性的刚度阵K和Householder变换矩阵Pi对应的向量vi以替代去除刚体模态后的结构的刚度矩阵Kp
6.根据权利要求4所述的基于Householder变换的无约束结构静力分析方法,其特征在于,采用修正的共轭梯度法求解方程(15)的方法为分步乘积的方式:
Kp*pk=PT*K*P*pk (21)
其中,pk为共轭梯度法中的修正方向;
公式(21)分解为三步进行:
sk=P*pk (22)
sk=K*sk (23)
sk=Kp*pk=PT*sk (24)。
7.根据权利要求6所述的基于Householder变换的无约束结构静力分析方法,其特征在于,计算公式(24)中PT*sk分解成6次乘法进行,每一次均将PT转换为Householder变换矩阵Pi对应的向量vi进行。
CN201410091321.7A 2014-03-12 2014-03-12 基于Householder变换的无约束结构静力分析方法 Expired - Fee Related CN103902764B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410091321.7A CN103902764B (zh) 2014-03-12 2014-03-12 基于Householder变换的无约束结构静力分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410091321.7A CN103902764B (zh) 2014-03-12 2014-03-12 基于Householder变换的无约束结构静力分析方法

Publications (2)

Publication Number Publication Date
CN103902764A CN103902764A (zh) 2014-07-02
CN103902764B true CN103902764B (zh) 2017-01-04

Family

ID=50994085

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410091321.7A Expired - Fee Related CN103902764B (zh) 2014-03-12 2014-03-12 基于Householder变换的无约束结构静力分析方法

Country Status (1)

Country Link
CN (1) CN103902764B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106919537A (zh) * 2017-03-07 2017-07-04 电子科技大学 一种基于FPGA的Jacobi变换的高效实现方法
CN108595893B (zh) * 2018-05-16 2021-06-01 电子科技大学 一种基于三层预处理子的三维力学模态仿真模拟方法
CN110704976B (zh) * 2019-09-30 2020-09-01 西北工业大学 航空航天薄壁壳体建模的五面体六节点体壳单元构造方法

Also Published As

Publication number Publication date
CN103902764A (zh) 2014-07-02

Similar Documents

Publication Publication Date Title
Kennedy et al. A parallel finite-element framework for large-scale gradient-based design optimization of high-performance structures
Liu et al. Smoothed point interpolation methods: G space theory and weakened weak forms
Damanpack et al. High-order free vibration analysis of sandwich beams with a flexible core using dynamic stiffness method
CN109902404B (zh) 不同阻尼形式的结构时程响应积分的统一递推计算方法
Vizotto Computational generation of free-form shells in architectural design and civil engineering
CN108986220B (zh) 一种加速有限元求解物体网格模型弹性变形的方法
CN101315649B (zh) 含大量输入端口的微机电系统降阶建模方法
Vieira et al. A higher order thin-walled beam model including warping and shear modes
CN106354954B (zh) 一种基于叠层基函数的三维力学模态仿真模拟方法
Sun et al. Topology optimization based on level set for a flexible multibody system modeled via ANCF
CN103902764B (zh) 基于Householder变换的无约束结构静力分析方法
Ruiz et al. Eigenvector sensitivity when tracking modes with repeated eigenvalues
CN106294975A (zh) 一种基于降阶模型的梁式结构自由振动分析方法
Woodgate et al. Fast prediction of transonic aeroelastic stability and limit cycles
CN109657408A (zh) 一种再生核粒子算法实现结构线性静力学仿真方法
Rezaiee-Pajand et al. An efficient mixed interpolated curved beam element for geometrically nonlinear analysis
Pargana et al. Fully integrated design and analysis of Tensioned Fabric Structures: Finite elements and case studies
Ren et al. An adaptive triangular element of absolute nodal coordinate formulation for thin plates and membranes
Ya Finite element models in stresses for plane elasticity problems
CN106126823A (zh) 一种基于提高迭代法稳定性和收敛性的位移求解方法
CN104008234B (zh) 密集模态含阻尼结构模型修正方法
Zhao et al. A new modeling method for flexible multibody systems
de la Fuente An efficient procedure to obtain exact solutions in random vibration analysis of linear structures
Shi et al. A new rotation-free shell formulation using exact corotational frame for dynamic analysis and applications
CN115146402A (zh) 高温合金的氧化过程模拟方法及装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C53 Correction of patent of invention or patent application
CB03 Change of inventor or designer information

Inventor after: Ding Hua

Inventor after: He Zhinian

Inventor after: Zhang Junbo

Inventor after: Li Tingting

Inventor after: Cai Zhihua

Inventor after: Xia Changqi

Inventor before: Ding Hua

Inventor before: Zhang Junbo

Inventor before: Li Tingting

Inventor before: Cai Zhihua

Inventor before: Xia Changqi

COR Change of bibliographic data

Free format text: CORRECT: INVENTOR; FROM: DING HUA ZHANG JUNBO LI TINGTING CAI ZHIHUA XIA CHANGQI TO: DING HUA HE ZHINIAN ZHANG JUNBO LI TINGTING CAI ZHIHUA XIA CHANGQI

C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170104