有限元方法讲义
有限元总结讲义
5 网格分界面和分界点 结构中的一些特殊界面和特殊点应分为网格边界
或节点，以便定义材料特性、物理特性、载荷和位移 约束条件。
常见的特殊界面和特殊点有材料分界面、 几何尺寸突变面、分布载荷分界线(点)、集中 载荷作用点和位移约束作用点等。
6 位移协调性 位移协调是指单元上的力和力矩能够通过节点传
递相邻单元。 为保证位移协调，一个单元的节点必须同时也是
划分网格原则
划分网格是建立有限元模型的一个重要环节，它要求 考虑的问题较多，需要的工作量较大，所划分的网格 形式对计算精度和计算规模将产生直接影响。为建立 正确、合理的有限元模型，这里介绍划分网格时应考
虑的一些基本原则（影响因素）。
1 网格数量 5 网格分界面和分界点 2 网格疏密 6 位移协调性 3 单元阶次 7 网格布局 4 网格质量 8 节点和单元编号
体内部趋近于边界的应力分量的关系。
Fsj ijni
位移边界条件 就是弹性体表面的变形协调，弹性体临近表
面的位移与已知边界位移相等
面（应）边界条件
给定面力分量 X ,Y ,Z 边界 —— 应力边界
cos( N ,x ) l cos( N , y ) m cos( N ,z ) n
N
Z
Y X
2 网格疏密 网格疏密是指在结构不同部位采用大小不
同的网格，以适应计算数据的分布特征。 在计算数据变化梯度较大的部位(如应力集 中处)，为了较好地反映数据变化规律，需要采 用比较密集的网格。 在计算数据变化梯度较小的部位，为减小模 型规模，则应划分相对稀疏的网格。
这样，整个结构便表现出疏密不同的网格 划分形式。
1、结构的离散化 2、单元特性分析 3、计算单元刚度矩阵 4、单元集成 5、施加边界条件 6、求解位移 7、求解应力应变等场量
有限元入门ppt课件
有限体积法 (Finite Volume Method)
其基本思路是：将计算区域划分为一系列不重复的控制体积，并使每个网格点周围有一个控制体积；将待解的微分方程对每一个控制体积积分，便得出一组离散方程。其中的未知数是网格点上的因变量的数值。为了求出控制体积的积分，必须假定值在网格点之间的变化规律，即假设值的分段的分布的分布剖面。
1-2 应力的概念
作用于弹性体的外力(或称荷载)可能有两种： 表面力，是分布于物体表面的力，如静水压力，一物体与另一物体之间的接触压力等。单位面积上的表面力通常分解为平行于座标轴的三个成分，用记号 来表示。 体力，是分布于物体体积内的外力，如重力、磁力、惯性力等。单位体积内的体力亦可分解为三个成分，用记号X、Y、Z表示。 弹性体受外力以后，其内部将产生应力。
边界元法 （Boundary Element Method）
边界元法是一种继有限元法之后发展起来的一种新的数值方法，与有限元法不同，边界元法仅在定义域的边界划分单元，用满足控制方程的函数去逼近边界条件。所以边界元与有限元相比具有单元和未知数少、数据准备简单等优点，但边界元法解非线性问题时，遇到同非线性项相对应的区域积分，这种积分奇异点处的强烈的奇异性，使求解遇到困难。边界元法在塑性问题中应用还比较少。
弹性力学 — 区别与联系 — 材料力学 弹性力学与材料力学既有联系又有区别。它们都同属于固体力学领域，但弹性力学研究的对象更普遍，分析的方法更严密，研究的结果更精确，因而应用的范围更广泛。 弹性力学 固有弱点： 由于研究对象的变形状态较复杂，处理的方法又较严谨，因而解算问题时，往往需要冗长的数学运算。但为了简化计算，便于数学处理，它仍然保留了材料力学中关于材料性质的假定：
塑性有限元常用软件
有限元经典PPT第4章
Pii Kiiui
Ki1u1 Ki2u2 Kiiui K u i,i1 i1
ui
n
Kiiui Kiiui
Kiju j
4.1.2 平面应力问题有限元的基本思想和瑞雷－里兹法
v3 f3y
3
u3
f3x
f1y v1 u1
1 f1x
v2 f2y u2
2 f2x
给定一个三角形单元和作用在角点上 的六个力，要求得六个角点的位移。 或者是要求三角形角点发生指定的位 移，在三角形三个角点如何加力？
很显然，问题的精确解很困难。采用 瑞雷－里兹法求近似式解
e号单元的三个节点I,j,k的力对应的 力的平衡方程是第2i-1,2i;2j-1,2j;2k1,2k个平衡方程
e号单元的三个节点I,j,k的位移是第 2i-1,2i;2j-1,2j;2k-1,2k个未知数
弹性模量：E 横截面积：A
1
1 L
2
2L
3
局部系单元刚度阵：
k
1
EA L
1 -1
-1
1
2 集成总刚：
0 1
解得：
ux uy
L EA
3.8284L
EA
i
j
第一类位移条件：
Ki1u1 Ki2u2 Kiiui Ki1ui1
ui 0
令： Kij 0 i j
m
vi 0
Kii 1
um 0
Pi 0
ui 0
第二类位移条件：um um
大数
充大数法： Kii Kii
第一步：求转换矩阵
k2
EA 1 2L -1
-1
1
P
cos 0
T sin
第五章 薄板弯曲问题有限元讲义
第五章薄板弯曲问题有限元法第一节薄板弯曲问题的有关概念一、基本概念1．薄板的定义：薄板是由上下两个平行的表面所构成的片状结构，其间距称为板厚。
同时，定义等分板厚的面为中面，当中面为平面时，称为平板，当中面为曲面时则称为壳体。
2．挠度; 板结构在承受横向载荷（弯矩、扭矩和横向剪力）作用下，发生弯扭而使薄板中面上各个点沿垂直中面方向发生的横向变形称为挠度，记为w。
3．薄板的两类问题：（1）平面应力板问题，载荷作用于板面内—（薄膜单元）；在拉、压力和面内切力作用下，板内将产生薄膜内力，从而使板产生面内变形。
（2）薄板弯曲问题：其特点为：a) 几何尺寸：板的厚度远较长与宽的几何尺寸为小(一般厚度与板面最小尺寸之比小于1/5-1/10)；（否则称为厚板）b) 载荷条件：结构仅承受垂直于板中面的横向载荷作用。
c) 小挠度条件；即挠度与板厚之比值较小，一般为w/t ≤1/5。
研究薄板弯曲问题时，通常以未变形的板的中面为xoy平面，厚度方向为z轴方向，3．板的一般问题：一般情况下，板既可承受横向载荷作用，也可同时承受平行于板中面的膜载荷作用。
(1) 薄板：在小挠度情况下，当两种载荷同时作用时，可认为两种变形互不影响，因此膜载荷的作用可按平面应力问题进行处理，而横向载荷的作用则按薄板弯曲问题来分析，两种问题引起的薄膜内力和弯曲内力的叠加便是一般载荷综合作用的结果。
（2）厚板：当1<w／t<5时为大挠度板,w／t≥5时为特大挠度板。
在大挠度情况下，薄板面内变形和弯扭变形之间将相互影响，即横向载荷也可能产生膜内力和面内变形，而膜载荷也可能产生弯曲内力和弯曲变形。
这时描述薄板变形的数学方程是非线性的，应采用更复杂的理论分析方法。
二．薄板弯曲问题求解的假设：（克希霍夫假设）1．法线假设垂直板中面的法线在板变形后仍垂直于弯曲的挠曲面，且法线线段没有伸缩，板的厚度无变化。
这样，垂直于中面的正应变便可忽略，即εz=0根据几何方程，可得因此挠度只是x,y的函数，表示为w=w(x,y)，也即薄板中面上法线的各点都有相同位移。
有限元讲义3-2
y y
A-17
第九节 有限元法分析的步骤
一、单元刚阵的推导 1 写出位移函数； 2 计算单元应变； 3 计算单元应力； 4 根据虚功方程，得出单元刚阵。 二、整体结构的分析 1 建立整体结构的静力平衡方程式； 2 进行边界条件处理； 3 解方程组，求节点位移； 4 根据节点位移求应变、应力。
u ( x, y) Niu i N j u j N k uk Nl ul
令ui 1, u j uk ul 0, 代如上ux, y 式
v( x, y) Ni vi N j v j N k vk Nl vl
u( x, y) Ni
综上对三种单元的分析，可以看出，形状函数是单元一些 可能位移的方程式。 • 二、形状函数的性质
性质1：任一形状函数在各节点处的值或为1或为0
1 Ni 0
在节点i处 在其它节点处
A-5
性质2：单元的各个形状函数之和总是等于1
Ni 1
这两个性质的意义是：第一，形状函数反应了相邻单元在共同节点 处位移的连续性；它反映了单元的刚体位移。 • 三、形状函数的设定的说明 形状函数既然是单元的一些可能产生的位移，因此它们与位移函数 具有相同的特性，可以用插值多项式来设定。设定时要满足上诉形 状函数性质以及连续性和常应变条件。即 1、形状函数应满足
A-15
u x x x v y 0 y xy u v y x y N i x 0 Ni y N j x 0 N j y
A-16
0 Ni y Ni x
ห้องสมุดไป่ตู้
0 N j y N j x
有限元讲义3-1
[
]
试中，[B]称为三角板单元的应变矩阵或几何矩阵。
四、根据物理方 由物理方程 程求应变
{σ} = [D]{ε} = [D][B]{q}
（4.3-8）
弹性矩阵[D]是常数矩阵，[B]也是常数矩阵,因此，当节点位移 求出后，就可以算出三角板单元的应力（常数值）。 在单元内，应变和应力均为常数值，一般是与实际情况不相符 合的，当单元划分相当小时，也只能说是近似的。 五、单元力的平衡方程 在弹性力学中，应力与体积力之间的平衡关系是由平衡微分方程来 体现；应力与表面力之间的平衡关系由静力边界条件来体现，以上可 统称为应力与外力之间的平衡方程。这种平衡关系在整个弹性体内是 逐点满足的。 在有限元法中，应力与外力之间的平衡关系不是逐点满足的，而是 在单元整体意义上满足平衡。通常用虚功方程代替平衡方程。
四、根据物理方程求应力
{σ} = [D][B]{q}
五、矩形板单元刚度矩阵的导出
− − [K] = ∫∫∫V [B]T [D][B]dV = t ⋅ ∫∫S [B]T [D][B]dxdy= t ⋅ ∫a a ∫b b[B]T [D][B]dxdy
（4.4-11）
vi = a5 − aa6 − ba7 + aba 8 v j = a5 + aa6 − ba7 − aba 8 vk = a5 + aa6 + ba7 + aba 8 vl = a5 − aa6 + ba7 − aba 8
（4.4-2）
写成矩阵形式并求出多项式系数，有
1 1 1 1 1 1 u 1 1 a1 − i − a a a a u a2 1 1 1 1 j = 1 − u a3 4 − b b b b k a4 1 1 1 1 ul − − ab ab ab ab
有限元法基础ppt课件
有限单元法
一、数值模拟方法概述 二、有限单元法简介 三、有限单元法分析步骤 四、利用有限元软件进行工程分析
一、数值模拟方法概述
工程技术领域中的许多力学问题和场问题，如固 体力学中的位移场、应力场分析、电磁学中的电磁 分析、振动特性分析、热力学中的温度场分析，流 体力学中的流场分析等，都可以归结为在给定边界 条件下求解其控制方程的问题。
结构矩阵分析方法认为：整体结构可以看作是由有限 个力学小单元相互连接而组成的集合体，每个单元的 力学特征可以看作建筑物的砖瓦，装配在一起就能提 供整体结构的力学特性。
结构矩阵分析方法分析的结构本身都明显地由杆件组 成，杆件的特征可通过经典的位移法分析建立。
虽然矩阵位移法整个分析方法和步骤都与有限单元法 相似，也是用矩阵来表达、用计算机来求解，但是它 与目前广泛应用的有限单元法是有本质区别的。
❖ 国际上早在20世纪50年代末、60年代初就投入大量的人力和 物力开发具有强大功能的有限元分析程序。其中最为著名的是 由美国国家宇航局（NASA）在1965年委托美国计算科学公司 和贝尔航空系统公司开发的NASTRAN有限元分析系统。该系 统发展至今已有几十个版本，是目前世界上规模最大、功能最 强的有限元分析系统。
有限元法
既可以分析杆系结构，又分析非杆系的连续 体结构。
三、有限单元法简介
有限单元法的常用术语：
有限元模型 是真实系统理想化的数学抽象。
定义
真实系统
有限元模型
自由度（DOFs－ degree of freedoms）
自由度(DOFs) 用于描述一个物理场的响应特性。
UY ROTY
ROTZ UZ
UX ROTX
目前在工程技术领域内常用的数值模拟方法有： 1、有限单元法FEM（ Finite Element Method） 2、边界元法BEM（Boundary Element Method ） 3、有限差分法FDM（ Finite Difference Method 4、离散单元法DEM（Discrete Element Method） 其中有限单元法是最具实用性和应用最广泛的。
有限元方法讲义
有限元方法讲义第1讲抛物问题有限元方法1、椭圆问题有限元方法考虑椭圆问题边值问题：（1）问题（1）的变分形式：求使满足（2）的性质，广义解的正则性结果。
区域的剖分，矩形剖分，三角剖分，剖分规则，正则剖分条件，拟一致剖分条件。
剖分区域上分片次多项式构成的有限元空间。
的逼近性质，逆性质：这里，为的插值逼近。
问题（2）的有限元近似：求使满足（3）（3）的解唯一存在，且满足。
（3）的解所满足的矩阵方程（离散方程组）形式：（4）刚度矩阵的由单元刚度矩阵组装而成。
模误差分析：由（2）－（3）可得（5）由（5）可首先得到则得到（6）－模误差分析设满足用与此方程做内积，由（5）式和插值逼近性质得到再利用模误差估计结果，得到（7）最优阶误差估计和超收敛估计概念。
当与时间相关时（为抛物问题准备），由（5）式得（8）利用（7），类似分析可得（9）2、抛物问题半离散有限元方法考虑抛物型方程初边值问题：（10）（10）的变分形式：求使满足（11）（11）的半离散有限元近似：求使满足（12）令，代入（12），依次取可导出常微分方程组：（13）其中为质量矩阵，K为刚度矩阵。
求解常微分方程组（13），得到代回的表达式，即得半离散有限元解。
定理1.问题（12）的解唯一存在且满足稳定性估计：（14）证明：在（12）中取得到整理为（注意是正定的）对此式积分，证毕。
误差分析。
引进解的椭圆投影逼近：满足（15）根据椭圆问题的有限元结果可知（16）分解误差：的估计由（16）式给出，只须估计。
由（11），（12）和（15）知，满足取，类似稳定性论证可得可取为的投影，插值逼近等。
由（17）式，三角不等式和（16），得到（18）3、抛物问题全离散有限元近似剖分时间区间：。
引进差分算子：规定，当为连续函数时，，则有由此得到（19）（20）定义问题（11）的全离散向后Euler有限元近似：求，使满足（21）将代入（21）可导出全离散方程组（22）其中。
有限元课件ppt
将所有单元的刚度矩阵依照一定的方式组合起来，形成整体的刚度 矩阵。
载荷向量与束缚条件
载荷向量
表示作用在结构上的外力，包括集中力和散布力。
束缚条件
表示结构在某些结点上的位移受到限制，常见的束缚有固定束缚、 弹性束缚等。
载荷向量和束缚条件的引入
在建立整体刚度矩阵后，需要将载荷向量和束缚条件引入到整体刚 度矩阵中，形成完全的线性方程组。
并行计算
采取并行计算技术，提高计算效率。
算法改进
优化算法，提高计算精度和效率。
06 有限元分析软件 介绍
ANSYS
01
功能特点
ANSYS是一款功能强大的有限元分析软件，广泛应用于结构、流体、
电磁等多种工程领域。它提供了丰富的建模工具和求解器，能够处理复
杂的工程问题。
02
优点
ANSYS具有友好的用户界面和强大的前后处理功能，使得建模和网格
有限元法的应用领域
结构分析
有限元法在结构分析中应用最 为广泛，可以用于分析各种类 型的结构，如桥梁、建筑、机
械零件等。
热传导
有限元法可以用于求解温度场 的问题，如热传导、热对流和 热辐射等问题。
流体动力学
有限元法在流体动力学领域也 有广泛应用，可以用于求解流 体流动和流体传热等问题。
其他领域
除了上述领域外，有限元法还 广泛应用于电磁场、声场、化
学反应等领域。
02 有限元的数学基 础
线性代数基础
向量与矩阵
01
介绍向量的基本概念、向量的运算、矩阵的表示和运算规则等
。
线性方程组
02
论述线性方程组的解法，包括高斯消元法、LU分解等。
特征值与特征向量
张年梅有限元方法讲义
张年梅有限元方法讲义全文共四篇示例，供读者参考第一篇示例：张年梅有限元方法讲义有限元方法是一种非常重要的数值计算方法，广泛应用于力学、电磁学、声学、地球物理学等领域。
张年梅是中国工程院院士、有限元方法的权威专家，他在有限元方法的研究和应用方面取得了很多成果。
他的有限元方法讲义成为了很多工程学子和研究人员学习的重要参考资料。
有限元方法是一种用数值方法解决复杂工程问题的工具。
它将实际工程问题抽象为有限个简单形状的单元，并通过适当的数学方法和计算机程序求解得到问题的近似解。
有限元方法的基本思想是将一个复杂的结构或领域分割成有限个简单的子结构或子域，然后在每一个子结构或子域上建立合适的数学模型，最后通过组合所有子结构或子域的模型获得整体结果。
张年梅有限元方法讲义详细介绍了有限元方法的基本原理、数学模型的建立和求解方法。
讲义先介绍了有限元方法的起源和发展历程，然后对基本概念和术语进行了解释，包括有限元模型、单元、节点、网格等。
接着讲义详细介绍了有限元方法的基本原理，包括离散化、变分原理、加权残差法、Galerkin法等。
有限元方法的数学模型的建立是有限元分析的关键步骤。
张年梅有限元方法讲义介绍了常见的结构、固体、流体、电磁等问题的有限元建模方法，包括线性弹性分析、非线性分析、热传导分析、流体动力学分析等。
在建立数学模型之后，有限元方法的求解方法也是十分重要的。
张年梅有限元方法讲义介绍了有限元方法的常用数值解法，包括直接法、迭代法、有限元展开法等。
有限元方法在实际工程问题中有着广泛的应用。
张年梅有限元方法讲义通过大量的案例和实例展示了有限元方法在结构分析、热力分析、电磁分析等领域的应用。
讲义还介绍了有限元方法在工程设计和优化中的应用，包括拓扑优化、材料优化、结构优化等。
张年梅有限元方法讲义是一部权威的、全面的有限元方法教材，受到了广大工程学子和研究者的欢迎和好评。
通过学习这本讲义，读者可以系统地了解有限元方法的基本原理和求解方法，掌握有限元方法在工程问题中的应用技能，为解决工程问题提供强有力的工具支持。
有限元法PPT课件
Motorola– Drop Test Fujitsu-Computers Intel –Chip Integrity
电子
Baxter - Equipment J&J – Stents Medtronic - Pacemakers
医疗
Principia-spain Arup-U.K. T.Y. Lin - Bridge
有限元法
左图所示,为分析齿轮上一个齿内的应力分布,可分析图中所示的一个平面截面内位移分布.作为近似解,可以先求出图中各三角形顶点的位移.这里的 三角形就是单元,其顶点就是节点。
从物理角度理解, 可把一个连续的齿形截面单元之间在节点处以铰链相链接,由单元组合而成的结构近似代替原连续结构,在一定的约束条件下,在给定的载荷作用下,就可以求出各节点的位移,进而求出应力.
一.Abaqus公司简介
公司
’00 ’01 ’02 ’03 ’04 ‘05 ’06 ‘07
18%
18%
20%
SIMULIA公司（原ABAQUS公司）成立于1978年，全球超过600名员工，100% 专注于有限元分析领域。 全球28个办事处和9个代表处 业务迅速稳定增长，是当前有限元软件行业中唯一保持两位数增长率的公司。 2005年5月ABAQUS加入DS集团,将共同成为全球PLM的领导者
Where :
Displacement interpolation functions (位移插值函数)
13.3 Approximating Functions for Two-Dimensional Linear Triangular Elements (二维线性三角形单元的近似函数)
node (节点)
element(单元)
有限元分析(FEA)方法PPT课件
(b)定义几何模型 应用实体建模
(c) 用P单元分网。 自适应网格对P方法是无效的
3.施加载荷、求解
应用实体模型加载，而不是有限元模型
求解：推荐采用条件共轭梯度法（PCG），但PCG对于壳体P单元无效
4.后处理 察看结果
有限元分析及应用讲义
举例： platep.dat
20 in
R=5 in
SEQV SMN=773.769 SMNB=708.94 SMX=4421 SMXB=4999
有限元分析及应用讲义
P方法及p单元的应用
P 单元的位移形函数
u=a1+a2x+a3y+a4x2+a5xy+a6y2
v=a7+a8x+a9y+ a10x2+a11xy+a12y2
P方法的优点：
如果使用 p-方法 进行结构分析，可以依靠p单元自动调整单元多项式阶数（2-
– 导出 MeshTool 工具, 划分方式设为自由划 分.
– 推荐使用智能网格划分 进行自由网格划分, 激活它并指定一个尺寸级别. 存储数据库.
– 按 Mesh 按钮开始划分网格. 按拾取器中 [Pick All] 选择所有实体 (推荐).
– 或使用命令 VMESH,ALL 或 AMESH,ALL.
savg = 1100
s = 1000 Elem 1
s = 1100
s = 1200 Elem 2
s = 1300
(节点的 ss 是积分点 的外插）
savg = 1200
7
有限元分析及应用讲义
ANSYS网格误差估计
误差估计作用条件：
• 线性静力结构分析及线性稳态热分析 • 大多数 2-D 或 3-D 实体或壳单元 • PowerGraphics off
有限元法基础讲稿-第8讲
青岛大学讲稿讲 授 内 容备 注 第8讲(第11周)2. 应变矩阵确定了单元位移后，可以很方便地利用几何方程和物理方程求得单元的应变和应力。
作为平面问题，单元内具有3个应变分量εx 、εy 、γxy （各符号的意义见附录1），用矩阵表示为⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡∂∂+∂∂∂∂∂∂=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=x v y u y v x u xy y x γεεε将（2.1.4）式代入上式中，得到e m mjjiim j i m j ib c b c b c c c c b b b A δε⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=00000021 或eδB ε = (2-1-7)式中B 称为应变矩阵，写为分块形式，即B =[B i B j B m ] (2-1-8)而其子阵为),,( 0021m j i b c c b A i ii ii ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=B (2-1-9)3节点三角形单元的B 是常量阵，所以称为常应变单元。
在应变梯度较大(也即应力梯度较大)的部位，单元划分应适当密集，否则将不能反映应变的真实变化而导致较大的误差。
上述应变中包括与应力有关的应变和与应力无关的应变两部分，无关的应变ε0又称为初应变⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=0000xy y x γεεεε0由温度变化、收缩、晶体生长等因素引起，对工程结构一般只考虑温度应变，无论线性和非线性温度，计算时可近似地采用平均温度33refT T T T T m j i -++=式中，T i 、T j 、T m 分别为节点i 、j 、m 的温度，T ref 为参考温度。
对于平面应力问题，温度T 引起的初始应变为⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=00TT ααε其中，α为线膨胀系数。
由于温度变化在各向同性介质中不引起剪切变形，所以γxy 0=0。
以后所述问题，除非特别说明，都指各向同性介质。
对平面应力问题，温度T 引起的初始应变为⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+=0)1(0TT ααμε 当不考虑温度的影响时，当前温度即为参考温度。
有限元ppt课件
因此有 y(x) (x)
试探函数中所取的项数越多，逼近的精度越高。
将试探函数代入式(1-9)，可以得到关于n个待定系数
的泛函表达式，简记为 I y(x) I(1,2,3, ,n)
根据多元函数有极值的必要条件，有
1
I (1,2 ,3,
2
I (1,2 ,3,
力,它反映了内力在截面上的分布密度。
z
y
o
zx
xz
z zy
yz
切应力互等定律 xy yx , xz zx , yz zy
y
应力矩阵
x xy
yx
T
x y z xy yz zx
y
x
z
微分体的应力分量
v y w z u v
0
0
yz
zx
y x y
v
w
0
y
0
x
0
z
u v
0
w
39
厚度为1的微分体,在水平方向拉
力F的作用下发生了位移 xdx
拉力表达式:
F xdy 1
x
x dy
拉力做的功:
dx
xdx
dW
1 2
F xdx
将F代入:
dW
1 2
x
x
dxdy
40
储存在微分体内的应变能:
张年梅有限元方法讲义
张年梅有限元方法讲义一、啥是有限元方法呀。
小伙伴们！今天咱来唠唠这有限元方法。
有限元方法呢，简单来说，就是一种数值计算方法啦。
它就像是一个超级厉害的数学工具，能帮咱们解决好多复杂的工程问题和物理问题呢。
比如说在结构力学里，要算一个复杂形状的结构受力情况，用传统方法可能会把人脑袋算大，这时候有限元方法就派上用场啦，轻松就能搞定。
想象一下哈，把一个复杂的大物体，划分成好多好多小的单元，就像把一个大蛋糕切成好多小块一样。
每个小单元呢，它的性质相对就比较简单啦，咱们就可以对这些小单元进行分析计算，最后再把这些小单元的结果组合起来，就得到整个大物体的结果啦，是不是挺神奇的！二、有限元方法的发展历程。
有限元方法的发展那可是经历了不少阶段呢。
早期的时候，它的理论和算法都还不太完善，计算能力也有限，只能解决一些相对简单的问题。
但是随着科技的不断进步，计算机的计算速度越来越快，存储能力也越来越强，有限元方法就像开了挂一样，发展得越来越厉害啦。
现在啊，有限元方法已经广泛应用到了好多领域。
像航空航天领域，设计飞机、火箭这些复杂的飞行器，就离不开有限元方法来分析它们的结构强度和气动性能；在汽车制造领域，设计汽车的车身、发动机等部件，也得靠有限元方法来优化设计，提高安全性和性能；还有土木工程领域，建高楼大厦、修桥梁隧道，都得用它来确保结构的安全可靠。
三、有限元方法的基本步骤。
这有限元方法要想用好，得按照一定的步骤来哦。
1. 建模。
首先得建立一个数学模型，把实际问题转化成数学问题。
比如说要分析一个汽车零部件的受力情况，就得把这个零部件的形状、尺寸、材料特性等信息都用数学的方式表示出来。
这就好比是给这个零部件画一幅精确的“画像”，让计算机能明白它长啥样，有啥特点。
2. 离散化。
接下来就是把这个模型划分成好多小单元啦，就像前面说的切蛋糕一样。
这个过程要根据问题的复杂程度和精度要求来确定单元的大小和形状。
一般来说，问题越复杂，要求精度越高，单元就得划分得越小越细。