ON RESIDUAL COMPRESSIVE STRENGTH PREDICTION OF COMPOSITE SANDWICH PANELS AFTER LOW-VELOCITY IMPACT DAMAGE

This paper introduces a nonlinear finite element analysis on damage propagation behavior of composite sandwich panels under in-plane uniaxial quasi-static compression after a low velocity impact. The major damage modes due to the impact, including the residual indentation on the impacted facesheet, the initially crushed core under the impacted area, and the delamination are incorporated into the model. A consequential core crushing mechanism is incorporated into the analysis by using an element deactivation technique. Damage propagation behavior, which corresponds to those observed in sandwich compression after impact (SCAI) tests, has been successfully captured in the numerical simulation. The critical far field stress corresponding to the onset of damage propagation at specified critical locations near the damage zone are captured successfully. They show a good correlation with experimental data. These values can be used to effectively predict the residual compressive strength of low-velocity impact damaged composite sandwich panels.     

A DAMAGE MECHANICS MODEL FOR TWISTED CARBON NANOTUBE FIBERS

Carbon nanotube fibers can be fabricated by the chemical vapor deposition spinning process. They are promising for a wide range of applications such as the building blocks of high-performance composite materials and micro-electrochemical sensors. Mechanical twisting is an effective means of enhancing the mechanical properties of carbon nanotube fibers during fabrication or by post processing. However, the effects of twisting on the mechanical properties remain an unsolved issue. In this paper, we present a two-scale damage mechanics model to quantitatively investigate the effects of twisting on the mechanical properties of carbon nanotube fibers. The numerical results demonstrate that the developed damage mechanics model can effectively describe the elastic and the plastic-like behaviors of carbon nanotube fibers during the tension process. A definite range of twisting which can effectively enhance the mechanical properties of carbon nanotube fiber is given. The results can be used to guide the mechanical twisting of carbon nanotube fibers to improve their properties and help optimize the mechanical performance of carbon nanotube-based materials.     

PREDICTION OF MODE Ⅰ CRACK PROPAGATION DIRECTION IN CARBON-FIBER REINFORCED COMPOSITE PLATE

A newly developed Z fracture criterion for the composite materials was introduced, the new concepts of in-plane average strain, in-plane dilatational strain energy density factor and reciprocal characteristic function were presented. Many experimental results show that the Z fracture criterion can be well used to predict the crack propagating direction for mode f crack in carbon-fiber reinforced composite laminates.     

THREE-DIMENSIONAL EXACT MODELING OF GEOMETRIC AND MECHANICAL PROPERTIES OF WOVEN COMPOSITES

A formulation for the prediction of the influence of various parameters on the elastic moduli of three-dimensional （3D） orthogonally woven composites has been given. These parameters can be classified into different groups according to their properties, such as input design and material parameters, structural parameters etc. Some, by their nature, can be well controlled during the design and manufacture of the composite. The composite is assumed to be homogeneous and orthotropic macroscopically. With a selected representative unit cell and the stiffness model developed by author in 2000, the influence of all of these parameters can be determined. Results showing the influence of the main design geometric parameters are presented. They demonstrate that an optimal design is possible for the through-the-thickness stiffness of the composites. The methodology used can be generalized to predict the behavior of other kinds of 3D woven structures.     

FINITE ELEMENT MODELLING OF COMPLEX 3D STATIC AND DYNAMIC CRACK PROPAGATION BY EMBEDDING COHESIVE ELEMENTS IN ABAQUS

<正>This study proposes an algorithm of embedding cohesive elements in Abaqus and develops the computer code to model 3D complex crack propagation in quasi-brittle materials in a relatively easy and efficient manner.The cohesive elements with softening traction-separation relations and damage initiation and evolution laws are embedded between solid elements in regions of interest in the initial mesh to model potential cracks.The initial mesh can consist of tetrahedrons,wedges,bricks or a mixture of these elements.Neither remeshing nor objective crack propagation criteria are needed.Four examples of concrete specimens,including a wedge-splitting test,a notched beam under torsion,a pull-out test of an anchored cylinder and a notched beam under impact,were modelled and analysed.The simulated crack propagation processes and load-displacement curves agreed well with test results or other numerical simulations for all the examples using initial meshes with reasonable densities.Making use of Abaqus's rich pre/post-processing functionalities and powerful standard/explicit solvers,the developed method offers a practical tool for engineering analysts to model complex 3D fracture problems.     

Study on spline wavelet finite-element method in multi-scale analysis for foundation

A new finite element method （FEM） of B-spline wavelet on the interval （BSWI） is proposed. Through analyzing the scaling functions of BSWI in one dimension, the basic formula for 2D FEM of BSWI is deduced. The 2D FEM of 7 nodes and 10 nodes are constructed based on the basic formula. Using these proposed elements, the multiscale numerical model for foundation subjected to harmonic periodic load, the foundation model excited by external and internal dynamic load are studied. The results show the pro- posed finite elements have higher precision than the tradi- tional elements with 4 nodes. The proposed finite elements can describe the propagation of stress waves well whenever the foundation model excited by extemal or intemal dynamic load. The proposed finite elements can be also used to con- nect the multi-scale elements. And the proposed finite elements also have high precision to make multi-scale analysis for structure.     

A PREDICTIVE APPROACH TO THE IN-PLANE MECHANICAL PROPERTIES OF STITCHED COMPOSITE LAMINATES

This contribution attempts to model the alteration of the in-plane elastic properties in laminates caused by stitching,and to predict the in-plane effective tensile strength of the stitched composite laminates.The distortion of in-plane fibers is considered to be the main cause that affects the in-plane mechanical properties.A fiber distortion model is proposed to characterize the fiber misalignment and the fiber content concentration due to stitching.The undistorted region,the fiber distortion region,the resin-rich pocket and the through-thickness reinforcement section are taken into account.The fiber misalignment and inhomogeneous fiber content due to stitching have been formulated by introducing two parameters,the distortion width and maximum misalignment.It has been found that the ply stress concentration in stitched laminates is influenced by the two concurrent factors,the stitch hole and inhomogeneous fiber content.The stitch hole brings about the stress concentration whereas the higher fiber content at the local region induced by stitching restrains the local deformation of the composite.The model is used to predict the tensile strength of the [0/45/0/-45/90/45/0/-45]_(2s) T300/QY9512 composite laminate stitched by Kevlar 29 yarn with different stitching configurations,showing an acceptable agreement with experimental data.     

A THEORETICAL MODEL FOR THERMAL-SENSITIVE MICROGEL WITH PNIPAM CORE AND ELASTIC SHELL

Poly(N-isopropylacrylamide)(PNIPAM) microgels are widely used in drug delivery due to their fast response to temperature.In order to get a better biocompatibility,PNIPAM microgels are typically coated with a layer of biocompatible material,resulting in composite microgels with core-shell structure.In a composite microgel prepared recently,for example,a microsphere of PNIPAM gel is enclosed by a phospholipid membrane,and the composite microgel exhibits a substantial volume transition in response to temperature changes.Here we develop a theoretical model to describe the thermal-responsive behavior of this composite microgel.In particular,we treat the phospholipid membrane as an elastic layer behaving like rubber-like elastomers and adopt the form of the free-energy function for nematic gels(which refer to anther species of thermalsensitive gels whose behavior has been intensively studied) as that for PNIPAM gels.We show that the thermal-responsive behavior of the composite microgel can be markedly influenced by the membrane.By investigating the state of stress on the interface,we further predict that when the coating membrane is stiff and thin,wrinkles are expected to occur on the outer surface of the composite microgel after the volume transition.     

FREE VIBRATION OF ANISOTROPIC RECTANGULAR PLATES BY GENERAL ANALYTICAL METHOD

According to the differential equation for transverse displacement function of anisotropic rectangular thin plates in free vibration, a general analytical solution is established. This general solution, composed of the composite solutions of trigonometric function and hyperbolic function, can satisfy the problem of arbitrary boundary conditions along four edges. The algebraic polynomial with double sine series solutions can also satisfy the problem of boundary conditions at four corners. Consequently, this general solution can be used to solve the vibration problem of anisotropic rectangular plates with arbitrary boundaries accurately. The integral constants can be determined by boundary conditions of four edges and four corners. Each natural frequency and vibration mode can be solved by the determinate of coefficient matrix from the homogeneous linear algebraic equations equal to zero. For example, a composite symmetric angle ply laminated plate with four edges clamped has been calculated and discussed.     

碳纳米管/碳纤维增强复合材料层合板低速冲击响应和破坏的数值模拟

碳纳米管/碳纤维增强复合材料（carbon nanotube/carbon fibre reinforced plastic,CNT/CFRP）是一种多尺度复合材料,比传统CFRP有更好的综合性能和更广阔的应用前景。对CNT/CFRP在低速冲击下的响应和破坏进行了数值模拟研究。首先,基于先前的研究通过引入基体增韧因子、残余强度因子并改进损伤耦合方程,建立了新的FRP动态渐进损伤模型;然后,利用新建立的本构模型并结合黏结层损伤模型,对4种碳纳米管含量的增韧碳纤维增强树脂基复合材料层合板在5个能量下的冲击实验进行了数值模拟;最后,将模拟结果与文献中的相关实验结果进行了比较,并讨论了冲击速度的影响。结果表明:新建立的FRP本构模型能够预测CNT/CFRP层合板在低速冲击载荷作用下的响应、破坏过程和分层形貌,模拟得到的载荷-位移曲线和破坏形貌与实验吻合较好;冲击速度会影响CNT/CFRP层合板拉伸和压缩破坏的比例,相同的冲击能量下,更大的冲击速度会造成更多的拉伸破坏。     

On the description of anisotropic damage in composite laminates

A general anisotropic damage theory of cracked laminates is formulated here. The deformation of composite laminates is composed of matrix elastic strains, pseudo-elastic damage strains due to cracking and permanent damage strains due to interlaminar slip. The surface of damage initiation is constructed according to the concept of linear elastic fracture mechanics for the virgin material. After the initial damage, a pesudo-elastic damage can be used to describe the damage behaviour if interlaminar slip is negligible. Damage evolution, load induced anisotropy and interlaminar intralaminar interaction for composite laminates are examined; the latter can perturb the normality structure of damage strain rate. Explicit expressions are given for pseudo-elastic (or secant) moduli of the damaging composite laminates, under a non-interacting assumption imposed on the cracks between different families.     

复合材料层合板低速冲击响应和损伤分析

通过在有限元软件Abaqus/Explicit中编写用户材料子程序VUMAT，建立了一种基于能量演化的复合材料低速冲击渐进损伤模型．该模型以三维Hashin准则来预测层内损伤的起始，以一种简化的损伤变量来模拟层内损伤的演化，将具有双线性牵引-分离本构的零厚度界面单元插入层间来模拟分层损伤．采用该模型对14.7 J冲击能量下的纤维增强复合材料低速冲击过程进行了仿真分析，计算得到的冲击力-时间曲线、能量-时间曲线和损伤分布与试验结果吻合较好．根据数值模拟结果，分析了纤维、基体和分层损伤的扩展规律，为低速冲击下的复合材料结构设计提供了理论依据．     

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

针对复合材料层合板低速冲击损伤问题,提出了一种各向异性材料连续介质损伤力学模型,模型涵盖损伤表征、损伤起始判定和损伤演化法则3 个方面. 通过材料断裂面坐标下的损伤状态变量矩阵完成损伤表征,并考虑断裂面角度的影响,建立了主轴坐标系下的材料损伤本构关系. 损伤起始由卜克(Puck) 失效准则预测,损伤演化由断裂面上的等效应变控制,服从基于材料应变能释放的线性软化行为. 模型区分了纤维损伤和基体损伤,并根据冲击载荷下层内产生多条基体裂纹继而扩展至界面形成层间裂纹(分层) 的试验观察,引入基体裂纹饱和密度参数表征层间分层. 以[0₃/45/-45]_S 和[45/0/-45/90]_4S 两种铺层的复合材料层合板为例,预测了不同冲击能量下复合材料层合板的低速冲击损伤响应参数,试验结果证明了连续介质损伤力学模型的有效性.模型在不同网格密度下的计算结果表明单元特征长度的引入可以在一定程度上降低损伤演化阶段对网格密度的依赖性.      

纤维/金属网增强层板低速冲击损伤机理和演化

本文首先通过落锤低速冲击实验测试了纯玻璃纤维增强环氧树脂复合材料和304不锈钢丝网(SSWM)/玻璃纤维混杂复合材料的力学性能，探究了SSWM嵌入数量对混杂复合材料抗冲击性能的影响．随后采用Abaqus有限元软件建立了混杂复合材料的低速冲击模型，分别采用三维Hashin失效准则和Jason-Cook破坏准则模拟了纤维/基体和SSWM的损伤；建立了基于表面接触的内聚力模型来模拟界面分层；编写了VUMAT用户子程序定义混杂复合材料层合板的渐进失效过程．结果表明：相较于纯玻璃纤维增强环氧树脂层合板，SSWM/玻璃纤维混杂增强环氧树脂层合板的抗冲击性能更优，其中铺层形式为铺层III的混杂复合材料抗冲击性能最佳．通过对比发现有限元仿真结果与实验结果吻合良好，表明建立的模型适用于SSWM/玻璃纤维混杂增强环氧树脂复合材料低速冲击损伤的评估．通过分析仿真结果发现混杂复合材料的低速冲击损伤主要是冲击区域的纤维断裂、基体破坏和层间分层；SSWM通过吸收和传递冲击能量从而提升了混杂复合材料的抗冲击性能．     

低速冲击下纤维/金属层合板抗分层性能研究

本文采用内聚力模型,对纤维/金属层合板(FMLs)在低速冲击载荷作用下抗分层性能进行研究。内聚力模型对裂纹的模拟具有它独特的优势：一是该模型不需要预先假设初始缺陷;二是在计算过程中随着裂纹的扩展,该方法不需要重新对结构进行网格划分。借助该模型,本文对低速冲击载荷作用下,纤维层合板(FRP)分层进行了模拟,并验证了该模型计算的有效性。在此基础上,本文研究了低速冲击载荷作用下,不同金属含量的纤维/金属层合板抗分层性能,并与纤维层合板进行了比较。最后从能量的角度讨论了金属含量与铺层结构对FMLs低速冲击性能的影响。     

粘结单元在模拟FRP层合板低速冲击响应中的应用

纤维增强树脂基复合材料层合板（fibre reinforced plastic composites，FRP）在航空、航天、交通、造船等诸多工程中得到了日益广泛的应用，而其在冲击载荷下的响应和破坏特别是分层一直为学术界所关注。本文中对FRP层合板在冲击载荷下的响应和破坏进行数值模拟，并通过引入粘结层重点研究其分层破坏。首先，介绍一种基于改进的粘结区域方法的粘结层损伤模型；其次，详细介绍了有限元模型建模过程和建模细节；最后，对有限元模型进行验证，并分析分层损伤发生的原因。模拟结果表明，该模型不仅能准确预测FRP层合板在低速冲击载荷下的载荷-时间曲线和载荷-位移曲线，还能成功地预测其分层破坏。     

含冲击损伤高强中模碳纤维复合材料层压板压缩剩余强度分析与试验验证

论文以碳纤维复合材料层压板为研究对象,发展了一种模拟复合材料层压板冲击及冲击后压缩的一体化数值分析方法．基于Puck 失效准则和粘聚区模型描述层内损伤与层间损伤,分别采用基于断裂能的双线性型、函数型以及直接折减型等不同损伤折减方法构建了层内损伤预测与演化模型;建立了碳纤维复合材料冲击后压缩数值仿真模型,通过开展不同能量冲击后压缩试验,验证了所发展的数值分析方法的有效性;研究结果表明,采用Puck 失效准则和基于断裂能的双线性损伤演化模型预测冲击后压缩强度时具有较高精度．     

材料强度对层板抗弹性能和损伤特性的影响

根据纤维增强复合材料宏细观结构特征,基于ABAQUS软件平台,建立了层合板高速冲击损伤三维有限元分析模型。该模型在复合材料层间引入界面单元模拟层间分层,采用三维粘弹性本构,结合Hashin失效准则模拟单层板面内损伤.利用该模型,深入研究了复合材料层板的抗弹性能和损伤特性,数值分析结果与实验结果吻合良好,证明了该方法的合理有效性。通过数值分析,详细探讨了材料强度参数对层板抗弹性能和损伤特性的影响规律,获得了一些有价值的结论。     

DAMAGE PROGRESSIVE MODEL OF COMPRESSION OF COMPOSITE LAMINATES AFTER LOW VELOCITY IMPACT

IntroductionCompressivepropertiesofcompositelaminatesafterlowvelocityimpactareoneofthemostseriouscircumstanceswhichmustbetakenintoaccountindamagetolerancedesignofcompositestructures[1].Impactdamagegenerallymanifestsintheformsofdelaminations,matrixcracksandfibrefracture[2 ,3].Sofaralotofresearcheshavebeentakeninlowvelocityimpactbehaviorofcompositelaminatesandtheirpost_impactcompressiveproperties.Compressivefailuremechanismsofcompositelaminatesafterlowvelocityimpacthavenotbeenknownclearly .Ther…