《论文平面问题有限元解法（公式推导讲解）课件》由会员分享，可在线阅读，更多相关《论文平面问题有限元解法（公式推导讲解）课件（79页珍藏版）》请在金锄头文库上搜索。
1、平面问题的有限单元解法,2018/9/27,有限单元法基本思想,有限单元法的思想是将物体（连续的求解域）离散成有限个且按一定方式相互联结在一起的单元组合，来模拟或逼近原来的物体，从而将一个连续的无限自由度问题简化为离散的有限自由度问题求解的一种数值分析法。物体被离散后，通过对其中各个单元进行单元分析，最终得到对整个物体的分析。有限单元法的分析步骤如下： 物体离散化 单元特性分析 单元组集，整体分析 求解未知节点的位移 由节点的位移求解各单元的位移和应力,2018/9/27,有限元单元模型中几个重要概念,单元 网格划分中每一个小的块体 节点 确定单元形状、单元之间相互联结的点 节点力 单元上节点
2、处的结构内力 载荷 作用在单元节点上的外力 （集中力、分布力） 约束 限制某些节点的某些自由度 弹性模量（杨式模量）E 泊松比（横向变形系数） 密度,单元,单元,载荷,节点,节点力,约束,2018/9/27,1. 研究内容,内容：弹性体在外力或温度作用下的应力、变形、位移等分布规律。,任务：解决弹性体的强度、刚度、稳定性问题。,弹性力学的内容及基本假定,2. 研究对象,一般弹性实体结构：,三维弹性固体、板状结构、杆件等,2018/9/27,弹性力学的内容及基本假定,3. 研究方法,由平衡方程、几何方程、物理方程三方面分析,4. 数学理论基础, 偏微分方程（高阶，二、三个变量）,数值解法：能量法
3、（变分法）、差分法、有限单元法等。,2018/9/27,弹性力学的内容及基本假定,5. 基本假定,（1）. 连续性假定,整个物体的体积都被组成物体的介质充满，不留下任何空隙。,作用：,使得、u 等量表示成坐标的连续函数。,2018/9/27,弹性力学的内容及基本假定,（2）. 完全弹性假定,假定物体完全服从虎克（Hooke）定律，应力与应变间成线性比例关系。,脆性材料 一直到破坏前，都可近似为线弹性的；,塑性材料 比例阶段，可视为线弹性的。,（3）. 均匀性假定,假定整个物体是由同一种材料组成的，各部分材料性质相同。,作用：,弹性常数（E、）等不随位置坐标而变化；,取微元体分析的结果可应用于整
4、个物体。,2018/9/27,弹性力学的内容及基本假定,（4）. 各向同性假定,（5）. 小变形假定,假定物体内一点的力学性质在所有各个方向都相同。,作用：,弹性常数（E、）不随坐标方向而变化；,假定位移和形变是微小的，即物体受力后物体内各点位移 远远小于物体的原来的尺寸。,作用：,建立方程时，可略去高阶微量；,可用变形前的尺寸代替变形后的尺寸。,使求解的方程线性化。,2018/9/27,基本概念：,外力、应力、形变、位移。,1. 外力：,体力、面力,(1) 体力, 分布在物体体积内的力, 体力分布集度,（矢量）,单位：,N/m3,kN/m3,说明：,f 是坐标的连续分布函数；,弹性力学中的几
5、个基本概念,2018/9/27,(2) 面力, 分布在物体表面的力, 面力分布集度（矢量）,单位：,1N/m2 =1Pa (帕),1MN/m2 = 106Pa = 1MPa (兆帕),说明：,弹性力学中的几个基本概念,是坐标的连续分布函数；,2018/9/27,2. 应力,(1) 一点应力的概念,内力,(1) 物体内部分子或原子间的相互作用力;,(2) 由于外力作用引起的相互作用力.,(不考虑),P,截面上P点的应力,应力矢量.,的极限方向,应力分量,应力的法向分量, 正应力,应力的切向分量, 切应力,单位:,MPa (兆帕),应力关于坐标连续分布,弹性力学中的几个基本概念,2018/9/27
6、,(2) 一点的应力状态,通过一点P 的各个面上应力状况的集合, 称为一点的应力状态,x面的应力：,y面的应力：,z面的应力：,弹性力学中的几个基本概念,2018/9/27,用矩阵表示：,其中，只有6个量独立。,切应力互等定理,应力正负号的规定：,正应力 拉为正，压为负。,切应力 坐标正面上，与坐标正向一致时为正；,坐标负面上，与坐标正向相反时为正。,弹性力学中的几个基本概念,2018/9/27,3. 形变,形变 物体的形状改变,（1）线段长度的改变,（2）两线段间夹角的改变。,P,B,C,A,用线（正）应变度量,切应变度量,（切应变两垂直线段夹角（直角）的改变量）,三个方向的线应变：,三个平
7、面内的切应变：,(1) 一点形变的度量,应变的正负：,线应变：,伸长时为正，缩短时为负；,切应变：,以直角变小时为正，变大时为负；,弹性力学中的几个基本概念,2018/9/27,(2) 一点应变状态,其中,应变无量纲；,4. 位移,注：,一点的位移 矢量S,应变分量均为位置坐标的函数,S,位移分量：,u x方向的位移 分量；,v y方向的位移 分量；,w z方向的位移 分量。,量纲：m 或 mm,弹性力学中的几个基本概念,2018/9/27,工程力学问题建立力学模型的过程中，一般从三方面进行简化：,结构简化,如空间问题向平面问题的简化，向轴对称问题的简化，实体结构向板、壳结构的简化。,受力简化
8、,如：根据圣维南原理，复杂力系简化为等效力系等。,材料简化,根据各向同性、连续、均匀等假设进行简化。,2018/9/27,基本变量 u （位移） （应变） （应力） 基本方程 力的平衡方程 几何方程 物理方程边界条件求解方法 经典解析 半解析 传统数值解法 现代数值解法（计算机硬件、规范化、标准化、规模化）,物体变形及受力情况的描述,三大方面,三大方程,即：, =E E 弹性模量,2018/9/27,平面问题的基本理论,任何一个实际的弹性力学问题都是空间问题，但是如果所考察的弹性体具有某种特殊的形状，并且承受的是某些特殊的外力和约束，就可以把空间问题简化为近似的平面问题。两种典型的平面问题 平
9、面应力问题平面应变问题,2018/9/27,平面应力问题,(1) 几何特征,一个方向的尺寸比另两个方向的尺寸小得多。, 平板,如：板式吊钩，旋转圆盘，工字形梁的腹板等,(2) 受力特征,外力（体力、面力）和约束，仅平行于板面作用，沿 z 方向不变化。,2018/9/27,(3) 应力特征,如图选取坐标系，以板的中面为xy 平面，垂直于中面的任一直线为 z 轴。,由于板面上不受力，有,因板很薄，且外力沿 z 轴方向不变。,可认为整个薄板的各点都有：,由切应力互等定理，有,结论：,平面应力问题只有三个应力分量：,应变分量、位移分量也仅为 x、y 的函数，与 z 无关。,2018/9/27,平面应变
10、问题,(1) 几何特征,水坝,滚柱,厚壁圆筒,一个方向的尺寸比另两个方向的尺寸大得多，且沿长度方向几何形状和尺寸不变化。, 近似认为无限长,(2) 外力特征,外力（体力、面力）平行于横截面作用，且沿长度 z 方向不变化。,约束 沿长度 z 方向不变化。,(3) 变形特征,如图建立坐标系：以任一横截面为 xy 面，任一纵线为 z 轴。,设 z方向为无限长，则,沿 z 方向都不变化，,仅为 x，y 的函数。,任一横截面均可视为对称面,2018/9/27,水坝,任一横截面均可视为对称面，则有,所有各点的位移矢量都平行于 x y 平面。, 平面位移问题, 平面应变问题,注：,平面应变问题中,但是，,2
11、018/9/27,如图所示三种情形，是否都属平面问题？是平面应力问题还是平面应变问题？,平面应力问题,平面应变问题,非平面问题,2018/9/27,三大基本方程,根据静力学、几何学和物理学三方面条件，建立三套方程。 平面问题中，根据微分体的平衡条件，建立平衡微分方程：(1-1)根据微分线段上形变与位移之间的几何关系，建立几何方程：(1-2) 根据应力与形变之间的物理关系，建立物理方程：(1-3)(1-3),2018/9/27,平衡微分方程,从弹性体中取出一个微分体，根据平衡条件导出应力分量与体力分量之间的关系式，也就是平面问题的平衡微分方程。 从弹性体中取出一个微小的正平行六面体，它在x和y方
12、向的尺寸分别为dx和dy，在z方向的尺寸为一个单位长度。,以x为投影轴，列出投影的平衡方程：,约简以后，两边除以dxdy，得：,同理，以y为投影轴，列出投影的平衡方程，化简得：,2018/9/27,几何方程,经过弹性体内的任意一点P，沿x轴和y轴的正方向取两个微小长度的线段PAdx和PBdy。假定弹性体受力后，P,A,B三点分别移动到P,A,B.,线段PA的线应变是：,注:由于位移微小，y方向的位移v引起的PA的伸缩，是高一阶微量，略去不计。,线段PB的线应变是：,线段PA与 PB之间的直角的改变，即切应变,线段PA的转角是：,线段PB的转角是：,2018/9/27,物理方程,在理想的弹性体中
13、，形变分量和应力分量之间的关系，在材料力学根据胡克定律导出如下：,在平面应力问题中，式变为：,在平面应变问题中，只要将上式中的E换为 ，换为 就得到平面应变问题的物理方程。,2018/9/27,假定已知任一点P处坐标面上的应力分量x， y ，x y = y x 。求经过该点的，平行于z轴而倾斜于x轴和 y轴的任何倾斜面上应力。 在P点附近取一个平面AB，它平行于上述斜面，并经过P点划出一个微小的三棱柱PAB。当AB无限小而趋于P点时，平面AB上的应力就成为斜面上的应力。,平面问题中一点的应力状态,设斜面AB 的长度为ds，则PB面及PA面的长度分别为 lds及mds，而PAB的面积为 ldsm
14、ds/2，棱柱的厚度设为1。 由x轴平衡条件，得：,其中，fx为体力分量。将上式除以ds，并令ds趋于0（斜面AB趋于P点），即得：,由y轴平衡条件，得：,用n表示斜面AB的外法线方向，其方向余弦为：,2018/9/27,边界条件,若在su部分边界上给定了约束位移分量 和 ，则对于此边界上的每一点，位移函数u和v应满足条件：,其中（u）s 和 （v）s 是位移的边界值， 和 在边界上是坐标的已知函数。,边界条件表示在边界上位移与约束，或应力与面力之间的关系式。它可以分为位移边界条件、应力边界条件和混合边界条件。,位移边界条件：,应力边界条件：,若在su部分边界上给定了面力 和 ，则由平衡条件得
15、出平面应力问题的应力(或面力）边界条件为：,其中，l,m是边界面外法线的方向余弦。,2018/9/27,圣维南原理,在求解弹性力学问题时，应力分量、形变分量和位移分量必须满足区域内的三套基本方程，还必须满足边界上的边界条件。但是，要使边界条件得到完全满足，往往遇到很大的困难。圣维南原理可为简化局部边界上的应力边界条件提供很大方便。圣维南原理表明，如果把物体的一小部分边界上的面力，变换为分布不同但静力等效的面力（主矢相同，对同一点的主矩也相同），那么，近处的应力分布将有显著的改变，但是远处所受的影响可以不计。,2018/9/27,圣维南原理的应用,例，设有柱形构件，在两端截面的形心受到大小相等而方向相反的拉力F（a）。如果把一端或两端的拉力变换为静力等效的力，则只有虚线划出的部分的应力分布有显著的改变，而其余部分所受影响是可以不计的。,由于（d）图中，面力连续分布，边界条件简单，应力容易求得。其它三种情况，应力难以求得。把d情况下的应力解答应用到其它三个情况，虽不能满足两端的应力边界条件，但仍然可以表明离杆端较远处的应力状态，没有显著的误差。 图e，构件右端有位移边界条件， ，d情况的解答，不能满足位移边界条件，但e图右端的面力，一定是合成为经过截面形心的力F。所以把图d情况的解答应用于图e时，仍然只是在靠近两端处有显著的误差，而在离两端较远之处，误差可以不计。,
