复合材料力学涉猎——单层的刚度和强度

1 绪论

1.1 复合材料

复合材料是一种多相固体材料,是指由两种或者两种以上具有不同的化学或物理性质的组分材料组成的一种与组分材料性质不同的新材料,且各组分之间具有明显的截面。

相对于金属材料,复合材料有自己的特点,比如比强度高、比模量高等。
在低应力的交变载荷作用下,复合材料构件从纤维或基体的薄弱环节开始形成裂纹,逐渐扩展到基体与纤维的结合面上,结合面具有阻止裂纹扩展或改变扩展方向的作用。
碳纤维-铝合金复合材料在400℃高温下,强度和弹性模量基本无变化。有的复合材料具有良好的烧蚀性能。
正是由于符合材料具有性能上的一系列特点,目前飞行器结构设计越来越多地使用复合材料。复合材料结构设计具有两个重要特性:可设计性、可成型性。

1.2 单层及模型假设

纤维增强符合材料是由两种基本原材料——基体和纤维组成的,构成复合材料的基本单元是单层板(简称单层、又名铺层)。纤维可以是连续的、非连续的、编制的、单向的、双向的或者随机分布的。

以力学角度分析复合材料,一般可分为宏观力学方法和细观力学方法。
前者以复合材料的单层、层合板或层合板结构作为研究对象,分析复合材料变现的力学性能,忽略两相材料各自性能差别及相互作用,将两相材料反映在平均的表现性能上。
后者考虑两相材料的各自性能及其相互作用,研究其如何反映在平均的表现性能上(宏观的力学性能上)。

用宏观力学分析方法分析是,假设单层为连续、均匀、正交各向异性的材料;而用细观力学方法时,除宏观假设与上述相同外,还需要细观假设各组分材料是均匀(细观单层是均匀的)、连续、各向同性的材料,并将这些分析限于线弹性与小变形的范围内。

本文将讨论单向纤维增强单层的刚度和强度,给出宏观力学分析方法的结果。

涉及的知识:材料力学、弹性力学、高等数学、线性代数。

2 单层的正轴刚度

单层的正轴刚度是指单层在正轴方向(即单层材料的弹性主方向,如图中x1​,x2​方向。例如碳纤维-环氧复合材料单层中的碳纤维铺层方向)所显示的刚度性能。
在任意坐标系Oxyz下的单层
假设讨论的复合材料限于在线弹性和小变形的情况下,材料力学中的应变叠加原理仍能适用于复合材料。

2.1 单层的正轴应力-应变关系

如图所示,单层在正轴下的平面应力状态只有三个应力分量σ1​、σ2​、τ12​。书中约定应力的符号规则为:正面正向或负面负向均为正,否则为负。(注意书中下标位置和弹性力学中相比是交换了的,即τ12​在书中表示在以2为法线的平面上沿1方向的剪应力;而在弹性力学中是指在以1为法线的平面上沿2方向的剪应力。由于τ12​=τ21​,这个区别在大多数情况下无需区分。)
在任意坐标系Oxyz下的单层
书中约定线应变伸长为正,缩短为负;剪应变与两个坐标方向一致的直角变小为正,变大为负。(与弹性力学一致)

因此,由σ1​引起的应变为:

ϵ1(1)​=EL​1​σ1​ϵ2(1)​=−EL​νL​​σ1​​}

由σ2​引起的应变为:

ϵ1(2)​=−ET​νT​​σ2​ϵ2(2)​=ET​1​σ2​​}

由τ12​引起的应变为

γ12​=GLT​1​τ12​

利用应变叠加原理可得单层在正轴方向的应变-应力关系式为:

ϵ1​=EL​1​σ1​−ET​νT​​σ2​ϵ2​=−EL​νL​​σ1​+ET​1​σ2​γ12​=GLT​1​τ12​​⎭⎪⎬⎪⎫​(2.1.1)

式中,EL​为纵向弹性模量;ET​为横向弹性模量;νL​为纵向泊松比;νT​为横向泊松比;GLT​为面内剪切弹性模量。
这些量称为弹性的正轴工程弹性常数,一共5个,采用功的互等原理,可以证明

νT​νL​​=ET​EL​​(2.1.2)

因此,单层板独立的工程弹性常数为4个。

将应变-应力关系式(2.1.1)写成矩阵的形式:

⎣⎢⎡​ϵ1​ϵ2​γ12​​⎦⎥⎤​=⎣⎢⎡​EL​1​−EL​νL​​0​−ET​νT​​ET​1​0​00GLT​1​​⎦⎥⎤​⎣⎢⎡​σ1​σ2​τ12​​⎦⎥⎤​

令:

S11​=1/EL​S22​=1/ET​S66​=1/GLT​S12​=−νT​/ET​S21​=−νL​/EL​

这些量被称为柔性分量。用柔性分量表示的单层板正轴应变-应力关系式为:

⎣⎢⎡​ϵ1​ϵ2​γ12​​⎦⎥⎤​=⎣⎢⎡​S11​S21​0​S12​S22​0​00S66​​⎦⎥⎤​⎣⎢⎡​σ1​σ2​τ12​​⎦⎥⎤​(2.1.3)

由上式可以解得应变表示的应力,得到单层板正轴应力-应变关系式:

⎣⎢⎡​σ1​σ2​τ12​​⎦⎥⎤​=⎣⎢⎡​Q11​Q21​0​Q12​Q22​0​00Q66​​⎦⎥⎤​⎣⎢⎡​ϵ1​ϵ2​γ12​​⎦⎥⎤​(2.1.4)

其中模量分量为:

Q11​=mEL​Q22​=mET​Q66​=GLT​Q12​=mνT​ET​Q21​=mνL​EL​m=(1−νL​νT​)−1

模量矩阵和柔量矩阵是互逆的。且Q12​=Q21​,S12​=S21​。可以看出无论是模量还是柔量其独立分量都为4个。

2.2 单层正轴刚度系数的取值范围

由能量守恒原理可以证明,单层的弹性模量、具有重复下表的柔量分量及模量分量均为正值:

EL​,ET​,GLT​>0(2.1.5)

S11​,S22​,S66​>0

Q11​,Q22​,Q66​>0

同时可以证明,有:

1−νL​νT​>0

νL2​<EL​/ET​νT2​<ET​/EL​​}(2.1.6)

式(2.1.2)、(2.1.5)和(2.1.6)被称为单层作为正交各向异性材料的工程弹性常数的限制条件。

这些工程弹性常数通常是由实验测得。

3 单层的偏轴刚度

偏轴是指选用的直角坐标系Oxyz和单层的正轴坐标系Ox1​x2​x3​的关系为:x3​轴与z轴方向相同,x1​、x2​轴相对于x、y轴的夹角为θ(铺层角,逆时针为正)。
设

m=cosθ,n=sinθ

3.1 应力转换与应变转换关系

两个不同直角坐标系之间的应力转换是静力平衡条件,应变转换是几何问题,均与材料的物理性质无关,材料力学和弹性力学中对各向同性材料所推导的转换公式完全适用于各项异性材料。

如图所示,分别用垂直x1​轴和x2​轴的斜截面在单元体上截出分离体,
在正轴坐标系123下的单层
根据静力平衡条件可得:

σ1​=σx​m2+σy​n2+2τxy​mnσ2​=σx​n2+σy​m2+2τxy​mnτ12​=(σy​−σx​)mn+τxy​(m2−n2)​⎭⎪⎬⎪⎫​(3.1.1)

写成矩阵的形式:

⎣⎢⎡​σ1​σ2​τ12​​⎦⎥⎤​=⎣⎢⎡​m2n2−mn​n2m2mn​2mn−2mnm2−n2​⎦⎥⎤​⎣⎢⎡​σx​σy​τxy​​⎦⎥⎤​或⎣⎢⎡​σx​σy​τxy​​⎦⎥⎤​=⎣⎢⎡​m2n2mn​n2m2−mn​−2mn2mnm2−n2​⎦⎥⎤​⎣⎢⎡​σ1​σ2​τ12​​⎦⎥⎤​(3.1.2)

如图所示,平面内一点O在外力作用下移动到O′点,它在Oxy坐标系和Ox′y′坐标系中的位移分别用u、v和u′、v′表示,
在正轴坐标系123下的单层
则存在下列关系:

u′=mu+nv,v′=−nu+mvx=mx′−ny′,y=nx′+my′

由链式求导法则:

∂x′∂u′​=m2∂x∂u​+n2∂y∂v​+mn(∂y∂u​+∂x∂v​)∂y′∂v′​=n2∂x∂u​+m2∂y∂v​−mn(∂y∂u​+∂x∂v​)∂y′∂u′​=−mn(∂x∂u​+∂y∂v​)+m2∂y∂u​−n2∂x∂v​∂x′∂v′​=−mn(∂x∂u​+∂y∂v​)+m2∂x∂v​−n2∂y∂u​

结合弹性力学几何方程:

ϵx​=∂x∂u​,ϵy​=∂y∂v​,γxy​=∂y∂u​+∂x∂v​ϵx′​=∂x′∂u′​,ϵy′​=∂y′∂v′​,γx′y′​=∂y′∂u′​+∂x′∂v′​

得:

ϵx′​=m2ϵx​+n2ϵy​+mnγxy​ϵy′​=m2ϵx​+n2ϵy​+mnγxy​γx′y′​=−2mn(ϵx​+ϵy​)+(m2−n2)γxy​​⎭⎪⎬⎪⎫​

将下标x′、y′换成1、2得:

ϵ1​=m2ϵx​+n2ϵy​+mnγxy​ϵ2​=m2ϵx​+n2ϵy​+mnγxy​γ12​=−2mn(ϵx​+ϵy​)+(m2−n2)γxy​​⎭⎪⎬⎪⎫​(3.1.3)

写成矩阵的形式:

⎣⎢⎡​ϵ1​ϵ2​γ12​​⎦⎥⎤​=⎣⎢⎡​m2n2−2mn​n2m22mn​mn−mnm2−n2​⎦⎥⎤​⎣⎢⎡​ϵx​ϵy​γxy​​⎦⎥⎤​或⎣⎢⎡​ϵx​ϵy​γxy​​⎦⎥⎤​=⎣⎢⎡​m2n22mn​n2m2−2mn​−mnmnm2−n2​⎦⎥⎤​⎣⎢⎡​ϵ1​ϵ2​γ12​​⎦⎥⎤​(3.1.4)

3.2 单层的偏轴应力-应变关系

将(3.1.2)和(3.1.4)代入(2.1.4)得:

⎣⎢⎡​σx​σy​τxy​​⎦⎥⎤​=⎣⎢⎡​m2n2mn​n2m2−mn​−2mn2mnm2−n2​⎦⎥⎤​⎣⎢⎡​Q11​Q21​0​Q12​Q22​0​00Q66​​⎦⎥⎤​⎣⎢⎡​m2n2−2mn​n2m22mn​mn−mnm2−n2​⎦⎥⎤​⎣⎢⎡​ϵx​ϵy​γxy​​⎦⎥⎤​

可简写成

⎣⎢⎡​σx​σy​τxy​​⎦⎥⎤​=⎣⎢⎡​Q​11​Q​21​Q​61​​Q​12​Q​22​Q​62​​Q​16​Q​26​Q​66​​⎦⎥⎤​⎣⎢⎡​ϵx​ϵy​γxy​​⎦⎥⎤​(3.2.1)

其中Q​ij​(i,j=1,2,6)被称为偏轴模量分量:

⎣⎢⎢⎢⎢⎢⎢⎢⎡​Q​11​Q​22​Q​12​Q​66​Q​16​Q​26​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​=⎣⎢⎢⎢⎢⎢⎢⎢⎡​m4n4m2n2m2n2m3nmn3​n4m4m2n2m2n2−mn3−m3n​2m2n22m2n2m4+n4−2m2n2mn3−m3nm3n−mn3​4m2n24m2n2−4m2n2(m2−n2)22(mn3−m3n)−2(mn3−m3n)​⎦⎥⎥⎥⎥⎥⎥⎥⎤​⎣⎢⎢⎢⎡​Q11​Q22​Q12​Q66​​⎦⎥⎥⎥⎤​Q​ij​=Q​ji​

同理可得:

⎣⎢⎡​ϵx​ϵy​γxy​​⎦⎥⎤​=⎣⎢⎡​S11​S21​S61​​S12​S22​S62​​S16​S26​S66​​⎦⎥⎤​⎣⎢⎡​σx​σy​τxy​​⎦⎥⎤​(3.2.2)

其中Sij​(i,j=1,2,6)被称为偏轴柔量分量:

⎣⎢⎢⎢⎢⎢⎢⎢⎡​S11​S22​S12​S66​S16​S26​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​=⎣⎢⎢⎢⎢⎢⎢⎢⎡​m4n4m2n24m2n22m3n2mn3​n4m4m2n24m2n2−2mn3−2m3n​2m2n22m2n2m4+n4−8m2n22(mn3−m3n)2(mn3−m3n)​m2n2m2n2−m2n2(m2−n2)2mn3−m3n−mn3+m3n​⎦⎥⎥⎥⎥⎥⎥⎥⎤​⎣⎢⎢⎢⎡​S11​S22​S12​S66​​⎦⎥⎥⎥⎤​Sij​=Sji​

偏轴模量矩阵和偏轴柔量矩阵互为逆矩阵。

3.3 单层的偏轴工程常数

3.3.1 偏轴工程常数定义

由于偏轴工程常数是表达偏轴下单轴应力或纯剪应力的刚度性能,因此可分别设:

σx​​=0,σy​=τxy​=0;σy​​=0,σx​=τxy​=0;τxy​​=0,σx​=σy​=0

来定义偏轴工程常数。

例如在xy平面内只有x向正应力作用,即σx​​=0,σy​=τxy​=0,由(3.2.2)得:

ϵx(x)​=S11​σx​ϵy(x)​=S21​σx​γxy(x)​=S61​σx​

定义x轴向的拉压弹性模量:

Ex​=ϵx(x)​σx​​

泊松耦合系数:

νx​≡νyx​=−ϵx(x)​ϵy(x)​​

拉剪耦合系数:

ηxy,x​=ϵx(x)​γxy(x)​​

这些弹性常数的前一个下标表示应力方向,后一个下标表示应变方向。

即可得到偏轴工程常数与柔量分量之间的关系:

Ex​=S11​1​,νx​=−S11​S21​​,ηxy,x​=S11​S61​​

同理可得:

Ey​=S22​1​,νy​=−S22​S12​​,ηxy,y​=S22​S62​​Gxy​=G66​1​,ηx,xy​=S66​S16​​,ηy,xy​=S66​S26​​

其中Gxy​为剪切弹性模量。

或写作:

S11​=Ex​1​,S12​=−Ey​νy​​,S16​=Gxy​ηx,xy​​S21​=−Ex​νx​​,S22​=Ey​1​,S26​=Gxy​ηy,xy​​S61​=Ex​ηxy,x​​,S62​=Ey​ηxy,y​​,S66​=Gxy​1​​⎭⎪⎪⎬⎪⎪⎫​

由于Sij​=Sji​,偏轴工程弹性常数由具有如下性质:

νy​νx​​=Ey​Ex​​ηx,xy​ηxy,x​​=Gxy​Ex​​ηy,xy​ηxy,y​​=Gxy​Ey​​​⎭⎪⎪⎬⎪⎪⎫​

3.3.2 偏轴工程弹性常数的转换关系

从上面的关系式可以看出工程弹性常数与柔量分量不能相互线性表示,所以偏轴工程弹性常数与正轴工程弹性常数之间不能得到矩阵的形式的转化公式。而只能用如下形式表示:

Ex​1​=EL​m4​+ET​n4​+m2n2(GLT​1​−EL​2νL​​)Ey​1​=ET​m4​+EL​n4​+m2n2(GLT​1​−EL​2νL​​)Gxy​1​=GLT​m4+n4​+2m2n2(EL​2​+ET​2​+EL​4νL2​​−GLT​1​)νx​=Ex​[(m4+n4)EL​νL​​−m2n2(EL​1​+ET​1​−GLT​1​)]νy​=Ey​[(m4+n4)EL​νL​​−m2n2(EL​1​+ET​1​−GLT​1​)]ηxy,x​=Ex​[m3n(EL​2​+EL​2νL​​−GLT​1​)−mn3(ET​2​+EL​2νL​​−GLT​1​)]ηxy,y​=Ex​[mn3(EL​2​+EL​2νL​​−GLT​1​)−m3n(ET​2​+EL​2νL​​−GLT​1​)]ηx,xy​=Gxy​[m3n(EL​2​+EL​2νL​​−GLT​1​)−mn3(ET​2​+EL​2νL​​−GLT​1​)]ηy,xy​=Gxy​[mn3(EL​2​+EL​2νL​​−GLT​1​)−m3n(ET​2​+EL​2νL​​−GLT​1​)]​⎭⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎬⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎫​

上述的偏轴工程弹性常数可以看作关于θ的函数。弹性系数Ex​取极值的条件为:

dθdEx​​=0

解得:

sin2θ=0cos2θ=EL​+ET​2ET​νL​−GLT​EL​ET​​EL​−ET​​

第一式表明θ=0时Ex​有极值,恰好等于EL​和ET​;第二式表明在非弹性主方向有极值的条件。由于∣cos2θ∣≤1和EL​>ET​,可得:

  1. 若

    GLT​>2(1+νL​)EL​​

    Ex​在非弹性主方向有极大值。

  2. 若

    GLT​<2(1+νT​)ET​​

    Ex​在非弹性主方向有极小值。

  3. 若

    2(1+νT​)ET​​<GLT​<2(1+νL​)EL​​

    Ex​在非弹性主方向无极值。

目前所知的单向纤维增强材料均属于后两种情况。

4 单层的三维应力-应变关系

前面讨论单层刚度都是基于单层为平面应力状态下的应力-应变关系。本节讨论单层三维-应力应变关系,及它与平面应力状态下应力-应变关系的联系。

4.1 单层的一般三维应力-应变关系

如图所示,在弹性小变形的情况下,单层在任意符合右手螺旋的坐标系Oxyz下,
在正轴坐标系123下的单层
可以仿照平面应力状态下利用叠加原理得到应变-应力关系,推广到具有三维应力状态的情况,得到单层的一般三维应变-应力关系(广义胡克定律)为:

⎣⎢⎢⎢⎢⎢⎢⎢⎡​ϵx​ϵy​ϵz​γyz​γzx​γxy​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​=⎣⎢⎢⎢⎢⎢⎢⎢⎡​S11​S21​S31​S41​S51​S61​​S12​S22​S32​S42​S52​S62​​S13​S23​S33​S43​S53​S63​​S14​S24​S34​S44​S54​S64​​S15​S25​S35​S45​S55​S65​​S16​S26​S36​S46​S56​S66​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​⎣⎢⎢⎢⎢⎢⎢⎢⎡​σx​σy​σz​τyz​τzx​τxy​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​

或者简写成

ϵ=Sσ(4.1.1)

式中,Sij​被称为三维柔量分量。

由(4.1.1)可得单层的一般三维应力-应变关系式

⎣⎢⎢⎢⎢⎢⎢⎢⎡​σx​σy​σz​τyz​τzx​τxy​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​=⎣⎢⎢⎢⎢⎢⎢⎢⎡​C11​C21​C31​C41​C51​C61​​C12​C22​C32​C42​C52​C62​​C13​C23​C33​C43​C53​C63​​C14​C24​C34​C44​C54​C64​​C15​C25​C35​C45​C55​C65​​C16​C26​C36​C46​C56​C66​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​⎣⎢⎢⎢⎢⎢⎢⎢⎡​ϵx​ϵy​ϵz​γyz​γzx​γxy​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​

或者简写成

σ=Cϵ(4.1.2)

式中,Cij​被称为三维模量分量。

且有

C−1=SCij​=Cji​Sij​=Sji​

可见独立的弹性系数实际为21个。

4.2 单层的正轴三维应力-应变关系

如图所示,将单层至于正轴坐标系Ox1​x2​x3​上。
在正轴坐标系123下的单层
由于一点处的线应变ϵ1​、ϵ2​和ϵ3​只与该点处的正应力σ1​、σ2​和σ3​有关,同时该点处的剪应变γ23​、γ31​和γ12​分别仅与剪应力τ23​、τ31​和τ12​有关。所以单层正轴三维应力-应变关系式为:

⎣⎢⎢⎢⎢⎢⎢⎢⎡​ϵ1​ϵ2​ϵ3​γ23​γ31​γ12​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​=⎣⎢⎢⎢⎢⎢⎢⎢⎡​S11​S21​S31​000​S12​S22​S32​000​S13​S23​S33​000​000S44​00​0000S55​0​00000S66​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​⎣⎢⎢⎢⎢⎢⎢⎢⎡​σ1​σ2​σ3​τ23​τ31​τ12​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​

或者简写成

ϵ=Sσ(4.2.1)

同样可得单层的正轴三维应力-应变关系式:

⎣⎢⎢⎢⎢⎢⎢⎢⎡​σ1​σ2​σ3​τ23​τ31​τ12​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​=⎣⎢⎢⎢⎢⎢⎢⎢⎡​C11​C21​C31​000​C12​C22​C32​000​C13​C23​C33​000​000C44​00​0000C55​0​00000C66​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​⎣⎢⎢⎢⎢⎢⎢⎢⎡​ϵ1​ϵ2​ϵ3​γ23​γ31​γ12​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​

或者简写成

σ=Cϵ(4.2.2)

且有

C−1=SCij​=Cji​Sij​=Sji​

可见独立的正轴弹性系数为9个。

4.3 单层的偏轴三维应力-应变关系

可以证明,偏轴三维应变-应力关系式为:

⎣⎢⎢⎢⎢⎢⎢⎢⎡​ϵx​ϵy​ϵz​γyz​γzx​γxy​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​=⎣⎢⎢⎢⎢⎢⎢⎢⎡​S11​S21​S31​00S61​​S12​S22​S32​00S62​​S13​S23​S33​00S63​​000S44​S54​0​000S45​S55​0​S16​S26​S36​00S66​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​⎣⎢⎢⎢⎢⎢⎢⎢⎡​σx​σy​σz​τyz​τzx​τxy​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​

偏轴三维应力-应变关系式为:

⎣⎢⎢⎢⎢⎢⎢⎢⎡​σx​σy​σz​τyz​τzx​τxy​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​=⎣⎢⎢⎢⎢⎢⎢⎢⎡​C11​C21​C31​00C61​​C12​C22​C32​00C62​​C13​C23​C33​00C63​​000C44​C54​0​000C45​C55​0​C16​C26​C36​00C66​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​⎣⎢⎢⎢⎢⎢⎢⎢⎡​ϵx​ϵy​ϵz​γyz​γzx​γxy​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​

且有

C−1=SCij​=Cji​Sij​=Sji​

可见正交各向异性的单层偏轴三维弹性系数为13个。可以证明,正交各向异性的单层三维独立的弹性系数与正交三维弹性系数一样,为9个。

4.4 与平面应力状态的关系

观察偏轴三维应变-应力关系式和偏轴平面应力状态应变-应力关系式:

⎣⎢⎢⎢⎢⎢⎢⎢⎡​ϵx​ϵy​ϵz​γyz​γzx​γxy​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​=⎣⎢⎢⎢⎢⎢⎢⎢⎡​S11​S21​S31​00S61​​S12​S22​S32​00S62​​S13​S23​S33​00S63​​000S44​S54​0​000S45​S55​0​S16​S26​S36​00S66​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​⎣⎢⎢⎢⎢⎢⎢⎢⎡​σx​σy​σz​τyz​τzx​τxy​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​⎣⎢⎡​ϵx​ϵy​γxy​​⎦⎥⎤​=⎣⎢⎡​S11​S21​S61​​S12​S22​S62​​S16​S26​S66​​⎦⎥⎤​⎣⎢⎡​σx​σy​τxy​​⎦⎥⎤​

前一个等式是后一个等式的推广。
对于平面应力状态,σz​=τyz​=τxz​=0。于是可删去前一个等式中的应力矢量σ的第三、四和五行,以及柔量矩阵S的第三、四和五列。同时,由于不关心ϵz​、γyz​和γzx​三个量,可以删去前一个等式中的应变矢量ϵ的第三、四和五行,以及柔量矩阵S的第三、四和五行。则第一式变为第二式。

可见平面应力状态的柔量分量和偏轴三维状态的柔量分量对应相等。对于一般三维应力状态和平面应力状态也有一样的关系。

观察偏轴三维应力-应变关系式和偏轴平面应力状态应力-应变关系式:

⎣⎢⎢⎢⎢⎢⎢⎢⎡​σx​σy​σz​τyz​τzx​τxy​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​=⎣⎢⎢⎢⎢⎢⎢⎢⎡​C11​C21​C31​00C61​​C12​C22​C32​00C62​​C13​C23​C33​00C63​​000C44​C54​0​000C45​C55​0​C16​C26​C36​00C66​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​⎣⎢⎢⎢⎢⎢⎢⎢⎡​ϵx​ϵy​ϵz​γyz​γzx​γxy​​⎦⎥⎥⎥⎥⎥⎥⎥⎤​⎣⎢⎡​σx​σy​τxy​​⎦⎥⎤​=⎣⎢⎡​Q​11​Q​21​Q​61​​Q​12​Q​22​Q​62​​Q​16​Q​26​Q​66​​⎦⎥⎤​⎣⎢⎡​ϵx​ϵy​γxy​​⎦⎥⎤​

对于平面应力状态,σz​=τyz​=τxz​=0。于是可删去前一个等式中的应力矢量σ的第三、四和五行,以及模量矩阵C的第三、四和五行。若想要删去应变矢量ϵ第三、四和五行,需要将应变矢量ϵ中的ϵz​、γyz​和γzx​三个量的影响折合到其它量中。
由于σz​=0,则由前一个矩阵方程得:

σz​=C31​ϵx​+C32​ϵy​+C33​ϵz​+C36​γxy​=0

所以:

ϵz​=−C33​1​(C31​ϵx​+C32​ϵy​+C36​γxy​)

代入到前一个矩阵方程中的σx​部分得:

σx​=(C11​−C33​C132​​)ϵx​+(C12​−C33​C13​C23​​)ϵy​+(C16​−C33​C13​C63​​)γxy​

由后一个矩阵方程的σx​部分得:

σx​=Q​11​ϵx​+Q​12​ϵy​+Q​16​γxy​

于是:

Q​11​=C11​−C33​C132​​,Q​12​=C12​−C33​C13​C23​​,Q​16​=C16​−C33​C13​C63​​

同理由前一个和后一个矩阵方程的σy​和τxy​部分,综合可得:

Q​ij​=Cij​−C33​Ci3​Cj3​​(i,j=1,2,6)

于是在前面的基础上,将前一个矩阵方程的应变矢量ϵ第三、四和五行和模量矩阵C的第三、四和五列删去。并利用上式将Cij​换成Q​ij​。则前一个矩阵方程变为后一个矩方程。

可见平面应力状态的模量分量和偏轴三维状态的模量分量是不同的。对于一般三维应力状态和平面应力状态也有类似的关系。

4.5 单层的三维工程弹性常数

单层的三维工程常数是单层在三维情况下,由单轴应力与纯剪应力确定的刚度性能参数。以单层的正轴情况为例,在(4.2.1)中分别设6个应力分量不为零,其余分量均为零来定义。

例如,设σ1​​=0,σ2​=σ3​=τ23​=τ13​=τ12​=0,可得:

ϵ1(1)​=S11​σ1​ϵ2(1)​=S12​σ1​ϵ3(1)​=S13​σ1​​⎭⎪⎪⎬⎪⎪⎫​

定义主方向1轴向的拉压弹性模量:

E1​=ϵ1(1)​σ1​​

泊松耦合系数:

ν21​=−ϵ1(1)​ϵ2(1)​​ν31​=−ϵ1(1)​ϵ3(1)​​​⎭⎪⎪⎬⎪⎪⎫​

于是可得三维正轴工程常数与柔量分量之间的关系:

E1​=S11​1​ν21​=−S11​S21​​ν31​=−S11​S31​​​⎭⎪⎬⎪⎫​

同理可得:

E2​=S22​1​ν12​=−S22​S12​​ν31​=−S22​S32​​​⎭⎪⎬⎪⎫​E3​=S33​1​ν13​=−S33​S13​​ν23​=−S33​S23​​​⎭⎪⎬⎪⎫​G23​=S44​1​G31​=S55​1​G12​=S66​1​​⎭⎪⎬⎪⎫​

由于Sij​=Sji​:

νji​νij​​=Ej​Ei​​(i,j=1,2,3)

5 单层板的强度

5.1 单层板的强度指标

各向异性的单层板强度指标需给出5个基本强度,即:
Xt​——单层板纵向拉伸强度;
Xc​——单层板纵向压缩强度;
Yt​——单层板横向拉伸强度;
Yc​——单层板横向压缩强度;
S——单层板面内剪切强度。
同弹性系数一样,强度指标也由单层板的拉伸、压缩和剪切试验确定。单层板5个强度指标与4个工程弹性系数一起合称为复合材料单层板的九参数。

5.2 单层板的失效准则

设单层板位于正轴坐标系Ox1​x2​x3​下。应力应变下标中的123表示其方向沿正轴坐标Ox1​x2​x3​的坐标轴。不要与材料力学中各向同性材料的主应力主应变方向混淆。

5.2.1 最大应力理论

单层的最大应力准则认为,复合材料在复杂应力状态下进入破坏是由于其中某个应力分量达到了材料的基本强度值。其判据式为:

σ1​=Xt​(或∣σ1​∣=Xc​)σ2​=Yt​(或∣σ2​∣=Yc​)∣τ12​∣=S​⎭⎪⎬⎪⎫​(5.2.1.1)

这三个表达式各自独立,未考虑各应力分量与材料强度的相互影响。只要有一个左边等于或大于右边的值,则认为材料已经失效。

5.2.2 最大应变理论

单层的最大应变准则认为,复合材料在复杂应力状态下进入破坏的主要原因,是材料正轴方向的应变值达到了各基本强度值所对应的应变值。其判据式为:

ϵ1​=ϵXt​​(或∣ϵ1​∣=ϵXc​​)ϵ2​=ϵYt​​(或∣ϵ2​∣=ϵYc​​)∣γ12​∣=γS​​⎭⎪⎬⎪⎫​

由于上式中的极限应变是与单轴应力或纯剪切应力状态下基本强度对应的,而材料失效是线弹性的,故

ϵXt​​=EL​Xt​​,ϵXc​​=EL​Xc​​ϵYt​​=ET​Yt​​,ϵYc​​=ET​Yc​​γS​=GLT​S​​⎭⎪⎬⎪⎫​

于是得:

σ1​−νL​σ2​=Xt​(或∣σ1​−νL​σ2​∣=Xc​)σ2​−νT​σ2​=Yt​(或∣σ2​−νT​σ2​∣=Yc​)∣τ12​∣=S​⎭⎪⎬⎪⎫​(5.2.2.1)

最大应变准则也是由3个互不影响、各自独立的表达式组成。只要满足上式中的一个,就意味着单层板被破环。

5.2.3 蔡-希尔(Tsai-Hill)理论与霍夫曼(Hoffman)理论

蔡-希尔从各向同性材料的畸变能强度理论出发,提出了单层板的蔡-希尔失效准则。其一般形式为:

F(σ2​−σ3​)2+G(σ3​−σ1​)2+H(σ1​−σ2​)2+2Lτ232​+2Mτ312​+2Nτ122​=1

系数F、G、H、L、M、N由单层板的拉、压和纯剪切实验确定。

如,单层板仅有τ12​作用时,(τ12​)max​=S,代入上式,得:

N=2S21​

同样,当单层板分别仅有σ1​、σ2​、σ3​作用时,可分别得到:

G+H=X21​F+H=Y21​F+G=Z21​​⎭⎪⎬⎪⎫​

X、Y、Z分别为主轴坐标系三个方向的强度(拉伸或者压缩),对于纤维增强板,Y≈Z,由上述三个式子,可以得到:

G=H=2X21​F=Y21​−2X21​​}

对于平面应力状态,σz​=τyz​=τxz​=0,蔡-希尔准则变为:

(Xσ1​​)2+(Yσ2​​)2−X2σ1​σ2​​+(Sτ12​​)2=1(5.2.3.1)

蔡-希尔理论未考虑拉伸、压缩强度不一致的影响。一般情况下,使用蔡-希尔理论,在结构承受拉伸载荷时,使用拉伸强度;压缩时,使用压缩强度。

Hoffman分析了单层板拉压强度不相同的这一因素,在蔡-希尔理论中补充了几个一次项,即:

C1​(σ1​−σ2​)2+C2​(σ3​−σ1​)2+C3​(σ2​−σ3​)2+C4​σ1​+C5​σ2​+C6​σ3​+C7​τ232​+C8​τ312​+C9​τ122​=1

其平面应力状态下的失效准则为:

Xt​Xc​σ12​−σ1​σ2​​+Yt​Yc​σ22​​+Xt​Xc​Xc​−Xt​​σ1​+Yt​Yc​Yc​−Yt​​σ2​+S2τ122​​=1(5.2.3.2)

5.2.4 蔡-胡失效准则

对于单层板,蔡-胡失效准则为:

F11​σ12​+2F12​σ1​σ2​+F22​σ22​+F66​σ62​+F1​σ1​+F2​σ2​=1(5.2.4.1)

其中:

F11​=Xt​Xc​1​,F22​=Yt​Yc​1​,F66​=S21​F1​=Xt​1​−Xc​1​,F2​=Yt​1​−Yc​1​,σ6​=τ12​​}

另外,F12​应包含σ1​、σ2​的双轴试验测得,但一般可采用:

\[F_{12}=-\frac{1}{2}\sqrt{F_{11}F_{22}}\]

若设2F12​=−F11​,则蔡-胡失效准则等同于Hoffman失效准则;若2F12​=−X21​,则蔡-胡失效准则等同于蔡-希尔失效准则。

5.3 单层板的强度比方程

为了定量说明单层板在不失效时的安全裕度,引进强度/应力比,简称强度比,使失效准则表达式变成强度比方程,对于给定的作用应力分量,能定量地给出它的安全裕度。

单层在作用应力作用下,极限应力的某一个分量与其对应的作用应力分量之比称为强度/应力比,简称强度比,记为R

R=σi​σi(a)​​(5.3.1)

式中,σi​为作用的应力分量,σi(a)​为对应于σi​的极限应力分量。
这里的“对应”的含义是基于假设σi​(i=1,2,3)为比例加载的。因此,对于任何应力分量,强度比值是相同的。

R的含义:
R=∞:作用力σi​=0;
∣R∣>1:结构安全,R−1为作用力达到单层板失效时尚可增加的应力倍数;
∣R∣=1:作用应力正好达到结构极限应力值;
∣R∣<1:结构已经失效,无实际意义。在设计计算中表明结构强度不够。

将强度比(5.3.1)代入强度理论可以得到强度比方程。如蔡-希尔失效判据的强度比方程为:

AR2=1A=(Xσ1​​)2+(Yσ2​​)2−X2σ1​σ2​​+(Sτ12​​)2

则:

∣R∣=A1​​⎩⎪⎪⎨⎪⎪⎧​>1安全,其值为安全裕度=1临界破坏状态<1材料已被破坏​

6 总结

  1. 纤维增强复合材料的单层的材料属性是各向异性的,所以描述其性能参数以及应力应变都必须以矩阵(矢量)表示。

  2. 正轴刚度是通过实验测得的,通常其他坐标系的参数都可以用正轴刚度参数表示出来。

  3. 复合材料的强度理论是各向同性材料力学强度理论的推广。

  4. 复合材料力学除了研究单层的性质,也会研究层合板的属性。


本文章使用limfx的vsocde插件快速发布