有限元发展现状学院建筑工程学院日期20161223硕士生.ppt

有限元发展现状 学 院：建筑工程学院 日 期：2016.12.23 硕士生：徐余锋 专 业：建筑与土木工程 二 一 1.概述 随着我国科技的迅速发展，人们更是不断的创建更快捷，更简便，规模更大的建筑物以及更加精密的设备。到今天有限元在工程上得到广泛应用，理论和算法都已经日趋完善。有限元法(Finite Element Method，简写为FEM)是求解微分方程的一种非常有效的数值计算方法，用这种方法进行波动数值模拟受到越来越多的重视。有限元法起源于固体力学，并逐步扩展到热传导、计算流体力学、电磁学等不同领域，已经成为数学物理中很重要的数值计算方法。  有限元在各个领域的广泛使用已使设计水平发生了质的飞跃，主要表现在以下几个方面： （1） 增加产品和工程的可靠性； （2）在产品的设计阶段发现潜在的问题； （3）经过分析计算，采用优化设计方案，降低原材料成本； （4）缩短产品投向市场的时间； （5）模拟试验方案，减少试验次数，从而减少试验经费。  2.主要表现 有限元法的基本思路是将计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。 3.基本思路 随着有限元技术应用的不断扩大，其发展呈现如下特点: （1）单一场计算向多物理耦合场问题的求解发展； （2）由求解线性问题发展到求解非线性问题； （3）与CAD/CAM等软件的集成； （4）提高自动化的网格处理能力； （5）软件面向专业用户的开放性； （6）软件开发强强联合。 4.发展趋势 二 一 二 1.基本原则 为正确、合理地建立有限元模型, 这里介绍划分网格时应考虑的一些基本原则: （1）网格数量：网格数量直接影响计算精度和计算时耗, 网格数量增加会提高计算精度, 但同时计算时耗也会增加。当网格数量较少时增加网格计算精度可明显提高, 但计算时耗不会有明显增加; 当网格数量增加到一定程度后, 再继续增加网格时精度提高就很小, 而计算时耗却大幅度增加。所以在确定网格数量时应权衡这两个因素综合考虑，如图 1 所示。  （2）网格密度：为了适应应力等计算数据的分布特点, 在结构不同部位需要采用大小不同的网格。从图 2 的有限元模型可以看出: 在孔的附近有集中应力, 因此网格需要加密; 周边应力梯度相对较小,网格划分较稀。该网格反映了疏密不同的网格划分原则: 在计算数据变化梯度较大的部位, 为了较好地反映数据变化规律, 需要采用比较密集的网格; 而在计算数据变化梯度较小的部位, 为减小模型规模, 网格则应相对稀疏。 （3）单元阶次： 单元阶次与有限元的计算精度有着密切的关联,单元一般具有线性、二次和三次等形式, 其中二次和传递给相邻单元。 （4）单元形状： 网格单元形状的好坏对计算精度有着很大的影响, 单元形状太差的网格甚至会中止计算。 （5）单元协调性 单元协调是指单元上的力和力矩能够通过节点传递给相邻单元。 2.网格生成通用方法 2.1 映射法： 映射法的基本思想是实际图形与标准图形的双向映射, 如图 3 所示, 具体有三个步骤:(1) 根据形体边界的参数方程, 利用适当的映射函数, 将待划分的物理域映射到参数空间, 形成规则参数域；(2) 对参数域进行网格划分；(3) 将参数空间内单元的网格反向映射到欧氏空间, 从而生成实际的网格。  2.2 基于栅格法： 基于栅格法 ( Grid- based Approach) 也叫空间分解法, 该算法的基本流程:首先用一组不相交的栅格 ( cells) 覆盖在物体之上, 既可在栅格的规则点处布置节点, 也可在栅格单元中随机布置节点; 再对栅格和物体进行相交检测, 保留完全或部分落在目标区域之内的栅格, 删除完全落在目标区域之外的栅格; 然后对与物体边界相交的栅格进行调整、剪裁、再分解等操作, 使其更准确地逼近目标区域; 最后对内部栅格和边界栅格进行栅格级的网格。 2.3 节点连元法： 最 早提出的节点生成法是随机布点法, 该方法不能保证布点的均匀性, 同时为了确定可用节点, 要进行点距检查和多次试探, 算法效率低。近年出现的长方形网格直接布点法、等距水平扫描线直接布点法等方法在一定程度上提高了节点连元法的计算效率, 但Delaunay 三角剖分方法 ( 简称 DT 法) 的出现才真正使得节点连元法成为目前主流的全自动网格划分方法之一。 2.