《有限元分析的数学原理》由会员分享,可在线阅读,更多相关《有限元分析的数学原理(77页珍藏版)》请在金锄头文库上搜索。
1、Chapter3Chapter3有限元分析的数学求解原理有限元分析的数学求解原理o 一般说来,求解方程的途径有两大类:1)直接针对原始方程进行求解2)间接针对原始方程进行求解直接解法解析法:o解析法从力平衡关系、几何关系以及物理关系出发,推导出一个 或一组关于应力或者关于应变、有时是同时含有应力、应变的微 分方程或偏微分方程,通过求解微分方程,解出应力、应变和变 形量。工程中,常采用的解析方法有材料力学中对杆件的分析, 弹性力学中平面问题的求解,板壳理论等。 o解析法的很多基本理论是建立在一些简化的假设基础之上的,经 过大量的工程实践,被证明能很好的符合构件实际工作情况,已 成为成熟的理论。解
2、析法得到的结果是未知量(应力、应变等) 的函数解,可直接得到结构中任意点的精确解。解析法在分析理 论问题以及一些工程问题时起着重要作用。但是解析法在应用到 一些形状复杂或应力分布复杂的结构时,往往由于数学上的问题 而显得无能为力,因而使解析法在应力分析中的应用受到限制。根据问题的性质,确定基本未知量和相应的基本方程,并且假设一组满足全部基本方程的应力函数或 位移函数。然后在确定的坐标系下,考察具有确定的 几何尺寸和形状的物体,其表面将受什么样的面力作 用或者将有什么样的位移。直接解法逆解法:对于给定的弹性力学问题,根据弹性体的几何形状,受力特征和变形特点,或已知简单结论,如材料力学 解,假设部
3、分应力分量或者部分位移分量的函数形式为 已知,由基本方程确定其他的未知量,然后根据边界条 件确定未知函数中的待定系数。直接解法半逆解法:逆解法和半逆解法的应用将在以后的章节中介绍,其求解过程带有 “试算”的性质。直接解法有限差分法:o有限差分法:微分方程和积分微分方程数值解的方法。其基本思 想是:有限差分方法(finite difference method )是计算机数值模拟最早采用的方法,至今仍被广泛运用。该方法将求解域划分为差分网格, 用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代 替进行离散,从而建立以网格节点上
4、的值为未知数的代数方程组 。该方法是一种直接将微分问题变为代数问题的近似数值解法, 数学概念直观,表达简单,是发展较早且比较成熟的数值方法。 有限差分格式o格式精度:一阶格式、二阶格式和高阶格式。o差分的空间形式:中心格式o时间因子:显格式、隐格式、显隐交替格式等。o构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方 法。其基本的差分表达式主要有三种形式:一阶向前差分、一阶 向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为 一阶计算精度,后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。 间接解法加权残值法:o 是一种应用广泛的求解
5、微分方程的方法.该方法先假定 一族带有待定参数的定义在全域上的近似函数,该近似 解不能精确满足微分方程和边界条件,即存在残差.在加 权平均的意义下消除残差,就得到加权残值法的方程.由 于试函数定义在全域上,所得方程的系数矩阵一般为满 阵.选取不同的权函数,可得到不同的加权参量法. 虚功原理定义:弹性体处于平衡状态,对于满足变形连续条件的 虚位移及其虚应变,外力在虚位移上所做的虚功,等于真实应力 分量在对应的虚应变上所做的虚功,即虚应变能。最小势能原理要求最后得间接解法虚功原理:应变 能 应变 余能 应变能应变余能间接解法最小势能原理:有限元上的应用(位移法):假设单元位移模式 单元刚度方程o
6、把一个物理学问题用变分法化为求泛函极值(或驻值) 的问题,后者就称为该物理问题 的变分原理。如果建 立了一个新的变分原理,它解除了原有的某问题变分原 理的某些约束条件,就称为该问题的广义变分原理;如 果解除了所有的约束条件,就称为无条件广义变分原理 ,或称为完全的广义变分原理。 o 在当代,变分原理已成为有限元法的理论基础,而广义 变分原理已成为混合和杂交有限元的理论基础。在实际 应用中,通常很少能求出精确的解析解,因此大多采用 近似计算方法。间接解法变分原理:1)假定2)将上式代入泛函 ,计算变分 。3)由极值条件,算出待定常数 ,使之满足基本微分方程。4)把得到的常数代回 ,得到所求问题的
7、解。与有限元方法比较:相同点:都是求解极值问题的方法,方法类似。不同点:求解问题区域不同:局部和整体 关系。对于泛函工程问题无论是几何形状、受力方式还是材料特性都是前变万化的,因此一种求解方法是否有优势,其判断标 准应该是p 具有良好的规范性(不需要太多的经验和个人技巧) p 具有良好的适应性(可以处理任何复杂的工程实际) p 具有良好的可靠性(计算结果收敛稳定,精度高) p 具有良好的求解可行性(计算工作量) 本章主要内容o 3.1简明问题的解析求解o 3.2弹性力学问题近似求解的加权残值法o 3.3弹性力学近似求解的虚功原理、最小势能原理及其变分原理o 3.4各种求解方法的特点及比较3.1
8、简明问题的解析求解一维拉压杆问题o 基本变量:ux(x),x(x),x(x)o 基本方程:几何方程物理方程平衡方程边界条件o 对三大方程直接进行求解得几何方程物理方程平衡方程根据边界条件可得,c=P/A, c1=0o 讨论1若用材料力学的经验方法求解,则需先作平面假设,即假设应力为均匀分布x=P/A由广义胡克定律得x=P/EA右端的伸长量为 u=xl=Pl/EAo 讨论2应变能动能势能有限元分析步骤-单元分析o 由于杆单元只有两个节 点位移,故可以设杆单 元的位移模式为之包含 两个待定常数的形式 u(x)=a1+a2x根据有限元法的基本思路,将弹性体离散成有限个单元体的组合,以 结点的位移作为
9、未知量。弹性体内实际的位移分布可以用单元内的位 移分布函数来分块近似地表示。在单元内的位移变化可以假定一个函 数来表示,这个函数称为单元位移函数、或单元位移模式。 o 回代得 写成矩阵形式为其中Ni,Nj是形函数。o 根据位移条件有u(0)=u0, u(l)=ul,从而得o根据几何方程得p根据物理方程得p从而,根据单元分析结果,进行整体分析,求解整体方程组, 进行结果分析由虚功原理可以推得 3.1简明问题的解析求解平面梁问题o 基本方程1)一般的建模及分析方法,取微单元体2)特征建模法,采用工程宏观特征量来进行问题描述简支梁的特征:o 梁为细长梁,因此可以只用x坐标来刻画;o 主要变形是垂直于
10、x的挠度,可以只用挠度来描述位移场。针对这两个特征,可以做出以下两个假定:p 直法线假定p 小变形和平面假设o 直法线假定:一垂直中面的直线(称为法线),变形时不 伸缩,并且仍为弹性曲面的法线。 o 平截面假定:平截面假定是材料力学中最基本的假定之 一。这个假定认为所有与杆轴线垂直的截面在杆件变形后仍保持为平面。这样截面上每一个点的变形趋势就可 以确定,如果知道了中性轴的位置和任意一点的应变( 变形),整个断面的应变就可以知道,这是建立该假设 的基础。实验也证明匀质弹性体根据此项假定所得的计 算结果是准确的。 o 基于以上假定,该问题的三类基本变量为位移:中性层的挠度v(x=y=0) 应力:x
11、方向的正应力x,其他应力分量很小忽略不 计,该变量对应于梁截面上的弯矩M应变:采用x ,满足直法线假定o 平衡方程:x方向y方向o 几何方程:ab的变形为:因此正应变为: 其中为曲率半径。o 物理方程:由广义虎克定律有整理得p 边界条件。弯矩ox方向平衡y方向平衡。o 求解方程 得其中c1,c2,c3是待定系数。最后可得o 讨论应变能外力功势能由于单元有四个位移 分量,可设梁单元的 位移模式v(x)为包含4个待定常数的三次多 项式:有限元分析步骤-单元分析根据边界条件可以确定待定系数,将其进一步回代,可以得到用节点位移表示的梁单元位移。 式中根据梁的平面假定可知梁单元的轴向应变为:这里利用平面
12、假设(变形后横截面仍保持平面,与纵线正交)如图: 从而可以由单向虎克定律得出单元的轴向应力:由虚功原理可以推得 3.2弹性力学问题近似求解的加权残值法直接针对原始三大类方程边界条件下求解三大类变量往往是非常困那的,尤其是当几何形状和边界条件比较复杂时,一般求不出相应的解析解。如果事先假定满足一定边界条件的试函数,再在此基础上进行近似求解,则可以大大降低求解难度。这种试函数方法可以使得求解过程比较规范和简单,并有一定的适应性,但是求解的精度有所降低。试函数方法的基本原理:先假定满足一定边界条件的试函数,然后将其带入需要求解的方程中(控制方程),通过使与原来的方程的误差残值最小来确定试函数中的待定
13、系数。为了提高解或逼近精度,可以采用较多项数的试函数来进行计算,这种方法叫做加权残值法。加权残值法WRM:Galerkin加权残值残值最小二乘法3.2.1梁弯曲问题近似求解的Galerkin加权残值法设满足以下方程和边界条件的位移场为公式中的L为微分算子。由于平面弯曲梁的平衡方程为故假设能找到事先满足式中的边界条件的一个试函数,将其带入到控制方程,则一定存在残差,记为对于更一般的情形,设有一组满足所有边界条件的试函数,将其线性组合为新的试函数其中c1,c2,c3cn为待定系数。将试函数代入原始方程组,则必有残差 ,真实的c1,c2,c3cn使得残值的积分为零,即其中w1,w2,w3wn为权函数
14、。以上为关于c1,c2,c3cn的方程组,由上式可以求出他们,最后由线性组合形式的试函数得到真实的 。如果将权函数w1,w2,w3wn取为1, 2, 3 n,则该方法称作伽辽金法。 受均布外载荷简支梁的Galerkin加权残值法求解代入控制方程得残差由Galerkin加权残值方程分析可得 ,求解上式可得受均布外载荷简支梁的Galerkin加权残值法求解代入控制方程得残差由Galerkin加权残值方程分析可得 ,求解上式可得受均布外载荷简支梁的Galerkin加权残值法求解同时满足面力边 界条件根据Galerkin法分析可得 几种函数结果比较 1、仅仅取1项试函数时,由伽辽金加权残数法得到的结果
15、与精确 解得相对误差为0.3861%。2、仅仅取2项试函数时,由伽辽金加权残数法得到的结果与精确 解得相对误差为-0.027%。由以上比较可以看出,此解法精度还是比较高的。3.2.2梁弯曲问题近似求解的最小二乘法同样,设有满足所有边界条件的试函数,若将试函数代入原始方程中,则必有残差值,真实地待定常数使得残差值平方的加权积分取极小值,即其中w为权函数,一般可以取1。将上式进一步具体化,即对上式取极值,有这是一组关于c1,c2,c3cn的方程组,由上式可以求出它们,从而得到 真实的试函数。受均布外载荷简支梁的残值最小二乘法求解将其带入到控制方程,则一定存在残差,取权函数为1,则残差平方的积分为由
16、最小二乘法有由上式可以解出与用伽辽金加权残值法相同的结果。3.2.3一般弹性问题近似求解的加权残值法平面问题有:三大应力分量,三个应变分量,以及两个位移分量。三个物理方程,三个几何方程,两个平衡方程,三类边界条件基于三大物理量之间的关系(见三大方程),可以得到基于位移表达的平衡方程(控制方程)如下:再加上边界条件,就构成了基于位移求解的平面应力问题的所有方程。设有位移的试函数为式中的cui和cvi是待定常数,ui和vi是满足所有边界条件的基底函数。将试函数代入控制方程,使得其残值在加权积分下为零,即可得到Garlerkin加权残值方程,得到一组关于待定系数的方程组。同样,也可以得到空间问题的加权残值法,包括伽