有限元分析的力学基础教学文案.ppt

有限元分析的力学基础本章要点o变形体的三大类基本变量o变形体的三大类基本方程及两类边界条件o弹性问题中的能量表示o平面应力、平面应变、刚体位移的特征及表达o应力及应变的分解2.1弹性力学同有限元分析的关系o弹性力学：弹性力学也称弹性理论，主要研究弹性体在外力作用或温度变化等外界因素下所产生的应力、应变和位移，从而解决结构或机械设计中所提出的强度和刚度问题。o是固体力学的重要分支，它研究弹性物体在外力和其它外界因素作用下产生的变形和内力。o研究对象：包括杆状构件在内的各种形状的弹性体。o弹性力学基本规律：变形连续规律、应力-应变关系和运动(或平衡)规律，它们有时被称为弹性力学三大基本规律。弹性力学中许多定理、公式和结论等，都可以从三大基本规律推导出来。2.1弹性力学同有限元分析的关系o弹性力学同材料力学的比较 1、研究内容：基本上没有什么区别。弹性力学也是研究弹性体在外力作用下的平衡和运动，以及由此产生的应力和变形。2、研究的对象：材料力学基本上只研究杆、梁、柱、轴等杆状构件，即长度远大于宽度和厚度的构件；弹性力学虽然也研究杆状构件，但还研究材料力学无法研究的板与壳及其它实体结构，即两个尺寸远大于第三个尺寸，或三个尺寸相当的构件。2.1弹性力学同有限元分析的关系o弹性力学同材料力学的比较 3、研究的方法：相同点：静力学、几何学与物理学三方面进行研究；不同点：材料力学：材料力学：对构件的整个截面建立分析方程，引用一些截面的变形状况或应力情况的假设，因而得出的结果往往是近似的，不精确。弹性力学：弹性力学：对构件采用无限小单元体来建立分析方程的，因而无须引用那些假设，分析的方法比较严密，得出的结论也比较精确。所以，可以用弹性力学的解答来估计材料力学解答的精确程度，并确定它们的适用范围。2.1弹性力学同有限元分析的关系 从几何形状复杂程度来考虑可以分为：1）简单形状变形体材料力学 2）任意形状变形体弹性力学 任意变形体是有限元方法有限元方法处理的对象，因而，弹性力学中有关变量和方程的描述是有限元方法的重要基础。弹性力学的弱点：弹性力学的弱点：由于研究对象的变形状态较复杂，处理的方法又较严谨，因而解算问题时，往往需要冗长的数学运算。但为了简化计算，便于数学处理，它仍然保留了材料力学中关于材料性质的假定。2.2 弹性力学中关于材料性质的假定o连续性：亦即物体整个体积内部被组成这种物体的介质填满，不留任何空隙。这样，物体内的一些物理量，如应力、应变、位 移等等才可以用座标的连续函数来表示。o完全弹性：亦即当使物体产生变形的外力被除去以后，物体能够完全恢复原形，而不留任何残余变形。这样，当温度不变时，物 体在任一瞬时的形状完全决定于它在这一瞬时所受的外力，与它过去的受力情况无关。服从虎克定律(应力应变成比例)o均匀性：也就是说整个物体是由同一种材料组成的。这样，整个物体的所有各部分才具有相同的物理性质，因而物体的弹性常 数(弹性模量和泊松系数)才不随位置座标而变。o各向同性：也就是说物体内每一点各个不同方向的物理性质和机械性质都是相同的。o物体的变形是微小的：亦即当物体受力以后，整个物体所有各点的位移都远小于物体的原有尺寸，因而应变和转角都远小于1，这样，在考虑物体变形以后的平衡状态时，可以用变形前的尺寸来代替变形后的尺寸，而不致有显著的误差；并且，在考虑物体的变形时，应变和转角的平方项或乘积项都可以略去不计，这就使得弹性力学中的微分方程都成为线性方程。2.2 弹性力学中关于材料性质的假定2.3弹性力学基本变量o基本变量2.3弹性力学基本变量o外力：指其他物体对研究对象（弹性体）的作用力。可以分为体积力和表面力 1、表面力：是分布于物体表面的力，如静水压力，一物体与另一物体之间的接触压力等。2、体力：是分布于物体体积内的外力，如重力、磁力、惯性力等。均为矢量。弹性体受外力以后，其内部将产生应力（内力）2.3弹性力学基本变量o内力：应力应力 -外力外力(或温度或温度)的作用的作用 内力内力 设作用于 上的内力为 ,则内力的平均集度,即平均应力,为 /这个极限矢量极限矢量S S,就是物体在截面mn上、P点的应力应力。应力就是弹性体内某一点作用于某截面单位面积上的内力2.3弹性力学基本变量oo正正应力应力 oo剪剪应力应力 每一个面上的应力分解为一个正应力和两个剪应力正应力下标表示作用在垂直于轴的面上同时也沿着轴方向作用的剪应力加上两个角码，前一个角码表明作用面垂直于哪一个坐标轴，后一个角码表明作用方向沿着哪一个坐标轴。2.3弹性力学基本变量o正面（外法线是沿着坐标轴的正方向）o负面（外法线是沿着坐标轴的负方向）o正面上的应力以沿坐标轴正方向为正，沿坐标轴负方向为负o负面上的应力以沿坐标轴负方向为正，沿坐标轴正方向为负正应力以拉应力为正，压应力为负2.3弹性力学基本变量o剪应力互等定律：作用在两个互相垂直的面上并且垂直于该两面交线的剪应力是互等的。(大小相等，正负号也相同)。因此剪应力记号的两个角码可以对调。不同的坐标表示应力张量一点的应力状态应变应变应变应变形状的改形状的改形状的改形状的改变变变变（形变形变形变形变）长度的改变和角度长度的改变和角度长度的改变和角度长度的改变和角度的改变的改变的改变的改变 应变和位移应变和位移 为了分析物体在其某一点 P 的形变状态,在这一点沿着坐标轴x,y,z 的正方向取三个微小的线段 PA,PB,PC。2.3弹性力学基本变量正正应变应变各线段的每单位长度的伸缩,即单位伸缩或相对伸缩。以伸长为正、缩短为负 剪剪应变应变各线段之间的直角的改变,用弧度表示。以直角减小为正、增大为负。2.3弹性力学基本变量位移位移位移位移就是位置的移就是位置的移就是位置的移就是位置的移动动动动。物体内任意一点的位移物体内任意一点的位移,用它在用它在x,y,z三三轴轴上的上的投影投影 ，来表示来表示以正以正标标向向为为正正。一般而论,弹性体内任意一点的体力分量、面力分量、应力分量、应变分量和位移分量,都是随着该点的位置而变的,因而都是位置坐标的函数位置坐标的函数位置坐标的函数位置坐标的函数。2.3弹性力学基本变量o位移与应变的关系2.3弹性力学基本变量应变应变位移位移刚体刚体位移位移位移位移刚体刚体转动转动strain-displacement relations.（几何方程（几何方程柯西方程）柯西方程）应力分量的矩阵表示称为应力列阵或应力向量。弹性体在载荷作用下，将产生位移和变形，即弹性体位置的移动和形状的改变。弹性体内任一点的位移可由沿直角坐标轴方向的3个位移分量 来表示。它的矩阵形式是：称作位移列阵或位移向量。o基本方程 受外部作用的任意形状变形体，在其微小体元dxdydz中，基于位移、应变和应力这三大类变量，可以建立以下三大类方程 平衡方程：外力和内力之间的平衡关系 几何方程：描述的是位移和应变之间关系 物理方程：应力和应变之间的关系2.4平面问题的基本力学方程o平衡方程：外力和内力之间的平衡关系o几何方程：描述的是位移和应变之间关系o物理方程：应力和应变之间的关系o边界条件：o平面(二维)平衡方程 平面问题的静力学平衡,设微小正六面体,在X,Y方向的尺寸dx,dy,Z方向的尺寸取一个单位长度.两个对面存在微小差量,通过中心点C,平行与Z轴的直线为轴,列出平衡方程上式两边除dxdy,可得:剪力互等关系以X轴为投影轴,满足平衡方程:上式两边除dxdy,可得:同理o平面(二维)几何方程 经过弹性体内任一点P,沿X轴和Y轴的方向取两个微小长度的线段PA=dx,PB=dy见图o变形协调条件 它的物理意义是：材料在变形过程中应该是整体连续的，不应该出现“撕裂”和“重叠”现象发生。写成矩阵形式为物理方程物理方程E称为杨氏模量反映材料对于拉伸或压缩变形的抵抗能力。是泊松系数，描写材料横向收缩或膨胀的特性。线应变（相对伸长或压缩）绝对伸长（或压缩）与原长之比称为相对伸长（或压缩）。公式：其中：设想直杆横截面是正方形每边长为 ,横向形变后为 。横向形变和纵向形变之比为泊松系数：当 时，为拉伸形变；时，为压缩形变，因而，它很好地反映形变程度。如直杆拉伸压缩时，还产生横向形变，则对应的应变(或形变)为：按照边界情况，弹性力学问题一般分为三类：位移边界问题：在边界面上全部给定位移，即全部是 Su 边界 应力边界问题：在边界面上全部给定表面力，即全部是应力边界。这时，外力（包括体力和面力）应是平衡力系。混合边界问题：既有Su 边界，又有应力边界。二者可以分别在边界表面不同的区域上，或同一区域不同的方向上。边界条件边界条件在 上弹性体的位移已知为 即有:用矩阵形式表示是弹性体V的全部边界为S,一部分边界上已知外力 称为力的边界条件，这部分边界用 表示;另一部分边界上弹性体的位移 已知，称为几何边界条件与位移边界条件，这部分边界用 表示。这两部分边界构成弹性体的全部边界，即:o几何边界条件作用在任意平面上该点的应力分量可以由下式表示为：其中 2.5空间问题的基本力学方程o平衡方程：外力和内力之间的平衡关系o几何方程：描述的是位移和应变之间关系o物理方程：应力和应变之间的关系o边界条件：平衡方程平衡方程oX方向负面oX方向正面oY方向负面oY方向正面oZ方向负面oZ方向正面oX方向力平衡o化简得oY方向力平衡o化简得oZ方向力平衡o化简得o如果这六个量在某点是已知的，就可以求得经过该点的任何面上的正应力和剪应力，因此，这六个量可以完全确定该点的应力状态，它们就称为在该点的应力分量。o一般说来，弹性体内各点的应力状态都不相同，因此，描述弹性体内应力状态的上述六个应力分量并不是常量，而是坐标x、y、z的函数。o六个应力分量的总体，可以用一个列矩阵来表示：几何方程几何方程工程应变o写成矩阵形式为o几何方程可见，当弹性体的位移分量完全确定时，应变分量是完全确定的。反过来，当应变分量完全确定时，位移分量却不完全确定；这是因为，具有确定形状的物体，可能发生不同的刚体位移。为了说明这一点，试命：式中的u0,v0,w0,x,y,z是积分常数。u0弹性体沿x方向的刚体移动v0 弹性体沿y方向的刚体移动 w0 弹性体沿z方向的刚体移动x 弹性体绕x轴的刚体转动y 弹性体绕y轴的刚体转动z 弹性体绕z轴的刚体转动为了完全确定弹性体的位移，必须有六个适当的约束条件来确定这六个刚体位移。o变形协调条件 当6个应变分量满足以上应变协调方程时，就能保证得到单值连续的位移函数。当沿x轴方向的两个对面受有均匀分布的正应力时，在满足先前假定的材料性质条件下，正应力不会引起角度的任何改变，而其在x方向的单位伸长则可表以方程 弹性体在x方向的伸长还伴随有侧向收缩，即在y和z方向的单位缩短可表示为：方程既可用于简单拉伸，也可用于简单压缩，且在弹性极限之内，两种情况下的弹性模量和波桑系数相同。应力分量与应变分量之间的关系-虎克定律物理方程物理方程设图中的弹性体在各面上都受有均匀分布的正应力，则合成应变的分量可用前面两式求得。实验证明，只须将三个应力中的每一应力所引起的应变分量叠加，就得到合成应变的分量。单位伸长与应力之间的关系完全由两个物理常数E及所确定。两个常数也可用来确定剪应力与剪应变之间的关系。公式的适用范围公式的适用范围 :在线弹性范围内,小 变形条件下,各向同性材料。如果弹性体的各面有剪应力作用任何两坐标轴的夹角的改变仅与平行于这两轴的剪应力分量有关，即得到：正应变与剪应变是各自独立的。因此，由三个正应力分量与三个剪应力分量引起的一般情形的应变，可用叠加法求得；即将六个关系式写在一起，得弹性方程或物理方程，这种空间状态的应力应变关系称为广义虎克定律。写成矩阵形式为 边界条件边界条件XN,YN,ZN分别为作用在某一任意平面上的沿三个坐标轴方向的分量。对于已知应力边界条件的情况，相应的应力边界条件为 o二维问题：2个位移分量，3个应力分量，3个应变分量 2个平衡方程，3个几何方程，3个物理方程o三维问题：3个位移分量，6个应力分量，6个应变分量 3个平衡方程，6个几何方程，6个物理方程 我们得到的变量和方程都是从任意变形体中所取出来的微单元体来建立的，因此无论对象的几何形状和边界条件如何不同，其基本变量和基本方程是完全相同，不同之处在于边界条件，所以求解的难度是如何处理边界条件（几何形状）。2.5弹性问题中的能量表示o能量分类1）施加外力在可能位移上所作的功(即外力在弹性变形过程中所做的功)。2）变形体由于变形而存储的能量(即由于变形而储存于弹性体内的能量)。o主要是研究泛函及其极值的求解方法o泛函：就是以函数为自变量的一类函数，简单地讲，泛函就是函数的函数。o弹性力学变分法中所研究的泛函，就是弹性体的能量，如形变势能、外力势能等。因此，弹性力学中的变分法又称为能量法。o取位移为基本未知函数2.5弹性问题中的能量表示o2.5.1外力功 施加外力在可能位移上所作的功，外力有两种，包括作用在物体上的面力和体力，这些力被假设为与变形无关的不变力系（保守力），则外力功包括这两部分力在可能位移上所作的功。2.5弹性问题中的能量表示o2.5.2应变能 以位移（或应变）为基本变量所表达的变形能叫做应变能（strain energy）。它也包括两部分 1）对应于正应力与正应变的应变能 2）对应于剪应力和剪应变的应变能2.5弹性问题中的能量表示1.1.单向拉伸杆单向拉伸杆单向拉伸杆单向拉伸杆外力做功外力做功外力做功外力做功弹性体应变能弹性体应变能弹性体应变能弹性体应变能单位体积应变能单位体积应变能应变能密度应变能密度应变能密度应变能密度静加载是线性静加载是线性的，没有动能的，没有动能与热能的变化与热能的变化2.5弹性问题中的能量表示对应于正应力与正应变的应变能，另外两个方向上的计算类似。对应于剪应力和剪应变的应变能（其它两个剪应力类似）2.受均匀剪应力时受均匀剪应力时受均匀剪应力时受均匀剪应力时应变能密度应变能密度应变能密度应变能密度3.受复杂应力状态受复杂应力状态受复杂应力状态受复杂应力状态最终弹性应变能与最终弹性应变能与变形过程变形过程无关，只取决于无关，只取决于变形的最终状态。变形的最终状态。2.5弹性问题中的能量表示在平面问题中在平面问题中,。在平面应力问题中还有在平面应力问题中还有 ;在平面应变问题中在平面应变问题中 ,还有还有 。因此。因此,在两种平在两种平面问题中面问题中,弹性体的形变势能密度的表达式都简化为弹性体的形变势能密度的表达式都简化为 应变能密度应变能密度应变能密度应变能密度2.5弹性问题中的能量表示 在一般的平面问题中在一般的平面问题中,弹性体各部分的受力并非均匀弹性体各部分的受力并非均匀,各个应力分量和形变分量都是坐标各个应力分量和形变分量都是坐标 x x 和和 y y 的函数的函数,因而因而形变势能密度形变势能密度U U l l 一般也是坐标一般也是坐标 x x 和和 y y 的函数。的函数。为了得为了得出整个弹性体具有的形变势能出整个弹性体具有的形变势能 U U,必须将形变势能密度必须将形变势能密度 u 在整个弹性体内积分。和以前一样在整个弹性体内积分。和以前一样,为了简便为了简便,在在 z z 方向取一个单位长度。这样就得到方向取一个单位长度。这样就得到(在平面区域在平面区域 A A 内内 )2.5弹性问题中的能量表示2.5弹性问题中的能量表示空间问题的能力密度考虑初始应力及应变o从而得系统势能*负号表示外力势能为负值图图表表示示一一平平衡衡的的杠杠杆杆，对对C点点写写力力矩平衡方程：矩平衡方程：图图表表示示杠杠杆杆绕绕支支点点C转转动动时时的的刚刚体体位移图：位移图：综合可得：综合可得：即：即：上上式式是是以以功功的的形形式式表表述述的的。表表明明：图图a的的平平衡衡力力系系在在图图b的的位位移移上上作作功功时时，功功的的总总和和必必须须等等于于零零。这这就叫做虚功原理。就叫做虚功原理。2.5弹性问题中的能量表示虚功原理虚功原理 进进一一步步分分析析。当当杠杠杆杆处处于于平平衡衡状状态态时时，和和 这这两两个个位位移移是是不不存存在在的的，但但是是如如果果某某种种原原因因，例例如如人人为为地地振振一一下下让让它它倾倾斜斜，一一定定满满足足的关系。的关系。将将这这个个客客观观存存在在的的关关系系抽抽象象成成一一个个普普遍遍的的原原理理，去去指指导导分分析析和和计计算算结构。结构。对对于于在在力力的的作作用用下下处处于于平平衡衡状状态态的的任任何何物物体体，不不用用考考虑虑它它是是否否真真正正发发生生了了位位移移，而而假假想想它它发发生生了了位位移移，(由由于于是是假假想想，故故称称为为虚虚位位移移)，那那么么，物物体体上上所所有有的的力力在在这这个个虚虚位位移移上上的的总总功功必必定定等等于于零零。这这就就叫叫做做虚虚位位移移原原理理，也也称称虚虚功功原原理理。在在图图a中中的的 和和 所所作作的的功功就就不不是是发发生生在在它它本本身身(状状态态a)的的位位移移上上，(因因为为它它本本身身是是平平衡衡的的，不不存存在在位位移移)，而而是是在在状状态态(b)的的位位移移上上作作的的功功。可可见见，这这个个位位移移对对于于状状态态(a)来来说就是虚位移，亦即是状态说就是虚位移，亦即是状态(a)假象的位移。假象的位移。虚功原理应用范围虚功原理应用范围 必必须须指指出出，虚虚功功原原理理的的应应用用范范围围是是有有条条件件的的，它它所所涉涉及及到到的的两两个个方方面面，力力和和位位移移并并不不是是随随意意的的。对对于于力力来来讲讲，它它必必须须是是在在位位移移过过程程中中处处于于平平衡衡的的力力系系；对对于于位位移移来来讲讲，虽虽然然是是虚虚位位移移，但但并并不不是是可可以以任任意意发发生的。它必须是和约束条件相符合的微小的刚体位移。生的。它必须是和约束条件相符合的微小的刚体位移。还还要要注注意意，当当位位移移是是在在某某个个约约束束条条件件下下发发生生时时，则则在在该该约约束束力力方方向向的的位位移移应应为为零零，因因而而该该约约束束力力所所作作的的虚虚功功也也应应为为零零。这这时时该该约约束束力力叫叫做做被被动动力力。(如如图图中中的的反反力力 由由于于支支点点C没没有有位位移移，故故 所所作作的的虚虚功功对对于于零零)。反反之之，如如图图中中的的 和和 是是在在位位移移过过程程中中作作功功的的力力，称称为为主主动动力力。因因此此，在在平平衡衡力力系系中中应应当当分分清清楚楚哪哪些些是是主主动动力力，哪哪些些是是被被动动力力，而而在在写写虚虚功功方方程程时时，只只有有主主动动力力作作虚虚功功，而而被被动动力力是是不不作作虚虚功功的。的。虚功原理与虚功方程虚功原理与虚功方程虚功原理表述如下：虚功原理表述如下：在在力力的的作作用用下下处处于于平平衡衡状状态态的的体体系系，当当发发生生与与约约束束条条件件相相符符合合的的任任意意微微小小的的刚刚体体位位移移时时，体体系系上上所所有有的的主主动动力力在在位位移移上上所所作作的的总总功功(各各力力所所作作的的功功的的代代数数和和)恒恒对对于于零。零。虚功原理用公式表示为：虚功原理用公式表示为：这就是虚功方程，其中这就是虚功方程，其中P和和 相应的代表力和虚位移。相应的代表力和虚位移。虚功原理虚功原理-用于弹性体的情况用于弹性体的情况 虚虚功功方方程程是是按按刚刚体体的的情情况况得得出出的的，即即假假设设图图的的杠杠杆杆是是绝绝对对刚刚性性，没没有任何的变形，因而在方程中没有有任何的变形，因而在方程中没有内功项内功项出现，而只有外功项。出现，而只有外功项。将将虚虚功功原原理理用用于于弹弹性性变变形形时时，总总功功W要要包包括括外外力力功功(T)和和内内力力功功(U)两两部部分分，即即：W=T-U ；内内力力功功(-U)前前面面有有一一负负号号，是是由由于于弹弹性性体体在在变变形形过过程程中中，内内力力是是克克服服变变形形而而产产生生的的，所所有有内内力力的的方方向向总是与变形的方向相反，所以内力功取负值。总是与变形的方向相反，所以内力功取负值。根据虚功原理，总功等于零得：根据虚功原理，总功等于零得：T -U=0 外力虚功外力虚功 T =内力虚功内力虚功 U 弹弹性性力力学学中中的的虚虚功功原原理理可可表表达达为为：在在外外力力作作用用下下处处于于平平衡衡状状态态的的弹弹性性体体，如如果果发发生生了了虚虚位位移移，那那么么所所有有的的外外力力在在虚虚位位移移上上的的虚虚功功(外外力力功功)等于整个弹性体内应力在虚应变上的虚功等于整个弹性体内应力在虚应变上的虚功(内力功内力功)。虚功原理虚功原理-用于弹性体的情况用于弹性体的情况虚应变能虚应变分量外力虚功内力虚功即应力在虚应变上做的的虚功，也称虚应变能外力虚功即作用于弹性体上的外力在虚位移上做的功由于虚位移是微小的，可认为在虚位移发生过程中外力保持为常量，则上式的变分符号可提到积分号外。最小势能原理 在有限元的理论中，最小势能原理是在所有满足给定边界条件的位移时，满足平衡微分方程的位移使得势能取得最小值。最小势能原理就是说当一个体系的势能最小时，系统会处于稳定平衡状态。或者说在所有几何可能位移中，真实位移使得总势能取最小值 最小势能原理：表明在满足位移边界条件的所有可能位移中，实际发生的位移使弹性体的势能最小。即对于稳定平衡状态，实际发生的位移使弹性体总势能取极小值。显然，最小势能原理与虚功原理完全等价。虚功原理的矩阵表示虚功原理的矩阵表示i点外力分量点外力分量j点外力分量点外力分量外外力力分分量量用用 表表示示；引起的应力分量用引起的应力分量用 表示表示虚功原理的矩阵表示虚功原理的矩阵表示假设发生了虚位移假设发生了虚位移虚位移分量为虚位移分量为用用 表表示示；引引起起的的虚虚应应变分量用变分量用 表示表示虚功原理的矩阵表示虚功原理的矩阵表示 在虚位移发生时，外力在虚位移上的虚功是：在虚位移发生时，外力在虚位移上的虚功是：式中式中 是是 的转置矩阵。的转置矩阵。同同样样，在在虚虚位位移移发发生生时时，在在弹弹性性体体单单位位体体积积内内，应应力力在在虚虚应应变变上上的的虚虚功是：功是：因此，在整个弹性体内，应力在虚应变上的虚功是：因此，在整个弹性体内，应力在虚应变上的虚功是：根据虚功原理得到：根据虚功原理得到：这这就就是是弹弹性性变变形形体体的的虚虚功功方方程程，它它通通过过虚虚位位移移和和虚虚应应变变表表明明外外力力与与应应力力之间的关系。之间的关系。这是以后推导有限元方程的基础。这是以后推导有限元方程的基础。2.6应用实例应用实例A1,l1A2,l2R31.离散化1232.位移函数A,lijF i uiF j ujxx2.位移函数A,lijF i uiF j uj（单元内位移线性分布）（单元内位移线性分布）形函数矩阵形函数矩阵 2.6应用实例应用实例763.单元刚度矩阵方程A 虚功原理 外力虚功虚应变能应变应力应变矩阵弹性矩阵2.6应用实例应用实例单元刚度矩阵单元的刚度方程单元刚度矩阵ElementElement2.6应用实例应用实例78B 最小势能定理 外力虚功虚应变能2.6应用实例应用实例由势能变分原理（势能最小原理）得势能变分，整理得平衡方程2.6应用实例应用实例4 整体分析整体分析 整体分析就是建立整个离散结构所有节点位移与整体分析就是建立整个离散结构所有节点位移与外力之间的关系，实现未知节点位移的求解外力之间的关系，实现未知节点位移的求解1）整体平衡方程）整体平衡方程 整体刚度方程可基于势能变分原理建立，也可根据节点的静力平衡来实现（即每个节点静力平衡）。节点i的平衡为三个节点三个自由度，即（A A）扩充单元刚度方程法）扩充单元刚度方程法位移协调性位移协调性载荷的叠加性载荷的叠加性A1,l1A2,l2R31232.6应用实例应用实例整体刚度方程2.6应用实例应用实例A1,l1A2,l2P123（B B）“对号入座对号入座”法法 （方便编程）（方便编程）i j 1 2 1 2 2 3 2 3Total1 2 31 2 31 1 2 2 3 32.6应用实例应用实例5.引入边界条件求解A1,l1A2,l2R3123边界条件支反力2.6应用实例应用实例o结构离散结构离散o单元分析单元分析o整体分析整体分析2.7基本步骤为三大步骤基本步骤为三大步骤1、结构离散:就是用假想的线或面将连续物体分割成有限个单元组成的集合体且单元之间仅在节点处连接，单元之间的作用仅由节点传递。（基本要求）注意的问题o单元：满足一定几何特性和物理特性的最小结构域o节点：单元与单元间的连接点o节点力：单元与单元间通过节点的相互作用力o节点载荷：作用于节点上的外载2.7基本步骤为三大步骤基本步骤为三大步骤2单元分析单元分析:1），选择插值（位移）函数；2），构造位移函数。插值函数：用以表示单元内物理量变化（如位移或位移场）的近似函数。由于该近似函数常由单元节点物理量值插值构成，故称为插值函数，如单元内物理量为位移，则该函数称为位移函数。选择位移函数的一般原则o位移函数在单元节点的值应等于节点位移（即单元内部是连续的）；o所选位移函数必须保证有限元的解收敛于真实解。2.7基本步骤为三大步骤基本步骤为三大步骤位移函数一般采用多项式形式，在单元内选适当阶次的多项式可得到与真实解接近的近似解构造位移函数：如平面问题位移函数的一般形式为1.多项式项数越多，则逼近真实位移的精度越高，项数的多少由单元的自由度数决定。2多项式选取应由低阶到高阶，尽量选择完全多项式以提高单元精度。3.选取多项式时，还应使所选取的多项式具有坐标的对称性，即按Pascal（帕斯卡）三角形来选择2.7基本步骤为三大步骤基本步骤为三大步骤位移函数构造方法：1.广义坐标法：2插值函数法：即将位移函数表示为各个节点位移与已知插值基函数积的和2.7基本步骤为三大步骤基本步骤为三大步骤3）单元特性分析单元特性分析：单元特性分析的基本任务就是建立单元的平衡方程，也称为刚度方程。在选择了单元类型和相应的位移函数后，即可按弹性力学的几何方程、物理方程导出单元应变与应力的表达式，最后利用虚位移原理或最小势能原理或直接法或加权残值法建立单元的平衡方程，即单元节点力与节点位移间的关系。2.7基本步骤为三大步骤基本步骤为三大步骤2.7基本步骤为三大步骤基本步骤为三大步骤2.7基本步骤为三大步骤基本步骤为三大步骤3.整体分析整体分析的基本任务包括建立整体平衡方程，引入边界条件，完成整体方程求解。整体平衡方程的建立有多种方法，可基于能量原理（势能变分或虚位移原理）推导，也可基于节点力平衡得到。在引入边界条件之前，整体平衡方程是奇异的，这意味着整体方程是不可解的。方程求解包括边界条件引入和数值计算，一旦利用适当的数值方法求出未知的节点位移，则可按前述的应力应变公式计算出各个单元的应变、应力等物理量o刚度由使其产生单位变形所需的外力值来量度，刚度是指零件在载荷作用下抵抗弹性变形的能力。o单元的刚度矩阵：单元刚度矩阵反应的是单元节点力与单元节点位移的关系；o总刚度矩阵反应的是整体的节点力与节点位移的关系；刚度矩阵将总体坐标下的节点位移与整个结构的总体力联系在一起。补充实例（刚度矩阵的理解）o单元刚度矩阵行数等于位移向量的分量个数，列数等于为位移的列向量的分量个数，由于两者相等所以单刚是个方阵。o结构的总体刚度矩阵即结构的原始刚度矩阵，每1个元素的物理意义就是当其所在列对应的节点位移分量等于单位位移（其余结点位移分量为0）时，其所在行对应的节点力的数值。表示由于第j个自由度的单位位移dj在第i个自由度需要的力补充实例（刚度矩阵的理解）3.2.5 单元的刚度矩阵的性质单元的刚度矩阵的性质 a.单元刚度矩阵仅与单元的几何特征和材料性质有关。仅与单元的横截面积A、惯性矩I、单元长度l、单元的弹性模量E有关。b.单元刚度矩阵是一个对称阵。在单元刚度矩阵对角线两侧对称位置上的两个元素数值相等，即，根据是反力互等定理。c.单元刚度矩阵是一个奇异阵。d.单元刚度矩阵可以分块矩阵的形式表示。具有确定的物理意义。整体刚度矩阵的性质整体刚度矩阵的性质 整体刚度矩阵 中位于主对角线上的子块 ，称为主子块，其余 为副子块。a.中主子块 由结点i的各相关单元的主子块扩展之后叠加求得，即 b.当结点i、j为单元e的相关结点时，中副子块 为该单元e相应的副子块，即 。c.当结点i、j为非相关结点时，中副子块 为零子块，即 。d.仅与各单元的几何特性、材料特性，即A、I、l、E等因素有关。e.为对称方阵，f.为奇异矩阵，其逆矩阵不存在，因为建立整体刚度矩阵时没有考虑结构的边界约束条件。g.为稀疏矩阵，整体刚度矩阵中的非零元素分布区域的宽度与结点编号有关，非零元素分布在以对角线为中心的带状区域内，称为带状分布规律，见图a。在包括对角线元素在内的区域中，每行所具有的元素个数叫做把半带宽，以d表示。最大半带宽等于相邻结点号的最大差值加 1 与结点自由度数的乘积，结点号差越大半带宽也就越大。计算机以半带宽方式存储，见图b。半带宽越窄，计算机的存储量就越少，而且可以大幅度减少求解方程所需的运算次数。其效果对大型结构显得尤为突出。图 整体刚度矩阵存储方法 h.整体刚度矩阵稀疏阵。故整体刚度矩阵不能求逆，必须作约束处理方能正确地将结点位移求出，进而求出结构的应力场。(a)带状分布规律(b)带状存储 约束处理及求解约束处理及求解 约束处理的必要性约束处理的必要性 建立结构原始平衡方程式 时，并未考虑支承条件（约束），也就是说，将原始结构处理成一个自由悬空的、存在刚体位移的几何可变结构。整体刚度矩阵是奇异矩阵，因此，无法求解。可以参照第 2 章的原则，结合实际工程结构引入支承条件，即对结构原始平衡方程式 做约束处理。约束处理后的方程称为基本平衡方程。统一记为 约束处理方法约束处理方法 约束处理常用方法有填0置1法和乘大数法。采用这两种方法不会破坏整体刚度矩阵的对称性、稀疏性及带状分布等特性。1 填0置1法 将整体刚度矩阵中与之相对应的主对角元素全部置换成1，相应行和列上的其它元素均改为0。同时，所在同一行上的载荷分量替换成0，则有1 填0置1法 也可简便地采用划行划列的办法。在整体刚度矩阵中将与约束位移为 0 的行和列划掉，包括相关的所在行的位移和载荷向量。2 乘大数法将整体刚度矩阵中与之相对应的主对角元素全部乘以一个大数N，一般取 。同时，将相应同一行上的载荷分量替换成 N 乘以其主对角刚度系数和给定的强迫位移（包括零位移）。补充实例（势能原理的应用）3根弹簧弹性系数分别为k1、k2、k3,在弹簧2末端施加一个F的力，利用势能原理推导系统的刚度矩阵。第三章 几个特殊问题o实际问题中，经常有一些比较典型的情况，需要有针对性的进行处理，如o厚度较薄的平面问题o厚度较厚的等截面平面应变问题o弹性力学可分为空间问题和平面问题，严格地说，任何一个弹性体都是空间物体，一般的外力都是空间力系，因而任何实际问题都是空间问题，都必须考虑所有的位移分量、应变分量和应力分量。但是，如果所考虑的弹性体具有特殊的形状，并且承受的是特殊外力，就有可能把空间问题简化为近似的平面问题，只考虑部分的位移分量、应变分量和应力分量即可。o平面应力问题o平面应变问题3.1平面应力问题o几何形状特征：物体在一个坐标方向（例如z）的几何尺寸远远小于其他两个坐标方向的几何尺寸，如图所示的薄板。o载荷特征：在薄板的两个侧表面上无表面载荷，作用于边缘的表面力平行于板面，且沿厚度不发生变化，或虽沿厚度变化但对称于板的中间平面，体积力亦平行于板面且沿厚度不变。几何方程物理方程平衡方程考察应力边界特点注意：注意：3.2平面应变问题o几何形状特征：物体沿一个坐标轴（例如z轴）方向的长度很长，且所有垂直于z轴的横截面都相同，即为一等直柱体；位移约束条件或支承条件沿z方向也相同。o 载荷特征：柱体侧表面承受的表面力以及体积力均垂直于z轴，且分布规律不随z变化。几何方程物理方程平衡方程特点：方向位移以及应变为零（但z方向的应力不为零）注意注意：由于z方向的伸缩被阻止，所以可以看出,平面应力和平面应变问题的物理方程可以通过以下变换互相得到因为在平面应变问题中也有和，所以有相应的剪应变为零。平面应力平面应变平面应变平面应力平面应力和平面应变的区别平平 面面 应应 力力平平 面面 应应 变变基本量基本量平衡方程平衡方程相 同几何方程几何方程相 同物理方程物理方程形式相同、内容不同。将代入平面应力问题的物理方程即可得到平面应变问题的物理方程。此课件下载可自行编辑修改，仅供参考！此课件下载可自行编辑修改，仅供参考！感谢您的支持，我们努力做得更好！谢谢感谢您的支持，我们努力做得更好！谢谢