弹性力学基础

1一、基本物理量应力张量：在直角坐标系中，过弹性体内任一点取分别平行于三个坐标平面的三个微平面，它们的外法线方向分别为三个坐标轴的方向，将三个剪应力平行于坐标轴的两个分量；由此共得九个应力分量，记为：zzzyzxyzyyyxxzxyxxτ；每个分量的第一下标表示应力分量所在平面的外法线方向，第二下标表示应力分量的方向。应力分量的正负号规定为：当应力分量所在平面的外法线方向与某坐标轴同向时，应力分量的方向也与相应坐标轴同向；当应力分量所在平面的外法线方向与某坐标轴反向时，应力分量的方向也与相应坐标轴反向。3、应变弹性体内某一点的正应变（线应变）：设P为弹性体内任意点，过P点某一微元线段变形前的长度为l，变形后的长度为'l，定义P点l方向的正应变为：llllll'lim0。即正应变表示单位长度线段的伸长或缩短。弹性体内某一点的剪应变（角应变）：设rl和sl为过P点的两微元线段，变形前两线段相互垂直，定义变形后两线段间夹角的改变量（弧度）为角应变，夹角减小则角应变为正。应变张量：在直角坐标系中，过弹性体内任一点取分别平行三个坐标轴的线段，按上述原则定义各应变分量，得：zzzyzxyzyyyxxzxyxxε；两个下标相同的分量为正应变，其它为剪应变。关于主应变和主应变方向的讨论与主应力基本相同，可以证明，主应变方向与主应力方向重合。4、外力体积力：作用于弹性体内部每一点上，如重力、电磁力、惯性力等。设V为包含P点的微元体，作用于该微元体上的体积力为VF，则定义P点的体积力为：TzyxVVfffVFf0lim。表面力：作用于弹性体表面，如压力，约束力等。设S为包含P点的微元面，作用于该微元面上的表面力为SF，则定义P点的表面力为：TzyxSSsssSFs0lim。二、基本方程1、平衡方程2应用牛顿第二定律，建立力学平衡方程，表达应力与位移之间的关系。设P为弹性体内任意点，由P点沿三坐标轴的正向分别取长度为dx、dy和dz的三条棱边，由此构成一个微长方体。微长方体共六个面，每个面上有一个正应力和两个剪应力。按前节应力定义，过P点的X平面（平面的法线方向与X轴平行，即平面与YOZ坐标面平行。）上的应力分量为：),,(zyxxxxx、),,(zyxxyxy和),,(zyxxzxz。（图3）与该平面平行而相距dx的X平面上的应力分量为：),,(''zydxxxxxx、),,(''zydxxxyxy和),,(''zydxxxzxz。将这三个应力分量在P点作幂级数展开，并略去二次以上小量，得：dxxxxxxxx'、dxxxyxyxy'和dxxxzxzxz'。同理可得其它四个面上的应力分量。微长方体平衡必须满足三个方向上的力的平衡和三个方向上的力矩平衡。在X方向上力的平衡，按牛顿第二定律有：dxdydztudxdydzfdxdydzzdxdydxdzdyydxdzdydzdxxdydzxzxzxzxyxyxyxxxxxxx22)()()(整理得：22tufzyxxzxyxxx，其中为弹性体体积密度。同理可得其它两个方向的平衡方程：22tvfzyxyzyyyxy；22twfzyxzzzyzxz。讨论力矩平衡时，为计算方便，将坐标原点移到微长方体的重心处。由于体积力、惯性力和正应力的合力都通过微长方体的重心，因此它们对三个坐标轴的力矩都为零，即力矩平衡方程中仅包含剪应力。由于剪应力yx和zx与X轴平行，因此对X轴的力矩为零。绕X轴的力矩平衡方程为：（图4）021)(2121)(21dzdxdydzzdzdxdydydxdzdyydydxdzzyzyzyyzyzyz；整理得：022dzzdyyzyzyyzyz；略去高阶小量，得：022zyyz；即：zyyz。同理可得对于其它两个坐标轴的力矩平衡方程：yxxy；xzzx。3由此导出剪应力互等定理，九个应力分量中只有六个是独立的。2、几何方程应用变形几何关系，表达应变与位移之间的关系。（图5）X方向正应变：变形前微线段长度为dxl；变形后长度为dxxuudxxuudxl)1()('。按正应变定义，X方向正应变为：xudxdxdxxullllxx)/1('lim0。同理可得：yvyy；zwzz。XY平面剪应变：在小变形的前提下，dx和dy线段变形后的转角分别为：xvxuxvudxdxxuuvdxxvvxyxy1)()(tan；yuyvyuyxyx1tan。剪应变：)(21)(21xvyuyxxyyxxy。同理可得：)(21ywzvzyyz；)(21zuxwxzzx。共六个独立应变分量。1、本构方程材料物理性质的数学表达，又称广义虎克定理，表达应变与应力之间的关系。设材料的弹性模量为E，剪切弹性模量为G，泊松比为，三个材料常数之间的关系为：)1(2EG。可以证明，各向同性材料只有两个独立的材料常数，常用的材料常数有五个，统称拉梅系数，用其中任意两个可表达其余三个系数。如X方向受到简单拉伸，按虎克定理和横向收缩系数的泊松比关系为：xxxxE；xxzzyy。当受到纯剪时，剪应力与剪应变的关系为：xyyxxyxyGG2)(。如所有应力都存在，则按迭加原理可得应力－应变关系式：Ezzyyxxxx/)]([；Exxzzyyyy/)]([；Eyyxxzzzz/)]([；Gxyxy2/；Gyzyz2/；Gzxzx2/。或以应变表达应力的关系式：4)]()1[()21)(1(zzyyxxxxE；xyxyG2；)]()1[()21)(1(xxzzyyyyE；yzyzG2；)]()1[()21)(1(yyxxzzzzE；zxzxG2。定义体积应变为zzyyxx，则应变表达应力的关系式又可表达为：xxxxG2；xyxyG2；yyyyG2；yzyzG2；zzzzG2；zxzxG2；其中拉梅系数)21)(1(E。2、边界条件和初始条件弹性力学问题的边界条件分成两类：力边界条件和位移边界条件。力边界条件：设P为弹性体表面受到表面力区域中的一点，取一包含P点的微四面体，四面体的三个界面平行于坐标平面，另一面则为包含P点的弹性体表面曲面（图6）。由于四面体很小，表面曲面可近视为三角形斜面。设表面在P点的外法线方向为N，方向余弦分别为xn、yn和zn，表面力强度为Tzyxssss，斜面面积为dS，三个界面的面积分别为xdS、ydS和zdS。各面积间有关系：dSndSxx，dSndSyy和dSndSzz。如体力为0，则微四面体在X方向的力平衡方程为：zzxyyxxxxxdSdSdSdSs；代入面积间有关系，得：zzxyyxxxxxnnns；同理可得Y、Z方向的力平衡方程：zzyyyyxxyynnns；zzzyyzxxzxnnns。位移边界条件：uu；vv；ww，或nunu；nvnv；nvnv对于动力学问题还需加上初始条件，0t时有：),,(0zyxuu；),,(0zyxvv；),,(0zyxww；),,(0zyxutut；),,(0zyxvtvt；),,(0zyxwtwt。三、基本解法［3，p49］5三组基本方程（本构方程、几何方程、平衡方程）共15个方程，三组基本变量（位移、应力、应变）共15个变量，加上边界（及初始）条件构成弹性力学问题的基本数学形式，理论上可求解。实际工作中将问题简化，一般以三个位移分量或六个应力分量为未知数求解，称为“位移法”或“力法”。1、位移法以三个位移分量为基本未知量，将其它未知量的方程和边界条件用位移来表示。将几何方程代入本构方程，消去应变分量，得：xuGzwyvxuxx2)(；)(xvyuGxy；yvGzwyvxuyy2)(；)(ywzvGyz；zwGzwyvxuzz2)(；)(zuxwGzx；将各应力代入平衡方程，X方向的平衡方程为：222222222222222222222)()()()()()()(2)(tufuGzwyvxuxGfzuyuxuGzxwyxvxuGzwyvxuxfzuzxwGyxvyuGxuGzwyvxuxfzyxxxxxzxyxxx同理：222)()(tvfvGzwyvxuyGy；222)()(twfwGzwyvxuzGz；2、力法四、圣维南（Saint-Venant）原理[1，p152]又称局部影响原理。如作用于弹性体边界局部表面S上的一组外力，被另一组与其等效（合力和合力矩相等）的力系所代替，则在弹性体内产生的应力的改变随着与S的距离的增加而迅速衰减。应力的改变的影响区域与S的尺度为同数量级。例：两端受力的杆（外力为集中力或分布于杆端表面）；虎钳夹铁丝。五、应变能和能量原理6弹性体内单位体积的应变能称应变能密度：)222(2121zxzxyzyzxyxyzzzzyyyyxxxxijij。弹性体应变势能：dVUijijV21；弹性体外力势能（变形过程中外力所做的功的负值）：SzyxVzyxdSwsvsusdVwfvfufV)()(；弹性体总势能：VU。最小势能原理（线弹性保守系统中的虚位移原理）：在一切变形可能的状态中，真实状态相对应的总势能最小。即弹性体总势能应取极小值。7第二章弹性力学各类经典专题一、平面应力问题平面应力问题以特定受力条件的薄板为研究对象。1、几何特征弹性体为均匀厚度的薄平板，其厚度远小于板平面内的最小尺寸。2、外力特征板的上下表面上无表面菏载；作用于板边缘的菏载平行于板面，且沿厚度不变，或虽沿厚度变化但对称于板的中间平面；体积力也平行于板面且沿厚度不变。取板的中间平面为XY坐标面，即Z轴方向垂直于板面。由于板的上下表面上无外力，即0)2/()2/()2/(hzhzhzzyzxzz。严格意义上，在板的内部此三个应力分量不为零，但由于板很薄，且在此外力特征下，板不受弯曲作用，因此板的内部此三个应力分量数值都很小，可近似为零。于是板内仅包含xx、yy和xy三个应力，由于此三个应力都属于XY平面内的应力，故称这类问题为平面应力问题。由于板很薄，故此三个应力沿厚度的变化较小，可假定它们沿厚度为均匀分布，或认为所计算的是此三个应力沿厚度方向的平均值，即它们的值与Z变量无关。由上述应力假定，从本构方程可得：02/Gyzyz；02/Gzxzx；EEyyxxyyxxzzzz/)(/