用于单元刚度计算的矩阵类程序的设计与开发.pdfVIP

于卜面一?用于单元刚度计算的矩阵类程序的设计与开发???黔整恻度艇，计篁翱蜷?????????汪继勇，兰箭’，华林，王德清配函数——???函数和内联函数，可有效地提高单元矩阵计算效率，从而提高有限元软件执行效率。本???????????，???，???辝??武汉理工大学材料加工系；湖北武汉???摘要：在有限元软件开发中，程序模块的重用和执行效率是影响软件性能的两个关键问题。运用面向对象语言开发有限元软件中的单元刚度矩阵计算程序，可有效解决其重用问题；运用引用变量、内存动态分文开发的单元刚度矩阵计算程序是用??镅员嘈础Ｓ氪?车拿嫦蚬?逃镅员嘀频挠邢拊3绦蛳啾龋?梅ū?制的程序执行速度快、重用性好，易于实现多种类型单元共存的情况，是解决大型有限元程序编制的有效方法。最后，文中通过应用实例比较，说明了此矩阵类的高效性。关键词：有限元；面向对象；矩阵类；???曲有限元方法是从应用电予计算机进行结构力学中的矩阵法计算而发展起来的【¨。经过?多年的发展，有限元法在工程问题的很多领域中得到广泛应用，成为目前最流行的、卓有成效的解决非线性问题的数值计算方法。有限元分析过程一般包括?鼋锥危河邢拊DＰ徒?⒑褪?菔淙虢锥?又称前处理??邢拊7?析计算阶段?蠼馄?和有限元分析结果的输出与评估阶段?殖坪蟠????Ｆ渲校?蠼馄饔址治H?霾糠郑????，???????，????韈??????．????????籓???猳????篗???见图????．???．，?????????．?????，?????????甎??????．???????，??????????????????甌?????．????????．???篺??????籆?图??有限元分析关系图图??矩阵类的结构????：?????，??????琧?????????．????．??????髓处理后处理单元蛔毒计冀求解线性方程姬?时??廿????，???‘ ピ8斩染卣罄喑绦蛏杓朴肟7?法跨越的缺陷——重用性和扩充性差，研究与开发基于面向对象的求解器程序是当前急需解决的。?岣呔卣罄喑绦蛑葱行?实姆椒?到目前为止，工程师和程序设计人员已经研究和开发了很多通用与专用的有限元软件，如????????等，其求解器都是传统的结构化?嫦蚬??有限元程序。该程序在应用中逐渐显露出其无由图?籰可知，单元计算是实现求解器求解的一个重要环节。在单元计算中，单元和节点是最基本的要素且数量较多，要计算出各个单元刚度矩阵?且桓龃罅康木卣蠹扑愎?獭？杉??岣叩ピ卣蠹扑阈?率是提高求解器运算效率的关键。??面向过程的程序设计方法传统的面向过程编程方法，如?????语言等，它们对问题域的认识和描述不是以问题域中的固有的事物作为基本单位，并保持它们的原貌，而是打破了各项事物之间的界限，在全局范围内以功能、数据为中心进行分析??Ｒ杂邢拊7治鑫@??嫦蚬?坛绦虬延邢拊7治龉?炭醋饕恍┕δ芎妥庸δ艿淖楹希?将数据流的加工和处理作为考虑重点。这样的分析结果不能直接反映实际情况，容易隐藏一些对问题域的理解偏差，与后续开发阶段的衔接就比较困难。在进行大型有限元分析时，由于数据量巨大，这种采用数据流为对象的编程思想，使得前面程序的一点小错误就可导致后面所有程序的崩溃，不利于程序的扩充。而且用此法编制的程序很不利于代码的重用，不利于在此基础上进行再次开发。例如，在单元计算中，将三角形膜单元的程序变换到四边形单元时，以前所编制的程序基本上都得重新编写，代码的可重用率太低，这将是程序员的一个沉重负担。??面向对象的程序设计方法面向对象程序设计方法????鏑?等，是二十世纪八十年代中期迅速发展起来的程序设计方法，它综合了结构化、数据封装等基本概念和方法，被誉为软件工程学的最新进展和现代软件的核心技术之一??Ｃ嫦蚨韵笊杓品椒ㄊ鞘辜扑慊?绦虻拇?砉?烫??嗣嵌钥凸凼澜缡挛锏拇?恚?渥钕灾?奶氐闶且?入了类的概念。它具有如下优越性：?程序可读性好，代码结构性更好、更合理：?软件的可重用性、可扩充性、可维护性及可靠性更好：?尤其适用于中等和大型的程序开发。面向对象程序设计方法优越性是由其本身的三个显著特点所决定的：封装、继承、多态与功能重载。??有限元单元刚度矩阵计算程序的开发目前，在有限元分析程序的开发过程中，通常的做法是采用??镅员嘀魄啊⒑笾么?聿糠值拇?耄???限元求解器模块的各个部分那则使用结构化语言????嘈础????语言虽然是世界上广泛流行的、最适用于数值计算的一种计算机语言，但是它最大的缺点是程序难读懂，可扩充性和可重用性较差。本文用??镅钥7⒌挠邢拊5ピ8斩染卣蠹扑愠绦颍?淅嘟峁谷缤?—?基类??????庾傲司卣蟮淖罨?疽K兀?K?芯卣笞永嗵峁┩骋唤涌凇Ｔ诨?嗟南乱徊悖?缮??了列向量?????⑿邢蛄縍????、满阵?????和按行存储的下三角对称矩阵