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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Influence of Strain-Gage Balance Characteristics on the Precision Error of the Drag Coefficient

Improvements and results of a new method are presented that computes a pre-test estimate of the precision error of the drag coefficient of a wind tunnel model. The error estimate is defined as the part of the drag coefficient's precision error that is primarily associated with the precision error of the angle of attack measurement and physical characteristics of the chosen strain-gage balance. The method indirectly describes the precision error of the angle of attack measurement by using an assumed balance gage output variation of one microV/V. The physical characteristics of the balance, on the other hand, are described by partial derivatives of the axial and normal forces with respect to the strain-gage outputs. These derivatives can directly be obtained from the data reduction matrix of the balance. The precision error estimate itself is calculated by applying a simple explicit equation that uses the model reference area, the dynamic pressure, the angle of attack, the coefficients of the linear terms of the data reduction matrix, and the electrical output variation of one microvolt per volt as input. Precision errors at constant angle of attack may be visualized as contour plots by plotting them, for example, versus the Mach number and the total pressure. Characteristics of NASA's MC60E balance are used in combination with the reference area of a generic wind tunnel model in order to demonstrate that error estimates are independent of both the balance load format and the units chosen for the description of balance loads, model reference area, and the dynamic pressure. Finally, experimental data from a wind tunnel test of the Ames Check Standard Model in the NASA Ames 11-foot Transonic Wind Tunnel illustrates the application of the method to real-world test data.

precision error↗
NASA NTRSDate not supplied

Multiscale Failure Analysis of a 3D Woven Unit Cell Containing Defects

Multiscale failure simulations have been performed for a Three-dimensional woven composite unit cell considering five, or more, length scales spanning the woven composite mesoscale to the sub-microscale voids. The multiscale recursive micromechanics approach, which enables recursive integration of general micromechanics theories over an arbitrary number of length scales, has been employed within the NASA Multiscale Analysis Tool. The multiscale model uses both the generalized method of cells and Mori-Tanaka micromechanics theories, and considers failure in the constituent materials using a simple damage model. Baseline results, containing distributed voids, are compared to uniaxial experimental data for an AS4 carbon fiber/ RTM6 epoxy matrix 3D orthogonal woven composite with good agreement in terms of global stiffness and global failure stress. The simulations demonstrate that the 3D woven composite exhibits damage tolerance through sustaining increasing axial load far beyond the first initiation of damage. The multiscale model is used to examine the nonlinear response of the material to other loading conditions. Case studies, motivated by X-ray computed tomography data, are presented on the effects of manufacturing induced voids and cracks.

3D woven↗
NASA NTRS2024

Summary of Model-based Attitude Control of LUVOIR in Modular Dynamic Analysis (MDA) Simulink Environment

This document overviews the work I completed in the second half of my summer 2018 internship experience in Code 591 at NASA Goddard Space Flight Center. Please see the former memorandum for additional context on the project. The take away point is that LUVOIR A has been modeled as three rigid bodies linked by 1-DOF rotary joints. An LQR attitude controller was designed for precision pointing of the spacecraft with the design requirement that steady state oscillations have amplitude less than a miliarcsecond. The performance of the algorithm was tested in the Modular Dynamic Analysis Simulink library developed by J. Roger Chen at NASA Goddard. While the initial simulations were of rigid bodies, subsequent simulations included flexible body modes on all three bodies of the model. The simulated structural deformations initially destabilized the LQR controller thus requiring a redesign of the K gain matrix. The final simulation demonstrated a stabilizing feedback gain with miliarcsecond precision within 400 minutes. This run used 20 modes on the spacecraft, 9 on the bus, and 100 on the payload. Future recommended work includes the development of a Vibration Isolation and Precision Pointing System (VIPPS) Simulink model. Armed with this model, the system should be modeled as four bodies whose associated flexible files must be generated from the spacecraft through Gimbal 1, Gimbal 1 through Gimbal 2, Gimbal 2 through the VIPPS interface, and the VIPPS interface through the payload.

William Bentz↗
NASA NTRS1994

Orientation of human optokinetic nystagmus to gravity: a model-based approach

Optokinetic nystagmus (OKN) was induced by having subjects watch a moving display in a binocular, head-fixed apparatus. The display was composed of 3.3 degrees stripes moving at 35 degrees/s for 45 s. It subtended 88 degrees horizontally by 72 degrees vertically of the central visual field and could be oriented to rotate about axes that were upright or tilted 45 degrees or 90 degrees. The head was held upright or was tilted 45 degrees left or right on the body during stimulation. Head-horizontal (yaw axis) and head-vertical (pitch axis) components of OKN were recorded with electro-oculography (EOG). Slow phase velocity vectors were determined and compared with the axis of stimulation and the spatial vertical (gravity axis). With the head upright, the axis of eye rotation during yaw axis OKN was coincident with the stimulus axis and the spatial vertical. With the head tilted, a significant vertical component of eye velocity appeared during yaw axis stimulation. As a result the axis of eye rotation shifted from the stimulus axis toward the spatial vertical. Vertical components developed within 1-2 s of stimulus onset and persisted until the end of stimulation. In the six subjects there was a mean shift of the axis of eye rotation during yaw axis stimulation of approximately 18 degrees with the head tilted 45 degrees on the body. Oblique optokinetic stimulation with the head upright was associated with a mean shift of the axis of eye rotation toward the spatial vertical of 9.2 degrees. When the head was tilted and the same oblique stimulation was given, the axis of eye rotation rotated to the other side of the spatial vertical by 5.4 degrees. This counterrotation of the axis of eye rotation is similar to the "Muller (E) effect," in which the perception of the upright is counterrotated to the opposite side of the spatial vertical when subjects are tilted in darkness. The data were simulated by a model of OKN with a "direct" and "indirect" pathway. It was assumed that the direct visual pathway is oriented in a body, not a spatial frame of reference. Despite the short optokinetic after-nystagmus time constants, strong horizontal to vertical cross-coupling could be produced if the horizontal and vertical time constants were in proper ratio and there were no suppression of nystagmus in directions orthogonal to the stimulus direction. The model demonstrates that the spatial orientation of OKN can be achieved by restructuring the system matrix of velocity storage. We conclude that an important function of velocity storage is to orient slow-phase velocity toward the spatial vertical during movement in a terrestrial environment.

Non-NASA Center↗
NASA NTRS1981

Creep and creep rupture of laminated graphite/epoxy composites

An incremental numerical procedure based on lamination theory is developed to predict creep and creep rupture of general laminates. Existing unidirectional creep compliance and delayed failure data is used to develop analytical models for lamina response. The compliance model is based on a procedure proposed by Findley which incorporates the power law for creep into a nonlinear constitutive relationship. The matrix octahedral shear stress is assumed to control the stress interaction effect. A modified superposition principle is used to account for the varying stress level effect on the creep strain. The lamina failure model is based on a modification of the Tsai-Hill theory which includes the time dependent creep rupture strength. A linear cumulative damage law is used to monitor the remaining lifetime in each ply.

Dillard, D. A.↗
NASA NTRS2014

Modeling Quasi-Static and Fatigue-Driven Delamination Migration

An approach was proposed and assessed for the high-fidelity modeling of progressive damage and failure in composite materials. It combines the Floating Node Method (FNM) and the Virtual Crack Closure Technique (VCCT) to represent multiple interacting failure mechanisms in a mesh-independent fashion. Delamination, matrix cracking, and migration were captured failure and migration criteria based on fracture mechanics. Quasi-static and fatigue loading were modeled within the same overall framework. The methodology proposed was illustrated by simulating the delamination migration test, showing good agreement with the available experimental data.

De Carvalho, N. V.↗
NASA NTRS1989

Structural robustness with suboptimal responses for linear state space model

A relationship between the closed-loop eigenvalues and the amount of perturbations in the open-loop matrix is addressed in the context of performance robustness. If the allowable perturbation ranges of elements of the open-loop matrix A and the desired tolerance of the closed-loop eigenvalues are given such that max(j) of the absolute value of Delta-lambda(j) (A+BF) should be less than some prescribed value, what is a state feedback controller F which satisfies the closed-loop eigenvalue perturbation-tolerance requirement for a class of given perturbation in A? The paper gives an algorithm to design such a controller. Numerical examples are included for illustration.

Keel, L. H.↗
NASA NTRS1991

Turbulence modeling in supersonic combusting flows

To support the National Aerospace Plane project, the RPLUS3D CFD code has been developed at NASA Lewis. The code has the ability to solve three-dimensional flowfields with finite rate combustion of hydrogen and air. The combustion processes of the hydrogen-air system are simulated by an 18-reaction path, 8-species chemical kinetic mechanism. The code uses a Lower-Upper (LU) decomposition numerical algorithm as its basis, making it a very efficient and robust code. Except for the Jacobian matrix for the implicit chemistry source terms, there is no inversion of a matrix even though it uses a fully implicit numerical algorithm. A k-epsilon (two equation) turbulence model is incorporated into the RPLUS3D code.

Chitsomboon, Tawit↗
NASA NTRS1988

Extension and validation of a method for locating damaged members in large space trusses

The damage location approach employs the control system capabilities for the structure to test the structure and measure the dynamic response. The measurements are then used in a system identification algorithm to produce a model of the damaged structure. The model is compared to one for the undamaged structure to find regions of reduced stiffness which indicate the location of damage. Kabe's 3,4 stiffness matrix adjustment method was the central identification algorithm. The strength of his method is that, with minimal data, it preserves the representation of the physical connectivity of the structure in the resulting model of the damaged truss. However, extensive storage and computational effort were required as a result. Extension of the damage location method to overcome these problems is the first part of the current work. The central system identification algorithm is replaced with the MSMT method of stiffness matrix adjustment which was previously derived by generalizing an optimal-update secant method form quasi-Newton approaches for nonlinear optimization. Validation of the extended damage location method is the second goal.

Smith, Suzanne Weaver↗
NASA NTRS1989

Thermal histories of the samples of two KOSI comet nucleus simulation experiments

Temperatures recorded during two KOSI comet nucleus simulation experiments strongly suggest that heat transport by vapor flow into the interior of the sample is very important. Two comet nucleus simulation experiments have been done by the KOSI team in a big space simulator. The thermal evolution of the sample during insolation and the results of simplified thermal evolution calculations are discussed. The observed thermal histories cannot be explained by a simple model with heat transferred by heat conduction at a constant conductivity, so a coupled heat and mass transfer problem was considered. The porous ice matrix was assumed to have a constant thermal conductivity and to be in thermal equilibrium with vapor in the pores, the internal pressure being the vapor pressure. The vapor was modelled as an ideal gas because, at the temperatures relevant to the problem, the mean free path length of the vapor molecules is large in comparison with the pore dimensions. The heat capacity at constant volume per unit mass of the two phase mixture was also assumed constant. The vapor was allowed to flow and transfer heat in response to an internal pressure gradient.

Spohn, T.↗
NASA NTRS1977

Crystallization kinetics, viscous flow, and thermal history of lunar breccia 67975

The maximum crystal growth rate for the 67975 matrix composition, about 0.00073 cm/min, was found to occur at about 990 C (undercooling 220 C). The Shaw (1972) model and, to a slightly lesser extent, the Bottinga and Weill (1972) model indicate the variation of viscosity with temperature. Time-temperature-transformation and continuous cooling curves were constructed from the measured crystal growth rates and viscosities and used to estimate the cooling rate required to form bodies with the degree of crystallinity observed in the matrix of the breccia. It is estimated on the basis of kinetic analysis that the 67975 breccia cooled at a rate of about 0.06 C/min over the temperature range below the liquidus. Heat flow analyses indicate that the 67975 breccia did not cool in the size of the lunar sample but as part of a much larger body. Viscous singering as well as cooling of a melt-clast mixture are investigated.

Uhlmann, D. R.↗
NASA NTRS1982

Fracture and crack growth in orthotropic laminates. Part 1: Analysis of a hybrid, unidirectional laminate with damage

The fracture behavior of unifirectional hybrid (buffer strip) composite laminates is studied. Three particular solutions are discussed: (1) broken fibers in a unidirectional half plane; (2) adjoined half planes of different fiber and matrix properties and (3) the solution of two half planes bounding a third distinct region of finite width. This finite width region represents a buffer strip and the potential of this strip to arrest a crack that originates in one of the half planes is investigated. The analysis is based on a materials modeling approach using the classical shear lag assumption to described the stress transfer between fibers. Explicit fiber and matrix properties of the three regions are retained and changes in the laminate behavior as a function of the relative material properties, buffer strip width and initial crack length are discussed.

Goree, J. G.↗
NASA NTRSDate not supplied

Multiscale Modeling of Thermoplastics Using Atomistic-informed Micromechanics

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

thermoplastics↗
NASA NTRS1979

Regressions by leaps and bounds and biased estimation techniques in yield modeling

The author has identified the following significant results. It was observed that OLS was not adequate as an estimation procedure when the independent or regressor variables were involved in multicollinearities. This was shown to cause the presence of small eigenvalues of the extended correlation matrix A'A. It was demonstrated that the biased estimation techniques and the all-possible subset regression could help in finding a suitable model for predicting yield. Latent root regression was an excellent tool that found how many predictive and nonpredictive multicollinearities there were.

Marquina, N. E.↗
NASA NTRS2018

Plasticity and Damage Modeling of Stress Asymmetry and Dynamic Behavior of AFS Additive Manufactured Aluminum Alloy 2219

The Solid State Additive Manufacturing (AM) process referred as MELD that fabricated the samples in this study, provides a new path for repairing, coating, joining and additive manufacturing metals and metal matrix composites. This research will be the first application of a physics-based microstructure dependent internal state variable (ISV) plasticity and damage material model to capture the mechanical response of an AM Aluminum Alloy (AA) 2219 via the MELD process. In this research, a microstructure-based internal state variable (ISV) plasticity-damage model was used to capture the mechanical behavior of AFS 2219 aluminum alloy. Aeroprobe Corporation, creator and patent holder for the MELD process, fabricated the material by pushing a solid filler rod of AA2219-T861 material through a hollow rotating tool onto an AA2219 T851 plate substrate. As feedstock, solid or powder precursor metals are pushed through a nonconsumable rotating cylindrical tool. Herein, added layers are deposited and metallurgically bonded to substrate material or previously deposited layers by the heat generated from the rotating tool through plastic deformation of the filler material. Once a layer has been added, the tool height increases, and starts the deposition of the next layer. This process results in beneficial properties such as grain refinement, homogenization and reduced porosity (fully dense). This process will experience temperatures similar to those in the weld nugget zone (WNZ) in friction stir welding (FSW), ranging from 0.6-0.9 Tm, with Tm being the melting point of the material. MELD is highly scalable with AA deposition rates reaching over 1000 cm3/hr, which allows for MELD being used for repairs, coatings, and building components. A motivating factor driving the research for physics-based history dependent material modeling of MELD components is the ability to accurately capture the stress-state and strain rate dependence in the material caused by variations in material microstructure from the MELD processing of new or repaired components. The ISV model incorporates microstructural content and is consistent with continuum level kinematics, kinetics, and thermodynamics. These features allow the ISV model to capture large deformations at the structural scale using the kinematic and isotropic hardening, while microscale damage is obtained from the microstructural features. The benefits of the ISV model arise from the inclusion of structure-property relationships identified from microstructural characterization and experimentation. The Bauschinger effect (BE) is an important concept, vital in the accurate prediction of cyclic stress-strain response of ductile materials such as metals. The ISV model has been successfully used to capture the behavior and damage, and the BE of different aluminum alloys and steels. The ISV model uses kinematic and isotropic hardening to help capture deformations of the material at the macro scale. To understand this hardening relationship, calculating the kinematic and isotropic hardening relationship in the material is warranted for a high-fidelity model. Electron Backscattered Diffraction (EBSD) was used to characterize the as-fabricated microstructure, where a fully-dense equiaxed grain morphology with average grain size of 2.5 μm was observed. Microhardness mapping of the as-built structures, monotonic tension and compression experiments at both quasi-static (0.001/s) strain rates, tension-followed-by-compression and compression-followed-by-tension experiments were performed to obtain the set of plasticity and damage constants necessary to capture strain rate and stress state behavior of this additive material. To calibrate the plasticity-damage model, a single set of constants were determined to capture the different stress states the MELD AA2219. One set of the constants was determined from experimental true stress-strain curves for the tension and compression data. Additionally, microstructural information and data from the open literature were used as the other model constants. This research is a first of its kind for AFS AA2219, includes correlating the ISV model to the monotonic experimental results that capture the isotropic and kinematic plasticity mechanical response.

Rivera, O. G.↗
NASA NTRS2009

Bounded Linear Stability Analysis - A Time Delay Margin Estimation Approach for Adaptive Control

This paper presents a method for estimating time delay margin for model-reference adaptive control of systems with almost linear structured uncertainty. The bounded linear stability analysis method seeks to represent the conventional model-reference adaptive law by a locally bounded linear approximation within a small time window using the comparison lemma. The locally bounded linear approximation of the combined adaptive system is cast in a form of an input-time-delay differential equation over a small time window. The time delay margin of this system represents a local stability measure and is computed analytically by a matrix measure method, which provides a simple analytical technique for estimating an upper bound of time delay margin. Based on simulation results for a scalar model-reference adaptive control system, both the bounded linear stability method and the matrix measure method are seen to provide a reasonably accurate and yet not too conservative time delay margin estimation.

Nguyen, Nhan T.↗
NASA NTRS2006

HFGMC Enhancement of MAC/GMC

Additional information about a mathematical model denoted the high-fidelity generalized method of cells (HFGMC) and implementation of the HFGMC within version 4.0 of the MAC/GMC software has become available. MAC/GMC (Micromechanics Analysis Code With Generalized Method of Cells) was a topic of several prior NASA Tech Briefs articles, version 4.0 having been described in "Comprehensive Micromechanics-Analysis Code - Version 4.0" (LEW-17495-1), NASA Tech Briefs, Vol. 29, No. 9 (September 2005), page 54. MAC/GMC predicts elastic and inelastic thermomechanical responses of composite materials. MAC/GMC utilizes the generalized method of cells (GMC) - a model of micromechanics that predicts macroscopic responses of a composite material as functions of the properties, sizes, shapes, and responses of its constituents (e.g., matrix and fibers). The accuracy of the GMC is limited by neglect of coupling between normal and shear stresses. The HFGMC was developed by combining elements of the GMC and a related model, denoted the higher-order theory for functionally graded materials (HOTFGM), that can account for this coupling. Hence, the HFGMC enables simulation of stress and strain with greater accuracy. Some alterations of the MAC/GMC data structure were necessitated by the greater computational complexity of the HFGMC.

Arnold, Steven M.↗
NASA NTRS2014

The Development of Environmental Barrier Coatings for SiCSiC Ceramic Matrix Composites: Challenges and Opportunities

Environmental barrier coatings (EBCs) and SiC/SiC ceramic matrix composites (CMCs) systems will play a crucial role in future turbine engines for hot-section component applications because of their ability to significantly increase engine operating temperatures, reduce engine weight and cooling requirements. The development of prime-reliant environmental barrier coatings is a key to enable the applications of the envisioned CMC components to help achieve next generation engine performance and durability goals. This paper will primarily address the performance requirements and design considerations of environmental barrier coatings for turbine engine applications. The emphasis is placed on current candidate environmental barrier coating systems for SiCSiC CMCs, their performance benefits and design limitations in long-term operation and combustion environments. Major technical barriers in developing advanced environmental barrier coating systems, the coating integrations with next generation CMC turbine components having improved environmental stability, cyclic durability and system performance will be described. The development trends for turbine environmental barrier coating systems by utilizing improved compositions, state-of-the-art processing methods, and simulated environment testing and durability modeling will be discussed.

Environmental Barrier Coatings↗