This project aimed to integrate, optimize, and validate the fiber orientation and length distribution models previously developed and implemented in the Autodesk® Simulation Moldflow® Insight (ASMI) software package for injection-molded long-carbon-fiber (LCF) thermoplastic composite structures. The project was organized into two phases. Phase 1 demonstrated the ability of the advanced ASMI package to predict fiber orientation and length distributions in LCF/polypropylene (PP) and LCF/polyamide-6, 6 (PA66) plaques within 15% of experimental results. Phase 2 validated the advanced ASMI package by predicting fiber orientation and length distributions within 15% of experimental results for a complex three-dimensional (3D) Toyota automotive part injection-molded from LCF/PP and LCF/PA66 materials. Work under Phase 2 also included estimate of weight savings and cost impacts for a vehicle system using ASMI and structural analyses of the complex part. The present report summarizes the completion of Phases 1 and 2 work activities and accomplishments achieved by the team comprising Pacific Northwest National Laboratory (PNNL); Purdue University (Purdue); Virginia Polytechnic Institute and State University (Virginia Tech); Autodesk, Inc. (Autodesk); PlastiComp, Inc. (PlastiComp); Toyota Research Institute North America (Toyota); Magna Exteriors and Interiors Corp. (Magna); and University of Illinois. Figure 1 illustrates the technical approach adopted in this project that progressed from compounding LCF/PP and LCF/PA66 materials, to process model improvement and implementation, to molding and modeling LCF/PP and LCF/PA66 plaques. The lessons learned from the plaque study and the successful validation of improved process models for fiber orientation and length distributions for these plaques enabled the project to go to Phase 2 to mold, model, and optimize the 3D complex part.

During the last quarter of FY 2015, the following technical progress has been made toward project milestones: 1) PlastiComp used the PlastiComp direct in-line (D-LFT) Pushtrusion system to injection mold 40 30wt% LCF/PP parts with ribs, 40 30wt% LCF/PP parts without ribs, 10 30wt% LCF/PA66 parts with ribs, and 35 30wt% LCF/PA66 parts without ribs. In addition, purge materials from the injection molding nozzle were obtained for fiber length analysis, and molding parameters were sent to PNNL for process modeling. 2) Magna cut samples at four selected locations (named A, B, C and D) from the non-ribbed Magna-molded parts based on a plan discussed with PNNL and the team and shipped these samples to Virginia Tech for fiber orientation and length measurements. 3) Virginia Tech started fiber orientation and length measurements for the samples taken from the complex parts using Virginia Tech's established procedure. 4) PNNL and Autodesk built ASMI models for the complex parts with and without ribs, reviewed process datasheets and performed preliminary analyses of these complex parts using the actual molding parameters received from Magna and PlastiComp to compare predicted to experimental mold filling patterns. 5) Autodesk assisted PNNL in developing the workflow to use Moldflow fiber orientation and length results in ABAQUS® simulations. 6) Autodesk advised the team on the practicality and difficulty of material viscosity characterization from the D-LFT process. 7) PNNL developed a procedure to import fiber orientation and length results from a 3D ASMI analysis to a 3D ABAQUS® model for structural analyses of the complex part for later weight reduction study. 8) In discussion with PNNL and Magna, Toyota developed mechanical test setups and built fixtures for three-point bending and torsion tests of the complex parts. 9) Toyota built a finite element model for the complex parts subjected to torsion loading. 10) PNNL built the 3D ABAQUS® model of the complex ribbed part subjected to 3-point bending. 11) University of Illinois (Prof. C.L. Tucker) advised the team on fiber orientation and fiber length measurement options, modeling issues as well as interpretation of data.

A set of edge-gated and center-gated plaques were injection molded with long carbon fiber-reinforced thermoplastic composites, and the fiber orientation was measured at different locations of the plaques. Autodesk Simulation Moldflow Insight (ASMI) software was used to simulate the injection molding of these plaques and to predict the fiber orientation, using the Anisotropic Rotary Diffusion and the Reduced Strain Closure models. The phenomenological parameters of the orientation models were carefully identified by fitting to the measured orientation data. The fiber orientation predictions show very good agreement with the experimental data.

This quarterly report summarizes the status for the project planning to initiate all the legal and contract documents required for establishing the subcontracts needed and a Cooperative Research and Development Agreement (CRADA) with Autodesk, Inc., Toyota Motor Engineering and Manufacturing North America (Toyota), and Magna Exterior and Interiors Corporation (Magna). During the first quarter (10/1/2012 to 12/31/2012), the statements of work (SOW) for the subcontracts to Purdue University, University of Illinois, and PlastiComp, Inc. were completed. A draft of the CRADA SOW was sent to Autodesk, Toyota, and Magna for technical and legal reviews. PNNL Legal Services contacted project partners’ Legal counterparts for preparing legal documents for the project. A non-disclosure agreement was drafted and sent to all the parties for reviews.

Long fiber filled polymer composite materials have attracted a great attention and usage in recent years. However, the injection and compression molded long fiber composite materials possess complex microstructures that include spatial variations in fiber orientation and length. This paper presents the recent implemented anisotropic rotary diffusion - reduced strain closure (ARD-RSC) model for predicting fiber orientation distribution[1] and a newly developed fiber breakage model[2] for predicting fiber length distribution in injection and compression molding simulation, and Eshelby-Mori-Tanaka model[3,4] with fiber-matrix de-bonding model[5] have been implemented to calculate the long fiber composite property distribution with predicted fiber orientation and fiber length distributions. A validation study on fiber orientation, fiber breakage and mechanical property distributions are given with injection molding process simulation.

In this work, the simulation results from a program developed for the three-dimensional analysis of compression molding of thermoset composites were compared to the experimental results of an actual molded thermoset sample. The program predicts the flow pattern, fiber orientation, fiber length distribution, curing of the resin and mechanical properties. The comparison between simulation and experiment includes short shots and fiber orientation. Actual fiber orientation images were obtained using CT scanning. The comparison shows reasonable agreement between simulation and experiment.

This report describes the work conducted under the CRADA Nr. PNNL/304 between Battelle PNNL and Autodesk whose objective is to validate the new process models developed under the previous CRADA for large injection-molded LFT composite structures. To this end, the ARD-RSC and fiber length attrition models implemented in the 2013 research version of Moldflow was used to simulate the injection molding of 600-mm x 600-mm x 3-mm plaques from 40% glass/polypropylene (Dow Chemical DLGF9411.00) and 40% glass/polyamide 6,6 (DuPont Zytel 75LG40HSL BK031) materials. The injection molding was performed by Injection Technologies, Inc. at Windsor, Ontario (under a subcontract by Oak Ridge National Laboratory, ORNL) using the mold offered by the Automotive Composite Consortium (ACC). Two fill speeds under the same back pressure were used to produce plaques under slow-fill and fast-fill conditions. Also, two gating options were used to achieve the following desired flow patterns: flows in edge-gated plaques and in center-gated plaques. After molding, ORNL performed measurements of fiber orientation and length distributions for process model validations. The structure of this report is as follows. After the Introduction (Section 1), Section 2 provides a summary of the ARD-RSC and fiber length attrition models. A summary of model implementations in the latest research version of Moldflow is given in Section 3. Section 4 provides the key processing conditions and parameters for molding of the ACC plaques. The validations of the ARD-RSC and fiber length attrition models are presented and discussed in Section 5. The conclusions will be drawn in Section 6.

Under the Predictive Engineering effort, PNNL developed linear and nonlinear property prediction models for long-fiber thermoplastics (LFTs). These models were implemented in PNNL’s EMTA and EMTA-NLA codes. While EMTA is a standalone software for the computation of the composites thermoelastic properties, EMTA-NLA presents a series of nonlinear models implemented in ABAQUS® via user subroutines for structural analyses. In all these models, it is assumed that the fibers are linear elastic while the matrix material can exhibit a linear or typical nonlinear behavior depending on the loading prescribed to the composite. The key idea is to model the constitutive behavior of the matrix material and then to use an Eshelby-Mori-Tanaka approach (EMTA) combined with numerical techniques for fiber length and orientation distributions to determine the behavior of the as-formed composite. The basic property prediction models of EMTA and EMTA-NLA have been subject for implementation in the Autodesk® Moldflow® software packages. These models are the elastic stiffness model accounting for fiber length and orientation distributions, the fiber/matrix interface debonding model, and the elastic-plastic models. The PNNL elastic-plastic models for LFTs describes the composite nonlinear stress-strain response up to failure by an elastic-plastic formulation associated with either a micromechanical criterion to predict failure or a continuum damage mechanics formulation coupling damage to plasticity. All the models account for fiber length and orientation distributions as well as fiber/matrix debonding that can occur at any stage of loading. In an effort to transfer the technologies developed under the Predictive Engineering project to the American automotive and plastics industries, PNNL has obtained the approval of the DOE Office of Vehicle Technologies to provide Autodesk, Inc. with the technical support for the implementation of the basic property prediction models of EMTA and EMTA-NLA in the Autodesk® Moldflow® packages. This report summarizes the recent results from Autodesk Simulation Moldlow Insight (ASMI) analyses using the EMTA models and EMTA-NLA/ABAQUS® analyses for further assessment of the EMTA-NLA models to support their implementation in Autodesk Moldflow Structural Alliance (AMSA). PNNL’s technical support to Autodesk, Inc. included (i) providing the theoretical property prediction models as described in published journal articles and reports, (ii) providing explanations of these models and computational procedure, (iii) providing the necessary LFT data for process simulations and property predictions, and (iv) performing ABAQUS/EMTA-NLA analyses to further assess and illustrate the models for selected LFT materials.

The reduced strain closure (RSC) and anisotropic rotary diffusion (ARD) models have been implemented in a research version of Autodesk Moldflow Insight software to improve the accuracy of fiber orientation predictions in simulations of injection-molded parts. The RSC model captures the slow fiber orientation kinetics, which are observed in experiments but over-predicted by the widely used Folgar-Tucker model. The ARD model accounts for fiber-fiber interactions using anisotropic diffusion to rectify the problem that for long-fiber materials, the Folgar-Tucker model does not match all aspects of measured fiber orientation data. Simulating the injection molding process in Autodesk Moldflow Insight software, the fiber orientation is predicted for a number of end-gated ISO plaques and center-gated disks of different thicknesses (from 1.5 to 6 mm) filled at different injection rates, with Midplane meshes as well as 3D meshes, using the RSC model for short-fiber materials and the ARD-RSC model for long-fiber materials. In Midplane models, with an appropriate inlet condition for fiber orientation applied at the gate, the RSC and ARD models show good agreement with measured fiber orientation data. In 3D models with random inlet orientation applied at the gate, the RSC model gives good agreement with one orientation component of the experimental data for short-fiber materials.

Recent experiments suggest that short fibers in concentrated suspensions align more slowly as a function of strain than models based on Jeffery’s equation predict. We develop an objective model that captures the slow orientation kinetics exhibited by short-fiber suspensions. The standard moment-tensor equation of fiber orientation is used to find equations for the change rates of the eigenvalues and eigenvectors of the orientation tensor. As a phenomenological assumption, the growth rates of the eigenvalues are reduced by a constant scalar factor, while the rotation rate expressions for the eigenvectors are unchanged. The eigenvalue/eigenvector equations are then reassembled as a tensor equation. An equivalent kinetic theory is also developed. The new model is tested in a variety of flows, and found to exhibit slower kinetics than the standard model but similar steady-state orientations. The model provides an excellent fit to the shear stress transient in a shear reversal experiment with a 30% glass fiber filled polybutylene terephthalate resin melt, and we show how this experiment can be used to determine the parameters of the model.

The objective of our work is to enable the optimum design of lightweight automotive structural components using injection-molded long fiber thermoplastics (LFTs). To this end, an integrated approach that links process modeling to structural analysis with experimental microstructural characterization and validation is developed. First, process models for LFTs are developed and implemented into processing codes (e.g. ORIENT, Moldflow) to predict the microstructure of the as-formed composite (i.e. fiber length and orientation distributions). In parallel, characterization and testing methods are developed to obtain necessary microstructural data to validate process modeling predictions. Second, the predicted LFT composite microstructure is imported into a structural finite element analysis by ABAQUS to determine the response of the as-formed composite to given boundary conditions. At this stage, constitutive models accounting for the composite microstructure are developed to predict various types of behaviors (i.e. thermoelastic, viscoelastic, elastic-plastic, damage, fatigue, and impact) of LFTs. Experimental methods are also developed to determine material parameters and to validate constitutive models. Such a process-linked-structural modeling approach allows an LFT composite structure to be designed with confidence through numerical simulations. Some recent results of our collaborative research will be illustrated to show the usefulness and applications of this integrated approach.