弹塑性力学总结.docx

弹塑性力学总结弹塑性力学的任务是分析各种结构物或其构件在弹性阶段和塑性阶段的应力 和位移，校核它们是否具有所需的强度、刚度和稳定性，并寻求或改进它们的计 算方法。并且弹塑性力学是以后有限元分析、解决具体工程问题的理论基础，这 就要求我们掌握其必要的基础知识和具有一定的计算能力。通过一学期的弹塑性 力学的学习，对其内容总结如下：一、弹性力学1、弹性力学的基本假定求解一个弹性力学问题，通常是已知物体的几何形状（即已知物体的边界）， 弹性常数，物体所受的外力，物体边界上所受的面力，以及边界上所受的约束； 需要求解的是物体内部的应力分量、应变分量与位移分量。求解问题的方法是通 过研究物体内部各点的应力与外力所满足的静力平衡关系，位移与应变的几何学 关系以及应力与应变的物理学关系，建立一系列的方程组；再建立物体表面上给 定面力的边界以及给定位移约束的边界上所给定的边界条件；最后化为求解一组 偏分方程的边值问题。在导出方程时，如果考虑所有各方面的因素，则导出的方程非常复杂，实际 上不可能求解。因此，通常必须按照研究对象的性质，联系求解问题的范围，做 出若干基本假定，从而略去一些暂不考虑的因素，使得方程的求解成为可能。（1）假设物体是连续的。就是说物体整个体积内，都被组成这种物体的物质 填满，不留任何空隙。这样，物体内的一些物理量，例如：应力、应变、位移等， 才可以用坐标的连续函数表示。（2）假设物体是线弹性的。就是说当使物体产生变形的外力被除去以后，物 体能够完全恢复原来形状，不留任何残余变形。而且，材料服从虎克定律，应力 与应变成正比。（3）假设物体是均匀的。就是说整个物体是由同一种质地均匀的材料组成的。 这样，整个物体的所有部分才具有相同的物理性质，因而物体的弹性模量和泊松 比才不随位置坐标而变。（4）假设物体是各向同性的。也就是物体内每一点各个不同方向的物理性质 和机械性质都是相同的。（5）假设物体的变形是微小的。即物体受力以后，整个物体所有各点的位移 都小于物体的原有尺寸，因而应变和转角都远小于1。这样，在考虑物体变形以 后的平衡状态时，可以用变形前的尺寸代替变形后尺寸，而不致有显著的误差； 并且，在考虑物体的变形时，应变和转角的平方项或乘积都可以略去不计，使得 弹性力学中的微分方程都成为线性方程。2、外力和应力的概念作用于弹性体的外力可以分为体（积）力和（表）面力。体力是分布在弹性 体体积内质量上的力，例如重力和惯性力、磁力等。在物体内任一点的体力，用 作用于其上的单位体积的体力沿坐标轴上的投影X、Y、Z来表示。它们的指向以 沿坐标轴正方向为正；反之为负。这三个投影称为该点的体力分量。面力是指作用于弹性体表面上的外力，例如流体压力和接触力等。可以是分 布力，也可以是集中力。在弹性表面上任一点的面力，用作用于其上的单位面积 上面力沿坐标轴上的投影X、Y、z来表示。它们的指向也以沿坐标轴正方向的 为正，反之为负。这三个投影称为该点的面力分量。弹性体在外力作用下变形，而在弹性体内部为了阻止其变形就产生了内力来 平衡外力。作用在单位面积上的内力称为应力。3、一点的应力状态为了研究弹性体内任一点p的应力，就在这一点设想从弹性体中取出一个微分体（无限小的平行六面体）如下图1：图1微小平行六面体的应力状态如果某一个截面上的外法线是沿着坐标轴的正方向，这个截面就称为一个正面，而这个面上的应力分量就以沿坐标轴正方向为正，沿坐标轴负方向为负。相反，如果某一个截面上的外法线是沿着坐标轴的负方向，这个截面就称为一个负面，而这个面上的应力分量就以沿坐标轴负方向为正，沿坐标轴正方向为负。图 上所示的应力分量全部都是正的。注意，虽然上述正负号规定对于正应力说来， 结果是和材料力学中的规定相同（拉应力为正而压应力为负），但是，对于剪应 力说来，结果却和材料力学中的规定不完全相同。剪应力的互等关系：作用在两个互相垂直的面上并且垂直于该两面交线的剪 应力，是互等的（大小相等，正负号也相同）。T =T ， T =T ， T =T（ 1）yz zy 女 yr4、斜截面应力公式，物体表面给定力的边界条件现在，假定物体在任一点P的六个应力分量。、。、。、T、T、T为已知， r y z yz zr ry试求经过P点的任一斜面上的应力。为此，在P点附近取一个平面ABC，平行于 这一斜面，并与经过P点而平行于坐标面的三个平面形成一个微小的四面体 PABC，如图2所示。当平面ABC趋近于P点时，平面ABC上的应力就成为该斜 面上的应力。图2物体内任意一点的应力状态设斜面ABC的向外法线为N，而N的方向余弦为:l = cos （n, r）m = cos （n, y ）n =cos（n, z）Z弓=0可得出与上式相似的两个方程。由平衡条件2fx= 0、&广0及简化后三个方程为：X = lb + mT + nTY = mb + nT + ItZ = nb . + It + mT设三角形ABC上的正应力为。,则由投影可得：Nb = 12b + m2。+ n2b + 2mn + 2nh + 2lm（ 4）设三角形ABC上的剪应力为t n，则由于：S 2 =b 2 +T 2 = X 2 + Y 2 + Z 2（ 5）而有：T2 = X 2 +Y2 + Z2 -b2（ 6）由公式（4）和（5）可见，在物体的任意一点，如果已知六个应力分量b、b、b、T、T、T ，就可以求得任一斜面上的正应力和剪应力。因此，可 工 y z 兴 吞 母以说，六个应力分量完全决定了一点的应力状态。在特殊情况下，如果ABC是物体的边界面，则Xn、'、Zn成为面力分量X、Y、Z，于是由公式（3）得出：lb + mT + nT = Xxyxzxmb + nT + lT=Y（7）yzyxynb + lT + mT=Zzxzyz这就是弹性体的应力边界条件，它表明应力分量的边界值与面力分量之间的 关系。5、应力分量的坐标转换关系若物体处在某一确定的应力状态，在某一组坐标系中，这个应力状态可以用 六个应力分量b表示，在另一组坐标系中，同一个应力状态却以另外一组不同的 ij应力分量b表示。两组应力分量之间应力满足一定的坐标转换关系。在物体上任ll2l3”1b 11n；21331T12b22T32T13；II m b23 h331l3m3n3孔L L3TbT212223匕T 32E上式也可以表示成抽象的矩阵乘式：b 疽 kb P jl（ 9）例如：若第一组坐标系为直角坐标系（X、y、乙），第二组坐标系为圆柱坐标系 （r、0、z），可知两组坐标系的转换矩阵为：cos 0 sin 0 0（10）虹=一 sin 0cos 0 0'0016、主应力、应力主方向、主剪应力若经过物体中一点P处的某一斜面上的剪应力等于零，则该斜面上的正应力 称为P点的一个主应力，该斜面称为P点的一个主应力面，而该斜面的垂线方向 称为p点的一个主应力方向。可以证明，在弹性体的任一点，一定存在三个相互垂直的主应力面及和它们 对应的三个主应力，通常用匕、。2、七。而且，任何一个斜面上的正应力都不会 大于三个主应力中最大的一个，也不会小于三个主应力中最小的一个。主应力与 主方向可以用以下的方法求得：假设N是P点应力状态a.的一个主方向，N与原始坐标系X、必z的夹角方 IJ向余弦为1、m、n，它们间总满足：（11）在垂直于N的截面上只有正应力a （某个主应力）作用，则由柯西公式知：（12）lb + mr + nr = la ma + nr + 1r = ma na + 1r + mr = na上式中1、m、n为待求的方向余弦，将上式移项可以得到求解的齐次线性方 程组:1 (a a)+ mr + nr = 0（13）1r + mG a)+ nr . = 0Ir . + m . + n(a . a)= 0方程（13）零解的条件是其系数行列式值为零，即：(ax-a)(r 虾)七一 rva -a Jr =a 3 -1 a 2 +1 a -1 = 0( 14)r r(a a)式（14）称为该应力状态的特征方程式，它是一个三次代数方程，可以证明 它有三个实根，称为特征根，就是应力状态a.所对应的主应力。可以证明，特征 ij方程（14）式的系数I，I，I是只与应力状态有关，与所选择的原始坐标系无关的 123量，分别称为该应力状态的第一、第二、第三不变量。即I =b +b +b1x yz"b T 一"bT"bT 一I =xyx+yzy+xzx2T bTbT b1- xyyu yzzxzzb TTxyxzxI =TbT3xyyzyT Tb1- xzyzz（15）（16）（17）7、叠加原理与圣维南原理在解决一个弹性力学问题时，我们常常利用叠加原理来有效地处理各种复杂 载荷作用的情况。叠加原理是：考虑同一物体受两组载荷作用，第一组为体力f ； （i = 1，2,3）和面力 F （i = 1，2,3）；第二组为体力f" （i = 1，2,3）和面力F （i = 1，2,3），它们引起的应力和 iii全移场分别为。和u以及C '和u （i，j = 1，2,3）。如果物体处于线弹性、小变形状 l iij i态，两组载荷同时作用时物体内的应力和位移场等于它们单独作用时相应的应力 与位移场之和。弹性理论要求在物体的每个边界点上都给定边界条件。实际工程问题却往往 只知道总的载荷量，只能提出等效的近似边界条件，给不出详细的载荷分布规律。 另外，解题时往往难于满足逐点给定的精确边界条件，因而也希望能找到一种边 界条件的简化方案。圣维南原理指出：由作用在物体局部表面上的自平衡力系（即合力与合力矩为 零的力系），所引起的应变，在远离作用区（距离远大于该局部作用区的线性尺寸） 的地方可以忽略不计。圣维南原理的另一种提法是：若把作用在物体局部表面上的 外力，用另一组与它静力等效的力系来代替。则这种等效处理对物体内部应力应 变状态的影响将随远离作用区的距离增加而迅速衰减。显然，上述两种提法是完 全等效的。8、平面问题的基本方程平衡微分方程:8b机八（18）云 F+x=08b机y + + Y = 08y8x几何方程:du& =xdx(19 )dvy d-yYxydvdu一 + 一dx dy物理方程:_pG +b )iyz8 y ELb _G +b)yz x -1(20 )1Y yz = G T yz1Y zx = G Lx1Y = 匚xy G xy9、平面应力问题与平面应变问题平面应力：只在平面内有应力，与该面垂直方向的应力可忽略，例如薄板拉 压问题。平面应变：只在平面内有应变，与该面垂直方向的应变可忽略，例如水坝侧 向水压问题。具体说来：平面应力是指所有的应力都在一个平面内，如果平面是oxy平面， 那么只有正应力b、b和剪应力t （它们都在一个平面内），没有b、T、匚。x yxyz yz zx平面应变是指所有的应变都在一个平面内，同样如果平面是oxy平面，则只有正 应变8、8和剪应变Y ，而没有8、Y、Y。x yxyz yz zx举例说来：平面应变问题比如压力管道、水坝等，这类弹性体是具有很长的 纵向轴的柱形物体，横截面大小和形状沿轴线长度不变；作用外力与纵向轴垂直， 并且沿长度不变；柱体的两端受固定约束。平面应力问题讨论的弹性体为薄板， 薄壁厚度远远小于结构另外两个方向的尺度。薄板的中面为平面，其所受外力， 包括体力均平行于中面面内，并沿厚度方向不变。而且薄板的两个表面不受外力 作用。10、弹性力学的基本方法在弹性力学里求解问题，主要有三种基本方法，分别是按位移求解、按应力 求解和混合求解。按位移求解时，以位移分量为基本未知函数，根据基本方程和边界条件求出 位移分量，从而求出其他分量。按应力求解一般有逆解法和半逆解法。所谓逆解法，就是先设定各种形式的、 满足相容方程的应力函数中，从而求出应力分量。然后根据应力边界条件来考察， 在各种形状的弹性体上，这些应力分量对应于什么样的面力，从而得知所设定的 应力函数可以解决什么问题。所谓半逆解法，就是针对所要解的问题，根据弹性 体的边界形状和受力情况，假设部分或全部应力分量为某种形式的函数，从而推 出应力函数中，然后来考察这个应力函数是否满足相容方程以及原来假设的应力 分量和由这个应力函数求出其他应力分量，是否满足应力边界条件和位移单值条 件。相容方程:七2 a2）（）（21）疝+寥r x y=二、塑性力学1、塑性力学的基本假设当作用在物体上的外力取消后，物体的变形不完全恢复，而产生一部分永久 变形时，这中变形为塑性变形。在实验的基础上，塑性力学一般采用以下假设:（1）材料是连续的，均匀的。（2）平均正应力（静水压力）不影响屈服条件和加载条件。（3）体积的变化是弹性的。（4）不考虑时间因素对材料性质的影响。2、变形体的模型对于不同的材料，不同的应用领域，我们可以采用不同的变形体的模型，这 种模型必须符合材料的实际性质。不同的材料有不同的拉伸曲线，但它们具有一 些共同性质。其拉伸曲线图如图3。3、屈服条件对于处于单向拉伸（或压缩）的物体，当应力达到屈服极限时，材料开始进 入塑性状态，对于处于复杂应力状态的物体，由弹性状态过渡到塑性状态的临界 条件称为屈服条件。在应力空间将初始屈服的应力点连成的弹性和塑性的分界面 称为屈服面。描述屈服面的数学表达式称为屈服函数。常用的各向同性金属材料 的屈服试验表明，屈服应力数据点介于屈雷斯卡（Tresca）屈服条件和密赛斯 （Mises）屈服条件之间，而更接近于密赛斯屈服条件。1）、屈雷斯卡屈服条件（最大切应力条件）屈雷斯卡屈服条件为：当最大切应力达到某一极限值时，材料开始进入塑性 状态，即b b =b, b >Q >Q(22)在主应力空间，当差值b bj、b2 b3|、b3 bj中任意一个达到2左时，材 料进入塑料性状态，即'b 2b J 京乍3b| 3b b <bVl 2 i s(23)因此用屈雷斯卡条件表示的屈服面为由下列六个平面组成的正六边形柱体。如图5所示:图5在主应力空间中Mises和Tresca屈服条件材料常数k由实验确定。在拉伸试验时，b广2k =气，即k =b s/2。在纯剪切 试验时，b1 b3 - 2k - 2t ,即k -T。如果屈雷斯卡条件成立，必有T = b /2。2)、密赛斯屈服条件密赛斯条件为：当切应力强度等于剪切屈服极限T I者当应力强度七等于拉伸屈服极限气时，材料开始屈服，(b b)2 +(b b)2 +(b b)2 - 2b 2(b b ) +(b b ) +(b b)2 + 6 ( 2 +T 2 +T时，材料开始屈服；或2)= 2b 2女s(24)对于密赛斯条件，c =b 。密赛斯条件与屈雷斯卡条件的最大差别不超过15%。在主应力空间，密赛斯屈服面为一外接于屈雷斯卡屈服面的圆柱面。在平面应力 状态，设七=0，则在。1、。2应力平面上，密赛斯条件为一椭圆，屈雷斯卡条件 为内接六边形（图6）。-1七/七-1图6当° 3 = 0时的Mi ses和Tresca屈服条件4、塑性应力应变关系塑性应力应变关系有增量（流动）理论和全量（形变）理论两种类型。1）、增量理论全量理论用应力和应变的瞬时值表示的塑性应力应变关系，是塑性应力应变 增量关系沿加载途径的积分形式。当满足小变形及简单加载（应力分量成比例增 长）条件，应力强度和应变强度之间存在单一的函数关系。2）、全量理论材料在塑性变形时，应力与应变之间一般不存在一一对应的关系。增量理论 假设在塑性流动的任一瞬时，塑性应变增量矢量与加载面正交。三、ANSYS求解实例为加深对理论的理解，将书上的习题用ANSYS来实现，可以将习题的理论 结果和ANSYS计算的数值结果进行对比。以下用弹性力学求解方法和ANSYS有 限元分析分析书上的一个习题。设图7中三角形悬臂梁中受重力的作用，而梁的密度为P，试用纯三次式的 应力函数求解。图71、运用弹性力学求解方法:设此应力函数为中=Ax3 + Bx2y + Cxy2 + Dy3b = 2Cx + 6 Dy可得 a = 6Ax + 2By-pgy > y t = -2 Bx - 2Cy xy考虑边界条件对于上界，有G） = 0 , G ） = 0 y y=0xy y=0得 A = 0，B = 0对于斜面AD面有-sin 以（2C + 6 Dtg）+ cos 以（一2Ctg以）=0cos 以（-p gtg以）+ sin 以（一2Ctg以）=0得 C = -p2gctga，D = -Pgctg 2ab x = P gxctga - 2 p gyctg 2a b =-P gyyt =-p gyctgaxy2、运用ANSYS有限元分析分析以下图8、图9、图10分别为应力分量七、a y、t xy分布图MLHirtL SOLUnCiB3TEF=1SUB =1TIME=1ANSYSJ辎 £5 £407 赤：蛔：罚图8应力分量。x分布图图9应力分量° 分布图图10应力分量T ”分布图3、两种方法结果对比。令AB =用，BC = 1, p = 1, g=9.8 ,对于悬臂梁上点E比分析。分析如下表。"，Z，B1 V3 7侦，0）处对bxbyTxy弹性力学求解方法ANSYS有限元分 析分析弹性力学求解方法ANSYS有限元分 析分析弹性力学求解方法ANSYS有限元分 析分析B S，0 )29.453.3010.20-10.5T招-16.974-16.125-5.658-5.533-9.800-9.610绝对误差B : 23.9E: 0.849B : 10.2E: 0.125B : 10.5E: 0.190相对误差B : 0.82E: 0.05B : 8E: 0.02B : 8E: 0.02E点处3个应力分量误差很小，两种分析方法所得结果相当吻合；B点处3 个应力分量误差相当大。从图8、9、10中可以看出在B点处进行ANSYS有限元 分析分析时出现很大的偏差。因为ANSYS有限元分析分析具体问题时与网格划 分、采用的具体算法以及建模时所考虑的边界条件等各种因素有一定的关系。以上是我对弹塑性力学的总结，还有诸多纸漏之处，理解不够透彻，望老师 批评指正。