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