有限单元法原理与应用
朱伯芳《有限单元法原理与应用》怎么样
朱伯芳《有限单元法原理与应用》怎么样《有限单元法原理与应用》是朱伯芳编写的一本关于有限单元法（Finite Element Method，简称FEM）的专业教材。
该书主要介绍了有限单元法在工程领域中的应用原理和实践技巧。
本文将从书籍的结构、内容、作者的资历和读者反馈等方面进行详细的分析，以期给出一个全面的评价。
首先，让我们来看看书籍的结构和内容。
《有限单元法原理与应用》共分为六章，涵盖了有限单元法的基础知识、离散法、结构分析、流体力学、固体力学和热传导，主要包括一维和二维有限元的建模方法、应力分析和位移分析等内容。
每章都以一个引言开始，然后逐步介绍理论和应用的细节，最后以一些实例或案例进行总结。
这种递进的结构使读者可以循序渐进地学习和理解有限单元法的原理和应用，使其成为一本适合初学者和专业人士的教材。
其次，让我们来看看作者朱伯芳的资历。
朱伯芳是国内有限单元法领域的知名专家，他在该领域拥有多年的研究和教学经验。
他曾在国内外多所知名大学进行过学术交流，并发表了大量有关有限单元法的论文。
他的丰富经验和深厚的理论基础使得他编写的教材具有一定的权威性和可信度。
另外，让我们来看看读者的反馈。
据我了解，该书在工程领域和学术界中广泛应用，并受到了读者的好评。
读者认为这本书的内容详实、通俗易懂，特别是在实际工程应用方面提供了很多实用的技巧和方法。
此外，该书还提供了大量的习题和案例分析，有助于读者巩固学到的知识并提升解决实际问题的能力。
然而，也有读者认为该书在某些地方没有深入或过于简单，希望在后续版本中能够进行改进和完善。
总的来说，《有限单元法原理与应用》是一本涵盖了有限单元法基本原理和应用技巧的优秀教材。
它的结构清晰，内容详实，适合初学者和专业人士学习和参考。
作者的资历和读者的反馈也增加了该书的可信度。
然而，该书也存在一些可以改进的地方，希望在后续版本中能够进行修改和完善。
总体来说，这本书可以说是有限单元法领域的一本经典之作，值得推荐给广大读者。
有限单元法基本思想,原理,数值计算过程
有限单元法学习报告在对力学问题分析求解过程中，方法可以概括为两种方法，一种为解析法，对具体问题具体分析，通过一定的推导用具体的表达式获得解答，由于实际工程中结构物的复杂性，此方法在处理工程问题是十分困难的；另一种是数值法，有限元法是其中一种方法，其数学逻辑严谨，物理概念清晰，又采用矩阵形式表达基本公式，便于计算机编程，因此在工程问题中获得广泛的应用。
有限元法基本原理是，将复杂的连续体划分为简单的单元体；将无限自由度问题化为有限自由度问题，因为单元体个数是有限的；将偏微分方程求解问题化为有限个代数方程组的求解问题。
通常以位移为基本未知量，通过虚功原理和最小势能原理来求解。
基本思想是先化整为零，即离散化整体结构，把整体结构看作是由若干个通过结点相连的单元体组成的整体；再积零为整，通过结点的平衡来建立代数方程组，最后计算出结果。
我将采用最简单的三结点三角形为基本单元体，解决弹性力学中的平面问题为例，解释有限单元法的基本原理、演示数值计算过程和一般性应用结论。
一、离散化解决平面问题时，主要单元类型包括三角形单元（三结点、六结点）和四边形单元（四结点矩形、四结点四边形、八结点四边形）等。
选用不同的单元会有不同的精度，划分的单元数越多，精度越高，但计算量也会越大。
因此在边界曲折，应力集中处单元的尺寸要小些，但最大与最小单元的尺寸倍数不宜过大。
在集中力作用点及分布力突变的点宜选为结点，不同厚度，不同材料不能划分在同一单元中。
三角形单元以内角接近60°为最好。
充分利用对称性与反对称性。
二、单元分析将一个单元上的所有未知量用结点位移表示，并将分布在单元上的外力等效到结点上。
1、位移函数选取：根据有限元法的基本思路，将连续体离散为有限的单元集合后，此时单元体满足连续性、均匀性、各向同性、完全线弹性假设。
单元与单元之间通过结点连接并传递力，位移法（应用最广）以结点位移δi=(u i v i)T为基本未知量，以离散位移场代替连续位移场。
有限单元法及工程应用
有限单元法及工程应用有限单元法（Finite Element Method，FEM）是一种数值计算方法，广泛应用于工程领域。
它是一种将复杂的连续体分割为有限个简单形状的小单元，并将偏微分方程转化为代数方程求解的方法。
有限单元法通过将计算领域离散化为一个有限的单元网络，然后通过求解每个单元上的方程来得到整个计算领域的解。
这种方法在解决复杂问题上具有很大的优势，并已经在工程应用中得到广泛应用。
有限单元法在工程应用中有许多不同的方面。
以下是其中一些主要的应用领域：1. 结构力学分析：有限单元法可以用于结构的形状、变形、应力和振动等问题的分析。
通过将结构离散为有限个单元，可以准确地计算结构的应力分布和变形情况，进而评估结构的稳定性和可靠性。
这在建筑、桥梁、飞机和船舶等领域中得到广泛应用。
2. 热传导分析：有限单元法可以用于热传导问题的分析，如温度分布、热流量和热应力等。
通过建立传导方程和边界条件，可以计算不同材料和结构的热行为，进而为热处理、热设备设计和热工艺优化提供指导。
3. 流体力学分析：有限单元法可以用于求解流体力学方程，如流体流动、湍流、传质和热传递等。
通过将流体域划分为有限个单元，可以计算流速、压力和流体力学特征等。
这在空气动力学、水力学和化工工艺等领域中得到广泛应用。
4. 电磁场分析：有限单元法可以用于求解静电场、磁场和电磁波等问题。
通过建立电磁方程和边界条件，可以计算电场、磁场和电磁波的分布和特性。
这在电力系统、电子器件和电磁辐射等领域中得到广泛应用。
5. 生物医学工程：有限单元法可以应用于生物医学领域的各种问题，如骨骼力学、组织力学、生物电流和生物传递等。
通过对生物体或医学设备建立有限元模型，可以模拟和预测生物体的行为和反应，为生物医学研究和医学工程设计提供指导。
以上只是有限单元法在工程应用中的一部分方面。
由于其灵活性和适用性，有限单元法被广泛应用于各种工程领域，为工程师提供了一种有效的工具来解决现实世界中的复杂问题。
有限单元法原理与应用
有限单元法原理与应用有限单元法（Finite Element Method，简称FEM）是一种数值计算方法，广泛应用于工程领域的结构分析、流体力学、热传导等问题的求解。
它将复杂的结构或物理现象分割成有限数量的简单单元，通过对每个单元进行数学建模和分析，最终得出整个系统的行为。
本文将介绍有限单元法的基本原理和其在工程领域中的应用。
有限单元法的基本原理是将连续的物理现象离散化为有限数量的单元，每个单元都可以通过简单的数学方程来描述。
这些单元相互连接，形成一个整体的系统，通过对每个单元的行为进行分析，最终得出整个系统的行为。
有限单元法的核心思想是将复杂的问题简化为简单的数学模型，通过数值计算方法求解这些模型，从而得到系统的行为。
有限单元法在工程领域有着广泛的应用。
在结构分析中，可以用有限单元法来模拟各种复杂的结构，如桥梁、建筑、飞机机翼等，通过对结构的受力、变形等进行分析，来评估结构的安全性和稳定性。
在流体力学中，有限单元法可以用来模拟流体的流动行为，如水流、气流等，通过对流体的速度、压力等进行分析，来优化流体系统的设计。
在热传导问题中，有限单元法可以用来模拟物体的温度分布和传热行为，如热传导、对流、辐射等，通过对热场的分析，来优化热传导系统的设计。
有限单元法的应用还不仅限于工程领域，它也被广泛应用于地质勘探、医学图像处理、材料科学等领域。
在地质勘探中，有限单元法可以用来模拟地下岩层的力学行为，来评估地下资源的分布和开采方案。
在医学图像处理中，有限单元法可以用来模拟人体组织的力学行为，来辅助医学诊断和手术设计。
在材料科学中，有限单元法可以用来模拟材料的力学性能和热物理性能，来指导新材料的设计和制备。
总的来说，有限单元法作为一种数值计算方法，具有广泛的应用前景和重要的理论意义。
通过对有限单元法的深入理解和应用，可以更好地解决工程领域中的复杂问题，推动工程技术的发展和进步。
希望本文对有限单元法的原理和应用有所帮助，也希望读者能够进一步深入研究和应用有限单元法，为工程领域的发展做出更大的贡献。
有限单元法原理与应用
yi yj ym
ai x j ym xm y j , bi y j ym , ci xm x j a j x m y i xi y m , b j y m y i , c j xi x m am xi y j x j yi , bm yi y j , cm x j xi
v 4 5 x 6 y
因此可以得到:
ui 1 2X i 3Yi ui 1 2X j 3Yj ui 1 2X m 3Y m v i 4 5X i 6Yi v i 4 5X j 6Yj v i 4 5X m 6Y m
2.3单元应变
•单元内的应变分量可用矩阵表示为：
u x x y v y xy u v y x
应变分 量是常 量
其子矩阵：
bi 1 Bi 0 2A c i
在（x，y）中，
， D , ,
,
,
D
,
,
T
,
D
,
T
整体坐标 系的弹性 矩阵
0 ci bi
bj 0 cj
0 cj bj
bm 0 cm
0 e e c m B bm
应力矩阵
S D 0
e
S D B S i
Sj
Sm
有限单元法原理及应用
目
录
(4) 形成整体系统的矩阵方程
(5) 约束处理，求解系统方程
(6) 其它参数计算
整理课件
6
第一章 概述
图1-2 工程问题有限单元法分析流程
整理课件
7
第一章 概述
1.3 工程实例
返 回 章 节 目 录
(a) 铲运机举升工况测试
(b) 铲运机工作装置插入工况有限元分析
图1-3 WJD-1.5型电动铲运机
结构是用来承受和传递载荷的。如果不计材料的
返 回
应变，在其受到任意载荷作用时其形状和位置没有发
章
生刚体位移时，称之为几何不变结构或几何稳定结构，
节 目
反之则称为几何可变结构或几何不稳定结构。几何可
录
变结构不能承受和传递载荷。对结构进行几何构造分
析也是能够对工程结构作有限单元法分析的必要条件。
整理课件
对称结构在正对称载荷下，对称轴截面上只能产生 正对称的位移，反对称的位移为零；对称结构在反对称 载荷下，对称轴截面上只有反对称的位移，正对称的位 移为零。 (1) 具有奇数跨的刚架
① 正对称载荷作用
(a) 对称刚架
（b) 变形状态分析
(c) 对称性利用
图2-22对称性利用示意图
整理课件
19
第二章 结构几何构造分析
整理课件
16
第二章 结构几何构造分析
(3) 按结构自由度分 ①静定结构——自由度为零的几何不变结构。其特征: a. 静定结构的内力及支座反力可全部由平衡方程式
求出，并且解答是唯一的。
b. 静定结构的内力及支座反力与材料的性质和截面 特征构上无外载荷作用时，其内力及支座反力 全为零。
感谢
大学力学论坛
有限单元法原理及应用
有限单元法原理及应用有限单元法（Finite Element Method，FEM）是一种数值分析方法，广泛应用于工程领域中结构力学、流体力学、热传导等问题的数值求解。
它的基本思想是将一个复杂的结构或物理现象分割成有限数量的简单单元，通过对单元的力学行为进行建模，最终得到整个系统的数值解。
本文将围绕有限单元法的原理及其在工程领域中的应用进行详细介绍。
有限单元法的原理。
有限单元法的原理基于力学原理和数学方法，其基本步骤包括，建立数学模型、离散化、单元划分、建立单元刚度矩阵和载荷向量、组装和求解方程、计算结果后处理等。
在建立数学模型时，需要根据实际问题选择合适的数学方程和边界条件，将问题转化为求解一组代数方程。
离散化是指将连续的物理问题划分成若干个小单元，每个单元内的物理行为可以用简单的数学方程描述。
单元划分是将整个结构或领域划分成若干个有限单元，通常采用三角形、四边形、四面体、六面体等几何形状。
建立单元刚度矩阵和载荷向量是对每个单元进行力学行为的建模，根据材料性质和几何形状计算单元的刚度矩阵和载荷向量。
组装和求解方程是将所有单元的刚度矩阵和载荷向量组装成整个系统的刚度矩阵和载荷向量，然后通过数值方法求解代数方程组。
最后，计算结果后处理是对数值解进行分析和可视化，评估结构的性能和稳定性。
有限单元法的应用。
有限单元法在工程领域中有着广泛的应用，包括结构力学、流体力学、热传导等方面。
在结构力学中，有限单元法可以用于分析和设计各种结构，如桥梁、建筑、机械零件等。
通过对结构的受力分析，可以评估结构的安全性和稳定性，指导工程设计和施工。
在流体力学中，有限单元法可以用于模拟流体的流动行为，如水力学、空气动力学等问题的数值模拟。
在热传导中，有限单元法可以用于分析材料的热传导性能，评估材料的热稳定性和散热效果。
总结。
有限单元法作为一种数值分析方法，在工程领域中有着重要的应用价值。
通过对结构、流体、热传导等问题的数值模拟，可以为工程设计和科学研究提供重要的参考和支持。
有限元法的基本原理
第二章有限单元法的基本原理作为一种比较成熟的数值计算方法，有限元的数学基础是变分原理。
经过半个过世纪的发展，它的数学基础已经比较完善。
从数学角度分析，有限元法是以变分原理和剖分插值为基础的数值计算方法。
它广泛的应用于解算各种类型的偏微分方程，特别对椭圆型方程，因为椭圆型方程的边值问题等价于适当的变分问题，即能量积分的级值问题。
通过变分，导出相应的泛涵，再把作用域从几何上剖分为足够小的单元，这样就能够用简单的图形去拟合复杂的边界，用简单的初等函数去模拟单元的性质。
在解算中先对每个单元进行分析，后在通过连接单元的节点对作用域的整体进行分析，就是对泛涵求极值，从而把一个复杂的偏微分方程求解问题，变成解线形代数方程组的问题。
尽管这样会出现大量的未知数，由于采用了矩阵分析的方法，总体上很有规律，适合编制程序用计算机完成。
通常的数学考虑包括这些：1）从古典变分方法原理去定义微分方程边值问题的广义解以及在古典变分方法的框架对有限元进行理论分析。
2）保证偏微分方程边值问题的提法正确，即要求解存在、唯一和稳定，即保证数值解法是可靠的。
3）有限元中重要的一点是采用了分块多项式插值函数，因此，有限元的误差估计转化为插值逼近的误差估计问题。
4）有限元的收敛性和误差估计。
由于本文是应用有限元的理论解决大地测量中的问题，因此，这里将不讨论上叙问题，而是从固体力学的基本方程出发，通过虚功原理建立起离散化的有限元方程。
另外，还以八节点六面体单元为例，简要叙述了实际中最常用的等参单元的概念及其数值变化的一些公式。
§2.1 弹性力学基本方程有限元法中经常要用到弹性力学的基本方程，这里写出这些方程的矩阵表达式。
2-1-1、平衡方程对任意一点的受力情况分析，沿坐标轴方向x, y ,z分解得到平衡方程0*00000000=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂z y xxz yz xy z y x F F F z yzz x y z y x τττσσσ 记为： 0=+F A σ其中A 是微分算子，F 是体积力向量。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
E2
2 2
n(1 n 2 ) n (1 ) 0 2 2 1 n 2(1 1 ) 1 12 0 0 0 m(1 1 ) 1 1 2n 22
当层面与x轴有一倾角时[D]必须进行变换。
则在（x’,y’）中，
Hale Waihona Puke 对于平面应力问题，温度引起的初应变为
0
T 其中：α为膨胀系 T 数 0 温度不引起剪切变形
对于平面应变问题
μ为泊松 比
0
T 1 T 0
对于层状各项异性材料： 由于线膨胀系数可能随方向的变化而变化，设X’,Y’为层状材料 的主应力方向，设相应的膨胀系数为α1，α2， X’与X夹角为φ，则在 平面应力问题的初应变为：
2.3单元应变
•单元内的应变分量可用矩阵表示为：
u x x y v y xy u v y x
应变分 量是常 量
其子矩阵：
bi 1 Bi 0 2A c i
3.各向异性体
设y轴垂直于层面，则有如下应力应变关 系：
x x E 1 2 y E 2 1 z E 1 , xy xy G 2 y 2 x E 2 y E 2 2 z E 2 , yz yz G 2 z 1 x E 1 2 y E 2 z E 1 , zx 2 1 ) zx E 1 (1
E
E
由矩阵表示为：
x 1 E y 2 1 xy 0
1 0
0 0 1 2
x x 0 y y 0 xy xy 0
解得
u
1 (ai bi x ci y)ui (a j b j x c j y)u j (am bm x cm y)u m 2A
v
1 (ai bi x ci y)vi (a j b j x c j y)v j (am bm x cm y)vm 2A
1 xi 1 A 1 x j 2 1 x m
yi yj ym
ai x j ym xm y j , bi y j ym , ci xm x j a j x m y i xi y m , b j y m y i , c j xi x m am xi y j x j yi , bm yi y j , cm x j xi
0 ci (i, j , m) bi
带入到位 移函数
bi 1 0 2A c i
0 ci bi
bj 0 cj
0 cj bj
bm 0 cm
ui v b 0 i u j 1 i cm 0 v j 2A c i bm u m v m
记
E1 E 2 n,G2 E 2 m
平面应力问题的弹性矩阵为：
n n 2 E2 D 1 (1 n 2 ) n 2 2 0 0
平面应变问题的弹性矩阵：
0 2 m(1 n 2 ) 0
D
(1 1 ) 1 1 2n
由以上可以解出：
x y xy
E
1
2
( x x 0 y y 0 ) ( x x 0 y y 0 ) 1 xy 0 ) ( xy xy 0） 2 2 1
E
1
2
( xy 2 ) (1
0
T 1 T 0
1 (1 ) E(1 ) D 1 (1 ) 1 1 2 0 0
0 1 2 2(1 ) 0
在（x，y）中，
， D , ,
,
,
D
,
,
T
,
D
,
T
整体坐标 系的弹性 矩阵
因此，单元位移：
u r N e INi v
IN j
INm
e
其中：
1 0 I 0 1
位移模式需满足以下三个条件： 1。位移模式必须反映单元的刚体位移 2。位移模式必须反映单元的常量应变 3。位移模式应尽可能反映位移的连续性
0
T , 1 x0 y ,0 2T 0 x ,y ,0
, 0
T
0
cos2 sin2 sin cos
为了简化位移函数的表达式，记：
ai bi x ci y N a j b j x c j y Ni j 2A 2A a m bm x c m y Nm 2A u N i ui N j u j N m u m
简单表达式：
v Ni vi N j v j N m vm
0 ci bi
bj 0 cj
0 cj bj
bm 0 cm
0 e e c m B bm
2.4初应变
初应变是指与应力无关，由温度变化、收缩、晶体生长等因素引起的应 变，即
0
x 0 yo xy 0
E S i D Bi 2 2 )A (1 1
bi bi
ci 2
ci ci 1 bi 2
2.各项同性体——平面应变
1 x ( x E 1 y ( y E (1 2 xy E 其中， x
z x ) T ) xy 1 ( z x y ) T 0 E y z ) T
再次得到
D ( 0 )
其中：
或
D ( 0 )
[D]为弹性矩阵，它决定于E和μ
bi 1 0 2A c i
0 ci bi
bj 0 cj
0 cj bj
bm 0 cm
ui v b 0 i u j 1 i cm 0 v j 2A c i bm u m v m
0 ci bi
bj 0 cj
0 cj bj
bm 0 cm
0 e e c m B bm
应力矩阵
S D 0
e
S D B S i
Sj
Sm
v 4 5 x 6 y
因此可以得到:
ui 1 2X i 3Yi ui 1 2X j 3Yj ui 1 2X m 3Y m v i 4 5X i 6Yi v i 4 5X j 6Yj v i 4 5X m 6Y m
D
,
D, D
2 sin cos 2 cos 2 sin cos sin cos cos2 sin 2 sin2
经过变换，在（x，y）中，γxy0可能不为0
2.5平面应力
1.各向同性体——平面应力
有广义胡克定律：
y x x x 0 E E y x y y 0 E E 2 ) (1 xy xy xy 0 E
2.2位移函数
每一结点具有两个位移分量：
i
ui v i
每个单元六个结点的位移分量可以表示为一个向量：
i e j m 假定单元内的位移分量是坐标的线性函数：
u 1 2 x 3 y
有限单元法原理与应用
——平面问题
2.1连续介质的离散化
连续介质的离散
连续介质的有限单元分析包含三个基本方面：介质的离散化、单元特性计算以及单元 组合体的结构分析。
对于二维连续介质，以图所示的建筑在岩石 基础上的支墩坝为例，用有限单元法进行分 析的步骤如下： (1)用虚拟的直线把原介质分割成有限个 三角形单元，这些直线是单元的边界，几条 直线的交点称为结点。 (2)假定各单元在结点上互相铰接，结点 位移是基本的未知量。 (3)选择位移函数。 (4)通过位移函数，用结点位移唯一地表 示单元内任一点的应变；再利用广义虎克定 律，用结点位移可唯一地表示单元内任一点 的应力。 (5)利用能量原理，找到与单元内部应力状态等效的结点力，再利用单元应力与结点位移 的关系，建立等效结点力与结点位移的关系。 (6)将每一单元所承受的荷载，按静力等效原则移置到结点上。 (7)在每一结点建立用结点位移表示的静力平衡方程，得到一个线性方程组：解出这个 方程组，求出结点位移，然后可求得每个单元的应力。