Search NASA⌕ Search

SEARCH · Search NASA

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Modeling of crack bridging in a unidirectional metal matrix composite

The effective fatigue crack driving force and crack opening profiles were determined analytically for fatigue tested unidirectional composite specimens exhibiting fiber bridging. The crack closure pressure due to bridging was modeled using two approaches: the fiber pressure model and the shear lag model. For both closure models, the Bueckner weight function method and the finite element method were used to calculate crack opening displacements and the crack driving force. The predicted near crack tip opening profile agreed well with the experimentally measured profiles for single edge notch SCS-6/Ti-15-3 metal matrix composite specimens. The numerically determined effective crack driving force, Delta K(eff), was calculated using both models to correlate the measure crack growth rate in the composite. The calculated Delta K(eff) from both models accounted for the crack bridging by showing a good agreement between the measured fatigue crack growth rates of the bridged composite and that of unreinforced, unbridged titanium matrix alloy specimens.

Ghosn, Louis J.↗
NASA NTRS1999

Improved Stratospheric Temperature Retrievals for Climate Reanalysis

The Data Assimilation Office (DAO) is embarking on plans to generate a twenty year reanalysis data set of climatic atmospheric variables. One of the focus points will be in the evaluation of the dynamics of the stratosphere. The Stratospheric Sounding Unit (SSU), flown as part of the TIROS Operational Vertical Sounder (TOVS), is one of the primary stratospheric temperature sensors flown consistently throughout the reanalysis period. Seven unique sensors made the measurements over time, with individual instrument characteristics that need to be addressed. The stratospheric temperatures being assimilated across satellite platforms will profoundly impact the reanalysis dynamical fields. To attempt to quantify aspects of instrument and retrieval bias we are carefully collecting and analyzing all available information on the sensors, their instrument anomalies, forward model errors and retrieval biases. For the retrieval of stratospheric temperatures, we adapted the minimum variance approach of Jazwinski (1970) and Rodgers (1976) and applied it to the SSU soundings. In our algorithm, the state vector contains an initial guess of temperature from a model six hour forecast provided by the Goddard EOS Data Assimilation System (GEOS/DAS). This is combined with an a priori covariance matrix, a forward model parameterization, and specifications of instrument noise characteristics. A quasi-Newtonian iteration is used to obtain convergence of the retrieved state to the measurement vector. This algorithm also enables us to analyze and address the systematic errors associated with the unique characteristics of the cell pressures on the individual SSU instruments and the resolving power of the instruments to vertical gradients in the stratosphere. The preliminary results of the improved retrievals and their assimilation as well as baseline calculations of bias and rms error between the NESDIS operational product and col-located ground measurements will be presented.

Rokke, L.↗
NASA NTRS2004

Theoretical Investigation of OCN(-) Charge Transfer Complexes in Condensed Phase Media: Spectroscopic Properties in Amorphous Ice

Density functional theory (DFT) calculations of cyanate (OCN(-)) charge-transfer complexes were performed to model the "XCN" feature observed in interstellar icy grain mantles. OCN(-) charge-transfer complexes were formed from precursor combinations of HNCO or HOCN with either NH3 or H2O. Three different solvation strategies for realistically modeling the ice matrix environment were explored, including (1) continuum solvation, (2) pure DFT cluster calculations, and (3) an ONIOM DFT/PM3 cluster calculation. The model complexes were evaluated by their ability to reproduce seven spectroscopic measurements associated with XCN: the band origin of the OCN(-) asymmetric stretching mode, shifts in that frequency due to isotopic substitutions of C, N, O, and H, plus two weak features. The continuum solvent field method produced results consistent with some of the experimental data but failed to account for other behavior due to its limited capacity to describe molecular interactions with solvent. DFT cluster calculations successfully reproduced the available spectroscopic measurements very well. In particular, the deuterium shift showed excellent agreement in complexes where OCN(-) was fully solvated. Detailed studies of representative complexes including from two to twelve water molecules allowed the exploration of various possible solvation structures and provided insights into solvation trends. Moreover, complexes arising from cyanic or isocyanic acid in pure water suggested an alternative mechanism for the formation of OCN(-) charge-transfer complexes without the need for a strong base such as NH3 to be present. An extended ONIOM (B3LYP/PM3) cluster calculation was also performed to assess the impact of a more realistic environment on HNCO dissociation in pure water.

Park, Jin-Young↗
NASA NTRS1994

A coupled/uncoupled deformation and fatigue damage algorithm utilizing the finite element method

A fatigue damage computational algorithm utilizing a multiaxial, isothermal, continuum based fatigue damage model for unidirectional metal matrix composites has been implemented into the commercial finite element code MARC using MARC user subroutines. Damage is introduced into the finite element solution through the concept of effective stress which fully couples the fatigue damage calculations with the finite element deformation solution. An axisymmetric stress analysis was performed on a circumferentially reinforced ring, wherein both the matrix cladding and the composite core were assumed to behave elastic-perfectly plastic. The composite core behavior was represented using Hill's anisotropic continuum based plasticity model, and similarly, the matrix cladding was represented by an isotropic plasticity model. Results are presented in the form of S-N curves and damage distribution plots.

Wilt, Thomas E.↗
NASA NTRS1991

Modeling of crack bridging in a unidirectional metal matrix composite

The effective fatigue crack driving force and crack opening profiles were determined analytically for fatigue tested unidirectional composite specimens exhibiting fiber bridging. The crack closure pressure due to bridging was modeled using two approaches; the fiber pressure model and the shear lag model. For both closure models, the Bueckner weight function method and the finite element method were used to calculate crack opening displacements and the crack driving force. The predicted near crack tip opening profile agreed well with the experimentally measured profiles for single edge notch SCS-6/Ti-15-3 metal matrix composite specimens. The numerically determined effective crack driving force, Delta K(sup eff), was calculated using both models to correlate the measure crack growth rate in the composite. The calculated Delta K(sup eff) from both models accounted for the crack bridging by showing a good agreement between the measured fatigue crack growth rates of the bridged composite and that of unreinforced, unbridged titanium matrix alloy specimens.

Ghosn, Louis J.↗
NASA NTRS1972

Structural model optimization using statistical evaluation

The results of research in applying statistical methods to the problem of structural dynamic system identification are presented. The study is in three parts: a review of previous approaches by other researchers, a development of various linear estimators which might find application, and the design and development of a computer program which uses a Bayesian estimator. The method is tried on two models and is successful where the predicted stiffness matrix is a proper model, e.g., a bending beam is represented by a bending model. Difficulties are encountered when the model concept varies. There is also evidence that nonlinearity must be handled properly to speed the convergence.

Collins, J. D.↗
NASA NTRS1973

Decomposition-aggregation stability analysis

This report presents the development and description of the decomposition aggregation approach to stability investigations of high dimension mathematical models of dynamic systems. The high dimension vector differential equation describing a large dynamic system is decomposed into a number of lower dimension vector differential equations which represent interconnected subsystems. Then a method is described by which the stability properties of each subsystem are aggregated into a single vector Liapunov function, representing the aggregate system model, consisting of subsystem Liapunov functions as components. A linear vector differential inequality is then formed in terms of the vector Liapunov function. The matrix of the model, which reflects the stability properties of the subsystems and the nature of their interconnections, is analyzed to conclude over-all system stability characteristics. The technique is applied in detail to investigate the stability characteristics of a dynamic model of a hypothetical spinning Skylab.

Siljak, D. D.↗
NASA NTRS1988

Nonlinear modeling, identification and estimation of slew induced structural deformation

This proposal discusses a new nonlinear, nonparametric method for off-line modeling and on-line estimation of the deformation of a flexible structure undergoing rapid retargeting maneuvers. In these circumstances, the structural stiffness and damping coefficients depend on the angular acceleration omega(dot), the angular rate omega, and the square of the angular rate omega. In the single axis case, the excitation of the structure is represented by the vector u(exp T) = (omega(dot), omega(exp 2), 2(omega)), to which the structural dynamics responds as a 'bilinear' (i.e., parametrically excited) system. A similar technique for multiaxial rotations yields a bilinear model with respect to matrix valued excitations. Three methods of estimation and modeling are described in this proposal to achieve deformation state determination: (1) a method based on a feedback linearized procedure which gives an estimate of the state by means of observers installed in the deformable body; (2) off-line modeling of the deformation state of the structure by means of topological interpolators; and (3) an on-line structural state estimation method based on a combination of the two previous techniques.

Dwyer, Thomas A., III↗
NASA NTRS2001

Ultrasonic Assessment of Impact-Induced Damage and Microcracking in Polymer Matrix Composites

The main objective of this NASA FAR project is to conduct ultrasonic assessment of impact-induced damage and microcracking in polymer matrix composites at various temperatures. It is believed that the proposed study of impact damage assessment on polymer matrix composites will benefit several NASA missions and current interests, such as ballistic impact testing of composite fan containment and high strain rate deformation modeling of polymer matrix composites. Impact-induced damage mechanisms in GLARE and ARALL fiber-metal laminates subject to instrumented drop-weight impacts at various temperatures were studied. GLARE and ARALL are hybrid composites made of alternating layers of aluminum and glass (for GLARE) and aramid- (for ARALL) fiber-reinforced epoxy. Damage in pure aluminum panels impacted by foreign objects was mainly characterized by large plastic deformation surrounding a deep penetration dent. On the other hand, plastic deformation in fiber-metal laminates was often not as severe although the penetration dent was still produced. The more stiff fiber-reinforced epoxy layers provided better bending rigidity; thus, enhancing impact damage tolerance. Severe cracking, however, occurred due to the use of these more brittle fiber-reinforced epoxy layers. Fracture patterns, e.g., crack length and delamination size, were greatly affected by the lay-up configuration rather than by the number of layers, which implies that thickness effect was not significant for the panels tested in this study. Immersion ultrasound techniques were then used to assess damages generated by instrumented drop-weight impacts onto these fiber-metal laminate panels as well as 6061-T6 aluminum/cast acrylic sandwich plates adhered by epoxy. Depending on several parameters, such as impact velocity, mass, temperature, laminate configuration, sandwich construction, etc., various types of impact damage were observed, including plastic deformation, radiating cracks emanating from the impact site, ring cracks surrounding the impact site, partial and full delamination, and combinations of these damages.

Liaw, Benjamin↗
NASA NTRS2006

An Initial Non-Equilibrium Porous-Media Model for CFD Simulation of Stirling Regenerators

The objective of this paper is to define empirical parameters for an initial thermal non-equilibrium porous-media model for use in Computational Fluid Dynamics (CFD) codes for simulation of Stirling regenerators. The two codes currently used at Glenn Research Center for Stirling modeling are Fluent and CFD-ACE. The codes porous-media models are equilibrium models, which assume solid matrix and fluid are in thermal equilibrium. This is believed to be a poor assumption for Stirling regenerators; Stirling 1-D regenerator models, used in Stirling design, use non-equilibrium regenerator models and suggest regenerator matrix and gas average temperatures can differ by several degrees at a given axial location and time during the cycle. Experimentally based information was used to define: hydrodynamic dispersion, permeability, inertial coefficient, fluid effective thermal conductivity, and fluid-solid heat transfer coefficient. Solid effective thermal conductivity was also estimated. Determination of model parameters was based on planned use in a CFD model of Infinia's Stirling Technology Demonstration Converter (TDC), which uses a random-fiber regenerator matrix. Emphasis is on use of available data to define empirical parameters needed in a thermal non-equilibrium porous media model for Stirling regenerator simulation. Such a model has not yet been implemented by the authors or their associates.

Tew, Roy C.↗
NASA NTRS2004

Modeling the Nonlinear, Strain Rate Dependent Deformation of Woven Ceramic Matrix Composites With Hydrostatic Stress Effects Included

An analysis method based on a deformation (as opposed to damage) approach has been developed to model the strain rate dependent, nonlinear deformation of woven ceramic matrix composites with a plain weave fiber architecture. In the developed model, the differences in the tension and compression response have also been considered. State variable based viscoplastic equations originally developed for metals have been modified to analyze the ceramic matrix composites. To account for the tension/compression asymmetry in the material, the effective stress and effective inelastic strain definitions have been modified. The equations have also been modified to account for the fact that in an orthotropic composite the in-plane shear stiffness is independent of the stiffness in the normal directions. The developed equations have been implemented into a commercially available transient dynamic finite element code, LS-DYNA, through the use of user defined subroutines (UMATs). The tensile, compressive, and shear deformation of a representative plain weave woven ceramic matrix composite are computed and compared to experimental results. The computed values correlate well to the experimental data, demonstrating the ability of the model to accurately compute the deformation response of woven ceramic matrix composites.

Goldberg, Robert K.↗
NASA NTRS2000

Ultrasonic Assessment of Impact-Induced Damage and Microcracking in Polymer Matrix Composites

The main objective of this NASA FAR project is to conduct ultrasonic assessment of impact-induced damage and microcracking in polymer matrix composites at various temperatures. It is believed that the proposed study of impact damage assessment on polymer matrix composites will benefit several NASA's missions and current interests, such as ballistic impact testing of composite fan containment and high strain rate deformation modeling of polymer matrix composites. Currently, impact-induced delamination and fracture in 6061-T6 aluminum/cast acrylic sandwich plates adhered by epoxy were generated in an instrumented drop-weight impact machine. Although only a small dent was produced on the aluminum side when a hemispherical penetrator tup was dropped onto it from a couple of inches, a large ring of delamination at the interface was observed. The delamination damage was often accompanied by severe shattering in the acrylic substratum. Damage patterns in the acrylic layer include radial and ring cracks and, together with delamination at the interface, may cause peeling-off of acrylic material from the sandwich plate. Theory of stress-wave propagation can be used to explain these damage patterns. The impact tests were conducted at various temperatures. The results also show clearly that temperature effect is very important in impact damage. For pure cast acrylic nil-ductile transition (NDT) occurs between 185-195 F Excessive impact energy was dissipated into fracture energy when tested at temperature below this range or through plastic deformation when tested at temperature above the NDT temperature. Results from this study will be used as baseline data for studying fiber-metal laminates, such as GLARE and ARALL for advanced aeronautical and astronautical applications.

Liaw, Benjamin↗
NASA NTRS2004

Bayesian Approach to the Joint Inversion of Gravity and Magnetic Data, with Application to the Ismenius Area of Mars

This viewgraph presentation reviews a Bayesian approach to the inversion of gravity and magnetic data with specific application to the Ismenius Area of Mars. Many inverse problems encountered in geophysics and planetary science are well known to be non-unique (i.e. inversion of gravity the density structure of a body). In hopes of reducing the non-uniqueness of solutions, there has been interest in the joint analysis of data. An example is the joint inversion of gravity and magnetic data, with the assumption that the same physical anomalies generate both the observed magnetic and gravitational anomalies. In this talk, we formulate the joint analysis of different types of data in a Bayesian framework and apply the formalism to the inference of the density and remanent magnetization structure for a local region in the Ismenius area of Mars. The Bayesian approach allows prior information or constraints in the solutions to be incorporated in the inversion, with the "best" solutions those whose forward predictions most closely match the data while remaining consistent with assumed constraints. The application of this framework to the inversion of gravity and magnetic data on Mars reveals two typical challenges - the forward predictions of the data have a linear dependence on some of the quantities of interest, and non-linear dependence on others (termed the "linear" and "non-linear" variables, respectively). For observations with Gaussian noise, a Bayesian approach to inversion for "linear" variables reduces to a linear filtering problem, with an explicitly computable "error" matrix. However, for models whose forward predictions have non-linear dependencies, inference is no longer given by such a simple linear problem, and moreover, the uncertainty in the solution is no longer completely specified by a computable "error matrix". It is therefore important to develop methods for sampling from the full Bayesian posterior to provide a complete and statistically consistent picture of model uncertainty, and what has been learned from observations. We will discuss advanced numerical techniques, including Monte Carlo Markov

data analysis↗
NASA NTRS2004

A Method for Incorporating Changing Structural Characteristics Due to Propellant Mass Usage in a Launch Vehicle Ascent Simulation

Launch vehicles consume large quantities of propellant quickly, causing the mass properties and structural dynamics of the vehicle to change dramatically. Currently, structural load assessments account for this change with a large collection of structural models representing various propellant fill levels. This creates a large database of models complicating the delivery of reduced models and requiring extensive work for model changes. Presented here is a method to account for these mass changes in a more efficient manner. The method allows for the subtraction of propellant mass as the propellant is used in the simulation. This subtraction is done in the modal domain of the vehicle generalized model. Additional computation required is primarily for constructing the used propellant mass matrix from an initial propellant model and further matrix multiplications and subtractions. An additional eigenvalue solution is required to uncouple the new equations of motion; however, this is a much simplier calculation starting from a system that is already substantially uncoupled. The method was successfully tested in a simulation of Saturn V loads. Results from the method are compared to results from separate structural models for several propellant levels, showing excellent agreement. Further development to encompass more complicated propellant models, including slosh dynamics, is possible.

McGhee, D. S.↗
NASA NTRS1985

Short-wavelength buckling and shear failures for compression-loaded composite laminates

The short-wavelength buckling (or the microbuckling) and the interlaminar and inplane shear failures of multi-directional composite laminates loaded in uniaxial compression are investigated. A laminate model is presented that idealizes each lamina. The fibers in the lamina are modeled as a plate, and the matrix in the lamina is modeled as an elastic foundation. The out-of-plane w displacement for each plate is expressed as a trigonometric series in the half-wavelength of the mode shape for laminate short-wavelength buckling. Nonlinear strain-displacement relations are used. The model is applied to symmetric laminates having linear material behavior. The laminates are loaded in uniform end shortening and are simply supported. A linear analysis is used to determine the laminate stress, strain, and mode shape when short-wavelength buckling occurs. The equations for the laminate compressive stress at short-wavelength buckling are dominated by matrix contributions.

Shuart, M. J.↗
NASA NTRS1993

Evolution of damage and plasticity in titanium-based, fiber-reinforced composites

The inelastic deformation mechanisms were evaluated for a model titanium-based, fiber-reinforced composite: a beta titanium alloy (Ti-15V-3Al-3Cr-3Sn) reinforced with SiC (SCS-6) fibers. The primary emphasis of this article is to illustrate the sequence in which damage and plasticity evolved for this system. The mechanical responses and the results of detailed microstructural evaluations for the 0(8), 90(8), and +/- 45(2s) line oriented laminates are provided. It is shown that the characteristics of the reaction zone around the fiber play a very important role in the way damage and plasticity evolve, particularly in the microyield regime of deformation, and must be included in any realistic constitutive model. Fiber-matrix debonding was a major damage mode for the off-axis systems. The tension test results are also compared with the predictions of a few constitutive models.

Majumdar, B. S.↗
NASA NTRS1987

Identification of large space structures: A state-of-practice report

An outline of this work is presented. It begins with a schematic flow diagram and a logical flow diagram of the identification process for large space structures (LSS). Next, the task is defined by a structure model definition. A matrix polynomial formulation with a node displacement equation and a state variable formulation with node displacement and velocities are outlined. Further outlined is the identification of LSS on orbit; modeling errors and uncertainties; verification and validation of model; and noise, computations, and data collection.

Source record↗
NASA NTRS2018

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

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

layer shifting↗