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.

SourceResearch connectionFrom yearThrough year

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

At least 793 records · Page 44

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Model reduction for discrete bilinear systems

A model reduction method for discrete bilinear systems is developed which matches q sets of Volterra and covariance parameters. These parameters are shown to represent both deterministic and stochastic attributes of the discrete bilinear system. A reduced order model which matches these q sets of parameters is defined to be a q-Volterra covariance equivalent realization (q-Volterra COVER). An algorithm is presented which constructs a class of q-Volterra COVERs parameterized by solutions to a Hermitian, quadratic, matrix equation. The algorithm is applied to a bilinear model of a robot manipulator.

King, A. M.↗

NASA NTRS1988

A micromechanics model of the stiffness and strength of laminates with fiber waviness

A mathematical model based on the Euler-Bernoulli beam theory is proposed for predicting the effective Young's moduli of piecewise isotropic composite laminates with local ply curvatures in the main load-carrying layers. Strains in corrugated layers, In-Phase layers, Out-of-Phase layers are predicted for various geometries and material configurations by assuming matrix layers as elastic foundations of different spring constants. The effective Young's moduli measured from corrugated aluminum specimens and aluminum/epoxy specimens with In-Phase and Out-of-Phase wavy patterns coincide very well with the model predictions. Moire fringe analysis of an In-Phase specimen and an Out-of-Phase specimen are also presented confirming the main assumption of the model related to the elastic constraint due to the matrix layers. The present model is also compared with the experimental results and other models, including the micro-buckling models, published in the literature. The results of the present study show that even a very small scale local ply curvature produces a noticeable effect on the mechanical constitutive behavior of a laminated composite.

Harris, Charles E.↗

NASA NTRS1990

A deformation-formulated micromechanics model of the effective Young's modulus and strength of laminated composites containing local ply curvature

A mathematical model based on the Euler-Bermoulli beam theory is proposed for predicting the effective Young's moduli of piecewise isotropic composite laminates with local ply curvatures in the main load-carrying layers. Strains in corrugated layers, in-phase layers, and out-of-phase layers are predicted for various geometries and material configurations by assuming matrix layers as elastic foundations of different spring constants. The effective Young's moduli measured from corrugated aluminum specimens and aluminum/epoxy specimens with in-phase and out-of-phase wavy patterns coincide very well with the model predictions. Moire fringe analysis of an in-phase specimen and an out-of-phase specimen are also presented, confirming the main assumption of the model related to the elastic constraint due to the matrix layers. The present model is also compared with the experimental results and other models, including the microbuckling models, published in the literature. The results of the present study show that even a very small-scale local ply curvature produces a noticeable effect on the mechanical constitutive behavior of a laminated composite.

Lee, Jong-Won↗

NASA NTRS1987

Material damping in aluminum and metal matrix composites

The material damping in beam-like specimens of aluminum and metal matrix composites was measured. A unique apparatus to determine damping by free decay while the specimens are in free fall in a vacuum was used. The specimens tested include 2024-T3 and 6061-T4 aluminum, and unidirectional graphite/metal matrix specimens with P55 and P100 fibers and 6061 Aluminum and AZ91C Magnesium as matrix materials. Tests were conducted to determine the dependence of damping on frequency and stress level. For the aluminum specimens, the material damping followed the Zener model at very low stress levels. Below the Zener relaxation frequency, a strong dependence of damping on stress was found for even moderate stress levels. Damping for the aluminum matrix materials was slightly above that predicted by the Zener model for a homogeneous bar of the matrix aluminum. For the magnesium matrix specimens, damping significantly above the Zener prediction for the homogeneous matrix material was observed.

Crawley, Edward F.↗

NASA NTRS1988

Micromechanical analysis of the time-dependent response of unidirectional composites

The micromechanical model of Aboudi (1982, 1984, and 1986) is used to predict the time-dependent creep behavior of viscoelastic-matrix composites with unidirectional elastic-fiber reinforcement. The derivation of the model and its numerical implementation are outlined, and predictions for T300/934 CFRP at 72 and 250 F are compared with experimental measurements in extensive tables and graphs. Good agreement is demonstrated, and the applicability of the present approach to composite laminates and hybrid materials such as ARALL is indicated.

Yancey, R. N.↗

NASA NTRSDate not supplied

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↗

NASA NTRS1997

Complexity of Kronecker Operations on Sparse Matrices with Applications to the Solution of Markov Models

We present a systematic discussion of algorithms to multiply a vector by a matrix expressed as the Kronecker product of sparse matrices, extending previous work in a unified notational framework. Then, we use our results to define new algorithms for the solution of large structured Markov models. In addition to a comprehensive overview of existing approaches, we give new results with respect to: (1) managing certain types of state-dependent behavior without incurring extra cost; (2) supporting both Jacobi-style and Gauss-Seidel-style methods by appropriate multiplication algorithms; (3) speeding up algorithms that consider probability vectors of size equal to the "actual" state space instead of the "potential" state space.

Buchholz, Peter↗

NASA NTRS1993

An analysis of fiber-matrix interface failure stresses for a range of ply stress states

A graphite/bismaleimide laminate was prepared without the usual fiber treatment and was tested over a wide range of stress states to measure its ply cracking strength. These tests were conducted using off-axis flexure specimens and produced fiber-matrix interface failure data over a correspondingly wide range of interface stress states. The absence of fiber treatment, weakened the fiber-matrix interfaces and allowed these tests to be conducted at load levels that did not yield the matrix. An elastic micromechanics computer code was used to calculate the fiber-matrix interface stresses at failure. Two different fiber-array models (square and diamond) were used in these calculations to analyze the effects of fiber arrangement as well as stress state on the critical interface stresses at failure. This study showed that both fiber-array models were needed to analyze interface stresses over the range of stress states. A linear equation provided a close fit to these critical stress combinations and, thereby, provided a fiber-matrix interface failure criterion. These results suggest that prediction procedures for laminate ply cracking can be based on micromechanics stress analyses and appropriate fiber-matrix interface failure criteria. However, typical structural laminates may require elastoplastic stress analysis procedures that account for matrix yielding, especially for shear-dominated ply stress states.

Crews, J. H.↗

NASA NTRS1991

A model for predicting high-temperature fatigue failure of a W/Cu composite

The material studied, a tungsten-fiber-reinforced, copper-matrix composite, is a candidate material for rocket nozzle liner applications. It was shown that at high temperatures, fatigue cracks initiate and propagate inside the copper matrix through a process of initiation, growth, and coalescence of grain boundary cavities. The ductile tungsten fibers neck and rupture locally after the surrounding matrix fails, and complete failure of the composite then ensues. A simple fatigue life prediction model is presented for the tungsten/copper composite system.

Verrilli, M. J.↗

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 NTRS1989

Thermoplastic matrix composites - Finite-element analysis of mode I and mode II failure specimens

A finite-element analysis was conducted to evaluate the stress distributions within mode I and mode II failure specimens, assuming both isotropic and orthotropic elastic material properties. The effects of anisotropy on both the magnitude and the location of the highest stress concentration at the vicinity of the crack tip are significant. The results from modeling realistic blunt crack tip geometry and resin-rich zones imply that local variations in the microstructure strongly influence the stress state near the crack tip and therefore the measured fracture properties. In addition, the features of a viscoelastic model for thermoplastic matrices are described. This model will be used in future investigations of matrix-dominated failure phenomena.

Bankert, Ray J.↗

NASA NTRS1991

Reliability analysis of ceramic matrix composite laminates

At a macroscopic level, a composite lamina may be considered as a homogeneous orthotropic solid whose directional strengths are random variables. Incorporation of these random variable strengths into failure models, either interactive or non-interactive, allows for the evaluation of the lamina reliability under a given stress state. Using a non-interactive criterion for demonstration purposes, laminate reliabilities are calculated assuming previously established load sharing rules for the redistribution of load as the failure of laminae occur. The matrix cracking predicted by ACK theory is modeled to allow a loss of stiffness in the fiber direction. The subsequent failure in the fiber direction is controlled by a modified bundle theory. Results using this modified bundle model are compared with previous models which did not permit separate consideration of matrix cracking, as well as to results obtained from experimental data.

Thomas, David J.↗

NASA NTRS2017

Incorporation of Failure Into an Orthotropic Three-Dimensional Model with Tabulated Input Suitable for Use in Composite 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. In the damage model, a semi-coupled approach is employed where the overall damage in a particular coordinate direction is assumed to be a multiplicative combination of the damage in that direction resulting from the applied loads in various coordinate directions. 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.

Impact↗

NASA NTRS1988

An approach to solving large reliability models

This paper describes a unified approach to the problem of solving large realistic reliability models. The methodology integrates behavioral decomposition, state trunction, and efficient sparse matrix-based numerical methods. The use of fault trees, together with ancillary information regarding dependencies to automatically generate the underlying Markov model state space is proposed. The effectiveness of this approach is illustrated by modeling a state-of-the-art flight control system and a multiprocessor system. Nonexponential distributions for times to failure of components are assumed in the latter example. The modeling tool used for most of this analysis is HARP (the Hybrid Automated Reliability Predictor).

Boyd, Mark A.↗

NASA NTRS2000

Application of a Reduced Order Kalman Filter to Initialize a Coupled Atmosphere-Ocean Model: Impact on the Prediction of El Nino

A reduced order Kalman Filter, based on a simplification of the Singular Evolutive Extended Kalman (SEEK) filter equations, is used to assimilate observed fields of the surface wind stress, sea surface temperature and sea level into the nonlinear coupled ocean-atmosphere model. The SEEK filter projects the Kalman Filter equations onto a subspace defined by the eigenvalue decomposition of the error forecast matrix, allowing its application to high dimensional systems. The Zebiak and Cane model couples a linear reduced gravity ocean model with a single vertical mode atmospheric model of Zebiak. The compatibility between the simplified physics of the model and each observed variable is studied separately and together. The results show the ability of the model to represent the simultaneous value of the wind stress, SST and sea level, when the fields are limited to the latitude band 10 deg S - 10 deg N. In this first application of the Kalman Filter to a coupled ocean-atmosphere prediction model, the sea level fields are assimilated in terms of the Kelvin and Rossby modes of the thermocline depth anomaly. An estimation of the error of these modes is derived from the projection of an estimation of the sea level error over such modes. This method gives a value of 12 for the error of the Kelvin amplitude, and 6 m of error for the Rossby component of the thermocline depth. The ability of the method to reconstruct the state of the equatorial Pacific and predict its time evolution is demonstrated. The method is shown to be quite robust for predictions I up to six months, and able to predict the onset of the 1997 warm event fifteen months before its occurrence.

Ballabrera-Poy, Joaquim↗

NASA NTRSDate not supplied

Unsteady Aerodynamic Modeling of Atmospheric Entry Vehicles in Subsonic and Incompressible flow: A Frequency Response Approach

The determination of time-varying lift force and pitch moment generated by a purely pitching Earth-entry capsule is investigated. Experiments were conducted in the 12-foot Low-Speed Tunnel at NASA Langley Research Center, testing a range of oscillation frequencies at a zero-mean angle of attack with a pitching amplitude of 10 degrees. In light of these measurements, a closed-form set of analytically derived equations for lift and moment was used to develop a semi-empirical formulation, incorporating empirically determined values from the experimental runs. The equations are grounded in potential flow theory, Theodorsen's classical theory of unsteady aerodynamics, and the Joukowski theorem of conformal mapping. The unsteady aerodynamics generated by the oscillating body are then modeled by constructing frequency response functions, with quasi-steady forces and moments serving as inputs and unsteady forces and moments as outputs. The experimentally determined gain and phase variations characterize the unsteady nature of the flow and the system's response and flow time-lag to input flow parameters for a blunt-body entry vehicle. The final semi-empirical model is validated with a set of parameters beyond the initial test matrix.

Entry Systems Modeling↗

NASA NTRS2019

Bayesian Inference and the Effects of Varying Uncertainty Models in Charring Ablator Calibration and Uncertainty Quantification Problems

The Mars Science Laboratory (MSL) vehicle utilized a heat shield constructed from NASA’s Phenolic-Impregnated Carbon Ablator (PICA) material to protect the main structure from the high enthalpy environment encountered during hypersonic atmospheric entry. During the vehicle’s descent through Martian atmosphere, multiple thermocouples embedded within the heat shield captured in-depth material temperature data that allow for studies to be conducted on current material response reconstruction tools. In the present work, material temperature data obtained from thermocouples within the MISP-4 plug (MEDLI (Mars Science Laboratory Entry, Descent, and Landing Instrument) Integrated Sensor Plug) are utilized in the calibration of Theoretical Ablative Composite for Open Testing (TACOT) model parameters in conjunction with NASA’s Porous material Analysis Toolbox (PATO) through Bayesian inference where uncertainty due to parametric, modeling, and experimental sources is simultaneously quantified. Prior to the study, a sensitivity analysis is performed through computation of the robust Sobol indices in an effort to study the relationship between input space and model response and to reduce the dimensionality of the statistical inverse problem. The Bayesian inference methodology necessitates an a-priori choice to be made for the uncertainty model for which numerous possibilities are available. Across most works, however, only basic additive or multiplicative models are utilized with pre-defined magnitudes of uncertainty based on a-priori knowledge or to-be-calibrated multipliers of static covariance matrix structures. The present effort explores the effects of informed uncertainty models, ones with temporal dependence that are simultaneously calibrated through Bayesian inference, on calibrated results for parameters that make up the uncertain input space.

Sensitivity Analysis↗

NASA NTRS1979

Analysis of a unidirectional composite containing broken fibers and matrix damage

An analytical solution is developed for the determination of the stresses and displacements in a unidirectional fiber-reinforced composite containing an arbitrary number of broken fibers as well as longitudinal yielding and splitting of the matrix. The solution is developed using a materials-modeling approach which is based on a shear-lag stress transfer mechanism. The equilibrium equation in the axial direction gives a pair of integral equations which are solved numerically. Excellent agreement is shown to exist between the solution and experimental results for notched unidirectional boron/aluminum laminates without splitting. For brittle matrix composites (i.e. epoxy) equally good results are indicated for both matrix yielding and splitting. For yielding without splitting the fracture strength depends on crack length while for large splitting it is crack length independent.

George, J. G.↗