弹塑性力学-第3章应变状态

第三章应变状态理论在外力、温度变化或其他因素作用下，物体内部各质点将产生位置的变化，即发生位移。如果物体内各点发生位移后仍保持各质点间初始状态的相对位置，则物体实际上只发生了刚体平移和转动，这种位移称为刚体位移。如果物体各质点发生位移后改变了各点间初始状态的相对位置，则物体同时也产生了形状的变化，其中包括体积改变和形状畸变，物体的这种变化称为物体的变形运动或简称为变形，它包括微元体的纯变形和整体运动。应变状态理论就是研究物变形后的几何特性。即给定物体内各点变形前后的位置，确定无限接近的任意两点之间所连矢量因物体变形所引起剧烈变化。这是一个单纯的几何问题，并不涉及物体变形的原因，也就是说并不涉及物体抵抗变形的物理规律。本章主要从物体变形前后的几何变化论述物体内一点的应变状态。位移与线元长度、方向的变化坐标与位移设变形前物体上各点的位置在笛卡尔坐标（Descartercoordinate）系的轴（X、、Y、Z）上的投影为（x,y,z），又设物体上各点得到一位移，并在同一坐标轴上的投影为（u、八w），这些位移分量可看作是坐标（x,y,z）的函数。于是物体上任点的最终位置由下述坐标值决定。即己=x+u（x,y,z）、n=y+v（x,y,z）＞c=z+w（x,y,z）,上式中函数u、v、w以及它们对坐标（x,y,z）的偏导数假设是连续的，则式确定了变量（x,y,z）与&,n工）之间的关系。因为物体中变形前各点对应看变形后的各点，因此式是单值的，所以式可看成是坐标的一个变换。如果在中，假设x=x,y=y，则由式可得如下三个方程00己=X+u((x,y,z)、0 00n=y+v(x,y,z)10 00£=z+同X,y,z)00式决定了一条曲线，曲线上各点M*,M*,A,在物体变形前为平行于z轴的直线12(X=X,y=y)上(图。由此可见，变形前物体上与坐标轴平行的坐标线，在变形00后的物体上一般将成为曲线。换句话说，如果用没有变形状态的坐标(X,y,z)末表征物体上各点的位置，到变形终了状态将是曲线坐标；反之，如果用&,n工)表示各点的坐标，则对巳变形物体是笛卡尔坐标，而对于变形前的物体将是曲线坐标。由以上可见，描述连续介质变形的方法有上述两种，分别称为Lagrange法Euler法。Lagrange描述法是用变形前的坐标(x,y,z)做自变量，而Euler法则是用变形后的坐标色内工)做自变量。在固体力学中，通常物体的初始形状、固定情况以及载荷是一定的，需要确定的是物体各点的位移u、v、放和应力o。对于小变形一般采用Lagrange坐标法；而对于大变形有时用Euler法。在数值计算中，通常采用矢量来表示，因为要计算变形前后两次应变的变化，所以用Euler法比较方便。在以后的讨论中，我们采用Lagrange坐标法。 图变形表示法变形体的应变设物体中变形前相距十分近的两点M,N，变形后移位至M*,N*。变形前M,N的坐标分别为M(x,y,z)，N(x+dX,y+dy,z+dz)，变形后M*,N*的坐标分别M*化,n<),N*化+dW,n+dn<+dq)。那么，矢量mn所表示的线元在物体变形后由矢量M*N*表示线元。那么，MN和M*N*的平方为MNF=dS2=dx2+dy2+dz2 (a)(b)根据式，点(b)根据式，点N*在％=dS2=d己2+d月2+dC2

*(c)此处du是因M,N两点所产生的增量，将其在(x,y,z)处展开为Taylor级数，即(d),du ,du ,du ,d2u d2u d2u/(d)du--dx+—dy+——dz+ (dx)2+ (dy)2+ (dz)2+ATOC\o"1-5"\h\zdxdydzdx2 dy2 dz2略去(d)式中的高阶微量(dx)2，…，并将(d)式代入(。)式，则可得eur.duI,duIC+dQ-x+dx+u+——dx+——dy+——dzdx dy dz由式知，Q-x+u，所以du du dudQ-(1+——)dx+——dy+——dzdx dy dz同理可得dv dv dvdn---dx+(1+—)dy+—dx dy dz"d*,dw dwdC-——dx+——dy+(1+——)dzdx dy dz式表示用物体的任意线元在变形前的投影表出它在变形后的投影。我们的目的是为了计算dS与dS之差，于是由(a)式和(e)式可得*dS2一dS2-2(8dx2+8dy2+8dz2+丫dxdy+ydydz+Ydzdx) (f)* x y z xy yz zx式中==--8x8y8zdu11,d==--8x8y8zdu11,du. ,dv. ,dw.-——+-I(——)2+(—)2+(——)2

dx21dx dx dxdv11zdu、dv/dw、1=——+，-I(——)2+(—)2+(——)21dy2Ldydy7dy，」dw1|\du.dv/dw.]=—+-I(—)2+(—)2+(—)21dz2Ldzdzdz_|dvdududvdvdwdw+ + + + 28dxdxdydxdydxdy xydwdududvdvdwdw+ + + + 28dydydzdydzdydz yzdwdudu+ + dx dxdzdvdvdwdw+ + dxdzdxdzzx>-28yx-28zy-28xz式实际上就是应变在各坐标方向的分量，它是非线性的。如果知道了变形体各点的位移u、V、攻,则可由该式求得各点的应变分量，式可采用张量表示为ij8x12712ij8x127127yxzx8y1万丫2zy12丫127xzryz线元的长度变化引入符号EMNdS—dS

^

EMNE是点M和N间由变形引起的距离的增加量对二者间变形前的距离的匕我MN们把这个量称作点M在点N方向的相对伸长度。根据式(a)和式(f)，并注意式，则可得伸长蜃的表达式为TOC\o"1-5"\h\z1 MNE(1+E)dS2=8dx2+8dy2+8dz2+7dxdy+7dydz+7dzdxMN2MN x y z xy yz zx=8l2+8m2+8n2+7lm+7mn+7nl

x y z xy yz zx式中l=竺，m=dy，n=喧是矢量MN的方向余弦。如果在（g）式中令dSdSdSl=1,m=n=0，可口么有E=,：1+28—1x x此处E表示M点在x方向的相对伸长度。类似有M点在y、z方向的相对伸x长度为E=1+28—1 E='1+28—1y' y z z因此，应变分量8、8、8描述了变形前平行于坐标轴的那些线元的伸长度，xyz它们称为正应变。线元方向的变化变形物体中的线段，在变形时不仅长度要改变，而且方向也会发生变化。矢量mN与坐标轴(X，Y，Z)形成的方向余弦分别为l、m、n；而矢量M*N：与坐标轴夹角的方向余弦分别为

d己

dsd己

dsd叩m*- dSd。

dSTOC\o"1-5"\h\z利用式解得dS=(1+E)dS，并注意到式可得* MN11+EMN11+EMN1MNn* 1+EMN(1+一)l+一m+一nd.x Sy dzdv Sv Sv—I+(1+——)m+—nSx Sy SzSw, Sw / S叭——l+——m+(1+——)nSx Sy Sz式表示任意线元在变形后的方向，即变形后M*N*的方向余弦可以用变形前MN的方向余弦表示。如果变形前线元dx与乂轴平行，则该线元的方向余弦为l-1，m-n-0，那么由式知，该线元变形后的方向余弦为rSu Sv Sw1+—— — ——s Sx e Sx “ Sxl*-m*-n*-1+E 1+E 1+Ex x x此处E是变形前与乂轴平行线元的伸长度。由上式可以看出，对于任意线元，因x各个方向的位移u、v、w不相同，因此方向要改变(图；同时各个方向的伸长度也不相同，方向也要改变。因为线元dx在变形后成为已变形物体上坐标曲线自上的线元，所以式实际上给出了点M*上坐标曲线己的切线方向的方向余弦。类似地可以由式得出已变形物体上坐标曲线y和z的切线的方向余弦。如果用i、i、i表示点M*在坐标白、小。xyz切线方向的三个单位矢量，那么该三个单位矢 图线元的方向余弦量相对于笛卡尔坐标的方向余弦可由式如同线元dx那样得到类似的式。具体列于表。类似于的方法也导出用m*N*的方向余弦表示变形前mN的方向余弦，读者可自行推导。表变形后相对于笛卡尔坐标的方向余弦

剪切度与切应变如图所示，设变形前物体中经过M点的两条任意纤维I和H，此两纤维在M点的切线的方向余弦分别为l、m、n和l、1 112m、n；变形后，物体中的M点移动到M*，22纤维I和II变成纤维I*和II*，纤维I*和II*的方向余弦也变为l*、m*、n*和l*、m*、n*。TOC\o"1-5"\h\z1 112 22由前面可知，变形后两纤维的方向余弦可用变形前的方向余弦表示，同时由解析几何知图剪切变形cos(I*,II*)=l*l*+m*m*+n*n*12 1 2 12图剪切变形则可求得变形后纤维I*和II*之间夹角的方向余弦。将式代入上式，并注意式，则可得cos(I*,II*)= 1(1+2e)ll+(1+2e)mm+(1+2e)nn(1+E)(1+E) /12、 /12 J12I II+y(ml+1m)+y(nm+nm)+y(nl+nl)]

xy12 1 2yz2 1 1 2zx12 21注意，式中纤维I和II的伸长度E和E由确定，但必须用变形前物体的纤维 I和IIIII的方向余弦l、m、n和l、m、n。1 112 22由显然可知，当知道了6个应变分量8、£、8、Y、Y、Y和变形前经xyzxyyzzx过物体中任意一点处的两纤维的方向余弦后，则可由式和求得该两纤维变形后的夹角。如果变形前物体中纤维和II分别平行于X轴和丫轴，则l=m=1,其余的12方向余弦为m=n=l=n=0，且变形后物体中纤维I*和II*的切线方向分别与1122

单位矢量1•、i重合，则根据式和表达式可知xyTOC\o"1-5"\h\z8 8cos(i,i)= xy = xy —xy(1+E)(1+E) ,(1+28)(1+28)X yV x y在变形前，纤维I和n的夹角为直角，令9为变形后纤维引起的夹角减少量，xy8xy、.(8xy、.(1+28)(1+28)x y兀cos(i,i)=cosj—9)=sin①xy 2xy xy类似可得分别与KZ轴和Z、X轴平行的两纤维夹角的减少量p和①为yzzx8 、sinp= yz-yzJ(1+28)(1+28)yyz>8sinp=— ―zx；(1+28)(1+28)zzxJ称角p、p和p为剪切度。由以上分析可知，y、y和丫zx表示了切应变，当xyyzzx xyyz它们均为零时，则纤维之间的夹角变形后保持不变。由以上的分析可知，在式中，应变I篇出现的量不外乎下面9个分量j称u为相对位移张量。因此，称u为相对位移张量。因此，当知道了位移对坐标的偏导数，则可根据式计算出应变分量I｝，从而也就知道了任何一根线段在任何方向的伸长E(由式)，同时TOC\o"1-5"\h\zij MN还可计算出原来与坐标轴平行的两线段角度的减少量方，更一般地还可计算cos(I*,『)，因此8｝充分地表示了应变。如果I20：则意味着没有变形，仅有ij ( ij刚体移动或转动；如果已知应变分量七土则不能求得一根线段的绝对角度变化，ij因为这时并不知道中的任何值，所以也无法由式求得l*、l*、…等，反过来也不知道1 2l、l…等，所以cos(I*,II*)无法求出。1 23.2应变张量与转动张量一般来说，物体中各点的变形由式中Ij｝的6个分量可完全确定，因为知道了这6个分量就等于知道了伸长度和剪切度。。

在变形理论分析中，通常还需引入9个参数，即SueySvSySvSwezSwSzSwSue=一SvxSxSue=e +—e=e +一e=e + >xy yxSySxyz, zuSzSyzx xzSxSz1SwSv、1SuSw1SvSu、3=—(—)3=—(■———)3=-(——■—)x2SySzy2SzSxz2SxSy这样，位移的所有一阶偏导数都于由这式的个参数表示为SuSu1SuSu1——=e，，=—e—3Sx xSy2xy加1Sv—=—e+3—=eSx2xy zSy ydw1Sw1—=—e—3—=—e+3Sx2xz， ySy2yzzxeg&加在1=—e+32xz1=—e—3微元体的转动为了研究微元体的转动，首先阐述3,3,3的几何意义。为此，设垂直于Z轴xyz的线元为MN，利用，并注意此时dz=0，则式成为d工=(1+e)dx+(―e—3)dyx 2xyzd"=(—e+3)dx+(1+e)dy2xyz ydC=(—e—3)dx+(—e+3)dy2xzz 2yzx在图中表示出平面xoy，线段MN是变形前线元的长度，而线段M*N*是变1形后在平面XOY上的投影。根据图显然有tan9=dy，tan9*二坐dx d&再根据式和分子分母均除MN原长dS，则有(_Le+3)cos9+(1+e)sin9tan9*=2^y_z 1 y (1+e)cos9+(—e—3)sin9x 2xyz图线元角度的变化

变形物体在变形过程中，由前节已经知道，线元不仅产生尺度变化，而且线元的方向也发生变化。但是在变形时起变化的不仅线元的相对方向，而且还有它的绝对方向，因为从初始状态的物体中割离出来的无限小微元体，到终了状态时，除了产生形变外，还有—些转动。把这术语应用到微元体上（它在产生位移过程中不仅位置要发生改变，而且还改变了大小和形状），意指所有属于微元体的许多个线元转动的平均值。同时，约定w作为绕z轴的转动角，此处z轴是变形前和线元MN垂直的轴，W是线元MN（在变形前的位置）和它在变形后M*N；在垂直z轴平面上投影之间的夹角（图.因此MN因变形产生绕z轴的转动角为W=8*—8z由、和三角函数关系式可得tan8+tanw（2,十^z）cos8+（1+1）sin8 -二 1tantanwz（1+e）cos8+（L-3）sin8x2xyz从上式可解得tanwz图线元绕从上式可解得tanwz3+—ecos28+—C—eIin28TOC\o"1-5"\h\zz2xy 2yx1+ecos28+esin28+eesin28x y 2xy取W从8=0到8=2兀间隔中的平均值（即取所有垂直于z轴的线元的tanw的平zz均值），tanw=——Jtanwd8=H+H

z 2兀 z 1 2d8d8其中1+ecos28+esin28+—esin28

x y 2xyecos28+(e—e)sin28

xy y x

1+ecos28+esin28+—esin28x y 2xycos28+esin28+—esin28,则积分H可化为x y 2xy 2

Inf|2k0而积分H可化为132F d0z~1兀2+e+e+-(e一e)2+e2sin(20+a)

0 xy xy xy3a于兀3a于兀2tT」ady■2+e+e+■'(e-e)2+e2siny

xy xy xy(a)式中—arcsine-式中—arcsine-ee—arccos( xy]j'e-e力+e2Vvx y xy由以上可得tanvz1 3 tanvz1 3 . z X2兀： 1■1+e+e+ee-e2xyxy4xyarctan入 y- (2+e+e)tan—+■'(e-e)+e2xy2 *xyxy1一2;1+e+e+ee——e2xyxy4xy2兀+a(b)因为所求得的arctan()且多值函数，所以结果不定，但这种不定性可以揭示出来。如果考虑到使e,e,e趋近于零时，从(a)式可知积分H必然趋近于3，据此在xyz 1 z(b)式中必须使(2+earctan (2+earctan ：y .一+e)tan—+：(ey2 %x一e)+e2

yxyT——e24xy2兀+a(c)于是得到而口的表达式为tanvtanvz1(1+e)(1+e)—e2

,xy4xy用类似的方法可以写出而正和而铲为x y

tanvtanvyTOC\o"1-5"\h\ztanvtanvy,,(1+e)(1+e)— e2।yz4yz3

y三个参数tanw,tanv，tanv表征出包围M点的无限小体积的转动，它们x y z分别和3,3,3成正比，而且在后者等于零时也等于零。因此，如果它们在任xyz一坐标系X,y,Z等于零，那么它们在任何其他坐标系中都等于零。因此可得出结论，在物体上任一点，如果有3=3=3=0xyz则表示通过该点的线元对于通过这点的任意轴平均来说没有转动。所以，等式是物体上M点周围任意无限小微元体没有转动的条件。应变张量与转动张量在直角坐标系下，式称为Green应变张量。虽然是在直角坐标系中导出的，但它们所描述的几何关系与坐标系的选择无关，因此适用于任意正交坐标系。从数学角度出发，Green应变张量属对称二阶张量。对于式，如果忽略高阶微量，则式中将成为due=—xd.x加e ysyswy=yy=ySuSv+ =2e=2exyyxSySxxyyxSvSwy=y+一=2e=2syzzySzSyyzzySuSwy=y+一=2e=2ezxxzSzSxzxxze= z Sz >8ij)J

四办加一很

++

加」ax8ij)J

四办加一很

++

加」ax史axeg一&史&而一小史氏1-21-2在直角坐标系中，称式为Cauchy应变张量，它也是二阶张量。由张量分析知，任何一个二阶张量都可唯一地分解成一个对称张量和一个反对称张量。因此式可写成1U=—(U+U)+—(U—U)=8+3i,j2i，jj，i2i，j j，i jij其中8的元素为式，而3的各元素为ijij1Sw1Sw(——2Sx666-61-21-2av-axeg一办o/(1-2SwSxSwSy根据式，上式也可写为TOC\o"1-5"\h\z0 —3 3zy3=3 0 —3ijz x—3 3 0y x由知，3反映了包围某点无限小微元体的转动，因此称3为转动张量。ijij物体变形的描述与简化以上的讨论阐述了已知位移U、v、w决定物体任意一点无限小区域的位移、转动和纯应变三个因素，这些因素确定了假想从物体中切割出来的无限小微元体受载后的终了位置和终了形状。但是必需指出，整体位移和微元体的转动并不是微元体变形的特征，变形是由应变8所决定的。但是，如果说整个物体的变形，则其具有特征的是物体上各ij点的位移和各纤维的转动角。如梁的变形通常是指梁的挠度（即位移），而轴的变形是指轴的一端相对于另一端的扭转角（即转动角）。因此从这个观点出发，位移和转动是整个物体变形的特征，而伸长度和剪切度是物体在无限小微元体变形的特征。这两个特征具有实际意义，前者决定承力构件或物体的刚度，后者决定承力构件或物体的强度。由上面的叙述可知，“小变形”这一术语就产生二种解析，“小形变”可以了解释为小伸长度和小剪切度（同1相比）；或者可以解释为小位移（同物体的尺度相比）和小转动角（同1相比）．必须指出，小变形的经典理论实际上是建立在小位移和小转动角的假设基础上的。但是这一情况很少用应有的明确程度加以说明。因此，习惯于将变形的概念联想到应变分量，以为是指小伸长度和小切应变。应孩强调的是关于小位移和小转动角的假投，比之小应变分量的假设，很大程度地限制了论述的普遍性、并且前一假设的结果也服从后一假设，但反过来则并不肯定。其次必须指出，在必须表明小位移时，常常没有预先说明它应和什么比较是微小的，其实这样的预先说明是十分必要的，因为位移是具有量纲的量。今后使用“小变形”这一术语时，始终是指小伸长度和小剪切度（同1相比）。此外，如果所研究的问题同时又有小位移和小转动，总是预先加以说明。直到现在，还没有对伸长度和剪切度的大小加以任何限制。但是在弹塑性力学的论述中这样的普遍性照例是没有必要的，因为只有很少的材料（例如橡皮）才在颇大的相对应变的情况下还保留它的弹性的性质，大多数应用在工程上的材科（例如所有的金属和合金）仅在同1相比很小的伸长度和剪切度下才是弹性的。例如钢的弹性变形区域大概在相对伸长度的数值是10-3到5x10一3这一量级，钢的弹性切应变的最大数值也在类似的量级上。有色金属（及其合金），数字稍有不一样，但它们的弹性变形区域也限制在很小的伸长度和切应变内。从上面可见，弹塑性力学应用于金属结构时，略去同1比较起来很小的伸长度和剪切度以简化公式是很自然的而且是合理的。由此，小应变理论提供了最大的实际兴趣。当忽略同1相比很小的伸长度和剪切度来简化前面己求得的公式，则对于实行简化，得到：E=£,E=8,E=8xx yy zy而根据式得到①2£, ①Eg,平28xyxy yzyz zxzx因此在小应变的情况下，应变分量8,8,8可以与相应的伸长度等同看待，而应变xyz分量8=1y,8=-y,8=-y可以与相应的剪切度等同但因它们仍需用式xy2xyyz2yzzx2zx计算，所以小应变仍属非线性。应当指出，如果转动很大，而剪切度却很小，那么在决定变形后线元的方向时，同转动相比可以略去剪切度（这里是指可一般略去剪切度；在有剪切度但不同

时有转动的式子中，就不能略去剪切度）。材料力学梁的理论中的平面假设、和板理论中的Kirchhoff假设，就是这种简化处理可变形物体线元方向的例子，两者都是是假定同转动相比可以略去剪切度而作出结论的。如果应变和转动角都很小，此时同1相比微小的不仅仅是应变分量，而且在中还可略去与转动角相关的平方项高，从而可以获得tanv»V23，tanv邓-3 tanv川-3xxx yyy zzz当伸长度、剪切度和转动角同1相比都很小时，利用式可将改写为e2+(，e+3)2+(1e—3)2x2xyz 2xzye2+(1e+3)2+(1e—3)2y2yzx 2xyze2+(1e+3)2+(1e—3)

z2xzy 2yzx8xy8yz8zx=e8xy8yz8zx=exy=eyz=ezxxy3+yzxze1-21-21-2yzxeee+++)))zxy333———xyyzxzeee1-21-21-2xyzeee++++3)+z+3)+—3)(-ey2k、/1—3)(ez2xz1—3)Gex2xy+3)=-Yx2xy1+3)=Yy2yz1+3)=-Yz2zx在应变分量公式中包括了以下项：（a）参数e的线性项；（b）参数e彼此相乘ij ij的项；（c3彼此相乘的项和（d）参数e与3的乘积项。k ijk当e和3同1相比很小时，则有两种可能性，则有两种可能性：（1）3与eijk kij同阶或更高阶的微量；（2）e与32同阶或更高阶的微量。ijk对于第一种情况，式中只需保留线性项，因此应变分量用式计算，即为Cauchy张量。对于第二种情况，在式中只需保留（a）和（c）形式的项，简化后可得应变分8-8-e—33xyxyxy8-e—33>yz yz y z8-e—33zx zx z x8-e+—(32+32),xx2yz1 18-e+—(32+32),yy2xz18-e+—(32+32),zz2xy式和式应用很广，式用于小应变分量和小转动分量，而且两者属同一量级时，则就是通常指的小变形情况。式用于小应变和小转动，但转动仍比应变大很多，因此适合柔性构件问题，如细长杆、板壳等，这种情况通常称为大变形小应变。

而式属大变形大应变问题。总之，Cauchy应变张量属于线性问题，其余均为非线性问题。主应变和应变不变量应变张量的坐标变换同一个变形可在不同的坐标系中研讨。在所有各种情况下，可以用前面所确定的六个应变分量把变形的特征充分地表示出来，但这六个应变分量的值却随坐标轴方向的选择而变更。设原有的坐标系为X、y、Z，另一坐标系为X'、F、Z、，它的各轴的方向对第一个坐标系各轴的方向余弦如表所示。表新旧坐标系之间的方向余弦表新旧坐标系之间的方向余弦因为二个坐标系均为直角坐标，12因为二个坐标系均为直角坐标，12+12+12=1123m2+m2+m2=1123n2+n2+n2=1123上式也可写为12+m2+n2=11 1112+m2+n2=12 2212+m2+n2=13 33因此表所列方向余弦之间有下面的关系：、11+mm+nn=02 12 12+mm+nn=0A13 1 3 1311+mm+nn=023 2 3 231m+1m+1m=011 22 331n+1n+1n=0 >11 22 33mn+mn+mn=011 22 33如果线段在第二个坐标系X'、y‛、Z，各轴上的投影是dx‛，dy，dz，，那么在第一个坐标系x、y、z各轴上的投影是：dx=ldx'+1dy'+1dz'123dy=mdx'+mdy'+mdz'>123dz=ndx'+ndy'+ndz'123注意节中的式(灯左边是表示点M和N之间距离的平方因变形而引起的变化，由于这两点的选择与坐标无关，该式左边也应与坐标选择无关，因此在坐标变换过程是应是不变量。于是将式右边的dx,dy,dz用矢量MN在新坐标上的投影dx'，dy'，dz'的代入，将有E(1+-E)dS2=8'(dx')2+8，(dy)2+8、(dz')2+2(8'dx'dy'+8'dy'dz'+8'dzdx')TOC\o"1-5"\h\zMN2MN x y z xy yz zx其中8'=812+8m2+8n2+2(81m+8mn+81n)xx1y1z1 xy1 1yz11zx118' =812+8m2+8n2+2(8 1m +8 mn+8 1n)yx2y2z2 xy2 2yz22zx228' =812+8m2+8n2+2(8 1m +8 mn+8 1n)zx3y3z3 xy3 3yz33zx338' = (8 1 1 +8 mm +8 nn)+8 (1m+1m)+8 (mn+nm)+8(1n+n1)xyx12y1 2z12xy1 2 2 1yz12 1 2zx12 128' = (8 1 1 +8 mm +8 nn)+8(1m+1m)+8 (mn+nm)+8(1n+n1 )yzx23y2 3z23xy2 3 3 2yz23 2 3zx23 238' = (8 1 1 +8 mm +8 nn)+8 (1m+1m)+8 (mn+nm)+8(1n+n1)zxx13y1 3z13xy3 1 1 3yz31 3 1zx31 31由以上可见，式与式在形式上相类似，因此式中所含的各系数8',8',8',8',8',8'在坐标系X'、Y'、Z，中的意义与系数8,8,8,8,8,8在坐xyzxyyzzx xyzxyyzzx标系X、Y、Z的意义一样。显然，式就是坐标轴变换时应变分量的变换规律。主应变和应变不变量现在讨论在那一个方向伸长度E会具有极值。设取X，轴平行这个方向，那MN么根据式有1TOC\o"1-5"\h\zE(1+E)=8′或E=.，1+28‛-1x2x x x x由上式可见，求E的极值归结为求8‛的极值，也即要确定1,m,n的值，使x' x 111得在该方向上使式中的第一式有极值。由式于知，1,m,n之间存在如下关系1 1112+m2+n2-1=01 11那么，假设一函数为f=8'-8(12+m2+n2-1)x1 11式中8为Lagrange乘子。现将上式分别对l,m,n求偏导数，并使其等于零，则得1 1 1如下线性方程组(8-8)l+8m+8n=0兀 1xy1zx18l+(8-8)m+8n=0>xy1y1yz18l+8m+(8-8)n=0TOC\o"1-5"\h\zzx1yz1z 1J由于条件式的存在，l,m,n不可能同时为零，因此是关于l,m,n的线性齐次1 1 1 1 1 1方程组。根据齐次方程组有非零解的条件，式中的系数行列式必为零，即8-8 8 8x xy zx8 8-8 8 =0xy y yz8 8 8-8zx yz z它至少有一个实根，将它记为。注意到中8'的表达式还可写为1 x8'=(8 l+8m+8n )l +(8 l+8m +8 n)m +(8 l+8m +8 n)nxx1xy1zx11xy1y1yz1 1zx1yz1z1 1将上式括号中的式子用中的值代入，并注意到式，将发现‛=8，x1即8'的极值就是8。x1当分别设丫‛,z，平行于伸长度e具有极值的方向时，采用类似的方法可分别MN得到关于l,m,n和关于l,m,n的类似于式的线性齐次方程组，且其系数行列式2 2 2 3 3 3与式完全一样。将式展开该式得83-I82+I18-I1=01 2 3I'=8+8+81xyz其中{I1=88+88+88-(82+82+82)2 xyyzzxxyyzzxI'=888+2888-(882+882+882)3xyz xyyzzxxyzyzxzxy这三个参数I,I,I分别称为第一、第二、第三应变不变量。123TOC\o"1-5"\h\z方程式有三个实根。设这三个实根分别为、8、8。则由根与系数关系，有1 2 3I'=8+8+81123I'=88+88+88〉2 12 23 31I1=8883 123J称8、8、8为主应变，所在方向称之为主方向。1 2 3另外，注意到应变分量8'和8'可以写为xyyz8'=(8l+8m+8n)l+(8l+8m+8n)m+(8l+8m+8n)nxyx1xy1zx1 2xy1y1yz1 2zx1yz1z1 28'=(8l+8m+8n)l+(8l+8m+8n)m+(8l+8m+8n)nzxx1xy1zx1 3xy1y1yz1 3zx1yz1z1 3

将上面两式中括号内的式子用代入，则有8'=8(11+mm+nn)xy112 12128'=8(11+mm+nn)

zx113 1313由可知，8'=8'=0。因此，如果在X'轴方向的伸长度是极值，那末应变分量xyzx8'=8'=0，也就是变形发生时，在X，方向和F方向之间的直角以及X，方向和Z、xyzx方向之间的直角没有变化。由此可见，不论在物体上任何点的变形怎样，总可以找出通过物体的三条纤维，它们在变形前是互相垂直的，而在变形后仍然还是互相垂直。TOC\o"1-5"\h\z将8代入式，则可求得1,m,n，从而可确定8的方向，即主方向。如果将式中1 1 1 1 1的1,m,n分别用1,m,n和1,m,n代替，以及别将8用8和8代替，则可求得1 1 1 2 2 2 3 3 3 2 3,m,n和1,m,n，从而确定主应力8和8的主方向。2 2 3 3 3 2 3完全类似地还可求得最大切应变为Y=±(8-8)]1 2 3Y=±(8-8)\3 1y=±(8-8)1 2 1以及八面体的切应变为=2(8—3 121(8—8)2+(83xyy—8)2+=2(8—3 121(8—8)2+(83xyy—8)2+(8z z—8)2+6(82+82+82)-T

x xyyzzx应变偏量及其不变量分别为TOC\o"1-5"\h\z」(28—8—8) 0 03 12 31e= 0 -(28—8—8) 0j 3 2 3 1八 八 1入0 0 -(28—8—8)3 3 12J1=0J'=ee+ee+ee>2 12 23 31J'=eee3 123J3.4应变率张量和应变增量张量应变率张量

在小变形条件下，应变张量可简写为。1..、£一一（U+U）ij2 i，j ji而当介质处在运动状态时，以V（九y,z,t）表示质点的速度，V表示速度的三个分量，i以时间t作为起点，则经过无限小时间段dt以后，位移为U=Vdt，由于dt很小，iiU及其对坐标的导数也很小，因此可以应用小变形公式，即£ij£iji,jj,i=—(U+U)=—(V+i,jj,i如果令J-..dt=£，则有jijTOC\o"1-5"\h\z. 1/ 1、£=—(V+V)ij2 i,j j,i£..称为应变率张量，上式定义不论£..大小都成立，但要求是对每一瞬时状态进行j j计算，不是按初始位置计算。因为，在一般情况下当按初始位置计算时.d£中——£ijdtij只有在小变形条件下才有d=d=—£dtija=—£atij由和式可知£ijd=£ijd=—£dtijdt_1/•・、——(u+u)2 1，j j，1于是应变对时间的变化率为au・aVau£=——y=——+——xat xy ax ayaV ・awaV£=—— y=——+——>y ay yz ay az・・・

・ aw ・ au aw£= y= + TOC\o"1-5"\h\zz az zx az ax将上式写为张量形式为ij£x1.5ij£x1.5y2xy1・5y2xz—y y2 xy 2 xz・ 1・£ 万yy 2 yz1・ ・5y£2yzz

应变增量张量应当指出，对于固体材料，当温度不变时或变形是缓慢的，则其力学行为与应变率关系不大，只有在受到动载荷时，因变形速率很快，材料的力学性质才会与应变速率有关，这类材料通常称为应变率敏感材料。因此，根据第一章中的基本假设，时间因素对物体的弹塑性力学行为不发生影响（即不考虑粘性效应），而且这里的dt并不代表真实的时间，仅仅代表加载变形的过程。于是，对于这里所讨论的问题主要关心的不是应变速率.而是应变增量d8j。于是采用应变增量d8j代替应变率8・更能表示不受时间参数选择的特点。 " "ij以du代表位移增量，则式成为id8d8ij=—(du+du)2 i,j j,id8d8ij=—(du+du)2 i,j j,i=d—(u+u)_2i，j j，i=d(8)

ij这说明在小变形时、按瞬时状态计算d8与按初始状态计算d8（近似地）没有什ijijij么区别。类似地、应变增量张量的应变增量偏量为de=d8-d86

ij ij ij注意，在求应变增量时，每一次都应从瞬时位置计起，而不是从初始位置算起。例如在简单拉伸时，轴向应变增量为dld8=-l此处l是拉伸时的瞬时长度（为了不与8相混淆，令dE=d）。一般情况下应变增l量的累计值Jd8的物理意义并不明显，但是当应变张量的主方向不变时，它们的ij积分才有明确的物理意义。对于简单拉伸问题有E=JdE=ln（1+8）这就是对数应变，又称为真应变。

小变形的应变协调方程对于一个连续的物体，按某一应变状态变形后必须既不出现开裂，又不会出现重叠，即保持其连续性。此时所给定的应变状态是协调的，否则是不协调的。这就要求位移函数u在所定义的域内为单值连续函数。一旦出现了开裂，位移函数就会出现间断；出现了重叠，位移函数就不可能为单值。因此、为保持物体变形后的连续性，各应变分量之间．必须有一定的关系。在小变形情况下，6个应变分量是通过6个几何方程与3个位移函数相联系。TOC\o"1-5"\h\z若已知位移分量u，则由求得各应变分量。若给定一组应变8，式是关于未知位i ij移函数u的微分方程组，它包含6个方程，仅三个未知函数，方程的个数超过了未知数的个数，若任意给定8，则方程不一定有解，仅当8满足某种可积条件，ij ij或称为应变协调关系时，才能由方程积分得到单值连续的位移场u。在小变形条件下，应变的计算式为。将式中的6个应变分量分为两组。第一组为式中的前三式，将该式中前两式分别对y和X求二阶偏导数，得d28 d3ud.xdy2将上面两式相加，得S28 将上面两式相加，得S28 S28 S2Su SvSy2 Sx2 (——十—)=SxSySySx这就是我们需要的应变之间的一个关系式．将上式内各字母循环替换，就得到另外二式，第一组中共有三个关系式。Sy S2攻S2Sy S2攻S2u小= + ， y^= +SS2vS小= + ， y^= +SzSxSzSySz SxSySxSzSx SySxSySySzSSy SSy Sy Sy S3w S28—(——y^+——田———x^)=2 =2SzSx Sy Sz SxSySz将上式各字母循环替换，就得到另外二式。第二组也共有三个关系式。于是第一组和第二组的6个关系如下TOC\o"1-5"\h\za28 a2e a2 + = yay2 ax2 axayxya28 a28a2—y-+ z= yaz2 ay2 ayaz yza28a28a2—z+—= yax2 az2 azax zxaay ay ay a28一(——z^+—————y^)=2*-axay az ax ayaza ay ay ay a28一(——iy+——y^———z^)=2 ay az ax ay axazaay ay ay a28—(——yz+——z^———xy)=2—zazax ay az axay上列应变分量之间的6个微分关系式，称为应变协调方程，又称变形连续方程（圣维南恒等式）。当弹塑性变形固体在外界因素影响下，物体中产生应力与应变，如能先求得位移u，匕攻，对于小变形问题则可由式可计算应变分量，这时应变协调方程自然满足，因应变协调方程本是由式所导得。但是，如先求出应力，然后再求应变，则所求的应变分量必须同时满足应变协调方程，否则，应变分量之间可能互不相容，因此也就不能用式求得正确的位移。式可以视为位移分量u,v,w的微分方程，如果应变分量和转动分已知，则求式的积分，就可求得位移分星u,v,w，进一步可证明，在求上述积分时、必须满足应变协调方程。正交曲线坐标中的应变几何方程在求解具有曲线或曲面边界的弹塑性力学问题时，一般选用正交曲线坐标系比笛卡儿直角坐标系更为方便。正交曲线坐标设以三个独立变量（x,y,z）来定义的三个独立标量函数（X,RZ）为X=X（x,y,z）, Y=Y（x,y,z）, z=Z（x,y,z）如果（X,Y,Z）表示笛卡尔坐标，则对任何一组（X,Y,Z），变量（x,y,z）都是空间坐标。那么在一规则区域内，独立标量函数与独立变量之间存在唯一解，即x=x（X,Y,Z）, y=y（X,Y,Z）, z=z（X,Y,Z）,如果（x,y,z）为常值（x,y,z），则方程给出000x(X,Y,Z)=X, y(X,Y,Z)=y, z(X,Y,Z)=z0 00方程x(X,Y,Z)=x定义一个坐标面。当x取不同的值，就得到与x对应的一族坐00 0标面。类似地，方程y(X,Y,Z)=y和z(X,Y,Z)=z给出另外两族坐标面。两个坐00标面的交线定义一坐标线。如y=y和z二z的交线定义一条坐标线，沿这条线只00有x在变化，该交线称为x坐标线。类似的面x=x0和z二z0的交线定义出y坐标，而x二x和y二y的交线定义出z坐标。一般来说,坐标线均为曲线。因此，变量00(x,y,z)称为曲绒坐标。通常三个坐标面在空间相交于一点，因此空间中的一点与三线的交汇点(x,y,z)有关。如果通过任何点(x,y,z)的曲线坐标线相互垂直，则称它为正交iii的，其曲线坐标系称为正交曲线坐标系。如柱坐标系(r,9,z)和球坐标系(r,9,⑺均属正交曲线坐标系。在直角坐标系(X,Y,Z)中，一点(x,y,z)的位置矢量可以写作r=X»+Yj+Zk，其中i,j,k分别为X,Y,Z坐标方向的单位矢量。因此，一个曲线坐标系(x,y,z)可以用矢量方程r=r(x,y,z)定义，如果令r,r,r分别是对(x,y,z)的偏导数(以下各xyz式中的下标x,y,z均分别表示对x,y,z的偏导数)，则它们分别是对(x,y,z)坐标线的切向矢量。对于正交坐标系有r•r=r•r=r•r=0xyyzzx注意，此处“•”为标量积符号。两相邻点之间的距离dS可以由dS2=dr•dr确定，而dr=rdx+rdy+rdz，根据方程可得xyzdS2=a2dx2+P2dy2+y2dz2其中a2=r•r，P2=r•r,y2=r•rxx yy zz因为r,r,r分别是对(x,y,z)坐标线的切向矢量，所以相对于坐标(x,y,z)的单位xyz切向矢量定义为rrre—■―,e—■/,e——z-1a2P3yTOC\o"1-5"\h\z因为r,r,r是对坐标线(x,y,z)的正交矢量，所以(e,e,e)也是对坐标线xyz 123(x,y,z)的正交矢量，于是其它任何矢量都可用它们线性地表示出。如r对x的二阶偏导数可以表示为r—ae+be+cexx1 2 3为了计算系数(a,b,c)，可以分别取r与e,e,e的标量积，结果得xx1 2 3a—e•r，b—e•r,c—e•r1xx 2xx 3xx

为了求得标量积e•r等，可以将方程和对x,y,z)求导，并利用方程。典型结果1xx为r•rr•r=aaxxx xr•r+r•r—0zxxxxzr•r=aaxxy yr•r=aaxzxz'r•r+r•r—0xxyyxx由方程、、可得a=a, b=—aa.P,c——aa.九x y, x类似地，r的其它二阶偏导数可由e,e,e的线性组合表示。完整的关系式如下;1 2 3Pze)

y3aaa

r—a(fe—甘e—fe)

xxPze)

y3r—P(上e+Pyeyya1P2yy.y、卜r—y(——e——ee+re)zza1P2y3r —ae+Pexy y1 y2r —Pe+yeyz z2 y3r—ae+yezxz1x3应变几何关系设(x,y,z)是相对于物体未变形状态时的正交曲线坐标，令(u,v,w)是物体内点(x,y,z)的位移矢量通过点(x,y,z)的坐标线切线上的投影。因与坐标线相切的单位矢量e,e,e是由方程定义，所以最初位于点(x,y,z)的一质点的位移矢量123可由位置矢量r―r(x,y,z)来定义，即,uvw△p——r+r+—raxPyyz物体变形后，最初在点(x,y,z)的质点移动到由位置矢量r*(x*,y*,z*)—r+Ap所确定的点(x*,y*,z*)。因此