收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

结构布局修改静力重分析的预条件共轭梯度法

李正光  
【摘要】:在结构设计和优化过程中,有时需要对结构进行几十次甚至上百次的修改,而每作一次修改都执行一次完全的有限元分析其计算成本是巨大的。为了减少计算成本,以不直接求解结构修改后的隐式分析方程,而利用初始分析的相关信息设计高效算法来求解修改后结构响应的重分析方法,日益受到人们的重视并得到飞速发展。 传统的结构修改主要侧重于参数修改。随着结构设计与优化技术的不断发展,人们发现结构布局的设计与优化比布局固定的设计与优化更具有节约潜力。 然而,由于在布局修改中,结构的有限元模型发生变化,其自由度数量的改变将导致分析方程组的结构和阶数发生变化。目前大部分重分析方法都是针对结构自由度没有变化的情形。 发展适用于结构布局修改的重分析方法是具有巨大挑战性和经济价值的。 预条件技术是科学与工程计算中复杂问题高效求解器的最关键组成部分。 所谓预条件就是将一个难处理的问题转化成另一个解能够被迅速逼近的问题。 本文主要应用预条件技术解决结构布局修改的静力重分析问题。 论文主要由五部分组成。第一章绪论,主要介绍了重分析问题的工程背景和国内外的研究状况。 第二章首先回顾了解结构修改后自由度未发生变化的重分析问题的组合逼近方法和预条件共轭梯度法,然后从计算复杂度和收敛精度等方面将二者做了比较,指出了预条件共轭梯度法的优越性。第三、四、五章是作者的主要研究工作(即摘要的1-3部分),分别建立了求解结构自由度减少、增加和一般布局修改的静力重分析问题的预条件共轭梯度法。 结点自由度减少情形的静力重分析方法 给定一个初始结点自由度为 的结构,相应的刚度矩阵记为 ,荷载向量记为 ,对应的位移 可由求解如下平衡方程 WP=93 (1) 来获得。由初始分析,刚度矩阵 可写成如下的Cholesky分解形式 (2) 式中 是上三角矩阵。 假设增加和删除一些单元并且从初始结构中删除一些结点后,结构有 个自由度。记 与 分别是结构修改后 阶刚度矩阵和荷载向量。则修改后的分析方程组可以写成如下形式: (3) 式中 是位移向量。我们选取 (4) 为预条件矩阵,再利用预条件共轭梯度法来求解。在方程(4)中, 是 阶投影矩阵,变换 通过除去m-维向量 中对应于被删除结点的分量而得到n-维向量 ; 是将n-维向量 转化为m-维向量 , , 在保留自由度上的分量与 的对应分量相同,而在被删除自由度上的分量为零。 本文建立了求解该问题的预条件共轭梯度法的有效辅助,研究了其计算复杂性。 结点自由度增加情形的静力重分析方法 在结点自由度为 的初始设计基础上,假设增加和删除一些单元并增加一些新的结点,这使得结构的自由度数目增加。修改后的分析方程可以写为如下形式: (5) WP=94 式中 表示在修改设计中新增加结点的自由度数目, 对应新增结点的子刚度矩阵, 为联系初始结点和新增结点的子刚度矩阵。方程(5)可以分解为以下两个线性系统: (6) (7) 由方程(7)解出 ,然后代入方程(6)可以得到如下方程: . (8) 式中 , , . (9) 这样就把一个 维的问题化为一个 维的问题。本文提出的第一个方法是以 为预条件矩阵应用预条件共轭梯度方法求解(8)。在求得 后代入(7)式即可求得 。 本文提出的第二个方法是取 (10) 为预条件矩阵,再利用预条件共轭梯度方法来求解方程(5)。 本文证明了 作为预条件矩阵的必要性质——对称正定性,建立了求解该问题的两种预条件共轭梯度法的有效辅助并研究了它们的计算复杂度。 结构一般布局修改的静力重分析方法 在结点自由度为 的初始设计基础上,假设增加和删除一些单元,删除了 个已有自由度并增加了 个新的自由度,修改后结构的分析方程为 (11) WP=95 式中 为 阶刚度矩阵。修改后的刚度矩阵可以分块的形式写为: (12) 下角标中小写字母表示子刚度矩阵的阶数, 其中 。 本文构造 (13) 为预条件矩阵,用预条件共轭梯度方法求解方程(11)。在方程(13)中, 是 阶投影矩阵,变换 通过除去m0-维向量 中对应于被删结点的分量而得到m-维向量 ; 将m-维向量 变换为m0-维向量 , 中对应保留自由度上分量与 向量的对应分量相同,而在被删除自由度上的分量为零。 本文证明了 作为预条件矩阵的必要性质——对称正定性,建立了求解该问题的预条件共轭梯度法的有效辅助并研究了该算法的计算复杂度。 在本文提出的预条件方法中,预条件算子的构造主要是基于修改前的刚度矩阵,这样,预条件算子的构造是廉价的且预条件后的系统很容易求解。预条件方法的主要优点是计算量小,需要的存储空间小,收敛速度快,算法容易辅助并能自适应地监视逼近解的精度,因而该方法能够大大地节省计算费用。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 李铮,邵新慧,李长军;矩阵B~TA~(-1)B的特征值估计及预条件处理[J];东北大学学报(自然科学版);2005年06期
2 刘喜武,刘洪;实现稀疏反褶积的预条件双共轭梯度法[J];物探化探计算技术;2003年03期
3 刘兴平;广义预条件迭代方法[J];计算数学;1990年01期
4 卫加宁,王伟沧,皮新明,章社生;一类虚拟边界预条件多重网格并行算法[J];武汉理工大学学报(交通科学与工程版);2001年01期
5 黄文彬,许传炬;Poisson方程谱元法的一个有限元预条件分析[J];厦门大学学报(自然科学版);2003年04期
6 孙丽英;线性方程组的AOR预条件迭代法的两个性质[J];广东教育学院学报;2003年02期
7 庄伟芬,卢琳璋;(I+S_(max))预条件Gauss-Seidel迭代法进一步探索[J];厦门大学学报(自然科学版);2004年S1期
8 雷刚;畅大为;王慧勤;;预条件[I+C(α)]加速2PPJ型方法的收敛性[J];江汉大学学报(自然科学版);2006年01期
9 卫加宁,章社生,Yakup Paker;区域分裂并行计算中预条件迭代若干性质[J];武汉理工大学学报(交通科学与工程版);2001年02期
10 雷刚;;两类预条件后SOR迭代法收敛性的比较[J];江南大学学报(自然科学版);2011年02期
11 李春光,游兆永;一类本身具有预条件功能的ABS型共轭梯度算法[J];数值计算与计算机应用;1999年02期
12 宋梅,张石峰;共轭梯度法在新井参数设计中的应用[J];新疆地质;2004年03期
13 黄维章;一些模型问题高阶离散格式的低阶预条件[J];科学通报;1992年07期
14 李荣;畅大为;;(I+S_(max))预条件的2PPJ迭代方法及比较定理[J];太原城市职业技术学院学报;2009年02期
15 黄文彬,许传炬;二维Poisson方程谱元法有限元预条件分析[J];福州大学学报(自然科学版);2004年04期
16 雷刚,王慧勤;一类预条件下2PPJ型方法收敛性的加速[J];宝鸡文理学院学报(自然科学版);2005年03期
17 李宏峰;王慧勤;杨亚强;;预条件下2PPJ型方法收敛性的加速[J];山西煤炭管理干部学院学报;2006年01期
18 雷刚;王慧勤;畅大为;;预条件下2PPJ型方法收敛性的加速[J];江西师范大学学报(自然科学版);2006年01期
19 雷刚;王慧勤;;一类新预条件下AOR迭代法收敛性的讨论[J];安徽大学学报(自然科学版);2007年03期
20 王慧勤;雷刚;;预条件P_c=(I+C)后SSOR迭代法收敛性的加速[J];宝鸡文理学院学报(自然科学版);2007年03期
中国重要会议论文全文数据库 前10条
1 李正光;吴柏生;;结构单元和结点删除后的一种重分析方法[A];第十二届全国结构工程学术会议论文集第Ⅰ册[C];2003年
2 梅金顺;刘洪;;预条件方程组及其应用[A];中国科学院地质与地球物理研究所二○○四学术论文汇编·第三卷(油气·矿产·水资源)[C];2004年
3 谷同祥;刘兴平;;预条件多搜索方向共轭梯度方法[A];中国工程物理研究院科技年报(2003)[C];2003年
4 李梅霞;王长钰;;线搜索下带误差项的Dai-Yuan共轭梯度算法(英文)[A];中国运筹学会第七届学术交流会论文集(中卷)[C];2004年
5 刘喜武;刘洪;;实现稀疏反褶积的预条件双共轭梯度法[A];中国科学院地质与地球物理研究所二○○三学术论文汇编·第四卷(油气资源)[C];2003年
6 龙桂华;李小凡;张美根;;粘弹性声波波形反演中的预条件最速下降法研究[A];中国地球物理学会第二十四届年会论文集[C];2008年
7 石磊;郭良辉;孟小红;;磁场源位置估计的扩展倾角梯度法[A];中国地球物理2010——中国地球物理学会第二十六届年会、中国地震学会第十三次学术大会论文集[C];2010年
8 李振华;王彦飞;杨长春;;一种能提高叠前偏移剖面分辨率的预条件正则化反演方法[A];中国地球物理·2009[C];2009年
9 连淑君;王长钰;;在Armijo型线搜索下共轭梯度法簇的全局收敛性[A];中国运筹学会第七届学术交流会论文集(上卷)[C];2004年
10 刘景昭;张玉忠;;变分不等式问题的带非负约束转化的一类信赖域法[A];中国运筹学会第七届学术交流会论文集(下卷)[C];2004年
中国博士学位论文全文数据库 前10条
1 李正光;结构布局修改静力重分析的预条件共轭梯度法[D];吉林大学;2004年
2 单润红;Toeplitz系统的高效预处理技术和分布式并行算法研究[D];国防科学技术大学;2004年
3 陈(日午);关于真实感服装生成算法的研究[D];西北工业大学;2000年
4 钟萍;Newton-PCG型算法的效率分析[D];中国农业大学;2002年
5 吴建平;稀疏线性代数方程组迭代法中的预处理技术研究[D];中国人民解放军国防科学技术大学;2002年
6 李永杰;调强放射治疗中的优化技术研究[D];电子科技大学;2004年
7 张振宇;数值线性代数中的若干问题[D];复旦大学;2003年
8 王振营;人口迁移的规律[D];中国人民大学;1993年
9 连淑君;共轭梯度算法的全局收敛性研究[D];大连理工大学;2004年
10 秦梅;一类奇异方程组的求解和扰动分析[D];复旦大学;2006年
中国硕士学位论文全文数据库 前10条
1 孟国艳;SAOR-AL预条件共轭梯度法[D];太原理工大学;2003年
2 肖映雄;一类晶格材料(Lattice Materials)的无结构代数多重网格法[D];湘潭大学;2001年
3 白文军;线性方程组的预条件迭代法和误差估计[D];陕西师范大学;2010年
4 李娟;一类共轭梯度算法的收敛性[D];首都师范大学;2005年
5 王炜;线性方程组新的预条件迭代法的研究[D];太原理工大学;2012年
6 黄文姬;线性约束优化的仿射内点共轭梯度路径方法及其应用[D];上海师范大学;2010年
7 刘翠;离散正弦变换与TOEPLITZ系统的快速算法与并行算法[D];国防科学技术大学;2004年
8 卞文良;一类新型共轭梯度算法[D];首都师范大学;2004年
9 蔡良弟;求解微波电磁场反问题的神经网络方法[D];大连理工大学;2005年
10 徐畅;一类共轭梯度法的研究[D];南京航空航天大学;2004年
中国重要报纸全文数据库 前10条
1 齐继成 毛冬蕾校;共轭雌激素依然领先[N];医药经济报;2001年
2 本报记者 李佳 王健玮;罗盘准星:科学·生态[N];中华建筑报;2005年
3 王如松;“摊大饼”有违生态城市建设理念[N];光明日报;2007年
4 张未名 作者单位:中国科学院上海微系统与信息技术研究院;环保要与发展阶段相适应[N];中国经济导报;2005年
5 江其勤;对南方构造提出新认识[N];中国石化报;2006年
6 天津社会科学院城市经济研究所 陈柳钦;共轭互补的风险投资与高新技术产业[N];中国高新技术产业导报;2007年
7 邵斌;国家重大科技专项关注百姓杯中奶[N];大众科技报;2005年
8 本报记者 肖琳;让科技释放经济发展潜能[N];中国高新技术产业导报;2003年
9 宁宕;巴西:研制成新型奶牛饲料[N];广东科技报;2003年
10 王新佳;奶业重大专项实现关键技术突破[N];中国高新技术产业导报;2005年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978