塑性力学-第二章应力状态

1、第二章应力状态 第二章应力状态理论2.1应力和应力张量在外力作用下，物体将产生应力和变形，即物体中诸元素之间的相对位置发生变化，由于这种变化，便产生了企图恢复其初始状态的附加相互作用力。用以描述物体在受力后任何部位的内力和变形的力学量是应力和应变。本章将讨论应力矢量和某一点处的应力状态。o,即为了说明应力的概念，假想把受组平衡力系作用的物体用一平面A分成A和B两部分（图2.1）。如将B部分移去，则B对A的作用应代之以B部分对A部分的作用力。这种力在B移去以前是物体内A与B之间在截面C的内力，且为分布力。如从C面上点P处取出一包括P点在内的微小面积元素AS，而AS上的内力矢量为AF，则内力的平均

2、集度为AF/AS，如令AS无限缩小而趋于点P,则在内力连续分布的条件下AF/AS趋于一定的极.AFlim=bAStOAS这个极限矢量b就是物体在过c面上点P处的应力。由于AS为标量，故,b的方向与AF的极限方向一致。内力矢量AF可分解为所在平面的外法线方向和切线方向两个分量AF和AF。ns图2.1应力矢量，AS面上的正应力和切应力同样，应力b可分解为所在平面的外法线方向和切线方向两个分量。沿应力所在平面的外法线方向n的应力分量称为正应力，记为b，沿切线方向的应力分量称为切应力，记为nt。此处脚注n标明其所在面的外法线方向，n分别为在上面的讨论中，过点P的平面C是任选的。显然，过点P可以做无穷多

3、个这样的平面C,也就是说，过点P有无穷多个连续变化的n方向。不同面上的应力是不同的。这样，就产生了如何描绘一点处的应力状态的问题。为了研究点P处的应力状态，在点P处沿坐标轴x,y，z方向取一个微小的平行六面体（图2.2），其六个面的外法线方向分别与三个坐标轴的正负方向重合,其边长分别为Ax,y，z。假定应力在各面上均匀分布，于是各面上的应力便可用作用在各面中心点的一个应力矢量来表示，每个面上的应力矢量又可分解关一个正应力和两个切应力分量，如图2.2所示。以后，对正应力只用一个字母的下标标记，对切应力则用两个字母标记*其中第一个字母表示应力所在面的外法线方向；第二个字母表示应力分量的指向。正应力

4、的正负号规定为：拉应力为正，压应力为负。切应力的正负早规定分为两种情况：当其所在面的外法线与坐标轴的正方向一致时，则以沿坐标轴正方向的切应力为正反之为负；当所在面的外法线与坐标袖的负方向一致时，则以沿坐标轴负方向的切应力为正，反之为负。图2.2中的各应力分量均为正。应力及其分量的单位为Pa。图2.2应力表示法由图2.2可知，当微小的平行六面体趋于无穷小时，六面体上的应力就代表一点处的应力。因此，一点处的应力分量共有9个，其中有3个正应力分量、6个切应力分量，由切应力互等定理可知，实际上独立的切应力分量只有3个。把这9个应力分量按一定规则排列，令其中每一行为过一点的一个面上的3个应力分量，即得如

5、下应力张量，在数学上称之为二阶张量。CTTxxyxzTOC o 1-5 h zC=TCTijyxyyzTTCzxzyz其中i,j=(x,y,z)，当i,j任取x,y,z时，则得到相应的应力分量,但c,c,c分别简写为c,C,c。xxyyzzxyz应当指出，物体内各点的应力状态，一般并不是相同的，即非均匀分布的，因此各点的应力分量是坐标z,y,z的函数。所以，应力张量C与给定点的空间ij位置有关，同时应力张量是针对物体中的某一确定点而言的，今后将会看到，应力张量完全确定了一点处的应力状态。张量符号与下标记号法使冗长的弹塑性力学公式变得简明醒目，在文献中已被广泛应用，今后将逐渐熟悉这种标记法。2.

6、2二维应力状态与平面问题的平衡微分方程式上节中讨论应力概念时，是从三维受力物体出发的，其中点P是从一个三维空间中取出来约点。为简单起见，首先讨论平面问题。掌提了平面问题以后再讨论空间问题就比较容易了。当受载物体所受的面力和体力以及其应力都与某一个坐标轴(例如z轴)无关。平面问题又分为平面应力问题与平面应变问题。1.平面应力问题如果考虑如图2.3所示物体是一个很薄的平板，荷载只作用在板边，且平行于板面，即xy平面，z方向的体力分量Z及面力分量F均z为零，则板面上(z=5/2处)应力分量为(G)广0zz=丘2(T)x=(T)x=0zxz=5zyz=b22因板的厚度很小，外荷载又沿厚度均匀分布，所以

7、可以近似地认为应力沿厚度均匀分布。由此，在垂直于z轴的任一微小面积上均有G=0，T=T=0图2.3平面应力问题zzxzyT=TyxxzG=ijGxTyx0TxyGy0根据切应力互等定理，即应力张量的对称性，必然有T=T=0。因而对于平面应力状态的应力张量为也可写为TxyGy如果Z方向的尺寸为有限量，仍假设G=0,zijGxTyxT=T=0，且认为g,Gzxzyxy和T(T)为沿厚度的平均值，则这类问题称为广义平面应力问题。2.平面应变问题如果物体纵轴方向(OZ坐标方向)的尺寸很长，外荷载及体力为沿z轴均匀分布地作用在垂直于oz方向，如图2.4所示的水坝是这类问题的典型例子。忽略端部效应，则因外

8、载沿Z轴方向为一常数，因而可以认为，沿纵轴方向各点的位移与所在Z方向的位置无关，即Z方向各点的位移均相同。令u、v、w分别表示一点在x、y、Z坐标方向的位移分量，则有w为常数。等于常数的位移w并不伴随产生任一xy平面的翘曲变形，故研究应力、应变问题时,可取w=0。此外，由于物体的变形只在xy平面内产生,因此w与z无关。故对于平面应变状态有u二u(x,y)v二v(x,y)w=0由对称条件可知，在xy平面内T(T)和T(T)xzzxyzzy恒等于零，但因Z方向对变形的约束，故C一般并不z为零，所以其应力张量为TOC o 1-5 h zQxQyQzQtQtQo+z0QxQyQz(2.2-3)而G当i

9、=x,y,z时，则分别i(2.2-4)式(2.2-4)即是不计体力时们三维平衡微分方程式。2.3一点的应力状态所谓一点的应力状态是指受力变形物体内一点的不同截面上的应力变化的状况。现以平面问题为例说明一点处应力状态。在受力物体中取一个如图2.6所示的微小三角形单元，其中AC，AB与坐标轴x,y重合，而BC的外法线与zz轴成0角。取坐标x,y,使BC的外法线方向与x方向重合(如图2.6)。如果o,o,txyxy已知，则BC面上的正应力o，和切应力T可用已知量表示。因0角的任意性,xxy若BC面趋于点A时，则可认为求得了描绘过点4处的应力状态的表达式。实际上，这里所讨论的问题是一点处不同方向的面上

10、的应力的转换，即BC面无限趋于点A时，该面上的应力如何用与原坐标相平行的面上的应力来表示。在这种问题的分析中，可不必引入应力增量和体力，因为它们与应力相比属于小量。假定BC的面积为1，则AB和AC的面积分别为cos0与sin0。于是，由力在坐标x,y的平衡条件图2.6一点的应力状态工F=0和工F=0,xy可得p=gcos0+tsin0 xxxy(a)p=tcos0+gsin0yxyy式中p,p为BC面上单位面积的力p在坐标轴x,y方向上的分力(图2.6)。将xyp,p投影到x,y坐标轴方向，有xyc=pcos0+psin0TOC o 1-5 h zxxy(b)t=pcos0一psin0 xyx

11、y将式(b)代入式(a)，并注意到2cos20=1+cos20，2sin20=1-cos20，cos20-sin20=cos20和2sin0cos0=sin20，可得c+cc-cc=x斗+斗cos20+tsin20(2.3la) HYPERLINK l bookmark52 x22xy(2.31b)c-ct=-x斗sin20+tcos20 HYPERLINK l bookmark54 xy2xy兀将式(23-1a)中的0换成0+-，则得cyc+cxy2c-cxy2cos20-tsin20 xy(2.31c)图2.7四面体的应力分布如果BC面趋近于A点，且已知A点的应力分量c,c,t时，则由式(

12、2.3-1)xyxy可求得过该点任意方向的平面上的应力分量。因此，对于平面问题，式(2.3-1)描述了该点的应力分布规律，即描述了该点的应力状态。对于三向应力状态，可以采用类似于二维应力状态分析的方法。现在研究从受力物体中取出的任一无穷小的四面体(图2.7)。斜面ABC的法线N与坐标轴间的夹角的方向余弦分别是1、m、n。四面体棱边的长度分别dx、dy和dz。设斜面的面积为1,则三角形OBC、OAC、OAB的面积分别为1xcos(N,x)=11xcos(N,y)=m1xcos(N,z)=n如果ABC面上单位面积上的力为p，沿坐标轴方向的分量p,p,p可由傲小四面体单元xyz的平衡条件得到p=cl

13、+tm+tnxxxyxzp=tl+cm+tn(2.32)yyxyyzp=tl+tm+cnzzxzyz式(2.32)是与坐标轴呈任意倾斜面止单位面积上的面力，该式也可按下标记号法和求和约定缩写为pn(i,j=x,y,z)(2.33)iijj式中n为斜面ABC外法线n与(j二x,y,z)轴间夹角的方向余弦1、m、n。j为了分析一点处应力的某些特征，现将坐标系oxyz变换到新坐标系oxyz,且新坐标系的ox轴与图2.7中的法线方向n重合，新旧坐标系间的方向余弦cos(x,x)二1,cos(x,y)二m,cos(x,z)二n,，如表2.1所示，则x方向的正应力111c为xc=pl+pm+pnxx1y1

14、z1表2.1新旧坐标系间的方向余弦xyzx1mn1i1y1mn222z1mn333将(2.3-2)代入上式，并注意到1、m、n分别等于1,m,n，则得111c=c12+cm2+cn2+2(t1m+tmn+tn1)xx1y1z1xy11yz11xz11类似地将p,p,p在y,z方向投影，可得到xyz=c12+cm2+cn2+2(t1m+tmn+tn1)yx2y2z2xy22yz22xz22c=c12+cm2+cn2+2(t1m+tmn+tn1)zx3y3z3xy33yz33xz33txytyz=c11+cmm+cnn+t(/m+/m)+t(mn+mn)+t(/n+/n)x12y12z12xy12

15、21yz1221zx1221=c11+cmm+cnn+t(1m+1m)+t(mn+mn)+t(1n+1n)x23y23z23xy2332yz2332zx2332第二章应力状态第二章应力状态 Tzx=cll+cmm+cnn+t(lm+Im)+t(mn+mn)+t(ln+In)x31y31z31xy3113yz3113zx3113采用张量的方法，可将以上各式统一表示为c=llc(2.3-4)IIII77ijiijjij式(2.3-4)则是c在坐标变换时所遵循的法则。凡是一组9个量c，在坐标变换ijij时遵从式(2.3-4)的法则就称为二阶张量。2.4边界条件当物体处于平衡状态时，除物体内部各点要满

16、足平衡微分方程式(2.2-4)外还应满走解条件。定解条件一般包括初始条件、边界条件或其它能确定唯一解答的补充条件。对于弹塑性静力学问题，定解条件主要是边界条件，所以弹塑性力学问题也就是数学物理方程中的边值问题。其它如约束条件、位移单值条件等也是常遇到的定解条件。在弹塑性力学中，给定面力的边界，用S表示，结定位移的过界，用S表示,cu如图2.8所示。本节主要讨论弹塑性力学平面问题的边界条件。a)b)图2.8平面问题边界条件位移边界条件所谓位移边界条件，就是在给定位移的边界上，物体的位移分量必须等于边界上的已知位移。设平面弹塑性体在S边界上给定X、y方向上的位移分别为U和V；，它们是u边界坐标的已

17、知函数；而位移分量u、v则是坐标的待求函数。当把它们代入S边u界的坐标时，则必等于该点所给定的位移，即u=u，v=V在S(2.4-1)u对于三维问题，在S边界的位移边界条件为uu=u(2.4-2)ii此处i=(x,y,z)，且对应于u、v、w。应力边界条件弹塑性体在外力作用下，处于平衡状态的条件，除物体内部各点的应力分量应满足平衡方程式（2.2-4）外，物体边界上各点也必须都是平衡的。由后者将导出应力边界条件。所谓应力边界条件就是在给定面力S的边界上应力分量与面力分量之间的关系。实质上，它是弹塑性体内部各点的平衡条件在其边界上的延续。因此，应力边界条件就是物体边界上点的平衡条件。设平面弹性体在

18、S上给定面力X、Y_，它们是边界坐标的已知函数；而应力c分量c、c、c则是坐标的待求函数。它们之间的关系可由边界上微元体的平xyz衡条件求出。不失一般性，在物体的边界上取一微元体（一般取为三角形微元，因为它可以描述任意曲线边界）.如图2.8b）所示，它在平面问题中显然是三角板（平面应力）或三棱柱（平面应变）。若令微元体边界面外法线N与x轴和y轴夹角的方向余弦分别为cos（N,x）=l,cos（N,y）=m；斜边长为ds，两直角边长分别为dx和dy，微元体的厚度仍取为1,则由图2.18b），根据力的平衡条件有、c/+tm=Xxxy卜（2.43）tl+cm=Yxyy如当边界平行于x轴时，有l=0,

19、m=1。这时，式（2.2-7）则为TOC o 1-5 h zc=Y，t=X（在S边界上）（a）yxyc而当边界平行y轴时，有l=1,m=0。这时，式（2.2-7）则为c=X,t=Y（在S边界上）（b）xxyc由此可见，当物体的边界线与某一坐标轴平行（或垂直）时，应力边界条件变得十分简单，即应力分量的边界值就等于对应的面力分量，应力分量的符号取决于边界面的外法线方向。当边界面的外法线方向与坐标正向一致时，等式右边取正号，否则取负号。但应注意，面力本身还有正负号。其规定与应力符号法则相同。对于三维问题，由力的平衡条件可得cl+tm+tn=Xxxyxztl+cm+tn=Y（2.4-4）xyyyztl

20、+tm+cn=Zxzyzz需要指出的是：在垂直X轴的边界面上，应力边界条件中不出现G，而在垂y直y轴的边界上不出现b。当作用在边界面上的面力不连续时，应分段或展开成x级数写出其边界条件；没有给定位移的自由边界，实际上是给定面力为零的应力边界，不能遗漏。3.混合边界条件在一般情况下，若用S表示整个物体的表面积，则往往在其中一部分面积S上b给出了面力，而在另一部分面积S上给定的是位移。如图2.9所示悬臂梁，固定u端部分属于S部分，它给定位移而末给定外力；其余边界均属S部分，它的外力ub已给定(包括外力等于零的部分)。显然，在S上各点应满足位u移边界条件式(2.4-1)，在S上各点应满足b应力边界条

21、件式(2.4-3)。对于混合边界条件，可以分别给在边界面的不同区域上，也可以给在同一区域的不同方向上。也即，对于边界上的一个点，在某一确定方向上，必须且只能给出S和S中ub的一种，既不能同时给定，也不能同时不给定；而同点在两个互相垂直方向止，可以是其中一个为S，另一个为S。bu例21如图2.9所示的一矩形截面悬臂梁，跨度为l，梁上表面作用均匀载荷q。试写出该问题的边界条件。并检查材料力学的应力公式是否满足力的边界条件。解：由材料力学所得的应力分量为bqx3y0qxsb=0,t=z(a)x2IyxyIzzh1)梁的上表面y=处2X=0，Y=q而/=cos(N,x)=0,m=cos(N,y)=一1

22、代入力的边界条件(2.4-3)，则解得t=0，b=qyxy由上式可知，因为材料力学作了纵向纤维无挤压的假设，无法算出b的分布规y律。因此，材料力学的应力计算公式(a)结果并不满足上表面b=q的边界条件。hy2)梁的下表面y二-处2X=0，l=cos(N,x)=-1m=cos(N,y)=0代入式(2.4-3)后解得T=0yxG=0y由上式可见，材料力学的应力计算公式(a)的结果满足该边界的力边界条件,其中G=0是由材料力学的假设得出的。y3)x=0的自由端处X=0，Y=0又l=cos(N,x)=-1，代入式(2.4-3)后解得m=cos(N,y)=0G=0y因此，在该边材料力学的应力计算公式(a

23、)的结果也满足该边界的力边界条件。4)x=l的固定端处因为固定端的外力分布没有具体给定我们只能求出该端面上的合力和合力矩的大小。且固定端限制了梁的移动和转动，所以该截面的位移边界条件是很重要的。位移边界条件可表示为T=0,xyu=0,v=0,色=0或竺=0dxdy(在x=l，y=0处)有关这方面的内容和处理方法将在后面的章节中详细介绍。2.5主应力、主切应力和八面体应力在受力物体内一点任意方向的微小面元上，一般都有正应又和切应力，不同方向的面元上这些应力有不同的数值。当此微小面元转动时，它的法线方向N随之改变，面元上的正应力g和切应力T的方向和它们的值也都要发生变化。在外法线方向不断改变过程中

24、，必然会出现面元上只有正应力，而切应力T等于零的情况。把这时面元的法线方向N称为主应力方向(主方向)，相应的正应力G称为主应N力，它所在的面称为主平面。以下将说明，物体中任一点都有3个主应力和相应的3个主方向。1.主应力在图2.7中，如令p,p,p为ABC面上单位面积面力的三个分量，则有xyzP2=p2+p2+p2(a)xyz将面元ABC上单位面积的三个分量p,p,p投影到面元的法线方向N,即得面元xyzABC的正应力为TOC o 1-5 h zb二pl+pm+pn(b)Nxyz将(2.3-2)式代入(b)式，并经整理后则得b=b12+bm2+bn2+2(tIm+Tmn+Tnl)(2.51)N

25、xyzxyyzzx式(2.5-1)即为任意法线方向N的斜面上正应力的表达式。该面上的切应力为T2=p2b2(2.5-2)NN将式(a)和式(2.5-1)代入上式(2.5-2)，可得法线方向为N的斜面上的切应力。注意到/2+m2+n2=1(2.5-3)因而三个方向余弦并不是独立的。现以l、m为独立变量，b和n看成是/和m的N函数，并求(2.5-1)式的极值。因此，其一阶偏导数应满足二=0dm即bl+tm+tn+(tl+tm+bn)dnjdl=0、(c)(d)(e)xxyzxzxyzz,tl+bm+tn+(tl+tm+bn)dn)dm=0 xyyyzzxzyz丿由式(2.5-3)可求得n对l和m的

26、两个偏导数为dndl=一l.n，dn.dm=一m.n将(d)式代入(c)式，并注意到(2.3-2)式，可得p.l=pm=pnTOC o 1-5 h zxyz令其比值为b，则有Np=bl、xNp=bmyNp=bnzN式(e)说明，在正应力取极值的斜平面上，全应力投影与斜平面的方向余弦成正比,比值b当然是正应力，正应力投影就是斜平面上全部应力的投影，而切应力不存N在，因此主应力(主平面)确实存在。将(2.3-2)式代入(e)式，经整理后得(2.5-4)(2.5-5)Q-c)l+tm+tn=0 xNxyxztl+Q-c)m+tn=0 xyyNyztl+tm+(c-c)n=0 xzyzzN或用张量符号

27、写为定义为6J1i=jij0i丰j(c-6c)l=0ijijNj此处6为Ktonecker-6,ij在(2.5-4)式中，共有4个未知数卩1、m、n和c，但由(2.5-3)式知，1、Nm、n这3个方向余弦不可能同时为零，因此，(2.5-4)可看成是关于1、m、n的线性齐次方程组，而且应有非零解存在，由线性齐次方程组有非零解的条件可得到c-cttxNxyxz=0tc-ctxyyNyzttc-cxzyzzN展开上式得c3-Ic2+Ic-I=0(2.5-6)N1N2N3其中I=c+c+c1xyzI=cc+cc+cc-(t2+t2+t2)2xyyzzxxyyzxzcttxxyxz2ttI=tct=cc

28、c+t-(ct2+ct2+ct2)3xyyyzxyzxyyzxzxyzyxzzxyttcxzyzz方程式(2.5-6)是一关于c的三次方程，它至少有一个实根。令其为c=c，N3z该上t=t=0。这样式(2.5-6)中剩下的应力分量只有c,c,t，可由平面应yzxzxyxy力状态理论求得其余两主应力c、c以12及它们作用的方向。这就简单地证明了，在物体内的任意一点，一定存在三个互相垂直的应力主平面，以及对应的三个主应力，它们的方向称为应力主方向。因为主应力c，c，c是方程(2.5-6)的根，按大小排列为ccc，它123123们分别位于三个互相垂直的主平面，且在主平面上切应力为零，所以式(2.5-

29、6)也可改写为C3Q+b+Q)Q2+Qq)(5=0N123N122331N123由代数学可知，为保证此方程和式(2.5-6)的解相同，其系数应相同，出此可得三个系数为I=5+5+5=5+5+51xyz123I=55+55+55(T2+T2+T2)=55+55+552xyyzzxxyyzxz1223315TTxxyxz+2TI=T5T=555TT(5T2+5T2+5T2)=5553xyyyzXyzxyyzxzxyzyxzzxy12TT5xzyzz由于在一定的应力状态下，物体内任一点的主应力不会随坐标系的改变而改变，所以式(2.5-6)所给出的系数I，I，I分别称为第一、第三、第三应力张量123不

30、变量，简称应力不变量。以主应力5，5，5的方向为坐标轴(分别记为1、2、3)的几何空间，称为123主向空间。在主向空间，(2.5-1)和(2.5-2)式则为5=5l2+5m2+5n2(2.5-7)N123t=C212+q2m2+52n2(q12+qm2+5n2)2(2.58)N1231232.主切应力当在主向空间讨论切应力t的变化时，(2.5-2)式可写为Nt2=C212+C2m2+C2n2(C12+Cm2+Cn2)2(2.5-9)N123123由(2.5-3)可知n2=112m2将n2用上式代替后，(2.5-9)式可得t2=(C2C2)12+(C2C2)m2+C2(CC)12+(CC)m2+

31、CTOC o 1-5 h zN1323313233为了求出T的极值，取t2对1和m的偏导数，并令它等于零，这时有NN1Q5)12+Qc)m2一一Q5)=0323213mQ5)12+(55)m2一一(55)=01323223满足上式的解有以下四种情况：I=0、m=0，由(2.5-3)式可得n=，由(2.5-7)式得工二0，这是N一主平面。l丰0、m=0，由式(f)的第一式得(Q-Q)(1212)=013因(Q-Q)丰0，故13l=由式(2.5-3)可知n=112并平分q、Q所夹再的平面，如图2.10a)所示。13a)该解表示通过Q,2图2.10主切应力平面用同样的方法可得(3)l=0，m=n=42(4)n=0，/=m=-J2解(3)代表通过Q，并平分q、Q所夹角的平面，见图2.10b)；而解123代表通过q并平2分Q、Q所夹角的平面，见图1.10c)。现将所有的解列于表2.2中。13表2.2切应力工有极值的平面方位l00101p21/J2m0101/V201迈n1001/J21迈0将以上所得到的I、m、n值代