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

基于并行的特大增量步算法在计算固体力学中的应用

龙丹冰  
【摘要】:广义逆力法是一种基于加载全过程的迭代算法,因解决材料非线性问题时增量步长不受计算精度限制只由加卸载状态决定,又称特大增量步算法(Large Increment Method,简称LIM)。作为一种有限元算法,LIM具有很强的空间上的可并行性,无需进行子结构划分就可以直接以单元为最小单位并行求解本构方程。由于引入了整体迭代优化算法,在分析材料非线性问题时,LIM具有很强的时间上的可并行性,这是这一算法独具的优点。在复杂的加载情况下,时间上的并行可大大提高计算速度。本文分析了本算法的时间并行性优势,找出需要进一步完善的算法并行化的前期工作和难点,提出了包含数值并行、空间并行及时间并行的算法并行化解决方案。但由于并行计算是应对大规模计算需求的产物,计算规模越大,并行越有效。而在本文的工作开展之前,LIM由于单元库贫乏,只能应对杆系结构的分析及简单的二维问题的计算,无法应对复杂结构的分析,计算规模无法增大。因而当务之急是丰满单元库。当算法无论从计算精度还是计算规模上都满足实用要求,再进行算法的并行化,将能充分发挥算法的并行特点,大大提高计算效率。 本文的主要工作着重于LIM在计算固体力学中的应用扩展——即开发二维和三维实体单元、构建LIM单元库。论文首先将杆系单元写入LIM的串行程序单元库,以其为例讨论并突出了LIM独有的计算特点。为丰满单元库,文中创新性地提出了两种建立二维和三维实体单元的通用办法,建立了两类单元,其一命名为“节点力型”单元,另一种命名为“应力型”单元。 “节点力型”单元是选取单元部分节点力作为单元广义内力的一类单元。单元广义内力由能够唯一确定单元内部应力分布的最少独立未知力组成,可以有多种选择。在综合考虑了各种二维、三维单元之后,本文提出了选取单元基本静定体系的特殊方式,确定了基本静定体系的局部坐标并选取了合理的基本力系作为单元广义内力系。在这样明确定义的基本体系下,文中通过坐标转换及平衡分析建立了单元的控制方程通式。通式没有限定单元的形状与节点个数,给定节点数目即可通过二维或三维通式建立多节点单元的控制方程。而且通式建立过程系统规范,不存在随机性,非常利于编程计算。本文已将四个二维及两个三维“节点力型”单元增加至LIM单元库,并在LIM串行程序中完成了单元控制方程通式及上述各单元的编写工作。数值算例表明,使用“节点力型”单元时,LIM的计算收敛速度及计算精度满足使用要求,精度受一定畸变度的不规则单元影响小,数值计算不受全局坐标系取向及节点编号的影响。 “应力型”单元分别采用不同的插值函数逼近应力场和位移场,再通过虚功原理建立单元的控制方程。由于将应力场插值函数的线性无关系数定义为单元广义内力,因而命名为“应力型”。本文首先通过假设应力函数得到应力各分量表达式,并建立了二维及三维单元的单元控制方程;而后,基于LIM的控制方程,文中提出了判别零能变形模式的充要条件,并分别从单元柔度矩阵是否病态以及是否存在零能变形模式两个方面考查了各单元。为抑制零能变形模式,论文提出了另一种构造“应力型”单元的方法,即基于基本位移模式的等函数法,并利用等函数法构造了一个新的无零能变形模式的六面体单元。基于以上两种方法,本文为LIM单元库增加了六个二维及三个三维“应力型”单元,并完成了所有单元的程序编写工作。数值算例表明,使用“应力型”单元时,LIM的计算收敛速度及计算精度均满足使用要求。 为了便于用户使用LIM求解弹塑性问题,文章在最后提出了适用于LIM的一维至三维一致弹塑性矩阵的推导,并以数值算例展示了LIM在求解弹塑性问题中的适用性及并行潜能。 本文共为LIM算法增加二维平面实体单元十个、三维空间实体单元五个。建立单元的方法具有可扩展性,用户可以根据需要、依据文中提出的方法增加适用单元。论文研究成果使得一直停滞不前的LIM的研究向“在计算固体力学中的应用”方面迈进了一大步。同时,论文作者独立完成了LIM的面向对象的串行程序的全部编写工作。在此基础上,对程序的进一步优化工作以及实现并行化计算将可以顺利开展。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 龚尧南;计算固体力学的发展及其在航空航天工程中的应用[J];计算力学学报;1993年02期
2 杨海天;有限体积法在固体力学中的应用[J];大连理工大学学报;1995年03期
3 温卫东,高德平;三维弹性边界元分析中面力不连续和1/r积分奇异性问题的处理[J];南京航空航天大学学报;1990年04期
4 杨忠清,钱祖钰,张霄文;应力极值分布、强度威布尔分布的可靠性[J];南京航空航天大学学报;1991年02期
5 毛军,薛琳,熊艳;在本科力学课程中引入计算力学内容的可行性分析[J];北方交通大学学报(社科版);2003年02期
6 詹世革;孟庆国;白坤朝;;从“十五”期间基金资助情况看力学学科发展现状[J];力学进展;2005年04期
7 谭礼乔;有限点法及其在固体力学中的应用[J];武汉大学学报(工学版);1988年05期
8 苏永华;何满潮;;利用矩点估计法简化响应面可靠度指标的计算[J];工程力学;2007年07期
9 滕弘飞,滕晓雷;机械动力优化设计的摄动逆求法及其应用[J];大连理工大学学报;1987年01期
10 郑贵明,王辉林,朱永强;冷藏集装箱有限元分析[J];山东工程学院学报;2002年04期
11 刘西拉,王开健,张春俊;空间刚架结构的广义逆力法[J];力学季刊;2004年01期
12 吴逸群;;多个约束下的高阶拉氏乘子法及其应用[J];解放军理工大学学报(自然科学版);1986年03期
13 蒋利顺;刘定生;;遥感图像K-Means并行算法研究[J];遥感信息;2008年01期
14 刘忠;周长顺;杨涛;;混合配筋预应力混凝土密肋楼盖三维有限元数值模拟[J];湘潭大学自然科学学报;2008年04期
15 A.Needleman,董务民;9.计算力学[J];力学进展;1987年03期
16 王思群,魏紫銮;线性规划内点法的并行计算[J];数值计算与计算机应用;1996年01期
17 宋贞霞;丁海平;胡颖平;;土-结构动力相互作用几种计算模型的比较[J];自然灾害学报;2009年05期
18 范峰;马会环;沈世钊;;半刚性螺栓球节点受力性能理论与试验研究[J];工程力学;2009年12期
19 茅春浦;上海交通大学工程力学系简介[J];水动力学研究与进展A辑;1986年02期
20 曹丹青;使用I—deas软件进行有限元分析[J];纺织机械;1995年06期
中国重要会议论文全文数据库 前10条
1 邢誉峰;;写教材的一点心得体会[A];探索的脚步——“十一五”北京高等教育教材建设论文集[C];2010年
2 隋永枫;孙义冈;丁旭东;刘象拯;张军辉;;不调频叶片设计改型研究[A];2009年中国动力工程学会透平专业委员会2009年学术研讨会论文集[C];2009年
3 王家林;;实体单元的转动位移分析及应用[A];中国计算力学大会'2010(CCCM2010)暨第八届南方计算力学学术会议(SCCM8)论文集[C];2010年
4 林少玲;陈铭年;徐健全;李进彬;;ANSYS结构分析单元的分类[A];福建省科协第四届学术年会——提升福建制造业竞争力的战略思考专题学术年会论文集[C];2004年
5 周凌远;李彤梅;;实体单元几何非线性CR列式有限元方法[A];第18届全国结构工程学术会议论文集第Ⅰ册[C];2009年
6 高艳平;焦亦凡;王余庆;;斜坡场地上建筑物对坡面地震反应的影响[A];中国土木工程学会第八届土力学及岩土工程学术会议论文集[C];1999年
7 饶舒;王克海;叶英华;;大跨径单箱单室刚构桥箱梁的空间有限元分析[A];第15届全国结构工程学术会议论文集(第Ⅲ册)[C];2006年
8 林少玲;陈铭年;徐健全;李进彬;;ANSYS结构分析单元的分类[A];福建省科协第四届学术年会提升福建制造业竞争力的战略思考专题学术年会论文集[C];2004年
9 占玉林;赵人达;毛学明;;钢-混凝土组合桥面板荷载试验的数值模拟[A];第十四届全国工程设计计算机应用学术会议论文集[C];2008年
10 王蓬波;刘振成;张崇亮;;钢板弹簧仿真设计[A];第五届中国CAE工程分析技术年会论文集[C];2009年
中国博士学位论文全文数据库 前10条
1 龙丹冰;基于并行的特大增量步算法在计算固体力学中的应用[D];上海交通大学;2012年
2 陈思;计算力学中高精度无网格法基础理论研究[D];上海交通大学;2007年
3 韩丛英;若干优化问题的并行算法研究[D];上海交通大学;2008年
4 周凌远;斜拉桥非线性理论及极限承载力研究[D];西南交通大学;2007年
5 田文;多视图图像的快速三维场景重建[D];华中科技大学;2010年
6 刘毅;面向人群的并行多目标疏散模型研究[D];武汉理工大学;2012年
7 范洪博;快速精确字符串匹配算法研究[D];哈尔滨工程大学;2011年
8 谢洪阳;弹性地基板动力问题的数值分析[D];华中科技大学;2006年
9 白洪涛;基于GPU的高性能并行算法研究[D];吉林大学;2010年
10 张文学;预应力混凝土连续箱梁局部应力分析及拉—压杆设计[D];同济大学;2007年
中国硕士学位论文全文数据库 前10条
1 陈爱平;基于Hadoop的聚类算法并行化分析及应用研究[D];电子科技大学;2012年
2 赵悦;子空间类超分辨算法并行化研究[D];哈尔滨工业大学;2010年
3 王璟瑶;基于Mobile CPU/GPU协同计算回声消除算法的优化与实现[D];电子科技大学;2011年
4 肖韬;基于MapReduce的信息检索相关算法并行化研究与实现[D];南京大学;2012年
5 邹彩辉;基于Hadoop平台的自适应局部超平面K近邻算法的研究[D];华南理工大学;2011年
6 党纪;现浇混凝土蜂窝式空心双向板的试验研究[D];东南大学;2006年
7 管晔;大型卵形消化池结构静力数值分析[D];天津大学;2006年
8 单玉乐;类电磁机制算法研究和改进[D];太原科技大学;2010年
9 张宇宸;人工颈椎间盘置换对多节段下颈椎活动影响的三维有限元分析[D];河北医科大学;2008年
10 文思;基于Hadoop的K近邻分类算法的研究[D];华南理工大学;2011年
中国重要报纸全文数据库 前10条
1 安世亚太 孟志华;ANSYS在航空薄壁构件结构设计的应用[N];中国航空报;2005年
2 安世亚态 马骏;航空CAE中的梁壳单元[N];中国航空报;2005年
3 ;亚光亚装饰:模块家居变革家装业[N];北京商报;2007年
4 张自然;企业家应该做到“四不”[N];中国电子报;2007年
5 张惟;优化指挥体制结构谋求指挥效能增值[N];解放军报;2007年
6 宋学春;青岛创立共建文明网络新机制[N];人民日报;2008年
7 ;ANSYS 9.0新功能[N];中国航空报;2005年
8 裘建;建数字国土迎信息时代[N];中国矿业报;2001年
9 ;ANSYS 9.0融合更多功能[N];中国计算机报;2005年
10 中航工业计算所 牛文生;计算机技术的一些新进展[N];中国航空报;2011年
中国知网广告投放
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978