EEP法超收敛计算的程序实现.pdfVIP

第四章 EEP 法超收敛计算的程序实现 第四章 EEP法超收敛计算的程序实现 4.1 引 言 为了将二维有限元EEP 法超收敛计算实现为通用计算机程序，我们在最初 的考虑中设计了两个方案。一是在现有的工程软件上进行二次开发。一是在原 有有限元线法通用程序的基础上自行编制有限元程序。 在现有的工程软件上进行二次开发是最初的想法。由于我们的超收敛计算 是一种后处理的办法，它的主要工作集中在得到常规的结点位移值以后，所以 我们希望能够近可能利用现有工程软件已有的模块。概括起来有以下几点考虑： 1．现有的有限元大型计算软件大多提供了较完善的计算功能和强大的前后处 理功能，特别是后者，为用户的建模和结果的考察提供了十分便利的条件， 图形界面的输入和输出尤为见长。在现有软件上进行二次开发，就可以充分 利用其这方面的优势，如本课题，就可以在其有限元位移结果的基础上，直 接修正其数据文件，利用其强大的后处理功能，实现修正后的结果的图形输 出和文档输出，而前处理的建模过程和中间的有限元计算过程可以原封不动 的利用起来，即减少了工作量，又达到了非常好的效果。 2 ．如果可以将研究成果顺利的应用于广泛使用的商业软件中，其实用价值是相 当可观的。如前所述，此算法的改进效果是明显的，工作量相对较少的同时 就提高了计算的效率，而设计人员往往最为关心的就是我们进行修正的应力 项。数值精度的提高使我们用相对少的单元就获得可以信赖的结果提供了可 能，大大提高了工作的效率，即有可能实现高效率、高效果、高效益的统一。 总之，基于有限元应力超收敛计算的工程软件二次开发是一项很有意义的工 作。 基于这种想法，我们首先选择了ANSYS作为我们的目标软件。经过一段时 间的调研，我们基本熟悉了在ANSYS上进行二次开发的工具，也了解了ANSYS 应力计算的处理。在证实其可行性的同时也发现了一些问题和困难。 首先，国内外在ANSYS上进行二次开发的例子不少，但多是针对具体的工 程或添加特殊的单元，而针对应力的二次计算而进行的二次开发的例子较少， 可以借鉴的经验和资料都很有限。另一方面，开发工具的专有性很强，方式途 - 57 - 第四章 EEP 法超收敛计算的程序实现 径多的同时也带来了操作复杂、上手较慢的特点，需要的准备期较长。另外， 与ANSYS二次开发相关的资料少，现有的资料多是集中在使用方面，特别是理 论资料方面的短缺，将会给本课题的进行带来一些阻力。其次，我们在进行算 法的理论研究是抽象到数学模型，而在ANSYS都是针对具体的物理模型使用特 定的单元，这个转换的过程也会给我们的研究带来一定的困难。再次，ANSYS 虽然提供了大量的结果数据接口，但作为中间的一些计算结果和过程被封装了 起来，而这其中的有些数据又恰恰是我们所需要的。如何去访问这些数据还是 有待解决的一个问题。 回到实现超收敛计算的目的，我们不仅仅是单纯的将理论“翻译”成代码， 我们希望通过这个程序来充分验证理论方面的结论。其中很重要的一部分内容 是单元边界上和内部任意点的超收敛性质，但是ANSYS没有提供这方面的接口。 通俗的说，大多数商业有限元软件多以h 法加密为基础，寻求在较密的网格下达 到计算精度的要求，所以一般只提供结点上的信息。这方面接口的欠缺使我们 难以验证超收敛理论中所有对内部点做的工作。 所以，为把我们的科研成果充分的体现出来，尤其是要体现单元内部点超 收敛的位移和应力值，我们最终选择了第二套方案，即在现有线法程序的基础 上自行编制有限元程序。有限元线法程序FEMOL 系列从90 年代初编制以来， 经过许多人的不断努力已初具规模，为自行编制有限元程序提供了坚实的基础。 另一方面，从最初的网格到整个超收敛理论的建立，我们都和有限元线法有着 密不可分的联系，在现有的线法程序进行开发从思路上合理连贯，也可以利用 到许多现有的资源。 4.2 系统设计和程序概述 整个程序的系统框架沿用线法程序的做法，即以VB 程序作为集成的总控制。 前处理程序包括数据输入机制的建立和网格的图形显示，这部分程序直接利用 VB 进行开发。核心计算程序由 Fortran90 语言来编制，它也分为前处理、求解 和后处理部分，分别进行相应的数据操作和计算。结果文档的输出返回到VB 的 界面上