
第
31
卷第
6
期
振动与冲击
JOURNAL
OF
VIBRATION
AND
SHOCK
Vo
l.
31
No.6
2012
基于多重多级动力子结构的
Lanczos
算法
张盛
1
,
2
方杰
l
气张洪武
I
气陈飘松
1
,
2
(1.大连理工大学工业装备结构分析国家重点实验室,大连
116024;
2.
大连理工大学运载工程与力学学部工程力学系,大连
116024
)
摘
要:提出利用多重多级子结构技术与
Lanczos
方法求解超大型复杂结构动力特性的子结构算法。该算法利用
子结构周游树技术,分别对每个子结构进行
Lanczos
迭代,通过累加各个子结构的正交化系数组成全局三对角矩阵,最后
求解得到整体结构的特征值。算法能够计算超大型结构特征值和特征向量,计算效率高;消耗计算机资源少,稳定性高。
由于考虑了各子结构内部自由度对整体求解的贡献,算法精度得到显著提高,并与不作凝聚的单一整体结构分析具有相
同的计算精度,计算结果不受复杂子结构划分方式的限制。数值算例验证了所提出算法的正确与有效性。
关键词:多重多级子结构;子结构树周游;
Lanczos
方法;特征值;特征向量
中图分类号
TB122
文献标识码
A
Lanczos
algorithm
based
on
multi-level
dynamic
substructures
ZHANG Sheng
1
,2 , FANG
lie
1
,2 , ZHANG
Hong-w
μ
1
,
2
CHEN
Bi
αo-songl
,
2
(1.
State
Key
Laboratory
of
Structural
Analysis
for
lndustrial
Equipment
,
Dalian
University
of
Technology
,
Dalian
116024
, China;
2.
Department
of
Engineering
Mechanics
,
Dalian
University
of
Technology
,
Dalian
116024
, China)
Abstract:
A combined method based on multi -level substructuring
technique
and
Lanczos algorithm was proposed for
solving dynamic characteristics of large-scale
and
complicate structures.
1n
implementation of the method , a total
tridiagonal matrix was given by accumulating the orthogonal coefficients of
each
substructure ,
and
the eigenvalues
corresponding to the whole structure were solved. The method was able to solve eigenvalue
and
eigenvector problems for
super large-scale structures
and
had
good computational efficiency and stability , and less computational cost. Since the
contribution of all
intemal
DOF of
each
substructure to the dynamic characteristics of the whole structure was considered ,
the accuracy was improved significantly
and
the same as the result achieved by directly solving the whole structure without
condensation process.
1t
was found that the partitioning scheme of substructure does not affect numerical results ,
and
the
finite
element
modeling has great flexibility. The numerical results showed that the solutions obtained by the proposed
multi-level substructuring method are effective and valid.
Key
words:
multi-level substructure; substructure tree travelling; Lanczos
method;
eigenvalue;
eigenvector
随着工程技术的飞速发展,结构系统越来越庞大,
越来越复杂,如飞机、大型轮船、高层建筑、大型机械和
各种航天器。在分析计算大型复杂结构的动力特性和
动力响应时,有限元离散化所得到的系统自由度是成
千上万阶,有时甚至高达几十万阶、几百万阶。对于这
种庞大的多自由度系统,用传统的求解特征值方法求
解是十分困难的,甚至是不可能的。子结构方法是计
算大型复杂结构动态特性十分有效的方法。
对于广义特征值问题
[KJI
ct>
f=λ
[M
川剧的求
基金项目:国家自然科学基金
(90715037
,
10872041
,
51021140004);
因
家基础性发展规划项目
(20
1O
CB832704)
;国家高技术研究发展计划
(2009AA044501
)
收稿
~IWJ:
2010
-12
-07
修改稿收到日期:
20
11
-
03
-
08
第一作者张盛男,博士,讲师,
1976
年生
解,经典的子结构方法如模态综合法
[IJ
计算精度的好
坏直接由选取的位移表达式即假设模态决定。界面位
移凝聚法
[2J
是通过对子结构刚度阵和质量阵凝聚处
理,达到降阶目的。胡海昌
[3
J
利用解析的模态分析方
法构造了约束界面模态综合法。邱吉宝等
[4
J
采用半解
析法提出了三类精确动态子结构方法。而对于一般广
义特征值问题的求解,里兹向量法、子空间迭代法川、
Lanczos
算法
[6J
都是很实用的近似解法。喻永声等
l
7;
提出利用动态子结构周游技术实现子空间迭代,求解
大型结构的广义特征值方程,并取得较好的计算精度,
但其需要迭代收敛判断;文献
[8J
提出利用子结构凝聚
实现
Lanczos
算法的反迭代过程,实际上其求解仍然是
在整体结构中进行。
Lanczos
算法被认为是求解大型矩阵特征值问题