第七讲_+弹性本构关系

第四讲弹性力学基础胡才博中国科学院大学地球科学学院中国科学院计算地球动力学重点实验室本构关系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一般表示()ijijf弹性的数学表达：如果材料σij=f(εij)呈单值连续关系（不一定线性），则称为柯西弹性材料，即为一般意义下的弹性。•对于线性弹性材料，应力与应变是线性关系x=c11x+c12y+c13z+c14xy+c15yz+c16zxy=c21x+c22y+c23z+c24xy+c25yz+c26zxz=c31x+c32y+c33z+c34xy+c35yz+c36zxxy=c41x+c42y+c43z+c44xy+c45yz+c46zxyz=c51x+c52y+c53z+c54xy+c55yz+c56zxzx=c61x+c62y+c63z+c64xy+c65yz+c66zx系数cmn共36个，取决于材料弹性性质，与坐标系选取有关。称为广义虎克定律的一般形式。•张量形式表示ij=Cijklkl其中Cijkl称为四阶弹性张量，共81个分量。同样也取决于坐标系，服从四阶张量的坐标变换定律。•弹性张量的对称性（1）根据应力张量的对称性Cijkl=Cjikl（σij=σji）（2）根据应变张量的对称性Cijkl=Cijlk（εij=εji）独立的分量也是36个。线弹性本构关系的一般形式：矩阵形式：[][][]Cx=c11x+c12y+c13z+c14xy+c15yz+c16zxy=c21x+c22y+c23z+c24xy+c25yz+c26zxz=c31x+c32y+c33z+c34xy+c35yz+c36zxxy=c41x+c42y+c43z+c44xy+c45yz+c46zxyz=c51x+c52y+c53z+c54xy+c55yz+c56zxzx=c61x+c62y+c63z+c64xy+c65yz+c66zx张量形式ijijklklC[σ]、[ε]分别为应力和应变列向量，[C]：弹性矩阵，其元素cmn为36个。cmn和Cijkl的下标关系：m、n123456ij、kl112233122331c22=C2222,c56=C2331根据应力张量和应变张量的对称性，Cijkl也只有36个独立参数：Cijkl=Cjikl，Cijkl=Cijlk两种形式是完全等效的。（3）应变能存在，则弹性张量关于ij和kl也应对称Cijkl=Cklij独立的弹性常数Cijkl共有21个。•弹性系数cmn也应具有对称性cmn＝cnm独立的弹性系数cmn共有21个。•弹性对称面该面对称的两个方向具有相同的弹性关系4.1.2材料对称性1.横观各向异性材料仅具有一个弹性对称面的材料称为横观各向异性材料。设Oxy平面为材料的弹性对称面，z轴为弹性主轴。体内一点P(x,y,z)的应力和应变为[σ][ε]，则有[σ]=[C][ε][C]为一般意义下的各向异性的弹性矩阵。现将z轴反向，考察其在新坐标系下的本构关系。由于弹性对称，应力应变关系应该保持不变。[σ']=[C][ε']，[σ]=[C][ε]应力张量的坐标变换：[σ']=[σxσyσzσxy-σyz-σxz]T,[ε']=[εxεyεzεxy-εyz-εxz]T联立以上各式，比较系数，c15=c16=c25=c25=c35=c36=c45=c46=0xyzx'100y'010z'00-1zzyzxyzyyxxzxyxzzyzxyzyyxxzxyx212121212121zzyzxyzyyxxzxyx---zzyzxyzyyxxzxyx---212121212121以最后一个方程为例zx反号，而x，y，z和xy不变，c61＝c62＝c63＝c64＝0x=c11x+c12y+c13z+c14xyy=c12x+c22y+c23z+c24xyz=c13x+c23y+c33z+c34xyxy=c14x+c24y+c34z+c44xyyz=c55yz+c56zxzx=c56yz+c66zx13个独立常数横观各向异性材料的广义虎克定律可以表示为：有13个独立的弹性常数。x=c11x+c12y+c13z+c14xyy=c12x+c22y+c23z+c24xyz=c13x+c23y+c33z+c34xyxy=c14x+c24y+c34z+c44xyyz=c55yz+c56zxzx=c56yz+c66zx正应变会产生切应力（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=c11x+c12y+c13zy=c12x+c22y+c23zz=c13x+c23y+c33zxy=c44xyyz=c55yzzx=c66zx正交各向异性材料的广义虎克定律可表示为：有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(c11c12)横观各向同性材料的广义虎克定律可表示为：有5个独立参数正应变只产生正应力；剪应变只产生剪应力。工程中，常用两个杨氏模量(Exy,Ez)，两个泊松比（vxy,vz），一个剪切模量（Gz）地层、层状岩体、复合板材等可简化为横观各向同性弹性体。x=c11x+c12y+c13zy=c12x+c11y+c13zz=c13x+c13y+c33zxy=0.5(c11c12)xyyz=c55yzzx=c55zx•广义Hooke定律在横观各向同性材料的基础上，将x轴与z轴互换，或将y轴与z轴互换时，材料弹性关系不变，这种材料称为各向同性材料。c11=c33，c12=c13，c55=c66=0.5(c11c12)于是，独立的弹性常数减少到2个。4.1.3各向同性弹性体x=c11x+c12y+c12zy=c12x+c11y+c12zz=c12x+c12y+c11zxy=0.5(c11c12)xyyz=0.5(c11c12)yzzx=0.5(c11c12)zx各向同性弹性体令c12=，0.5(c11c12)=G、G称为Lame（拉梅）弹性常数x=2Gx+xy=Gxyy=2Gy+yz=Gyzz=2Gz+zx=Gzx=x+y+z是体积应变x=c11x+c12y+c12zy=c12x+c11y+c12zz=c12x+c12y+c11zxy=0.5(c11c12)xyyz=0.5(c11c12)yzzx=0.5(c11c12)zx•广义Hooke定律的张量形式ij=kkij+2Gijij=CijklklCijkl=ijkl+G(ikjl+iljk)某个面上的剪切应力为零时，剪应变也为零，因此应力的主方向与应变的主方向重合。ij=kkij+2Gij由上式反解应变，令此式即为虎克定律的工程形式，其中常数E，G，μ为熟知的杨氏模量、剪切模量和泊松比，仅两个是独立的。2(1)EG•应变用应力表示kk=(3+2G)kkkkijijijGGG)23(221•体积应力与体积应变关系（体积虎克定律）将等式对应相加，可得平均应力与体积应变的关系：30=(2G+3)式中0=(x+y+z)/3是平均应力。0=K式中K=(3+2G)/3是体积变形模量。•偏应力与偏应变关系x=2Gx+sx+0=2G(ex+)+将体应力与体应变关系代入：sx=2Gex同理可得：sy=2Geysz=2Gez31张量形式表示为sij=2Geij在线弹性范围内，偏应力只产生偏应变，即只产生形状改变，体积应力只产生体应变，即只产生体积改变。sij=2Geij0=K广义虎克定律的偏量形式：此形式便于塑性分析。静水压缩实验体积模量GGKkkkk323233/kk313322114.1.4弹性常数的测定•单轴拉伸实验使用物理关系，有弹性模量和泊松比：相反，有00000000xijG)3G2(GxxE)G(2xy)(12EG))((211E纯剪实验使用物理方程，xy=2Gxy，因此，G也是剪切模量。0000000yxxyijxyxyxyxyG2•单轴应变实验•有唯一应变分量•约束模量：11GM21111•各向同性弹性本构关系用其他参数表示：•正应力只产生正应变；剪应力只产生剪应变。•每个应变等于各个应力单独作用时产生的应变之和。ijkkijijEE)21)(1(1ijkkijijEE1•弹性常数的限制•实验结果表明，E、G、K总为正值，有•大多数材料为正值，而，有•即材料弹性不可压缩，如橡胶。5.015.00/1,3KEGEEzyxxxyxyE12EExzyyyzyzE12EEyxzzzxzxE12][C4.1.5矩阵形式表达)1(20)1(200)1(20001000100011][称对EC][D2/)21(02/)21(002/)21(000100010001)21)(1(][称对ED•平面应力情况0zxyzzxyyxxyyxE2/)1(00010112)(1)(yxyxzE•平面应变情况（重力坝）0zxyzz