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,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Controllability of inherently damped large flexible space structures

Graph theoretic techniques are used to study controllability of linear systems which represent large flexible orbiting space systems with inherent damping. The controllability of the pair of matrices representing the system state and control influence matrices is assured when all states in the model are reachable in a digraph sense from at least one input and also when the term rank of a Boolean matrix whose non trivial components are based on the state and control influence matrices has a term rank of the order of the state vector. The damping matrix does not influence the required number of actuators but gives flexibility to the possibility locations of the actuators for which the system is controllable.

Reddy, A. S. S. R.↗
NASA NTRS1983

A tire contact solution technique

An efficient method for calculating the contact boundary and interfacial pressure distribution was developed. This solution technique utilizes the discrete Fourier transform to establish an influence coefficient matrix for the portion of the pressurized tire surface that may be in the contact region. This matrix is used in a linear algebra algorithm to determine the contact boundary and the array of forces within the boundary that are necessary to hold the tire in equilibrium against a specified contact surface. The algorithm also determines the normal and tangential displacements of those points on the tire surface that are included in the influence coefficient matrix. Displacements within and outside the contact region are calculated. The solution technique is implemented with a finite-element tire model that is based on orthotropic, nonlinear shell of revolution elements which can respond to nonaxisymmetric loads. A sample contact solution is presented.

Tielking, J. T.↗
NASA NTRS1984

Radiative transfer in spherical atmospheres

A method for defining spherical model atmospheres in radiative/convective and hydrostatic equilibrium is presented. A finite difference form is found for the transfer equation and a matrix operator is developed as the discrete space analog (in curvilinear coordinates) of a formal integral in plane geometry. Pressure is treated as a function of temperature. Flux conservation is maintained within the energy equation, although the correct luminosity transport must be assigned for any given level of the atmosphere. A perturbed integral operator is used in a complete linearization of the transfer and constraint equations. Finally, techniques for generating stable solutions in economical computer time are discussed.

Kalkofen, W.↗
NASA NTRS1985

Iterative methods for mixed finite element equations

Iterative strategies for the solution of indefinite system of equations arising from the mixed finite element method are investigated in this paper with application to linear and nonlinear problems in solid and structural mechanics. The augmented Hu-Washizu form is derived, which is then utilized to construct a family of iterative algorithms using the displacement method as the preconditioner. Two types of iterative algorithms are implemented. Those are: constant metric iterations which does not involve the update of preconditioner; variable metric iterations, in which the inverse of the preconditioning matrix is updated. A series of numerical experiments is conducted to evaluate the numerical performance with application to linear and nonlinear model problems.

Nakazawa, S.↗
NASA NTRS1987

Nonlinear programming extensions to rational function approximations of unsteady aerodynamics

This paper deals with approximating unsteady generalized aerodynamic forces in the equations of motion of a flexible aircraft. Two methods of formulating these approximations are extended to include both the same flexibility in constraining them and the same methodology in optimizing nonlinear parameters as another currently used 'extended least-squares' method. Optimal selection of 'nonlinear' parameters is made in each of the three methods by use of the same nonlinear (nongradient) optimizer. The objective of the nonlinear optimization is to obtain rational approximations to the unsteady aerodynamics whose state-space realization is of lower order than that required when no optimization of the nonlinear terms is performed. The free 'linear' parameters are determined using least-squares matrix techniques on a Lagrange multiplier formulation of an objective function which incorporates selected linear equality constraints. State-space mathematical models resulting from the different approaches are described, and results are presented which show comparative evaluations from application of each of the extended methods to a numerical example. The results obtained for the example problem show a significant (up to 63 percent) reduction in the number of differential equations used to represent the unsteady aerodynamic forces in linear time-invariant equations of motion as compared to a conventional method in which nonlinear terms are not optimized.

Tiffany, Sherwood H.↗
NASA NTRS1988

Nonlinear programming extensions to rational function approximation methods for unsteady aerodynamic forces

The approximation of unsteady generalized aerodynamic forces in the equations of motion of a flexible aircraft are discussed. Two methods of formulating these approximations are extended to include the same flexibility in constraining the approximations and the same methodology in optimizing nonlinear parameters as another currently used extended least-squares method. Optimal selection of nonlinear parameters is made in each of the three methods by use of the same nonlinear, nongradient optimizer. The objective of the nonlinear optimization is to obtain rational approximations to the unsteady aerodynamics whose state-space realization is lower order than that required when no optimization of the nonlinear terms is performed. The free linear parameters are determined using the least-squares matrix techniques of a Lagrange multiplier formulation of an objective function which incorporates selected linear equality constraints. State-space mathematical models resulting from different approaches are described and results are presented that show comparative evaluations from application of each of the extended methods to a numerical example.

Tiffany, Sherwood H.↗
NASA NTRS1987

A simple method to model truss-beams as equivalent continua

A new method is presented for obtaining stiffness constants for structures modeled as three-dimensional linear elastic beams. The three-dimensional lattice is comprised of multiple two-dimensional lattices. Matrix manipulation is used to determine the contribution of each two-dimensional substructure to the complete three dimensional lattice.

Barton, O., Jr.↗
NASA NTRS1988

Tension fatigue analysis and life prediction for composite laminates

A tension fatigue life prediction methodology for composite laminates is presented. Tension fatigue tests were conducted on quasi-isotropic and orthotropic glass epoxy, graphite epoxy, and glass/graphite epoxy hybrid laminates. Edge delamination onset data were used to generate plots of strain energy release rate as a function of cycles to delamination onset. These plots were then used along with strain energy release rate analyses of delaminations initiating at matrix cracks to predict local delamination onset. Stiffness loss was measured experimentally to account for the accumulation of matrix cracks and for delamination growth. Fatigue failure was predicted by comparing the increase in global strain resulting from stiffness loss to the decrease in laminate failure strain resulting from delaminations forming at matrix cracks through the laminate thickness. Good agreement between measured and predicted lives indicated that the through-thickness damage accumulation model can accurately describe fatigue failure for laminates where the delamination onset behavior in fatigue is well characterized, and stiffness loss can be monitored in real time to account for damage growth.

Obrien, T. K.↗
NASA NTRS1989

Numerical methods in Markov chain modeling

Several methods for computing stationary probability distributions of Markov chains are described and compared. The main linear algebra problem consists of computing an eigenvector of a sparse, usually nonsymmetric, matrix associated with a known eigenvalue. It can also be cast as a problem of solving a homogeneous singular linear system. Several methods based on combinations of Krylov subspace techniques are presented. The performance of these methods on some realistic problems are compared.

Philippe, Bernard↗
NASA NTRS1992

Recursive formulation of operational space control

A recently developed spatial operator algebra approach to modeling and analysis of multibody robotic systems is used to develop O(n) recursive algorithms that compute the operational space mass matrix and the operational space coriolis/centrifugal and gravity terms of an n-link serial manipulator. These algorithms enable an O(n) recursive implementation of operational space control.

Kreutz-Delgado, K.↗
NASA NTRS1991

Regenerator performance in a Vuilleumier refrigerator compared with a third-order numerical model

A three-stage Vuilleumier refrigerator was used to measure the performance of various third stage regenerators. The refrigerator operates between 2.5 and 5.0 Hz and, depending on the material used in the third stage regenerator, achieves temperatures of 8 to 20 K at the cold end of the third stage. This paper presents a comparison of regenerator performance for four regenerator materials: 229 micron diameter spheres of Pb(+)5 pct Sb, 229 micron diameter spheres of brass, 216 micron irregularly-shaped GdRh powder, and a mixture of 229 microns and 762 microns diameter spheres of Pb(+5) pct Sb. The experimental results are compared with a first-order model that neglects the void volume within the regenerator and with a third-order model that considers the effect of pressure oscillations in the regenerator void volume. Experimental results indicate that regenerator losses are dominated by the pressure oscillation in the void volume rather than the mass flow through the temperature gradient in the regenerator. These results are consistent with the third-order numerical model. This model shows that the heat capacity of the gas in the void space as well as the heat capacity of the matrix influences the regenerator performance.

Bradley, P. E.↗
NASA NTRS1999

Inversion of Gravity Fields From the Spacecraft Orbital Data Using an Adjoint Operator Approach

In perturbation approximation, the forward problem of orbital dynamics (equations with initial conditions) is linear with respect to variations of coordinates and/or velocities of the spacecraft and to corresponding variations of the gravity field in the models used. The linear operator adjoint to the linear operator of such forward problem turns out to be instrumental in inversion of differences between observed and predicted coordinates/velocities in terms of the updates of harmonics in the initial gravity field model. Based on this approach, the solution of resulting adjoint problem of orbital dynamics can be used to directly evaluate the matrix of partial derivatives of observable differences with respect to the gravity field harmonics. General discussion of the adjoint problem of orbital dynamics is given and an example of a mathematical formalism for the practical retrieval algorithm is presented.

Ustinov, E. A.↗
NASA NTRS2000

Planned Axial Reorientation Investigation on Sloshsat

This paper details the design and logic of an experimental investigation to study axial reorientation in low gravity. The Sloshsat free-flyer is described. The planned axial reorientation experiments and test matrixes are presented. Existing analytical tools are discussed. Estimates for settling range from 64 to 1127 seconds. The planned experiments are modelled using computational fluid dynamics. These models show promise in reducing settling estimates and demonstrate the ability of pulsed high thrust settling to emulate lower thrust continuous firing.

Chato, David J.↗
NASA NTRS2001

Characterization of Composites Response at High Rates of Deformation

The objective of the proposed research is to experimentally study the effect of strain rate on mechanical response (deformation and failure) of IM-7/977-2 carbon fiber/epoxy matrix composites. The experimental data will provide the information needed for the development of a nonlinear, rate dependent deformation and strength model for this material that can subsequently be used in design.

Gilat, Amos↗
NASA NTRS2004

Analytical Solution for the Critical Velocity of Pushing/Engulfment Transition

The distribution of ceramic particles in a metal matrix composite material depends primarily on the interaction of the particles with the solid/liquid interface during the solidification process. A numerical model that describes the evolution of the shape of the solid/liquid interface in the proximity of a foreign particle will presented in this paper. The model accounts for the influence of the temperature gradient and the Gibbs-Thomson and disjoining pressure effects. It shows that for the systems characterized by k(sub p) < k(sub L) the disjoining pressure causes the interface curvature to change its sign in the close-contact particle/interface region. It also shows that the increase of the temperature gradient diminishes the effect of the disjoining pressure. The analysis of the numerical results obtained for a large range of processing conditions and materials parameters has led to the development of an analytical solution for the critical velocity of pushing/engulfinent transition. The theoretical results will be discussed and compared with the experimental measurements performed under microgravity conditions.

Catalina, Adrian V.↗
NASA NTRS2006

Multiscale Analysis of Delamination of Carbon Fiber-Epoxy Laminates with Carbon Nanotubes

A multi-scale analysis is presented to parametrically describe the Mode I delamination of a carbon fiber/epoxy laminate. In the midplane of the laminate, carbon nanotubes are included for the purposes of selectively enhancing the fracture toughness of the laminate. To analyze carbon fiber epoxy carbon nanotube laminate, the multi-scale methodology presented here links a series of parameterizations taken at various length scales ranging from the atomistic through the micromechanical to the structural level. At the atomistic scale molecular dynamics simulations are performed in conjunction with an equivalent continuum approach to develop constitutive properties for representative volume elements of the molecular structure of components of the laminate. The molecular-level constitutive results are then used in the Mori-Tanaka micromechanics to develop bulk properties for the epoxy-carbon nanotube matrix system. In order to demonstrate a possible application of this multi-scale methodology, a double cantilever beam specimen is modeled. An existing analysis is employed which uses discrete springs to model the fiber bridging affect during delamination propagation. In the absence of empirical data or a damage mechanics model describing the effect of CNTs on fracture toughness, several tractions laws are postulated, linking CNT volume fraction to fiber bridging in a DCB specimen. Results from this demonstration are presented in terms of DCB specimen load-displacement responses.

Riddick, Jaret C.↗
NASA NTRS2008

Further Developments in Modeling Creep Effects Within Structural SiC/SiC Components

Anticipating the implementation of advanced SiC/SiC composites into turbine section components of future aero-propulsion engines, the primary objective of this on-going study is to develop physics-based analytical and finite-element modeling tools to predict the effects of constituent creep on SiC/SiC component service life. A second objective is to understand how to possibly manipulate constituent materials and processes in order to minimize these effects. Focusing on SiC/SiC components experiencing through-thickness stress gradients (e.g., airfoil leading edge), prior NASA creep modeling studies showed that detrimental residual stress effects can develop globally within the component walls which can increase the risk of matrix cracking. These studies assumed that the SiC/SiC composites behaved as isotropic viscoelastic continuum materials with creep behavior that was linear and symmetric with stress and that the creep parameters could be obtained from creep data as experimentally measured in-plane in the fiber direction of advanced thin-walled 2D SiC/SiC panels. The present study expands on those prior efforts by including constituent behavior with non-linear stress dependencies in order to predict such key creep-related SiC/SiC properties as time-dependent matrix stress, constituent creep and content effects on composite creep rates and rupture times, and stresses on fiber and matrix during and after creep.

Lang, Jerry↗
NASA NTRS2014

Line Interference Effects Using a Refined Robert-Bonamy Formalism: the Test Case of the Isotropic Raman Spectra of Autoperturbed N2

A symmetrized version of the recently developed refined Robert-Bonamy formalism [Q. Ma, C. Boulet, and R. H. Tipping, J. Chem. Phys. 139, 034305 (2013)] is proposed. This model takes into account line coupling effects and hence allows the calculation of the off-diagonal elements of the relaxation matrix, without neglecting the rotational structure of the perturbing molecule. The formalism is applied to the isotropic Raman spectra of autoperturbed N2 for which a benchmark quantum relaxation matrix has recently been proposed. The consequences of the classical path approximation are carefully analyzed. Methods correcting for effects of inelasticity are considered. While in the right direction, these corrections appear to be too crude to provide off diagonal elements which would yield, via the sum rule, diagonal elements in good agreement with the quantum results. In order to overcome this difficulty, a re-normalization procedure is applied, which ensures that the off-diagonal elements do lead to the exact quantum diagonal elements. The agreement between the (re-normalized) semi-classical and quantum relaxation matrices is excellent, at least for the Raman spectra of N2, opening the way to the analysis of more complex molecular systems.

nitrogen↗