有限元方法与ANSYS应用二讲课件.ppt

1、l弹性力学中的五个基本假定弹性力学中的五个基本假定。关于材料性质的假定及其在建立弹力理论中的作用：关于材料性质的假定及其在建立弹力理论中的作用：（1）连续性连续性 假定物体是连续的。各物理量可用连续函数表示。假定物体是连续的。各物理量可用连续函数表示。（2）完全弹性完全弹性 假定物体是假定物体是,a.完全弹性完全弹性外力取消外力取消,变形恢复变形恢复,无残余变形。无残余变形。b.线性弹性线性弹性应力与应变成正比。应力与应变成正比。即应力与应变关系可用胡克定律表示（物理线性）。即应力与应变关系可用胡克定律表示（物理线性）。（3）均匀性均匀性 假定物体由同种材料组成。假定物体由同种材料组成。E、等

2、与位置等与位置 无关无关（4）各向同性各向同性 假定物体各向同性。假定物体各向同性。E、等与方向无关。等与方向无关。符合（符合（1）-（4）假定的称为理想弹性体。）假定的称为理想弹性体。（5）小变形假定小变形假定 假定位移和形变为很小。假定位移和形变为很小。a.位移物体尺寸位移物体尺寸,例：梁的挠度例：梁的挠度v梁高梁高h.a.简化平衡条件：考虑微分体的平衡条件时，可以用变形前简化平衡条件：考虑微分体的平衡条件时，可以用变形前的尺寸代替变形后的尺寸。的尺寸代替变形后的尺寸。b.简化几何方程：在几何方程中，由于简化几何方程：在几何方程中，由于 可略去可略去 等项等项,使几何方程成为线性方程。使几

3、何方程成为线性方程。.1 ,.b,),(),(),(32 2),(弹性力学中几个基本概念外力外力其他物体对研究对象（弹性体）的作用力其他物体对研究对象（弹性体）的作用力 体力体力作用于物体体积内的力。如重力、惯性力和电磁力。以作用于物体体积内的力。如重力、惯性力和电磁力。以单位体积内所受的力来量度，单位体积内所受的力来量度，面力面力作用于物体表面上的力。如流体压力和接触力。以单位作用于物体表面上的力。如流体压力和接触力。以单位面积所受的力来量度面积所受的力来量度 集中力集中力如牵引力如牵引力 弹性体受外力以后，其内部将产生应力。弹性体受外力以后，其内部将产生应力。.,zyxfff.,zyxff

4、fv 0FVlimf222xyzffffxfyfzf22L MToxyzfyfzPVfFfxfF xozyS 0FSlimf222xyzffffxfyfzf12L MT PSFzfyfxffFf内力内力假想切开物体，截面两边互相作用的力（合力和合力矩，假想切开物体，截面两边互相作用的力（合力和合力矩，称为内力。称为内力。应力应力截面上某一点处，单位截面面积上的内力值。截面上某一点处，单位截面面积上的内力值。正应力是作用在垂直于正应力是作用在垂直于x轴的面上同时也沿着轴的面上同时也沿着X轴轴方向作用的。方向作用的。剪应力是作用在垂直于剪应力是作用在垂直于X轴的面上而沿着轴的面上而沿着y轴方向轴方

5、向作用的。作用的。应力成对出现，坐标面上的应力以正面正向，负面负向应力成对出现，坐标面上的应力以正面正向，负面负向为正。为正。(如果某一个面上的外法线是沿着坐标轴的正方向，这个如果某一个面上的外法线是沿着坐标轴的正方向，这个面上的应力就以沿坐标轴正方向为正，沿坐标轴负方向为负。面上的应力就以沿坐标轴正方向为正，沿坐标轴负方向为负。相反，如果某一个面上的外法线是沿着坐标轴的负方向，这个相反，如果某一个面上的外法线是沿着坐标轴的负方向，这个面上的应力就以沿坐标轴的负方向为正，沿坐标轴正方向为负面上的应力就以沿坐标轴的负方向为正，沿坐标轴正方向为负。)xxyoxyzPAFpIIImnA 0FAlim

6、ppF12L MT应力分量的空间方位，注意每个分量的作用面与作用方向！剪应力互等定律剪应力互等定律 作用在两个互相垂直的面上并且垂直于该两面交线的剪应力是作用在两个互相垂直的面上并且垂直于该两面交线的剪应力是互等的。互等的。(大小相等，正负号也相同大小相等，正负号也相同)。因此剪应力记号的两个角。因此剪应力记号的两个角码可以对调。码可以对调。应力分量应力分量一般说来，弹性体内各点的应力状态都不相同，因此，一般说来，弹性体内各点的应力状态都不相同，因此，描述弹性体内应力状态的上述六个应力分量并不是常描述弹性体内应力状态的上述六个应力分量并不是常量，而是坐标量，而是坐标x、y、z的函数。的函数。六

7、个应力分量的总体，可以用一个列矩阵六个应力分量的总体，可以用一个列矩阵 来表示：来表示：位移、应变、刚体位移 弹性体在受外力以后，还将发生变形。物体的变形状态，一般有两种方式来描述：1、给出各点的位移；2、给出各微元体的变形 弹性体内任一点的位移，用此位移在x、y、z三个坐标轴上的投影u、v、w来表示。这三个投影称为位移分量。以沿坐标轴正方向为正，沿坐标轴负方向为负.一般情况下，弹性体受力以后，各点的位移并不是定值，而是坐标的函数。微元体的变形可以分为两类：1、长度的变化，2、角度的变化。任一线素的长度的变化与原有长度的比值称为线应变(或称正应变)，用符号 来表示。沿坐标轴的线应变，则加上相应

8、的角码，分别用 来表示。当线素伸长时，其线应变为正。反之，线素缩短时，其线应变为负。这与正应力的正负号规定相对应。任意两个原来彼此正交的线素，在变形后其夹角的变化值称为角应变或剪应变，用符号 来表示。两坐标轴之间的角应变，则加上相应的角码，分别用 来表示。规定当夹角变小时为正，变大时为负，与剪应力的正负号规定相对应(正的 引起正的 ，等等)。描述变形体的基本方程描述变形体的基本方程 弹性力学基本方程 考虑微元体各个面上的考虑微元体各个面上的法向应力和剪应力法向应力和剪应力与其与其体力平衡体力平衡，注意，注意应力从一个面到对面是变化的，即有增量，将作用于微元体各个应力从一个面到对面是变化的，即有

9、增量，将作用于微元体各个方向的力求和，略去高阶项，可得平衡方程（受力状态的描述）方向的力求和，略去高阶项，可得平衡方程（受力状态的描述）：平衡微分方程的推导过程平衡微分方程的推导过程二、几何方程二、几何方程空间独立的应变分量空间独立的应变分量正应变分量：切应变分量：应变分量列阵：如何计算正应变和切应变？正应变切应变如何描述位移分量？如何描述位移分量？几何方程几何方程几何方程描述几何方程描述应变分量应变分量与与位移分量位移分量之间的关系！之间的关系！如何推导几何方程？如何推导几何方程？几何坐标：几何位移：因此，正应变切应变切应变采用同样的思路，分析微单元体在xoz与yoz上的投影及之后的几何变形

10、，可获得其他三个分量（1个正应变分量和2个切应变分量)。几何方程应变分量 三、物理方程三、物理方程对均匀、连续、各向同性的弹性材料，存在虎克定律对均匀、连续、各向同性的弹性材料，存在虎克定律l 未知数未知数 位移位移3个个+应力应力6个个+应变应变6个个=15l 方程个数方程个数 平衡方程平衡方程3个个+几何方程几何方程6个个+物理方程物理方程6个个=15原则上原则上15个方程可以求解个方程可以求解15个未知量，但实际上先求出一部分个未知量，但实际上先求出一部分，再通过方程求剩下的。再通过方程求剩下的。四、边界条件四、边界条件4.1 位移边界条件给出边界面上所有点的位移给定边界上位移4.2 应

11、力边界条件应力边界条件给出边界面给出边界面abc上表面力上表面力!目前有限元法采用位移法，以三个位移分量为基目前有限元法采用位移法，以三个位移分量为基本未知量。本未知量。位移产生变形，变形产生应变，应变产生应力，位移产生变形，变形产生应变，应变产生应力，应力于外力平衡。应力于外力平衡。协调的含义：不重叠、不开裂！图a示一平衡的杠杆，对C点写力矩平衡方程：图b表示杠杆绕支点C转动时的刚体位移图：综合可得：即：式是以功的形式表述的。表明：图a的平衡力系在图b的位移上作功时，功的总和必须等于零。这就叫做虚功原理。五、虚功原理五、虚功原理进一步分析。当杠杆处于平衡状态时，进一步分析。当杠杆处于平衡状态

12、时，和和 这两个位移是这两个位移是不存在的，但是如果某种原因，例如人为地振一下让它倾斜，一不存在的，但是如果某种原因，例如人为地振一下让它倾斜，一定满足上式的关系。将这个客观存在的关系抽象成一个普遍的原定满足上式的关系。将这个客观存在的关系抽象成一个普遍的原理，去指导分析和计算结构。理，去指导分析和计算结构。对于在力的作用下处于平衡状态的任何物体，不用考虑它是否真对于在力的作用下处于平衡状态的任何物体，不用考虑它是否真正发生了位移，而假想它发生了位移，正发生了位移，而假想它发生了位移，(由于是假想，故称为虚位由于是假想，故称为虚位移移)，那么，物体上所有的力在这个虚位移上的总功必定等于零。，那

13、么，物体上所有的力在这个虚位移上的总功必定等于零。这就叫做这就叫做虚位移原理虚位移原理，也称，也称虚功原理虚功原理。在图。在图a中的中的 和和 所作所作的功就不是发生在它本身的功就不是发生在它本身(状态状态a)的位移上，的位移上，(因为它本身是平衡因为它本身是平衡的，不存在位移的，不存在位移)，而是在状态，而是在状态(b)的位移上作的功。可见，这个的位移上作的功。可见，这个位移对于状态位移对于状态(a)来说就是虚位移，亦即是状态来说就是虚位移，亦即是状态(a)假象的位移。假象的位移。必须指出，虚功原理的应用范围是有条件的，它所涉及必须指出，虚功原理的应用范围是有条件的，它所涉及到的两个方面，力

14、和位移并不是随意的。对于到的两个方面，力和位移并不是随意的。对于力力来讲，它必须是来讲，它必须是在位移过程中处于在位移过程中处于平衡的力系平衡的力系；对于；对于位移位移来讲，虽然是虚位移，来讲，虽然是虚位移，但并不是可以任意发生的。它必须是和约束条件相符合的微小的但并不是可以任意发生的。它必须是和约束条件相符合的微小的刚体位移。刚体位移。还要注意，当位移是在某个约束条件下发生时，则在该还要注意，当位移是在某个约束条件下发生时，则在该约束力方向的位移应为零，因而该约束力所作的虚功也应为零。约束力方向的位移应为零，因而该约束力所作的虚功也应为零。这时该约束力叫做这时该约束力叫做被动力被动力。(如上

15、图中的反力如上图中的反力 ，由于支点，由于支点C没没有位移，故有位移，故 所作的虚功等于零所作的虚功等于零)。反之，如上图中的。反之，如上图中的 和和 是是在位移过程中作功的力，称为主动力。因此，在平衡力系中应当在位移过程中作功的力，称为主动力。因此，在平衡力系中应当分清楚哪些是主动力，哪些是被动力，而在写虚功方程时，只有分清楚哪些是主动力，哪些是被动力，而在写虚功方程时，只有主动力作虚功，而被动力是不作虚功的。主动力作虚功，而被动力是不作虚功的。虚功原理虚功原理表述如下：表述如下：在力的作用下处于平衡状态的体系，当发生与约束条件相在力的作用下处于平衡状态的体系，当发生与约束条件相符合的任意微

16、小的虚刚体位移时，体系上的主动力在虚位移上所符合的任意微小的虚刚体位移时，体系上的主动力在虚位移上所作的总功作的总功(各力所作的功的代数和各力所作的功的代数和)恒对于零。恒对于零。虚功原理用公式表示为：虚功原理用公式表示为：这就是这就是虚功方程虚功方程，其中，其中P和和 相应的代表力和虚位移。相应的代表力和虚位移。虚功方程是按刚体的情况得出的，即假设上图的杠杆是绝虚功方程是按刚体的情况得出的，即假设上图的杠杆是绝对刚性，没有任何的变形，因而在方程中没有内功项出现，而只对刚性，没有任何的变形，因而在方程中没有内功项出现，而只有外功项。有外功项。将虚功原理用于弹性变形时，总虚功将虚功原理用于弹性变

17、形时，总虚功 要包括外力要包括外力虚功虚功 和内力虚功和内力虚功 两部分，即：两部分，即：内力功前面有一负号，由于弹性体在变形过程中，内力内力功前面有一负号，由于弹性体在变形过程中，内力是克服变形而产生的，所有内力的方向总是与变形的方向相反，是克服变形而产生的，所有内力的方向总是与变形的方向相反，所以内力功取负值所以内力功取负值 根据虚功原理，总功等于零得：根据虚功原理，总功等于零得：外力虚功外力虚功 =内力虚功内力虚功 弹性力学中的弹性力学中的虚功原理虚功原理可表达为：在外力作用下处于平衡状可表达为：在外力作用下处于平衡状态的弹性体，如果发生了虚位移，那么所有的外力在虚位移上的态的弹性体，如

18、果发生了虚位移，那么所有的外力在虚位移上的虚功虚功(外力功外力功)等于整个弹性体内应力在虚应变上的虚功等于整个弹性体内应力在虚应变上的虚功(内力功内力功)。虚功原理虚功原理-用于弹性体的情况用于弹性体的情况虚位移虚位移：满足物体内变形连续条件，边界上位移约束条件：满足物体内变形连续条件，边界上位移约束条件的任何可能的无限小位移。的任何可能的无限小位移。虚功虚功：真实外力在虚位移上所做的功。：真实外力在虚位移上所做的功。虚应变虚应变：对可变形的弹性体，虚位移也必将导致虚应变，：对可变形的弹性体，虚位移也必将导致虚应变，虚应变和虚位移之间满足弹性体几何方程。虚应变和虚位移之间满足弹性体几何方程。虚

19、功原理虚功原理：外力作用下处于平衡状态的弹性体，外力在虚：外力作用下处于平衡状态的弹性体，外力在虚位移上做的总虚功等于弹性体内真实应力在虚应变上做的位移上做的总虚功等于弹性体内真实应力在虚应变上做的总虚变形功。总虚变形功。刚体虚位移原理刚体虚位移原理：对刚体而言，如果它处于平衡状态，则其上作用力在虚：对刚体而言，如果它处于平衡状态，则其上作用力在虚位移上做的总虚功等于零。位移上做的总虚功等于零。刚体刚体：刚体受力后，仅发生刚体位移，无变形和应变发生。：刚体受力后，仅发生刚体位移，无变形和应变发生。如何导出弹性体虚功原理？如何导出弹性体虚功原理？外力虚功外力虚功；作用在弹性体上的外力在虚位移上做

20、的功作用在弹性体上的外力在虚位移上做的功体积力体积力 表面力表面力外力虚功=体积力虚功+表面力虚功 内力虚功即应力在虚应变上做的虚功，也称虚应变能。内力虚功即应力在虚应变上做的虚功，也称虚应变能。正应力虚功切应力虚功外力虚功=总变形虚功弹性力学平面问题基础一、平面应力问题以薄板的中面为xy面，以垂直于中面的任一直线为Z轴。由于薄板两表面没有垂直和平行于板面的外力，所以板面上各点均有：另外由于平板很薄，外力又不沿厚度变化，可认为在整个薄板内各点均有：于是，在六个应力分量中，只需要研究剩下的平行于XOY平面的三个应力分量，即 ，所以称为平面应力问题。xyyx,xyyx,),(),(yxvvyxuu

21、0zyzxz)0(z二、平面应变问题0zxzyz0zxyyx,xyyx,),(),(yxvvyxuu三、平面问题基本方程和边界条件xvyuyvxuxyyx vuxyyxxyyx00物理方程中后两式可见，这时的剪应变：由第三式可见：一般，并不一定等于零，但可由 及 求得，在分析问题时不必考虑。于是只需要考虑三个应变分量即可，于是应变矩阵转化成应力分量用应变分量表示的形式：转化成应力分量用应变分量表示的形式：用矩阵方程表示：用矩阵方程表示：21EEu100yyxyxxyxfyxfyx为体力分量。，yxffuStSvvuu,uSyxyyxxyxtlmtml为边界上面力分量。，。为边界外法线方向余弦，yxttml