Search

Fiber reinforced composite materials are a heavily sought after material for next generation vehicles for light-weighting components due to their high specific strength and stiffness. However, when these composite materials are manufactured into laminates, they have relatively weak interlaminar strength and are prone to delamination. This is especially the case when a delamination crack already exists. Quasi-3D (Q3D) braided composites seek to solve this issue by weaving the bias tows into the adjacent (above and below) plies. The plies are physically connected through fiber tows as opposed to being bonded simply by the epoxy, and the composite will achieve a higher interlaminar strength due to fiber failure being required for crack propagation as opposed to simply matrix failure. The [0◦/60◦/-60◦] UD and Q3D carbon composites are investigated in this study for their better in-plane isotropy. Both laminates were manufactured using the vacuum assisted resin transfer molding (VARTM) method with API SC-15 toughened epoxy. The tensile property of the Q3D composite has been found to be competitive with UD composite. Mode I and Mode II interlaminar fracture toughness tests were conducted on UD and Q3D samples. In Mode I experiments, the samples were continuously loaded to full beam split using the double cantilever beam method to obtain the fracture toughness throughout the sample. The Q3D composite shows a large increase in fracture toughness once the crack stabilized and the interlaminar tows become engaged. The Q3D composites also shows a larger energy absorption when the area under the load-extension curve is taken into account. In Mode II testing under end-notch flexural test conditions, the Q3D composite shows

This study investigates the effects of process conditions on the inherent variabilities in fused filament fabrication (FFF) of short carbon-fiber-reinforced Nylon-6 composites, where the sources of uncertainty and their adverse effects on microstructures and Young’s modulus are quantified. Microstructural characteristics such as fiber volume fraction, void volume fraction, and their spatial distributions are first extracted via image-based data analytics, and then their uncertainties are quantified by the analysis of variance. A Monte Carlo sampling method is introduced to enrich the datasets for analyzing uncertainty propagation from micro-level (microstructures) to macro-level (mechanical property). A modified Halpin-Tsai model with the consideration of fiber and void distributions is developed to quantify the propagated uncertainties on Young’s modulus, which are further validated through quasi-static tensile tests. This study examined the process-structure-property relationship of FFF samples and quantified the underlying variations in both micro- and macro-levels.

The material behavior of chopped carbon fiber chip-reinforced Sheet Molding Compound (SMC) composite is investigated under both quasi-static and cyclic fatigue loading conditions in the present study. Quasi-static tensile and compressive strength data is collected, followed by a series of fatigue tests under three different cyclic loading conditions, i.e., Tension-Tension (T-T), Compression-Compression (C-C), and Tension-Compression (T-C). The experimental results show a considerable amount of variation for each loading condition, which is found to be induced by the spatial distribution of fiber orientation within the tested specimens. An analytical approach is developed to correlate the fatigue performance of the SMC material with the local modulus, which is determined by the fiber orientation distribution. Additionally, the similar failure modes are observed under different loading conditions. Based on these findings, a new form of Constant Fatigue Life (CFL) diagram is established in which the impact of fiber orientation tensor to the fatigue behavior of SMC is considered. This developed CFL diagram is examined by comparing the measured fatigue performance with prediction where local fiber orientation distribution is quantified from microscopic images.

Injection molded plastics are widely used for mass production in auto industry. However, it is challenge to produce parts that do not warp. It is helpful to have an accurate way to measure the warpage when investigating and solving warpage issues. In this paper, a new technology is introduced to measure the part warpage quickly and accurately with the help of digital image correlation (DIC) method. The polypropylene plaque were fabricated by injection molding based on An L9 (3ˆ4) Taguchi design of experiments. The warpage of plaque are characterized for all the conditions and at different time after being ejected. Different kinds of warpage results are obtained by DIC method, such as maximum deflection, deflection distribution along desired path, warpage contour on the surface. The plaque warpage increases with the time. The maximum difference can reach as high as 64.9%. The cooling time has the largest influence on the warpage, while the packing pressure has the minor effect. Nevertheless, the warpage difference due to time varying depends on melt temperature, mold temperature, cooling time and packing pressure, in which mold temperature has the highest effects.

Failure analysis is a critical process for the lifetime evaluation of high density polyethylene (HDPE) under cyclic loading since HDPE may encounter different failure modes under different service environments. However, there are only a few investigations on the failure mechanism and failure mode of HDPE in cyclic-loading tests. The objective of this paper is to understand the failure mechanism and propose an approach to predict the failure mode of HDPE under cyclic loading. On the mechanism analysis, our results demonstrated that there could be three failure modes including necking, buckling and fatigue cracking for HDPE under cyclic loading. Both necking and buckling were strongly dependent on the increase of the specimen temperature, which led to the reduction of both yield stress defining necking and tangent modulus governing buckling. On the mode prediction, the yield stress and the tangent modulus could be computed as the functions of specimen temperature which was calculated by FEM heat transfer analysis considering heat generated by inelastic deformation energy and heat transferred with surroundings. Then, necking was predicted by comparing the decreased yield stress with the applied peak stress, and buckling was predicted by comparing the declined buckling stress proportional to tangent modulus with the absolute of the applied valley stress. The proposed methodology was validated as the predicted results were well matched with the experimental data for the tests of HDPE under different cyclic loadings.

The sink mark is one common defect of the injection-molded thermoplastics, which is harmful to the part appearance. A procedure for the sink mark prediction was proposed to build the finite element method (FEM) model in Abaqus. The temperature and pressure distributions were imported as boundary conditions at a different time from the injection molding simulation software Moldflow. The user subroutines UEXPAN and UMAT were coded and implemented in Abaqus to describe the pressure-volume-temperature (PVT) behaviors under actual cooling conditions and the polymer constitutive relationship. The PVT curves of the thermoplastic were obtained by the piston-die method at a cooling rate of 5 °C/min. The crystallizing behaviors were measured by a differential scanning calorimetry and adapted to modify the PVT curves at different cooling rates. The tension-compression tests at different temperatures were performed to obtain the stress-strain curves. The simulation of injection molding was conducted in the software Moldflow. The deflection tests were conducted on the coordinate measuring machine (CMM) to validate the predicted sink mark. It was found that the subroutines of UEXPAN described the volume variations under different cooling rates successfully. Furthermore, the subroutine UMAT predicted stress-strain curves successfully at different temperatures and strain rates. The predicted sink mark in Abaqus agreed well with the measured one. The sink mark depth from Abaqus reached one third of the experimental one and had obvious improvement than that from Moldflow.

In present study, the fatigue damage behavior of chopped carbon fiber chip-reinforced composite has been experimentally and numerically investigated. A large scatter is exhibited in the S-N diagram due to the random distribution of carbon fiber chips, and therefore a new analysis procedure is proposed to link the local microstructure to the fatigue behavior of the bulk material. Interrupted fatigue tests are also performed so that microstructure characterization can be conducted on the tested samples to analyze the crack initiation and propagation. Inspired from the experimental findings, a multi-scale progressive damage fatigue model is proposed to predict the fatigue behavior. The model incorporates a new stochastic chip-packing algorithm for microstructure reconstruction along with continuum damage models into the representative volume element (RVE) model in ABAQUS/Explicit. The simulation results based on the proposed model are in good agreement with the experimental data in terms of cracking modes and predicted life.

This paper presents a computational framework to predict the tensile failure of chopped carbon fiber SMC composites based on the Representative Volume Element (RVE) model. In the framework, a conforming mesh-based reconstruction algorithm is developed to generate mesostructure RVE models of SMC based on the statistical characteristics. Moreover, constitutive models are proposed for the three phases in SMC composites: fiber chip, matrix, and interface. The model predictions are validated by coupon-level experimental tests on SMC samples of various chip orientation distributions. The meso-scale failure mechanism and failure strength predicted by the RVE models are in good agreement with the experimental results. The presented study demonstrates the capability of the present modeling framework in the design of composites materials with complex micro- and meso- structure.

It is important to investigate the constitutive behavior of high density polyethylene (HDPE) under cyclic loading since it is widely used for industrial applications. However, there are only a few preliminary investigations on the constitutive behavior of HDPE under cyclic loading. This work aims at investigating the deformation mechanism, developing the constitutive model especially for the cyclic-loading behavior, and proposing an efficient approach for model calibration. Firstly, the deformation mechanism was successfully explored by conducting a special designed relaxation-unloading test by which the time-independent elastic-plastic behavior was decoupled from the time-dependent overall deformation behavior. Then, a nonlinear elastic-viscoplastic constitutive model was developed based on the parallel rheological framework (PRF) composed of an elastic-plastic network in parallel with multiple nonlinear viscoelastic networks. Finally, an approach was proposed to calibrate the constitutive model by means of relaxation-unloading tests and material parameter optimizations. The agreement between the prediction and the test data successfully demonstrated that the model could be used to precisely simulate the constitutive behavior of HDPE under static, quasi-static and dynamic cyclic loading. The developed constitutive model is beneficial to the long-term durability evaluation of HDPE. The methodology for exploring the deformation mechanism and calibrating the constitutive model may also be applicable to other polymers.

This study proposed a spatiotemporal algorithm to quantitatively characterize the in-situ 3D fracture behavior of carbon-fiber-reinforced polymer (CFRP) composites at microscale . In-situ micro X-ray computed tomography (µXCT) integrated with a tensile stage was applied to capture the 3D fracture evolution of the CFRP composites, where the initiation and propagation of fracture features (e.g., fiber tip-end crack and fiber/matrix debonding) were identified. After the reconstruction of the 3D material microstructure , the proposed spatiotemporal algorithm thereafter extracted the fracture features by employing multiple image processing techniques for quantitative analysis. A similar distribution of the 3D strain obtained from the volumetric digital image correlation demonstrated the feasibility of the developed spatiotemporal algorithm. Moreover, this algorithm provided in-depth and quantitative analysis of fracture features, which provided insights into the microscale failure mechanism and thus shed light on the improvement of failure criteria for CFRP composites with complex microstructures .

This paper presents a discontinuous Galerkin weak form for bond-based peridynamic models to predict the damage of fiber-reinforced composite laminates. To represent the anisotropy of a laminate in a peridynamic model, a lamina is simplified as a transversely isotropic medium under a plane stress condition. The laminated structure is modeled by stacking the surface mesh layers along the thickness direction according to the laminate sequence. To avoid a mesh dependence on either the fiber orientation or the discretization, the spherical harmonic expansion theory is employed to construct a function for the micro-elastic modulus in terms of the bond-fiber angle. The laminate material is decomposed into an isotropic matrix material part and a transversely isotropic material part. The bond stiffness can be evaluated using the engineering material constants, based on the equivalence between the elastic energy density in the peridynamic theory and the elastic energy density in the classic continuum mechanics theory. Benchmark tests are conducted to verify the proposed model. Numerical results illustrate that the convergence of simulations with different horizon sizes and meshes can be achieved. In terms of damage analysis, the proposed model can capture the dynamic process of the complex coupling of the inner-layer and delamination damage modes.

The mechanical properties and damage behaviors of carbon/epoxy woven fabric composite under in-plane tension and compression are studied at the meso-scale level through experiment and simulation. An efficient representative volume element (RVE) modeling method with consistent mesh, high yarn volume fraction and realistic geometry is proposed. The material constitutive laws with plasticity, tension-compression asymmetry and damage evolution are established for the three components - yarn, matrix and interface, respectively. Significantly different mechanical properties and damage evolutions are observed depending on loading conditions and initial geometry characteristics. It shows a non-linear stress-strain curve with clear transition region and intensive damage in tension, while a quasi-linear behavior up to facture is observed in compression with little damage prior to final fracture. Moreover, compared to the constant Poisson's ratio with straining in compression, a dramatic increase in Poisson's ratio appears in tension. Simulation shows damage mechanisms including transverse damage, matrix damage and delamination, which all play critical roles in the property evolution. In particular, the rapid damage accumulation after elastic deformation destroys the strong bonds and causes the easy deformation of transverse yarns which results in the transition region and large Poisson's ratio in tension. All the mechanical behaviors and damage evolutions are well captured and explained with the current RVE model.

Woven fiber composites have been increasingly employed as light-weight materials in aerospace, construction, and transportation industries due to their superior properties. These materials possess a hierarchical structure that necessitates the use of multiscale simulations in their modeling. To account for the inherent uncertainty in materials, such simulations must be integrated with statistical uncertainty quantification (UQ) and propagation (UP) methods. However, limited advancement has been made in this regard due to the significant computational costs and complexities in modeling spatially correlated structural variations coupled at different scales. In this work, a non-intrusive approach is proposed for multiscale UQ and UP to address these limitations. We introduce the top-down sampling method that allows to model non-stationary and continuous (but not differentiable) spatial variations of uncertainty sources by creating nested random fields (RFs) where the hyperparameters of an ensemble of RFs is characterized by yet another RF. We employ multi-response Gaussian RFs in top-down sampling and leverage statistical techniques (such as metamodeling and dimensionality reduction) to address the considerable computational costs of multiscale simulations. We apply our approach to quantify the uncertainty in a cured woven composite due to spatial variations of yarn angle, fiber volume fraction, and fiber misalignment angle. Our results indicate that, even in linear analysis, the effect of uncertainty sources on the material’s response could be significant.

Sheet metals such as advanced high strength steels often exhibit varying degrees of anisotropy due to the existence of preferred orientation of their grain structures from the rolling process. While extensive research has been devoted to characterize and model anisotropic plasticity behavior over the past decades, only recently people began to tackle their anisotropic fracture behavior. The focus of these studies is almost always on the in-plane anisotropy, and ignores out-of-plane fracture behavior. The reasons are two-folds: (1). Sheet metal deformation during a stamping or crash application is mostly thought as in a plane stress state and modeled as such, so there is no need to consider out-of-plane fracture behavior; and (2). It is very challenging to develop experimental techniques to reliably measure out-of-plane fracture strengths for sheet metals with a thickness only about 1 mm. In this paper, we theorize, based on strong empirical evidence, that the fracture strength of a thin sheet metal in out-of-plane shear is inherently weaker than that under in-plane shear, analogous to the interlaminar properties of composite laminates. This weakness is responsible for the so-called “shear fracture” behavior, often observed during the stamping operations of advanced high strength steels when the sheet metal flows around a tight radius. A new double-notched shear specimen with shear plane perpendicular to the thickness direction is carefully designed and tested to accurately measure the out-of-plane shear fracture strength. Simulations are executed to validate the rationality of the out-of-plane shear specimen design. Test results of a DP980 sheet show that the out-of-plane shear fracture strength is indeed about 15% lower than that measured under the in-plane shear condition. A new anisotropic ductile fracture model based on linear transformation of stresses is proposed to account for the out-of-plane fracture strength for AHSS sheets. The model is calibrated with experimental data, and the fracture surface is compared to the isotropic Mohr-Coulomb (M-C) fracture model to illustrate its effectiveness.

Predicting the mechanical behavior of the chopped carbon fiber Sheet Molding Compound (SMC) due to spatial variations in local material properties is critical for the structural performance analysis but is computationally challenging. Such spatial variations are induced by the material flow in the compression molding process. In this work, a new multiscale SMC modeling framework and the associated computational techniques are developed to provide accurate and efficient predictions of SMC mechanical performance. The proposed multiscale modeling framework contains three modules. First, a stochastic algorithm for 3D chip-packing reconstruction is developed to efficiently generate the SMC mesoscale Representative Volume Element (RVE) model for Finite Element Analysis (FEA). A new fiber orientation tensor recovery function is embedded in the reconstruction algorithm to match reconstructions with the target characteristics of fiber orientation distribution. Second, a metamodeling module is established to improve the computational efficiency by creating the surrogates of mesoscale analyses. Third, the macroscale behaviors are predicted by an efficient multiscale model, in which the spatially varying material properties are obtained based on the local fiber orientation tensors. Our approach is validated through experiments at both meso- and macro-scales, such as tensile tests assisted by Digital Image Correlation (DIC) and mesostructure imaging.

In current work, the crashworthiness of woven carbon fiber reinforced plastic (CFRP) hat section component is investigated under two types of loading experimentally and numerically, i.e., axial crushing and dynamic three-point bending. The crushing characteristics of woven CFRP hat section sample are evaluated at various impact velocities and layup orientations which exhibit different effects under two types of loading. With impact velocity/energy increasing, the bending behaviors illustrate different evolutions with distinct stages which are affected by the damage development. Layup orientation shows a little effect on the bending behaviors, while it exerts remarkable effects on the axial crushing which leads to different crushing modes and thus different energy absorption capabilities. The finite element model with multiple layers of thick shell element and cohesive element is shown to capture the crushing behaviors reasonably well in both axial crushing and dynamic bending. The corresponding damage evolutions are also analyzed with simulations to explain the crushing behaviors.

Woven composites are considered as a promising material choice for lightweight applications. An improved non-orthogonal material model that can decouple the strong tension and weak shear behavior of the woven composite under large shear deformation is proposed for simulating the preforming of woven composites. The tension, shear and compression moduli in the model are calibrated using the tension, bias-extension and bending experiments, respectively. The interaction between the composite layers is characterized by a sliding test. The newly developed material model is implemented in the commercial finite element software LS-DYNA® and validated by a double dome study.

Carbon fiber reinforced plastics (CFRPs) are attracting growing attention in industry because of their enhanced properties. Preforming of thermoset carbon fiber prepregs is one of the most common production techniques of CFRPs. To simulate preforming, several computational methods have been developed. Most of these methods, however, obtain the material properties directly from experiments such as uniaxial tension and bias-extension where the coupling effect between tension and shear is not considered. Neglecting this coupling effect deteriorates the prediction accuracy of simulations. To address this issue, we develop a Bayesian model calibration and material characterization approach in a multiscale finite element preforming simulation framework that utilizes mesoscopic representative volume element (RVE) to account for the tension-shear coupling. A new geometric modeling technique is first proposed to generate the RVE corresponding to the close-packed uncured prepreg. This RVE model is then calibrated with a modular Bayesian approach to estimate the yarn properties, test its potential biases against the experiments, and fit a stress emulator. The predictive capability of this multiscale approach is further demonstrated by employing the stress emulator in the macroscale preforming simulation which shows that this approach can provide accurate predictions.

We present a data-driven nonlinear uncertainty quantification and propagation framework to study the microstructure-induced stochastic performance of unidirectional (UD) carbon fiber reinforced polymer (CFRP) composites. The proposed approach integrates (1) microscopic image characterization, (2) stochastic microstructure reconstruction, and (3) efficient multiscale finite element simulations enabled by self-consistent clustering (SCA) analysis. To model the complex microstructural variability, the proposed UQ methods take the non-Gaussian uncertainty sources into account through a distribution-free sampling approach leveraging nonparametric and asymptotic statistical tools. A hierarchical conditional sampling strategy enables the simultaneous sampling of multiple sources of uncertainties. Our approach provides insights into the impact of microstructural variabilities, which are shown to have an increasing impact on the nonlinear responses of UD CFRP parts under progressive compression loading and ultimately on the failure rate over time. We discover that before CFRP parts start to fail, a characteristic time period emerges with distinctive uncertainty distributions specific to the microstructure variability and the probability of failure. Identifying the failure time period is crucial to the reliability prediction, which is an essential component of CFRP design.

The mechanical behaviors and damage evolutions of carbon/epoxy woven fabric composites with three different geometries, i.e., one plain weave and two twill weave patterns with different areal densities, are studied under tensile loading. The effects of weave patterns on mechanical properties are investigated by monotonic and cyclic tension tests. Remarkable variations in stress–strain curve, Poisson’s ratio, residual strain and strain map exist in the three composites. Crimp ratio is found to be a critical factor to govern the mechanical properties. With smaller crimp ratio, a quasi-linear stress–strain curve with higher elastic modulus and strength is observed. The stress–strain curves of composites with higher crimp ratio contain transition stages with significant tangent modulus degradation. Elastic modulus, strength and damage initiation are all correlated with the crimp ratio linearly regardless of the fabric pattern. Dramatic nonlinear evolution in Poisson’s ratio occurs in the composite with higher crimp ratio. Cyclic tension results indicate that the residual strain is a more appropriate damage indicator than the unloading elastic modulus. Microstructure examination shows that damage developments are essentially related to the fabric geometry, and result in various mechanical behaviors. This study provides important insights into the geometry-deformation mechanism-mechanical property relationship of the woven composites.

To develop a composite with high in-plane properties and sufficient delamination resistance, a quasi-three-dimensional (Q3D) woven composite has been designed. In Q3D, some yarns in each layer are woven into those in the adjacent layers and form a 3D interlocking network. This study examines the performance of a carbon fiber/SC-15 epoxy triaxial Q3D [ (0°/60°/−60°)4 ] composite. For comparison, a [0°/60°/ − 60°]4 composite laminate with unidirectional (UD) plies and a [ 0°/60°/−60° ]4 laminate with two-dimensional woven (2DW) layers were fabricated using the same carbon fiber and epoxy system. The in-plane tensile, interlaminar fracture, and low-velocity impact tests were performed for all three types of composites. The results show that the Q3D composite has the in-plane tensile modulus and strength comparable to that of UD and 2DW composites, but offers a higher resistance to interlaminar crack propagation and impact damage. This work provides insights into the influences of fiber structures on composite properties and performances.

The mechanical properties and damage evolutions of carbon/epoxy woven fabric composites with three different fabric architectures, including one plain weave and two twill weave patterns, are experimentally investigated under multiaxial stress states. In particular, the effects of weave patterns are investigated by monotonic and cyclic off-axis tension tests. Both elastic modulus and strength degrade remarkably with increasing off-axis loading angle, while Poisson's ratio is much higher than that measured from on-axis tests and increases with loading strain gradually. Different fabric architectures show limited effects on the modulus and strength under multiaxial stress states, and they are well predicted by transformation equation and Tsai-Wu failure criteria, respectively. However, significantly different failure behaviors are observed in three fabric composites, and microstructure observation shows that fabric architecture affects the stress concentration and the damage development. Smaller crimp ratio and compacted structure postpone the damage development but result in more abrupt failure under multiaxial stress states.

In this work, an RVE model with embedded solid elements (RVE-ESE) is developed to predict the effective elastic constants of discontinuous fiber reinforced composites (DFRC). Unlike the traditional finite element (FE) modeling, the matrix and fibers can be meshed separately and independently due to the implementation of the embedded element technique, which enables the ease of RVE meshing without strictly satisfying the mesh conformity at the interfaces. Since the structured meshes can be used for the matrix, it becomes quite easy to impose the periodic boundary conditions on the RVE. In addition, a modified random sequential adsorption (RSA) algorithm based on the Boolean operation is proposed to generate the RVE by using the Python script in the ABAQUS/CAE. This algorithm can determine whether a new generated fiber overlaps with the pre-existing ones without looping calculation. Some important issues related to the modeling accuracy, such as the RVE size, the fiber separation effect and the mesh sensitivity of both matrix and fibers, are also investigated. It is found that the proposed model agrees very well with the traditional FE direct meshing approach when the element sizes of the matrix and fibers in the RVE-ESE model are sufficiently fine. Compared with commercial software Digimat, the current model yields relatively higher prediction on the elastic moduli of the composites with randomly distributed fibers. However, for the composites with aligned fibers, Digimat shows lower result for the transverse Young's modulus while a bit higher prediction for the longitudinal one.

The dynamic crushing characteristics of unidirectional carbon fiber reinforced plastic composites under two loading types, dynamic three-point bending and axial crushing, are investigated by experiment and finite element simulation. Hat section samples with two different layup orientations are tested at various impact velocities to investigate its effects. The experiment results show that delamination plays a critical role in dynamic bending deformation. Different layup orientations lead to similar crushing bending behaviors, but remarkable variations in delamination pattern evolutions. In contrast, layup orientation exerts more significant effects on failure modes and energy absorption in axial crushing. A finite element model with multiple layers of thick shell element and cohesive element is established, and simulations are performed for both dynamic three-point bending and axial crushing. The model can capture the crushing behaviors in the two loading types reasonably well. The damage evolution and energy absorption through various mechanisms at different loading conditions are discussed.

This paper addresses the challenge of reconstructing nonuniformly orientated fiber-reinforced polymer composites (FRPs) with three-dimensional (3D) geometric complexity, especially for fibers with curvatures, and proposes a framework using micro X-ray computed tomography (μXCT) images to quantify the fiber characteristics in 3D space for elastic modulus prediction. The FRP microstructure is first obtained from the μXCT images. Then, the fiber centerlines are efficiently extracted with the proposed fiber reconstruction algorithm, i.e., iterative template matching, and the 3D coordinates of the fiber centerlines are adopted for quantitative characterization of the fiber morphology. Finally, Young's modulus is predicted using the Halpin-Tsai model and laminate analogy approach, and the fiber configuration averaging method with the consideration of the fiber morphology. The new framework is demonstrated on both injection-molded short and long carbon fiber-reinforced polymer composites, whose fiber morphology and predicted mechanical properties are validated through previous pyrolysis and quasi-static tensile tests, respectively.

This paper presents a predictive multiscale modeling scheme for Unidirectional (UD) Carbon Fiber Reinforced Polymers (CFRP). A bottom-up modeling procedure is discussed for predicting the performance of UD structures. UD material responses are computed from high fidelity Representative Volume Elements (RVEs). A data-driven Reduced Order Modeling approach compresses RVEs into Reduced Order Models so material responses can be computed in a concurrent fashion along with the structural level simulation. The approach presented in this paper is validated against experimental data and is believed to provide design guidance for future fiber reinforced polymers development.

Carbon fiber reinforced composites have received growing attention because of their superior performance and high potential for lightweight systems. An economic method to manufacture the parts made of these composites is a sequence of forming followed by a compression molding. The first step in this sequence is called preforming that forms the prepreg, which is the fabric impregnated with the uncured resin, to the product geometry, while the molding process cures the resin. Slip between different prepreg layers is observed in the preforming step, and it is believed to have a non-negligible impact on the resulting geometry. This paper reports a method to characterize the interaction between different prepreg layers, which should be valuable for future predictive modeling and design optimization. An experimental device was built to evaluate the interactions with respect to various industrial production conditions. The experimental results were analyzed for an in-depth understanding about how temperature, relative sliding speed, and fiber orientation affect the tangential interaction between two prepreg layers. Moreover, a hydro-lubricant model was introduced to study the relative motion mechanism of this fabric-resin-fabric system, and the results agreed well with the experiment data. The interaction factors obtained from this research will be implemented in a preforming process finite element simulation model.

In most current works related to the structure analysis of cured woven composites, the influence of yarn angle variation on the local material properties is neglected. In this paper, its influence was studied, and an integrated preforming-performance simulation model, accounting for the influence of yarn angle variation, was proposed. The multi-scale modeling approach was adopted to predict the material properties of cured woven composites with different yarn angles. The results were compared with experiments to validate its effectiveness. In the integrated preforming-performance simulation model, non-orthogonal constitutive law was used in the preforming simulation to compute yarn orientations and yarn angle distribution. To demonstrate the capabilities of the integrated simulation model, simulations were conducted and compared with experiments. The results show that the proposed simulation model has more accurate prediction than the simulation model without considering the influence of yarn angle, and the impact of yarn angle shouldn’t be neglected.

Compared with the conventional steels, "shear fracture" is one of the main issues for advanced high strength steels (AHSS). Due to rolling, anisotropy is an intrinsic property for sheet metals. Not only the plastic responses of sheet metals but also the fracture strengths are orientation dependent. In the small radius forming process, for example, the stretch-bending deformation of sheet metals under small radius condition, the normal stress cannot be neglected. Three-dimensional loading condition constructs complex shear stress states of sheet metals especially the out-of-plane shear stress. The out-of-plane performance must be considered in order to better understand the "shear fracture" phenomenon of AHSS. Compared to in-plane shear test, the out-of-plane shear test is more difficult to carry out due to the severe restriction of the dimensions in the thickness direction. In this paper, a new specimen is presented for out-of-plane shear test. Failure of the specimen occurs in shear between two centrally located notches machined halfway through its thickness from opposing sides. Meanwhile, the finite element (FE) model and possible failure modes of this specimen are investigated in detail. At last, brief experimental results between out-of-plane shear fracture strength and the in-plane shear fracture strength are compared for DP980 sheets.