拓宏亮, 吴涛, 卢智先, 马晓平

(1.长安大学 理学院, 陕西 西安 710064; 2.长安大学 建筑工程学院, 陕西 西安 710064;3.西北工业大学 航空学院, 陕西 西安 710072; 4.中国科学院 工程热物理研究所, 北京 100190)

与金属材料相比，复合材料有相对较好的抗疲劳特性，在当前的飞机结构设计中，通常以静力覆盖疲劳(即认为在满足静力设计要求情况下，疲劳性能自然满足)为准则，但在交变载荷下，疲劳问题一定程度上必然存在。疲劳循环载荷作用下，复合材料的宏观力学性能(诸如刚度和强度)将会发生退化，细观组分纤维和基体的疲劳性能迥异，疲劳损伤起始和演化规律又各不相同，互相诱发和耦合，导致了复杂的疲劳损伤机理，严重影响复合材料的结构完整性。复合材料的脆性、损伤隐蔽性和突然失效(sudden death)等对复合材料结构安全构成严重威胁。目前飞机结构复合材料设计许用值较小，使用较为保守，未能充分发挥其潜在的减重性能。随着复合材料开始应用于复杂的主承力结构，其疲劳和耐久性问题会更加突出。

研究复合材料疲劳问题最直接的方法就是通过疲劳试验获取特定试件在特定载荷下的疲劳寿命，虽然该方法数据可靠，但复合材料具有良好的可设计性，铺层和结构形式千变万化，分散性大，通过疲劳试验确定疲劳寿命将耗费大量的人力和物力，且试验周期长。

目前关于复合材料的疲劳损伤分析方法主要归纳为3类：疲劳寿命模型、宏观唯象模型和渐进疲劳损伤模型。疲劳寿命模型主要是指对S-N曲线或Goodman疲劳寿命曲线进行插值，结合疲劳失效准则，预测特定载荷下的疲劳寿命。此类模型不考虑疲劳损伤的微观机理、材料性能退化、损伤之间的耦合和损伤的累积，但需要进行大量的基础试验来确定不同结构形式、载荷水平和应力比下的S-N曲线。

Hashin和Rotem[1]于20世纪70年代建立了区分纤维与基体失效模式的复合材料疲劳失效准则。Philippidis和Vassilopoulos[2]将Tsai-Hill静力失效准则扩展为疲劳失效准则。 Ellyin和El-Kadi[3]利用幂函数关系式建立了应变能密度和疲劳寿命之间的关系。

宏观唯象模型是利用宏观力学性能在疲劳载荷下的退化现象来表征疲劳损伤，主要包含剩余强度模型和剩余刚度模型。Broutman和Sahu[4]提出线性剩余强度模型，该模型认为复合材料的剩余强度随疲劳循环次数的增加线性下降，Schaff和Davidson[5]对线性退化模型进行改进。Yao和Himmel[6-7]基于新的强度退化和损伤累积理论，建立了变幅疲劳载荷作用下的复合材料层合板疲劳寿命预测方法。Stojkovic等[8]提出等幅谱加载下的双参数强度退化模型，该模型仅需要进行较少的疲劳试验，可以明显减少试验数量。Philippidis和Vassilopoulos[9]根据疲劳试验结果得出复合材料层合板的模量E(n)与加载循环次数呈线性关系。Liu等[10]提出了基于刚度退化的双参数复合材料疲劳损伤模型，该模型能够阐述疲劳损伤演化的三阶段特性。

渐进疲劳损伤分析模型主要涵盖疲劳加载过程中应力分析、疲劳失效准则和材料性能退化模型三部分。Shokrieh和Lessard[11]提出了单向板正则化材料性能退化模型，采用三维非线性Hashin疲劳失效准则进行疲劳失效判定，可有效预测螺栓连接件的疲劳损伤演化和寿命。Noll等[12]考察了面内剪切非线性应力应变关系对碳纤维增强树脂基复合材料(CFRP)疲劳寿命的影响，采用断裂面失效准则评估层内损伤，基于连续介质力学与刚度退化方法处理纤维间裂纹和层间分层损伤。王丹勇[13]基于复合材料单向板的材料性能退化，发展了剩余强度模型和剩余刚度模型，并采用正则化疲劳寿命分析方法，建立了复合材料螺栓接头疲劳载荷下三维渐进疲劳损伤分析方法，研究了不同损伤模式之间的耦合作用，疲劳数值结果与疲劳试验结果吻合较好。Lian和Yao[14]利用单向板纵向、横向和面内剪切方向的疲劳特性，采用修正的Hashin准则预测疲劳过程中的损伤起始，提出新的刚度退化方法，对E-glass/epoxy层合板的寿命进行预测。Naderi和Maligno[15]采用渐进疲劳损伤分析方法，考虑了材料性能的概率特性，结合应力分析、疲劳失效分析和材料性能退化方法，分析了AS4/3501-6层合板的疲劳性能。周银华[16]以复合材料弹塑性本构为基础，引入正则化疲劳寿命模型和剩余刚度/强度模型，结合修正的金属Lemaitre疲劳模型，研究了复合材料/金属混合连接结构的疲劳损伤特性。张文姣[17]结合连续损伤理论和Ladevèze理论，建立了复合材料疲劳损伤模型，研究平均应力、铺层和加载方式对单向板疲劳性能的影响，并通过疲劳试验验证了模型准确性。刘建明等[18]采用了基于断裂韧性的疲劳失效准则，包含纤维、基体和分层3种失效模式，结合剩余强度理论、非线性损伤累积理论和正则化疲劳寿命模型，对螺栓连接结构进行疲劳寿命预测，并和王丹勇的试验数据进行对比，寿命和失效模式吻合较好。Zhao等[19]提出了剩余应变模型，并与剩余刚度模型、疲劳加载下材料性能逐渐退化方法相结合，建立了复合材料双搭接三钉连接件疲劳性能的渐进疲劳损伤模型，分析了试验件分别在三级应力水平下的疲劳寿命，试验和数值计算结果吻合。Zhang等[20]通过有限元和试验方法研究了T800大型复合材料结构件的疲劳性能。

基于复合材料疲劳损伤和寿命预测的复杂性及重要性，展开相关研究工作对于优化复合材料疲劳性能，进一步扩大复合材料的应用范围具有重要意义。本文以复合材料层合板为研究对象，展开以下研究：①将最大应力准则和基于物理失效机制的Puck准则扩展为疲劳失效准则；②建立能够表征复合材料疲劳损伤演化并能有效预测疲劳寿命的分析方法；③通过疲劳试验验证模型的准确性，结合试验和数值方法探究层合板疲劳加载下疲劳损伤演化过程和失效机理。

1 层合板疲劳试验

层合板材料为环氧树脂基碳纤维增强复合材料，试验件单层的名义厚度为0.125 mm，采用的疲劳试验机为MTS-810，具体铺层为[-45/90/45/0/45/0/0/-45/0/-45]s，试验件尺寸为135 mm×25 mm×2.5 mm，加载方式为等幅正弦波加载，应力比R=-1。考虑到层合板压缩强度小于拉伸强度，故选取应力水平时根据压缩载荷来确定，静力试验包含3件试验件，表1给出了试验结果。根据静压缩强度，选取55%,60%和65%三级应力水平进行疲劳试验，疲劳加载频率为6 Hz，疲劳过程中试验件表面温度均未发生明显升高现象。

表1 静强度试验结果

从表2的疲劳试验结果看出，层合板三级应力水平下的分散性均较小，数据可靠。65%应力水平下，疲劳寿命仅为33 699，分散性最小，仅为1.51%；当应力水平下降至55%，循环次数增加至20万～35万次之间，分散性稍有变大趋势，也验证了疲劳寿命“高应力区分散性小、低应力区分散性大”这一普遍规律。

表2 疲劳试验结果

图1为层合板疲劳断口图，按照试验标准层合板试验件两端有加强片，大部分试验件疲劳断裂发生在中间段考核区，个别试验件的断口位于加强片和考核区的倒角处，疲劳损伤从两侧自由边向内侧发展，随着损伤扩展，损伤覆盖整个横截面，最终出现大量的纤维断裂，破坏断面参差不齐，出现了大面积的基体损伤。

图1 层合板典型疲劳断口

2 层合板疲劳损伤模型

2.1 疲劳失效准则及损伤演化准则

复合材料在疲劳加载下，随着循环次数的增加，材料性能会发生退化，当强度退化至所施加应力时，复合材料会突然失效，因此可以认为复合材料的疲劳失效是某种意义上单次循环载荷加载下的静力失效。复合材料有多种疲劳失效模式，如纤维拉伸断裂、纤维屈曲和基体开裂等失效模式，基于这一假设，结合最大应力失效准则和基于物理失效机制的Puck准则，并将其扩展至疲劳失效准则。

对于纤维疲劳损伤，将最大应力失效准则扩展为疲劳失效准则，根据该准则判定纤维疲劳损伤的起始，失效准则如下所示：

纤维拉伸模式

(1)

纤维压缩模式

(2)

疲劳分析过程中发生纤维损伤后,纤维损伤因子按照(3)式计算。

(3)

(4)

对于基体疲劳损伤,将Puck静力失效准则扩展为疲劳失效准则:

(5)

(6)

针对Puck准则中断裂面角度的求解,采用黄金搜索算法和反二次插值法来确定断裂角度,进而求得应力危险系数。当应力危险系数最大值大于或等于1时,表明基体损伤起始,此时的潜在断裂面即为实际发生断裂的作用面。

损伤起始后,分别定义断裂面内的等效应力和等效应变

建立混合模式下基于能量的损伤判据

(9)

当基体损伤完全失效时，满足

(10)

(11)

根据双线性本构关系,基体损伤变量由(12)式计算所得。

(12)

2.2 疲劳循环加载下材料性能逐渐退化模型

本文采用基于宏观唯象的剩余强度和剩余刚度模型来表征复合材料疲劳损伤的扩展,即复合材料的强度和刚度在疲劳循环载荷下会发生逐渐退化现象。

复合材料在单轴疲劳载荷下,随着循环次数的增加,剩余强度下降,当下降至施加应力水平σ时,试验件彻底失效。复合材料结构在实际使用期间需承受不同载荷情况,因此需要研究不同应力水平下的强度退化规律。图2给出了不同应力水平下强度退化规律,在高应力水平下,退化规律呈现出“突然失效”(sudden death)现象,加载期间剩余强度几乎为常数,临近失效时,强度迅速下降,试验件彻底破坏;低应力水平下,试验件的强度退化呈现出“缓慢失效”(wear out)规律。

图2 不同应力水平下强度退化规律

Shokrieh和Lessard[11,24]提出了一个能够表征不同应力水平下强度退化特征的模型,如(13)式所示,其中,α和β为无量纲的参数,独立于应力水平,通过调整其具体数值来实现“缓慢失效”和“突然失效”强度退化特征之间的变换。

S(n,σ,R)=

(13)

式中：Ss为静强度;n为疲劳循环次数;σ为施加疲劳应力;Nf为特定应力加载下的疲劳寿命,α和β为试验拟合参数。

和强度退化一样,复合材料刚度退化也是衡量损伤的重要宏观力学性能参数。在一定应力水平下,随着循环次数的增加,剩余刚度下降,当刚度下降至临界刚度值时,试验件彻底失效。不同应力水平下刚度退化的规律不一致,和强度退化规律类似,在高应力水平下,剩余刚度几乎为一常数,临近失效时,刚度迅速下降,试验件彻底破坏,低应力水平下,刚度呈缓慢退化规律。借鉴复合材料剩余强度模型,采用不同应力比和应力水平下的正则化剩余刚度模型

E(n,σ,R)=

(14)

式中：E(n,σ,R)为剩余刚度；Es为初始无损刚度；σ为施加应力；εf为失效应变；Nf为当前施加应力σ下的疲劳寿命；λ和γ为拟合参数。

2.3 正则化疲劳寿命模型

复合材料疲劳损伤分析中,等寿命曲线的建立和插值需要大量的试验数据和计算工作,因此,有学者提出了一些解析模型并结合适量的试验数据来表征平均应力对疲劳寿命的影响。Adam等[25]对等寿命曲线的数据进行整理和分析,提出了两参数的正则化疲劳寿命模型,并通过大量层合板疲劳试验进行验证,表明模型能够合理描述任意应力比下的疲劳寿命。

a/f=(1-q)u(c+q)v

(15)

式中,a=σa/σt,q=σm/σt,c=σc/σt,σa=(σmax-σmin)/2为疲劳载荷作用下的应力幅值,σm=(σmax+σmin)/2为疲劳平均应力,f,u和v为材料拟合常数。 Gathercole等[26]通过大量复合材料疲劳寿命试验,发现复合材料的寿命曲线呈钟状分布(见图3),指数u和v分别决定钟状曲线左侧和右侧的形状,从图中可看出曲线两侧基本对称,假设u=v=A+BlgNf,且u和v都是疲劳寿命的线性函数，则可得

(16)

式中,f取1.06,A和B为试验拟合参数。

图3 正则化等疲劳寿命曲线[27]

2.4 疲劳损伤累积理论

根据损伤力学理论,复合材料在疲劳载荷下的损伤扩展规律可以通过不可逆损伤变量D来表征。

Miner理论假设复合材料在疲劳加载过程中损伤累积满足线性关系。

(17)

式中,di为第i个载荷谱造成的疲劳损伤,ni为当前加载谱下的循环次数,Ni为当前加载谱下对应的疲劳寿命。根据(17)式,如果D达到某一临界值,材料便达到疲劳寿命,临界值取1。在有限元实现过程中,定义状态变量Nratio。

neff=Nf·Nratio

(18)

式中,Nratio在每一增量步开始前进行更新。

(19)

3 层合板疲劳损伤分析和寿命预测

根据复合材料试验件的尺寸,取名义值建立数值模型,根据疲劳试验夹持情况,将夹持端固支,另一端施加疲劳载荷,不考虑加强片区损伤,把试验件两端加强片区域单元设置为线弹性单元。整个模型单元尺寸为1 mm×1 mm,单元采用实体单元(C3D8R),沿厚度方向各单层设置一个单元。根据试验选取三级应力水平,分别对其进行数值计算,材料性能参数及疲劳参数由表3～5给出。

表3 单向板材料性能参数

表4 临界能量释放率参数 J/m2

表5 疲劳性能参数

根据建立的疲劳损伤分析方法,复合材料层合板的疲劳损伤模型分析流程如图4所示。

图4 复合材料疲劳损伤模型分析流程图

为了提高疲劳数值分析模型的计算效率,可以采用如下的疲劳循环增量策略:在103疲劳循环次数内,疲劳循环增量Δn设置为50;在103～104疲劳循环次数内,疲劳循环增量Δn选取为500;在104～106疲劳循环次数内,疲劳循环增量Δn选取为5×103;当达到106疲劳循环次数时终止分析程序。

本文疲劳数值模型计算所得的疲劳寿命与疲劳试验结果的比较情况见表6和图5，可以看出本文发展的复合材料疲劳损伤分析模型预测结果和疲劳试验结果较好吻合。

表6 试验结果和计算结果对比

图5 疲劳寿命试验和数值结果对比

通过本文建立的疲劳数值模型可以揭示复合材料层合板的疲劳损伤演化过程和失效机理，根据数值计算结果，三级应力水平下损伤累积呈现的规律基本相同，本文选取65%的应力水平进行细致分析，分析各铺层疲劳损伤出现的顺序及演化规律，45°和-45°铺层呈对称分布，其规律相似，只选取了45°铺层为代表进行分析。后续损伤分析中，为了分析整体损伤，损伤云图均为相同铺层的损伤叠加云图。

图6 90°铺层基体损伤演化过程

图7 45°铺层基体疲劳损伤演化过程 图8 45°铺层纤维损伤演化过程

层合板90°铺层纤维方向和加载方向垂直，承载能力较弱，主要依靠基体承受拉伸和压缩损伤，如图6所示，在n/N=0.27时，90°铺层两自由边首先出现基体损伤，但基体损伤程度较轻(d值较小)。随后，损伤沿着两自由边不断发展，发展至一定程度后，沿中央横截面由两侧向中心区域扩展，损伤程度较为严重，90°铺层的基体损伤导致了自由边的应力集中，诱发其他铺层损伤的起始，当n/N=0.6，45°铺层出现了基体损伤(见图7)，45°铺层的基体疲劳损伤会进一步导致45°铺层出现纤维损伤(见图8)，随着疲劳载荷继续加载，n/N=0.8时，0°铺层中央位置两自由边边缘出现纤维损伤，损伤逐渐向中心扩展(见图9)，n/N=0.85时，90°铺层基体疲劳损伤贯穿整个中央面，损伤面积宽，此时45°铺层纤维和基体疲劳损伤也贯穿了层合板中央区域，相比90°铺层的损伤，45°铺层损伤面积较小， 此时0°铺层尚有承载能力，当n/N=0.95，纤维损伤几乎贯穿截面，此后纤维损伤迅速扩展并贯穿整个层合板，试验件彻底丧失刚度，数值模型中以0°铺层纤维损伤贯穿层合板为试验件彻底失效的判断标准，数值结果得到的失效模式和试验结果良好吻合。

图9 0°铺层纤维损伤演化过程

4 结 论

本文以复合材料层合板为研究对象，针对其疲劳损伤和寿命问题，展开试验和数值研究，通过系统研究得出以下结论：

1) 将复合材料的最大应力失效准则和Puck静力失效准则扩展为疲劳失效准则，结合复合材料性能退化模型，并采用复合材料正则化疲劳寿命模型和Miner线性损伤累积理论，建立了分析复合材料层合板疲劳损伤演化的数值模型。

2) 本文建立的复合材料层合板疲劳数值模型能够准确预测复合材料层合板的疲劳损伤演化过程、失效机理和疲劳寿命，并通过55%,60%和65%三级应力水平下的疲劳试验对该模型进行了试验验证。

3) 疲劳载荷作用下，复合材料层合板疲劳损伤从两侧自由边向内侧发展，90°铺层两自由边首先出现基体疲劳损伤，随后诱发45°铺层的基体和纤维疲劳损伤，0°铺层纤维损伤出现较晚，迅速向中心扩展，最终损伤覆盖整个横截面。