第七讲_+弹性本构关系

第四讲弹性力学基础,胡才博中国科学院大学地球科学学院中国科学院计算地球动力学重点实验室,本构关系,4.1弹性应力应变关系,4.1.1一般表示4.1.2材料对称性4.1.3各向同性弹性体4.1.4弹性常数的测定4.1.5矩阵形式表达4.1.6弹性应变能,应力只取决于应变状态，与达到该状态的过程无关x=x(x，y，z，xy，yz，zx)y=y(x，y，z，xy，yz，zx).zx=zx(x，y，z，xy，yz，zx)呈线性单值连续关系的材料性质称为线弹性。,4.1.1一般表示,弹性的数学表达：,如果材料ij=f(ij)呈单值连续关系（不一定线性），则称为柯西弹性材料，即为一般意义下的弹性。,对于线性弹性材料，应力与应变是线性关系x=c11x+c12y+c13z+c14xy+c15yz+c16zxy=c21x+c22y+c23z+c24xy+c25yz+c26zxz=c31x+c32y+c33z+c34xy+c35yz+c36zxxy=c41x+c42y+c43z+c44xy+c45yz+c46zxyz=c51x+c52y+c53z+c54xy+c55yz+c56zxzx=c61x+c62y+c63z+c64xy+c65yz+c66zx系数cmn共36个，取决于材料弹性性质，与坐标系选取有关。称为广义虎克定律的一般形式。,张量形式表示ij=Cijklkl其中Cijkl称为四阶弹性张量，共81个分量。同样也取决于坐标系，服从四阶张量的坐标变换定律。,弹性张量的对称性（1）根据应力张量的对称性Cijkl=Cjikl（ij=ji）（2）根据应变张量的对称性Cijkl=Cijlk（ij=ji）独立的分量也是36个。,线弹性本构关系的一般形式：,矩阵形式：,x=c11x+c12y+c13z+c14xy+c15yz+c16zxy=c21x+c22y+c23z+c24xy+c25yz+c26zxz=c31x+c32y+c33z+c34xy+c35yz+c36zxxy=c41x+c42y+c43z+c44xy+c45yz+c46zxyz=c51x+c52y+c53z+c54xy+c55yz+c56zxzx=c61x+c62y+c63z+c64xy+c65yz+c66zx,张量形式,、分别为应力和应变列向量，C：弹性矩阵，其元素cmn为36个。,cmn和Cijkl的下标关系：,c22=C2222,c56=C2331,根据应力张量和应变张量的对称性，Cijkl也只有36个独立参数：Cijkl=Cjikl，Cijkl=Cijlk,两种形式是完全等效的。,（3）应变能存在，则弹性张量关于ij和kl也应对称Cijkl=Cklij独立的弹性常数Cijkl共有21个。弹性系数cmn也应具有对称性cmncnm独立的弹性系数cmn共有21个。,弹性对称面该面对称的两个方向具有相同的弹性关系,4.1.2材料对称性,横观各向异性材料仅具有一个弹性对称面的材料称为横观各向异性材料。设Oxy平面为材料的弹性对称面，z轴为弹性主轴。体内一点P(x,y,z)的应力和应变为，则有=CC为一般意义下的各向异性的弹性矩阵。现将z轴反向，考察其在新坐标系下的本构关系。由于弹性对称，应力应变关系应该保持不变。=C，=C应力张量的坐标变换：=xyzxy-yz-xzT,=xyzxy-yz-xzT联立以上各式，比较系数，c15=c16=c25=c25=c35=c36=c45=c46=0,以最后一个方程为例zx反号，而x，y，z和xy不变，c61c62c63c640x=c11x+c12y+c13z+c14xyy=c12x+c22y+c23z+c24xyz=c13x+c23y+c33z+c34xyxy=c14x+c24y+c34z+c44xyyz=c55yz+c56zxzx=c56yz+c66zx13个独立常数,横观各向异性材料的广义虎克定律可以表示为：,有13个独立的弹性常数。,x=c11x+c12y+c13z+c14xyy=c12x+c22y+c23z+c24xyz=c13x+c23y+c33z+c34xyxy=c14x+c24y+c34z+c44xyyz=c55yz+c56zxzx=c56yz+c66zx,正应变会产生切应力（xy）剪应变会产生正应力（x、y、z）,单斜晶体（如正长石）可简化为横观各向异性弹性体。,2.正交各向异性材料具有三个相互正交的弹性对称面的材料称为正交各向异性材料。将z轴反向，由横观各向异性材料的特点可知：c15=c16=c25=c25=c35=c36=c45=c46=0同理，将x轴反向，c14=c16=c24=c26=c34=c36=c46=c56=0将y轴反向，不产生新的结果x=c11x+c12y+c13zy=c12x+c22y+c23zz=c13x+c23y+c33zxy=c44xyyz=c55yzzx=c66zx,正交各向异性材料的广义虎克定律可表示为：,有9个独立的弹性参数，正应变仅产生正应力；剪应变仅产生剪应力。工程中，一般用三个弹性模量（E1,E2,E3），三个泊松比（v1,v2,v3），三个剪切模量（Gxy,Gyz,Gzx）表示。煤、木材、增强纤维复合材料等可简化为正交各向异性弹性体。,3.横观各向同性材料存在一个弹性对称轴（z轴），在垂直该轴的平面内材料各向同性（在此平面内所有射线方向的弹性性质均相同）。横观各向同性是正交各向异性的进一步特殊化。取两个特殊的变换：将x，y轴互换时，材料弹性关系不变c11=c22，c13=c23，c55=c66将坐标系绕z轴旋转450，剪切应力应变关系不变，得c44=0.5(c11c12),横观各向同性材料的广义虎克定律可表示为：,有5个独立参数正应变只产生正应力；剪应变只产生剪应力。工程中，常用两个杨氏模量(Exy,Ez)，两个泊松比（vxy,vz），一个剪切模量（Gz）,地层、层状岩体、复合板材等可简化为横观各向同性弹性体。,x=c11x+c12y+c13zy=c12x+c11y+c13zz=c13x+c13y+c33zxy=0.5(c11c12)xyyz=c55yzzx=c55zx,广义Hooke定律在横观各向同性材料的基础上，将x轴与z轴互换，或将y轴与z轴互换时，材料弹性关系不变，这种材料称为各向同性材料。c11=c33，c12=c13，c55=c66=0.5(c11c12)于是，独立的弹性常数减少到2个。,4.1.3各向同性弹性体,x=c11x+c12y+c12zy=c12x+c11y+c12zz=c12x+c12y+c11zxy=0.5(c11c12)xyyz=0.5(c11c12)yzzx=0.5(c11c12)zx,各向同性弹性体,令c12=，0.5(c11c12)=G、G称为Lame（拉梅）弹性常数x=2Gx+xy=Gxyy=2Gy+yz=Gyzz=2Gz+zx=Gzx=x+y+z是体积应变,x=c11x+c12y+c12zy=c12x+c11y+c12zz=c12x+c12y+c11zxy=0.5(c11c12)xyyz=0.5(c11c12)yzzx=0.5(c11c12)zx,广义Hooke定律的张量形式ij=kkij+2Gijij=CijklklCijkl=ijkl+G(ikjl+iljk)某个面上的剪切应力为零时，剪应变也为零，因此应力的主方向与应变的主方向重合。,ij=kkij+2Gij,由上式反解应变，令,此式即为虎克定律的工程形式，其中常数E，G，为熟知的杨氏模量、剪切模量和泊松比，仅两个是独立的。,应变用应力表示kk=(3+2G)kk,体积应力与体积应变关系（体积虎克定律）将等式对应相加，可得平均应力与体积应变的关系：30=(2G+3)式中0=(x+y+z)/3是平均应力。0=K式中K=(3+2G)/3是体积变形模量。,偏应力与偏应变关系x=2Gx+sx+0=2G(ex+)+将体应力与体应变关系代入：sx=2Gex同理可得：sy=2Geysz=2Gez,张量形式表示为sij=2Geij在线弹性范围内，偏应力只产生偏应变，即只产生形状改变，体积应力只产生体应变，即只产生体积改变。,sij=2Geij,0=K,广义虎克定律的偏量形式：,此形式便于塑性分析。,静水压缩实验体积模量,4.1.4弹性常数的测定,单轴拉伸实验使用物理关系，有弹性模量和泊松比：相反，有,纯剪实验使用物理方程，xy=2Gxy，因此，G也是剪切模量。,单轴应变实验有唯一应变分量约束模量：,各向同性弹性本构关系用其他参数表示：正应力只产生正应变；剪应力只产生剪应变。每个应变等于各个应力单独作用时产生的应变之和。,弹性常数的限制实验结果表明，E、G、K总为正值，有大多数材料为正值，而，有即材料弹性不可压缩，如橡胶。,4.1.5矩阵形式表达,平面应力情况,平面应变情况（重力坝）,一维情况一细长杆，长度L，横截面积S，两端受拉力P作用，伸长量为L，外力功为由于应力x=P/S，x=L/L，上式可写成,4.1.6弹性应变能,单位体积的应变能W为求应变能相对应变的偏导,三维情况考察微小六面体，应力分量ij产生的应变分量ij，各应力分量ij都只在与它相同的应变分量ij上做功，,根据能量平衡，单位体积的应变能应是所以dW=ijdij对于弹性体，应变能只取决于状态，是应变状态的单值函数W=W(ij)，应变能增量dW必须是全微分,于是对于任意的应变增量dij都应成立：这是从能量角度出发建立的弹性物体的应力-应变关系可导出如下对称性Cijkl=Cklij,将物理方程ij=Cijklkl代入dW=ijdij，考虑对称性，则W=Cijklijkl=ijij,应变能分解应变能