SEARCH · Engineering Papers

Results for “Matrix elements”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,693 records · Page 94

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Advanced cast aluminum alloys for automotive engine application with superior high-temperature properties

A high fatigue strength aluminum alloy comprises in weight percent copper 3.0-3.5%, iron 0-1.3%, magnesium 0.24-0.35%, manganese 0-0.8%, silicon 6.5-12.0%, strontium 0-0.025%, titanium 0.05-0.2%, vanadium 0.20-0.35%, zinc 0-3.0%, zirconium 0.2-0.4%, a maximum of 0.5% other elements and balance aluminum plus impurities. The alloy defines a microstructure having an aluminum matrix with the Zr and the V in solid solution after solidification. The matrix has solid solution Zr of at least 0.16% after heat treatment and solid solution V of at least 0.20% after heat treatment, and both Cu and Mg are dissolved into the aluminum matrix during the heat treatment and subsequently precipitated during the heat treatment. A process for heat treating an Al—Si—Cu—Mg—Fe—Zn—Mn—Sr-TMs alloy comprises heat treating the alloy to produce a microstructure having a matrix with Zr and V in solid solution after solidification.

36 MATERIALS SCIENCE↗

NASA NTRS1991

Effect of fiber-matrix debonding on notched strength of titanium metal matrix composites

Two specimen configuration of a (0/90)2s SCS-6/Ti-15-3 laminate were tested and analyzed: a center hole (CH) specimen and a double edge notch (DEN) specimen. The two specimen configurations failed at similar stress levels. Two analytical techniques, a 3-D finite-element analysis and a macro-micromechanical analysis were used to predict the overall stress-deformation behavior and the notch-tip fiber-matrix interface stresses in both configurations.

Bigelow, C. A.↗

NASA NTRS1987

Elastic-plastic stress concentrations around crack-like notches in continuous fiber reinforced metal matrix composites

Continuous fiber silicon-carbide/aluminum composite laminates with slits were tested statically to failure. Five different layups were examined: (0) sub 8, (0 sub 2/ + or - 45) sub s, (0/90) sub 2s), (0/ + or - 45/90 sub s), and (+ or - 45) sub 2s. Either a 9.5 or a 19 mm slit was machined in the center of each specimen. The strain distribution ahead of the slit tip was found experimentally with a series of strain gages bonded ahead of the slit tip. A three-dimensional finite element program (PAFAC) was used to predict the strain distribution ahead of the slit tip for several layups. For all layups, except the (0) sub 8, the yielding of the metal matrix caused the fiber stress concentration factor to increase with increasing load. This is contrary to the behavior seen in homogeneous materials where yielding causes the stress concentration to drop. For the (0) sub 8 laminate, yielding of the matrix caused a decrease in the fiber stress concentration. The finite element analysis predicted these trends correctly.

Johnson, W. S.↗

NASA NTRS1989

Elastic-plastic stress concentrations around crack-like notches in continuous fiber reinforced metal matrix composites

Johnson, W. S.↗

NASA NTRS1996

Crack Development in Cross-Ply Laminates Under Uniaxial Tension

This study addresses matrix-dominated failures in carbon fiber/polymer matrix composite laminates in a cross-ply lay-up. The events of interest are interlaminar fracture in the form of transverse cracks in the 90' plies and longitudinal splitting in the 0 deg plies and interlaminar fracture in the form of 0/90 delamination. These events were observed using various nondestructive evaluation (NDE) techniques during static tensile tests. Acoustic emission (AE) x radiography, and edge view microscopy were the principal ones utilized in a real-time environment. A comparison of the NDE results with an analytical model based on the classical linear fracture mechanics concept of strain energy release rate as a criterion for crack growth was performed. The virtual crack closure theory was incorporated with a finite element model to generate strain energy release rate curves for the analytical case. Celion carbon fiber/polyimide matrix (G30-500/PMR-15) was the material tested with cross-ply lay-ups of (0(2)/90(6))s and (0(4)/90(4))s. The test specimens contained thermally induced cracks caused by the high-temperature processing. The analytical model was updated to compensate for the initial damage and to study further accumulation by taking into account the crack interactions. By correlating the experimental and analytical data, the critical energy release rates were found for the observable events of interest.

Gyekenyesi, Andrew L.↗

NASA NTRS1993

Convection equation modeling: A non-iterative direct matrix solution algorithm for use with SINDA

The determination of the boundary conditions for a component-level analysis, applying discrete finite element and finite difference modeling techniques often requires an analysis of complex coupled phenomenon that cannot be described algebraically. For example, an analysis of the temperature field of a coldplate surface with an integral fluid loop requires a solution to the parabolic heat equation and also requires the boundary conditions that describe the local fluid temperature. However, the local fluid temperature is described by a convection equation that can only be solved with the knowledge of the locally-coupled coldplate temperatures. Generally speaking, it is not computationally efficient, and sometimes, not even possible to perform a direct, coupled phenomenon analysis of the component-level and boundary condition models within a single analysis code. An alternative is to perform a disjoint analysis, but transmit the necessary information between models during the simulation to provide an indirect coupling. For this approach to be effective, the component-level model retains full detail while the boundary condition model is simplified to provide a fast, first-order prediction of the phenomenon in question. Specifically for the present study, the coldplate structure is analyzed with a discrete, numerical model (SINDA) while the fluid loop convection equation is analyzed with a discrete, analytical model (direct matrix solution). This indirect coupling allows a satisfactory prediction of the boundary condition, while not subjugating the overall computational efficiency of the component-level analysis. In the present study a discussion of the complete analysis of the derivation and direct matrix solution algorithm of the convection equation is presented. Discretization is analyzed and discussed to extend of solution accuracy, stability and computation speed. Case studies considering a pulsed and harmonic inlet disturbance to the fluid loop are analyzed to assist in the discussion of numerical dissipation and accuracy. In addition, the issues of code melding or integration with standard class solvers such as SINDA are discussed to advise the user of the potential problems to be encountered.

Schrage, Dean S.↗

NASA NTRS2021

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↗

DOE OSTI2023

Characterization of actinide abundances and isotopic compositions by HR-ICP-MS. Part 2: Results from actinide doping studies

Previously, we presented actinide isotopic and elemental data from fallout melt glass that was measured using high-resolution inductively coupled plasma mass spectrometry (HR-ICP-MS). Direct comparison between these measurements and ‘gold standard’ data obtained by multi-collector ICP-MS showed broad overlap, indicating that HR-ICP-MS is a potentially valuable technique for producing actinide elemental and isotopic data on a relatively rapid timescale. To test the usefulness of this technique further, we doped varying amounts of uranium certified reference materials (CRMs) into a rhyolitic rock standard to establish the effects of uranium concentration and isotopic composition on the accuracy of uranium isotopic analyses by HR-ICP-MS. This also enabled us to quantify peak tailing effects from 238 U on the measurement of 239 Pu and 237 Np, which in turn allows us to constrain correction factors based on measured 237 Np/ 238 U and 239 Pu/ 238 U ratios. A second doping study involved the addition of a mixed actinide standard into three samples with different matrix compositions (termed soil, city, and seawater) to assess whether sample chemistry affects the accuracy and precision of these analyses. Results suggest that this is not the case. Systematic offsets were not observed in elemental or isotopic data derived from the three matrix samples. Results indicate that useful actinide isotopic data can be obtained from whole rock solutions by HR-ICP-MS. Our findings also have implications for solid sampling techniques such as laser ablation ICP-MS, which do not require sample dissolution.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2024

A Low-Rank QTT-based Finite Element Method for Elasticity Problems

We present an efficient and robust numerical algorithm for solving the linear elasticity problem that combines the Quantized Tensor Train format and a domain partitioning strategy. This approach makes it possible to solve the linear elasticity problem on a computational domain that is more general than a square. By integrating Z-ordering and subdomain concatenation, our method substantially decreases memory usage and achieves a notable reduction in rank compared to established Finite Element implementations like the FEniCS platform. This efficiency is maintained while still guaranteeing exponential convergence with respect to the number of degrees of freedom. This performance gain, however, requires a fundamental rethinking of how core finite element operations are implemented. This includes changes to mesh discretization, node and degree of freedom ordering, stiffness matrix and internal nodal force assembly, and the execution of algebraic matrix-vector operations. In this work, we discuss all these aspects in detail and assess the method’s performance in the numerical approximation of three representative test cases.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2015

Development and Characterization of a Rate-Dependent Three-Dimensional Macroscopic Plasticity Model Suitable for Use in Composite Impact Problems

Several key capabilities have been identified by the aerospace community as lacking in the material/models for composite materials currently available within commercial transient dynamic finite element codes such as LS-DYNA. Some of the specific desired features that have been identified include the incorporation of both plasticity and damage within the material model, the capability of using the material model to analyze the response of both three-dimensional solid elements and two dimensional shell elements, and the ability to simulate the response of composites composed with a variety of composite architectures, including laminates, weaves and braids. In addition, a need has been expressed to have a material model that utilizes tabulated experimentally based input to define the evolution of plasticity and damage as opposed to utilizing discrete input parameters (such as modulus and strength) and analytical functions based on curve fitting. To begin to address these needs, an orthotropic macroscopic plasticity based model suitable for implementation within LS-DYNA has been developed. Specifically, the Tsai-Wu composite failure model has been generalized and extended to a strain-hardening based orthotropic plasticity model with a non-associative flow rule. The coefficients in the yield function are determined based on tabulated stress-strain curves in the various normal and shear directions, along with selected off-axis curves. Incorporating rate dependence into the yield function is achieved by using a series of tabluated input curves, each at a different constant strain rate. The non-associative flow-rule is used to compute the evolution of the effective plastic strain. Systematic procedures have been developed to determine the values of the various coefficients in the yield function and the flow rule based on the tabulated input data. An algorithm based on the radial return method has been developed to facilitate the numerical implementation of the material model. The presented paper will present in detail the development of the orthotropic plasticity model and the procedures used to obtain the required material parameters. Methods in which a combination of actual testing and selective numerical testing can be combined to yield the appropriate input data for the model will be described. A specific laminated polymer matrix composite will be examined to demonstrate the application of the model.

Impact↗

DOE OSTI2022

Using θ' interfaces as templates for planar L1 2 precipitation in AlCuMnZr alloys

Controlled Mn and Zr additions to Al-Cu alloys have allowed for the improved retention of mechanical properties after extended 350°C exposures by stabilizing the main strengthening θ' (Al 2 Cu) phase. Ultimately, θ'/L1 2 (Al 3 Zr) co-precipitate formation stabilizes θ' most effectively; however, Zr diffuses sluggishly and has low solubility in aluminum castings. Increasing the Zr segregation rate would allow for faster and more effective θ'/L1 2 co-precipitation. It is demonstrated that the Zr segregation rate is faster when the Zr matrix content is higher. A much higher Zr matrix content was achieved by rapid cooling during additive manufacturing (AM) that produces θ'/L1 2 co-precipitation faster, which is shown by scanning transmission electron microscopy and atom probe tomography experiments. It was also found that Zr continuously segregates to θ' interfaces up to the most aggressive heat treatment studied such that planar L1 2 precipitates remain after the metastable θ' dissolves. In this manner, we demonstrate that θ' coherent interfaces serve as perfect templates to form stable planar L1 2 precipitates that can provide strength at higher temperatures than traditional θ' strengthened AlCu alloys. This work introduces an alloy design strategy that uses metastable precipitates to quickly nucleate and grow co-precipitates with a desired geometry that contain slow diffusing elements. These ideas can be applied to engineer more heat resistant alloys by taking advantage of high solute matrix contents enabled by rapid cooling during additive manufacturing.

36 MATERIALS SCIENCE↗

NASA NTRS1996

Computational Simulation of Continuous Fiber-Reinforced Ceramic Matrix Composites Behavior

This report describes a methodology which predicts the behavior of ceramic matrix composites and has been incorporated in the computational tool CEMCAN (CEramic Matrix Composite ANalyzer). The approach combines micromechanics with a unique fiber substructuring concept. In this new concept, the conventional unit cell (the smallest representative volume element of the composite) of the micromechanics approach is modified by substructuring it into several slices and developing the micromechanics-based equations at the slice level. The methodology also takes into account nonlinear ceramic matrix composite (CMC) behavior due to temperature and the fracture initiation and progression. Important features of the approach and its effectiveness are described by using selected examples. Comparisons of predictions and limited experimental data are also provided.

Murthy, Pappu L. N.↗

DOE OSTI2021

A Multiscale Modeling Approach to Cryocompressed Hydrogen Storage Pressure Vessels – Part II: Constitutive Modeling and Finite Element Analysis

Subjected to thermomechanical cycling, composite cryogenic hydrogen (H2) storage pressure vessels experience high stresses in the carbon fiber (CF)/epoxy overwrap which can lead to vessel failure due to a combination of degradation mechanisms such as matrix cracking, fiber/matrix debonding, delamination, and fiber rupture. The present paper is the second part of a set of two articles on this subject that addresses analysis and design of these pressure vessels. Predictive finite element (FE) modeling capabilities have been used to support a material acceptance process that aims at providing guidance to evaluate specialty resins, vessel liner options, and CF composites through thermomechanical testing. We have applied a multiscale modeling approach recently developed (B.N. Nguyen et al., International Journal of Hydrogen Energy, 2019, https://doi.org/10.1016/j.ijhydene.2019.09.200) to perform three-dimensional (3D) FE analysis and design of a cryo-compressed H2 storage pressure vessel for its operation in a large temperature range from room to cryogenic temperatures. The developed approach termed the micro-meso-macro approach determines the laminar stresses in different layers of the vessel in addition to constituent (i.e., fiber and matrix) stresses that are important to the design of the composite overwrap for its layup and material combination to reduce the risk of vessel failure.

Nguyen, Ba Nghiep↗

NASA NTRS1989

Improved performance in NASTRAN (R)

Three areas of improvement in COSMIC/NASTRAN, 1989 release, were incorporated recently that make the analysis program run faster on large problems. Actual log files and actual timings on a few test samples that were run on IBM, CDC, VAX, and CRAY computers were compiled. The speed improvement is proportional to the problem size and number of continuation cards. Vectorizing certain operations in BANDIT, makes BANDIT run twice as fast in some large problems using structural elements with many node points. BANDIT is a built-in NASTRAN processor that optimizes the structural matrix bandwidth. The VAX matrix packing routine BLDPK was modified so that it is now packing a column of a matrix 3 to 9 times faster. The denser and bigger the matrix, the greater is the speed improvement. This improvement makes a host of routines and modules that involve matrix operation run significantly faster, and saves disc space for dense matrices. A UNIX version, converted from 1988 COSMIC/NASTRAN, was tested successfully on a Silicon Graphics computer using the UNIX V Operating System, with Berkeley 4.3 Extensions. The Utility Modules INPUTT5 and OUTPUT5 were expanded to handle table data, as well as matrices. Both INPUTT5 and OUTPUT5 are general input/output modules that read and write FORTRAN files with or without format. More user informative messages are echoed from PARAMR, PARAMD, and SCALAR modules to ensure proper data values and data types being handled. Two new Utility Modules, GINOFILE and DATABASE, were written for the 1989 release. Seven rigid elements are added to COSMIC/NASTRAN. They are: CRROD, CRBAR, CRTRPLT, CRBE1, CRBE2, CRBE3, and CRSPLINE.

Chan, Gordon C.↗

NASA NTRS1990

Analytical and experimental studies on creep behavior of polymeric matrix composites

The creep behavior of graphite/epoxy composites is studied both analytically and experimentally. In the analytical study, a special finite element procedure was developed for the accurate and efficient analysis of creep response in anisotropic materials. This procedure was used to study the stress and strain distributions and histories in composite laminates containing circular holes. In the experimental study, creep tests were performed to investigate the time-dependent response of graphite/epoxy composites at elevated temperature. The moire interferometry technique was employed to determine the deformation histories in notched composites. The experimental results were compared with the analytical predictions and good agreement was observed.

Lin, K. Y.↗

NASA NTRS2017

Implementation of a Tabulated Failure Model Into a Generalized Composite Material Model Suitable for Use in Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased use in the aerospace and automotive communities. The aerospace community has identified several key capabilities which are currently lacking in the available material models in commercial transient dynamic finite element codes. To attempt to improve the predictive capability of composite impact simulations, a next generation material model is being developed for incorporation within the commercial transient dynamic finite element code LS-DYNA. The material model, which incorporates plasticity, damage and failure, utilizes experimentally based tabulated input to define the evolution of plasticity and damage and the initiation of failure as opposed to specifying discrete input parameters such as modulus and strength. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. For the damage model, a strain equivalent formulation is used to allow for the uncoupling of the deformation and damage analyses. For the failure model, a tabulated approach is utilized in which a stress or strain based invariant is defined as a function of the location of the current stress state in stress space to define the initiation of failure. Failure surfaces can be defined with any arbitrary shape, unlike traditional failure models where the mathematical functions used to define the failure surface impose a specific shape on the failure surface. In the current paper, the complete development of the failure model is described and the generation of a tabulated failure surface for a representative composite material is discussed.

polymer matrix composites↗

NASA NTRS2019

Effect of Phenolic Matrix Microcracking on the Structural Response of a 3-D Woven Thermal Protection System

The effect of microcracking in the phenolic matrix of a three-dimensional woven thermal protection system (TPS) and the resulting material stiffness reduction was studied via a comparison of finite element analysis results from a linear analysis and an iterative linear analysis. A TPS is necessary to protect space vehicles from the aerodynamic heating of planetary entry. The Heatshield for Extreme Entry Environment Technology (HEEET) project has developed a TPS for use in high heat-flux and pressure missions. The material is a dual-layer continuous dry weave, which is then infiltrated with a low-density phenolic resin matrix to form a rigid ablator. The phenolic resin matrix does not have structural load transfer requirements, and testing has shown that the phenolic resin can fully satisfy thermal requirements when the matrix contains microcracks. Due to high stresses in the through-the-thickness direction of the material, phenolic microcracks may form in the matrix material, which would result in a reduction of stiffness. An exploratory study was conducted to determine if reduction in material stiffness would change the load paths and/or decrease the structural margins. A comparison was performed between a linear finite element analysis that did not take into account phenolic microcracking and an iterative linear finite element analysis that accounted for propagation of phenolic microcracking. Four subcases using varying assumptions were analyzed and the results indicate that the assumed strength at which the phenolic microcracking propagates was the critical parameter for determining the extent of microcracking in the phenolic matrix. Phenolic microcracking does not have an adverse effect on the structural response of the test article and is not a critical failure.

Langston, Sarah L.↗

NASA NTRS2011

Finite Element Model Characterization Of Nano-Composite Thermal And Environmental Barrier Coatings

Thermal and environmental barrier coatings have been applied for protecting Si based ceramic matrix composite components from high temperature environment in advanced gas turbine engines. It has been found that the delamination and lifetime of T/EBC systems generally depend on the initiation and propagation of surface cracks induced by the axial mechanical load in addition to severe thermal loads. In order to prevent T/EBC systems from surface cracking and subsequent delamination due to mechanical and thermal stresses, T/EBC systems reinforced with nano-composite architectures have showed promise to improve mechanical properties and provide a potential crack shielding mechanism such as crack bridging. In this study, a finite element model (FEM) was established to understand the potential beneficial effects of nano-composites systems such as SiC nanotube-reinforced oxide T/EBC systems.

Yamada, Yoshiki↗