工程中的计算机模拟技术-有限元法与数学new

《工程中的计算机模拟技术-有限元法与数学new》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、工程中的计算机模拟技术-有限元法与数学机械系贺永近年来随着计算机技术的普及和计算速度的不断提高，有限元(finiteelement)分析在工程设计和分析中得到了越来越广泛的重视，已经成为解决复杂的工程分析计算问题的有效途径，现在从汽车到航天飞机几乎所有的设计制造都已离不开有限元分析计算，其在机械制造、材料加工、航空航天、汽车、土木建筑、电子电器，国防军工，船舶，铁道，石化，能源，科学研究等各个领域的广泛使用已使设计水平发生了质的飞跃，其主要表现在以下几个方面：增加产品和工程的可靠性；在产品的设计阶段

2、发现潜在的问题；经过分析计算，采用优化设计方案，降低原材料成本缩短产品投向市场的时间；模拟试验方案，减少试验次数，从而减少试验经费。图1、图2和图3为有限元的实际应用实例。图1、眼球运动特征分析（很多视网膜手术后眼角膜会受到损伤，而为了更好地了解其原因，美国波士顿医学院利用有限元软件对眼球进行建模分析，给出了应力的分布特征，从而精确地了解眼球在运动过程中的受力情况）图2、轮胎在零气压条件下的最大主应力云图（对于汽车来说，良好膨胀性能的轮胎不但可以节省汽油的消耗，而且可以提高汽车的可操作性以及乘客的舒

3、适性。现在一种新型的轮胎正被研究开发，这种轮胎可以在没有空气压力的情况下使汽车以正常的速度行驶一段有限的距离。为了更好地研究其性能，使用采用有限元分析软件对轮胎在行驶过程中的受力进行了详细的分析）图3牙齿的受力分析过程（建立CAD模型，划分网格，求解出牙齿受力的应力分布）有限元法的基本思想是将求解域空间离散成一系列的子空间，每个子空间称为单元，单元上有一定的节点。单元与单元之间通过单元上的节点来传递物理量（如力、热量、变形等），单元内部的物理量通过节点间的插值求出。根据胡克定律，可得出每个单元受力和

4、变形的关系式，将这些关系式按一定的规则整合，可将整个空间上问题的求解化为P=Kd这一基本形式，P为结构的总刚度矩阵，d为结构节点的位移向量，P为结构节点的载荷向量。由于划分单元的个数有限，节点的数目也有限，故此方法称为有限元法。用有限元方法解决问题时采用的是物理模型的近似。这种方法概念清晰，通用性与灵活性兼备，能妥善处理各种复杂情况，只要改变单元的数目就可以使解的精确度改变，得到与真实情况无限接近的解。有限元法的求解步骤如下：(1)定义求解域：根据实际问题近似确定求解域的物理性质和几何区域。(2)求

5、解域离散化：结构离散是有限元法的基础，也就是将分析对象按一定的规则划分成有限个具有不同大小和形状单元体的集合，使相邻单元在节点处连接，单元之间的载荷也仅由节点来传递，这一步骤习惯上称为有限元网格划分。(3)单元推导：在结构离散完成之后就可以对单元进行特性分析，建立各单元节点位移与节点力之间的关系，从而求出单元的刚度矩阵。(4)等效节点载荷计算。(5)总装求解。(6)联立方程组求解和结果解释。有限元法的理论基础可以追溯到变分法原理的提出。Bernoulli，Euler和Lagrange创立变分法，指出

6、从变分原理可以推出Newton力学。以后Hamilton提出著名的Hamilton原理：设T是体系的动能，D是体系的势能．L=T-D是体系的Lagrange函数，则从时刻t0到时刻t1的的真实运动使积分取最小值。热传导理论、电动力学、量子力学、广义相对论和量子场论都可以广泛应用这个原理。而有限元方法也是这个原理的产物。到了二十个世纪，几何学发展了无穷维空间的几何，在这种新的几何里，将变分问题转换成为几何问题。以Riemann，Hilbert等人为工作为基础，Friedrichs(1901-1982)

7、给出变分原理一个漂亮的理论，在这个理论发展的过程中，产生了Sobolev空间．Sobolev空间的几何是一种以能量为长度的几何，这种几何成为研究偏微分方程和有限元方法的重要工具。Courant等数学家在研究使用差分的方法解决椭圆边值问题时，将差分和变分结合起来，将区域剖分为许多三角形，利用函数在在各三角形顶点的区数值在每个三角形上插出一个函数，试图用这样的方法构造坐标函数，从而于1943年构造出了Courant三角形，该方法是最早的有限元模型。但由于计算量太大，缺乏有效快速的计算工具，该方法未能得到

8、推广。尽管Friedrichs，Sobolev、Ritz及Mikhlin等数学家的工作为有限元的数学理论做好了准备工作，随着后来电子计算机的设计与制造为有限元的实现提供了载体。但有限元方法的应用发展是由工程师们来完成的。从1955-1958年间，阿吉里斯(J．H．Argyris)等人研究复杂杆系结构，而进一步发展了结构分析的矩阵方法，在此基础上，他导出了平面应力状态下矩形板格的单元刚度矩阵。1956年特约(M．J．Turner)和克拉夫(B．w．Clough)等将刚架