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

《有限元方法与ANSYS应用二讲.ppt》由会员分享，可在线阅读，更多相关《有限元方法与ANSYS应用二讲.ppt（59页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

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

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

3、何方程成为线性方程。弹性力学中几个基本概念外力外力其他物体对研究对象（弹性体）的作用力其他物体对研究对象（弹性体）的作用力 体力体力作用于物体体积内的力。如重力、惯性力和电磁力。以作用于物体体积内的力。如重力、惯性力和电磁力。以单位体积内所受的力来量度，单位体积内所受的力来量度，面力面力作用于物体表面上的力。如流体压力和接触力。以单位作用于物体表面上的力。如流体压力和接触力。以单位面积所受的力来量度面积所受的力来量度 集中力集中力如牵引力如牵引力 弹性体受外力以后，其内部将产生应力。弹性体受外力以后，其内部将产生应力。v正负号规定：正负号规定：、沿坐标正沿坐标正沿坐标正沿坐标正向为正，负向为负

4、。向为正，负向为负。向为正，负向为负。向为正，负向为负。量纲（因次）：量纲（因次）：oxyzfyfzPffx矢量矢量 方向沿方向沿 的极的极限方向限方向 xozy 、沿坐标沿坐标沿坐标沿坐标轴正向为正、负向为负轴正向为正、负向为负轴正向为正、负向为负轴正向为正、负向为负 量纲：量纲：P 方向沿方向沿方向沿方向沿 极限方向极限方向极限方向极限方向正负号规定：正负号规定：内力内力假想切开物体，截面两边互相作用的力（合力和合力矩，假想切开物体，截面两边互相作用的力（合力和合力矩，称为内力。称为内力。应力应力截面上某一点处，单位截面面积上的内力值。截面上某一点处，单位截面面积上的内力值。正应力是作用在

5、垂直于正应力是作用在垂直于x轴的面上同时也沿着轴的面上同时也沿着X轴轴方向作用的。方向作用的。剪应力是作用在垂直于剪应力是作用在垂直于X轴的面上而沿着轴的面上而沿着y轴方向作轴方向作用的。用的。应力成对出现，坐标面上的应力以正面正向，负面负向应力成对出现，坐标面上的应力以正面正向，负面负向为正。为正。(如果某一个面上的外法线是沿着坐标轴的正方向，这个如果某一个面上的外法线是沿着坐标轴的正方向，这个面上的应力就以沿坐标轴正方向为正，沿坐标轴负方向为负。面上的应力就以沿坐标轴正方向为正，沿坐标轴负方向为负。相反，如果某一个面上的外法线是沿着坐标轴的负方向，这个相反，如果某一个面上的外法线是沿着坐标

6、轴的负方向，这个面上的应力就以沿坐标轴的负方向为正，沿坐标轴正方向为负。面上的应力就以沿坐标轴的负方向为正，沿坐标轴正方向为负。)oxyzPmn矢量矢量 方向沿方向沿 的极限方向的极限方向量纲：量纲：沿截面切向和法沿截面切向和法向分解为向分解为 和和应力分量的空间方位，注意每个分量的作用面与作用方向！剪应力互等定律剪应力互等定律 作用在两个互相垂直的面上并且垂直于该两面交线的剪应力是作用在两个互相垂直的面上并且垂直于该两面交线的剪应力是互等的。互等的。(大小相等，正负号也相同大小相等，正负号也相同)。因此剪应力记号的两个角。因此剪应力记号的两个角码可以对调。码可以对调。应力分量应力分量一般说来

7、，弹性体内各点的应力状态都不相同，因此，一般说来，弹性体内各点的应力状态都不相同，因此，描述弹性体内应力状态的上述六个应力分量并不是常描述弹性体内应力状态的上述六个应力分量并不是常量，而是坐标量，而是坐标x、y、z的函数。的函数。六个应力分量的总体，可以用一个列矩阵六个应力分量的总体，可以用一个列矩阵 来表示：来表示：位移、应变、刚体位移弹性体在受外力以后，还将发生变形。物体的变形状态，一般有两种方式来描述：1、给出各点的位移；2、给出各微元体的变形 弹性体内任一点的位移，用此位移在x、y、z三个坐标轴上的投影u、v、w来表示。这三个投影称为位移分量。以沿坐标轴正方向为正，沿坐标轴负方向为负.

8、一般情况下，弹性体受力以后，各点的位移并不是定值，而是坐标的函数。微元体的变形可以分为两类：1、长度的变化，2、角度的变化。任一线素的长度的变化与原有长度的比值称为线应变(或称正应变)，用符号 来表示。沿坐标轴的线应变，则加上相应的角码，分别用 来表示。当线素伸长时，其线应变为正。反之，线素缩短时，其线应变为负。这与正应力的正负号规定相对应。任意两个原来彼此正交的线素，在变形后其夹角的变化值称为角应变或剪应变，用符号 来表示。两坐标轴之间的角应变，则加上相应的角码，分别用 来表示。规定当夹角变小时为正，变大时为负，与剪应力的正负号规定相对应(正的 引起正的 ，等等)。描述变形体的基本方程描述变

9、形体的基本方程 已知量：物体的形状和大小、材料性质、体力、边界上的面已知量：物体的形状和大小、材料性质、体力、边界上的面已知量：物体的形状和大小、材料性质、体力、边界上的面已知量：物体的形状和大小、材料性质、体力、边界上的面力或约束。力或约束。力或约束。力或约束。待求量：应力、形变和位移。待求量：应力、形变和位移。待求量：应力、形变和位移。待求量：应力、形变和位移。弹性力学的研究方法弹性力学的研究方法弹性力学的研究方法弹性力学的研究方法解法：解法：解法：解法：根据微分体上力的平衡条件，建立根据微分体上力的平衡条件，建立根据微分体上力的平衡条件，建立根据微分体上力的平衡条件，建立平衡微分方程平衡

10、微分方程平衡微分方程平衡微分方程；根据微分线段上应变位移的几何条件，建立根据微分线段上应变位移的几何条件，建立根据微分线段上应变位移的几何条件，建立根据微分线段上应变位移的几何条件，建立几何方程几何方程几何方程几何方程；根据应力应变间的物理条件，建立根据应力应变间的物理条件，建立根据应力应变间的物理条件，建立根据应力应变间的物理条件，建立物理方程物理方程物理方程物理方程在弹性体边界上：在弹性体边界上：在弹性体边界上：在弹性体边界上：根据面力条件，建立应力边界条件；根据面力条件，建立应力边界条件；根据面力条件，建立应力边界条件；根据面力条件，建立应力边界条件；根据约束条件，建立位移边界条件。根据

11、约束条件，建立位移边界条件。根据约束条件，建立位移边界条件。根据约束条件，建立位移边界条件。在弹性体内部：在弹性体内部：在弹性体内部：在弹性体内部：弹性力学基本方程 考虑微元体各个面上的考虑微元体各个面上的法向应力和剪应力法向应力和剪应力与其与其体力平衡体力平衡，注意，注意应力从一个面到对面是变化的，即有增量，将作用于微元体各个应力从一个面到对面是变化的，即有增量，将作用于微元体各个方向的力求和，略去高阶项，可得平衡方程（受力状态的描述）：方向的力求和，略去高阶项，可得平衡方程（受力状态的描述）：平衡微分方程的推导过程平衡微分方程的推导过程二、几何方程二、几何方程空间独立的应变分量空间独立的应

12、变分量正应变分量：切应变分量：应变分量列阵：如何计算正应变和切应变？正应变切应变如何描述位移分量？如何描述位移分量？几何方程几何方程几何方程描述几何方程描述应变分量应变分量与与位移分量位移分量之间的关系！之间的关系！如何推导几何方程？如何推导几何方程？几何坐标：几何位移：因此，正应变切应变切应变采用同样的思路，分析微单元体在xoz与yoz上的投影及之后的几何变形，可获得其他三个分量（1个正应变分量和2个切应变分量)。几何方程应变分量三、物理方程三、物理方程对均匀、连续、各向同性的弹性材料，存在虎克定律对均匀、连续、各向同性的弹性材料，存在虎克定律l未知数未知数 位移位移3个个+应力应力6个个+

13、应变应变6个个=15l方程个数方程个数 平衡方程平衡方程3个个+几何方程几何方程6个个+物理方程物理方程6个个=15原则上原则上15个方程可以求解个方程可以求解15个未知量，但实际上先求出一部分个未知量，但实际上先求出一部分，再通过方程求剩下的。再通过方程求剩下的。四、边界条件四、边界条件4.1 位移边界条件给出边界面上所有点的位移给定边界上位移4.2 应力边界条件应力边界条件给出边界面给出边界面abc上表面力上表面力!目前有限元法采用位移法，以三个位移分量为基目前有限元法采用位移法，以三个位移分量为基本未知量。本未知量。位移产生变形，变形产生应变，应变产生应力，位移产生变形，变形产生应变，应

14、变产生应力，应力于外力平衡。应力于外力平衡。协调的含义：不重叠、不开裂！图a示一平衡的杠杆，对C点写力矩平衡方程：图b表示杠杆绕支点C转动时的刚体位移图：综合可得：即：式是以功的形式表述的。表明：图a的平衡力系在图b的位移上作功时，功的总和必须等于零。这就叫做虚功原理。五、虚功原理五、虚功原理进一步分析。当杠杆处于平衡状态时，进一步分析。当杠杆处于平衡状态时，和和 这两个位移是这两个位移是不存在的，但是如果某种原因，例如人为地振一下让它倾斜，一定不存在的，但是如果某种原因，例如人为地振一下让它倾斜，一定满足上式的关系。将这个客观存在的关系抽象成一个普遍的原理，满足上式的关系。将这个客观存在的关

15、系抽象成一个普遍的原理，去指导分析和计算结构。去指导分析和计算结构。对于在力的作用下处于平衡状态的任何物体，不用考虑它是否真对于在力的作用下处于平衡状态的任何物体，不用考虑它是否真正发生了位移，而假想它发生了位移，正发生了位移，而假想它发生了位移，(由于是假想，故称为虚位由于是假想，故称为虚位移移)，那么，物体上所有的力在这个虚位移上的总功必定等于零。，那么，物体上所有的力在这个虚位移上的总功必定等于零。这就叫做这就叫做虚位移原理虚位移原理，也称，也称虚功原理虚功原理。在图。在图a中的中的 和和 所作所作的功就不是发生在它本身的功就不是发生在它本身(状态状态a)的位移上，的位移上，(因为它本身

16、是平衡的，因为它本身是平衡的，不存在位移不存在位移)，而是在状态，而是在状态(b)的位移上作的功。可见，这个位移对的位移上作的功。可见，这个位移对于状态于状态(a)来说就是虚位移，亦即是状态来说就是虚位移，亦即是状态(a)假象的位移。假象的位移。必须指出，虚功原理的应用范围是有条件的，它所涉及必须指出，虚功原理的应用范围是有条件的，它所涉及到的两个方面，力和位移并不是随意的。对于到的两个方面，力和位移并不是随意的。对于力力来讲，它必须是来讲，它必须是在位移过程中处于在位移过程中处于平衡的力系平衡的力系；对于；对于位移位移来讲，虽然是虚位移，来讲，虽然是虚位移，但并不是可以任意发生的。它必须是和

17、约束条件相符合的微小的但并不是可以任意发生的。它必须是和约束条件相符合的微小的刚体位移。刚体位移。还要注意，当位移是在某个约束条件下发生时，则在该还要注意，当位移是在某个约束条件下发生时，则在该约束力方向的位移应为零，因而该约束力所作的虚功也应为零。约束力方向的位移应为零，因而该约束力所作的虚功也应为零。这时该约束力叫做这时该约束力叫做被动力被动力。(如上图中的反力如上图中的反力 ，由于支点，由于支点C没没有位移，故有位移，故 所作的虚功等于零所作的虚功等于零)。反之，如上图中的。反之，如上图中的 和和 是是在位移过程中作功的力，称为主动力。因此，在平衡力系中应当在位移过程中作功的力，称为主动

18、力。因此，在平衡力系中应当分清楚哪些是主动力，哪些是被动力，而在写虚功方程时，只有分清楚哪些是主动力，哪些是被动力，而在写虚功方程时，只有主动力作虚功，而被动力是不作虚功的。主动力作虚功，而被动力是不作虚功的。虚功原理虚功原理表述如下：表述如下：在力的作用下处于平衡状态的体系，当发生与约束条件相在力的作用下处于平衡状态的体系，当发生与约束条件相符合的任意微小的虚刚体位移时，体系上的主动力在虚位移上所符合的任意微小的虚刚体位移时，体系上的主动力在虚位移上所作的总功作的总功(各力所作的功的代数和各力所作的功的代数和)恒对于零。恒对于零。虚功原理用公式表示为：虚功原理用公式表示为：这就是这就是虚功方

19、程虚功方程，其中，其中P和和 相应的代表力和虚位移。相应的代表力和虚位移。虚功方程是按刚体的情况得出的，即假设上图的杠杆是绝虚功方程是按刚体的情况得出的，即假设上图的杠杆是绝对刚性，没有任何的变形，因而在方程中没有内功项出现，而只对刚性，没有任何的变形，因而在方程中没有内功项出现，而只有外功项。有外功项。将虚功原理用于弹性变形时，总虚功将虚功原理用于弹性变形时，总虚功 要包括外力要包括外力虚功虚功 和内力虚功和内力虚功 两部分，即：两部分，即：内力功前面有一负号，由于弹性体在变形过程中，内力内力功前面有一负号，由于弹性体在变形过程中，内力是克服变形而产生的，所有内力的方向总是与变形的方向相反，

20、是克服变形而产生的，所有内力的方向总是与变形的方向相反，所以内力功取负值所以内力功取负值 根据虚功原理，总功等于零得：根据虚功原理，总功等于零得：外力虚功外力虚功 =内力虚功内力虚功 弹性力学中的弹性力学中的虚功原理虚功原理可表达为：在外力作用下处于平衡状可表达为：在外力作用下处于平衡状态的弹性体，如果发生了虚位移，那么所有的外力在虚位移上的态的弹性体，如果发生了虚位移，那么所有的外力在虚位移上的虚功虚功(外力功外力功)等于整个弹性体内应力在虚应变上的虚功等于整个弹性体内应力在虚应变上的虚功(内力功内力功)。虚功原理虚功原理-用于弹性体的情况用于弹性体的情况虚位移虚位移：满足物体内变形连续条件

21、，边界上位移约束条件：满足物体内变形连续条件，边界上位移约束条件的任何可能的无限小位移。的任何可能的无限小位移。虚功虚功：真实外力在虚位移上所做的功。：真实外力在虚位移上所做的功。虚应变虚应变：对可变形的弹性体，虚位移也必将导致虚应变，：对可变形的弹性体，虚位移也必将导致虚应变，虚应变和虚位移之间满足弹性体几何方程。虚应变和虚位移之间满足弹性体几何方程。虚功原理虚功原理：外力作用下处于平衡状态的弹性体，外力在虚：外力作用下处于平衡状态的弹性体，外力在虚位移上做的总虚功等于弹性体内真实应力在虚应变上做的位移上做的总虚功等于弹性体内真实应力在虚应变上做的总虚变形功。总虚变形功。刚体虚位移原理刚体虚

22、位移原理：对刚体而言，如果它处于平衡状态，则其上作用力在虚：对刚体而言，如果它处于平衡状态，则其上作用力在虚位移上做的总虚功等于零。位移上做的总虚功等于零。刚体刚体：刚体受力后，仅发生刚体位移，无变形和应变发生。：刚体受力后，仅发生刚体位移，无变形和应变发生。如何导出弹性体虚功原理？如何导出弹性体虚功原理？外力虚功外力虚功；作用在弹性体上的外力在虚位移上做的功作用在弹性体上的外力在虚位移上做的功体积力体积力 表面力表面力外力虚功=体积力虚功+表面力虚功 内力虚功即应力在虚应变上做的虚功，也称虚应变能。内力虚功即应力在虚应变上做的虚功，也称虚应变能。正应力虚功切应力虚功外力虚功=总变形虚功弹性力

23、学平面问题基础 任何实际变形体的力学问题都是空间问题（三维问题），任何实际变形体的力学问题都是空间问题（三维问题），所受的外力一般都是空间力系。所受的外力一般都是空间力系。在某些特殊情况下，比如物体具有特殊形状，受特殊的外力，特在某些特殊情况下，比如物体具有特殊形状，受特殊的外力，特殊的位移约束时，空间问题就可以简化成平面问题。此时，问题殊的位移约束时，空间问题就可以简化成平面问题。此时，问题的几何和力学量仅仅是二维坐标的函数。所求未知力学量只是二的几何和力学量仅仅是二维坐标的函数。所求未知力学量只是二维空间内的分量。维空间内的分量。这种平面问题模型下，所得到的结果能满足工程上的精度要求，这种

24、平面问题模型下，所得到的结果能满足工程上的精度要求，而分析计算工作量大大减少。而分析计算工作量大大减少。大量的固体力学问题都可以简化为平面问题。大量的固体力学问题都可以简化为平面问题。平面问题包括：平面应力问题和平面应变问题。平面问题包括：平面应力问题和平面应变问题。一、平面应力问题 平面应力问题的基本特征：平面应力问题的基本特征：平平面面应应力力问问题题模模型型1 1）几何特征）几何特征2 2）受力特征）受力特征vv薄板的两个侧面上无载荷作用；薄板的两个侧面上无载荷作用；vv边缘上受到平行于板面且沿板厚均匀分布的面力作用；边缘上受到平行于板面且沿板厚均匀分布的面力作用；vv体力平行于板面且不

25、沿板厚变化（体力平行于板面且不沿板厚变化（x,yx,y的函数）。的函数）。vv物体在一个方向（物体在一个方向（z)z)的尺寸远远小于其它两个方向（的尺寸远远小于其它两个方向（x,y)x,y)的尺寸。的尺寸。几何为均匀薄板。几何为均匀薄板。以薄板的中面为xy面，以垂直于中面的任一直线为Z轴。由于薄板两表面没有垂直和平行于板面的外力，所以板面上各点均有：另外由于平板很薄，外力又不沿厚度变化，可认为在整个薄板内各点均有：于是，在六个应力分量中，只需要研究剩下的平行于XOY平面的三个应力分量，即 ，所以称为平面应力问题。vv应力分量：应力分量：x,yx,y的函数的函数vv应变分量：应变分量：x,yx,

26、y的函数的函数vv位移分量：位移分量：（非独立）（非独立）平面应力问题的应力、应变、位移分量平面应力问题的应力、应变、位移分量平面应力问题的例子平面应力问题的例子由基本特征推出：由基本特征推出：二、平面应变问题平面应变问题的例子平面应变问题的例子 平面应变问题的基本特征：平面应变问题的基本特征：1 1）几何特征）几何特征vv一个方向（一个方向（z)z)尺寸远远大于其它两个方向（尺寸远远大于其它两个方向（x,y)x,y)的尺寸的尺寸,呈现为无呈现为无限长等截面柱体。限长等截面柱体。2 2）受力特征）受力特征vv外力（体力、面力）平行于横截面作用，且沿纵向不变化。外力（体力、面力）平行于横截面作用

27、，且沿纵向不变化。平面应变问题的应力、应变、位移分量平面应变问题的应力、应变、位移分量（非独立）（非独立）vv应变分量：应变分量：x,yx,y的函数的函数x,yx,y的函数的函数vv应力分量：应力分量：vv位移分量：位移分量：平面应变问题的例子平面应变问题的例子由基本特征推出：由基本特征推出：三、平面问题基本方程和边界条件 平面应力和平面应变问题都归结为求解平面内的平面应力和平面应变问题都归结为求解平面内的3 3个应力分量、个应力分量、3 3个应变分量、个应变分量、2 2个位移分量。个位移分量。要求解这些未知力学量，需要通过研究弹性体的平衡、几何、要求解这些未知力学量，需要通过研究弹性体的平衡

28、、几何、物理关系得到足够的方程。物理关系得到足够的方程。一）平面问题几何方程一）平面问题几何方程应变应变 位移关系位移关系算子矩阵算子矩阵vv几何方程对于平面应力和平面应变问题相几何方程对于平面应力和平面应变问题相同同物理方程中后两式可见，这时的剪应变：由第三式可见：一般，并不一定等于零，但可由 及 求得，在分析问题时不必考虑。于是只需要考虑三个应变分量即可，于是应变矩阵二）平面问题的物理方程二）平面问题的物理方程应力应力 应变关系应变关系转化成应力分量用应变分量表示的形式：转化成应力分量用应变分量表示的形式：用矩阵方程表示：用矩阵方程表示：平面应力弹性矩平面应力弹性矩阵，对称方阵。阵，对称方

29、阵。vv平面应变问题弹性矩阵按如下办法得到。平面应变问题弹性矩阵按如下办法得到。vv因此，对于平面问题的推导和编程，因此，对于平面问题的推导和编程，只按平面应力问题处理。只按平面应力问题处理。将平面应力弹性矩阵中的弹性常数作如下变换：将平面应力弹性矩阵中的弹性常数作如下变换：通过研究微分单元体平衡，得到下列平衡微分方程：通过研究微分单元体平衡，得到下列平衡微分方程：三）平面问题的平衡方程三）平面问题的平衡方程四）平面问题边界条件四）平面问题边界条件vv平面弹性体的边界分为位移边界平面弹性体的边界分为位移边界 和应力边界和应力边界vv通过前面的基本方程求解弹性力通过前面的基本方程求解弹性力学问题时，必须考虑上述边界学问题时，必须考虑上述边界上位移的协调和力的平衡。上位移的协调和力的平衡。vv边界条件描述如下：边界条件描述如下：在在 上上2 2）应力边界条件）应力边界条件1 1）位移边界条件）位移边界条件