Engineering topics

Bednarcyk, Brett A.

Publications and source records attributed to Bednarcyk, Brett A..

At least 19 records

NASA NTRS2020

Solution of the Nonlinear High-Fidelity Generalized Method of Cells Micromechanics Relations via Order-Reduction Techniques

The High-Fidelity Generalized Method of Cells (HFGMC) is one technique, distinct from traditional finite element approaches, for accurately simulating nonlinear composite material behavior. In this work, the HFGMC global system of equations for doubly-periodic repeating unit cells with nonlinear constituents has been reduced in size through the novel application of a Petrov-Galerkin Proper Orthogonal Decomposition order-reduction scheme in order to improve its computational efficiency. Order-reduced models of an E-glass/Nylon 12 composite led to a 4.8 to 6.3x speedup in the equation assembly/solution runtime while maintaining model accuracy. This corresponded to a 21 to 38 percent reduction in total runtime. The significant difference in assembly/solution and total runtimes was attributed to the evaluation of integration point inelastic field quantities; this step was identical between the unreduced and order-reduced models. Nonetheless, order-reduced techniques offer the potential to significantly improve the computational efficiency of multiscale calculations.

Ricks, Trenton M.↗

NASA NTRS2020

Stochastic Simulation of Mudcrack Damage Formation in an Environmental Barrier Coating

The FEAMAC/CARES program, which integrates finite element analysis (FEA) with the MAC/GMC (Micromechanics Analysis Code with Generalized Method of Cells) and the CARES/Life (Ceramics Analysis and Reliability Evaluation of Structures / Life Prediction) programs, was used to simulate the formation of mudcracks during the cooling of a multilayered environmental barrier coating (EBC) deposited on a silicon carbide substrate. FEAMAC/CARES combines the MAC/GMC multiscale micromechanics analysis capability (primarily developed for composite materials) with the CARES/Life probabilistic multiaxial failure criteria (developed for brittle ceramic materials) and Abaqus (Dassault Systèmes) FEA. In this report, elastic modulus reduction of randomly damaged finite elements was used to represent discrete cracking events. The use of many small-sized low-aspect-ratio elements enabled the formation of crack boundaries, leading to development of mudcrack-patterned damage. Finite element models of a disk-shaped three-dimensional specimen and a twodimensional model of a through-the-thickness cross section subjected to progressive cooling from 1,300 °C to an ambient temperature of 23 °C were made. Mudcrack damage in the coating resulted from the buildup of residual tensile stresses between the individual material constituents because of thermal expansion mismatches between coating layers and the substrate. A two-parameter Weibull distribution characterized the coating layer stochastic strength response and allowed the effect of the Weibull modulus on the formation of damage and crack segmentation lengths to be studied. The spontaneous initiation of cracking and crack coalescence resulted in progressively smaller mudcrack cells as cooling progressed, consistent with a fractal-behaved fracture pattern. Other failure modes such as delamination, and possibly spallation, could also be reproduced. The physical basis assumed and the heuristic approach employed, which involves a simple stochastic cellular automaton methodology to approximate the crack growth process, are described. The results ultimately show that a selforganizing mudcrack formation can derive from a Weibull distribution that is used to describe the stochastic strength response of the bulk brittle ceramic material layers of an EBC.

Nemeth, Noel N.↗

NASA NTRS2020

Computationally-efficient Structural Models for Analysis of Woven Composites

The paper presents a novel approach to model woven composite using the computationally efficient one-dimensional models. The framework is built within the scheme of the Carrera Unified Formulation (CUF), a generalized hierarchical formulation that generates variable kinematic structural theories. Various components of the woven composite unit cell are modeled using a combination of straight and curved one-dimensional CUF models. By employing a component-wise approach, a modeling technique within CUF, the complex geometry of the woven composite components is modeled precisely. The ability of CUF models to accurately resolve stress and strain fields are exploited to capture complex deformation within a woven composite unit cell. Numerical results include analyses of a non-crimped textile composite, a curved tow under tension, and a dry woven textile unit cell.

Kaleel, Ibrahim↗

NASA NTRS2019

Comparison of Design Tools for Stress Analysis of Adhesively Bonded Joints

Analytical or semi-analytical models for stress analysis have long been a part of initial adhesive bond design and sizing. Even with the rise of general finite element software and methods, these design models have still remained a preferred method for fast and simple joint analysis. While these methods can yield fairly representative results, the models are constructed on the foundation of geometrical and material simplifications or assumptions that allow closed form or semi-closed form solutions. This study outlines the major differences in basic assumptions for three common design software packages under use at NASA, and shows the ramifications these assumptions in a few exemplar bonded joints.

Composites↗

NASA NTRS2019

A Comparison of Different Modeling Strategies for Predicting Effective Properties of 3D Woven Composites

Three-dimensional (3D) woven composites are an attractive means of achieving superior mechanical performance in aerospace structures. Limited analysis capability currently exists to predict both effective elastic and strength properties for these complex composites. In this study, a comparison of three modeling strategies was performed to assess the ability of the different methods to predict the effective elastic properties of four distinct 3D orthogonal woven composites. Two finite element techniques (in-plane and triply-periodic boundary conditions) and one method of cells technique, the Multiscale Generalized Method of Cells, were considered.

Ricks, Trenton M.↗

NASA NTRS2019

Integration of MAC/GMC into CalculiX, an Open Source Finite Element Code

An analysis framework is presented that makes available multiscale analysis of composite structures using the open-source FEA solver package CalculiX CrunchiX (CCX). At the center of this framework is the coupling of the Finite Element Analysis - Micromechanics Analysis Code (FEAMAC) library from NASA's Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC) to allow micromechanics analysis. The results show that the proposed coupling can be used with appropriate care for multiscale FEA simulations of composite materials. The largest error reported in this validation was in a four-point bend test specimen with an error of less than 1% difference in the maximum deflection of the beam.

Yapor Genao, Francisco A.↗

NASA NTRS2018

Modeling-Driven Damage Tolerant Design of Graphene Nanoplatelet/Carbon Fiber/Epoxy Hybrid Composite Panels for Full-Scale Aerospace Structures

The objective of this study is to design a new nano graphenecarbon fiberpolymer hybrid composite that can be used for the NASA SLS Composite Exploration Upper Stage (CEUS) forward skirt structure. The new material will improve the resistance to open-hole compression failure of the structure relative to traditional polymer fiber composites. The material is designed rapidly and with little cost using the Integrated Computational Materials Engineering (ICME) approach. Multiscale modeling and experiments are used to synergistically optimize the material design to yield improved properties and performance by controlling key processing parameters for manufacturing nano-enhanced materials. Specifically, the nanocomposite panel showed a 22 reduction in mass relative to the traditional composite panel, while both designs are equal in terms of ease of manufacture. This potential mass savings corresponds to an estimated 45 savings in materials and manufacturing costs. The multiscale ICME workflow developed for this project can be readily applied to the development of nano-enhanced composite materials and large aerospace structures. In addition, all key aspects of ICME were employed to complete this project including multiscale modeling, experimental characterization and visualization, data management, visualization, error and uncertainty quantification, and education. The results presented herein indicate a dramatic level of success, as well as the power and potential of ICME approach and multiscale modeling for composite materials.

computational mechanics↗

NASA NTRS2018

Analysis of the Influence of Layer Shifting on the Elastic Response and Damage Nucleation and Growth in Woven Composite Laminates

The influence of relative layer shifting on the elastic and damage response of plain weave composite laminates is analyzed using a continuum damage mechanics approach in combination with the finite element method. First, the homogenized properties of the woven composite as a function of the number of layers and of layer shifting are presented. Next, the damage development in various shifting configurations is studied using different damage constitutive models for the matrix and the fiber bundles. It is shown that the impact of layer shifting on both the elastic response and the nonlinear damage response is significant. Most notably, the model captures changes in the damage mechanisms within the woven composite that occur due to layer shifting, resulting in stiffer, more brittle behavior, which has been shown experimentally in the literature. Model results in the linear and nonlinear regimes are shown to be consistent with both an independent analytical model and reported experiments.

layer shifting↗

NASA NTRS2018

Modeling of Time-Dependent Strength Degradation of SiC/SiC Ceramic Matrix Composites via MAC/GMC Computer Code

Silicon carbide fiber reinforced silicon carbide (SiC/SiC) ceramic matrix composites (CMCs) display time-dependent strength degradation at intermediate temperatures (600 to 900 degrees Centigrade). This is generally believed to be an oxidation induced phenomenon. The understanding of the effect of temperature with environment (oxidation) is key towards development of SiC/SiC CMCs with a reliable load carrying capacity. Various theories have been proposed to explain the strength degradation. One suggests that the boron nitride (BN) coating deposited on the fibers oxidizes causing fusion of fibers. Another theory proposes that the SiC fibers are oxidized forming a silica scale leading to premature fiber failure. A more recent theory suggests that SiC fiber strength is intrinsically time-dependent due to slow crack growth of flaws in the fibers. An empirical model, termed as a "fiber classic model," which is based on a standard slow crack growth type power-law, has been implemented within NASA's micromechanics-based MAC/GMC computer code as a user routine. Model parameters for this "classic model" were calibrated from stress-rupture data of Hi-Nicalon (TM) monofilaments using the maximum likelihood estimation (MLE) technique. This new capability in the MAC/GMC computer code was then used to predict the stress-rupture behavior of Hi-Nicalon (TM) tows as well as 2-D SiC/SiC composites reinforced with Hi-Nicalon (TM) fibers. Results demonstrate that the MAC/GMC with this new capability successfully predicts the time-to-failure vs. applied stress within the intermediate temperature range at various scales as well as laminated composites in an oxidizing environment.

Environmental degradation↗

NASA NTRS2018

Multiscale Modeling of PEEK Using Reactive Molecular Dynamics Modeling and Micromechanics

Polyether ether ketone (PEEK) is a high-performance, semi-crystalline thermoplastic that is used in a wide range of engineering applications, including some structural components of aircraft. The design of new PEEK-based materials requires a precise understanding of the multiscale structure and behavior of semi-crystalline PEEK. Molecular Dynamics (MD) modeling can efficiently predict bulk-level properties of single phase polymers, and micromechanics can be used to homogenize those phases based on the overall polymer microstructure. In this study, MD modeling was used to predict the mechanical properties of the amorphous and crystalline phases of PEEK. The hierarchical microstructure of PEEK, which combines the aforementioned phases, was modeled using a multiscale modeling approach facilitated by NASA's MSGMC. The bulk mechanical properties of semi-crystalline PEEK predicted using MD modeling and MSGMC agree well with vendor data, thus validating the multiscale modeling approach.

Pisani, William A.↗

NASA NTRS2018

Development and Application of a Tool for Optimizing Composite Matrix Viscoplastic Material Parameters

This document describes a recently developed analysis tool that enhances the resident capabilities of the Micromechanics Analysis Code with the Generalized Method of Cells (MAC/GMC) and its application. MAC/GMC is a composite material and laminate analysis software package developed at NASA Glenn Research Center. The primary focus of the current effort is to provide a graphical user interface (GUI) capability that helps users optimize highly nonlinear viscoplastic constitutive law parameters by fitting experimentally observed/measured stress-strain responses under various thermo-mechanical conditions for braided composites. The tool has been developed utilizing the MATrix LABoratory (MATLAB) (The Mathworks, Inc., Natick, MA) programming language. Illustrative examples shown are for a specific braided composite system wherein the matrix viscoplastic behavior is represented by a constitutive law described by seven parameters. The tool is general enough to fit any number of experimentally observed stress-strain responses of the material. The number of parameters to be optimized, as well as the importance given to each stress-strain response, are user choice. Three different optimization algorithms are included: (1) Optimization based on gradient method, (2) Genetic algorithm (GA) based optimization and (3) Particle Swarm Optimization (PSO). The user can mix and match the three algorithms. For example, one can start optimization with either 2 or 3 and then use the optimized solution to further fine tune with approach 1. The secondary focus of this paper is to demonstrate the application of this tool to optimize/calibrate parameters for a nonlinear viscoplastic matrix to predict stress-strain curves (for constituent and composite levels) at different rates, temperatures and/or loading conditions utilizing the Generalized Method of Cells. After preliminary validation of the tool through comparison with experimental results, a detailed virtual parametric study is presented wherein the combined effects of temperature and loading rate on the predicted response of a braided composite is investigated.

thermo-mechanical; MATLAB; GUI-Based Too↗

NASA NTRS2017

Order-Reduced Solution of the Nonlinear High-Fidelity Generalized Method of Cells Micromechanics Relations

The High-Fidelity Generalized Method of Cells (HFGMC) is one technique for accurately simulating nonlinear composite material behavior. The HFGMC uses a higher-order approximation for the subcell displacement field that allows for a more accurate determination of the subcell stressstrain fields at the cost of some computational efficiency. In order to reduce computational costs associated with the solution of the ensuing system of simultaneous equations, the HFGMC global system of equations for doubly-periodic repeating unit cells with nonlinear constituents was reduced in size through the use of a Petrov-Galerkin-based Proper Orthogonal Decomposition order-reduction scheme. A number of cases were presented that address the computational feasibility of using order-reduction techniques to solve solid mechanics problems involving complex microstructures.

order-reduction↗

NASA NTRS2017

Constitutive Modeling of the Facesheet to Core Interface in Honeycomb Sandwich Panels Subject to Mode I Delamination

A new cohesive zone traction-separation law, which includes the effects of fiber bridging, has been developed, implemented with a finite element (FE) model, and applied to simulate the delamination between the facesheet and core of a composite honeycomb sandwich panel. The proposed traction-separation law includes a standard initial cohesive component, which accounts for the initial interfacial stiffness and energy release rate, along with a new component to account for the fiber bridging contribution to the delamination process. Single cantilever beam tests on aluminum honeycomb sandwich panels with carbon fiber reinforced polymer facesheets were used to characterize and evaluate the new formulation and its finite element implementation. These tests, designed to evaluate the mode I toughness of the facesheet to core interface, exhibited significant fiber bridging and large crack process zones, giving rise to a concave downward concave upward pre-peak shape in the load-displacement curve. Unlike standard cohesive formulations, the proposed formulation captures this observed shape, and its results have been shown to be in excellent quantitative agreement with experimental load-displacement and apparent critical energy release rate results, representative of a payload fairing structure, as well as local strain fields measured with digital image correlation.

composites↗

NASA NTRS2017

Comparison of Multiscale Method of Cells-Based Models for Predicting Elastic Properties of Filament Wound C/C-SiC

Three different multiscale models, based on the method of cells (generalized and high fidelity) micromechanics models were developed and used to predict the elastic properties of C/C-SiC composites. In particular, the following multiscale modeling strategies were employed: Concurrent multiscale modeling of all phases using the generalized method of cells, synergistic (two-way coupling in space) multiscale modeling with the generalized method of cells, and hierarchical (one-way coupling in space) multiscale modeling with the high fidelity generalized method of cells. The three models are validated against data from a hierarchical multiscale finite element model in the literature for a repeating unit cell of C/C-SiC. Furthermore, the multiscale models are used in conjunction with classical lamination theory to predict the stiffness of C/C-SiC plates manufactured via a wet filament winding and liquid silicon infiltration process recently developed by the German Aerospace Institute.

multiscale modeling↗

NASA NTRS2017

A Compilation of MATLAB Scripts and Functions for MACGMC Analyses

The primary aim of the current effort is to provide scripts that automate many of the repetitive pre- and post-processing tasks associated with composite materials analyses using the Micromechanics Analysis Code with the Generalized Method of Cells. This document consists of a compilation of hundreds of scripts that were developed in MATLAB (The Mathworks, Inc., Natick, MA) programming language and consolidated into 16 MATLAB functions. (MACGMC). MACGMC is a composite material and laminate analysis software code developed at NASA Glenn Research Center. The software package has been built around the generalized method of cells (GMC) family of micromechanics theories. The computer code is developed with a user-friendly framework, along with a library of local inelastic, damage, and failure models. Further, application of simulated thermo-mechanical loading, generation of output results, and selection of architectures to represent the composite material have been automated to increase the user friendliness, as well as to make it more robust in terms of input preparation and code execution. Finally, classical lamination theory has been implemented within the software, wherein GMC is used to model the composite material response of each ply. Thus, the full range of GMC composite material capabilities is available for analysis of arbitrary laminate configurations as well. The pre-processing tasks include generation of a multitude of different repeating unit cells (RUCs) for CMCs and PMCs, visualization of RUCs from MACGMC input and output files and generation of the RUC section of a MACGMC input file. The post-processing tasks include visualization of the predicted composite response, such as local stress and strain contours, damage initiation and progression, stress-strain behavior, and fatigue response. In addition to the above, several miscellaneous scripts have been developed that can be used to perform repeated Monte-Carlo simulations to enable probabilistic simulations with minimal manual intervention. This document is formatted to provide MATLAB source files and descriptions of how to utilize them. It is assumed that the user has a basic understanding of how MATLAB scripts work and some MATLAB programming experience.

Materials Engineering↗

NASA NTRS2017

Micromechanics of Composite Materials Governed by Vector Constitutive Laws

The high-fidelity generalized method of cells micromechanics theory has been extended for the prediction of the effective property tensor and the corresponding local field distributions for composites whose constituents are governed by vector constitutive laws. As shown, the shear analogy, which can predict effective transverse properties, is not valid in the general three-dimensional case. Consequently, a general derivation is presented that is applicable to both continuously and discontinuously reinforced composites with arbitrary vector constitutive laws and periodic microstructures. Results are given for thermal and electric problems, effective properties and local field distributions, ordered and random microstructures, as well as complex geometries including woven composites. Comparisons of the theory's predictions are made to test data, numerical analysis, and classical expressions from the literature. Further, classical methods cannot provide the local field distributions in the composite, and it is demonstrated that, as the percolation threshold is approached, their predictions are increasingly unreliable. XXXX It has been observed that the bonding between the fibers and matrix in composite materials can be imperfect. In the context of thermal conductivity, such imperfect interfaces have been investigated in micromechanical models by Dunn and Taya (1993), Duan and Karihaloo (2007), Nan et al. (1997) and Hashin (2001). The present HFGMC micromechanical method, derived for perfectly bonded composite materials governed by vector constitutive laws, can be easily generalized to include the effects of weak bonding between the constituents. Such generalizations, in the context of the mechanical micromechanics problem, involve introduction of a traction-separation law at the fiber/matrix interface and have been presented by Aboudi (1987), Bednarcyk and Arnold (2002), Bednarcyk et al. (2004) and Aboudi et al. (2013) and will be addressed in the future.

micromechanics↗

NASA NTRS2016

Mechanical Model Development for Composite Structural Supercapacitors

Novel composite structural supercapacitor concepts have recently been developed as a means both to store electrical charge and to provide modest mechanical load carrying capability. Double-layer composite supercapacitors are often fabricated by impregnating a woven carbon fiber fabric, which serves as the electrodes, with a structural polymer electrolyte. Polypropylene or a glass fabric is often used as the separator material. Recent research has been primarily limited to evaluating these composites experimentally. In this study, mechanical models based on the Multiscale Generalized Method of Cells (MSGMC) were developed and used to calculate the shear and tensile properties and response of two composite structural supercapacitors from the literature. The modeling approach was first validated against traditional composite laminate data. MSGMC models for composite supercapacitors were developed, and accurate elastic shear/tensile properties were obtained. It is envisioned that further development of the models presented in this work will facilitate the design of composite components for aerospace and automotive applications and can be used to screen candidate constituent materials for inclusion in future composite structural supercapacitor concepts.

Micromechanics↗

NASA NTRS2016

Stochastic-Strength-Based Damage Simulation Tool for Ceramic Matrix and Polymer Matrix Composite Structures

Stochastic-based, discrete-event progressive damage simulations of ceramic-matrix composite and polymer matrix composite material structures have been enabled through the development of a unique multiscale modeling tool. This effort involves coupling three independently developed software programs: (1) the Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC), (2) the Ceramics Analysis and Reliability Evaluation of Structures Life Prediction Program (CARES/ Life), and (3) the Abaqus finite element analysis (FEA) program. MAC/GMC contributes multiscale modeling capabilities and micromechanics relations to determine stresses and deformations at the microscale of the composite material repeating unit cell (RUC). CARES/Life contributes statistical multiaxial failure criteria that can be applied to the individual brittle-material constituents of the RUC. Abaqus is used at the global scale to model the overall composite structure. An Abaqus user-defined material (UMAT) interface, referred to here as "FEAMAC/CARES," was developed that enables MAC/GMC and CARES/Life to operate seamlessly with the Abaqus FEA code. For each FEAMAC/CARES simulation trial, the stochastic nature of brittle material strength results in random, discrete damage events, which incrementally progress and lead to ultimate structural failure. This report describes the FEAMAC/CARES methodology and discusses examples that illustrate the performance of the tool. A comprehensive example problem, simulating the progressive damage of laminated ceramic matrix composites under various off-axis loading conditions and including a double notched tensile specimen geometry, is described in a separate report.

composite structures↗