第6期
机电技术
135
作者简介：徐建全(1977－)，男，讲师，硕士，湖南大学机械与运载工程学院博士生；研究方向：计算机辅助工程分析，汽车车身结构与设
计理论，汽车技术与产业发展研究。
《有限元法》本科课程综合改革与实践
徐建全 陈铭年 林大同
(福建农林大学，福建 福州350002)
摘 要：论述了农林院校本科有限元法课程进行改革的必要性，介绍了综合改革的指导思想、目标和内容。实践表明，通过教学改革，学生的综合能力得到提高。
关键词：有限元法；教学改革；课程
中图分类号：G642.0 文献标识码：A 文章编号：1672-4801(2011)06-135-03
1 《有限元法》课程改革的必要性
有限元法是计算机辅助工程分析的主要方法之一，是工程结构(包括机械结构，汽车结构，建筑结构，航天航空结构等)设计不可缺少的先进方法。是否掌握有限元法已成为工程设计人员业务水平的重要标志之一。高等学校《有限元法》课程是机械、土木等工程专业的重要选修课之一(个别学校已作为必修课)，对于培养高素质、高质量的高级专门人才有着重要作用。有限元法涉及的理论较深，对数学和力学基础要求较高，对本科生而言，特别是农林院校的本科生而言，有一定难度，并不是人人通过自学都可以学会的。因此，有限元法课程的教学好坏就显得特别有意义。
目前各校讲授的《有限元法》课程及所编教材，大多以大量篇幅和课时讲授有限元法和各种单元的力学原理，在课堂花费很多时间进行数学力学推导。而用很少时间讲授应用，并一般不讲授实验验证方法。实践表明，这样的教学效果使得多数学生感觉深奥难懂，枯燥乏味且不懂应用。极少数学会有限元法的学生也不懂实验验证。
造成这种现状的原因是多方面的：首先是教师对本科生有限元法课程的教学应以应用为主这一点缺乏足够的认识；其次是选修课学时少，一般只有30学时。在力学学时普遍压缩的情况下，加上一些学生的数学和力学基础较薄弱，这么少的学时使学生学会有限元法有一定难度；最后是缺少适用的应用教材或实验条件。
因此，改革《有限元法》课程的教学，在有限的学时中，减少繁琐的推导，加强应用教学是很有必要的。而且，由于计算机使用的便利，有限元软件的具备，以及实验室条件的改善，已为有限元法的教学改革提供了可能。
2 改革的指导思想、目标和内容
几年来学校在车辆工程专业进行了有限元法课程综合改革的尝试，改革的指导思想是，农林院校车辆工程本科有限元法课程的教学应以应用为主[1]。改革的目标是在有限的学时内，使学生学会有限元法的基本理论，懂得有限元的应用和验证。改革的内容包括教学内容、教学方法，教学材料和实验教学[2]。
2.1 教学内容改革
改革后本科《有限元法》选修课为30学时，教学内容和学时安排如下：
第一章：结构力学基础(2学时)，主要介绍结构的分类，结构几何不变性，超静定，自由度，约束，对称结构的概念和性质等。
第二章：平面杆系结构分析(6学时)，主要介绍有限元法的基本原理，单元力学特性分析，坐标变换，结构整体分析。这一章是教学的重点，通过平面杆系教授有限元法的基本原理，学生学会基本原理后，对后续空间杆系的学习才会理解。因此对平面杆系结构分析进行详尽讲解是很有必要的。
第三章：空间杆系结构分析(3学时)，是将二维杆系结构分析的结果推广到三维。这一章着重介绍梁的空间坐标，空间梁单元的分析和应用。
第四章：平面问题的有限元分析(6学时)，主要介绍平面应力、应变问题，弹性体的离散化，位移模式的选择，收敛准则和实施步骤，平面单元的应用。其中位移模式是本章重点讲解内容。
第五章：轴对称单元、等参单元、板壳单元和实体单元等(3学时)。
机电技术2011年12月136
第六章：有限元软件ANSYS上机练习(8学时)。
第七章：实验验证(2学时)。
很多院校《有限元法》的授课首先是从弹性力学入手开始讲解，采用虚位移原理讲解单元刚度方程。而对于农林院校，太多弹性力学知识让学生普遍产生了畏难情绪，所以笔者采取从平面杆系结构分析入手，用直接法推导平面桁架的单元刚度矩阵，它所用的基础知识是在大二刚上过的《材料力学》，学生可以轻松理解，教学实践也证明了这一点。另外，增加了上机学时数，还特别增加了实验验证环节，让学生懂得如何通过实验来验证有限元模型的正确性。
2.2 教学方法
采用案例教学和启发式教学[3,4]，各章的教学方法略有不同。第一章结构力学基础，学时少，以介绍为主。第二章平面杆系结构分析，以讲有限元法基本原理为主，本章授课结束安排一个大作业加以巩固。第三章空间杆系结构分析，特别注意将二维概念推及到三维时，要给学生形象地说明。第四章对平面应力问题，位移模式选择要加以重点讲解，让学生深刻理解有限元的精髓之一插值思想，其它则以说明概念为主。第五章只进行简单介绍，重点讲解单元的基本特征及如何应用。理论教学采取多媒体和板书(一些简单推导)相结合的方法，充分发挥两者的优点，并紧密结合实际工程分析实例，大大调动了学生的积极性。每一章理论授课结束，安排一次上机，通过上机操作加深对理论知识的理解。
2.3 教学材料
编写了《有限元法应用》作为补充教材。补充教材共分六章，约6万字。第一章概述，包括前言，机械与汽车结构分析的过程，有限元程序，前、后处理；第二章机械与汽车结构的有限元计算模型，包括对计算模型的要求，机械结构的计算模型，支承与载荷；第三章有限元分析简单实例，包括分析步骤，平面刚架计算分析和空间刚架计算分析等；第四章ANSYS软件上机指南；第五章实验验证；第六章为车架综合实例。
2.4 实验教学
包括两部分，一部分是上机实验，一部分为实验验证。上机实验采用自编的ANSYS软件上机指南，包括桁架有限元分析、超静定粱有限元分析、平面支架有限元分析和连杆有限元分析等。其中有部分实例作为选作题，以便有兴趣的学生做更深入的实践。实验验证包括悬臂梁弯曲实验，L型结构的扭转实验如图1所示。国内高校的有限元法课程教学中普遍仅有上机实验，而缺乏实验验证内容。我们增加了实验验证内容，虽然仅有2学时的一次实验，但却使学生开拓了眼界，学到有限元法的实验验证方法和误差分析，这也是我们教改工作的主要特色。
为了更好的让学生真正掌握实际工程分析能力，作为拓展训练，老师给了一个实车车架(如图2所示)，让学生进行车架建模和有限元分析。通过真正的工程实例分析，增强了学生对有限元法的学习兴趣，提高了学生应用理论知识解决实际问题的能力，并为今后理论的深入学习和应用能
力的进一步提高打下基础。
图1 L型结构的扭转实验图
图2 车架实例
3 教学效果
从学生答卷、网上评教和课堂教学的效果来看，改革后的有限元法教学取得了良好的效果。主要表现在以下几个方面：
1) 原理听得懂。在平面杆系结构分析和位移模式选择中进行详尽的理论指导，都以讲清概念为主，减少了过多繁杂的推导，学生听得懂。
2) 方法会应用。由于树立了车辆工程本科生
第6期徐建全等：《有限元法》本科课程综合改革与实践137
的有限元法教学应以应用为主的指导思想，编写了应用讲义的辅助教学材料，在教学中注重应用，学生的应用能力得到了培养。也因为有应用，使得学生的兴趣大大增加。
3) 软件感兴趣。教学中所使用的ANSYS软件具有良好的人工界面，增加了学生学习兴趣。部分学生参加全国见习机械设计工程师资格认证考试，有限元上机操作部分平均通过率为98％。
4) 验证很直观。对计算结果的验证采用理论公式验证和实验验证两种方法。国内有限元法课程教学一般仅有上机实验而无实验验证内容，我们增设的实验验证内容在国内同类课程中属于领先。
总的来说，学生对有限元法课程教改是欢迎的，普遍反映听得懂、有兴趣、用得上。
4 结束语
通过几年来有限元法的综合改革与实践，学生初步掌握了有限元法的概念、基本原理和软件的操作，对有限元法有了更深刻的认识，学生的数值计算应用能力有了很大提高，为将来从事CAE相关研究工作打下了坚实的基础。
参考文献：
[1] 韦志林,尹辉俊,黄贵东.车辆专业有限元法课程的教学改革[J].中国科技信息,2008(3):225-226.
[2] 向家伟.机械类工程有限元法课程新体系的建设与实践[J].桂林电子科技大学学报,2008,28(2):150-152.
[3] 葛藤.《有限元法》本科教学改革探讨[J].科技信息,2010(13):116.
[4] 章泳健.面向应用的《有限元分析》课程的导入及教学实践[J].常熟理工学院学报,2009,23(6):116-118.
考试科目:《现代设计方法》 （总分100分） 时间：90分钟 __________学习中心（教学点） 批次： 层次： 专业： 学号： 身份证号： 姓名： 得分： 一、单项选择题（每小题1.5分，共27分） 1.试判别矩阵1111???? ? ?，它是（ ） A 、单位矩阵 B 、正定矩阵 C 、负定矩阵 D 、不定矩阵 2.约束极值点的库恩——塔克条件为：-?=?=∑F X g X i i q i ()()* * λ1 ，当约束函数是g i (X)≤0和 λi >0时，则q 应为（ ） A 、等式约束数目 B 、不等式约束数目 C 、起作用的等式约束数目 D 、起作用的不等式约束数目 3.在图示极小化的约束优化问题中，最优点为（ ） A 、A B 、B C 、C D 、D 4.下列优化方法中，不需计算迭代点一阶导数和二阶导数的是（ ） A 、可行方向法 B 、复合形法 C 、DFP 法 D 、BFGS 法 5.内点罚函数Φ(X,r (k) )=F(X)-r (k) 1 01g X g X u u u m () ,(())≤=∑，在其无约束极值点X ·(r (k))逼近原 目标函数的约束最优点时，惩罚项中（ ） A 、r (k) 趋向零， 11 g X u u m ()=∑ 不趋向零 B 、r (k) 趋向零，11g X u u m ()=∑ 趋向零 C 、r (k) 不趋向零， 11 g X u u m ()=∑ 趋向零 D 、④r (k) 不趋向零，11g X u u m ()=∑ 不趋向零 6.0.618法在迭代运算的过程中，区间的缩短率是（ ）
A 、不变的 B 、任意变化的 C 、逐渐变大 D 、逐渐变小 7.对于目标函数F(X)受约束于g u (X)≥0(u=1,2,…，m)的最优化设计问题，外点法惩罚函数的表 达式是（ ） A 、Φ(X,M (k) )=F(X)+M (k) {max[(),]},() g X M u u m k 012=∑为递增正数序列 B 、Φ(X,M (k))=F(X)+M (k){max[(),]},() g X M u u m k 012 =∑为递减正数序列 C 、Φ(X,M (k))=F(X)+M (k){min[(),]},()g x M u u m k 01 2 =∑为递增正数序列 D 、Φ(X,M (k))=F(X)+M (k){min[(),]},() g x M u u m k 01 2 =∑为递减正数序列 8.标准正态分布的均值和标准离差为（ ） A 、μ=1,σ=0 B 、μ=1,σ=1 C 、μ=0,σ=0 D 、μ=0,σ=1 9.在约束优化方法中，容易处理含等式约束条件的优化设计方法是（ ） A 、可行方向法 B 、复合形法 C 、内点罚函数法 D 、外点罚函数法 10.若组成系统的诸零件的失效相互独立，但只有某一个零件处于工作状态，当它出现故障后， 其它处于待命状态的零件立即转入工作状态。这种系统称为（ ） A 、串联系统 B 、工作冗余系统 C 、非工作冗余系统 D 、r/n 表决系统 11.对于二次函数F(X)=1 2 X T AX+b T X+c,若X *为其驻点，则▽F(X *)为（ ） A 、零 B 、无穷大 C 、正值 D 、负值 12.平面应力问题中(Z 轴与该平面垂直)，所有非零应力分量均位于（ ） A 、XY 平面内 B 、XZ 平面内 C 、YZ 平面内 D 、XYZ 空间内 13当选线长度l ，弹性模量E 及密度ρ为三个基本量时，用量纲分析法求出包含振幅A 在内的 相似判据为(E 的量纲为（ ）［ML -1T -2 ］ A 、A=l E 1 1212- ρ B 、A=l E -- 1 121 2 ρ C 、A=l E 100ρ D 、A l E =-11 12ρ 14.平面三角形单元内任意点的位移可表示为三个节点位移的（ ） A 、算术平均值 B 、代数和车员 C 、矢量和 D 、线性组合 15.已知F(X)=(x 1-2)2+x 22,则在点X (0)=00???? ??处的梯度为（ ） A 、?=?????? F X ()()000 B 、?=-?????? F X ()() 020
有限单元法与有限元分析 1.有限单元法 在数学中，有限元法（FEM，Finite Element Method）是一种为求解偏微分方程边值问题近似解的数值技术。求解时对整个问题区域进行分解，每个子区域都成为简单的部分，这种简单部分就称作有限元。它通过变分方法，使得误差函数达到最小值并产生稳定解。类比于连接多段微小直线逼近圆的思想，有限元法包含了一切可能的方法，这些方法将许多被称为有限元的小区域上的简单方程联系起来，并用其去估计更大区域上的复杂方程。它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的工程分析手段。 随着电子计算机的发展，有限单元法是迅速发展成一种现代计算方法。它是50年代首先在连续体力学领域--飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。 1.1.有限元法分析本质 有限元法分析计算的本质是将物体离散化。即将某个工程结构离散为由各种单元组成的计算模型，这一步称作单元剖分。离散后单元与单元之间利用单元的节点相互连接起来；单元节点的设置、性质、数目等应视问题的性质，描述变形形态的需要和计算精度而定（一般情况单元划分越细则描述变形情况越精确，即越接近实际变形，但计算量越大）。所以有限元中分析的结构已不是原有的物体或结构物，而是同新材料的由众多单元以一定方式连接成的离散物体。这样，用有限元分析计算所获得的结果只是近似的。如果划分单元数目非常多而又合理，则所获得的结果就与实际情况相符合。 1.2.特性分析 1）选择位移模式： 在有限单元法中，选择节点位移作为基本未知量时称为位移法；选择节点力作为基本未知量时称为力法；取一部分节点力和一部分节点位移作为基本未知量时称为混合法。位移法易于实现计算自动化，所以，在有限单元法中位移法应用范围最广。 当采用位移法时，物体或结构物离散化之后，就可把单元总的一些物理量如
目录 摘要 (1) 1引言 (1) 2 简要介绍有限元和极限平衡方法 (1) 3影响边坡稳定性的因素 (2) 3.1水位下降速度的影响 (2) 3.2 不排水粘性土对边坡失稳的影响 (5) 3.3 裂缝位置的影响 (9) 4 总结和结论 (12)
基于有限元法和极限平衡法的边坡稳定性分析 摘要：相较于有限元分析法，极限平衡法是一种常用的更为简单的边坡稳定性分析方法。这两种方法都可用于分析均质和不均质的边坡，同时考虑了水位骤降，饱和粘土和存在张力裂缝的条件。使用PLAXIS8.0(有限元法)和SAS-MCT4.0(极限平衡方法)进行了分析，并对两种方法获得的临界滑动面的安全系数和位置进行了比较。 关键词：边坡稳定；极限平衡法；有限元法；PLAXIS；SAS-MCT 1.引言 近年来，计算方法，软件设计和高速低耗硬件领域都得到快速发展，特别是相关的边坡稳定性分析的极限平衡法和有限元方法。但是，使用极限平衡方法来分析边坡，可能会在定位临界滑动面(取决于地质)时出现几个计算困难和前后数值不一致，因此要建立一个安全系数。尽管极限平衡法存在这些固有的局限性，但由于其简单，它仍然是最常用的方法。然而，由于个人电脑变得更容易获得，有限元方法已越来越多地应用于边坡稳定性分析。有限元法的优势之一是，不需要假设临界破坏面的形状或位置。此外，该方法可以很容易地用于计算压力，位移，路堤空隙压力，渗水引起的故障，以及监测渐进破坏。 邓肯(1996年)介绍了一个综合观点，用极限平衡和有限元两种方法对边坡进行分析。他比较了实地测量和有限元分析的结果，并且发现一种倾向，即计算变形大于实测变形。Yu 等人(1998年)比较了极限平衡法和严格的上、下界限法对于简单土质边坡的稳定性分析的结果，同时，他们也将采用毕肖普法和利用塑性力学上、下限原理的界限法得到的结果进行了比较。Kim等人(1999年)同时使用极限平衡法和极限分析法对边坡进行分析，发现对于均质土边坡，得自两种方法的结果大体是一致的，但是对于非均质土边坡还需要进行进一步分析工作。Zaki(1999年)认为有限元相对于极限平衡法更显优势。Lane和Griffiths (2000年) 提出一个看法，用有限元方法在水位骤降条件下评价边坡的稳定性，应绘制出适用于实际结构的操作图表。Rocscience有限公司(2001年)提出了一个文件，概述了有限元分析方法的能力，并通过与各种极限平衡方法的结果比较，提出了有限元方法更为实用。Kim等人(2002年)用上、下界限法和极限平衡法分析了几处非均质土体且几何不规则边坡的剖面。这两种方法给出了类似有限元分析法产生的安全系数，临界滑动面位置。 2.简要介绍有限元和极限平衡方法 有限元法(FEM)是一个应用于科学和工程中，求解微分方程和边值问题的数值方法。进一步的细节，读者可参考Clough和Woodward(1967年)，Strang和Fix(1973年)，Hughes(1987年)，Zienkiewicz和Taylor(1989年)所做的研究工作。 PLAXIS 8版(Brinkgreve 2002年)是一个有限元软件包，应用于岩土工程二维的变形和 折稳定性分析。该程序可以分析自然成型或人为制造的斜坡问题。安全系数的确定使用c
有限差分法、有限单元法和有限体积法的简介 1.有限差分方法 有限差分方法(Finite Difference Method，FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。 对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑，可分为中心格式和逆风格式。考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 2.有限元方法 有限元方法(Finite Element Method，FEM)的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。 有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。 在数值模拟中，常见的有限元计算方法是由变分法和加权余量法发展而来的
有限元分析概念 有限元法：把求解区域看作由许多小的在节点处相互连接的单元（子域）所构成，其模型给出基本方程的分片（子域）近似解，由于单元（子域）可以被分割成各种形状和大小不同的尺寸，所以它能很好地适应复杂的几何形状、复杂的材料特性和复杂的边界条件 有限元模型：它是真实系统理想化的数学抽象。由一些简单形状的单元组成，单元之间通过节点连接，并承受一定载荷。 有限元分析：是利用数学近似的方法对真实物理系统（几何和载荷工况）进行模拟。并利用简单而又相互作用的元素，即单元，就可以用有限数量的未知量去逼近无限未知量的真实系统。 线弹性有限元是以理想弹性体为研究对象的，所考虑的变形建立在小变形假设的基础上。在这类问题中，材料的应力与应变呈线性关系，满足广义胡克定律；应力与应变也是线性关系，线弹性问题可归结为求解线性方程问题，所以只需要较少的计算时间。如果采用高效的代数方程组求解方法，也有助于降低有限元分析的时间。 线弹性有限元一般包括线弹性静力学分析与线弹性动力学分析两方面。 非线性问题与线弹性问题的区别： 1）非线性问题的方程是非线性的，一般需要迭代求解； 2）非线性问题不能采用叠加原理； 3）非线性问题不总有一致解，有时甚至没有解。 有限元求解非线性问题可分为以下三类：
1）材料非线性问题 材料的应力和应变是非线性的，但应力与应变却很微小，此时应变与位移呈线性关系，这类问题属于材料的非线性问题。由于从理论上还不能提供能普遍接受的本构关系，所以，一般材料的应力与应变之间的非线性关系要基于试验数据，有时非线性材料特性可用数学模型进行模拟，尽管这些模型总有他们的局限性。在工程实际中较为重要的材料非线性问题有：非线性弹性（包括分段线弹性）、弹塑性、粘塑性及蠕变等。 2）几何非线性问题 几何非线性问题是由于位移之间存在非线性关系引起的。 当物体的位移较大时，应变与位移的关系是非线性关系。研究这类问题一般都是假定材料的应力和应变呈线性关系。它包括大位移大应变及大位移小应变问题。如结构的弹性屈曲问题属于大位移小应变问题，橡胶部件形成过程为大应变问题。 3）非线性边界问题 在加工、密封、撞击等问题中，接触和摩擦的作用不可忽视，接触边界属于高度非线性边界。 平时遇到的一些接触问题，如齿轮传动、冲压成型、轧制成型、橡胶减振器、紧配合装配等，当一个结构与另一个结构或外部边界相接触时通常要考虑非线性边界条件。 实际的非线性可能同时出现上述两种或三种非线性问题。
1.1 概念 有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。 1.2 差分格式 （1）从格式的精度来划分，有一阶格式、二阶格式和高阶格式。 （2）从差分的空间形式来考虑，可分为中心格式和逆风格式。 （3）考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。 目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。 1.3 构造差分的方法 构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 2. FEM 2.1 概述 有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。 2.2 原理 有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学、土力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。在河道数值模拟中，常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。 根据所采用的权函数和插值函数的不同，有限元方法也分为多种计算格式。 （1）从权函数的选择来说，有配置法、矩量法、最小二乘法和伽辽金法； （2）从计算单元网格的形状来划分，有三角形网格、四边形网格和多边形网格； （3）从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。 不同的组合同样构成不同的有限元计算格式。
中国地质大学研究生课程论文封面 课程名称现代设计方法 教师姓名 研究生姓名 研究生学号 研究生专业机械工程 所在院系机电学院 日期: 2013 年 1 月 8 日
评语 注：1、无评阅人签名成绩无效； 2、必须用钢笔或圆珠笔批阅，用铅笔阅卷无效； 3、如有平时成绩，必须在上面评分表中标出，并计算入总成绩。
有限元分析简介 摘要: ANSYS 软件具有建模简单、快速、方便的特点, 因而成为大型通用有限元程序的代表。对有限元作了一个总体的介绍, 并着重介绍了ANSYS 软件, 简要地叙述了ANSYS 软件的主要技术特点和各部分构成以及其主要的分析功能,从其构成及功能中可以看到,ANSYS 软件的确是工程应用分析的有效工具。 1、有限元分析的基本概念和计算步骤 1.1、有限元分析的基本概念 有人将CAE技术称为当今“科学与技术的完美结合”。这句话说得比较夸张，但不可否认，CAE技术的确是现代产品研发的重要基础技术，其理论性和需要的学科知识厚重而宽广。有限元软件是目前CAE的主流分析软件之一，在全球拥有最大的用户群。有限元分析（FEA，Finite Element Analysis）利用数学近似的方法对真实物理系统（几何和载荷工况）进行模拟。还利用简单而又相互作用的元素，即单元，就可以用有限数量的未知量去逼近无限未知量的真实系统。 有限元分析的基本概念是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的工程分析手段。 有限元是那些集合在一起能够表示实际连续域的离散单元。有限元的概念早在几个世纪前就已产生并得到了应用，例如用多边形（有限个直线单元）逼近圆来求得圆的周长，但作为一种方法而被提出，则是最近的事。有限元法最初被称为矩阵近似方法，应用于航空器的结构强度计算，并由于其方便性、实用性和有效性而引起从事力学研究的科学家的浓厚兴趣。经过短短数十年的努力，随着计算机技术的快速发展和普及，有限元方法迅速从结构工程强度分析计算扩展到几乎所有的科学技术领域，成为一种丰富多彩、应用广泛并且实用高效的数值分析方法。
有限差分和有限体积的有限元等 有限元法、有限差分法和有限体积法的区别 标签：函数有限元插值差分格式 有限差分方法(Finite Differential Method)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以泰勒级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。 对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑，可分为中心格式和逆风格式。考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。 构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 有限元法(Finite Element Method)的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。 根据所采用的权函数和插值函数的不同，有限元方法也分为多种计算格式。从权函数的选择来说，有配置法、矩量法、最小二乘法和伽辽金法。从计算单元网格的形状来划分，有三角形网格、四边形网格和多边形网格，从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。不同的组合同样构成不同的有限元计算格式。对于权函数，伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数；最小二乘法是令权函数等于余量本身，而内积的极小值则为对代求系数的平方误差最小；在配置法中，先在计算域内选取N个配置点。令近似解在选定的N个配置点上严格满足微分方程，即在配置点上令方程余量为0。插值函数一般由不同次幂的多项式组成，但也有采用三角函数或指数函数组成的乘积表示，但最常用的多项式插值函数。有限元插值函数分为两大类，一类只要求插值多项式本身在插值点取已知值，称为拉格朗日(Lagrange)多项式插值；另一种不仅要求插值多项式本身，还要求它的导数值在插值点取已知值，称为哈密特(Hermite)多项式插值。单元坐标有笛卡尔直角坐标系和无因次自然坐标，有对称和不对称等。常采用的无因次坐标是一种局部坐标系，它的定义取决于单元的几何形状，一维看作长度比，二维看作面积比，三维看作体积比。在二维有限元中，三角形单元应用的最早，近来四边形等参元的应用也越来越广。对于二维三角形和四边形电源单元，常采用的插值函数为有Lagrange插值直角坐标系中的线性插值函
有限差分方法(Finite Differential Method)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以泰勒级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。 对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑，可分为中心格式和逆风格式。考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。 构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 有限元法(Finite Element Method)的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。 根据所采用的权函数和插值函数的不同，有限元方法也分为多种计算格式。从权函数的选择来说，有配置法、矩量法、最小二乘法和伽辽金法。从计算单元网格的形状来划分，有三角形网格、四边形网格和多边形网格，从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。不同的组合同样构成不同的有限元计算格式。对于权函数，伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数；最小二乘法是令权函数等于余量本身，而内积的极小值则为对代求系数的平方误差最小；在配置法中，先在计算域内选取N个配置点。令近似解在选定的N个配置点上严格满足微分方程，即在配置点上令方程余量为0。插值函数一般由不同次幂的多项式组成，但也有采用三角函数或指数函数组成的乘积表示，但最常用的多项式插值函数。有限元插值函数分为两大类，一类只要求插值多项式本身在插值点取已知值，称为拉格朗日(Lagrange)多项式插值；另一种不仅要求插值多项式本身，还要求它的导数值在插值点取已知值，称为哈密特(Hermite)多项式插值。单元坐标有笛卡尔直角坐标系和无因次自然坐标，有对称和不对称等。常采用的无因次坐标是一种局部坐标系，它的定义取决于单元的几何形状，一维看作长度比，二维看作面积比，三维看作体积比。在二维有限元中，三角形单元应用的最早，近来四边形等参元的应用也越来越广。对于二维三角形和四边形电源单元，常采用的插值函数为有Lagrange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线性插值函数、二阶或更高阶插值函数等。 有限体积法(Finite V olume Method)又称为控制体积法。其基本思路是：将计算区域划分为一系列不重复的控制体积，并使每个网格点周围有一个控制体积；将待解的微分方程对每一个控制体积积分，便得出一组离散方程。其中的未知数是网格点上的因变量的数值。为了求出控制体积的积分，必须假定值在网格点之间的变化规律，即假设值的分段的分布的分布剖面。从积分区域的选取方法看来，有限体积法属于加权剩余法中的子区域法；从未知解的近似方法看来，有限体积法属于采用局部近似的离散方法。简言之，子区域法属于有限体积发的基本方法。
有限元法与有限差分法的主要区别 有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Ｔaylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观，表达简单，是发展较早且比较成熟的数值方法。对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑,可分为中心格式和逆风格式．考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等.目前常见的差分格式，主要是上述几种形式的组合,不同的组合构成不同的差分格式。差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内，选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成．在河道数值模拟中,常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等.根据所采用的权函数和插值函数的不同,有限元方法也分为多种计算格式。从权函数的选择来说,有配置法、矩量法、最小二乘法和伽辽金法,从计算单元网格的形状来划分,有三角形网格、四边形网格和多边形网格，从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。不同的组合同样构成不同的有限元计算格式。对于权函数，伽辽金(Gaｌeｒkiｎ）法是将权函数取为逼近函数中的基函数；最小二乘法是令权函数等于余量本身，而内积的极小值则为对代求系数的平方误差最小;在配置法中，先在计算域内选取N个配置点。令近似解在选定的Ｎ个配置点上严格满足微分方程，即在配置点上令方程余量为０．插值函数一般由不同次幂的多项式组成,但也有采用三角函数或指数函数组成的乘积表示,但最常用的多项式插值函数。有限元插值函数分为两大类,一类只要求插值多项式本身在插值点取已知值,称为拉格朗日(Lａgranｇe)多项式插值;另一种不仅要求插值多项式本身，还要求它的导数值在插值点取已知值,称为哈密特(Ｈｅrmite)多项式插值。单元坐标有笛卡尔直角坐标系和无因次自然坐标，有对称和不对称等。常采用的无因次坐标是一种局部坐标系，它的定义取决于单元的几何形状,一维看作长度比，二维看作面积比，三维看作体积比。在二维有限元中，三角形单元应用的最早，近来四边形等参元的应用也越来越广。对于二维三角形和四边形电源单元,常采用的插值函数为有Ｌa g ｒange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线性插值函数、二阶或更高阶插值函数等．对于有限元方法，其基本思路和解题步骤可归纳为(１）建立积分方程,根据变分原理或方程余量与权函数正交化原理,建立与微分方程初边值问题等价的积分表达式,这是有限元法的出发点。(2）区域单元剖分，根据求解区域的形状及实际问题的物理特点，将区域剖分为若干相互连接、不重叠的单元。区域单元划分是采用有限元方法的前期准备工作,这部分工作量比较大，除了给计算单元和节点进行编号和确定相互之间的关系之外,还要表示节点的位置坐标，同时还需要列出自然边界和本质边界的节点序号和相应的边界值。(3）确定单元基函数，根据单元中节点数目及对近似解精度的要求，选择满足一定插值条件的插值函
有限差分，有限元，有限体积等离散方法的区别介绍 1 有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑，可分为中心格式和逆风格式。考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。 构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 2 有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。 在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。 根据所采用的权函数和插值函数的不同，有限元方法也分为多种计算格式。从权函数的选择来说，有配置法、矩量法、最小二乘法和伽辽金法，从计算单元网格的形状来划分，有三角形网格、四边形网格和多边形网格，从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。不同的组合同样构成不同的有限元计算格式。 对于权函数，伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数；最小二乘法是令权函数等于余量本身，而内积的极小值则为对代求系数的平方误差最小；在配置法中，先在计算域内选取N个配置点。令近似解在选定的N个配置点上严格满足微分方程，即在配置点上令方程余量为0。插值函数一般由不同次幂的多项式组成，但也有采用三角函数或指数函数组成的乘积表示，但最常用的多项式插值函数。有限元插值函数分为两大类，一类只要求插值多项式本身在插值点取已知值，称为拉格朗日(Lagrange)多项式插值；另一种不仅要求插值多项式本身，还要求它的导数值在插值点取已知值，称为哈密特(Hermite)多项式插值。单元坐标有笛卡尔直角坐标系和无因次自然坐标，有对称和不对称等。常采用的无因次坐标是一种局部坐标系，它的定义取决于单元的几何形状，一维看作长度比，二维看作面积比，三维看作体积比。在二维有限元中，三角形单元应用的最早，近来四边形等参元的应用也越来越广。对于二维三角形和四边形电源单元，常采用的插值函数为有La grange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线性插值函数、二阶或更高
https://www.doczj.com/doc/2f4500105.html,/s/blog_501a61220100f9rs.html 有限差分，有限元，有限体积等离散方法的区别介绍 (2009-10-25 22:07:18) 转载 以下介绍是本人从网络上搜集的，供计算数学虫子参考。也许小木虫论坛有，我没搜索到。欢迎大家补充内容。 转自https://www.doczj.com/doc/2f4500105.html,/bbs/viewthread.php?tid=1618917&pid=16196206&page=1#pid16196206 1 有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑，可分为中心格式和逆风格式。考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。 构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 2 有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。 在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。 根据所采用的权函数和插值函数的不同，有限元方法也分为多种计算格式。从权函数的选择来说，有配置法、矩量法、最小二乘法和伽辽金法，从计算单元网格的形状来划分，有三角形网格、四边形网格和多边形网格，从插值函数的精度来划分，又分
《现代设计方法》作业关于有限元法的研究 学院：机械工程学院 专业：机械制造及其自动化
0.有限元法 有限元法分析起源于50年代初杆系结构矩阵的分析。随后，Clough于1960年第一次提出了“有限元法”的概念。其基本思想是利用结构离散化的概念，将连续介质体或复杂结构体划分成许多有限大小的子区域的集合体，每一个子区域称为单元（或元素），单元的集合称为网格，实际的连续介质体（或结构体）可以看成是这些单元在它们的节点上相互连接而组成的等效集合体；通过对每个单元力学特性的分析，再将各个单元的特性矩阵组集成可以建立整体结构的力学方程式，即力学计算模型；按照所选用计算程序的要求，输入所需的数据和信息，运用计算机进行求解。 当前，有限元方法/理论已经发展的相当成熟和完善，而计算机技术的不断革新，又在很大程度上推进了有限元法分析在工程技术领域的应用。然而，如此快速地推广和应用使得人们很容易忽视一个前提，即有限元分析软件提供的计算结果是否可靠、满足使用精度的前提，是合理地使用软件和专业的工程分析。有限元法分析一般包括四个步骤：物理模型的简化、数学模型的程序化、计算模型的数值化和计算结果的分析。每一个步骤在操作过程中都或多或少地引入了误差，这些误差的累积最终可能会对计算结果造成灾难性的影响，进而蒙蔽我们的认识和判断。 1.受内压空心圆筒的轴对称有限元分析 例图1.1所示为一无限长的受内压的轴对称圆筒，该圆筒置于内径为120mm的刚性圆孔中，试求圆筒内径处的位移。结构的材料参数
为：200 =，0.3 E GPa μ=。 图1 结构图 对该问题进行有限元分析的过程如下。 （1）结构的离散化与编号 由于该圆筒为无限长，取出中间一段（20mm高），采用两个三角形轴对称单元，如图1.2所示。对该系统进行离散，单元编号及结点编号如图1.3所示，有关结点和单元的信息见表1.1。 图1.2 有限元模型
有限元极限载荷分析法在压力容器分析设计中的应用2010-07-15 10:39:54| 分类：分析设计| 标签：极限分析分析设计asme规范先进设计方法经验分享|字号大 中 小订阅 在某炼化一体化项目中，几个加氢反应器均采用分析法设计。详细设计时，国内计算后，反应器的主要受压元件厚度均要比专利商建议的厚度多出10~30mm不等。这其中有国内设计出于保守的考虑，另一个原因：同是采用分析设计，ASME的非线性分析相对先进一点。参与国际竞争时，先进的设计方法值得我们研究。 1.背景 随着中国加入WTO，国内各工程公司正积极走向海外。随之进入国际市场的压力容器产品也面临着严峻的挑战，为了在国际舞台上获得竞争优势，各工程公司必须采用先进的技术设计出更安全和更低成本的产品。压力容器分析设计是力学与工程紧密结合产物，解决了常规设计无法解决的问题，代表了近代设计的先进水平[1]。过去，国内分析设计通常采用弹性应力分析法，通过路径分析，应力线性化处理获得路径上的一次应力和二次应力，进而进行强度评定。该方法主要存在以下问题：⑴对大多数情况是安全可靠的，但对某些结果可能出现安全裕度不足的情况（如球壳开打孔）；⑵如何对有限元法求解获得的总应力分解并正确分类遇到了困难。假如把一次应力误判为二次，则设计的结果将非常危险，反之，把二次应力误判为一次，则又非常保守。文[2]5.2.1.2节明确提到：应力分类需特殊的知识和识别力，应力分类方法可能产生模棱两可的结果。国内专家亦也认为对应力进行正确的分类存在一定困难[3-6]。 以弹性分析代替塑性分析,是一种工程近似方法。实际结构的破坏往往是一个渐进过程，随着载荷的增加，高应力区首先进入屈服，载荷继续增加时塑性区不断夸大，同时出现应力重新分布。当载荷增大到某一值时，结构变为几何可变机构，此时即使载荷不在增加，变形也会无限增大，发生总体塑性变形（overall plastic deformation），此时的载荷称为“极限载荷（limit load）”。 极限载荷分析法（下文简称极限分析）的目的是求出结构的极限载荷。在防止塑性垮塌失效时，极限分析相比弹性应力分析更接近工程实际，同时避免了应力分类，对防止塑性垮塌有比较精确的评定。 2．极限载荷的求解方法 塑性力学提出极限分析法由来已久。经典的极限分析方法有如下3种[8]：(1)广义内力与广义变形法；(2)上限定理与下限定理法；(3)静力法和机动法。经典解法的分析与计算均很复杂，只能应用于少数结构简单的压力容器元件，从而使极限分析的工程应用受到了限制。 上世纪七十年代出现三维有限元计算后，有限元的应用大大扩展。为了适应工程需要，有限元极限分析应运而生，形成了分析设计中的一个重要分支，它使得复杂的塑性极限分析可以通过计算机数值计算得以解决。在不久的将来，极限分析必与弹性应力分析法、弹-塑性应力分析法一同形成三足鼎立之势。极限分析的模型精度和计算成本居后两者之间。