Search NASA⌕ Search

Engineering topics

Evan J. Pineda

Publications and source records attributed to Evan J. Pineda.

At least 19 records

Out of Plane Stresses Using a Beam-Based Elastic Foundation Bonded Joint Model

Efficient, semi-analytical models of bonded joints are utilized for initial joint design and sizing due to their convenience and speed. Many of these models are based on beam and plate theories, which are most applicable when in-plane stresses/strains dominate. While this can be a sensible assumption for obtaining the deformations of the joint, out of plane (transverse) stresses are often of critical importance for failure in bonded joints composed of fiber reinforced composite materials. Since through-thickness strains are neglected in these semi-analytical theories, the present study shows a method to back-calculate the out of plane stresses based on stress balance equations, which use the in-plane stresses and displacements that have already been obtained. The method was implemented within the Joint Element Designer (JED) software, and verification and validation examples have been shown. As is typical for these types of models, corner stress singularities are not captured by the JED model, but the general trends of the stress and deformation fields agreed with a higher fidelity finite element model. As such, this method for obtaining out of plane stresses, when coupled with the JED code, is very useful for the complete stress analysis of bonded composite joints.

composites↗

Multiscale Progressive Failure Analysis of 3D Woven Composites

Application of three-dimensional (3D) woven composites is growing as an alternative to the use of ply-based composite materials. However, the design, analysis, modeling, and optimization of these materials is more challenging due to their complex and inherently multiscale geometries. Herein, a multiscale modeling procedure, based on efficient, semi-analytical micromechanical theories rather than the traditional finite element approach, is presented and applied to a 3D woven carbon–epoxy composite. A crack-band progressive damage model was employed for the matrix constituent to capture the globally observed nonlinear response. Realistic microstructural dimensions and tow-fiber volume fractions were determined from detailed X-ray computed tomography (CT) and scanning electron microscopy data. Pre-existing binder-tow disbonds and weft-tow waviness, observed in X-ray CT scans of the composite, were also included in the model. The results were compared with experimental data for the in-plane tensile and shear behavior of the composite. The tensile predictions exhibited good correlations with the test data. While the model was able to capture the less brittle nature of the in-plane shear response, quantitative measures were underpredicted to some degree.

3D woven composites↗

Residual Stresses Induced due to Curing of the Bulk Matrix in a Simplified Three-Dimensional (3D) Woven Repeating Unit Cell

It has been observed that 3D woven materials are highly susceptible to processes induced defects including tow misalignment, intratow and intertow cracking and/or voids. Adequate process modeling of 3D woven Polymer Matrix Composites (PMCs) is necessary to predict and estimate the effects of the manufacturing process on these defects. A preliminary step towards this is achieved with the commercial Finite Element Analysis (FEA) tool Abaqus complemented with written user subroutines to account for the effect of shrinkage and thermoelastic properties evolution as a function of the degree of cure of the bulk matrix phase in a 3D woven composite repeating unit cell (RUC). As a pathfinder, the process modeling framework is demonstrated on a simplified 3D woven RUC of an AS4/RIMR 135 system. Curing simulations are performed on the 3D woven RUC to observe the processing induced residual stresses. The results of the analysis show that high stress concentrations, present within the 3D intertow matrix bordering the binder tow, coincide with X-Ray computed tomography (X-Ray-CT) scan data depicting cracking in the same location

3D Woven↗

Relevence of 3D Simulations and Sandwich Core Topology for the Modeling of Honeycomb Core Sandwich Panels Undergoing Mode I Delamination

A recently developed cohesive zone traction-separation law, which includes the effects of fiber bridging in a novel way, is extended from 2D to 3D. The proposed cohesive model is applied to low fidelity (i.e. homogenized core representation) and high fidelity (i.e. directly accounting for the core topology) finite element models of a composite panel comprised of carbon fiber reinforced plastic facesheets and a honeycomb sandwich core. This enables the investigation of 2D to 3D parameter transferability, width-dependent effects such as thumbnail-shaped crack growth, and the verification of plane strain / plane stress assumptions. A pronounced curvature of the initial interface-related crack front is observed, while the bridging-related crack front is straight. Furthermore, it is found that the cohesive parameters can easily be transferred from 2D to 3D under plane stress assumptions, but not under plane strain assumptions. The numerical predictions are compared to experimental load-displacement and R-curves.

Composites↗

Application of the NASA Multiscale Analysis Tool: Multiscale Integration and Interoperability

In order to demonstrate NASMAT’s multiscale operability, a series of illustrative examples will be presented that focus on the application of NASMAT to practical problems. First, the multiscale integration and data recursion is demonstrated by performing a multiscale analysis using only built-in micromechanics methods. NASMAT’s integration is then highlighted by running a multiscale analysis where an external finite element software calls NASMAT. In this case, at each integration point within the finite element model, a local NASMAT analysis is performed to account for failure behavior at the constituent scale. In a similar example, an external program is called from within NASMAT. This case would be relevant for a user wanting to implement an outside micromechanics technique. A combination of these examples is then presented to further illustrate the code’s flexibility when interfacing with outside codes in a multiscale framework. For all examples, data is presented using a custom-developed visualization tool. Additional potential use cases are also addressed. Finally, the plan for upcoming features and added capabilities is discussed.

NASMAT↗

Application of the NASA Multiscale Analysis Tool: Multiscale Integration and Interoperability

The NASA Multiscale Analysis Tool (NASMAT) was developed recently to allow a wide variety of multiscale analysis problems to be effectively and efficiently solved. The architecture of NASMAT was established specifically to enable parallelized, “plug-and-play” functionality to reduce the complexity associated with adding new features to the code in the future and to allow end users to rapidly implement and evaluate user-defined capabilities. Additionally, the tool utilizes recursive data structures and subroutines to allow for an arbitrary number of length scales when performing multiscale analyses of heterogeneous materials. These features permit the rapid integration of user-defined capabilities (e.g., a material model, micromechanics approach, or failure theory) at all stages within a NASMAT calculation while leveraging built-in techniques where needed. Additionally, these features allow NASMAT to both be called from an external program as well as call an external program. This paper specifically focuses on the multiscale integration and interoperability of NASMAT with other analysis techniques through an illustrative, multiscale analysis of a 3D woven polymer matrix composite (PMC).

NASMAT↗

Software Architecture and Hierarchy of the Nasa Multiscale Analysis Tool

The NASA Multiscale Analysis Tool (NASMAT) serves as a state-of-the-art, “plug and play,” software package which utilizes multiscale recursive micromechanics as a platform for massively multiscale modeling for hierarchical materials and structures subjected to thermomechanical loads on high performance computing systems. This paper is intended to give an overview of the design of NASMAT and how the design supports modularity, upgradability and maintainability, interoperability, and utility. First, the software architecture and hierarchy will be explored. Details on each of the 11 NASMAT procedures and the arrangement of NASMAT data will be presented. Finally, application program interfaces (APIs) that were developed to facilitate the communication of NASMAT with other programs will be described.

multiscale modeling↗

Macro-scale Testing and Micromechanics Modeling of Fracture Behaviors for Boron Carbide Composites with Hierarchical Microstructures

With a unique combination of properties including high hardness, low density, wear and corrosion resistance, thermal stability, high neutron absorption, and semi-conductivity, boron carbide (B4C) is a candidate material for various engineering applications that involve extreme environment. The current applications of boron carbide, however, are limited by its intrinsic brittleness due to its strong covalent bonding. To toughen boron carbide, in this work hierarchical microstructure designs was used to provide multiple toughening mechanisms including crack deflection/bridging, micro-crack toughening, etc. Using field assisted sintering, B4C composites with hierarchical microstructure features including graphite platelets, micron and sub-micron sized TiB2 reinforcements were fabricated. The fracture toughness of fabricated B4C composites were previously measured at micro-scale using micro-indentation followed by post-testing microstructure inspection. However, questions including whether the fracture toughness enhancement measured at micro-scale can translate to macro-scale mechanical properties, and what are the fundamental mechanisms behind observed fracture toughness enhancement, remain to be answered. In this study, the fabricated B4C composites were tested using standardized four-point bending method to obtain fracture toughness at macro-scale. In addition, micromechanics modeling was conducted using MAC/GMC code and crack-band model to study the effect of residual stress and weak interphases on fracture behaviors of B4C composites reinforced with TiB2 particles. Through standardized four-point bending tests, fracture toughness enhancements up to 2.85, 3.32, and 3.65 MPa∙m(exp 1/2) (from 2.38 MPa∙m(exp 1/2)) were achieved for B4C composites with graphite platelets addition (micro/nano B4C), with TiB2 formation (micro B4C-TiB2), and with both graphite and TiB2 addition (micro/nano B4C-TiB2) respectively. Micromechanics modeling indicated that introduction of thermal residual stress and weak interphases caused enhanced micro-cracking behavior and resulted in the observed fracture toughness enhancement. These results furthered understanding about the mechanical behaviors at macro-scale and the mechanisms behind observed fracture toughness enhancement for B4C composites with hierarchical microstructures and can provide reference data for the future design of B4C composites with optimized microstructures for further fracture toughness enhancement.

Boron Carbide↗

Geometric Sensitivity of Residual Fields in Metal Additive Manufacturing

The sensitivity of geometrical parameters on the residual stress and distortion as computed by a metal additive manufacturing simulation is presented. A nozzle-like geometry was chosen with its length, wall thickness, and outer radius varied. A Latin hypercube sampling of 300 specimens produced the combinations of the geometrical parameters. Full thermomechanical simulations predicted the residual fields. Statistical analysis including regression curves shows that there is a strong correlation between these residual fields and the geometry, hinting that it may be possible to predictably influence the residual displacements and residual stresses by varying the overall geometry.

Multiscale Modeling↗

Artificial Generation of 2-D Fiber Reinforced Composite Microstructures with Statistically Equivalent Features

Fiber reinforced composites are used widely for their high strength and low weight advantages in various aerospace and automotive applications. While their use may be sought after, modeling of these material requires increasing fidelity at the lower scales to capture accurate material behavior under loading. The first steps in creating statistically equivalent models to real life cases is developing a method of rapid evaluation and artificial microstructure generation. The outlined work is capable of tracking microscale fiber positions and determining regions of localized volume fraction extrema (high and low end). Groupings of high and low volume fraction regions are called clusters and their geometry is used to characterize the microstructure. These cluster features can be evaluated for both artificial models and actual scans, allowing correlation to be established which can ultimately be used to regenerate statistically equivalent models. The results of this work show that if one feature is to be correlated, a model can be generated which matches almost exactly. But once more features are equally taken into account, the regeneration loses accuracy.

micromechanics↗

Thermal Conductivity of 3D Woven Composite Thermal Protection System Materials via Multiscale Recursive Micromechanics

Multiscale micromechanics predictions have been made for the effective thermal conductivities and local thermal fields for a novel 3D woven composite thermal protection system material. The Multiscale Recursive Micromechanics approach, which enables micromechanics models to call other micromechanics models (or themselves recursively) to consider finer and finer length scales, has been employed. The multiscale model uses a recently developed version of the High-Fidelity Generalized Method of Cells micromechanics theory at each of three length scales. Results focus on the impact of the microstructural geometry representation at each length scale on the material’s effective thermal conductivity, along with the local thermal flux and temperature fields induced in the microstructures.

Composites↗

Multiscale Thermal Conductivity Modeling of 3D Woven Composite Thermal Protection System Materials

Novel thermal protection systems (TPS) for re-entry and hypersonic vehicles utilize 3D woven composite materials with blended carbon filament tows (yarns) and a phenolic resin matrix. The presence of continuous carbon fibers provides superior stiffness and strength compared to legacy TPS, while the 3D weave pattern provides a great deal of design flexibility for both in-plane and through thickness thermomechanical behavior. However, modeling of 3D woven composites is notoriously complex and numerically intensive. The present investigation utilizes a unique ultra-efficient approach, known as Multiscale Recursive Micromechanics (MsRM), wherein recursive semi-analytical micromechanics methods are employed at various length scales within the composite. The MsRM approach has recently been extended to solve for the effective thermal conductivity as well as the local temperature and thermal flux fields throughout the composite. Presented model results focus on the impact of the microstructural geometry representation at each length scale on the material’s effective thermal conductivity, along with the local thermal flux and temperature fields induced in the microstructures, for a novel 3D woven TPS.

Composite↗

Estimation of 3D Woven Design Sensitivities Using a Rapid Multiscale Analysis Technique

Highly-refined finite element models of three-dimension (3D) woven composite systems currently require excessive computational demands that limit their use in sensitivity analysis, uncertainty quantification, and optimization. An alternative analysis methodology was developed using the NASA Multiscale Analysis Tool (NASMAT) where multiscale models of a 3D woven composite (including inter-tow matrix voids and constituent failure) can be completed on a single central processing unit(CPU)on the order of ~30 s. To develop inputs and validation data for the NASMAT model, coupon and acid-digesting testing and x-ray computed tomography were performed. The NASMAT inputs were parameterized using a set of 25 input variables and distributions. These inputs were randomly sampled to generate a total of 100,000 NASMAT analyses that could be used to understand the influence of different material and geometric properties on the warp and weft-direction stiffness and strength. These analyses (including pre/post-processing) were performed in less than eight hours on a 120 CPU cluster. The computational efficiency of the NASMAT model enabled a sensitivity analysis to be performed, and dominant input variables were able to be identified. Key results were consistent with theoretical and experimental observations for the specific 3D woven system studied in this work.

NASMAT↗

A Partitioned -Task Parallel Implementation of the NASA Multiscale Analysis Tool for High Performance Computing

The NASA Multiscale Analysis Tool (NASMAT) is a “plug and play” software package that allows users to conduct massively multiscale modeling of hierarchical and nonlinear materials. This work extends the scalability and improves the High Performance Computing friendliness of NASMAT by adopting a Partitioned Task-Parallel approach. Interoperability of NASMAT with external software is enhanced through preCICE, a open source library for multiphysics coupling in a partitioned manner. Enhancement through preCICE allows for easy integration of NASMAT to other macro solvers and dissociates the parallelization strategy adopted within NASMAT from the macro solver. The task-parallel framework based on Master-Worker approach is implemented as the parallelization scheme. The scheme accounts for hierarchy of multiple scales (task-dependence) and heterogeneous nature (dynamic load balancing) of computations. The applicability and scalability of the framework will be evaluated by analyzing large scale engineering problems through massively multiscale methods.

NASMAT↗

Modeling and Characterization of Crushable Composite Structures

The crushing response and the energy absorption capacity of a quasi-isotropic C-shaped composite panel is presented to simulate the progressive damage and failure mechanisms involved in the component testing. Crushable composite structures are widely used for impact applications, as because of their great high-energy absorption capability and lower weight. Simulating the crush characteristics of composite structures is challenging, as the conventional tools fail to accurately represent the failure mechanisms at the component level. Here, a cost-effective predictive tool ‘CZone’ is applied for simulating the crushing response of composite materials. CZone for Abaqus (CZA) is an add-on capability for Abaqus/Explicit and combines with the commercially available impact modeling capabilities included within the Abaqus finite element software package.

crush zone↗

Multiscale Modeling of Thermoplastics Using Atomistic-informed Micromechanics

A multiscale model was developed for predicting the thermoelastic behavior of semi-crystalline thermoplastic materials for composite aerospace applications. At the highest scale containing the semi-crystalline spherulite in an amorphous matrix, the generalized method of cells, or high fidelity method of cells, was used to perform the homogenization calculations to obtain the effective properties. Models were developed assuming a cubic, or spherical shape, for the spherulite to understand if the morphology of the spherulite affects the effective thermoelastic properties. The generalized method of cells was used to model at the repeating unit cells at the subscales of the microstructure including the lamellae stacks and granular crystal blocks. The scales are integrated using the multiscale micromechanics method in the NASA Multiscale Analysis Tool. Data from molecular dynamics simulations were used as inputs for the amorphous and crystalline constituents. Convergence studies were performed to determine the best level of discretization for the repeating unit cell at the highest scale. Effective Young’s modulus, shear modulus, Poisson’s ratio, coefficient of thermal expansion, and thermal conductivity were predicted for polyether ether ketone and polyether ketone ketone, and very good agreement between the model utilizing the cubic spherulite and the experimental data, where available, was observed for polyether ketone ketone. Normalization of the data for the bulk polyether ketone ketone, against amorphous data, improved the predictions as compared to experimental data. Overall, the high fidelity method of cells predicted a stiffer response then the generalized method of cells as the crystallinity was increased. The shape of the spherulite had a minimal effect on the predicted bulk properties of the polymers.

thermoplastics↗

Scaling Electric Machines to a Megawatt and Material Options

Megawatt (MW) electric aircraft propulsion (EAP) is seen as a significant contributor toward achieving the goals set forth by the Sustainable Flight National Partnership. A large part of enabling MW EAP is developing specific-power-dense electric machines. As specific-power-dense electric machines are scaled up from kW to MW power levels, the thermal stresses on the machines increase in both magnitude and performance-affecting characteristics. This is particularly true for the stators of these machines. Analysis via thermal resistance network modeling and multiscale modeling reveals that increasing amounts of heat will be trapped in the stator windings as the power levels increase. The challenges this presents can be addressed through material advancements whereby materials gain multifunctionality. Specifically, the electrical insulation and potting materials, along with the electrical conductor, that compose the stator slot must work together (gain multifunctionality) to relieve the increased thermal stress. Materials research at the NASA Glenn Research Center points to some useful solutions in this trade space.

Electric Machine↗