有限元与数值方法-讲稿5-2 有限元的基本概念

有限元与数值方法第五讲

有限元法的一般原理与基本格式

----有限元的基本概念授课教师：刘书田Tel：84706149；Email：stliu@教室：综合教学楼351时间：2013年4月12日：8：00—10：201弹性力学问题的有限元法有限元法的基本思想杆系结构的直接刚度法静定桁架的内力可以通过节点的平衡方程求得，由内力和杆件断面积可求得杆件应力、应变，再求得节点位移PP静不定桁架的内力无法简单通过节点平衡方程求得，需要采用力法或位移法求得。采用位移法求解时，假定每个节点的位移为未知量，然而可以将杆件伸长、杆件应变、杆件应力，杆件内力用节点位移表示，根据节点的平衡要求可以得到节点位移满足的平衡方程。2有限元法的基本思想杆系结构的直接刚度法静不定桁架的内力无法简单通过节点平衡方程求得，需要采用力法或位移法求得。采用位移法求解时，假定每个节点的位移为未知量，然而可以将杆件伸长、杆件应变、杆件应力直至杆件内力用节点位移表示，根据节点的平衡要求可以得到节点位移满足的平衡方程。由节点的平衡方程就可求得节点位移；这一平衡方程的系数矩阵就是结构刚度矩阵；结构刚度矩阵是由每个杆件的单元刚度矩阵适当地组装得到。F2x,,u2xF2y,,u2yF1x,,u1xF1y,,u1y12P3杆系有限元方法以桁架结构为例，介绍有限元的基本思想4杆单元的有限元分析一维线性杆单元基本假定：只能承受拉压内力（各杆两端的约束条件使得弯曲、扭转、剪切不能传递）轴线为直线材料满足胡克定律自由转动12桁架结构5位移插值建立轴线方向的坐标系记任一点轴向位移为并将节点位移表示为建立杆件位移与节点位移的插值关系其中，形函数必须满足112121节点位移协调关系满足6可简单地将形函数取为一次多项式的形式：考虑到边界条件，可得到因此位移插值杆上无分布力时，一次多项式可精确描述杆件变形7位移及应变小位移假设下，应变为位移模式为位移模式包括刚体位移和常应变模式N形函数矩阵B应变矩阵8单元刚度阵利用胡克定律，得到杆件应力和内力分别为则节点力为其矩阵形式表示为

单元刚度矩阵S应力矩阵9XYxyXYxyi坐标变换矩阵设OXY为结构坐标，oxy为单元坐标。为任意单元

i

端的任一矢量。它在结构坐标系中的分量为X、Y；在单元坐标系中的分量为x、y。X、Y

在单元坐标x轴上投影的代数和给出x

。同理，X、Y

在单元坐标y轴上投影的代数和给出y

10即坐标变换矩阵令表示两个端点的位移矢量在单元局部坐标系的分量，表示两个端点的位移矢量在全局坐标系的分量，则11上式可写成坐标变换矩阵[R]的具体内容为：用节点坐标描述方向余弦：坐标变换矩阵（Xi，Yi）和（Xj，Yj）分别为节点i

和节点j

在全局坐标系中的坐标值12平面内任意方向的杆单元记为而节点力列阵满足(或)由单元局部坐标系下的关系可得到或写成其中131.整体节点位移单元节点位移：总体控制方程：单元集成分析桁架结构扩充矩阵2.整体节点力14边界条件全局平衡方程如不考虑约束条件，总刚度阵是奇异的零位移约束条件15边界条件处理零位移约束条件代人平衡方程，得到约束反力外载荷未知位移16对于一般的指定位移约束，可将方程分块为其中，是指定位移，是主动位移边界条件即17在单元局部坐标系中的单元节点位移分量为根据位移插值关系单元应变和应力可给出单元轴向应变为由胡克定律可进一步给出单元轴向应力为18单元应变和应力而由可得到由总体坐标系位移分量表示的单元应变和单元应力19连续体问题的有限元方法的基本思想

20有限元法(FEM)是求解偏微分方程边值问题近似解的数值方法边值问题未知量是由控制方程（椭圆、双曲、抛物型）描述的场变量（如位移、温度、流体速度等）边界条件是给定的场变量值或者其偏导数有限元法的基本概念21有限元法的基本概念有限元分析的基本思想是将求解域场分成小的子区域，通常称为“单元”或“有限元”。对每一单元假定一个分片近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。有限元法方程的系数矩阵通常是稀疏的，便于求解。有限元法不仅计算精度高，而且能适应各种复杂形状，不同物理特性、多变的边界条件和任何承载情况的工程结构分析问题。有限元法应用于场（力场、电场、磁场、温度场、流体场等）分析、热传导、非线形材料的弹塑性蠕变分析等22(a)二维问题的几何域(b)三角形单元(c)有限元网格的一部分单元有限元网格有限元法中的离散各种几何形状的有限元单元23三角形的顶点称为节点（node）节点处的场变量（这里是温度）将作为自变量被直接求解node热传导问题的三角形单元node有限元法中的场变量表示以平面热传导问题的三角形单元为例24除了节点外的其他各位置的点对应的场变量如何确定？

单元内部点的场变量值由单元节点的插值(interpolation)给出:T=?有限元法中的场变量表示

,,和是插值函数，称为位移函数或形函数。插值函数所包括的多项式阶数越高，越能精确表示位移分布。25常见平面单元形状与节点数三节点三角形单元CST(常应变单元)六节点三角形单元二次插值八节点四边形单元二次插值CST三角形单元网格划分简单，但对于弯曲过刚；线性应变三角形元描写弯曲性能远优于CST单元四边形单元剖分有时比较困难，但性能较好四节点四边形单元双线性插值26四节点四面体单元线性插值（常应变）十节点四面体单元二次插值（线性应变）八节点四面体单元Lagrange单元非完全三次插值二十节点Serendipity单元四面体单元网格剖分简单，但四节点四面体精度较差八面体单元精度较好，但网格剖分比较困难常见三维单元形状与节点数27一维单元(x)(x)(x)不同形式的单元插值28二维单元不同形式的单元插值29三维单元不同形式的单元插值(x,y,z)(x,y,z)30有限元法总体思路有限元法通过加权余量法（或变分法、最小势能原理、虚功原理等）将偏微分方程转变为代数方程，便于计算机处理将求解域剖分为网格，对节点变量进行单元分片插值值单元矩阵形成规范化，形函数只与坐标有关，便于计算机计算单元刚度阵组装为整体系数矩阵后，考虑边界条件，求解矩阵方程即可得到节点未知量。31有限元法的理论基础概述将微分方程转化为等效积分弱形式变分原理加权余量法采用单元上的分片假设近似函数，将积分方程转化为代数方程组32FiniteDifference(FD)Method:FD对微分算子进行近似在节点上建立方程FiniteElement(FE)Method:FE采用精确的算子，但利用基函数对场变量进行近似在问题域上建立方程（需要积分形式）