复合材料是一种多相固体材料,是指由两种或者两种以上具有不同的化学或物理性质的组分材料组成的一种与组分材料性质不同的新材料,且各组分之间具有明显的截面。
相对于金属材料,复合材料有自己的特点,比如比强度高、比模量高等。
在低应力的交变载荷作用下,复合材料构件从纤维或基体的薄弱环节开始形成裂纹,逐渐扩展到基体与纤维的结合面上,结合面具有阻止裂纹扩展或改变扩展方向的作用。
碳纤维-铝合金复合材料在400℃高温下,强度和弹性模量基本无变化。有的复合材料具有良好的烧蚀性能。
正是由于符合材料具有性能上的一系列特点,目前飞行器结构设计越来越多地使用复合材料。复合材料结构设计具有两个重要特性:可设计性、可成型性。
纤维增强符合材料是由两种基本原材料——基体和纤维组成的,构成复合材料的基本单元是单层板(简称单层、又名铺层)。纤维可以是连续的、非连续的、编制的、单向的、双向的或者随机分布的。
以力学角度分析复合材料,一般可分为宏观力学方法和细观力学方法。
前者以复合材料的单层、层合板或层合板结构作为研究对象,分析复合材料变现的力学性能,忽略两相材料各自性能差别及相互作用,将两相材料反映在平均的表现性能上。
后者考虑两相材料的各自性能及其相互作用,研究其如何反映在平均的表现性能上(宏观的力学性能上)。
用宏观力学分析方法分析是,假设单层为连续、均匀、正交各向异性的材料;而用细观力学方法时,除宏观假设与上述相同外,还需要细观假设各组分材料是均匀(细观单层是均匀的)、连续、各向同性的材料,并将这些分析限于线弹性与小变形的范围内。
本文将讨论单向纤维增强单层的刚度和强度,给出宏观力学分析方法的结果。
涉及的知识:材料力学、弹性力学、高等数学、线性代数。
单层的正轴刚度是指单层在正轴方向(即单层材料的弹性主方向,如图中x1,x2方向。例如碳纤维-环氧复合材料单层中的碳纤维铺层方向)所显示的刚度性能。
假设讨论的复合材料限于在线弹性和小变形的情况下,材料力学中的应变叠加原理仍能适用于复合材料。
如图所示,单层在正轴下的平面应力状态只有三个应力分量σ1、σ2、τ12。书中约定应力的符号规则为:正面正向或负面负向均为正,否则为负。(注意书中下标位置和弹性力学中相比是交换了的,即τ12在书中表示在以2为法线的平面上沿1方向的剪应力;而在弹性力学中是指在以1为法线的平面上沿2方向的剪应力。由于τ12=τ21,这个区别在大多数情况下无需区分。)
书中约定线应变伸长为正,缩短为负;剪应变与两个坐标方向一致的直角变小为正,变大为负。(与弹性力学一致)
因此,由σ1引起的应变为:
ϵ1(1)=EL1σ1ϵ2(1)=−ELνLσ1}
由σ2引起的应变为:
ϵ1(2)=−ETνTσ2ϵ2(2)=ET1σ2}
由τ12引起的应变为
γ12=GLT1τ12
利用应变叠加原理可得单层在正轴方向的应变-应力关系式为:
ϵ1=EL1σ1−ETνTσ2ϵ2=−ELνLσ1+ET1σ2γ12=GLT1τ12⎭⎪⎬⎪⎫(2.1.1)
式中,EL为纵向弹性模量;ET为横向弹性模量;νL为纵向泊松比;νT为横向泊松比;GLT为面内剪切弹性模量。
这些量称为弹性的正轴工程弹性常数,一共5个,采用功的互等原理,可以证明
νTνL=ETEL(2.1.2)
因此,单层板独立的工程弹性常数为4个。
将应变-应力关系式(2.1.1)写成矩阵的形式:
⎣⎢⎡ϵ1ϵ2γ12⎦⎥⎤=⎣⎢⎡EL1−ELνL0−ETνTET1000GLT1⎦⎥⎤⎣⎢⎡σ1σ2τ12⎦⎥⎤
令:
S11=1/ELS22=1/ETS66=1/GLTS12=−νT/ETS21=−νL/EL
这些量被称为柔性分量。用柔性分量表示的单层板正轴应变-应力关系式为:
⎣⎢⎡ϵ1ϵ2γ12⎦⎥⎤=⎣⎢⎡S11S210S12S22000S66⎦⎥⎤⎣⎢⎡σ1σ2τ12⎦⎥⎤(2.1.3)
由上式可以解得应变表示的应力,得到单层板正轴应力-应变关系式:
⎣⎢⎡σ1σ2τ12⎦⎥⎤=⎣⎢⎡Q11Q210Q12Q22000Q66⎦⎥⎤⎣⎢⎡ϵ1ϵ2γ12⎦⎥⎤(2.1.4)
其中模量分量为:
Q11=mELQ22=mETQ66=GLTQ12=mνTETQ21=mνLELm=(1−νLνT)−1
模量矩阵和柔量矩阵是互逆的。且Q12=Q21,S12=S21。可以看出无论是模量还是柔量其独立分量都为4个。
由能量守恒原理可以证明,单层的弹性模量、具有重复下表的柔量分量及模量分量均为正值:
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)被称为单层作为正交各向异性材料的工程弹性常数的限制条件。
这些工程弹性常数通常是由实验测得。
偏轴是指选用的直角坐标系Oxyz和单层的正轴坐标系Ox1x2x3的关系为:x3轴与z轴方向相同,x1、x2轴相对于x、y轴的夹角为θ(铺层角,逆时针为正)。
设
m=cosθ,n=sinθ
两个不同直角坐标系之间的应力转换是静力平衡条件,应变转换是几何问题,均与材料的物理性质无关,材料力学和弹性力学中对各向同性材料所推导的转换公式完全适用于各项异性材料。
如图所示,分别用垂直x1轴和x2轴的斜截面在单元体上截出分离体,
根据静力平衡条件可得:
σ1=σxm2+σyn2+2τxymnσ2=σxn2+σym2+2τxymnτ12=(σy−σx)mn+τxy(m2−n2)⎭⎪⎬⎪⎫(3.1.1)
写成矩阵的形式:
⎣⎢⎡σ1σ2τ12⎦⎥⎤=⎣⎢⎡m2n2−mnn2m2mn2mn−2mnm2−n2⎦⎥⎤⎣⎢⎡σxσyτxy⎦⎥⎤或⎣⎢⎡σxσyτxy⎦⎥⎤=⎣⎢⎡m2n2mnn2m2−mn−2mn2mnm2−n2⎦⎥⎤⎣⎢⎡σ1σ2τ12⎦⎥⎤(3.1.2)
如图所示,平面内一点O在外力作用下移动到O′点,它在Oxy坐标系和Ox′y′坐标系中的位移分别用u、v和u′、v′表示,
则存在下列关系:
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−2mnn2m22mnmn−mnm2−n2⎦⎥⎤⎣⎢⎡ϵxϵyγxy⎦⎥⎤或⎣⎢⎡ϵxϵyγxy⎦⎥⎤=⎣⎢⎡m2n22mnn2m2−2mn−mnmnm2−n2⎦⎥⎤⎣⎢⎡ϵ1ϵ2γ12⎦⎥⎤(3.1.4)
将(3.1.2)和(3.1.4)代入(2.1.4)得:
⎣⎢⎡σxσyτxy⎦⎥⎤=⎣⎢⎡m2n2mnn2m2−mn−2mn2mnm2−n2⎦⎥⎤⎣⎢⎡Q11Q210Q12Q22000Q66⎦⎥⎤⎣⎢⎡m2n2−2mnn2m22mnmn−mnm2−n2⎦⎥⎤⎣⎢⎡ϵxϵyγxy⎦⎥⎤
可简写成
⎣⎢⎡σxσyτxy⎦⎥⎤=⎣⎢⎡Q11Q21Q61Q12Q22Q62Q16Q26Q66⎦⎥⎤⎣⎢⎡ϵxϵyγxy⎦⎥⎤(3.2.1)
其中Qij(i,j=1,2,6)被称为偏轴模量分量:
⎣⎢⎢⎢⎢⎢⎢⎢⎡Q11Q22Q12Q66Q16Q26⎦⎥⎥⎥⎥⎥⎥⎥⎤=⎣⎢⎢⎢⎢⎢⎢⎢⎡m4n4m2n2m2n2m3nmn3n4m4m2n2m2n2−mn3−m3n2m2n22m2n2m4+n4−2m2n2mn3−m3nm3n−mn34m2n24m2n2−4m2n2(m2−n2)22(mn3−m3n)−2(mn3−m3n)⎦⎥⎥⎥⎥⎥⎥⎥⎤⎣⎢⎢⎢⎡Q11Q22Q12Q66⎦⎥⎥⎥⎤Qij=Qji
同理可得:
⎣⎢⎡ϵxϵyγxy⎦⎥⎤=⎣⎢⎡S11S21S61S12S22S62S16S26S66⎦⎥⎤⎣⎢⎡σxσyτxy⎦⎥⎤(3.2.2)
其中Sij(i,j=1,2,6)被称为偏轴柔量分量:
⎣⎢⎢⎢⎢⎢⎢⎢⎡S11S22S12S66S16S26⎦⎥⎥⎥⎥⎥⎥⎥⎤=⎣⎢⎢⎢⎢⎢⎢⎢⎡m4n4m2n24m2n22m3n2mn3n4m4m2n24m2n2−2mn3−2m3n2m2n22m2n2m4+n4−8m2n22(mn3−m3n)2(mn3−m3n)m2n2m2n2−m2n2(m2−n2)2mn3−m3n−mn3+m3n⎦⎥⎥⎥⎥⎥⎥⎥⎤⎣⎢⎢⎢⎡S11S22S12S66⎦⎥⎥⎥⎤Sij=Sji
偏轴模量矩阵和偏轴柔量矩阵互为逆矩阵。
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=S111,νx=−S11S21,ηxy,x=S11S61
同理可得:
Ey=S221,νy=−S22S12,ηxy,y=S22S62Gxy=G661,ηx,xy=S66S16,ηy,xy=S66S26
其中Gxy为剪切弹性模量。
或写作:
S11=Ex1,S12=−Eyνy,S16=Gxyηx,xyS21=−Exνx,S22=Ey1,S26=Gxyηy,xyS61=Exηxy,x,S62=Eyηxy,y,S66=Gxy1⎭⎪⎪⎬⎪⎪⎫
由于Sij=Sji,偏轴工程弹性常数由具有如下性质:
νyνx=EyExηx,xyηxy,x=GxyExηy,xyηxy,y=GxyEy⎭⎪⎪⎬⎪⎪⎫
3.3.2 偏轴工程弹性常数的转换关系
从上面的关系式可以看出工程弹性常数与柔量分量不能相互线性表示,所以偏轴工程弹性常数与正轴工程弹性常数之间不能得到矩阵的形式的转化公式。而只能用如下形式表示:
Ex1=ELm4+ETn4+m2n2(GLT1−EL2νL)Ey1=ETm4+ELn4+m2n2(GLT1−EL2νL)Gxy1=GLTm4+n4+2m2n2(EL2+ET2+EL4νL2−GLT1)νx=Ex[(m4+n4)ELνL−m2n2(EL1+ET1−GLT1)]νy=Ey[(m4+n4)ELνL−m2n2(EL1+ET1−GLT1)]ηxy,x=Ex[m3n(EL2+EL2νL−GLT1)−mn3(ET2+EL2νL−GLT1)]ηxy,y=Ex[mn3(EL2+EL2νL−GLT1)−m3n(ET2+EL2νL−GLT1)]ηx,xy=Gxy[m3n(EL2+EL2νL−GLT1)−mn3(ET2+EL2νL−GLT1)]ηy,xy=Gxy[mn3(EL2+EL2νL−GLT1)−m3n(ET2+EL2νL−GLT1)]⎭⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎬⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎫
上述的偏轴工程弹性常数可以看作关于θ的函数。弹性系数Ex取极值的条件为:
dθdEx=0
解得:
sin2θ=0cos2θ=EL+ET2ETνL−GLTELETEL−ET
第一式表明θ=0时Ex有极值,恰好等于EL和ET;第二式表明在非弹性主方向有极值的条件。由于∣cos2θ∣≤1和EL>ET,可得:
若
GLT>2(1+νL)EL
Ex在非弹性主方向有极大值。
若
GLT<2(1+νT)ET
Ex在非弹性主方向有极小值。
若
2(1+νT)ET<GLT<2(1+νL)EL
Ex在非弹性主方向无极值。
目前所知的单向纤维增强材料均属于后两种情况。
前面讨论单层刚度都是基于单层为平面应力状态下的应力-应变关系。本节讨论单层三维-应力应变关系,及它与平面应力状态下应力-应变关系的联系。
如图所示,在弹性小变形的情况下,单层在任意符合右手螺旋的坐标系Oxyz下,
可以仿照平面应力状态下利用叠加原理得到应变-应力关系,推广到具有三维应力状态的情况,得到单层的一般三维应变-应力关系(广义胡克定律)为:
⎣⎢⎢⎢⎢⎢⎢⎢⎡ϵxϵyϵzγyzγzxγxy⎦⎥⎥⎥⎥⎥⎥⎥⎤=⎣⎢⎢⎢⎢⎢⎢⎢⎡S11S21S31S41S51S61S12S22S32S42S52S62S13S23S33S43S53S63S14S24S34S44S54S64S15S25S35S45S55S65S16S26S36S46S56S66⎦⎥⎥⎥⎥⎥⎥⎥⎤⎣⎢⎢⎢⎢⎢⎢⎢⎡σxσyσzτyzτzxτxy⎦⎥⎥⎥⎥⎥⎥⎥⎤
或者简写成
ϵ=Sσ(4.1.1)
式中,Sij被称为三维柔量分量。
由(4.1.1)可得单层的一般三维应力-应变关系式
⎣⎢⎢⎢⎢⎢⎢⎢⎡σxσyσzτyzτzxτxy⎦⎥⎥⎥⎥⎥⎥⎥⎤=⎣⎢⎢⎢⎢⎢⎢⎢⎡C11C21C31C41C51C61C12C22C32C42C52C62C13C23C33C43C53C63C14C24C34C44C54C64C15C25C35C45C55C65C16C26C36C46C56C66⎦⎥⎥⎥⎥⎥⎥⎥⎤⎣⎢⎢⎢⎢⎢⎢⎢⎡ϵxϵyϵzγyzγzxγxy⎦⎥⎥⎥⎥⎥⎥⎥⎤
或者简写成
σ=Cϵ(4.1.2)
式中,Cij被称为三维模量分量。
且有
C−1=SCij=CjiSij=Sji
可见独立的弹性系数实际为21个。
如图所示,将单层至于正轴坐标系Ox1x2x3上。
由于一点处的线应变ϵ1、ϵ2和ϵ3只与该点处的正应力σ1、σ2和σ3有关,同时该点处的剪应变γ23、γ31和γ12分别仅与剪应力τ23、τ31和τ12有关。所以单层正轴三维应力-应变关系式为:
⎣⎢⎢⎢⎢⎢⎢⎢⎡ϵ1ϵ2ϵ3γ23γ31γ12⎦⎥⎥⎥⎥⎥⎥⎥⎤=⎣⎢⎢⎢⎢⎢⎢⎢⎡S11S21S31000S12S22S32000S13S23S33000000S44000000S55000000S66⎦⎥⎥⎥⎥⎥⎥⎥⎤⎣⎢⎢⎢⎢⎢⎢⎢⎡σ1σ2σ3τ23τ31τ12⎦⎥⎥⎥⎥⎥⎥⎥⎤
或者简写成
ϵ=Sσ(4.2.1)
同样可得单层的正轴三维应力-应变关系式:
⎣⎢⎢⎢⎢⎢⎢⎢⎡σ1σ2σ3τ23τ31τ12⎦⎥⎥⎥⎥⎥⎥⎥⎤=⎣⎢⎢⎢⎢⎢⎢⎢⎡C11C21C31000C12C22C32000C13C23C33000000C44000000C55000000C66⎦⎥⎥⎥⎥⎥⎥⎥⎤⎣⎢⎢⎢⎢⎢⎢⎢⎡ϵ1ϵ2ϵ3γ23γ31γ12⎦⎥⎥⎥⎥⎥⎥⎥⎤
或者简写成
σ=Cϵ(4.2.2)
且有
C−1=SCij=CjiSij=Sji
可见独立的正轴弹性系数为9个。
可以证明,偏轴三维应变-应力关系式为:
⎣⎢⎢⎢⎢⎢⎢⎢⎡ϵxϵyϵzγyzγzxγxy⎦⎥⎥⎥⎥⎥⎥⎥⎤=⎣⎢⎢⎢⎢⎢⎢⎢⎡S11S21S3100S61S12S22S3200S62S13S23S3300S63000S44S540000S45S550S16S26S3600S66⎦⎥⎥⎥⎥⎥⎥⎥⎤⎣⎢⎢⎢⎢⎢⎢⎢⎡σxσyσzτyzτzxτxy⎦⎥⎥⎥⎥⎥⎥⎥⎤
偏轴三维应力-应变关系式为:
⎣⎢⎢⎢⎢⎢⎢⎢⎡σxσyσzτyzτzxτxy⎦⎥⎥⎥⎥⎥⎥⎥⎤=⎣⎢⎢⎢⎢⎢⎢⎢⎡C11C21C3100C61C12C22C3200C62C13C23C3300C63000C44C540000C45C550C16C26C3600C66⎦⎥⎥⎥⎥⎥⎥⎥⎤⎣⎢⎢⎢⎢⎢⎢⎢⎡ϵxϵyϵzγyzγzxγxy⎦⎥⎥⎥⎥⎥⎥⎥⎤
且有
C−1=SCij=CjiSij=Sji
可见正交各向异性的单层偏轴三维弹性系数为13个。可以证明,正交各向异性的单层三维独立的弹性系数与正交三维弹性系数一样,为9个。
观察偏轴三维应变-应力关系式和偏轴平面应力状态应变-应力关系式:
⎣⎢⎢⎢⎢⎢⎢⎢⎡ϵxϵyϵzγyzγzxγxy⎦⎥⎥⎥⎥⎥⎥⎥⎤=⎣⎢⎢⎢⎢⎢⎢⎢⎡S11S21S3100S61S12S22S3200S62S13S23S3300S63000S44S540000S45S550S16S26S3600S66⎦⎥⎥⎥⎥⎥⎥⎥⎤⎣⎢⎢⎢⎢⎢⎢⎢⎡σxσyσzτyzτzxτxy⎦⎥⎥⎥⎥⎥⎥⎥⎤⎣⎢⎡ϵxϵyγxy⎦⎥⎤=⎣⎢⎡S11S21S61S12S22S62S16S26S66⎦⎥⎤⎣⎢⎡σxσyτxy⎦⎥⎤
前一个等式是后一个等式的推广。
对于平面应力状态,σz=τyz=τxz=0。于是可删去前一个等式中的应力矢量σ的第三、四和五行,以及柔量矩阵S的第三、四和五列。同时,由于不关心ϵz、γyz和γzx三个量,可以删去前一个等式中的应变矢量ϵ的第三、四和五行,以及柔量矩阵S的第三、四和五行。则第一式变为第二式。
可见平面应力状态的柔量分量和偏轴三维状态的柔量分量对应相等。对于一般三维应力状态和平面应力状态也有一样的关系。
观察偏轴三维应力-应变关系式和偏轴平面应力状态应力-应变关系式:
⎣⎢⎢⎢⎢⎢⎢⎢⎡σxσyσzτyzτzxτxy⎦⎥⎥⎥⎥⎥⎥⎥⎤=⎣⎢⎢⎢⎢⎢⎢⎢⎡C11C21C3100C61C12C22C3200C62C13C23C3300C63000C44C540000C45C550C16C26C3600C66⎦⎥⎥⎥⎥⎥⎥⎥⎤⎣⎢⎢⎢⎢⎢⎢⎢⎡ϵxϵyϵzγyzγzxγxy⎦⎥⎥⎥⎥⎥⎥⎥⎤⎣⎢⎡σxσyτxy⎦⎥⎤=⎣⎢⎡Q11Q21Q61Q12Q22Q62Q16Q26Q66⎦⎥⎤⎣⎢⎡ϵxϵyγxy⎦⎥⎤
对于平面应力状态,σz=τyz=τxz=0。于是可删去前一个等式中的应力矢量σ的第三、四和五行,以及模量矩阵C的第三、四和五行。若想要删去应变矢量ϵ第三、四和五行,需要将应变矢量ϵ中的ϵz、γyz和γzx三个量的影响折合到其它量中。
由于σz=0,则由前一个矩阵方程得:
σz=C31ϵx+C32ϵy+C33ϵz+C36γxy=0
所以:
ϵz=−C331(C31ϵx+C32ϵy+C36γxy)
代入到前一个矩阵方程中的σx部分得:
σx=(C11−C33C132)ϵx+(C12−C33C13C23)ϵy+(C16−C33C13C63)γxy
由后一个矩阵方程的σx部分得:
σx=Q11ϵx+Q12ϵy+Q16γxy
于是:
Q11=C11−C33C132,Q12=C12−C33C13C23,Q16=C16−C33C13C63
同理由前一个和后一个矩阵方程的σy和τxy部分,综合可得:
Qij=Cij−C33Ci3Cj3(i,j=1,2,6)
于是在前面的基础上,将前一个矩阵方程的应变矢量ϵ第三、四和五行和模量矩阵C的第三、四和五列删去。并利用上式将Cij换成Qij。则前一个矩阵方程变为后一个矩方程。
可见平面应力状态的模量分量和偏轴三维状态的模量分量是不同的。对于一般三维应力状态和平面应力状态也有类似的关系。
单层的三维工程常数是单层在三维情况下,由单轴应力与纯剪应力确定的刚度性能参数。以单层的正轴情况为例,在(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=S111ν21=−S11S21ν31=−S11S31⎭⎪⎬⎪⎫
同理可得:
E2=S221ν12=−S22S12ν31=−S22S32⎭⎪⎬⎪⎫E3=S331ν13=−S33S13ν23=−S33S23⎭⎪⎬⎪⎫G23=S441G31=S551G12=S661⎭⎪⎬⎪⎫
由于Sij=Sji:
νjiνij=EjEi(i,j=1,2,3)
各向异性的单层板强度指标需给出5个基本强度,即:
Xt——单层板纵向拉伸强度;
Xc——单层板纵向压缩强度;
Yt——单层板横向拉伸强度;
Yc——单层板横向压缩强度;
S——单层板面内剪切强度。
同弹性系数一样,强度指标也由单层板的拉伸、压缩和剪切试验确定。单层板5个强度指标与4个工程弹性系数一起合称为复合材料单层板的九参数。
设单层板位于正轴坐标系Ox1x2x3下。应力应变下标中的123表示其方向沿正轴坐标Ox1x2x3的坐标轴。不要与材料力学中各向同性材料的主应力主应变方向混淆。
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=ELXt,ϵXc=ELXcϵYt=ETYt,ϵYc=ETYcγS=GLTS⎭⎪⎬⎪⎫
于是得:
σ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=X21F+H=Y21F+G=Z21⎭⎪⎬⎪⎫
X、Y、Z分别为主轴坐标系三个方向的强度(拉伸或者压缩),对于纤维增强板,Y≈Z,由上述三个式子,可以得到:
G=H=2X21F=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
其平面应力状态下的失效准则为:
XtXcσ12−σ1σ2+YtYcσ22+XtXcXc−Xtσ1+YtYcYc−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=XtXc1,F22=YtYc1,F66=S21F1=Xt1−Xc1,F2=Yt1−Yc1,σ6=τ12}
另外,F12应包含σ1、σ2的双轴试验测得,但一般可采用:
若设2F12=−F11,则蔡-胡失效准则等同于Hoffman失效准则;若2F12=−X21,则蔡-胡失效准则等同于蔡-希尔失效准则。
为了定量说明单层板在不失效时的安全裕度,引进强度/应力比,简称强度比,使失效准则表达式变成强度比方程,对于给定的作用应力分量,能定量地给出它的安全裕度。
单层在作用应力作用下,极限应力的某一个分量与其对应的作用应力分量之比称为强度/应力比,简称强度比,记为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材料已被破坏
纤维增强复合材料的单层的材料属性是各向异性的,所以描述其性能参数以及应力应变都必须以矩阵(矢量)表示。
正轴刚度是通过实验测得的,通常其他坐标系的参数都可以用正轴刚度参数表示出来。
复合材料的强度理论是各向同性材料力学强度理论的推广。
复合材料力学除了研究单层的性质,也会研究层合板的属性。
本文章使用limfx的vsocde插件快速发布