Non-linearity in piezo-fiber reinforced composites: an asymptotic approach

An asymptotically correct analysis is developed for Macro Fiber Composite unit cell using Variational Asymptotic Method (VAM). VAM splits the 3D nonlinear problem into two parts: A 1D nonlinear problem along the length of the fiber and a linear 2D cross-sectional problem. Closed form solutions are obtained for the 2D problem which are in terms of 1D parameters.

Modeling and mesoscopic damage constitutive relation of brittle short-fiber-reinforced composites

Aimed at brittle composites reinforced by randomly distributed short-fibers with a relatively large aspect ratio, stiffness modulus and strength, a mesoscopic material model was proposed. Based on the statistical description, damage mechanisms, damage-induced anisotropy, damage rate effect and stress redistribution, the constitutive relation were derived. By taking glass fiber reinforced polypropylene polymers as an example, the effect of initial orientation distribution of fibers, damage-induced anisotropy, and damage-rate effect on macro-behaviors of composites were quantitatively analyzed. The theoretical predictions compared favorably with the experimental results.

Theoretical and experimental researches of post-micro-buckling for fiber-reinforced composites

Post-microbuckling is a fundamental feature of compressive failure process for the unidirectional-fiber-reinforced composites and laminated composites. The post-microbuckling behavior of composites under compression in the light of the Kevlar49-reinforced 648/BF3.400 (brittle epoxy) and EP (flexible epoxy) is studied, theoretically and experimentally. Analytical results of compressive strength are in good agreement with experimental results, qualitatively and quantitatively. By the experimental research, the post-microbuckling feature of the advancing kink band model is clearly displayed.

An eigenfunction expansion-variational method in prediction of the transverse thermal conductivity of fiber reinforced composites considering interfacial characteristics

An eigenfunction expansion-variational method based on a unit cell is developed to deal with the steady-state heat conduction problem of doubly-periodic fiber reinforced composites with interfacial thermal contact resistance or coating. The numerical results show a rapid convergence of the present method. The present solution provides a unified first-order approximation formula of the effective thermal conductivity for different interfacial characteristics and fiber distributions. A comparison with the present high-order results, available experimental data and micromechanical estimations demonstrates that the first-order approximation formula is a good engineering closed-form formula. An engineering equivalent parameter reflecting the overall influence of the thermal conductivities of the matrix and fibers and the interfacial characteristic on the effective thermal conductivity, is found. The equivalent parameter can greatly simplify the complicated relation of the effective thermal conductivity to the internal structure of a composite. (c) 2010 Elsevier Ltd. All rights reserved.

Estimation of Transverse Thermal Conductivity of Doubly-periodic Fiber Reinforced Composites

For steady-state heat conduction a new variational functional for a unit cell of composites with periodic microstructures is constructed by considering the quasi-periodicity of the temperature field and in the periodicity of the heat flux fields. Then by combining with the eigenfunction expansion of complex potential which satisfies the fiber-matrix interface conditions, an eigenfunction expansion-variational method (EEVM) based on a unit cell is developed. The effective transverse thermal conductivities of doubly-periodic fiber reinforced composites are calculated, and the first-order approximation formula for the square and hexagonal arrays is presented,which is convenient for engineering application. The numerical results show a good convergency of the presented method, even through the fiber volume fraction is relatively high. Comparisons with the existing analytical and experimental results are made to demonstrate the accuracy and validity of the first-order approximation formula for the hexagonal array.

Explanation for Micromorphologies Around Broken Fibers in Fiber-Reinforced Composites

Experimental observations on micromorphologies around broken fibers in glass-fiber-reinforced epoxy matrix composites reveal different kinds of highly oriented patches at the circumambience of broken fibers, whereas the bulk of the matrix has been observed to be largely isotropic. These patches are interpreted to correlated areas where the stress gradients of the matrix are formed after fiber breaking, but the underlying cause for the orientation is still unknown. The authors have modified an embedded cell model to explain the experimental phenomena. The finite element simulation indicates that the surfaces around broken fibers display a change from an extension micromorphology to a mixed tension and shear micromorphology with the increase of applied strain.

Numerical investigation of the effects of drill geometry on drilling induced delamination of carbon fiber reinforced composites

Drilling is a major process in the manufacturing of holes required for the assemblies of composite laminates in aerospace industry. Simulation of drilling process is an effective method in optimizing the drill geometry and process parameters in order to improve hole quality and to reduce the drill wear. In this research we have developed three-dimensional (3D) FE model for drilling CFRP. A 3D progressive intra-laminar failure model based on the Hashin's theory is considered. Also an inter-laminar delamination model which includes the onset and growth of delamination by using cohesive contact zone is developed. The developed model with inclusion of the improved delamination model and real drill geometry is used to make comparison between the step drill of different stage ratio and twist drill. Thrust force, torque and work piece stress distributions are estimated to decrease by the use of step drill with high stage ratio. The model indicates that delamination and other workpiece defects could be controlled by selection of suitable step drill geometry. Hence the 3D model could be used as a design tool for drill geometry for minimization of delamination in CFRP drilling. © 2013 Elsevier Ltd.

Prediction of failure and fracture mechanisms of polymeric composites using finite element analysis. Part 2:Fiber reinforced composites

A robust finite element scheme for the micro-mechanical modeling of the behavior of fiber reinforced polymeric composites under external loads is developed. The developed model is used to simulate stress distribution throughout the composite domain and to identify the locations where maximum stress concentrations occur. This information is used as a guide to predict dominant failure and crack growth mechanisms in fiber reinforced composites. The differences between continuous fibers, which are susceptible to unidirectional transverse fracture, and short fibers have been demonstrated. To assess the validity and range of applicability of the developed scheme, numerical results obtained by the model are compared with the available experimental data and also with the values found using other methods reported in the literature. These comparisons show that the present finite element scheme can generate meaningful results in the analysis of fiber reinforced composites.

Fracture and resistance-curve behavior in hybrid natural fiber and polypropylene fiber reinforced composites

This article presents the results of a combined experimental and theoretical study of fracture and resistance-curve behavior of hybrid natural fiber- and synthetic polymer fiber-reinforced composites that are being developed for potential applications in affordable housing. Fracture and resistance-curve behavior are studied using single-edge notched bend specimens. The sisal fibers used were examined using atomic force microscopy for fiber bundle structures. The underlying crack/microstructure interactions and fracture mechanisms are elucidated via in situ optical microscopy and ex-situ environmental scanning microscopy techniques. The observed crack bridging mechanisms are modeled using small and large scale bridging concepts. The implications of the results are then discussed for the design of eco-friendly building materials that are reinforced with natural and polypropylene fibers.

Curing, defects and mechanical performance of fiber-reinforced composites

Tradicionalmente, la fabricación de materiales compuestos de altas prestaciones se lleva a cabo en autoclave mediante la consolidación de preimpregnados a través de la aplicación simultánea de altas presiones y temperatura. Las elevadas presiones empleadas en autoclave reducen la porosidad de los componentes garantizando unas buenas propiedades mecánicas. Sin embargo, este sistema de fabricación conlleva tiempos de producción largos y grandes inversiones en equipamiento lo que restringe su aplicación a otros sectores alejados del sector aeronáutico. Este hecho ha generado una creciente demanda de sistemas de fabricación alternativos al autoclave. Aunque estos sistemas son capaces de reducir los tiempos de producción y el gasto energético, por lo general, dan lugar a materiales con menores prestaciones mecánicas debido a que se reduce la compactación del material al aplicar presiones mas bajas y, por tanto, la fracción volumétrica de fibras, y disminuye el control de la porosidad durante el proceso. Los modelos numéricos existentes permiten conocer los fundamentos de los mecanismos de crecimiento de poros durante la fabricación de materiales compuestos de matriz polimérica mediante autoclave. Dichos modelos analizan el comportamiento de pequeños poros esféricos embebidos en una resina viscosa. Su validez no ha sido probada, sin embargo, para la morfología típica observada en materiales compuestos fabricados fuera de autoclave, consistente en poros cilíndricos y alargados embebidos en resina y rodeados de fibras continuas. Por otro lado, aunque existe una clara evidencia experimental del efecto pernicioso de la porosidad en las prestaciones mecánicas de los materiales compuestos, no existe información detallada sobre la influencia de las condiciones de procesado en la forma, fracción volumétrica y distribución espacial de los poros en los materiales compuestos. Las técnicas de análisis convencionales para la caracterización microestructural de los materiales compuestos proporcionan información en dos dimensiones (2D) (microscopía óptica y electrónica, radiografía de rayos X, ultrasonidos, emisión acústica) y sólo algunas son adecuadas para el análisis de la porosidad. En esta tesis, se ha analizado el efecto de ciclo de curado en el desarrollo de los poros durante la consolidación de preimpregnados Hexply AS4/8552 a bajas presiones mediante moldeo por compresión, en paneles unidireccionales y multiaxiales utilizando tres ciclos de curado diferentes. Dichos ciclos fueron cuidadosamente diseñados de acuerdo a la caracterización térmica y reológica de los preimpregnados. La fracción volumétrica de poros, su forma y distribución espacial se analizaron en detalle mediante tomografía de rayos X. Esta técnica no destructiva ha demostrado su capacidad para analizar la microestructura de materiales compuestos. Se observó, que la porosidad depende en gran medida de la evolución de la viscosidad dinámica a lo largo del ciclo y que la mayoría de la porosidad inicial procedía del aire atrapado durante el apilamiento de las láminas de preimpregnado. En el caso de los laminados multiaxiales, la porosidad también se vio afectada por la secuencia de apilamiento. En general, los poros tenían forma cilíndrica y se estaban orientados en la dirección de las fibras. Además, la proyección de la población de poros a lo largo de la dirección de la fibra reveló la existencia de una estructura celular de un diámetro aproximado de 1 mm. Las paredes de las celdas correspondían con regiones con mayor densidad de fibra mientras que los poros se concentraban en el interior de las celdas. Esta distribución de la porosidad es el resultado de una consolidación no homogenea. Toda esta información es crítica a la hora de optimizar las condiciones de procesado y proporcionar datos de partida para desarrollar herramientas de simulación de los procesos de fabricación de materiales compuestos fuera de autoclave. Adicionalmente, se determinaron ciertas propiedades mecánicas dependientes de la matriz termoestable con objeto de establecer la relación entre condiciones de procesado y las prestaciones mecánicas. En el caso de los laminados unidireccionales, la resistencia interlaminar depende de la porosidad para fracciones volumétricas de poros superiores 1%. Las mismas tendencias se observaron en el caso de GIIc mientras GIc no se vio afectada por la porosidad. En el caso de los laminados multiaxiales se evaluó la influencia de la porosidad en la resistencia a compresión, la resistencia a impacto a baja velocidad y la resistencia a copresión después de impacto. La resistencia a compresión se redujo con el contenido en poros, pero éste no influyó significativamente en la resistencia a compresión despues de impacto ya que quedó enmascarada por otros factores como la secuencia de apilamiento o la magnitud del daño generado tras el impacto. Finalmente, el efecto de las condiciones de fabricación en el proceso de compactación mediante moldeo por compresión en laminados unidireccionales fue simulado mediante el método de los elementos finitos en una primera aproximación para simular la fabricación de materiales compuestos fuera de autoclave. Los parámetros del modelo se obtuvieron mediante experimentos térmicos y reológicos del preimpregnado Hexply AS4/8552. Los resultados obtenidos en la predicción de la reducción de espesor durante el proceso de consolidación concordaron razonablemente con los resultados experimentales. Manufacturing of high performance polymer-matrix composites is normally carried out by means of autoclave using prepreg tapes stacked and consolidated under the simultaneous application of pressure and temperature. High autoclave pressures reduce the porosity in the laminate and ensure excellent mechanical properties. However, this manufacturing route is expensive in terms of capital investment and processing time, hindering its application in many industrial sectors. This fact has driven the demand of alternative out-of-autoclave processing routes. These techniques claim to produce composite parts faster and at lower cost but the mechanical performance is also reduced due to the lower fiber content and to the higher porosity. Corrient numerical models are able to simulate the mechanisms of void growth in polymer-matrix composites processed in autoclave. However these models are restricted to small spherical voids surrounded by a viscous resin. Their validity is not proved for long cylindrical voids in a viscous matrix surrounded by aligned fibers, the standard morphology observed in out-of-autoclave composites. In addition, there is an experimental evidence of the detrimental effect of voids on the mechanical performance of composites but, there is detailed information regarding the influence of curing conditions on the actual volume fraction, shape and spatial distribution of voids within the laminate. The standard techniques of microstructural characterization of composites (optical or electron microscopy, X-ray radiography, ultrasonics) provide information in two dimensions and are not always suitable to determine the porosity or void population. Moreover, they can not provide 3D information. The effect of curing cycle on the development of voids during consolidation of AS4/8552 prepregs at low pressure by compression molding was studied in unidirectional and multiaxial panels. They were manufactured using three different curing cycles carefully designed following the rheological and thermal analysis of the raw prepregs. The void volume fraction, shape and spatial distribution were analyzed in detail by means of X-ray computed microtomography, which has demonstrated its potential for analyzing the microstructural features of composites. It was demonstrated that the final void volume fraction depended on the evolution of the dynamic viscosity throughout the cycle. Most of the initial voids were the result of air entrapment and wrinkles created during lay-up. Differences in the final void volume fraction depended on the processing conditions for unidirectional and multiaxial panels. Voids were rod-like shaped and were oriented parallel to the fibers and concentrated in channels along the fiber orientation. X-ray computer tomography analysis of voids along the fiber direction showed a cellular structure with an approximate cell diameter of 1 mm. The cell walls were fiber-rich regions and porosity was localized at the center of the cells. This porosity distribution within the laminate was the result of inhomogeneous consolidation. This information is critical to optimize processing parameters and to provide inputs for virtual testing and virtual processing tools. In addition, the matrix-controlled mechanical properties of the panels were measured in order to establish the relationship between processing conditions and mechanical performance. The interlaminar shear strength (ILSS) and the interlaminar toughness (GIc and GIIc) were selected to evaluate the effect of porosity on the mechanical performance of unidirectional panels. The ILSS was strongly affected by the porosity when the void contents was higher than 1%. The same trends were observed in the case of GIIc while GIc was insensitive to the void volume fraction. Additionally, the mechanical performance of multiaxial panels in compression, low velocity impact and compression after impact (CAI) was measured to address the effect of processing conditions. The compressive strength decreased with porosity and ply-clustering. However, the porosity did not influence the impact resistance and the coompression after impact strength because the effect of porosity was masked by other factors as the damage due to impact or the laminate lay-up. Finally, the effect of the processing conditions on the compaction behavior of unidirectional AS4/8552 panels manufactured by compression moulding was simulated using the finite element method, as a first approximation to more complex and accurate models for out-of autoclave curing and consolidation of composite laminates. The model parameters were obtained from rheological and thermo-mechanical experiments carried out in raw prepreg samples. The predictions of the thickness change during consolidation were in reasonable agreement with the experimental results.

Intraply fracture of fiber-reinforced composites: microscopic mechanisms and modeling

The fracture behavior parallel to the fibers of an E-glass/epoxy unidirectional laminate was studied by means of three-point tests on notched beams. Selected tests were carried out within a scanning electron microscope to ascertain the damage and fracture micromechanisms upon loading. The mechanical behavior of the notched beam was simulated within the framework of the embedded cell model, in which the actual composite microstructure was resolved in front of the notch tip. In addition, matrix and interface properties were independently measured in situ using a nanoindentor. The numerical simulations very accurately predicted the macroscopic response of the composite as well as the damage development and crack growth in front of the notch tip, demonstrating the ability of the embedded cell approach to simulate the fracture behavior of heterogeneous materials. Finally, this methodology was exploited to ascertain the influence of matrix and interface properties on the intraply toughness.

Application of digital image correlation at the microscale in fiber-reinforced composites

Digital image correlation (DIC) is applied to analyzing the deformation mechanisms under transverse compression in a fiber-reinforced composite. To this end, compression tests in a direction perpendicular to the fibers were carried out inside a scanning electron microscope and secondary electron images obtained at different magnifications during the test. Optimum DIC parameters to resolve the displacement and strain field were computed from numerical simulations of a model composite and they were applied to micrographs obtained at different magnifications (250_, 2000_, and 6000_). It is shown that DIC of low-magnification micrographs was able to capture the long range fluctuations in strain due to the presence of matrix-rich and fiber-rich zones, responsible for the onset of damage. At higher magnification, the strain fields obtained with DIC qualitatively reproduce the non-homogeneous deformation pattern due to the presence of stiff fibers dispersed in a compliant matrix and provide accurate results of the average composite strain. However, comparison with finite element simulations revealed that DIC was not able to accurately capture the average strain in each phase.

High fidelity simulations of failure in fiber-reinforced composites

Los ensayos virtuales de materiales compuestos han aparecido como un nuevo concepto dentro de la industria aeroespacial, y disponen de un vasto potencial para reducir los enormes costes de certificación y desarrollo asociados con las tediosas campañas experimentales, que incluyen un gran número de paneles, subcomponentes y componentes. El objetivo de los ensayos virtuales es sustituir algunos ensayos por simulaciones computacionales con alta fidelidad. Esta tesis es una contribución a la aproximación multiescala desarrollada en el Instituto IMDEA Materiales para predecir el comportamiento mecánico de un laminado de material compuesto dadas las propiedades de la lámina y la intercara. La mecánica de daño continuo (CDM) formula el daño intralaminar a nivel constitutivo de material. El modelo de daño intralaminar se combina con elementos cohesivos para representar daño interlaminar. Se desarrolló e implementó un modelo de daño continuo, y se aplicó a configuraciones simples de ensayos en laminados: impactos de baja y alta velocidad, ensayos de tracción, tests a cortadura. El análisis del método y la correlación con experimentos sugiere que los métodos son razonablemente adecuados para los test de impacto, pero insuficientes para el resto de ensayos. Para superar estas limitaciones de CDM, se ha mejorado la aproximación discreta de elementos finitos enriqueciendo la cinemática para incluir discontinuidades embebidas: el método extendido de los elementos finitos (X-FEM). Se adaptó X-FEM para un esquema explícito de integración temporal. El método es capaz de representar cualitativamente los mecanismos de fallo detallados en laminados. Sin embargo, los resultados muestran inconsistencias en la formulación que producen resultados cuantitativos erróneos. Por último, se ha revisado el método tradicional de X-FEM, y se ha desarrollado un nuevo método para superar sus limitaciones: el método cohesivo X-FEM estable. Las propiedades del nuevo método se estudiaron en detalle, y se concluyó que el método es robusto para implementación en códigos explícitos dinámicos escalables, resultando una nueva herramienta útil para la simulación de daño en composites. Virtual testing of composite materials has emerged as a new concept within the aerospace industry. It presents a very large potential to reduce the large certification costs and the long development times associated with the experimental campaigns, involving the testing of a large number of panels, sub-components and components. The aim of virtual testing is to replace some experimental tests by high-fidelity numerical simulations. This work is a contribution to the multiscale approach developed in Institute IMDEA Materials to predict the mechanical behavior of a composite laminate from the properties of the ply and the interply. Continuum Damage Mechanics (CDM) formulates intraply damage at the the material constitutive level. Intraply CDM is combined with cohesive elements to model interply damage. A CDM model was developed, implemented, and applied to simple mechanical tests of laminates: low and high velocity impact, tension of coupons, and shear deformation. The analysis of the results and the comparison with experiments indicated that the performance was reasonably good for the impact tests, but insuficient in the other cases. To overcome the limitations of CDM, the kinematics of the discrete finite element approximation was enhanced to include mesh embedded discontinuities, the eXtended Finite Element Method (X-FEM). The X-FEM was adapted to an explicit time integration scheme and was able to reproduce qualitatively the physical failure mechanisms in a composite laminate. However, the results revealed an inconsistency in the formulation that leads to erroneous quantitative results. Finally, the traditional X-FEM was reviewed, and a new method was developed to overcome its limitations, the stable cohesive X-FEM. The properties of the new method were studied in detail, and it was demonstrated that the new method was robust and can be implemented in a explicit finite element formulation, providing a new tool for damage simulation in composite materials.

Influence of Inter-Fiber Spacing and Interfacial Adhesion on Failure of Multi-Fiber Model Composites: Experiment and Numerical Analysis

The uniaxial tension experiments on glass-fiber-reinforced epoxy matrix composites reveal that the fragmentations of fibers display vertically aligned fracture, clustered fracture, coordinated fracture, and random fracture with the increase of inter-fiber spacing. The finite element analysis indicates that the fragmentations of fibers displaying different phenomena are due to the stress concentration as well as the inherent randomness of fiber defects, which is the dominant factor. The experimental results show that matrices adjacent to the fiber breakpoints all exhibit birefringent-whitening patterns for the composites with different interfacial adhesion strengths. The larger the extent of the interfacial debonding, the less the domain of the birefringent-whitening patterns. The numerical analysis indicates that the orientation of the matrix adjacent to a fiber breakpoint is caused by the interfacial shear stress, resulting in the birefringent-whitening patterns. The area of shear stress concentrations decides on the domain of the birefringent-whitening patterns.

Thermal Degradation of Short Nylon -6 Fiber-Reinforced Styrene Butadiene Rubber Composite

The thermal properties of short Nylon-6 fiber-reinforced Styrene butadiene rubber (SBR) composites were studied by Thermogravimetric Analysis (TGA). The effect of epoxy-based bonding agent on thermal degradation of the gum and the composites was also studied. The thermal stability of the SBR was enhanced in the presence of Nylon-6 fibers and the stability of the composites increased in the presence of bonding agent. The epoxy resin did not significantly change the thermal stability of SBR gum vulcanizate. Results of kinetic studies showed that the degradation of SBR and the short nylon fiber-reinforced composites with and without bonding agents followed first-order kinetics.