Fatigue Damage of Composite Structures Applying a Micromechanical Approach

Fatigue damage calculations of unidirectional polymer composites is presented applying micromechanics theory. An orthotropic micromechanical damage model is integrated with an isotropic fatigue evolution model to predict the micromechanical fatigue damage of the composite structure. The orthotropic micromechanical damage model is used to predict the orthotropic damage evolution within a single cycle. The isotropic fatigue model is used to predict the magnitude of fatigue damage accumulated as a function of the number of cycles. The advantage of using this approach is the cheap determination of model parameters since the orthotropic damage model parameters can be determined using available data from quasi-static loading tests. Decomposition of the state variables down to the constituent scale is accomplished by micromechanics theory. Phenomenological damage evolution models are then postulated for each constituent and for interphase among them. Comparison between model predictions and experimental data is presented.

A micromechanics approach for damage modeling of polymer matrix composites

A new model for damage evolution in polymer matrix composites is presented. The model is based on a combination of two constituent-level models and an interphase model. This approach reduces the number of empirical parameters since the two constituent- level models are formulated for isotropic materials, namely fiber and matrix. Decomposition of the state variables down to the micro-scale is accomplished by micromechanics. Phenomenological damage evolution models are then postulated for each constituent. Determination of material parameters is made from available experimental data. The required experimental data can be obtained with standard tests. Comparison between model predictions and additional experimental data is presented.

A micro-mechanics damage approach for fatigue of composite materials

A new model for fatigue damage evolution of polymer matrix composites (PMC) is presented. The model is based on a combination of an orthotropic damage model and an isotropic fatigue evolution model. The orthotropic damage model is used to predict the orthotropic damage evolution within a single cycle. The isotropic fatigue model is used to predict the magnitude of fatigue damage accumulated as a function of the number of cycles. This approach facilitates the determination of model parameters since the orthotropic damage model parameters can be determined from available data from quasi-static-loading tests. Then, limited amount of fatigue data is needed to adjust the fatigue evolution model. The combination of these two models provides a compromise between efficiency and accuracy. Decomposition of the state variables down to the constituent scale is accomplished by micro-mechanics. Phenomenological damage evolution models are then postulated for each constituent and for the micro-structural interaction among them. Model parameters are determined from available experimental data. Comparison between model predictions and additional experimental data is presented.

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.

Distribution-based adaptive bounding genetic algorithm for continuous optimisation problems

A new search-space-updating technique for genetic algorithms is proposed for continuous optimisation problems. Other than gradually reducing the search space during the evolution process with a fixed reduction rate set ‘a priori’, the upper and the lower boundaries for each variable in the objective function are dynamically adjusted based on its distribution statistics. To test the effectiveness, the technique is applied to a number of benchmark optimisation problems in comparison with three other techniques, namely the genetic algorithms with parameter space size adjustment (GAPSSA) technique [A.B. Djurišic, Elite genetic algorithms with adaptive mutations for solving continuous optimization problems – application to modeling of the optical constants of solids, Optics Communications 151 (1998) 147–159], successive zooming genetic algorithm (SZGA) [Y. Kwon, S. Kwon, S. Jin, J. Kim, Convergence enhanced genetic algorithm with successive zooming method for solving continuous optimization problems, Computers and Structures 81 (2003) 1715–1725] and a simple GA. The tests show that for well-posed problems, existing search space updating techniques perform well in terms of convergence speed and solution precision however, for some ill-posed problems these techniques are statistically inferior to a simple GA. All the tests show that the proposed new search space update technique is statistically superior to its counterparts.

A pseudo-transient solution strategy for the analysis of delamination by means of interface elements

Recent efforts in the finite element modelling of delamination have concentrated on the development of cohesive interface elements. These are characterised by a bilinear constitutive law, where there is an initial high positive stiffness until a threshold stress level is reached, followed by a negative tangent stiffness representing softening (or damage evolution). Complete decohesion occurs when the amount of work done per unit area of crack surface is equal to a critical strain energy release rate. It is difficult to achieve a stable, oscillation-free solution beyond the onset of damage, using standard implicit quasi-static methods, unless a very refined mesh is used. In the present paper, a new solution strategy is proposed based on a pseudo-transient formulation and demonstrated through the modelling of a double cantilever beam undergoing Mode I delamination. A detailed analysis into the sensitivity of the user-defined parameters is also presented. Comparisons with other published solutions using a quasi-static formulation show that the pseudo-transient formulation gives improved accuracy and oscillation-free results with coarser meshes

Microstructure and microtexture evolution of aluminum alloy 3003 under ultrasonic welding process for embedding SiC fibre

Ultrasonic welding process can be used for bonding metal foils which is the fundament of ultrasonic consolidation (UC). UC process can be used to embed reinforcement fibres such as SiC fibres within an aluminum matrix materials. In this research we are investigating the phenomena occurring in the microstructure of the parts during ultrasonic welding process to obtain better understanding about how and why the process works. High-resolution electron backscatter diffraction (EBSD) is used to study the effects of the vibration on the evolution of microstructure in AA3003. The inverse pole figures (IPF) and the correlated misorientation angle distribution of the mentioned samples are obtained. The characteristics of the crystallographic orientation, the grain structure and the grain boundary are analyzed to find the effect of ultrasonic vibration on the microstructure and microtexture of the bond. The ultrasonic vibration will lead to exceptional refinement of grains to a micron level along the bond area and affect the crystallographic orientation. Ultrasonic vibration results in a very weak texture. Plastic flow occurs in the grain after welding process and there is additional plastic flow around the fibre which leads to the fibre embedding. © 2009 Editorial Board of CHINA WELDING.

3D FEM simulation of helical milling hole process for titanium alloy Ti-6Al-4V

As an emerging hole-machining methodology, helical milling process has become increasingly popular in aeromaterials manufacturing research, especially in areas of aircraft structural parts, dies, and molds manufacturing. Helical milling process is highly demanding due to its complex tool geometry and the progressive material failure on the workpiece. This paper outlines the development of a 3D finite element model for helical milling hole of titanium alloy Ti-6Al-4V using commercial FE code ABAQUS/Explicit. The proposed model simulates the helical milling hole process by taking into account the damage initiation and evolution in the workpiece material. A contact model at the interface between end-mill bit and workpiece has been established and the process parameters specified. Furthermore, a simulation procedure is proposed to simulate different cutting processes with the same failure parameters. With this finite element model, a series of FEAs for machined titanium alloy have been carried out and results compared with laboratory experimental data. The effects of machining parameters on helical milling have been elucidated, and the capability and advantage of FE simulation on helical milling process have been well presented.

Modelling the nonlinear behaviour and fracture process of AS4/PEKK thermoplastic composite under shear loading

The accurate determination of non-linear shear behaviour and fracture toughness of continuous carbon-fibre/polymer composites remains a considerable challenge. These measurements are often necessary to generate material parameters for advanced computational damage models. In particular, there is a dearth of detailed shear fracture toughness characterisation for thermoplastic composites which are increasingly generating renewed interest within the aerospace and automotive sectors. In this work, carbon fibre (AS4)/ thermoplastic Polyetherketoneketone (PEKK) composite V-notched cross-ply specimens were manufactured to investigate their non-linear response under pure shear loading. Both monotonic and cyclic loading were applied to study the shear modulus degradation and progressive failure. For the first time in the reported literature, we use the essential work of fracture approach to measure the shear fracture toughness of continuous fibre reinforced composite laminates. Excellent geometric similarity in the load-displacement curves was observed for ligament-scaled specimens. The laminate fracture toughness was determined by linear regression, of the specific work of fracture values, to zero ligament thickness, and verified with computational models. The matrix intralaminar fracture toughness (ply level fracture toughness), associated with shear loading was determined by the area method. This paper also details the numerical implementation of a new three-dimensional phenomenological model for carbon fibre thermoplastic composites using the measured values, which is able to accurately represent the full non-linear mechanical response and fracture process. The constitutive model includes a new non-linear shear profile, shear modulus degradation and load reversal. It is combined with a smeared crack model for representing ply-level damage initiation and propagation. The model is shown to accurately predict the constitutive response in terms of permanent plastic strain, degraded modulus as well as load reversal. Predictions are also shown to compare favourably with the evolution of damage leading to final fracture.

Predicting impact damage, residual strength and crashworthiness of composite structures

The development of the latest generation of wide-body carbon-fibre composite passenger aircraft has heralded a new era in the utilisation of these materials. The premise of superior specific strength and stiffness, corrosion and fatigue resistance, is tempered by high development costs, slow production rates and lengthy and expensive certification programmes. Substantial effort is currently being directed towards the development of new modelling and simulation tools, at all levels of the development cycle, to mitigate these shortcomings. One of the primary challenges is to reduce the extent of physical testing, in the certification process, by adopting a ‘certification by simulation’ approach. In essence, this aspirational objective requires the ability to reliably predict the evolution and progression of damage in composites. The aerospace industry has been at the forefront of developing advanced composites modelling tools. As the automotive industry transitions towards the increased use of composites in mass-produced vehicles, similar challenges in the modelling of composites will need to be addressed, particularly in the reliable prediction of crashworthiness. While thermoset composites have dominated the aerospace industry, thermoplastics composites are likely to emerge as the preferred solution for meeting the high-volume production demands of passenger road vehicles. This keynote presentation will outline recent progress and current challenges in the development of finite-element-based predictive modelling tools for capturing impact damage, residual strength and energy absorption capacity of thermoset and thermoplastic composites for crashworthiness assessments.

Assortative mating as a mechanism for rapid evolution of a migratory divide.

There have been numerous recent observations of changes in the behavior and dynamics of migratory bird populations, but the plasticity of the migratory trait and our inability to track small animals over large distances have hindered investigation of the mechanisms behind migratory change. We used habitat-specific stable isotope signatures to show that recently evolved allopatric wintering populations of European blackcaps Sylvia atricapilla pair assortatively on their sympatric breeding grounds. Birds wintering further north also produce larger clutches and fledge more young. These findings describe an important process in the evolution of migratory divides, new migration routes, and wintering quarters. Temporal segregation of breeding is a way in which subpopulations of vertebrates may become isolated in sympatry.