2.5D C/SiC复合材料连续损伤本构模型

基于连续损伤力学建立了一种包含拉伸与剪切损伤变量的2.5D C/SiC复合材料连续损伤本构模型。分别开展了拉伸和剪切试验,获得应力-应变曲线,并通过拟合试验曲线获得各损伤变量的演化参数。采用子程序技术将本构模型嵌入商用有限元软件ANSYS,应用有限元法计算了材料的应力-应变曲线。考虑了拉剪损伤耦合效应,计算了偏轴拉伸情况下的应力-应变曲线。结果表明:沿经纱拉伸、沿纬纱拉伸以及面内剪切的应力-应变曲线与试验结果吻合,最大偏差依次为4.30%、3.09%及3.73%;偏轴拉伸计算与试验应力-应变曲线也吻合较好。      

纤维金属层板低速冲击试验和数值仿真

为开展纤维金属层板（FML）低速冲击有限元数值仿真研究,改进了传统的连续损伤力学（CDM）模型,然后对FML落锤低速冲击试验进行数值仿真,并与实验结果进行对比验证。分别采用5.11 J 和10.33 J冲击能量对FML进行落锤低速冲击试验,得到冲击载荷、位移和能量时程曲线,分析FML的动态响应和失效模式。建立了考虑塑性应变、压缩刚度衰减特征和纤维拉伸断裂损伤的新CDM模型,描述S2-玻璃纤维/环氧树脂（S2-galss/epoxy）复合材料的损伤本构,并编写VUMAT子程序,通过ABAQUS/Explicit求解器对FML落锤冲击试验进行数值仿真。研究结果表明：低能量冲击条件下,FML背面主要为鼓包和裂纹等失效模式,位移峰值随冲击能量的提高而增加,冲击载荷峰值在穿透前也随冲击能量的提高而增加;采用改进的CDM模型描述FML中S2-galss/epoxy复合材料铺层后,有限元数值计算可以较好地预测FML低速冲击载荷下的动态响应;有限元数值仿真结果表明,FML中第2层复合材料铺层发生的纤维断裂损伤比第1层的更严重。     

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

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

基于应变能耗散的复合材料层合板面内缺口强度分析CDM模型

基于伴随能量释放的渐进损伤演化思想,建立了复合材料层合板面内失效分析的连续介质损伤力学（CDM）分析模型,该模型包含损伤表征、损伤起始判定和损伤演化法则3个方面。基于CDM模型,通过引入损伤状态变量表征损伤,建立了平面应力状态下的材料损伤本构模型。采用损伤参量 f_E改写Hashin准则,以判定损伤的起始。损伤演化由特征长度内的应变能释放密度控制,建立了损伤状态变量关于等效应变的渐进损伤演化法则。模型中还同时考虑了面内剪切非线性和网格敏感性,并进行了对比分析。对含缺口的[90/0/±45]_3s和[（±θ）₄]_s 2类典型复合材料层合板的面内拉伸失效进行了分析,结果表明,本文中的模型能有效预测复合材料层合板的面内拉伸强度。     

连续石墨纤维增强铝基复合材料准静态拉伸损伤演化与断裂力学行为

针对连续石墨纤维增强铝基(CF/Al)复合材料,采用三种纤维排布方式的代表体积单元(RVE)建立了其细观力学有限元模型,采用准静态拉伸试验与数值模拟结合的方法,研究了其在轴向拉伸载荷下的渐进损伤与断裂力学行为。结果表明,采用基体合金和纤维原位力学性能建立的细观力学有限元模型,对轴向拉伸弹性模量和极限强度的计算结果与实验结果吻合良好,而断裂应变计算值较实验结果偏低。轴向拉伸变形中首先出现界面和基体合金损伤现象,随应变增加界面发生失效并诱发基体合金的局部失效,最后复合材料因纤维发生失效而破坏,从而出现界面脱粘后纤维拔出与基体合金撕裂共存的微观形貌。细观力学有限元分析结果表明,在复合材料制备后纤维性能衰减而强度较低条件下,改变界面强度和刚度对复合材料轴向拉伸弹塑性力学行为的影响较小,复合材料中纤维强度水平是决定该复合材料轴向拉伸力学性能的主要因素。     

形状记忆合金的力学及界面参数和体积分数对大块金属玻璃基复合材料增韧的影响

采用考虑塑性的超弹性材料模型和基于损伤塑性的准脆性材料模型,建立了三维单胞有限元模型,模拟了形状记忆合金颗粒增韧大块金属玻璃基复合材料的单调拉伸行为。讨论了形状记忆合金的力学参数、体积分数、界面厚度和界面材料参数对金属玻璃增韧效果的影响。结果表明:提高形状记忆合金的相变应变和马氏体塑性屈服应力将显著提高形状记忆合金颗粒增韧大块金属玻璃基复合材料的拉伸失效应变;形状记忆合金弹性模量超过50.0GPa、马氏体塑性屈服应力超过1.8GPa后,复合材料的拉伸失效应变变化不大。能同时兼顾失效应变和失效应力的形状记忆合金体积分数为15%左右。复合材料界面弹性模量和界面屈服应力的增加将提高复合材料的失效应力,但对失效应变影响不大;复合材料界面厚度的增加在提高失效应变的同时,也降低了复合材料的失效应力。     

针刺C/C复合材料面内拉伸强度预测

为研究针刺碳纤维增强碳基体复合材料(针刺C/C复合材料)面内拉伸强度与渐进损伤,建立了针刺C/C复合材料代表性体积单元有限元模型。模型包含无纬布层、网胎层、针刺纤维束、界面4类子区域,并考虑了孔隙的影响。采用基于应变的破坏准则及指数型损伤演化规律研究无纬布层及针刺纤维束损伤,采用弹塑性本构研究网胎层损伤,采用内聚力牵引分离定律和二次应力破坏准则分析界面损伤。通过两步法计算了孔隙对材料性能的折减效果,并得到上述4个子区域的力学性能,通过ABAQUS UMAT预测了材料的面内拉伸应力-应变曲线及各子区域损伤起始、演化与失效过程,非线性趋势及拉伸强度数值与试验数值吻合较好,验证了该模型有效性。      

SiC/Ti-6Al-4V复合材料在横向拉伸载荷下界面失效行为的数值模拟

采用有限元模拟了SiC/Ti-6Al-4V复合材料冷却过程和横向拉伸试验过程, 横向拉伸试样采用十字形试样。分别建立了平面应力和轴对称有限元模型, 采用平面应力有限元模型计算环绕纤维圆周的界面微区应力分布, 预测界面失效机制。采用轴对称有限元模型分析复合材料界面脱粘过程以及残余应力对界面径向应力分布的影响。结果表明: 对于SiC/Ti-6Al-4V复合材料十字形试样,在横向拉伸载荷下的界面失效由径向应力导致,界面失效模式为法向失效, 剪切失效模式未发生; 十字形试样在横向拉伸载荷下界面初始脱粘位置处于界面中间; 随横向拉伸应力增加, 十字形试样的界面脱粘对称向两边扩展; 界面径向应力随残余应力降低而升高。     

复合材料层合板缺口强度的CDM三维数值模型

针对复合材料层合结构缺口强度问题,基于连续损伤力学（CDM）提出了一种三维损伤数值模型。模型区分了层内损伤（纤维失效、纤维间失效）和层间分层损伤的不同失效模式。采用三维Puck准则与Aymerich准则对上述2类损伤进行判定,材料失效后基于CDM中线性软化模型对材料损伤进行演化。模型考虑了复合材料层合板子层的就位效应和剪切非线性行为。对Carlsson的AS4/3501-6缺口拉伸强度试验进行数值模拟。结果表明:分析结果与试验结果吻合良好,证明了该模型能够准确地预测含缺口复合材料层合板面内拉伸强度。      

基于弹塑性损伤本构模型的复合材料层合板破坏荷载预测

在连续损伤力学和塑性力学框架内,建立一个同时考虑塑性效应和损伤累积导致材料属性退化的复合材料弹塑性损伤本构模型。基于最近点投影回映算法,开发本构模型的应变驱动隐式积分算法以更新应力及与解答相关的状态变量,并推导与所开发算法相应的数值一致性切线刚度矩阵,保证有限元分析采用NewtonRaphson迭代法解答非线性问题的计算效率。采用断裂带模型对已开发的本构模型软化段进行规则化,以减轻有限元分析结果的网格相关性问题。对损伤变量进行粘滞规则化,并推导出相应的粘滞规则化数值一致性切线刚度张量,解决了在有限元隐式计算程序中采用含应变软化段本构关系的数值分析由于计算困难而提前终止的问题。开发包含数值积分算法的用户材料子程序UMAT,并嵌于有限元程序Abaqus v6.14中。通过对力学行为展现显著塑性效应的AS4/3501-6V型开口复合材料层合板的渐进失效分析,验证本文提出的材料本构模型的有效性。结果显示,预测结果与已报道的试验结果吻合良好,并且预测精度高于其他已有弹性损伤模型。表明已建立的弹塑性损伤本构模型能够准确预测力学行为,展现显著塑性效应的复合材料层合板的破坏荷载,为其构件和结构设计提供一种有效的分析方法。     

Fatigue life prediction of composite laminates by FEA simulation method

This paper is to simulate the fatigue damage evolution in composite laminates and predict fatigue life of the laminates with different lay-up sequences on the basis of the fatigue characteristics of longitudinal, transverse and in-plane shear directions by finite element analysis (FEA) method. In FEA model, considering the scatter of the material’s properties, each element was assigned with different material’s properties. The stress analysis was carried out in MSC Patran/Nastran, and a modified Hashin’s failure criterion was applied to predict the failure of the elements. A new stiffness degradation model was proposed and applied in the simulation and then a strength degradation model was deduced, which is coupled with the presented stiffness degradation model. The reduced or discounted elastic constants were determined based on the failure mechanism of the laminates and the restrictive conditions of orthotropic property. The fatigue behavior and fatigue life of six kinds of E-glass/epoxy composite laminates with different lay-up sequences were experimentally studied, and the S–N curves and stiffness degradation models in longitudinal, transverse and in-plane shear direction were obtained. These fatigue data were adopted in the simulation to simulate fatigue behavior and estimate life of the laminates. The simulation results, including the fatigue life predicted and the residual stiffness, were coincident with the experimental results well except for the quasi-isotropic laminate for the lack of consideration of the out-of-plane fatigue character in the simulation.     

Effects of Shear‐Transverse Coupling and Plasticity in the Formulation of an Elementary Ply Composites Damage Model,Part I: Model Formulation and Validation

Abstract: Continuum damage mechanics (CDM) models have been employed successfully in the literature to predict the response of laminated composite materials. Some sophisticated models can include the effects of non‐linear shear and transverse damage progression, plasticity and shear‐transverse damage coupling. However, these models require non‐standard test data for calibration that may not always be available to a modeller. In this two‐part study, we examine the effect of neglecting plasticity parameters, and also the effect of neglecting both plasticity and shear‐transverse coupling parameters in simplified CDM models for predicting monotonic tensile strength. In part I, we develop simplified versions of the CDM model and test their ability to accurately predict the failure response of angle‐ply laminates. In part II, we provide details of the experimental test series carried out to determine the input parameters for the models. It was found that neglecting plasticity requires some approximations in the damage development laws, but the resulting model can predict well the response of the angle‐ply laminates tested under monotonic loading to failure. Neglecting shear‐transverse coupling is acceptable for the some materials.     

Intra/inter-ply shear behaviors of continuous fiber reinforced thermoplastic composites in thermoforming processes

The thermoforming of continuous fiber reinforced thermoplastic (CFRTP) composite panels generally involves significant in-plane shear deformation. In the present work, the in-plane shear behavior of woven thermoplastic composites (Carbon/Polyphenylene Sulfide) over a range of processing temperatures is studied by bias-test experiments at different velocities. The experimental data of force versus displacement and force versus shear strain are presented for different extension velocities and temperatures. A thermo-visco-elastic model for numerical simulations of woven thermoplastic composite forming is proposed considering the influences of temperature and of strain rate. We applied a large displacement three-dimensional cohesive element with eight nodes which has been used for crack analysis in fracture mechanics by other authors, to investigate the inter-ply shear mechanism of woven thermoplastic composites. Applying three-dimensional cohesive elements, multi-plies forming simulations are performed to show inter-ply slip behaviors at different temperatures. The proposed models can be useful to predict from the properties of reinforcement and resin the intra/inter-ply shear behaviors of woven thermoplastic composites at high temperatures if experimental characterization of composite laminate behaviors is difficult to conduct.     

A family of simple and robust finite elements for linear and geometrically nonlinear analysis of laminated composite plates

A family of simple, displacement-based and shear-flexible triangular and quadrilateral flat plate/shell elements for linear and geometrically nonlinear analysis of thin to moderately thick laminate composite plates are introduced and summarized in this paper.

The developed elements are based on the first-order shear deformation theory (FSDT) and von-Karman’s large deflection theory, and total Lagrangian approach is employed to formulate the element for geometrically nonlinear analysis. The deflection and rotation functions of the element boundary are obtained from Timoshenko’s laminated composite beam functions, thus convergence can be ensured theoretically for very thin laminates and shear-locking problem is avoided naturally.

The flat triangular plate/shell element is of 3-node, 18-degree-of-freedom, and the plane displacement interpolation functions of the Allman’s triangular membrane element with drilling degrees of freedom are taken as the in-plane displacements of the element. The flat quadrilateral plate/shell element is of 4-node, 24-degree-of-freedom, and the linear displacement interpolation functions of a quadrilateral plane element with drilling degrees of freedom are taken as the in-plane displacements.

The developed elements are simple in formulation, free from shear-locking, and include conventional engineering degrees of freedom. Numerical examples demonstrate that the elements are convergent, not sensitive to mesh distortion, accurate and efficient for linear and geometric nonlinear analysis of thin to moderately thick laminates.     

Finite element model for compression after impact behaviour of stitched composites

A finite element (FE) model using coupling continuum shell elements and cohesive elements is proposed to simulate the compression after impact (CAI) behaviour and predict the CAI strength of stitched composites. Continuum shell elements with Hashin failure criterion exhibit the composite laminate damage behaviour; whilst cohesive elements using traction-separation law characterise the laminate interfaces. Impact-induced delamination is explicitly modelled by reducing material properties of damaged cohesive elements. Computational results have demonstrated the trend of increasing CAI strength with decreasing impact-induced delamination area. Spring elements are introduced into the model to represent through-thickness stitch thread in the composite laminates. Results in this study validate experimental finding that CAI strength is improved when stitching is incorporated into the composite structure. The proposed FE model reveals good CAI strength predictions and indicates good agreement with experimental results, making it a valuable tool for CAI strength prediction of stitched composites.     

Characterization of Nonlinear Behavior in Woven Composite Laminates

Nonlinear stress–strain behavior in woven glass/epoxy laminates under off-axis tension has been investigated experimentally. The validity of an orthotropic plasticity model of such behavior, with three parameters, is discussed. The parameters have been determined from the experimental results. An attempt is also made to describe the nonlinear behavior of the woven composite as a cross-ply laminate using assumed unidirectional composite properties. The nonlinear behavior of the unidirectional laminate is assumed to be described by the one-parameter plasticity model. It is shown that there is a possibility that the one-parameter plasticity model can be used to predict the nonlinear behavior of woven composites.     

Damage behaviour of laminated composites during fatigue loading

This study deals with the modelling of damage evolution in the carbon/epoxy laminated composites under static and fatigue loading. A cumulative damage model is developed on the basis of damage evolution due to static and fatigue during cyclic loading. A continuum damage mechanics (CDM)‐based damage model coupling with the micromechanics has been utilized to predict the fatigue behaviour of laminate composites. A multicriterion approach has been introduced to predict the damage behaviour in the longitudinal, transverse, and shear direction at the ply scale. Extensive experimental results on T300/EPL1012 carbon/epoxy laminates are prepared to characterize under static and fatigue loading and to evaluate the proposed model in different conditions. The obtained results show that at the beginning of the cyclic loading, the damage grows suddenly and increases until final failure, which justifies the proposed method is able to predict the evolution of the damage due to static and fatigue loading separately during cyclic loading. The obtained results show that considering damage due to static loading leads to more accurate results, particularly in low‐cycle fatigue.     

Effects of Process-Induced Voids on the Properties of Fibre Reinforced Composites

It is well known that voids have detrimental effects on the performance of composites. This study aims to provide a practical method for predicting the effects of process induced voids on the properties of composites. Representative volume elements(RVE) for carbon ?bre/epoxy composites of various ?bre volume fractions and void contents are created, and the moduli and strengths are derived by ?nite element analysis(FEA). Regression models are ?tted to the FEA data for predicting composite properties including tensile, compressive and shear. The strengths of composite laminates including tensile strength and interlaminar shear strength(ILSS) are calculated with the aid of the developed models. The model predictions are compared with various experimental data and good agreement is found. The outcome from this study provides a useful optimisation and robust design tool for realising affordable composite products when process induced voids are taken into account.