有限元
百科名片
有限元法（FEA，Finite Element Analysis）的基本概念是用较简单的问题代替复杂问题后
再求解。它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的工程分析手段。
目录
简介
1）物体离散化
2）单元特性分析
3）单元组集
4）求解未知节点位移
5）有限元的未来是多物理场耦合
编辑本段简介
英文：Finite Element
有限单元法是随着电子计算机的发展而迅速发展起来的一种现代计算方法。它是50年代首先在连续体力学领域--飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。
有限元法分析计算的思路和做法可归纳如下：
编辑本段1）物体离散化
将某个工程结构离散为由各种单元组成的计算模型，这一步称作单元剖分。离散后单元与单元之间利用单元的节点相互连接起来；单元节点的设置、性质、数目等应视问题的性质，描述变形形态的需要和计算进度而定（一般情况单元划分越细则描述变形情况越精确，即越接近实际变形，但计算量越大）。所以有限元中分析的结构已不是原有的物体或结构物，而是同新材料的由众多单元以一定方式连接成的离散物体。这样，用有限元分析计算所获得的结果只是近似的。如果划分单元数目非常多而又合理，则所获得的结果就与实际情况相符合。
编辑本段2）单元特性分析
A、选择位移模式
在有限单元法中，选择节点位移作为基本未知量时称为位移法；选择节点力作为基本未知量时称为力法；取一部分节点力和一部分节点位移作为基本未知量时称为混合法。位移法易于实现计算自动化，所以，在有限单元法中位移法应用范围最广。
当采用位移法时，物体或结构物离散化之后，就可把单元总的一些物理量如位移，应变和应力等由节点位移来表示。这时可以对单元中位移的分布采用一些能逼近原函数的近似函数予以描述。通常，有限元法我们就将位移表示为坐标变量的简单函数。这种函数称为位移模式或位移函数。
B、分析单元的力学性质
根据单元的材料性质、形状、尺寸、节点数目、位置及其含义等，找出单元节点力和节点位移的关系式，这是单元分析中的关键一步。此时需要应用弹性力学中的几何方程和物理方程来建立力和位移的方程式，从而导出单元刚度矩阵，这是有限元法的基本步骤之一。
C、计算等效节点力
物体离散化后，假定力是通过节点从一个单元传递到另一个单元。但是，对于实际的连续体，力是从单元的公共边传递到另一个单元中去的。因而，这种作用在单元边界上的表面力、体积力和集中力都需要等效的移到节点上去，也就是用等效的节点力来代替所有作用在单元上的力。
编辑本段3）单元组集
利用结构力的平衡条件和边界条件把各个单元按原来的结构重新连接起来，形成整体的有限元方程
（1-1）
式中，K是整体结构的刚度矩阵；q是节点位移列阵；f是载荷列阵。
编辑本段4）求解未知节点位移
解有限元方程式（1-1）得出位移。这里，可以根据方程组的具体特点来选择合适的计算方法。
通过上述分析，可以看出，有限单元法的基本思想是"一分一合"，分是为了就进行单元分析，合则为了对整体结构进行综合分析。
有限元的发展概况
1943年courant在论文中取定义在三角形域上分片连续函数，利用最小势能原理研究St.Venant的扭转问题。
1960年clough的平面弹性论文中用“有限元法”这个名称。
1965年冯康发表了论文“基于变分原理的差分格式”，这篇论文是国际学术界承认我国独立发展有限元方法的主要依据。
1970年随着计算机和软件的发展，有限元发展起来。
涉及的内容：有限元所依据的理论，单元的划分原则，形状函数的选取及协调性。
有限元法涉及：数值计算方法及其误差、收敛性和稳定性。
应用范围：固体力学、流体力学、热传导、电磁学、声学、生物力学
求解的情况：杆、梁、板、壳、块体等各类单元构成的弹性（线性和非线性）、弹塑性或塑性问题（包括静力和动力问题）。能求解各类场分布问题（流体场、温度场、电磁场等的稳态和瞬态问题），水流管路、电路、润滑、噪声以及固体、流体、温度相互作用的问题。
编辑本段5）有限元的未来是多物理场耦合
5）有限元的未来是多物理场耦合
随着计算机技术的迅速发展，在工程领域中，有限元分析（FEA）越来越多地用于仿真模拟，来求解真实的工程问题。这些年来，越来越多的工程师、应用数学家和物理学家已经证明这种采用求解偏微分方程（PDE）的方法可以求解许多物理现象，这些偏微分方程可以用来描述流动、电磁场以及结构力学等等。有限元方法用来将这些众所周知的数学方程转化为近似的数字式图象。
早期的有限元主要关注于某个专业领域，比如应力或疲劳，但是，一般来说，物理现象都不是单独存在的。例如，只要运动就会产生热，而热反过来又影响一些材料属性，如电导率、化学反应速率、流体的粘性等等。这种物理系统的耦合就是我们所说的多物理场，分析起来比我们单独去分析一个物理场要复杂得多。很明显，我们现在需要一个多物理场分析工具。
在上个世纪90年代以前，由于计算机资源的缺乏，多物理场模拟仅仅停留在理论阶段，有限元建模也局限于对单个物理场的模拟，最常见的也就是对力学、传热、流体以及电磁场的模拟。看起来有限元仿真的命运好像也就是对单个物理场的模拟。
现在这种情况已经开始改变。经过数十年的努力，计算科学的发展为我们提供了更灵巧简洁而又快速的算法，更强劲的硬件配置，使得对多物理场的有限元模拟成为可能。新兴的有限元方法为多物理场分析提供了一个新的机遇，满足了工程师对真实物理系统的求解需要。有限元的未来在于多物理场求解。
千言万语道不尽，下面只能通过几个例子来展示多物理场的有限元分析在未来的一些潜在应用。
压电扩音器（Piezoacoustic transducer）可以将电流转换为声学压力场，或者反过来，将声场转换为电流场。这种装置一般用在空气或者液体中的声源装置上，比如相控阵麦克风，超声生物成像仪，声纳传感器，声学生物治疗仪等，也可用在一些机械装置比如喷墨机和压电马达等。
压电扩音器涉及到三个不同的物理场：结构场，电场以及流体中的声场。只有具有多物理场分析能力的软件才能求解这个模型。
压电材料选用PZT5-H晶体，这种材料在压电传感器中用得比较广泛。在空气和晶体的交界面处，将声场边界条件设置为压力等于结构场的法向加速度，这样可以将压力传到空气中去。另外，晶体域中又会因为空气压力对其的影响而产生变形。仿真研究了在施加一个幅值200V，震荡频率为300 KHz的电流后，晶体产生的声波传播。
这个模型的描述及其完美的结果表明在任何复杂的模型下，我们都可以用一系列的数学模型进行表达，进而求解。
多物理场建模的另外一个优势就是在学校里，学生们直观地获取了以前无法见到的一些现象，而简单易懂的表达方式也获得了学生们的好感。这只是Krishan Kumar Bhatia博士在纽约Glassboro的Rowan 大学给高年级的毕业生讲授传热方程课程时介绍建模及分析工具所感受到的，他的学生的课题是如何冷却一个摩托车的发动机箱。Bhatia博士教他们如何利用“设计－制造－检测”的理念来判断问题、找出问题、解决问题。如果没有计算机仿真的应用，这种方法在课堂上推广是不可想象的，因为所需费用实在是太大了。
COMSOL Multiphysics拥有优秀的用户界面，可以使学生方便地设置传热问题，并很快得到所需要的结果。“我的目标是使每个学生都能了解偏微分方程，当下次再
遇到这样的问题时，他们不会再担心，” Bhatia博士说，“这不需要了解太多的分析工具，总的来说，学生都反映…这个建模工具太棒了?”。
很多优秀的高科技工程公司已经看到多物理场建模可以帮助他们保持竞争力。多物理场建模工具可以让工程师进行更多的虚拟分析而不是每次都需要进行实物测试。这样，他们就可以快速而经济地优化产品。在印度尼西亚的Medrad Innovations Group中，由John Kalafut博士带领着一个研究小组，采用多物理场分析工具来研究细长的注射器中血细胞的注射过程，这是一种非牛顿流体，而且具有很高的剪切速率。
通过这项研究，Medrad的工程师制造了一个新颖的装置称为先锋型血管造影导
管（Vanguard Dx Angiographic Catheter）。同采用尖喷嘴的传统导管相比，采用扩散型喷嘴的新导管使得造影剂分布得更加均匀。造影剂就是在进行X光拍照时，将病变的器官显示得更加清楚的特殊材料。
另外一个问题就是传统导管在使用过程中可能会使得造影剂产生很大的速度，进而可能会损伤血管。先锋型血管造影导管降低了造影剂对血管产生的冲击力，将血管损伤的可能性降至最低。
关键的问题就是如何去设计导管的喷嘴形状，使其既能优化流体速度又能减少结构变形。Kalafut的研究小组利用多物理场建模方法将层流产生的力耦合到应力应变
分析中去，进而对各种不同喷嘴的形状、布局进行流固耦合分析。“我们的一个实习生针对不同的流体区域建立不同的喷嘴布局，并进行了分析，” Kalafut博士说，“我
们利用这些分析结果来评估这些新想法的可行性，进而降低实体模型制造次数”。
摩擦搅拌焊接（FSW），自从1991年被申请专利以来，已经广泛应用于铝合金的焊接。航空工业最先开始采用这些技术，现在正在研究如何利用它来降低制造成本。在摩擦搅拌焊接的过程中，一个圆柱状具有轴肩和搅拌头的刀具旋转插入两片金属的连接处。旋转的轴肩和搅拌头用来生热，但是这个热还不足以融化金属。反之，软化呈塑性的金属会形成一道坚实的屏障，会阻止氧气氧化金属和气泡的形成。粉碎，搅拌和挤压的动作可以使焊缝处的结构比原先的金属结构还要好，强度甚至可以到原来的两倍。这种焊接装置甚至可以用于不同类型的铝合金焊接。
空中客车（AirBus）资助了很多关于摩擦搅拌焊接的研究。在制造商大规模投资和重组生产线之前，Cranfield大学的Paul Colegrove博士利用多物理场分析工具帮
助他们理解了加工过程。
第一个研究成果是一个摩擦搅拌焊接的数学模型，这让空客的工程师“透视”到焊缝中来检查温度分布和微结构的变化。Colegrove博士和他的研究小组还编写了一个带有图形界面的仿真工具，这样空客的工程师可以直接提取材料的热力属性以及焊缝极限强度。
在这个摩擦搅拌焊接的模拟过程中，将三维的传热分析和二维轴对称的涡流模拟耦合起来。传热分析计算在刀具表面施加热流密度后，结构的热分布。可以提取出刀具的位移，热边界条件，以及焊接处材料的热学属性。接下来将刀具表面处的三维热分布映射到二维模型上。耦合起来的模型就可以计算在加工过程中热和流体之间的相互作用。
将基片的电磁、电阻以及传热行为耦合起来需要一个真正的多物理场分析工具。一个典型的应用是在半导体的加工和退火的工艺中，有一种利用感应加热的热壁熔炉，它用来让半导体晶圆生长，这是电子行业中的一项关键技术。
例如，金刚砂在2,000°C的高温环境下可以取代石墨接收器，接收器由功率接近10KW的射频装置加热。在如此高温下要保持炉内温度的均匀，炉腔的设计至关重要。经过多物理场分析工具的分析，发现热量主要是通过辐射的方式进行传播的。在模型内不仅可以看到晶圆表面温度的分布，还可以看到熔炉的石英管上的温度分布。
在电路设计中，影响材料选择的重要方面是材料的耐久性和使用寿命。电器小型化的趋势使得可在电路板上安装的电子元件发展迅猛。众所周知，安装在电路板上的电阻以及其他一些元件会产生大量的热，进而可能使得元件的焊脚处产生裂缝，最后导致整个电路板报废。
多物理场分析工具可以分析出整个电路板上热量的转移，结构的应力变化以及由于温度的上升导致的变形。这样做可以用来提升电路板设计的合理性以及材料选择的合理性。
计算机能力的提升使得有限元分析由单场分析到多场分析变成现实，未来的几年内，多物理场分析工具将会给学术界和工程界带来震惊。单调的“设计－校验”的设计方法将会慢慢被淘汰，虚拟造型技术将让你的思想走得更远，通过模拟仿真将会点燃创新的火花。
`
MATLAB: MATLAB是美国MathWorks公司出品的商业数学软件，用于数据分析、无线通信、深度学习、图像处理与计算机视觉、信号处理、量化金融与风险管理、机器人，控制系统等领域。 MATLAB是matrix&laboratory两个词的组合，意为矩阵工厂（矩阵实验室），软件主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式。 MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等。MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点，使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C++，JAVA的支持。 MATLAB有限元分析与应用:
《MATLAB有限元分析与应用》是2004年4月清华大学出版社出版的图书，作者是卡坦，译者是韩来彬。 内容简介： 《MATLAB有限元分析与应用》特别强调对MATLAB的交互应用，书中的每个示例都以交互的方式求解，使读者很容易就能把MATLAB用于有限分析和应用。另外，《MATLAB有限元分析与应用》还提供了大量免费资源。 《MATLAB有限元分析与应用》采用当今在工程和工程教育方面非常流行的数学软件MATLAB来进行有限元的分析和应用。《MATLAB有限元分析与应用》由简单到复杂，循序渐进地介绍了各种有限元及其分析与应用方法。书中提供了大量取自机械工程、土木工程、航空航天工程和材料科学的示例和习题，具有很高的工程应用价值。
西安市新城区某公司科研办公楼结构设计 有限元分析报告 撰写人：王平 班级：工程力学1203 学号：121010321 指导教师：张卫喜 2016年6月15日
目录 1 工程概况 (2) 2 分析依据 (3) 3 荷载与计算工况 (4) 3.1荷载简化及荷载组合 (4) 3.2 边界条件 (4) 3.3 工况 (5) 4 有限元模型 (6) 4.1 基本假定 (6) 4.2 力学模型 (6) 4.3 主要物理参数取值 (6) 4.4单元选取 (7) 4.5分网与有限元模型 (8) 5 静力分析 (9) 5.1模态结果 (9) 5.2静力分析结果 (13) 5.3 强度校核 (16) 6基于ANSYS、PKPM、手算的误差分析 (18) 6.1计算原理的不同 (18) 6.2 研究对象的复杂性 (19)
1工程概况 工程名称：西安市新城区某公司科研办公楼； 建筑所在地：西安市； 建设规模：总建筑面积约4700m2，主体结构6层，无地下室。结构总高度22.5m，底层结构高度4.5m，其余层结构高度为3.6m，几何模型图如图1所示； 抗震设防烈度：抗震设防烈度为8度，设计基本地震加速度值0.2g,第一组。场地类别为Ⅱ类，特征周期为0.35s。周期折减系数为0.75。 建筑设计使用年限：50年。 结构重要性等级：二级。 图1框架几何模型图
2分析依据 框架结构是由梁、板、柱以刚接相连接而成，构成承重体系的结构，即由梁、板、柱组成框架共同抵抗使用过程中出现的水平荷载和竖直荷载。本设计报告采用ANSYS有限元软件分析。 根据框架结构体系特点，本结构分析主要依据以下国家规范： [1]国家标准：《建筑结构荷载规范》（GB50009-2012).北京：中国建筑工业出版社.2012； [2]国家标准：《建筑抗震设计规范》（GB50011-2010).北京：中国建筑工业出版社.2010； [3]国家标准：《混凝土结构设计规范》（GB50010-2010).北京：中国建筑工业出版社.2010； [4]建筑、勘察等技术文件。
Ansys有限元分析实例[教学] 有限元分析案例:打点喷枪模组(用于手机平板电脑等电子元件粘接)，该产品主要是使用压缩空气推动模组内的顶针作高频上下往复运动，从而将高粘度的胶水从喷嘴中打出(喷嘴尺寸,0.007”)。顶针是这个产品中的核心零件，设计使用材料是:AISI 4140 最高工作频率是160HZ(一个周期中3ms开3ms关)，压缩空气压力3-8bar, 直接作用在顶针活塞面上，用Ansys仿真模拟分析零件的强度是否符合要求。 1. 零件外形设计图:
2. 简化模型特征后在Ansys14.0 中完成有限元几何模型创建:
3. 选择有限元实体单元并设定，单元类型是SOILD185，由于几何建模时使用的长度单位是mm, Ansys采用单位是长度:mm 压强: 3Mpa 密度:Ton/M。根据题目中的材料特性设置该计算模型使用的材料属性:杨氏模量 2.1E5; 泊松比:0.29; 4. 几何模型进行切割分成可以进行六面体网格划分的规则几何形状后对各个实体进行六面体网格划分，网格结果: 5. 依据使用工况条件要求对有限元单元元素施加约束和作用载荷:
说明: 约束在顶针底端球面位移全约束; 分别模拟当滑块顶断面分别以8Bar，5Bar，4Bar和3Bar时分析顶针的内应力分布，根据计算结果确定该产品允许最大工作压力范围。 6. 分析结果及讨论: 当压缩空气压力是8Bar时: 当压缩空气压力是5Bar时:
当压缩空气压力是4Bar时: 结论: 通过比较在不同压力载荷下最大内应力的变化发现，顶针工作在8Bar时最大应力达到250Mpa，考虑到零件是在160HZ高频率在做往返运动，疲劳寿命要求50百万次以上，因此采用允许其最大工作压力在5Mpa，此时内应力为156Mpa，按线性累积损伤理论[3 ]进行疲劳寿命L-N疲劳计算，进一部验证产品的设计寿命和可靠性。
Abaqus-基础与应用-第一章概述
第1章概述 有限元分析是使用有限元方法来分析静态或动态的物体或系统。在这种方法中一个物体或系统被分解为由多个相互联结的、简单、独立的点所组成的几何模型。在这种方法中这些独立的点的数量是有限的，因此被称为有限元。 1.1有限元分析简介 本节首先简要介绍有限元分析的基本概念，然后简要阐述其发展和应用概况。 1.1.1有限元分析的基本概念 在工程技术领域内，有许多问题归结为场问题的分析和求解，如位移场、应力场、应变场、流场和温度场等。这些场问题虽然已经得出应遵循的基本规律（微分方程）和相应的限制条件（边界条件），但因实际问题的复杂性而无法用解析方法求出精确解。 由于这些场问题的解是工程中迫切所需要的，人们从不同角度去寻找满足工程实际要求的近似解，有限元方法就是随着计算机技术的发展和应用而出现的一种求解数理方程的非常有效的数值方法。 有限元分析的基本思想是用离散近似的概念，把连续的整体结构离散为有限多个单元，单元构成的网格就代表了整个连续介质或结构。这种离散化的网格即为真实结构的等效计算模型，与真实结构的区别主要在于单元与单元之间除了在分割线的交点(节点)上相互连接外，再无任何连接，且这种连接要满足变形协调条件，单元间的相互作用只通过节点传递。这种离散网格结构的节点和单元数目都是有限的，所以称为有限单元法。 在单元内，假设一个函数用来近似地表示所求场问题的分布规律。这种近似函数一般用所求场问题未知分布函数在单元各节点上的值及其插值函数表示。这样就将一个连续的有无限自由度的问题，变成了离散的有限自由度的问题。根据实际问题的约束条件，解出各个节点上的未知量后，就可以用假设的近似函数确定单元内各点场问题的分布规律。 有限元方法进行结构分析主要涉及三个问题： （1）网格剖分和近似函数的选取
有限元 百科名片 有限元法（FEA，Finite Element Analysis）的基本概念是用较简单的问题代替复杂问题后 再求解。它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的工程分析手段。 目录 简介 1）物体离散化 2）单元特性分析 3）单元组集 4）求解未知节点位移 5）有限元的未来是多物理场耦合 编辑本段简介 英文：Finite Element 有限单元法是随着电子计算机的发展而迅速发展起来的一种现代计算方法。它是50年代首先在连续体力学领域--飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。 有限元法分析计算的思路和做法可归纳如下： 编辑本段1）物体离散化 将某个工程结构离散为由各种单元组成的计算模型，这一步称作单元剖分。离散后单元与单元之间利用单元的节点相互连接起来；单元节点的设置、性质、数目等应视问题的性质，描述变形形态的需要和计算进度而定（一般情况单元划分越细则描述变形情况越精确，即越接近实际变形，但计算量越大）。所以有限元中分析的结构已不是原有的物体或结构物，而是同新材料的由众多单元以一定方式连接成的离散物体。这样，用有限元分析计算所获得的结果只是近似的。如果划分单元数目非常多而又合理，则所获得的结果就与实际情况相符合。 编辑本段2）单元特性分析 A、选择位移模式
有限元分析案例 图1 钢铸件及其砂模的横截面尺寸 砂模的热物理性能如下表所示: 铸钢的热物理性能如下表所示: 一、初始条件:铸钢的温度为2875o F，砂模的温度为80o F；砂模外边界的对流边界条件：对流系数0.014Btu/hr.in2.o F，空气温度80o F；求3个小时后铸钢及砂模的温度分布。 二、菜单操作: 1.Utility Menu>File>Change Title, 输入Casting Solidification; 2.定义单元类型:Main Menu>Preprocessor>Element Type>Add/Edit/Delete, Add, Quad 4node 55； 3.定义砂模热性能:Main Menu>Preprocessor>Material Props>Isotropic,默认材料编号1, 在Density(DENS)框中输入0.054,在Thermal conductivity (KXX)框中输入0.025,在S pecific heat(C)框中输入0.28; 4.定义铸钢热性能温度表:Main Menu>Preprocessor>Material Props>-Temp Dependent->Temp Table,输入T1=0,T2=2643, T3=2750, T4=2875; 5.定义铸钢热性能:Main Menu>Preprocessor>Material Props>-Temp Dependent ->Prop Table, 选择Th Conductivity,选择KXX, 输入材料编号2,输入C1=1.44, C2=1.54, C3=1.22, C4=1.22,选择Apply,选择Enthalpy,输入C1=0, C2=128.1, C3=163.8, C4=174.2; 6.创建关键点:Main Menu>Preprocessor>-Modeling->Create>Keypoints>In Active
板结构有限元分析实例详解1：带孔平板结构静力分析本节介绍带孔平板结构静力分析问题，同时介绍布尔操作的基本用法。 8.3.1 问题描述与分析 有孔的矩形平板，左侧边缘固定，长400mm，宽200 mm，厚度为10 mm，圆孔在板的正中心，半径为40 mm，左侧全约束，右侧边缘均布应力1MPa，如图8.7所示。求板的变形、位移及应力变化情况。（材料的材料属性为：弹性模量为300000 MPa，剪切模量为0.31。） 图8.7 带孔的矩形平板 由于小孔处边缘不规则，本文采用PLANE82高阶平面单元进行分析。 8.3.2 求解过程 8.3.2.1 定义工作目录及文件名 启动ANSYS Mechanical APDL Product Launcher窗口，如图8.8所示。在License下 拉选框中选择ANSYS Multiphysics产品，在Working Directory输入栏中输入工作目 录：C:\ANSYS12.0 Structural Finite Elements Analysis and Practice\Chapter 8\8-1，在Job Name一栏中输入工作文件名：Chapter8-1。以上参数设置完毕后，单 击Run按钮运行ANSYS。
图8.8 ANSYS设置窗口菜单 可以先在目标文件位置建立工作目录，然后单击Browse按钮选择工作目录；也 可以通过单击Browse按钮选择工作文件名。 8.3.2.2 定义单元类型和材料属性 选择Main Menu>Preferences命令，出现Preferences for GUI Filtering对话框， 如图8.9所示，在Individual discipline(s) to show in the GUI中勾选Structural，过滤掉ANSYS GUI菜单中与结构分析无关的选项，单击OK按钮关闭该对话框。 图8.9 Preferences for GUI Filtering对话框
有限元分析 一个厚度为20mm的带孔矩形板受平面内张力，如下图所示。左边固定，右边受载荷p=20N/mm作用，求其变形情况 P 一个典型的ANSYS分析过程可分为以下6个步骤： ①定义参数 ②创建几何模型 ③划分网格 ④加载数据 ⑤求解 ⑥结果分析 1定义参数 1.1指定工程名和分析标题 (1)启动ANSYS软件，选择File→Change Jobname命令，弹出如图所示的[Change Jobname]对话框。 (2)在[Enter new jobname]文本框中输入“plane”，同时把[New log and error files]中的复选框选为Yes，单击确定 (3)选择File→Change Title菜单命令，弹出如图所示的[Change Title]对话框。 (4)在[Enter new title]文本框中输入“2D Plane Stress Bracket”，单击确定。 1.2定义单位
在ANSYS软件操作主界面的输入窗口中输入“/UNIT,SI” 1.3定义单元类型 (1)选择Main Menu→Preprocessor→Element Type→Add/Edit/Delete命令，弹出如图所示[Element Types]对话框。 (2)单击[Element Types]对话框中的[Add]按钮，在弹出的如下所示[Library of Element Types]对话框。 (3)选择左边文本框中的[Solid]选项，右边文本框中的[8node 82]选项，单击确定，。 (4)返回[Element Types]对话框，如下所示 (5)单击[Options]按钮，弹出如下所示[PLANE82 element type options]对话框。
SolidWorks Simulation概述 SolidWorks Simulation是一款基于有限元（即FEA数值）技术的设计分析软件，是SRAC 开发的工程分析软件产品之一。SRAC是DS SolidWorks@公司的子公司，成立于1982年，是将有限元分析带入微型电脑的先驱。1995年，SRAC开始与DS SolidWorks@公司合作开发了COSMOSWorks软件，从而进入了工程界主流有限元分析软件的市场，成为了DS SolidWorks@公司的金牌产品之一。同时，它作为嵌入式分析软件与SolidWorks无缝集成，迅速成为顶级销售产品。整合了SolidWorks CAD软件的COSMOS-Works软件在商业上取得了成功，并于2001年获得了Dassault Systemes（DS SolidWorks@母公司）的认可。2003年，SRAC公司与DS SolidWorks⑤公司合并。COSMOSWorks推出的2009版被重命名为Solid-Works Simulation。 SolidWorks是一款基于特征的参数化CAD系统软件。和许多最初在UNIX环境中开发，后来才向Windows系统开放的CAD系统不同，SolidWorks与SolidWorks Simulation茌一开始就是专为Windows操作系统开发的，所以相互整合是完全可行的。 SolidWorks Simulation有不同的程序包或应用软件以适应不同用户的需要。除了SolidWorks Simula-tionXpress程序包是SolidWorks的集成部分之外，所有的SolidWorks Simulation软件程序包都是插件式的。不同程序包的主要功能如下： ·SolidWorks SimulationXpress：能对带有简单载荷和支撑的零件进行静态分析。 ·SolidWorks Simulation：能对零件和装配体进行静力分析。 ·SolidWorks Simulation Professional：能进行零件和装配体的静态、热传导、扭曲、频率、跌落测试、优化和疲劳分析。 ·SolidWorks Simulation Premium：具有SolidWorks Simulation Professional的所有功能，外加非线性功能和动力学分析。 有限元分析概述 在数学术语中，FEA也称之为有限单元法，是一种求解关于场问题的一系列偏微分方程的数值方法。这种类型的问题涉及许多工程学科，如机械设计、声学、电磁学、岩土力学、流体动力学等。在机械工程中，有限元分析被广泛地应用在结构、振动和传热问题上。 FEA不是唯一的数值分析工具，在工程领域还有其他的数值方法，如有限差分法、边界元法和有限体积法。然而由于FEA的多功能性和高数值性能，它占据了绝大多数工程分析的软件市场，而其他方法则被归入小规模应用。使用FEA，通过不同方法理想化几何体，能够分析任何形状的模型，并得到预期的精度。当使用现代的商业软件，例如SolidWorks Simulation时，FEA理论、数值问题公式和求解方法对用户是完全透明的。 作为一个强有力的工程分析工具，FEA可以解决从简单到复杂的各种问题。一方面，设计工程师使用FEA在产品研发过程中分析设计改进，由于时间和可用的产品数据的限制，需要对所分析的模型作许多简化。另一方面，专家们使用FEA来解决一些非常深奥的问题，如车辆碰撞动力学、金属成形和生物结构分析。 不管项目多复杂或是应用领域多广，无论是结构、热传导或是声学分析，所有FEA的第一步总是相同的，都是从几何模型开始。在本课程中，即为SolidWorks的零件和装配件。我们给这些模型分配材料属性，定义载荷和约束，再使用数值近似方法，将模型离散化以便分析。 离散化过程也就是网格划分过程，即将几何体剖分成相对小且形状简单的实体，这些实体称为有限单元。单元称为“有限”的，是为了强调这样一个事实：它们不是无隈的小，而是与整个模型的尺寸相比之下适度的小。 当使用有限单元工作时，FEA求解器将把单个单元的简单解综合成对整个模型的近似解
有 限 元 分 析 作 业 作业名称 输气管道有限元建模分析 姓 名 陈腾飞 学 号 3070611062 班 级 07机制（2）班 宁波理工学院
题目描述： 输气管道的有限元建模与分析 计算分析模型如图1所示 承受内压：1.0e8 Pa R1=0.3 R2=0.5 管道材料参数：弹性模量E=200Gpa；泊松比v=0.26。 图1受均匀内压的输气管道计算分析模型（截面图） 题目分析： 由于管道沿长度方向的尺寸远远大于管道的直径，在计算过程中忽略管道的断面效应，认为在其方向上无应变产生。然后根据结构的对称性，只要分析其中1/4即可。此外，需注意分析过程中的单位统一。 操作步骤 1.定义工作文件名和工作标题 1.定义工作文件名。执行Utility Menu-File→Chang Jobname-3070611062，单击OK按钮。 2.定义工作标题。执行Utility Menu-File→Change Tile-chentengfei3070611062，单击OK按钮。 3.更改目录。执行Utility Menu-File→change the working directory –D/chen 2.定义单元类型和材料属性 1.设置计算类型 ANSYS Main Menu: Preferences →select Structural →OK
2.选择单元类型。执行ANSYS Main Menu→Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 8node 82 →apply Add/Edit/Delete →Add →select Solid Brick 8node 185 →OK Options…→select K3: Plane strain →OK→Close如图2所示，选择OK接受单元类型并关闭对话框。 图2 3.设置材料属性。执行Main Menu→Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic，在EX框中输入2e11，在PRXY框中输入0.26，如图3所示，选择OK并关闭对话框。 图3 3.创建几何模型 1. 选择ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS →依次输入四个点的坐标：input:1(0.3,0),2(0.5,0),3(0,0.5),4(0,0.3) →OK
有限元网格剖分方法概述 在采用有限元法进行结构分析时，首先必须对结构进行离散，形成有限元网格，并给出与此网格相应的各种信息，如单元信息、节点坐标、材料信息、约束信息和荷载信息等等，是一项十分复杂、艰巨的工作。如果采用人工方法离散对象和处理计算结果，势必费力、费时且极易出错，尤其当分析模型复杂时，采用人工方法甚至很难进行，这将严重影响高级有限元分析程序的推广和使用。因此，开展自动离散对象及结果的计算机可视化显示的研究是一项重要而紧迫的任务。 有限元网格生成技术发展到现在, 已经出现了大量的不同实现方法，列举如下： 映射法 映射法是一种半自动网格生成方法，根据映射函数的不同，主要可分为超限映射和等参映射。因前一种映射在几何逼近精度上比后一种高，故被广泛采用。映射法的基本思想是：在简单区域内采用某种映射函数构造简单区域的边界点和内点，并按某种规则连接结点构成网格单元。也就是根据形体边界的参数方程，利用映射函数，把参数空间内单元正方形或单元三角形(对于三维问题是单元立方体或单元四面体)的网格映射到欧氏空间，从而生成实际的网格。这种方法的主要步骤是，首先人为地把分析域分成一个个简单可映射的子域，每个子域为三角形或四边形，然后根据网格密度的需要，定义每个子域边界上的节点数，再根据这些信息，利用映射函数划分网格。 这种网格控制机理有以下几个缺点： （1）它不是完全面向几何特征的，很难完成自动化，尤其是对于3D区域。 （2）它是通过低维点来生成高维单元。例如，在2D问题中，先定义映射边界上的点数，然后形成平面单元。这对于单元的定位，尤其是对于远离映射边界的单元的定位，是十分困难的，使得对局部的控制能力下降。 （3）各映射块之间的网格密度相互影响程度很大。也就是说，改变某一映射块的网格密度，其它各映射块的网格都要做相应的调整。 其优点是：由于概念明确，方法简单，单元性能较好，对规则均一的区域，适用性很强，因此得到了较大的发展，并在一些商用软件如ANSYS等得到应用。 2 。拓扑分解法 拓扑分解法较其它方法发展较晚, 它首先是由Wordenwaber提出来的。该方法假设最后网格顶点全部由目标边界顶点组成, 那么可以用一种三角化算法将目标用尽量少的三角形完全分割覆盖。这些三角形主要是由目标的拓扑结构决定, 这样目标的复杂拓扑结构被分解成简单的三角形拓扑结构。该方法生成的网格一般相当粗糙, 必须与其它方法相结合, 通过网格加密等过程, 才能生成合适的网格。该方法后来被发展为普遍使用的目标初始三角化算法, 用来实现从实体表述到初始三角化表述的自动化转换。 单一的拓扑分解法因只依赖于几何体的拓扑结构使网格剖分不理想，有时甚至很差。 3．连接节点法 这类方法一般包括二步：区域内布点及其三角化。早期的方法通常是先在区域内布点, 然后再将它们联成三角形或四面体, 在三角化过程中, 对所生成的单元形状难于控制。随着Delaunay三角化(简称为DT ) 方法的出现, 该类方法已成为目前三大最流行的全自动网格生成方法之一。 DT法的基本原理：任意给定N个平面点Pi(i=1,2,…,N)构成的点集为S，称满足下列条件的点集Vi为Voronoi多边形。其中，Vi满足下列条件： Vi ={ X：|X- Pi|(|X- Pj|,X(R2,i(j,j=1,2,…,N }Vi为凸多边形，称{ Vi}mi=1为Dirichlet Tesselation
基于ANSYS软件的有限元分析报告 机制1205班杜星宇U201210671 一、概述 本次大作业主要利用ANSYS软件对桌子的应力和应变进行分析,计算出桌子的最大应力和应变。然后与实际情况进行比较,证明分析的正确性,从而为桌子的优化分析提供了充分的理论依据，并且通过对ANSYS软件的实际操作深刻体会有限元分析方法的基本思想，对有限元分析方法的实际应用有一个大致的认识。 二、问题分析 已知：桌子几何尺寸如图所示，单位为mm。假设桌子的四只脚同地面完全固定，桌子上存放物品，物品产生的均匀分布压力作用在桌面，压力大小等于300Pa，其中弹性模量E=9.3GPa，泊松比μ=0.35，密度ρ=560kg/m3，分析桌子的变形和应力。
将桌脚固定在地面，然后在桌面施加均匀分布的压力，可以看作对进行平面应力分析，桌脚类似于梁单元。由于所分析的结构比较规整且为实体，所以可以将单元类型设为八节点六面体单元。 操作步骤如下： 1、定义工作文件名和工作标题 （1）定义工作文件名：执行Utility Menu/ File/Change Jobname，在弹出Change Jobname 对话框修改文件名为Table。选择New log and error files复选框。 （2）定义工作标题：Utility Menu/File/ Change Title，将弹出Change Title对话框修改工作标题名为The analysis of table。 （3）点击：Plot/Replot。 2、设置计算类型 （1）点击：Main Menu/Preferences,选择Structural,点击OK。
有限元分析基础 第一章有限元法概述 在机械设计中，人们常常运用材料力学、结构力学等理论知识分析机械零构件的强度、刚度和稳定性问题。但对一些复杂的零构件，这种分析常常就必须对其受力状态和边界条件进行简化。否则力学分析将无法进行。但这种简化的处理常常导致计算结果与实际相差甚远，有时甚至失去了分析的意义。所以过去设计经验和类比占有较大比重。因为这个原因，人们也常常在设计中选择较大的安全系数。如此也就造成所设计的机械结构整体尺寸和重量偏大，而局部薄弱环节强度和刚度又不足的设计缺陷。 近年来，数值计算机在工程分析上的成功运用，产生了一门全新、高效的工程计算分析学科——有限元分析方法。该方法彻底改变了传统工程分析中的做法。使计算精度和计算领域大大改善。 §1.1 有限元方法的发展历史、现状和将来 一，历史 有限元法的起源应追溯到上世纪40年代（20世纪40年代）。1943年R.Courant从数学的角度提出了有限元法的基本观点。50年代中期在对飞机结构的分析中，诞生了结构分析的矩阵方法。1960年R.W.Clough在分析弹性力学平面问题时引入了“Finite Element Method”这一术语，从而标志着有限元法的思想在力学分析中的广泛推广。 60、70年代计算机技术的发展，极大地促进了有限元法的发展。具体表现在： 1）由弹性力学的平面问题扩展到空间、板壳问题。 2）由静力平衡问题——稳定性和动力学分析问题。 3）由弹性问题——弹塑性、粘弹性等问题。 二，现状 现在有限元分析法的应用领域已经由开始时的固体力学，扩展到流体力学、传热学和电磁力学等多个传统的领域。已经形成了一种非常成熟的数值分析计算方法。大型的商业化有限元分析软件也是层出不穷，如： SAP系列的代表SAP2000（Structure Analysis Program） 美国安世软件公司的ANSYS大型综合有限元分析软件 美国航天航空局的NASTRAN系列软件 除此以外，还有MASTER、ALGO、ABIQUES、ADINA、COSMOS等。 三，将来 有限元的发展方向最终将和CAD的发展相结合。运用“四个化”可以概括其今后的发展趋势。那就是：可视化、集成化、自动化和网络化。 §1.2 有限元法的特点 机械零构件的受力分析方法总体说来分为解析法和数值法两大类。如大家学过的材料力学、结构力学等就是经典的解析力学分析方法。在这些解析力学方法中，弹性力学的分析方法在数学理论上是最为严谨的一种分析方法。 其解题思路是：从静力、几何和物理三个方面综合考虑，建立描述弹性体的平衡、应力、应变和位移三者之间的微分方程，然后考虑边界条件，从而求出微分方程的解析解。其最大的有点就是，严密精确。缺点就是微分方程的求解困难，很多情况下，无法求解。 数值方法是一种近似的计算方法。具体又分为“有限差分法”和“有限元法”。 “有限差分法”是将得到的微分方程离散成近似的差分方程。通过对一系列离散的差分
1有限元法简介 1.1有限单法的形成 在工程技术领域内，经常会遇到两类典型的问题。其中的第一类问题，可以归结为有限个已知单元体的组合。例如，材料力学中的连续梁、建筑结构框架和桁架结构。我们把这类问题，称为离散系统。如图1-1所示平面桁架结构，是由6个承受轴向力的“杆单元”组成。尽管离散系统是可解的，但是求解图1-2所示这类复杂的离散系统，要依靠计算机技术。 图1-1 平面桁架系统
图1-2 大型编钟“中华和钟”的振动分析及优化设计（曾攀教授） 第二类问题，通常可以建立它们应遵循的基本方程，即微分方程和相应的边界条件。例如弹性力学问题，热传导问题，电磁场问题等。由于建立基本方程所研究的对象通常是无限小的单元，这类问题称为连续系统。 图1-3 V6引擎的局部 下面是热传导问题的控制方程与换热边界条件： t T c Q z T z y T y x T x ??=+??? ??????+??? ? ??????+??? ??????ρλλλ （1- 1） 初始温度场也可以是不均匀的，但各点温度值是已知的： () 00 x,y,z T T t == （1- 2） 通常的热边界有三种，第三类边界条件如下形式： ()f T-T h n T λ=??- （1- 3） 尽管我们已经建立了连续系统的基本方程，由于边界条件的限制，通常只能得到少数简单问题的精确解答。对于许多实际的工程问题，还无法给出精确的解答，例如，图1-3所示V6引擎在工作中的温度分布。这为解决这个困难，工程师们和数学家们提出了许多近似方法。 在寻找连续系统求解方法的过程中，工程师和数学家从两个不同的路线得到了相同的结果，即有限元法。有限元法的形成可以回顾到二十世纪50年代，来源于固体力学中矩阵结构法的发展和工程师对结构相似性的直觉判断。从固体力学的角度来看，桁架结构等标准离散系统与人为地分割成有限个分区后的连续系统在结构上存在相似性。 1956年M..J.Turner, R.W.Clough, H.C.Martin, L.J.Topp 在纽约举行的航空学会年会上介
《有限元分析》课程教学大纲 【课程编号】××××× 【课程名称】有限元分析/ Finite Element Analysis 【课程性质】专业核心课 【学时】144学时【实验/上机学时】144学时 【考核方式】试卷考【开课单位】XX学院 【授课对象】本科、机械设计制造及其自动化学生 一、课程的性质、目的和任务 有限元法作为边值问题的近似计算方法，随着计算机和计算技术的迅猛发展，其应用已从固体力学发展到流体力学、热力学、电磁学、声学、光学、生物学等多耦合场问题。《有限元分析基础》是材料成型类专业的一门专业基础课，主要介绍固体力学有限单元法的基本理论和应用。在对有限单元法的原理、方法进行讲授的同时配以相应的计算算例及大型工程软件的使用示例，加深学生的理解和消化。 课程教学所要达到的目的是：1、有限单元法的基本理论和实施方法；2、掌握工程结构和设备的受力及变形分析技能并最终提高他们的工程设计能力和解决实际问题的能力；3、利用ANSYS软件上机实践完成两个上机练习：刚架结构有限元分析和三维固体有限元分析；4、掌握利用有限元的加权残值法求解场问题的概念，重点介绍1维和2维热传导问。 题有限元分析。 二、教学内容、基本要求和学、课时分配 第一章：ANSYS概论（13学时） （一）基本要求：了解有限元法的分析过程，ANSYS 15.0的安装与启动，前处理、加载并求解、后处理。 （二）教学内容和课时分配： 1、有限元法的分析过程，ANSYS 15.0的安装与启动（2学时）
2、系统要求、设置运行参数（1学时） 3、ANSYS分析的基本过程（1学时） 4、实验内容（9学时） 实验1 梁的有限元建模与变形分析（1学时） 实验目的和要求： 1)要求选择不同形状的截面分别进行计算； 2) 梁截面分别采用以下三种截面； 3) 设置计算类型； 重点：有限元法的分析过程，ANSYS 15.0的安装与启动； 难点：ANSYS分析的基本过程； 第二章：图形用户界面（13学时） （一）基本要求：了解ANSYS软件界面下各窗口的功能，具体包括应用命令菜单、主菜单、工具栏、输入窗口、图形窗口和输出窗口。ANSYS架构及命令，具体包括简单模型的建立、材料属性输入、单元的选择和划分、求解处理和后置处理。 （二）教学内容和课时分配： 1、ANSYS 15.0图形用户界面的组成（1学时） 具体包括应用命令菜单、主菜单、工具栏、输入窗口、图形窗口和输出窗口。ANSYS 架构及命令，具体包括简单模型的建立、材料属性输入、单元的选择和划分、求解处理和后置处理2、对话框及其组件、通用菜单，输入窗口 2、主菜单，输出窗口，图形窗口的功能（1学时） 3、个性化界面（1学时） 4、实验内容（10学时） 实验1 超静定桁架的有限元建模与分析 实验目的和要求：上机熟悉ANSYS软件的命令，并对简单的例题进行有限元静、动态分析。 重点（黑体，小四号字）：ANSYS 15.0图形用户界面的组成；
有限元法的基本思想及计算步骤 有限元法是把要分析的连续体假想地分割成有限个单元所组成的组合体，简称离散化。这些单元仅在顶角处相互联接，称这些联接点为结点。离散化的组合体与真实弹性体的区别在于：组合体中单元与单元之间的联接除了结点之外再无任何关联。但是这种联接要满足变形协调条件，即不能出现裂缝，也不允许发生重叠。显然，单元之间只能通过结点来传递内力。通过结点来传递的内力称为结点力，作用在结点上的荷载称为结点荷载。当连续体受到外力作用发生变形时，组成它的各个单元也将发生变形，因而各个结点要产生不同程度的位移，这种位移称为结点位移。在有限元中，常以结点位移作为基本未知量。并对每个单元根据分块近似的思想，假设一个简单的函数近似地表示单元内位移的分布规律，再利用力学理论中的变分原理或其他方法，建立结点力与位移之间的力学特性关系，得到一组以结点位移为未知量的代数方程，从而求解结点的位移分量。然后利用插值函数确定单元集合体上的场函数。显然，如果单元满足问题的收敛性要求，那么随着缩小单元的尺寸，增加求解区域内单元的数目，解的近似程度将不断改进，近似解最终将收敛于精确解。 用有限元法求解问题的计算步骤比较繁多，其中最主要的计算步骤为： 1）连续体离散化。首先，应根据连续体的形状选择最能完满地描述连续体形状的单元。常见的单元有：杆单元，梁单元，三角形单元，矩形单元，四边形单元，曲边四边形单元，四面体单元，六面体单元以及曲面六面体单元等等。其次，进行单元划分，单元划分完毕后，要将全部单元和结点按一定顺序编号，每个单元所受的荷载均按静力等效原理移植到结点上，并在位移受约束的结点上根据实际情况设置约束条件。 2）单元分析。所谓单元分析，就是建立各个单元的结点位移和结点力之间的关系式。现以三角形单元为例说明单元分析的过程。如图1所示，三角形有三个结点i，j，m。在平面问题中每个结点有两个位移分量u，v和两个结点力分量F x，F y。三个结点共六个结点位移分量可用列阵(δ)e表示： {δ}e=[u i v i u j v j u m v m]T 同样，可把作用于结点处的六个结点力用列阵{F}e表示： {F}e=[F ix F iy F jx F jy F mx F my]T 应用弹性力学理论和虚功原理可得出结点位移与结点力之间的关系
有限元法的概述 有限元方法(Finite Element Method)是力学，数学物理学，计算方法，计算机技术等多种学科综合发展和结合的产物。在人类研究自然界的三大科学研究方法(理论分析，科学试验，科学计算)中，对于大多数新型领域，由于科学理论和科学实践的局限性，科学计算成为一种最重要的研究手段。在大多数工程研究领域，有限元方法是进行科学计算的重要方法之一;利用有限元方法几乎可以对任意复杂的工程结构进行分析，获取结构的各种机械性能信息，对工程结构进行评判，对工程事故进行分析。有限元法在设计过程中有极为关键的作用。 人们对各种力学问题进行分析求解，其方法归结起来可以分为解析法(Analytical Method)和数值法(Numeric Method).如果给定一个问题，通过一定的推导可以用具体的表达式来获得问题的解答，这样的求解方法就称为解析法。但是由于实际结构物的复杂性，除了少数极其简单的问题外，绝大多数科学研究和工程计算问题用解析法求解式极其困难的。因此，数值法求解便成为了一种不可替代的广泛应用的方法，并取得了不断的发展，如有限元法，有限差分法，边界元方法等都是属于数值求解方法。其中有限元法式 20 世纪中期伴随着计算机技术的发展而迅速发展起来的一种数值分析方法，它的数学逻辑严谨，物理概念清晰，应用非常广泛，能活灵活现处理和求解各种复杂的问题。有限元方法采用矩阵式来表达基本公式，便于计算机编程，这些优点赋予了它强大的生命力。 有限元方法的实质是将复杂的连续体划分成为有限多个简单的单元体，化无限自由度问题为优先自由度问题，将连续场函数的(偏)微分方程的求解问题转化为有限个参数的代数方程组的求解问题。用有限元方法分析工程结构的问题时，将一个理想体离散化后，如何保证其数值的收敛性和稳定性是有限元理论讨论的主要内容之一，而
1-1有限单元法的分析过程 土木工程、岩土工程等学科中的弹塑性、粘弹性、粘塑性力学，水利、码头工程等的流体力学和流固耦合作用，交通工程等学科中的层状介质路面力学等都是力学学科的重要分支，其研究结果最终归结为求解数学物理方程边值或初值问题。但这些学科传统的研究成果只对较为简单、规则的问题才能获得解析解答，大量的实际科学、工程计算问题，由于数学上的困难无法得到解决。 从有限单元法命名至今已经历了40年左右的发展，用有限单元法来解决问题，从理论上讲，无论是简单的一维杆系结构，还是受复杂荷载和不规则边界情况的二维平面、轴对称问题、三维空间块体等问题的静力、动力和稳定性分析，考虑材料具有非线性力学行为和有限变形的分析，温度场、电磁场，流体、液-固体、结构与土壤相互作用等工程复杂问题的分析都可得到满意的解决，而且其基本思路和分析过程是基本相同的。作为本科教材，本书将只介绍应用最广泛的、以结点位移作为基本未知量的“位移型有限元（简称位移元）”，讨论范围仅限于最基本的线弹性问题。 一、结构离散化 应用有限元法来分析工作问题的第一步，首先是将结构进行离散化。其过程就是将待分析的结构（或更数学化一点也可称为求解域）用一些假想的线或面进行切割，使其成为具有选定切割形状的有限个单元体（注意单元体跟材料力学中的微元体是根本不同的，它的尺度是有限值而不是微量）。这些单元体被认为仅仅在单元的一些指定点
处相互连接，这些单元上的点则称为单元的结点。这一步的实质也就是用单元的集合体来代替原来待分析的结构。 为了便于理论推导和用计算程序进行分析，一般来说结构离散化的具体步骤是：建立单元和整体坐标系、对单元和结点进行合理的编号，为后续有限元分析准备出所必需的数据化信息。目前市面上有各种类型的有限元分析软件，一般都具有友好的用户图形接口和图形直观输入、输出计算信息的强大功能，使用者应用这些软件越来越方便。即便如此，使用这些大型软件的第一步仍需“建模”工作，即建立离散化模型的准备所需的数据。 二、确定单元位移模式 结构离散化后，接下来的工作就是对结构离散化所得的任一典型单元进行所谓单元特性分析。为此，首先必须对该单元中任意一点的位移分布做出假设，即在单元内用只具有有限自由度的简单位移代替真实位移。对位移元来说，就是将单元中任意一点的位移近似地表示成该单元结点位移的函数，该位移称为单元的位移模式或位移函数。位移函数的假设合理与否，将直接影响到有限元分析的计算精度、效率和可靠性。目前比较常用的方法是以多项式作为位移模式，这主要是因为多项式的微积分运算比较简单，而且从泰勒级数展开的意义来说，任何光滑函数都可以用无限项的泰勒级数多项式来展开。位移模式的合理选择，是有限单元法的最重要内容之一，所谓创建一种新型的单元，位移模式的确定是其核心内容。本书后续各章将结合具体的单元进行较详细的讨论。
有限元综述 蔡璟、吕丹丹、李川 摘要：有限元法（Finite Element Method）是一种高效能、常用的数值计算方法。1965年“有限元”这个名词第一次出现，经历了三十多年的发展历史，理论和算法都已经日趋完善。如今，有限元在工程上得到广泛应用。本文首先介绍了有限元的研究背景和意义，其次从它的诞生、主要特点以及解题步骤三方面阐述相关概念，再讨论传统有限元算法及优化算法、有限元与其他算法结合得到的混合算法两个方面来分类阐述各自的研究现状与特点，最后总结有限元算法的应用以及发展趋势。 关键词：有限元法，FEM，经典算法，优化算法，网格优化，Herrmann算法，时域有限元，混合算法，矩量法，时域有限差分，应用研究，边界元法，光滑粒子法，发展趋势
前言 有限元法（Finite Element Method）是一种高效能、常用的数值计算方法，其基本思想是由解给定的泊松方程化为求解泛函的极值问题。有限元法在早期是以变分原理为基础发展起来的，所以它广泛地应用于以拉普拉斯方程和泊松方程所描述的各类物理场中（这类场与泛函的极值问题有着紧密的联系）。自从1969年以来，某些学者在流体力学中应用加权余数法中的迦辽金法(Galerkin)或最小二乘法等同样获得了有限元方程，解决了物理场应用中的限制。经历几十年的发展，有限元法已经被广泛用于各个领域。 1.研究背景和意义 有限元法的思想首先由 R. Courant 在 1943 年提出，十九世纪六十年代数值分析科学家认识了有限元基本思想，建立了有限元方法的数学基础。其中，我国数学家冯康独立地提出了有限元方法，将其命名为“基于变分原理的差分格式”，对有限元方法的创始及奠基工作做出了重要贡献。 以变分原理为基础建立起来的有限元法,因其理论依据的普遍性,不仅广泛地被应用于各种结构工程,而且作为一种声誉很高的数值分析方法已被普遍推广并成功地用来解决其他工程领域中的问题,例如热传导!渗流!流体力学、空气动力学、土壤力学、机械零件强度分析、电磁场工程问题等等。 有限元法由于可以模拟任意几何模型和各种特性的复杂材料而且具有的适应性强、程序较为通用等优势而得到了长足的发展。同时，结合其他方法和理论呈现出广阔的应用前景，如自适应网格剖分、三维场建模求解、耦合问题、开放域问题等领域取得较多成果。现阶段，为了进一步拓宽求解问题的广泛性以及适应求解问题对高精度，高复杂程度的要求，有限元还需要进行突破性的工作。2.有限元研究概况 2.1有限元的诞生 1943年，数学家库朗德第一次提出了可在定义域内分片地使用展开函数来表达其上的未知函数，最早提出有限元法基本思想。20世纪50年代，飞机设计师们发现无法用传统的力学方法分析飞机的应力、应变等问题。波音公司的一个技术小组，首先将连续体的机翼离散为三角形板块的集合来进行应力分析，经过一番波折后获得前述的两个离散的成功。20世纪50年代，大型电子计算机投入了解算大型代数方程组的工作，这为实现有限元技术准备好了物质条件。1960年前后，美国的R.W.Clough教授及我国的冯康教授分别独立地在论文中提出了“有限单元”这样的名词。此后，这样的叫法被大家接受，有限元技术从此正式