GIS平台视角下水利水电工程组织有限元前处理分析.docVIP

　　GIS平台视角下水利水电工程组织有限元前处理分析 1 绪论 1.1问题的提出 在 1943 年，Couraut 在其论文中表述了应用定义在三角形区域上的分片函数的最小势能原理来求解 Saint-Venant 扭转问题。在 20 世纪 50 年代，Argis、Clough等对结构力学有限元分析的计算机实现做出了巨大推动，并形象地描述：有限元法=Rayleigh Ritz 法+分片函数。此后，越来越多的数学家涉足这个领域，使得有限元的发展纳入数学的轨道，从逐渐完善的数学描述，到建立牢固的数学基础，以及提出各种新型非标准有限元等。应该说除了微积分之外没有一种数学方法会像有限元方法一样在科学、技术、工业工程以及社会科学中有如此广泛的应用，并获得如此巨大的成功。不断涌现的商业有限元软件已涉及到科学技术的广泛领域，每年有上万用户、花费数亿美元来不断研发新的有限元软件[3]。 根据有限元法求解工程问题的过程可以将它分为三个阶段：1）前处理过程：主要对分析对象进行有限元网格剖分（或结构离散）以及材料、边界、荷载等条件的定义；2）有限元分析：主要对有限元模型进行单元分析，有限元方程的建立以及求解；3）后处理过程：根据有限元分析的结果，以图形等方式将结果提供给用户，并辅助判断计算结果等。利用有限元方法分析工程问题的过程中，前处理过程约占所用时间的 40%~45%，计算分析约占 5%左右的时间，用于分析结果及评判的后处理过程约占50%~55%左右的时间。 有限元法在水利工程中目前已得到了广泛的应用，主要体现在水工结构、岩土工程、水力学和热力学有关问题的研究等方面， 特别是在水工结构的静动力分析及抗震设计、波动理论、温度场及温度应力分析、流固耦合作用、结构的安全监控及反馈分析、高边坡稳定及加固分析、围岩稳定分析、渗流场等问题上都有着重要的应用。 1.1.2 前处理技术的发展 有限元前处理主要是对分析对象进行有限元网格剖分、材料属性定义以及边界条件的确定等。其中，有限元网格剖分尤为重要，它将直接影响有限元分析结果的精度。在有限元法与计算机结合之前，网格离散工作必须由手工完成，这样不仅效率低下而且容易出错，通常要分析一个结构需要花费大量的时间在网格剖分工作上。在有限元前处理过程中，大部分时间都将用于有限元模型的建立以及有限元网格剖分[4]。有限元模型的简易程度也将影响有限元网格剖分的进度。网格剖分是有限元法的关键步骤，自有限元法诞生以来就有许多科学工致力于这方面的研究，特别是随着计算机辅助工程（puter Aided Engineering，CAE）的发展，这项技术已经趋于成熟[5]。 目前，国内外对有限元前处理的研究主要集中在有限元网格剖分方法和新的单元形式研究上。有限元网格剖分是有限元前处理过程中的重要组成部分。K.Ho-Le 将网格生产算法分成五类：映射法、栅格法、拓扑分解法、节点连接法和几何分解法[6~21,4]。目前主要是这五种方法的混合使用及现代技术的综合应用。 （1）映射法 映射法即是结构化网格生成方法，有事非结构化网格生成方法，其基本步骤是：通过适当的映射函数将待剖分物理域映射到参数空间中形成规则参数域；对规则参数域进行网格剖分；将参数域网格反向映射回物理空间，从而得到物理域的有限元网格。映射法可分为保角映射法、基于偏微分方程法及代数插值法三类。映射法一般可直接处理单连通域问题，对于复杂多连通域问题，需要先将区域分解称若干单连通区域，再对每个单连通区域运用映射法进行网格剖分[1]。 （2）栅格法 栅格法的基本步骤是：首先用一组不相交的栅格覆盖在目标范围上，保留在目标范围内的栅格，删除目标范围外的栅格；然后对与物体边界相交的栅格进行调整或再分解，使得更加逼近目标对象；最后对内部栅格和边界栅格根据栅格分布进行网格剖分，从而得到整个目标区域的有限元网格。栅格法可以分为正则栅格法和有限四（八）叉树法。 2有限元法的基本理论与方法 由于偏微分方程描述客观物理世界的重要性，其数值解法长期以来吸引着数学家、物理学家以及工程师们的关注[3]。有限元法作为求解偏微分方程的一个重要方法，随着计算机科学的发展，迅速地在工程中得到了广泛的应用并不断发展。有限单元法（The Finite Element Method，简称 FEM）是以变分原理为基础，将要求解的数学微分方程型数学模型转换成相应的变分问题，即泛函求极值问题，然后利用分片插值原理将变分问题离散为普通多元函数的极值问题，最终归结为一组多元代数方程组，求解该方程组即可获得该问题的数值解。 2.1有限元法的构造 约束条件主要分为两类，即强加约束条件和周期性约束条件。强加约束条件在力学问题及连续场问题中经常出现，如第一类边界条件就是常见的强加约束条件，对于强加约束条件的处理常采用近似处理方法（置大数法）和