Multi-scale ballistic material modeling of cross-plied compliant composites

The open-literature material properties for fiber and polymeric matrix, unit-cell microstructural characteristics, atomic-level simulations and unit-cell based finite-element analyses are all used to construct a new continuum-type ballistic material model for 0°/90° cross-plied highly-oriented polyethylene fiber-based armor-grade composite laminates. The material model is formulated in such a way that it can be readily implemented into commercial finite-element programs like ANSYS/Autodyn [ANSYS/Autodyn version 11.0, User Documentation, Century Dynamics Inc. a subsidiary of ANSYS Inc. (2007)] and ABAQUS/Explicit [ABAQUS version 6.7, User Documentation, Dessault Systems, 2007] as a User Material Subroutine. Model validation included a series of transient non-linear dynamics simulations of the transverse impact of armor-grade composite laminates with two types of projectiles, which are next compared with their experimental counterparts. This comparison revealed that a reasonably good agreement is obtained between the experimental and the computational analyses with respect to: (a) the composite laminates’ capability, at different areal densities, to defeat the bullets with different impact velocities; (b) post-mortem spatial distribution of damage within the laminates; (c) the temporal evolution of composite armor laminate back-face bulging and delamination; and (d) the existence of three distinct penetration stages (i.e. an initial filament shearing/cutting dominated stage, an intermediate stage characterized by pronounced filament/matrix de-bonding/decohesion and the final stage associated with the extensive back-face delamination and bulging of the armor panel).     

Acoustic emission analysis of composite pressure vessels under constant and cyclic pressure

The use of acoustic emission (AE) for the detection of damage in carbon fibre composite pressure vessels was evaluated for constant and cyclic internal gas pressure loading conditions. AE was capable of monitoring the initiation and accumulation of damage events in a composite pressure vessel (CPVs), although it was not possible to reliably distinguish carbon fibre breakage from other microscopic damage events (e.g. matrix cracks, fibre/matrix interfacial cracks). AE tests performed on the carbon fibre laminate used as the skin of pressure vessels revealed that the development of damage is highly variable under constant pressure, with large differences in the rupture life and acoustic emission events at final failure. Numerical analysis of the skin laminate under constant tensile stress revealed that the high variability in the stress rupture life is due mainly to the stochastic behaviour of the carbon fibre rupture process.     

Finite element modeling of post-fire flexural modulus of fiber reinforced polymer composites under constant heat flux

In this study, a simple 1D finite element model was developed to predict the temperature evolution and post-fire mechanical degradation of glass fiber reinforced polymers (FRPs) subjected to constant heat fluxes, including 35 kW/m², 50 kW/m², 75 kW/m², and 100 kW/m². A temperature-dependent post-fire mechanical property model was proposed and implemented. The calculated temperature and residual mechanical moduli showed good agreement with the experimental data. By properly selecting the parameters of the model, an effective strategy was demonstrated to design FRP structure with enhanced durability.     

A mechanism-based multi-scale model for predicting thermo-oxidative degradation in high temperature polymer matrix composites

This paper describes a mechanism-based multi-scale model for life prediction of high temperature polymer matrix composites (HTPMC) under thermo-oxidative aging conditions. The multi-scale model incorporates molecular level damage such as inter-crosslink chain scission in a thermoset polymer due to thermo-oxidative aging of the polymer resin. The degradation of inter-laminar stress depends on remaining inter-crosslink density of thermo-set polymer in fiber/matrix interface region subjected to thermo-oxidative aging environment. The degradation of inter-laminar shear stress of thermo-oxidatively aged unidirectional IM-7/PETI-5 composite specimens at 300 °C was modeled using an in-house test-bed FEA code (NOVA-3D). A micromechanics based viscoelastic cohesive layer model was used to model delamination. The model is fully rate dependent and does not require a pre-assigned traction-separation law. Viscoelastic regularization of the constitutive equations of the cohesive layer used in this model not only mitigates numerical instability, but also enables the analysis to follow load-deflection behavior beyond peak failure load. The model was able to successfully simulate delamination failure in thermo-oxidatively aged unidirectional IM-7/PETI-5 composite, and the model predictions were verified using test data.     

Predicting dynamic in-plane compressive properties of multi-axial multi-layer warp-knitted composites with a meso-model

Proper prediction of material microstructure from known processing conditions and constituent material properties is a critical step to determine the bulk properties of the composite. This paper reports a meso-structure model of multi-axial multi-layer warp-knitted (MMWK) composites from an elastic–plastic material model considering the strain rate effect for the components of the MMWK composite. The representative unit cell (RUC) of fiber tow is created to obtain the elastic–plastic parameters of the fiber tow. The 3D meso-structure model of the MMWK composite is based on an idealized geometrical model according to the preform structure of the MMWK fabric. The model is used to investigate the effect of the volume fraction of the knitting yarn on the dynamic in-plane compressive properties. The results show that the fiber tow failure at large extent is mainly caused by the micro cracking of the matrix, and the effects of the knitting yarn on the mechanical properties of MMWK composite are very limited. Particularly, MMWK composites could be considered as laminates when the volume fraction of the knitting yarn is low, such as below 1.5%. Experiments were also conducted to validate the results from the simplified meso-structure model of the MMWK composite. The material is found to be strain rate sensitive, and the experimental and predicted results agree well with respect to the compressive strength and modulus of the composite. This confirms that the meso-structure MMWK composite model proposed is capable of capturing the essential features for the response of the composite under different strain rate conditions at the meso-level.     

Effect of adhesive fillet geometry on bond strength between microelectronic components and composite circuit boards

Printed circuit boards (PCBs) assembled with ball grid array (BGA) microelectronics packages were tested in a double cantilever beam (DCB) configuration. The results were compared for a filled and an unfilled underfill epoxy adhesive as well as a cyanoacrylate adhesive. The original fillet, formed in the underfilling process, was modified to create fillets of different sizes. Regardless of the underfill thermal and mechanical properties as well as its curing profile, the crack initiation load and the failure mode were solely a function of the size of the underfill fillet, and the failure always initiated within the PCB. Moreover, the strength of the underfilled solder joints was increased significantly (approximately 100%) by the presence of a relatively large fillet. This effect of the underfill fillet on the crack path and the fracture load was then examined in terms of differences in the stress states using a finite element model.     

A combined optical-thermal model for near-infrared laser heating of thermoplastic composites in an automated tape placement process

A combined optical-thermal model is presented for near-infrared laser heating of carbon fibre reinforced thermoplastic composites in an automated tape placement process. For the first time a three dimensional ray tracing model is presented for a near-infrared laser tape placement process which captures the unique anisotropic scattering behaviour of the composite. Predicted irradiance distributions on the composite are subsequently applied to a 2D non-linear finite element thermal model. It is shown that a shadow is present in the process and causes a significant drop in temperature prior to the consolidation zone. The effect of various source and surface model simplifications was also studied. The modelled temperature profiles agree well with experimental data. Substrate fibre orientation was investigated and found to have only a small influence on the temperature history.     

A progressive damage simulation algorithm for GFRP composites under cyclic loading. Part I: Material constitutive model

A life prediction algorithm and its implementation for a thick-shell finite element formulation for GFRP composites under constant or variable amplitude loading is introduced in this work. It is a distributed damage model in the sense that constitutive material response is defined in terms of meso-mechanics for the unidirectional ply. The algorithm modules for non-linear material behaviour, pseudo-static loading-unloading-reloading response, Constant Life Diagrams and strength and stiffness degradation due to cyclic loading were implemented on a robust and comprehensive experimental database for a unidirectional glass/epoxy ply. The model, based on property definition in the principal coordinate system of the constitutive ply, can be used, besides life prediction, to assess strength and stiffness of any multidirectional laminate after arbitrary, constant or variable amplitude multi-axial cyclic loading. Numerical predictions were corroborated satisfactorily by test data from constant amplitude fatigue of glass/epoxy laminates of various stacking sequences.     

Nondestructive evaluation of forced delamination in glass fiber-reinforced composites by terahertz and ultrasonic waves

Glass fiber-reinforced composite laminates in polyetherimide resin have been studied via terahertz imaging and ultrasonic C-scans. The forced delamination is created by inserting Teflon film between various layers inside the samples prior to consolidating the laminates. Using reflective pulsed terahertz imaging, we find high-resolution, low-artifact terahertz C-scan and B-scan images locating and sizing the delamination in three dimensions. Furthermore, terahertz imaging enables us to determine the thicknesses of the delamination and of the layers constituting the laminate. Ultrasonic C-scan images are also successfully obtained; however, in our samples with small thickness-to-wavelength ratio, detailed ultrasonic B-scan images providing quantitative information in depth cannot be obtained by 5 MHz or 10 MHz focused transducers. Comparative analysis between terahertz imaging and ultrasonic C-scans with regard to spatial resolution is carried out demonstrating that terahertz imaging provides higher spatial resolution for imaging, and can be regarded as an alternative or complementary modality to ultrasonic C-scans for this class of glass fiber-reinforced composites.     

Deformation gradient tensor decomposition for representing matrix cracks in fiber-reinforced materials

A new method is presented for the representation of matrix cracks in continuum damage mechanics (CDM) models for fiber-reinforced materials. The method is based on the additive decomposition of the deformation gradient tensor into ‘crack’ and ‘bulk material’ components, analogous to the additive strain decomposition of the smeared-crack approach. The potential improvements to the accuracy of CDM models that utilize the presented method are demonstrated for a single element subjected to simple shear deformation and for a unidirectional open-hole tension specimen. The presented method avoids load transfer across matrix cracks and eliminates the prediction of spurious secondary failure modes that occurs when conventional strain-based CDM models are used in geometrically nonlinear finite element analyses involving large shear deformations.     

Hemp fiber composites for the design of a Naca cowling for ultra-light aviation

The present study focuses on the design of a Naca cowling of an acrobatic ultra-light airplane, where the traditional woven glass/epoxy laminate utilized for the production has been replaced by woven hemp reinforced epoxy composite. Specialized software (Fluent and ANSYS) was used for configuration, design and analysis. The results showed comparable mechanical performance, about same weight, but easier disposal and better eco-friendly characteristics as compared to their synthetic counterpart. The engine cover, produced by the use of hemp/epoxy composites, demonstrates the effective possibility to produce semi-structural aeronautical components using natural fiber composites in substitution of glass ones.     

Bending strength of adhesive joints in microelectronic components: Comparison of edge-bonding and underfilling

The bending strength of underfilled and edge-bonded ball grid array (BGA) microelectronic packages assembled on printed circuit boards (PCBs) was compared using double cantilever beam (DCB) specimens. All specimens with fillets of the same size and shape failed at the same load, with cracks initiating and propagating within the PCB. This was consistent with measurements of the crack initiation strain energy release rate for PCB interfacial failure, which was significantly smaller than that of cohesive failure within the adhesives. Finite element analysis (FEA) indicated that the stress state in the PCB near the PCB-fillet interface in both underfilled and edge-bonded specimens was only a function of the adhesive fillet size and shape, and independent of the extent of the adhesive layer between the PCB and the BGA, and independent of the adhesive mechanical and thermal properties over the broad range of properties of the tested adhesives. This explained why decreasing the fillet curvature in edge-bonded specimens produced a significant increase in the joint strength. The crack path in the PCB of the edge-bonded specimens was found to change with the adhesive cure temperature; however, this had a negligible effect on the failure load.     

Experimental and finite element studies on hot sizing process for L-shaped composite beams

This work aims at developing a hot sizing process on composite materials to correct the profiles of composite structures during manufacture. Hot sizing experiments were carried out at 150 °C with different sizing loads and hot sizing periods for L-shaped composite beams made of carbon fiber plain-weave fabric and epoxy resin. To predict the springback in hot sizing process, a corresponding finite element simulation method was developed using stress relaxation equations determined at the same temperature. Excellent agreements between the predicted and observed results were obtained. The effects of the component thickness and 45° ply percentage on the springback rate were investigated by simulation. Springback rate in hot sizing process on composite materials ranges from 60% to 95%. In conclusion hot sizing process is proved to be a valid method for compensation for the process-induced deformation (PID) of L-shaped composite beams.     

Mechanical,flow and electrical properties of thermoplastic polyurethane/fullerene composites: Effect of surface modification of fullerene

Thermoplastic polyurethane (TPU) composites with fullerene loadings varying from 0.5 to 2 weight% were prepared by melt-mixing method. Nitric acid oxidation and silanization were applied to fullerene surface to improve interfacial interactions with TPU matrix. The influence of surface modifications of fullerene on mechanical, melt flow and electrical properties of TPU based composites were investigated. Incorporation of fullerene leads to nearly twofold increase in tensile strength and Young's modulus of the composites in addition to enhancing the flexibility. The best results are obtained in nitric acid and silane modified fullerene containing composites at the lowest concentration (0.5%). Higher MFI values were observed for composites loaded with surface treated fullerenes compared to pristine fullerene because of their better dispersion in TPU. Electrical properties of TPU also improved by the addition of surface modified fullerene particles. Surface oxidation and silanization gave rise to dispersion homogeneity which may be the reason of both tensile strength and strain improvements at the same time.     

Continuous resistance welding of thermoplastic composites: Modelling of heat generation and heat transfer

A process model composed of electrical and heat transfer models was developed to simulate continuous resistance welding of thermoplastic composites. Glass fabric reinforced polyphenylenesulfide welded in a lap-shear configuration with a stainless steel mesh as the heating element was considered for modelling and experimental validation of the numerical results. The welding temperatures predicted by the model showed good agreement with the experimental results. Welding input power and welding speed were found to be the two most important parameters influencing the welding temperature. The contact quality between the electrical connectors and the heating element was found to influence the distribution of the welding temperature transverse to the welding direction. Moreover, the size of the electrical connectors was found to influence the achievable welding speed and required power input for a certain welding temperature.     

A progressive damage simulation algorithm for GFRP composites under cyclic loading. Part II: FE implementation and model validation

The FE implementation of FADAS, a material constitutive model capable of simulating the mechanical behaviour of GFRP composites under variable amplitude multiaxial cyclic loading, was presented. The discretization of the problem domain by means of FE is necessary for predicting the damage progression in real structures, as failure initiates at the vicinity of a stress concentrator, causing stress redistribution and the gradual spread of damage until the global failure of the structure. The implementation of the stiffness and strength degradation models in the principal material directions of the unidirectional ply was thoroughly discussed. Details were also presented on the FE models developed, the computational effort needed and the definition of final failure considered. Numerical predictions were corroborated satisfactorily by experimental data from constant amplitude uniaxial fatigue of multidirectional glass/epoxy laminates under various stress ratios. The validation of predictions included fatigue strength, stiffness degradation and residual static strength after cyclic loading.     

Numerical method for optimizing design variables of carbon-fiber-reinforced epoxy composite coil springs

To successfully reduce a vehicle's weight by replacing steel with composite materials, it is essential to optimize the material parameters and design variables of the structure. In this study, we investigated numerical and experimental methods for determining the ply angles and wire diameters of carbon fiber/epoxy composite coil springs to attain a spring rate equal to that of an equivalent steel component. First, the shear modulus ratio for two materials was calculated as a function of the ply angles and compared with the experimental results. Then, by using the equation of the spring rate with respect to the shear modulus and design variables, normalized spring rates were obtained for specific ply angles and wire diameters. Finally, a finite element model for an optimal composite coil spring was constructed and analyzed to obtain the static spring rate, which was then compared with the experimental results.     

Non-linear analysis of fiber-reinforced open conical shell panels considering variation of thickness and fiber orientation under thermo-mechanical loadings

Non-linear bending analysis of moderately thick laminated conical panels under various thermo-mechanical loadings and boundary conditions is presented using the generalized differential quadrature (GDQ) method together with the Newton–Raphson iterative scheme. The stiffness coefficients are assumed to be functions of the meridional and circumferential coordinates in panels for the realistic applications. In the first case of orthotropic open conical shell panels, the orientation of fibers are assumed to be in the meridional and circumferential directions. The stiffness coefficients of this type of fiber-reinforced panel are usually assumed to be constant. It is shown that due to the geometry of the conical surface, thickness of laminate will be changed along the meridional direction. The effect of stiffness variation on the non-linear response of panel is considered for the first time. In the second type, open conical shell panel can be made by cutting from a filament wound circular conical shell. In this case, thickness and ply orientation are functions of the shell coordinates. In this paper, different path definitions for variable stiffness filament wound shells are considered. The inclusion of this geometric complicating effect in large deformation analysis will add considerably to the complication and cost of a solution scheme. Paper presents some results to show when these assumptions have a significant effect on the end result. Assuming the effects of shear deformation and initial curvature, based on the first-order shear deformation theory (FSDT) and von Kármán-type of geometric non-linearity, the governing system of equations is obtained. Comparisons of the predictions with those available in the literature and finite element analyses show very good agreement. More results for panels with particular boundary conditions and thermo-mechanical load are presented for future references. For the sake of brevity, numerical results which presented in this paper are limited to deflection responses only.     

Elastic constants of metal/ceramic composites with lamellar microstructures: Finite element modelling and ultrasonic experiments

The elastic properties of single domains of lamellar AlSi12/Al₂O₃ composites produced by metal infiltration of freeze cast preforms have been examined. The anisotropic elastic constants determined from ultrasonic phase spectroscopy (UPS) experiments have been compared to microstructure based FE-models created with the program OOF2, micromechanical models (Mori–Tanaka and inverse Mori–Tanaka) and an analytical model for an ideal laminate. The influences of lamellae orientations and ceramic contents on the elastic constants have been investigated. Along the lamellae directions the microstructure based FE model and the inverse Mori–Tanaka model are in good agreement with experimental results. Perpendicular to the lamellae the experiment shows a stiffer than expected elastic behaviour.     

Establishing a new Forming Limit Curve for a flax fibre reinforced polypropylene composite through stretch forming experiments

In developing an understanding of the failure in natural fibre reinforced polymer composites, the failure limits of this class of the material system are required. It is found that the conventional Forming Limit Curve is not suitable to predict the failure initiated in the natural fibre composite as principal strains cannot differentiate the strain on the flax fibres and the polypropylene matrix. This study proposes a new Forming Limit Curve for the composite which expresses limiting fibre strain as a function of forming mode depicted by the ratio of minor strain to major strain. The new Forming Limit Curve, along with the Maximum Strain failure criterion have been successfully implemented in FEA simulations, and numerical simulations suggest that the former is more accurate. The current work provides an innovative method to predict the onset of failure in natural fibre composites, which can be applied in composite forming and structural design.