碳纤维复合材料层合板低速冲击损伤应力分析

目的 为掌握碳纤维复合材料板在低速冲击载荷作用下的损伤规律,延缓失效破坏,对其冲击损伤的应力状态进行研究。方法 基于ABAQUS平台,建立碳纤维复合材料层合板低速冲击有限元模型,采用Hashin失效准则和VUMAT用户子程序,对碳纤维复合材料层合板的冲击过程进行数值模拟,同时考虑层合板层内与层间失效,以此来研究低速冲击条件下复合材料的损伤机理,分析冲击损伤过程中的应力变化趋势,讨论应力的分布状态。重点研究铺层角度及铺层距离冲头远近对应力的影响。结果 不同角度铺层的应力传播轨迹均沿着纤维方向和垂直于纤维方向同时扩展,应力均先增加至极限值而后迅速下降;铺层角度越大,板料的承载能力越弱,0°铺层的极限应力为1 432 MPa,而90°铺层的极限应力降至1 206 MPa;离冲头越远的铺层应力越小,达到峰值的时间更早且率先下降,说明远离冲头的铺层更早发生失效。结论 揭示了碳纤维层合板在低速冲击载荷作用下的应力状态及其对损伤的影响规律,能够为复合材料层合板零件设计提供参考。     

T300/QY8911层合板低速冲击试验及有限元模拟

对T300/QY8911复合材料层合板进行了低速冲击试验研究及数值仿真模拟。通过自由落体装置对层板进行冲击,并使用超声C扫描技术检测了层板冲击后的损伤状态,获得了不同能量下层板内部的损伤面积。建立了用于预测复合材料层合板在低速冲击作用下损伤演化的3D有限元模型,模型包含了用于模拟分层损伤的界面元和用于模拟纤维断裂、纤维挤压、基体开裂、基体挤裂等面内损伤形式的3D实体单元。该模型考虑了面内基体损伤对层间强度的影响。本文中的数值仿真结果和试验结果的对比验证了模型的合理性和有效性,文中还分析了影响低速冲击后层板内部分层面积的主要因素。     

基于连续介质损伤力学的复合材料层合板低速冲击损伤模型

基于连续介质损伤力学（CDM）方法,建立了分析复合材料层合板低速冲击问题的三维数值模型。该模型考虑了层内损伤（纤维和基体损伤）、层间分层损伤和剪切非线性行为,采用最大应变失效准则预测纤维损伤的萌生,双线性损伤本构模型表征纤维损伤演化,基于物理失效机制的三维Puck准则判断基体损伤的起始,根据断裂面内等效应变建立混合模式下基体损伤扩展准则。横向基体拉伸强度和面内剪切强度采用基于断裂力学假设的就地强度（in-situ strength）。纤维和基体损伤本构关系中引入单元特征长度,有效降低模型对网格密度的依赖性。层间分层损伤情况由内聚力单元（cohesive element）预测,以二次应力准则为分层损伤的起始准则,B-K准则表征分层损伤演化。分别通过数值分析方法和试验研究方法对复合材料典型铺层层合板四级能量低速冲击下的冲击损伤和冲击响应规律进行分析,数值计算和试验测量的接触力-时间曲线、分层损伤的形状和面积较好吻合,表明该模型能够准确地预测层合板低速冲击损伤和冲击响应。     

复合材料层合板低速冲击逐渐累积损伤预测方法

针对复合材料层板在冲击载荷下,各种损伤的产生和扩展是一个随载荷、时间和空间而演变的过程,发展了复合材料层合板低速冲击逐渐累积损伤预测方法.采用刚度退化技术和改进的Chang-Chang失效准则、显式有限元法来模拟复合材料层合板受到低速冲击下逐渐损伤过程.使用所发展的方法分析了[0m/90n/0m]铺层的复合材料层合板在低速冲击过程中的逐渐损伤扩展,结果表明本文的方法能较好地模拟复合材料层板在低速冲击下的损伤扩展及变形过程,计算结果与实验结果吻合较好;对不同冲击能量下层合板损伤扩展研究表明,冲击能量与分层损伤面积成线性关系.     

纤维/镁合金混杂层合板低速冲击响应及损伤模拟

为了研究新型纤维增强镁合金混杂层合板在低速冲击下的力学响应,分别对由玻璃纤维、碳纤维和二者混杂增强的AZ31B镁合金层合板在不同冲击能量下的落锤低速冲击试验进行了数值模拟。基于镁合金各向异性塑性本构和指数关系界面脱粘内聚力本构模型,同时纤维复合材料层采用三维Hashin失效准则且引入刚度折减,编写了复合材料层板损伤的VUMAT子程序,并将该子程序嵌入ABAQUS/Explicit中实现对层合板冲击过程的模拟。研究了该纤维层合板在不同冲击能量下的动态冲击响应以及脱粘与损伤演化规律,分析了冲击载荷、形变和能量吸收随时间的变化规律。模拟结果表明:在冲击能较小时,首先在冲击背面出现基体开裂,随着冲击能的增加,层合板受冲击面出现由无明显损伤到出现基体开裂和纤维断裂的现象;与单一碳纤维增强的镁合金层合板复合材料相比,单一玻璃纤维增强的镁合金层合板在冲击载荷作用时能够吸收更多的能量,碳纤维层内混杂合适的玻璃纤维铺层能够提高碳纤维增强镁合金层合板的抗冲击性能。     

碳纤维-玻璃纤维混杂增强环氧树脂复合材料低速冲击性能及其模拟

本文基于实验和数值模拟方法研究了碳纤维-玻璃纤维混杂增强环氧树脂复合材料低速冲击性能。采用商业有限元软件ABAQUS建立了层间/层内两类混杂复合材料低速冲击模型,采用基于应变形式的Hashin失效准则模拟面内损伤;零厚度Cohesive内聚力单元预测层间分层;编写VUMAT子程序定义渐进失效过程,并结合C扫和Micro-CT扫描,分析了复合材料内部微观损伤形貌及损伤分布情况。结果表明,层间混杂结构复合材料的抗冲击性能更优,其中铺层形式为I-C的混杂复合材料抗冲击性能最佳,冲击面为玻璃纤维时混杂结构复合材料对冲击响应区别不明显,CN-1层内混杂结构复合材料抗冲击性能优于CN-2层内混杂结构复合材料。低速冲击损伤主要为冲击处纤维断裂、基体破坏及界面分层,混杂结构可有效降低冲击破坏,层间混杂结构中玻璃纤维层损伤较大,层内混杂结构损伤受混杂界面影响,碳纤维束对临近的玻璃纤维束具有保护作用。      

z-pin增韧复合材料层合板低速冲击损伤过程研究

从z-pin增韧复合材料层合板低速冲击实验出发,考虑了复合材料层合板在受低速冲击时的分层及基体开裂等四种损伤形式,引入适当的损伤判据,对T300/3234碳纤维树脂基复合材料层合板的低速冲击损伤过程进行有限元模拟。结果表明,采用碳纤维钉z-pin增韧使得冲击后层间分层区域面积减小50%左右,有限元结果与实验数据吻合。     

低速冲击作用下碳纤维复合材料铺层板的损伤分析

建立了一个有效计算模型, 以分析碳纤维复合材料层合板在低速冲击作用下的层内和层间失效行为。针对铺层板的层内损伤, 在基于应变描述的Hashin 失效准则的基础上, 建立了单层板的逐渐累积损伤分析模型;针对铺层板的脱层损伤, 建立了各向同性脱层损伤模型, 通过结合传统的应力失效准则和断裂力学中的能量释放率准则定义了界面损伤演化规律, 并在潜在产生脱层的区域模拟为粘结接触, 并将脱层损伤模型作为界面的接触行为。该计算模型通过商用有限元软件ABAQUS/ Explicit 的用户子程序实现。使用该计算模型对碳纤维增强环氧树脂复合材料层合板在横向低速冲击作用下的损伤和变形行为进行预测分析。数值仿真的结果与试验结果进行了比较, 取得了满意的结果, 验证了该模型的正确性。      

基于超弹性特性的SMA复合材料层合板低速冲击损伤数值模拟方法

基于ABAQUS有限元软件结合VC++6.0程序设计,建立了含不同铺层角度、不同排列密度形状记忆合金(SMA)纤维的复合材料层合板有限元模型。将基于Brinson本构模型的SMA分段线性超弹性模型以及判断复合材料层内失效的三维HASHIN失效准则编译至ABAQUS/VUMAT子程序,使用界面单元模拟复合材料层间区域,建立了SMA复合材料层合板的低速冲击损伤及冲击后剩余强度数值模拟方法。对比了不含SMA纤维层合板、含SMA纤维层合板、含普通金属丝层合板在不同冲击能量下的损伤响应。进一步分析了SMA纤维体积分数和直径变化对冲击响应的影响。冲击后剩余压缩强度模拟结果表明:冲击能量为16J时,含体积分数25%、直径0.5mm的SMA纤维层合板的冲击后剩余压缩强度相比不含SMA纤维层合板提高5.78%、相比含普通金属丝层合板提高4.69%。随着SMA纤维体积分数提高,层合板的抗低速冲击能力增强,当体积分数一定时,较细的(0.3mm)SMA纤维比粗的(0.6mm)SMA纤维对层合板的抗低速冲击能力增强效果更好。     

复合材料层合板分层损伤数值模拟方法研究现状

复合材料层合板在航空航天等领域受到广泛应用,分层损伤作为复合材料层合板主要的损伤形式,对复合材料结构的强度和刚度有显著的影响,是限制其重大工程应用的热点问题之一。通过实验的方法对复合材料结构进行研究,往往需要耗费大量的时间和成本,成熟的有限元数值模拟技术可以较低成本实现复合材料结构的分层行为模拟,成为分层损伤研究的重要手段。本文从数值模拟角度总结了国内外在纤维增强复合材料分层损伤取得的研究成果,并对目前主要的方法虚拟裂纹闭合技术(VCCT)、内聚力模型(CZM)及扩展有限元法(XFEM)进行阐述。最后,对其发展方向进行了展望。      

复合材料加筋板低速冲击损伤的数值模拟

建立了复合材料加筋板在横向低速冲击载荷作用下的渐进损伤有限元模型。该模型考虑了复合材料加筋板受低速冲击时的纤维断裂、基体开裂及分层脱粘等五种典型的损伤形式, 在层内采用应变描述的失效判据, 结合相应的材料性能退化方案, 通过编写VUMAT用户自定义子程序以实现相应损伤类型的判断和演化。在层间以及筋条与层板间加入界面元, 模拟层间区域的情况, 结合传统的应力失效判据和断裂力学中的能量释放率准则来判断分层损伤的起始和演化规律。通过对数值模拟结果与实验数据的比较, 验证了模型的合理性和有效性。同时探讨了不同位置、不同冲击能量以及含初始损伤(脱粘)等因素对复合材料加筋板低速冲击性能的影响。     

压缩预应力对复合材料层板抗冲击损伤性能的影响

为确定压缩预应力对复合材料层板抗冲击损伤性能的影响,首先对不同压缩预应力下的碳纤维/双马树脂CCF300/5428层板进行了低速冲击和准静态压痕试验,然后通过热揭层和冲击后压缩试验分别得到了层板分层面积和剩余强度。结果表明:压缩预应力会大幅降低层板的接触刚度和弯曲刚度,从而导致相同冲击能量下层板凹坑深度和背部基体开裂长度增大;对于准静态压痕过程和相同冲击能量下的冲击过程,分层起始载荷和峰值载荷均随压缩预应力的增大而减小;在相同冲击能量下,随着压缩预应力的增大,层板内部分层总面积及冲击能量吸收比不断增大,剩余压缩强度不断降低。因此,压缩预应力会降低复合材料层板的冲击损伤阻抗,对损伤容限性能不利,在对承受压缩载荷结构的试验验证过程中应考虑压缩预应力对抗冲击损伤性能的影响。      

Predicting low velocity impact damage and Compression-After-Impact (CAI) behaviour of composite laminates

Low-velocity impact damage can drastically reduce the residual strength of a composite structure even when the damage is barely visible. The ability to computationally predict the extent of damage and compression-after-impact (CAI) strength of a composite structure can potentially lead to the exploration of a larger design space without incurring significant time and cost penalties. A high-fidelity three-dimensional composite damage model, to predict both low-velocity impact damage and CAI strength of composite laminates, has been developed and implemented as a user material subroutine in the commercial finite element package, ABAQUS/Explicit. The intralaminar damage model component accounts for physically-based tensile and compressive failure mechanisms, of the fibres and matrix, when subjected to a three-dimensional stress state. Cohesive behaviour was employed to model the interlaminar failure between plies with a bi-linear traction–separation law for capturing damage onset and subsequent damage evolution. The virtual tests, set up in ABAQUS/Explicit, were executed in three steps, one to capture the impact damage, the second to stabilize the specimen by imposing new boundary conditions required for compression testing, and the third to predict the CAI strength. The observed intralaminar damage features, delamination damage area as well as residual strength are discussed. It is shown that the predicted results for impact damage and CAI strength correlated well with experimental testing without the need of model calibration which is often required with other damage models.     

搭接长度和铺层方式对CFRP复合材料层合板胶接结构连接性能和损伤行为的影响

针对不同搭接长度和铺层方式的碳纤维增强树脂（CFRP）复合材料层合板单搭胶接结构进行了拉伸试验,观察了试件的受力过程和失效形态,获得了载荷-位移曲线;同时基于连续损伤力学模型和三维Hashin失效准则模拟了CFRP复合材料层合板的层内损伤形成和演化,并利用内聚力模型来模拟层间及胶层的失效损伤,对CFRP复合材料层合板单搭胶接结构在拉伸作用下的失效强度和损伤机制进行了预测,通过对比验证了该数值方法的有效性;通过数值试验比较不同搭接长度和铺层方式的单搭胶接结构及双搭胶接结构的连接强度和损伤行为,并提出了一种优化的CFRP复合材料层合板胶接结构。结果表明:CFRP复合材料层合板胶接结构的极限失效载荷随着搭接长度的增大逐渐增加并趋于稳定值,且结构的失效形式逐渐从胶层自身剪切失效过渡到邻近胶层的层合板层间分层失效;CFRP复合材料层合板胶接结构的连接强度和损伤行为随着铺层方式的不同而改变,通过对3种铺层方式的对比和分析,得到性能最好的铺层方式是[0₃/90₃]_2S;在搭接长度为5~20 mm时,通过对搭接长度进行优化,得到单搭胶接结构的最优搭接长度是17 mm,双搭胶接结构的最优搭接长度是19.3 mm,与搭接长度为20 mm相比,单搭胶接结构和双搭胶接结构的连接强度分别提高了13.26%和0.43%。      

Predicting low-velocity impact damage on a stiffened composite panel

An intralaminar damage model, based on a continuum damage mechanics approach, is presented to model the damage mechanisms occurring in carbon fibre composite structures incorporating fibre tensile and compressive breakage, matrix tensile and compressive fracture, and shear failure. The damage model, together with interface elements for capturing interlaminar failure, is implemented in a finite element package and used in a detailed finite element model to simulate the response of a stiffened composite panel to low-velocity impact. Contact algorithms and friction between delaminated plies were included, to better simulate the impact event. Analyses were executed on a high performance computer (HPC) cluster to reduce the actual time required for this detailed numerical analysis. Numerical results relating to the various observed interlaminar damage mechanisms, delamination initiation and propagation, as well as the model’s ability to capture post-impact permanent indentation in the panel are discussed. Very good agreement was achieved with experimentally obtained data of energy absorbed and impactor force versus time. The extent of damage predicted around the impact site also corresponded well with the damage detected by non destructive evaluation of the tested panel.     

Element-Free Galerkin modelling of composite damage

Delamination and matrix cracking are routine damage mechanisms, observed by post-mortem analysis of laminated structures containing geometrical features such as notches or bolts. Current finite element tools cannot explicitly model an intralaminar matrix microcrack, except if the location of the damage is specified a priori. In this work, a meshless technique, the Element-Free Galerkin (EFG) method, is utilized for the first time to simulate delamination (interlaminar) and intralaminar matrix microcracking in composite laminates.     

Low-velocity impact-induced damage of continuous fiber-reinforced composite laminates. Part I. An FEM numerical model

A numerical model for simulating the process of low-velocity impact damage in composite laminates using the finite element method is presented in this paper, i.e. Part I of this two part series on the study of impact. In this model, the 9-node Lagrangian element of the Mindlin plate with consideration of large deformation analysis is employed. To analyze the transient response of the laminated plates, a modified Newmark time integration algorithm previously proposed by the authors is adopted here. We also proved that the impact process between a rigid ball and laminated plates is a stiff system, therefore a kind of A(α) stable method has been advocated here to solve the motion equation of the rigid ball. Furthermore, various types of damages including delamination, matrix cracking and fiber breakage, etc. and their mutual influences are modeled and investigated in detail. To overcome the difficulty of numerical oscillation or instability in the analysis of the dynamic contact problem between delaminated layers using the traditional penalty methods, we have employed dynamic spring constraints to simulate the contact effect, which are added to the numerical model by a kind of continuous penalty function. Moreover, an effective technique to calculate the strain energy release rate based on the Mindlin plate model is proposed, which can attain high precision. Finally, some techniques of adaptive analyses have been realized for improving the computational efficiency. Based on this model, a program has been developed for numerically simulating the damage process of cross-ply fiber-reinforced carbon/epoxy composite laminates under low-velocity impact load. In Part II, this numerical model will be verified by comparing with the experimental results. Also the impact damage will be investigated in detail using this numerical approach.