有限元应用及我的认识

 
我在网上找了好久，甚至到维普中文文献数据库去找，终于找到一些关于有限元应用及其发展现状和前景的文章。感觉我的搜索能力还有待提高啊，但是，我总算大致对有限元在实际中的应用有了初步的了解。我下了很多别人用有限元处理实际问题的论文，从中我发现有限元的应用范围也是相当的广的。他涉及到工程结构、传热、流体运动、电磁等连续介质的力学分析中，并在气象、地球物理、医学等领域得到应用和发展。电子计算机的出现和发展是有限元法在许多实际问题中的应用变为现实，并具有广阔的前景。
人类认识客观世界的第一任务就是获取复杂对象的各类信息，这是人们从事科学研究、进行工程设计的基础。有限元分析（finite element analysisi）是进行科学计算的极为重要的方法之一，利用有限元分析可以获得几乎任意复杂工程结构的各种机械性能信息，还可以直接就工程设计进行各种评判，可以就各种工程事故进行技术分析。1990年10月美国波音公司开始在计算机上对新型客机B－777进行“无纸设计”，仅用了三年半时间，于1994年4月第一架B－777就试飞成功，这是制造技术史上划时代的成就，其中在结构设计和评判中就大量采用有限元分析这一手段。
在有限元分析的发展初期，由于其基本思想和原理的“简单”和“朴素”，以至于许多学术权威都对其学术价值有所鄙视，国际著名刊物Journal of Applied Mechanics 许多年来都拒绝刊登有关于有限元分析的文章。然而现在，有限元分析已经成为数值计算的主流，不但国际上存在如ANSYS等数种通用有限元分析软件，而且涉及到有限元分析的杂志也有几十种之多。
有限元法（FEA，Finite Element Analysis）的基本概念是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的工程分析手段。
有限元是那些集合在一起能够表示实际连续域的离散单元。有限元的概念早在几个世纪前就已产生并得到了应用，例如用多边形（有限个直线单元）逼近圆来求得圆的周长，但作为一种方法而被提出，则是最近的事。有限元法最初被称为矩阵近似方法，应用于航空器的结构强度计算，并由于其方便性、实用性和有效性而引起从事力学研究的科学家的浓厚兴趣。经过短短数十年的努力，随着计算机技术的快速发展和普及，有限元方法迅速从结构工程强度分析计算扩展到几乎所有的科学技术领域，成为一种丰富多彩、应用广泛并且实用高效的数值分析方法。
有限元方法与其他求解边值问题近似方法的根本区别在于它的近似性仅限于相对小的子域中。20世纪60年代初首次提出结构力学计算有限元概念的克拉夫（Clough）教授形象地将其描绘为：“有限元法=Rayleigh Ritz法＋分片函数”，即有限元法是Rayleigh Ritz法的一种局部化情况。不同于求解（往往是困难的）满足整个定义域边界条件的允许函数的Rayleigh Ritz法，有限元法将函数定义在简单几何形状（如二维问题中的三角形或任意四边形）的单元域上（分片函数），且不考虑整个定义域的复杂边界条件，这是有限元法优于其他近似方法的原因之一。
对于不同物理性质和数学模型的问题，有限元求解法的基本步骤是相同的，只是具体公式推导和运算求解不同。有限元求解问题的基本步骤通常为：
第一步：问题及求解域定义：根据实际问题近似确定求解域的物理性质和几何区域。
第二步：求解域离散化：将求解域近似为具有不同有限大小和形状且彼此相连的有限个单元组成的离散域，习惯上称为有限元网络划分。显然单元越小（网络越细）则离散域的近似程度越好，计算结果也越精确，但计算量及误差都将增大，因此求解域的离散化是有限元法的核心技术之一。
第三步：确定状态变量及控制方法：一个具体的物理问题通常可以用一组包含问题状态变量边界条件的微分方程式表示，为适合有限元求解，通常将微分方程化为等价的泛函形式。
第四步：单元推导：对单元构造一个适合的近似解，即推导有限单元的列式，其中包括选择合理的单元坐标系，建立单元试函数，以某种方法给出单元各状态变量的离散关系，从而形成单元矩阵（结构力学中称刚度阵或柔度阵）。
为保证问题求解的收敛性，单元推导有许多原则要遵循。对工程应用而言，重要的是应注意每一种单元的解题性能与约束。例如，单元形状应以规则为好，畸形时不仅精度低，而且有缺秩的危险，将导致无法求解。
第五步：总装求解：将单元总装形成离散域的总矩阵方程（联合方程组），反映对近似求解域的离散域的要求，即单元函数的连续性要满足一定的连续条件。总装是在相邻单元结点进行，状态变量及其导数（可能的话）连续性建立在结点处。
第六步：联立方程组求解和结果解释：有限元法最终导致联立方程组。联立方程组的求解可用直接法、选代法和随机法。求解结果是单元结点处状态变量的近似值。对于计算结果的质量，将通过与设计准则提供的允许值比较来评价并确定是否需要重复计算。
简言之，有限元分析可分成三个阶段，前处理、处理和后处理。前处理是建立有限元模型，完成单元网格划分；后处理则是采集处理分析结果，使用户能简便提取信息，了解计算结果。
除了找到了解到有限元的这些理论知识与应用外，我还找到了很多有限元相关的论坛，如SinWE仿真论坛、仿真在线、仿真科技论坛等等。从中我学到了很多以前从没听过的东西，下载了很多有关有限元的资料，让我受益匪浅。具体的与实际结合的我就不一一详细说明，我相信，随着对有限元的进一步学习，我对有限元的认识会进一步加深，将来我一定能够在实际中用好它。它的功能如此强大，前景是很美好的，是值得我们好好用心学习和研究的。
 

ps：翻出以前选修有限元课写的一篇小文章，粗看貌似比较粗糙啊，不过贴出来，呵呵。

原创博文，转载请注明出处：卓越替身 —http://www.chaself.com/
文章地址：http://www.chaself.com/post/173.html

标签:

发布:替身 | 分类:机械工程 | 评论:3 | 引用:0 | 浏览: | TrackBack引用地址