Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MATRICES”

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 613 records · Page 34

Property changes induced by the space environment in composite materials on LDEF: Solar array materials passive LDEF experiment A0171 (SAMPLE)

Surface modifications to composite materials induced by long term exposure in low earth orbit (LEO) were dominated by atomic oxygen erosion and micrometeoroid and space debris impacts. As expected, calculated erosion rates were peculiar to material type and within the predicted order of magnitude. Generally, about one ply of the carbon fiber composites was eroded during the 70 month LDEF experiment. Matrix erosion was greater than fiber erosion and was more evident for a polysulfone matrix than for epoxy matrices. Micrometeoroid and space debris impacts resulted in small (less than 1mm) craters and splattered contaminants on all samples. Surfaces became more diffuse and darker with small increases in emissivity and absorption. Tensile strength decreased roughly with thickness loss, and epoxy matrices apparently became slightly embrittled, probably as a result of continued curing under UV and/or electron bombardment. However, changes in the ultimate yield stress of the carbon reinforced epoxy composites correlate neither with weave direction nor fiber type. Unexpected developments were the discovery of new synergistic effects of the space environment in the interaction of atomic oxygen and copious amounts of contamination and in the induced luminescence of many materials.

Source record↗

Discrete Kalman filtering equations of second-order form for control-structure interaction simulations

A second-order form of discrete Kalman filtering equations is proposed as a candidate state estimator for efficient simulations of control-structure interactions in coupled physical coordinate configurations as opposed to decoupled modal coordinates. The resulting matrix equation of the present state estimator consists of the same symmetric, sparse N x N coupled matrices of the governing structural dynamics equations as opposed to unsymmetric 2N x 2N state space-based estimators. Thus, in addition to substantial computational efficiency improvement, the present estimator can be applied to control-structure design optimization for which the physical coordinates associated with the mass, damping and stiffness matrices of the structure are needed instead of modal coordinates.

Park, K. C.↗

Formally biorthogonal polynomials and a look-ahead Levinson algorithm for general Toeplitz systems

Systems of linear equations with Toeplitz coefficient matrices arise in many important applications. The classical Levinson algorithm computes solutions of Toeplitz systems with only O(n(sub 2)) arithmetic operations, as compared to O(n(sub 3)) operations that are needed for solving general linear systems. However, the Levinson algorithm in its original form requires that all leading principal submatrices are nonsingular. An extension of the Levinson algorithm to general Toeplitz systems is presented. The algorithm uses look-ahead to skip over exactly singular, as well as ill-conditioned leading submatrices, and, at the same time, it still fully exploits the Toeplitz structure. In our derivation of this algorithm, we make use of the intimate connection of Toeplitz matrices with formally biorthogonal polynomials.

Freund, Roland W.↗

Efficient, massively parallel eigenvalue computation

In numerical simulations of disordered electronic systems, one of the most common approaches is to diagonalize random Hamiltonian matrices and to study the eigenvalues and eigenfunctions of a single electron in the presence of a random potential. An effort to implement a matrix diagonalization routine for real symmetric dense matrices on massively parallel SIMD computers, the Maspar MP-1 and MP-2 systems, is described. Results of numerical tests and timings are also presented.

Huo, Yan↗

Parallel O(log n) algorithms for open- and closed-chain rigid multibody systems based on a new mass matrix factorization technique

In this paper, parallel O(log n) algorithms for computation of rigid multibody dynamics are developed. These parallel algorithms are derived by parallelization of new O(n) algorithms for the problem. The underlying feature of these O(n) algorithms is a drastically different strategy for decomposition of interbody force which leads to a new factorization of the mass matrix (M). Specifically, it is shown that a factorization of the inverse of the mass matrix in the form of the Schur Complement is derived as M(exp -1) = C - B(exp *)A(exp -1)B, wherein matrices C, A, and B are block tridiagonal matrices. The new O(n) algorithm is then derived as a recursive implementation of this factorization of M(exp -1). For the closed-chain systems, similar factorizations and O(n) algorithms for computation of Operational Space Mass Matrix lambda and its inverse lambda(exp -1) are also derived. It is shown that these O(n) algorithms are strictly parallel, that is, they are less efficient than other algorithms for serial computation of the problem. But, to our knowledge, they are the only known algorithms that can be parallelized and that lead to both time- and processor-optimal parallel algorithms for the problem, i.e., parallel O(log n) algorithms with O(n) processors. The developed parallel algorithms, in addition to their theoretical significance, are also practical from an implementation point of view due to their simple architectural requirements.

Fijany, Amir↗

Application of the target decomposition theorem to a polarimetric random media model

With advances in polarimetric radar measurements of land surfaces, the need for understanding the underlying scattering mechanisms and dominant target features has become the focus of many studies. In particular, the maximum use of the polarimetric information to identify and/or separate parameters related to the surface features such as vegetation thickness, structure, water content, and soil surface characteristics will enhance the possibility of using polarimetric radars for monitoring the earth's surface from space. In this paper, Cloude's decomposition theorem is applied to a polarimetric random media model to simulate the radar measurements of vegetated canopies. The vegetated canopies are modeled as a three layer discrete random medium with leaves and branches in the first layer, tree trunks in the second layer, and a half space of homogeneous ground with rough interface as the bottom layer. The distorted born approximation (DBA) has been used to compute full Mueller matrix of the canopy, using canonical dielectric objects such as thin discs and cylinders as leaves, branches, and trunks, respectively. The Mueller matrix and the derived covariance matrix contain information on the second order statistics of radar signals at various polarizations from the canopy. To decompose the covariance matrix to its constituent targets, the eigenvalues and eigenvectors of the covariance matrix are computed in terms of the physical parameters of the canopy. In addition, each eigenvector explicitly shows the scattering mechanisms such as odd and even reflections in the canopy. Cloude's decomposition theorem is applied using the Pauli spin matrices as a basis and an expression for the degree of disorder or the entropy for the vegetated surface is found. Then, the physical parameters estimated from in situ measurements are used in the random media to obtain realistic covariance matrices. As a result, the sensitivity of the eigenvalue spectrums and the coefficients resulting from the target decomposition theorem to the physical parameters of the canopy are examined and the possible use of Cloude's theorem to estimate vegetation and soil parameters is discussed.

Saatchi, Sasan S.↗

Steady induction effects in geomagnetism. Part 1C: Geomagnetic estimation of steady surficial core motions: Application to the definitive geomagnetic reference field models

In the source-free mantle/frozen-flux core magnetic earth model, the non-linear inverse steady motional induction problem was solved using the method presented in Part 1B. How that method was applied to estimate steady, broad-scale fluid velocity fields near the top of Earth's core that induce the secular change indicated by the Definitive Geomagnetic Reference Field (DGRF) models from 1945 to 1980 are described. Special attention is given to the derivation of weight matrices for the DGRF models because the weights determine the apparent significance of the residual secular change. The derived weight matrices also enable estimation of the secular change signal-to-noise ratio characterizing the DGRF models. Two types of weights were derived in 1987-88: radial field weights for fitting the evolution of the broad-scale portion of the radial geomagnetic field component at Earth's surface implied by the DGRF's, and general weights for fitting the evolution of the broad-scale portion of the scalar potential specified by these models. The difference is non-trivial because not all the geomagnetic data represented by the DGRF's constrain the radial field component. For radial field weights (or general weights), a quantitatively acceptable explication of broad-scale secular change relative to the 1980 Magsat epoch must account for 99.94271 percent (or 99.98784 percent) of the total weighted variance accumulated therein. Tolerable normalized root-mean-square weighted residuals of 2.394 percent (or 1.103 percent) are less than the 7 percent errors expected in the source-free mantle/frozen-flux core approximation.

Voorhies, Coerte V.↗

Micromechanical characterization of nonlinear behavior of advanced polymer matrix composites

Due to the presence of curing stresses and oriented crystalline structures in the matrix of polymer matrix fiber composites, the in situ nonlinear properties of the matrix are expected to be rather different from those of the bulk resin. A plane stress micromechanical model was developed to retrieve the in situ elastic-plastic properties of Narmco 5260 and Amoco 8320 matrices from measured elastic-plastic properties of IM7/5260 and IM7/8320 advance composites. In the micromechanical model, the fiber was assumed to be orthotropically elastic and the matrix to be orthotropic in elastic and plastic properties. The results indicate that both in situ elastic and plastic properties of the matrices are orthotropic.

Gates, Thomas S.↗

Distributed parameter modeling of repeated truss structures

A new approach to find homogeneous models for beam-like repeated flexible structures is proposed which conceptually involves two steps. The first step involves the approximation of 3-D non-homogeneous model by a 1-D periodic beam model. The structure is modeled as a 3-D non-homogeneous continuum. The displacement field is approximated by Taylor series expansion. Then, the cross sectional mass and stiffness matrices are obtained by energy equivalence using their additive properties. Due to the repeated nature of the flexible bodies, the mass, and stiffness matrices are also periodic. This procedure is systematic and requires less dynamics detail. The first step involves the homogenization from a 1-D periodic beam model to a 1-D homogeneous beam model. The periodic beam model is homogenized into an equivalent homogeneous beam model using the additive property of compliance along the generic axis. The major departure from previous approaches in literature is using compliance instead of stiffness in homogenization. An obvious justification is that the stiffness is additive at each cross section but not along the generic axis. The homogenized model preserves many properties of the original periodic model.

Wang, Han-Ching↗

Component modes damping assignment methodology for articulated, multiflexible body structures

To simulate the dynamical motion of articulated, multiflexible body structures, one can use multibody simulation packages such as DISCOS. To this end, one must supply appropriate reduced-order models for all of the flexible components involved. The component modes projection and assembly model reduction (COMPARE) methodology is one way to construct these reduced-order component models, which when reassembled capture important system input-to-output mapping of the full-order model at multiple system configurations of interest. In conjunction, we must also supply component damping matrices which when reassembled generate a system damping matrix that has certain desirable properties. The problem of determining the damping factors of components' modes to achieve a given system damping matrix is addressed here. To this end, we must establish from first principles a matrix-algebraic relation between the system's modal damping matrix and the components' modal damping matrices. An unconstrained/constrained optimization problem can then be formulated to determine the component modes' damping factors that best satisfy that matrix-algebraic relation. The effectiveness of the developed methodology, called ModeDamp, has been successfully demonstrated on a high-order, finite element model of the Galileo spacecraft.

Lee, Allan Y.↗

Performance of customized DCT quantization tables on scientific data

We show that it is desirable to use data-specific or customized quantization tables for scaling the spatial frequency coefficients obtained using the Discrete Cosine Transform (DCT). DCT is widely used for image and video compression (MP89, PM93) but applications typically use default quantization matrices. Using actual scientific data gathered from divers sources such as spacecrafts and electron-microscopes, we show that the default compression/quality tradeoffs can be significantly improved upon by using customized tables. We also show that significant improvements are possible for the standard test images Lena and Baboon. This work is part of an effort to develop a practical scheme for optimizing quantization matrices for any given image or video stream, under any given quality or compression constraints.

Ratnakar, Viresh↗

SEQ-REVIEW: A tool for reviewing and checking spacecraft sequences

A key component of JPL's strategy to make space missions faster, better and cheaper is the Advanced Multi-Mission Operations System (AMMOS), a ground software intensive system currently in use and in further development. AMMOS intends to eliminate the cost of re-engineering a ground system for each new JPL mission. This paper discusses SEQ-REVIEW, a component of AMMOS that was designed to facilitate and automate the task of reviewing and checking spacecraft sequences. SEQ-REVIEW is a smart browser for inspecting files created by other sequence generation tools in the AMMOS system. It can parse sequence-related files according to a computer-readable version of a 'Software Interface Specification' (SIS), which is a standard document for defining file formats. It lets users display one or several linked files and check simple constraints using a Basic-like 'Little Language'. SEQ-REVIEW represents the first application of the Quality Function Development (QFD) method to sequence software development at JPL. The paper will show how the requirements for SEQ-REVIEW were defined and converted into a design based on object-oriented principles. The process starts with interviews of potential users, a small but diverse group that spans multiple disciplines and 'cultures'. It continues with the development of QFD matrices that related product functions and characteristics to user-demanded qualities. These matrices are then turned into a formal Software Requirements Document (SRD). The process concludes with the design phase, in which the CRC (Class, Responsibility, Collaboration) approach was used to convert requirements into a blueprint for the final product.

Maldague, Pierre F.↗

A Fully Implicit Time Accurate Method for Hypersonic Combustion: Application to Shock-induced Combustion Instability

A new fully implicit, time accurate algorithm suitable for chemically reacting, viscous flows in the transonic-to-hypersonic regime is described. The method is based on a class of Total Variation Diminishing (TVD) schemes and uses successive Gauss-Siedel relaxation sweeps. The inversion of large matrices is avoided by partitioning the system into reacting and nonreacting parts, but still maintaining a fully coupled interaction. As a result, the matrices that have to be inverted are of the same size as those obtained with the commonly used point implicit methods. In this paper we illustrate the applicability of the new algorithm to hypervelocity unsteady combustion applications. We present a series of numerical simulations of the periodic combustion instabilities observed in ballistic-range experiments of blunt projectiles flying at subdetonative speeds through hydrogen-air mixtures. The computed frequencies of oscillation are in excellent agreement with experimental data.

Yungster, Shaye↗

The effect of mixed mode precracking on the mode 1 fracture toughness of composite laminates

We subjected double cantilever beam specimens from four different composite materials to mixed-mode precracking. Three different precracking mode 1 to mode 2 ratios were used--1 to 4, 1 to 1, and 4 to 1. Following precracking the specimens were tested for mode I fracture toughness. The mixed-mode precracking often influenced the mode 1 toughness and its influence persisted for as much as 60 mm of mode 1 crack growth. We tested composites with untoughened matrices, composites with rubber-toughened matrices, and composites with interlayer toughening. Depending on material type and precracking mode ratio, the precracking could cause either a significant increase or a significant decrease in the mode 1 fracture toughness.

Shankar, Prashanth↗

The effect of interface properties on nickel base alloy composites

This program was performed to assess the extent to which mechanical behavior models can predict the properties of sapphire fiber/nickel aluminide matrix composites and help guide their development by defining improved combinations of matrix and interface coating. The program consisted of four tasks: 1) selection of the matrices and interface coating constituents using a modeling-based approach; 2) fabrication of the selected materials; 3) testing and evaluation of the materials; and 4) evaluation of the behavior models to develop recommendations. Ni-50Al and Ni-20AI-30Fe (a/o) matrices were selected which gave brittle and ductile behavior, respectively, and an interface coating of PVD YSZ was selected which provided strong bonding to the sapphire fiber. Significant fiber damage and strength loss was observed in the composites which made straightforward comparison of properties with models difficult. Nevertheless, the models selected generally provided property predictions which agreed well with results when fiber degradation was incorporated. The presence of a strong interface bond was felt to be detrimental in the NiAI MMC system where low toughness and low strength were observed.

Groves, M.↗

Prioritization methodology for chemical replacement

Since United States of America federal legislation has required ozone depleting chemicals (class 1 & 2) to be banned from production, The National Aeronautics and Space Administration (NASA) and industry have been required to find other chemicals and methods to replace these target chemicals. This project was initiated as a development of a prioritization methodology suitable for assessing and ranking existing processes for replacement 'urgency.' The methodology was produced in the form of a workbook (NASA Technical Paper 3421). The final workbook contains two tools, one for evaluation and one for prioritization. The two tools are interconnected in that they were developed from one central theme - chemical replacement due to imposed laws and regulations. This workbook provides matrices, detailed explanations of how to use them, and a detailed methodology for prioritization of replacement technology. The main objective is to provide a GUIDELINE to help direct the research for replacement technology. The approach for prioritization called for a system which would result in a numerical rating for the chemicals and processes being assessed. A Quality Function Deployment (QFD) technique was used in order to determine numerical values which would correspond to the concerns raised and their respective importance to the process. This workbook defines the approach and the application of the QFD matrix. This technique: (1) provides a standard database for technology that can be easily reviewed, and (2) provides a standard format for information when requesting resources for further research for chemical replacement technology. Originally, this workbook was to be used for Class 1 and Class 2 chemicals, but it was specifically designed to be flexible enough to be used for any chemical used in a process (if the chemical and/or process needs to be replaced). The methodology consists of comparison matrices (and the smaller comparison components) which allow replacement technology to be quantitatively compared in several categories, and a QFD matrix which allows process/chemical pairs to be rated against one another for importance (using consistent categories). Depending on the need for application, one can choose the part(s) needed or have the methodology completed in its entirety. For example, if a program needs to show the risk of changing a process/chemical one may choose to use part of Matrix A and Matrix C. If a chemical is being used, and the process must be changed; one might use the Process Concerns part of Matrix D for the existing process and all possible replacement processes. If an overall analysis of a program is needed, one may request the QFD to be completed.

Cruit, Wendy↗

Al2O3 fiber strength degradation in metal and intermetallic matrix composites

The mechanisms for fiber damage in single crystal Al2O3 fiber-reinforced composites were investigated. Both fiber fragmentation and fiber strength degradation were observed in composites with a variety of matrix compositions. Four mechanisms that may be contributing to the fiber strength loss have been proposed and include matrix reaction, reaction with binders, residual stress-induced damage, and pressure from hot pressing. The effect of matrix reaction was separated from the other three effects by sputter-coating the matrices on cleaned fibers and annealing with a temperature profile that simulates processing conditions. These experiments revealed that Y and Cr in FeCrAlY base alloys and Zr in NiAl alloys reacted with the fiber, and grooves and adherent particles were formed on the fiber surface which were responsible for the strength loss. The effects of the matrix reaction appeared to dominate over the other possible mechanisms, although evidence for reaction with binders was also found. Ridges on the fiber surface, which reflected the grain boundaries of the matrix, were also observed. In order for single-crystal Al2O3 to be used as a fiber in MMC's and IMC's, a matrix or protective coating which minimizes matrix reaction during processing will be necessary. Of the matrices investigated, the Thermo-span(sup TM) alloy was the least damaging to fiber properties.

ALUMINA/IRON CHROMIUM ALUMINUM↗

Development of Finite Elements for Two-Dimensional Structural Analysis Using the Integrated Force Method

The Integrated Force Method has been developed in recent years for the analysis of structural mechanics problems. This method treats all independent internal forces as unknown variables that can be calculated by simultaneously imposing equations of equilibrium and compatibility conditions. In this paper a finite element library for analyzing two-dimensional problems by the Integrated Force Method is presented. Triangular- and quadrilateral-shaped elements capable of modeling arbitrary domain configurations are presented. The element equilibrium and flexibility matrices are derived by discretizing the expressions for potential and complementary energies, respectively. The displacement and stress fields within the finite elements are independently approximated. The displacement field is interpolated as it is in the standard displacement method, and the stress field is approximated by using complete polynomials of the correct order. A procedure that uses the definitions of stress components in terms of an Airy stress function is developed to derive the stress interpolation polynomials. Such derived stress fields identically satisfy the equations of equilibrium. Moreover, the resulting element matrices are insensitive to the orientation of local coordinate systems. A method is devised to calculate the number of rigid body modes, and the present elements are shown to be free of spurious zero-energy modes. A number of example problems are solved by using the present library, and the results are compared with corresponding analytical solutions and with results from the standard displacement finite element method. The Integrated Force Method not only gives results that agree well with analytical and displacement method results but also outperforms the displacement method in stress calculations.

Kaljevic, Igor↗