09-2有限单元法的程序设计.docVIP

第九章 有限单元法实用软件 4. 更为强大的网格处理能力 有限元法求解问题的基本过程主要包括：分析对象的离散化、有限元求解、计算结果的后处理三部分。 由于结构离散后的网格质量直接影响到求解时间及求解结果的正确性与否，近年来各软件开发商都加大了其在网格处理方面的投入，使网格生成的质量和效率都有了很大的提高，但在有些方面却一直没有得到改进，如对三维实体模型进行自动六面体网格划分和根据求解结果对模型进行自适应网格划分，除了个别商业软件做得较好外，大多数分析软件仍然没有此功能。 自动六面体网格划分是指对三维实体模型程序能自动的划分出六面体网格单元，现在大多数软件都能采用映射、拖拉、扫略等功能生成六面体单元，但这些功能都只能对简单规则模型适用，对于复杂的三维模型则只能采用自动四面体网格划分技术生成四面体单元。 对于四面体单元，如果不使用中间节点，在很多问题中将会产生不正确的结果，如果使用中间节点将会引起求解时间、收敛速度等方面的一系列问题，因此人们迫切的希望自动六面体网格功能的出现。 自适应性网格划分是指在现有网格基础上，根据有限元计算结果估计计算误差、重新划分网格和再计算的一个循环过程。对于许多工程实际问题，在整个求解过程中，模型的某些区域将会产生很大的应变，引起单元畸变，从而导致求解不能进行下去或求解结果不正确，因此必须进行网格自动重划分。自适应网格往往是许多工程问题如裂纹扩展、薄板成形等大应变分析的必要条件。 §9-2 有限单元法的程序设计 程序设计原则 比较研究著名的通用程序，可以总结出有限元程序开发一些比较成功的设计原则和基本经验。 1. 程序商品化 各程序在第一次发布后（相当我国鉴定会后）随着用户的要求与科学技术发展，要不断增添新的功能，改进原有内容。例如最著名的NASTRAN程序，第一次发表时仅有15万语句，完全是个线性程序，不断开发的结果到1991年67版本时，综合起来已有近100万条语句，包括各种动力、非线性、场分析功能，同时更新了单元库，编入了反映最新科技水平的理论与算法，可以说是结构分析技术的滚动积累。再如ANSYS程序，经过不断发展，现在已经成为独一无二的CAD/CAE协同环境，独一无二的全新耦合场解决方案，并且具有良好的二次开发功能，极大地方便了用户。 所以鉴定（初版本）后的维护、推广、销售、服务、扩充与改进应是程序整个开发一个极重要的阶段。一个软件由出场到成熟约需十年时间，要付出艰辛劳动，投人大量资金。 2. 用户界面友善 当前各有限元开发公司都着力于改进和开发前、后处理程序功能，有的几乎付出与开发程序新功能相同的投资，力求让用户使用方便，减少纸面工作量，降低数据输人错误，输人、输出直观，并尽可能把计算机语言变成人类的自然语言。为达此目的采用了CAD新的技术，结构造型功能很强，同时又应用计算机新的图形功能。另一友善设计是把程序设计得尽可能开放。用户在工作中很可能发现，必须使用一些新的单元、材料，增添新的力学功能，即使通用性很强的程序也难于满足这些要求，而必须增添之。如果以已有通用程序为开发平台在其上增添这些功能，则可免去大部分具有共性的工作，达到事半功倍的效果。因此，通用程序设计得有良好开放性，对用户有很大帮助。现各程序均有用户子程序来达到这一目的。但是，都带有很大的局限性，这是值得注意的。 3. 技术支持全面 通用程序几乎是力学大全，功能很多，必须由用户选择，因此讲课、辅导及技术咨询很必要，如有可能开用户讨论会还可反馈信息改进程序。国外各公司均很注意技术服务与支援工作，往往从事这些工作的人员人数相当于程序开发人员人数。完整的程序文档是对用户支援不可缺少的条件，一般应包括理论手册、用户手册、算例手册。理论手册讲明本程序各功能的理论依据、算法说明以及单元库、材料库，使用户明确程序功能、解题范围及可选用的方法、单元等。用户手册告诉用户如何使用程序、输人数据以及输出形式、图形操作等。算例手册应列出各种功能计算的例题，它们的计算模型，输人数据的格式，便于用户直观地理解或曰照猫画虎。对于一些高级用户程序员手册也是必要的，他们在程序中增添程序模块时可参考这些手册。 4. 程序有通用性 有限元的特点是可将千姿百态的各种结构离散成一些简单的单元，正如人类用砖、瓦、预制件等形状简单的基本件建成形状万干的建筑物，因而有限元能编成通用程序适用非常广泛的结构，解决广阔的力学问题。但在设计思想上不明确也会设计出局限性很大的程序。起码的功能应具有内存与内--外存交互求解功能，可以当解题规模小时单纯内存求解，以提高效率，题目规模大时，内--外交流求解。较丰富的单元库是这一要求的另一必备条件，单元是有限元法及其程序的细胞，只有丰富单元库才谈得上通用性。 5. 分析高效率 这是每个程序永远追求的目标，也是一个程序的重要指标，当今结构分析的规