Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix Models”

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.

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

At least 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Analysis of damage in MMC components using an internal state variable model

A metal-matrix composite (MMC) model was developed which includes the concept of damage evolution. The evolution of damage is assumed to be governed by a Kachanov-type equation. This viscoplastic damage model was implemented in the finite element code, MARC. Both uniaxial (creep) and multiaxial (an internally pressurized thick-walled cylinder) problems were analyzed using this implementation. Some preliminary results are presented which consider monotonic (constant) loadings. The creep curves including damage for four fiber orientations are presented. As expected, the minimum creep occurs when load is applied in a direction parallel to the fibers. The tangential strains at the inner radius of a thick-walled MMC-cylinder for four fiber orientations are shown with damage included. The cylinder exhibits the maximum creep resistance when the fibers are oriented in the circumferential direction, perpendicular to the axis of the cylinder. Time-to-failure for the thick-walled cylinder for the same fiber orientation angles is also shown. As expected, the life of the cylinder can be increased by orientating the fibers in the circumferential direction, perpendicular to the axis of the cylinder. The results, although qualitative, indicate that significant benefits in creep-resistance and service life can be achieved by using MMC materials as structural materials for high-temperature design.

Arya, V. K.↗
DOE OSTI2020

Fast matrix algebra for Bayesian model calibration

In Bayesian model calibration, evaluation of the likelihood function usually involves finding the inverse and determinant of a covariance matrix. When Markov Chain Monte Carlo (MCMC) methods are used to sample from the posterior, hundreds of thousands of likelihood evaluations may be required. In this paper, we demonstrate that the structure of the covariance matrix can be exploited, leading to substantial time savings in practice. Here, we also derive two simple equations for approximating the inverse of the covariance matrix in this setting, which can be computed in near-quadratic time. The practical implications of these strategies are demonstrated using a simple numerical case study and the "quack" R package. For a covariance matrix with 1000 rows, application of these strategies for a million likelihood evaluations leads to a speedup of roughly 4000 compared to the naive implementation

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1994

Prediction of stress-strain response of SCS-6/Timetal-21S subjected to a hypersonic flight profile

Thermomechanical response of a cross-ply SCS-6/Timetal-21S composite subjected to a generic hypersonic flight profile with the temperature ranging from -130 C to 816 C was evaluated experimentally and analytically. A two dimensional micromechanical anlaysis, VISCOPLY, was used to predict the stress-strain response of the laminate and of the constituents in each ply during thermomechanical loading conditions. In the analysis, the fiber was modeled as elastic with transverse orthotropic and temperature dependent properties and the matrix was modeled using a thermoviscoplastic constitutive relation. The fiber transverse modulus was reduced in the analysis to simulate fiber-matrix interface failure. Reasonable agreement was found between measured and predicted laminate stress-strain response when fiber-matrix debonding was modeled.

Mirdamadi, Massoud↗
NASA NTRS2017

Fracture-Based Mesh Size Requirements for Matrix Cracks in Continuum Damage Mechanics Models

This paper evaluates the ability of progressive damage analysis (PDA) finite element (FE) models to predict transverse matrix cracks in unidirectional composites. The results of the analyses are compared to closed-form linear elastic fracture mechanics (LEFM) solutions. Matrix cracks in fiber-reinforced composite materials subjected to mode I and mode II loading are studied using continuum damage mechanics and zero-thickness cohesive zone modeling approaches. The FE models used in this study are built parametrically so as to investigate several model input variables and the limits associated with matching the upper-bound LEFM solutions. Specifically, the sensitivity of the PDA FE model results to changes in strength and element size are investigated.

Leone, Frank A.↗
NASA NTRS1999

Strain Rate Dependent Deformation and Strength Modeling of a Polymer Matrix Composite Utilizing a Micromechanics Approach

Potential gas turbine applications will expose polymer matrix composites to very high strain rate loading conditions, requiring an ability to understand and predict the material behavior under extreme conditions. Specifically, analytical methods designed for these applications must have the capability of properly capturing the strain rate sensitivities and nonlinearities that are present in the material response. The Ramaswamy-Stouffer constitutive equations, originally developed to analyze the viscoplastic deformation of metals, have been modified to simulate the nonlinear deformation response of ductile, crystalline polymers. The constitutive model is characterized and correlated for two representative ductile polymers. Fiberite 977-2 and PEEK, and the computed results correlate well with experimental values. The polymer constitutive equations are implemented in a mechanics of materials based composite micromechanics model to predict the nonlinear, rate dependent deformation response of a composite ply. Uniform stress and uniform strain assumptions are applied to compute the effective stresses of a composite unit cell from the applied strains. The micromechanics equations are successfully verified for two polymer matrix composites. IM7/977-2 and AS4/PEEK. The ultimate strength of a composite ply is predicted with the Hashin failure criteria that were implemented in the composite micromechanics model. The failure stresses of the two composite material systems are accurately predicted for a variety of fiber orientations and strain rates. The composite deformation model is implemented in LS-DYNA, a commercially available transient dynamic explicit finite element code. The matrix constitutive equations are converted into an incremental form, and the model is implemented into LS-DYNA through the use of a user defined material subroutine. The deformation response of a bulk polymer and a polymer matrix composite are predicted by finite element analyses. The results compare reasonably well to experimental values, with some discrepancies. The discrepancies are at least partially caused by the method used to integrate the rate equations in the polymer constitutive model.

Goldberg, Robert K.↗
NASA NTRS1993

A detailed study of mean-flow solutions for stability analysis of transitional flows

A finite-volume upwind-difference parabolized Navier-Stokes code is utilized to obtain laminar mean-flow solutions at Mach 3.5 on a half-angle cone of 5 deg at an angle-of-attack of 2 deg. A detailed study is conducted on this configuration; the main focus is the velocity profiles in the leeward and windward symmetry planes at various axial locations. Comparisons of the solution profiles are made with both a central-difference code that incorporates scalar and matrix dissipation models and another state-of-the-art upwind-difference finitevolume code. The results obtained emphasize the importance of using matrix dissipation models for schemes that require explicit artificial dissipation. These results also illustrate the accuracy and efficiency of the planeby-plane marching procedure for computing mean-flow solutions for predicting the onset of transition with linear instability.

Ramakrishnan, R.↗
NASA NTRS1994

Modeling of outgassing and matrix decomposition in carbon-phenolic composites

Work done in the period Jan. - June 1994 is summarized. Two threads of research have been followed. First, the thermodynamics approach was used to model the chemical and mechanical responses of composites exposed to high temperatures. The thermodynamics approach lends itself easily to the usage of variational principles. This thermodynamic-variational approach has been applied to the transpiration cooling problem. The second thread is the development of a better algorithm to solve the governing equations resulting from the modeling. Explicit finite difference method is explored for solving the governing nonlinear, partial differential equations. The method allows detailed material models to be included and solution on massively parallel supercomputers. To demonstrate the feasibility of the explicit scheme in solving nonlinear partial differential equations, a transpiration cooling problem was solved. Some interesting transient behaviors were captured such as stress waves and small spatial oscillations of transient pressure distribution.

Mcmanus, Hugh L.↗
NASA NTRS1985

Report of the LSPI/NASA Workshop on Lunar Base Methodology Development

Groundwork was laid for computer models which will assist in the design of a manned lunar base. The models, herein described, will provide the following functions for the successful conclusion of that task: strategic planning; sensitivity analyses; impact analyses; and documentation. Topics addressed include: upper level model description; interrelationship matrix; user community; model features; model descriptions; system implementation; model management; and plans for future action.

Nozette, Stewart↗
NASA NTRS1995

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 1: Theoretical developments and applications

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present first part of the report, we focus on the theoretical developments, and discussions of the results of numerical-performance studies using the integration schemes for GVIPS and NAV models.

Saleeb, Atef F.↗
NASA NTRS1995

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 2: Algorithmic developments and implementation

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present second part of the report, we focus on the specific details of the numerical schemes, and associated computer algorithms, for the finite-element implementation of GVIPS and NAV models.

Li, Wei↗
NASA NTRS1999

Development of GENOA Progressive Failure Parallel Processing Software Systems

A capability consisting of software development and experimental techniques has been developed and is described. The capability is integrated into GENOA-PFA to model polymer matrix composite (PMC) structures. The capability considers the physics and mechanics of composite materials and structure by integration of a hierarchical multilevel macro-scale (lamina, laminate, and structure) and micro scale (fiber, matrix, and interface) simulation analyses. The modeling involves (1) ply layering methodology utilizing FEM elements with through-the-thickness representation, (2) simulation of effects of material defects and conditions (e.g., voids, fiber waviness, and residual stress) on global static and cyclic fatigue strengths, (3) including material nonlinearities (by updating properties periodically) and geometrical nonlinearities (by Lagrangian updating), (4) simulating crack initiation. and growth to failure under static, cyclic, creep, and impact loads. (5) progressive fracture analysis to determine durability and damage tolerance. (6) identifying the percent contribution of various possible composite failure modes involved in critical damage events. and (7) determining sensitivities of failure modes to design parameters (e.g., fiber volume fraction, ply thickness, fiber orientation. and adhesive-bond thickness). GENOA-PFA progressive failure analysis is now ready for use to investigate the effects on structural responses to PMC material degradation from damage induced by static, cyclic (fatigue). creep, and impact loading in 2D/3D PMC structures subjected to hygrothermal environments. Its use will significantly facilitate targeting design parameter changes that will be most effective in reducing the probability of a given failure mode occurring.

Abdi, Frank↗
NASA NTRS1988

Unified constitutive model development for metal matrix composites at high temperature

Structural alloys used in high temperature applications exhibit complex thermomechanical behavior that is time dependent and hereditary. Recent attention is being focused on metal matrix composite materials for high temperature applications where they exhibit all the complexities of conventional alloys and their strong anisotropy adds further complexities. Here, a proven constitutive model for isotropic materials in which the inelastic strain rate and internal state are expressible as gradients of a dissipation potential is taken to depend on invariants that reflect local transverse isotropy. Applications illustrate the capability of the theory of representing the time dependent, hereditary, anisotropic behavior typical of these materials at high temperature.

Robinson, D. N.↗
NASA NTRS2010

Comparison of Damage Models for Predicting the Non-Linear Response of Laminates Under Matrix Dominated Loading Conditions

Five models for matrix damage in fiber reinforced laminates are evaluated for matrix-dominated loading conditions under plane stress and are compared both qualitatively and quantitatively. The emphasis of this study is on a comparison of the response of embedded plies subjected to a homogeneous stress state. Three of the models are specifically designed for modeling the non-linear response due to distributed matrix cracking under homogeneous loading, and also account for non-linear (shear) behavior prior to the onset of cracking. The remaining two models are localized damage models intended for predicting local failure at stress concentrations. The modeling approaches of distributed vs. localized cracking as well as the different formulations of damage initiation and damage progression are compared and discussed.

Schuecker, Clara↗
DOE OSTI2022

Model predictive control for matrix converter operating in current control mode with load current estimation

A matrix converter system operating in current control mode is provided. The matrix converter includes a switching matrix coupled between a low voltage side and a high voltage side, wherein the matrix converter is coupled at its low voltage side to a generator for receiving an input power including an input current and transforming the input power into an output power at its high voltage side, wherein a load is coupled to the high voltage side. The control system includes a load observer and a model predictive controller (MPC). The load observer is configured to estimate a load current that flows to the load from the high voltage side of the matrix converter as a function of a switching state of the switching matrix, an output voltage output at the high voltage side of the matrix converter, and the input current. The MPC is configured to select the switching state of the switching matrix to meet one or more control objectives defined by minimization of a multi-objective function that tracks the output voltage and the input current using the estimated load current.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1992

A differential CDM model for fatigue of unidirectional metal matrix composites

A multiaxial, isothermal, continuum damage mechanics (CDM) model for fatigue of a unidirectional metal matrix composite volume element is presented. The model is phenomenological, stress based, and assumes a single scalar internal damage variable, the evolution of which is anisotropic. The development of the fatigue damage model, (i.e., evolutionary law) is based on the definition of an initially transversely isotropic fatigue limit surface, a static fracture surface, and a normalized stress amplitude function. The anisotropy of these surfaces and function, and therefore the model, is defined through physically meaningful invariants reflecting the local stress and material orientation. This transversely isotropic model is shown, when taken to it's isotropic limit, to directly simplify to a previously developed and validated isotropic fatigue continuum damage model. Results of a nondimensional parametric study illustrate (1) the flexibility of the present formulation in attempting to characterize a class of composite materials, and (2) the capability of the formulation in predicting anticipated qualitative trends in the fatigue behavior of unidirectional metal matrix composites. Also, specific material parameters representing an initial characterization of the composite system SiC/Ti 15-3 and the matrix material (Ti 15-3) are reported.

Arnold, S. M.↗
NASA NTRS2013

Improved Subcell Model for the Prediction of Braided Composite Response

In this work, the modeling of triaxially braided composites was explored through a semi-analytical discretization. Four unique subcells, each approximated by a "mosaic" stacking of unidirectional composite plies, were modeled through the use of layered-shell elements within the explicit finite element code LS-DYNA. Two subcell discretizations were investigated: a model explicitly capturing pure matrix regions, and a novel model which absorbed pure matrix pockets into neighboring tow plies. The in-plane stiffness properties of both models, computed using bottom-up micromechanics, correlated well to experimental data. The absorbed matrix model, however, was found to best capture out-of- plane flexural properties by comparing numerical simulations of the out-of-plane displacements from single-ply tension tests to experimental full field data. This strong correlation of out-of-plane characteristics supports the current modeling approach as a viable candidate for future work involving impact simulations.

Cater, Christopher R.↗
DOE OSTI2023

An ICME Modeling Framework for Metal Matrix Composites Focusing on Ultrahigh Temperature Matrix Material and Tungsten Carbide Reinforcement Particulate (Final CRADA Report)

Ultra High Temperature Metal Matrix Composites (UHT-MMCs) are of interest to the due to their potential defense and aerospace applications. These materials consist of a metal matrix reinforced with a stiff ceramic or cermet. Based on a preliminary literature survey, this project will focused on a pure titanium (Ti) matrix reinforced with tungsten carbide-cobalt (WC-Co). UHT-MMC technology is still at a nascent stage of development and requires extensive experimental work to optimize the structure and resulting material properties of a particular UHT-MMC material system. An accurate and complete high-performance computing (HPC) model of the Ti/WC-Co material system could greatly accelerate the development and deployment Ti/WC-Co materials by connecting processing parameters to material properties. This type of model could then be used in an integrated computational materials engineering approach to optimize the material structure to achieve a targeted level of performance.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Phase-Field Modeling of Mechanical Damages in Ceramic Matrix Composites

Developed a phase-field model for mechanical damages in CMCs, which incorporates the CMC microstructures, matrix cracking, fiber breakage, and interfacial sliding. Two types of mechanisms, fiber bridging and fiber pull-out, are considered. The obtained simulation results agree with experimental observations and an analytical solution. Simulation results suggest that the performance of CMCs would be enhanced with thicker fibers, longer fibers, and higher fiber density. Opposite trends of interfacial sliding resistances are suggested for the two types of situations. In reality, a mixture of the two situations may exist, and then an intermediate interfacial sliding resistance may be optimal.

Xue, Fei↗