基于ANSYS的波动有限元程序前处理研究.pdfVIP

基于ANSYS的波动有限元程序前处理研究 刘红帅 薄景山 刘德东 (中国地震局工程力学研究所黑龙江哈尔滨 150080) 摘要 波动有限元能够真实反映地震动传播特性，因而逐渐成为地震工程和岩土工程重要的数值模拟 方法之一，但现有波动有限元程序建模难度大，缺乏方便而快捷的前处理程序，使其在科研工作和实 际工程中的应用受到了很大的限制。为此，提出了一套快速易学的前处理思路和方法，即以ANSYS 有限元程序完成的复杂地质体网格划分为基础，采用Visual Basic语言编写了波动有限元程序的前处理 数据转换程序和波动有限元程序控制性参数的可视化输入界面，使波动有限元程序前处理变得轻松自 如，简单易用。经实例验证，该方法快速高效，对其他程序的前处理具有重要的借鉴意义。 关键词 波动有限元：复杂模型；前处理程序；自动生成 引言 我国学者廖振鹏教授将时域集中质量有限单元法和人工透射边界结合起来给出一种波动问题 的时域有限元解法【l’2】，其方法解波动过程最为直接，可以很好地模拟地震动的传播过程，因而 逐渐受到了科研人员和工程人员的高度重视和广泛应用。同时该方法在实现的过程中采用了显式 解法，降低了对计算机硬件的要求，因而受到大家的欢迎。杨柏坡主持编写的“二维显式有限元 波动程序包LAND”于1992年通过部级鉴定，为科研和工程应用奠定了坚实的基础f3]，周正华【4】、 李山有【5】等也对该程序的功能进行了扩展和改进。需要指出的是，原有的计算程序均采用 FORTRAN77语言，数组采用固定数组，为了能够充分利用计算机内存，不得不结合解题规模进 行编译，使用起来比较麻烦。为了能够使波动有限元程序自动结合解题规模充分利用内存，作者 采用具有先进的FORTRAN95对上述程序进行了改造升级，采用了动态数组，使其能够高效利用 Wave 计算机内存，避免了多次编译，同时将程序重新命名为WPA2D(2-Dimension Propagation Analysis)，并重新编写两侧边界不等高的边界处理程序、增加了弹粘塑性分析、粘弹性非线性等 功能，作者将另文介绍。尽管主分析程序的计算稳定性和功能已经达到了实用化程序，但现有程 序在模型建立及数据准备等前处理问题上却存在以下不足，使其在科研工作和实际工程中的应用 受到了很大的限制： (1)模型的建立只能靠数据文件来实现，不是很直观，不能像ANSYS或ALGOR等有限元软 件，可以直接进行图形的处理； (2)对于比较复杂的工程模型，在建模时需要各控制点详细的数据，容易出错，检查起来也 不是很容易； ．， (3)建模工作量大，花费时间长，直接造成了模拟计算的周期长、难度大。 _ 为解决显式波动有限元程序建模的不足，杨柏坡等用Fortran语言编写了显式波动有限元的前 处理程序[3】，对于水平成层的场地实现了快速、便捷的建模。吴兆营编写了基于ANSYS专门针 ··-一157·__—— 对土石坝的前处理程序旧】，该程序只能针对两侧边界等高的情况，但是对于具有复杂的地质结构 等地质体，其显式波动有限元模型的建立及网格划分仍然非常不便。 然而。许多有限元程序对于复杂工程地质体数值模拟的建立有明显优势。由于研究对象的单 元节点信息，与大型有限元软件完全一致，即将研究对象划分为具有一定形状且有一定编制规则 的节点组成的单元面。因此，作者利用通用的建模与网格划分功能强大的ANSYS软件对复杂工 程地质体建立相应数值模型(包括网格划分)，再通过数据转换实现WPA2D前处理的自动生成。 1 WPA2D程序的前处理 型，剖分时，直接选取该单元类型。根据波动分析要求，计算出单元所需的尺寸界限，以保证波 动分析的有效性和稳定性，并完成复杂地质体的