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 739 records · Page 41

Nonlinear response of resin matrix laminates using endochronic theory

The nonlinear response of laminated, resin matrix fibrous composites is modeled using orthotropic endochronic theory. The theory is formulated in terms of elastic constants and endochronic parameters characterizing the linear and nonlinear response, respectively, of unidirectional composites. All constants and parameters can be determined from normal (tension and/or compression) and shear tests on unidirectional and off-axis specimens. The nonlinear constitutive response relations for the unidirectional lamina are presented and the procedure for determining constants and parameters from test is described. The results are then used to predict the nonlinear response of unidirectional laminae and angle-ply laminates. Comparison between theory and experiment for compression loading of AS4/3502 graphite-epoxy, angle-ply laminates shows excellent correlation.

Mathison, S. R.↗

Ultrasonic and micromechanical study of damage and elastic properties of SiC/RBSN ceramic composites

Ultrasonic techniques are employed to develop methods for nondestructive evaluation of elastic properties and damage in SiC/RBSN composites. To incorporate imperfect boundary conditions between fibers and matrix into a micromechanical model, a model of fibers having effective anisotropic properties is introduced. By inverting Hashin's (1979) microstructural model for a composite material with microscopic constituents the effective fiber properties were found from ultrasonic measurements. Ultrasonic measurements indicate that damage due to thermal shock is located near the surface, so the surface wave is most appropriate for estimation of the ultimate strength reduction and critical temperature of thermal shock. It is concluded that bonding between laminates of SiC/RBSN composites is severely weakened by thermal oxidation. Generally, nondestructive evaluation of thermal oxidation effects and thermal shock shows good correlation with measurements previously performed by destructive methods.

Chu, Y. C.↗

Linear analysis of opto-mechanical systems

A general framework is presented for the matrix-form, linear optical model analysis of controlled optomechanical systems; the models are conjoined with linear models of structures and controls to compute system performance as a function of optics, structures, and control parameters. Covariance analysis, optimization, and estimation/simulation are used. Attention is given to a tolerancing example for the Hubble Space Telescope's Wide Field and Planetary Camera, which involves the creation of a linear model of residual pupil shear.

Redding, David C.↗

Progressive Failure And Life Prediction of Ceramic and Textile Composites

An engineering approach to predict the fatigue life and progressive failure of multilayered composite and textile laminates is presented. Analytical models which account for matrix cracking, statistical fiber failures and nonlinear stress-strain behavior have been developed for both composites and textiles. The analysis method is based on a combined micromechanics, fracture mechanics and failure statistics analysis. Experimentally derived empirical coefficients are used to account for the interface of fiber and matrix, fiber strength, and fiber-matrix stiffness reductions. Similar approaches were applied to textiles using Repeating Unit Cells. In composite fatigue analysis, Walker's equation is applied for matrix fatigue cracking and Heywood's formulation is used for fiber strength fatigue degradation. The analysis has been compared with experiment with good agreement. Comparisons were made with Graphite-Epoxy, C/SiC and Nicalon/CAS composite materials. For textile materials, comparisons were made with triaxial braided and plain weave materials under biaxial or uniaxial tension. Fatigue predictions were compared with test data obtained from plain weave C/SiC materials tested at AS&M. Computer codes were developed to perform the analysis. Composite Progressive Failure Analysis for Laminates is contained in the code CPFail. Micromechanics Analysis for Textile Composites is contained in the code MicroTex. Both codes were adapted to run as subroutines for the finite element code ABAQUS and CPFail-ABAQUS and MicroTex-ABAQUS. Graphic user interface (GUI) was developed to connect CPFail and MicroTex with ABAQUS.

Xue, David Y.↗

Simulation of Non-Acoustic Combustion Instability in a Hybrid Rocket Motor

A transient model of a hybrid motor was formulated to study the cause and elimination of non-acoustic combustion instability. The transient model was used to simulate four key tests out of a series of seventeen hybrid motor tests conducted by Thiokol, Rocketdyne and Martin Marietta at NASA/Marshall Space Flight Center (NASA/MSFC). These tests were performed under the Hybrid Propulsion Technology for Launch Vehicle Boosters (HPTLVB) program. The first test resulted in stable combustion. The second test resulted in large-amplitude, 6.5 Hz chamber pressure oscillations that gradually damped away by the end of the test. The third test resulted in large-amplitude, 7.5 Hz chamber pressure oscillations that were sustained throughout the test. The seventh test resulted in the elimination of combustion instability with the installation of an orifice immediately upstream of the injector. The formulation and implementation of the model are the scope of this presentation. The current model is an independent continuation of modeling presented previously by joint Thiokol-Rocketdyne collaborators Boardman, Hawkins, Wassom, and Claflin. The previous model simulated an unstable IR&D hybrid motor test performed by Thiokol. There was very good agreement between the model and the test data. Like the previous model, the current model was developed using Matrix-x simulation software. However, the tests performed at NASA/MSFC under the HPTLVB program were actually simulated. In the current model, the hybrid motor consisting of the liquid oxygen (LOX) injector, the multi-port solid fuel grain and the nozzle was simulated. Also, simulated in the model was the LOX feed system consisting of the tank, venturi, valve and feed lines. All components of the hybrid motor and LOX feed system are treated by a lumped-parameter approach. Agreement between the results of the transient model and the actual test data was very good. This agreement between simulated and actual test data indicated that the combustion instability in the hybrid motor was due to two causes. The first cause was a LOX feed system of insufficient stiffness. The second cause was a LOX injector with an impedance or pressure drop that was too low to provide damping against the feed system oscillations. Also, it was discovered that testing with a new grain of solid fuel sustained the combustion instability. However, testing with a used grain of solid fuel caused the combustion instability to gradually decay.

Rocker, Marvin↗

Simulation of Non-Acoustic Combustion Instability in a Hybrid Rocket Motor

A transient model of a hybrid motor was formulated to study the cause and elimination of non-acoustic combustion instability. The transient model was used to simulate four key tests out of a series of seventeen hybrid motor tests conducted by Thiokol, Rocketdyne and Martin Marietta at NASA/Marshall Space Flight Center (NASAIMSFC). These tests were performed under the Hybrid Propulsion Technology for Launch Vehicle Boosters (HPTLVB) program. The first test resulted in stable combustion. The second test resulted in large-amplitude, 6.5 Hz chamber pressure oscillations that gradually damped away by the end of the test. The third test resulted in large-amplitude, 7.5 Hz chamber pressure oscillations that were sustained throughout the test. The seventh test resulted in the elimination of combustion instability with the installation of an orifice immediately upstream of the injector. The formulation and implementation of the model are the scope of this presentation. The current model is an independent continuation of modeling presented previously by joint Thiokol-Rocketdyne collaborators Boardman, Hawkins, Wassom, and Claflin. The previous model simulated an unstable IR&D hybrid motor test performed by Thiokol. There was very good agreement between the model and the test data. Like the previous model, the current model was developed using Matrix-x simulation software. However, the tests performed at NASA/MSFC under the HPTLVB program were actually simulated. In the current model, the hybrid motor consisting of the liquid oxygen (LOX) injector, the multi-port solid fuel grain and the nozzle was simulated. Also, simulated in the model was the LOX feed system consisting of the tank, venturi, valve and feed lines. All components of the hybrid motor and LOX feed system are treated by a lumped-parameter approach. Agreement between the results of the transient model and the actual test data was very good. This agreement between simulated and actual test data indicated that the combustion instability in the hybrid motor was due to two causes. The first cause was a LOX feed system of insufficient stiffness. The second cause was a LOX injector with an impedance or pressure drop that was too low to provide damping against the feed system oscillations. Also, it was discovered that testing with a new grain of solid fuel sustained the combustion instability. However, testing with a used grain of solid fuel caused the combustion instability to gradually decay.

Rocker, Marvin↗

Structural Modeling of a Five-Meter Thin Film Inflatable Antenna/Concentrator With Rigidized Support Struts

Dynamic characterization of a non-rigidized thin film inflatable antenna/solar concentrator structure with rigidized composite support struts is described in detail. A two-step finite element modeling approach in MSC/NASTRAN is utilized, consisting of: (1) a nonlinear static pressurization procedure used to obtain the updated stiffness matrix, and (2) a modal "restart" eigen solution that uses the modified stiffness matrix. Unique problems encountered in modeling of this large 5-m lightweight inflatable are identified, including considerable difficulty in obtaining convergence in the nonlinear pressurization solution. It was found that the extremely thin polyimide film material (.001 in or I mil) presents tremendous problems in obtaining a converged solution when internal pressure loading is applied. It was concluded that the ratios of film thickness to other geometric dimensions such as torus cross-sectional and ring diameter and lenticular diameter are the critical parameters for convergence of the pressurization procedure. Comparison of finite element predictions for frequency and mode shapes with experimental results indicated reasonable agreement considering the complexity of the structure, the film-to-air interaction, and the nonlinear material properties of the film. It was also concluded that analysis should be done using different finite element to codes to determine if a more robust and stable solution can be obtained.

Smalley, Kurt B.↗

The Role of Magnetograms in Testing Models of Coronal Mass Ejections

A general review of current and future magnetographic capabilities is presented with emphasis on how they address observational predictions of CME models (e.g., the model/observation matrix developed in the SHINE 2000 workshop). Existing line-of-sight magnetographs can observe certain structural and kinematic properties of CME source regions when they are on the visible disk. These include polarity configuration, location and structure of the magnetic inversion line, proper motions of magnetic elements, and time variation of flux. Additional chromospheric magnetograms provide information on how the field spreads with height and, if properly sequenced with photospheric magnetograms, can determine whether flux is emerging or submerging. Vector magnetographs can provide information about currents and magnetic shear but will only begin to attain the spatial and temporal coverage appropriate to CMEs with the advent of the SOLIS vector spectromagnetograph and, hopefully, its copies at several widely spaced longitudes. Important additions to the SHINE 2000 matrix are required to define an observational program. These include spatial and temporal resolutions and extents as well as needed sensitivity.

Jones, Harrison P.↗

Dynamics of a 4x6-Meter Thin Film Elliptical Inflated Membrane for Space Applications

Dynamic characterization of a thin film inflatable elliptical structure is described in detail. A two-step finite element modeling approach in MSC/NASTRAN is utilized, consisting of (1) a nonlinear static pressurization procedure used to obtain the updated stiffness matrix, and (2) a modal "restart" eigen solution that uses the modified stiffness matrix. Unique problems encountered in modeling of this large Hexameter lightweight inflatable arc identified, including considerable difficulty in obtaining convergence in the nonlinear finite element pressurization solution. It was found that the extremely thin polyimide film material (.001 in or 1 mil) presents tremendous problems in obtaining a converged solution when internal pressure loading is applied. Approaches utilized to overcome these difficulties are described. Comparison of finite element predictions for frequency and mode shapes of the inflated structure with closed-form solutions for a flat pre-tensioned membrane indicate reasonable agreement.

Casiano, Matthew J.↗

Effects of Solute Concentrations on Kinetic Pathways in Ni-Al-Cr Alloys

The kinetic pathways resulting from the formation of coherent gamma'-precipitates from the gamma-matrix are studied for two Ni-Al-Cr alloys with similar gamma'-precipitate volume fractions at 873 K. The details of the phase decompositions of Ni-7.5Al-8.5Cr at.% and Ni-5.2Al-14.2Cr at.% for aging times from 1/6 to 1024 h are investigated by atom-probe tomography, and are found to differ significantly from a mean-field description of coarsening. The morphologies of the gamma'-precipitates of the alloys are similar, though the degrees of gamma'-precipitate coagulation and coalescence differ. Quantification within the framework of classical nucleation theory reveals that differences in the chemical driving forces for phase decomposition result in differences in the nucleation behavior of the two alloys. The temporal evolution of the gamma'-precipitate average radii and the gamma-matrix supersaturations follow the predictions of classical coarsening models. The compositional trajectories of the gamma-matrix phases of the alloys are found to follow approximately the equilibrium tie-lines, while the trajectories of the gamma'-precipitates do not, resulting in significant differences in the partitioning ratios of the solute elements.

Booth-Morrison, Christopher↗

Measuring Massive Black Hole Binaries with LISA

The coalescence of two massive black holes produces gravitational waves (GWs) which can be detected by the space-based detector LISA. By measuring these waves, LISA can determine the various parameters which characterize the source. Measurements of the black hole masses and spins will provide information about the growth of black holes and their host galaxies over time. Measurements of a source's sky position and distance may help astronomers identify an electromagnetic counterpart to the GW event. The counterpart's redshift, combined with the GW-measured luminosity distance, can then be used to measure the Hubble constant and the dark energy parameter $w$. Because the potential science output is so high, it is useful to know in advance how well LISA can measure source parameters for a wide range of binaries. We calculate expected parameter estimation errors using the well-known Fisher matrix method. Our waveform model includes the physics of spin precession, as well as subleading harmonics. When these higher-order effects are not included, strong degeneracies between some parameters cause them to be poorly determined by a GW measurement. When precession and subleading harmonics are properly included, the degeneracies are broken, reducing parameter errors by one to several orders of magnitude.

Lang, Ryan N.↗

Comparison of Multiscale Method of Cells-Based Models for Predicting Elastic Properties of Filament Wound C/C-SiC

Three different multiscale models, based on the method of cells (generalized and high fidelity) micromechanics models were developed and used to predict the elastic properties of C/C-SiC composites. In particular, the following multiscale modeling strategies were employed: Concurrent multiscale modeling of all phases using the generalized method of cells, synergistic (two-way coupling in space) multiscale modeling with the generalized method of cells, and hierarchical (one-way coupling in space) multiscale modeling with the high fidelity generalized method of cells. The three models are validated against data from a hierarchical multiscale finite element model in the literature for a repeating unit cell of C/C-SiC. Furthermore, the multiscale models are used in conjunction with classical lamination theory to predict the stiffness of C/C-SiC plates manufactured via a wet filament winding and liquid silicon infiltration process recently developed by the German Aerospace Institute.

multiscale modeling↗

Critique of Macro Flow/Damage Surface Representations for Metal Matrix Composites Using Micromechanics

Guidance for the formulation of robust, multiaxial, constitutive models for advanced materials is provided by addressing theoretical and experimental issues using micromechanics. The multiaxial response of metal matrix composites, depicted in terms of macro flow/damage surfaces, is predicted at room and elevated temperatures using an analytical micromechanical model that includes viscoplastic matrix response as well as fiber-matrix debonding. Macro flow/damage surfaces (i.e., debonding envelopes, matrix threshold surfaces, macro 'yield' surfaces, surfaces of constant inelastic strain rate, and surfaces of constant dissipation rate) are determined for silicon carbide/titanium in three stress spaces. Residual stresses are shown to offset the centers of the flow/damage surfaces from the origin and their shape is significantly altered by debonding. The results indicate which type of flow/damage surfaces should be characterized and what loadings applied to provide the most meaningful experimental data for guiding theoretical model development and verification.

Lissenden, Cliff J.↗

Krylov vector methods for model reduction and control of flexible structures

Krylov vectors and the concept of parameter matching are combined here to develop model-reduction algorithms for structural dynamics systems. The method is derived for a structural dynamics system described by a second-order matrix differential equation. The reduced models are shown to have a promising application in the control of flexible structures. It can eliminate control and observation spillovers while requiring only the dynamic spillover terms to be considered. A model-order reduction example and a flexible structure control example are provided to show the efficacy of the method.

Su, Tzu-Jeng↗

Modeling stiffness loss in boron/aluminum below the fatigue limit

Boron/aluminum can develop significant internal matrix cracking when fatigued. These matrix cracks can result in a 40 percent secant modulus loss in some laminates, even when fatigued below the fatigue limit. It is shown that the same amount of fatigue damage will develop during stress or strain-controlled tests. Stacking sequence has little influence on secant modulus loss. The secant modulus loss in unidirectional composites is small, whereas the losses are substantial in laminates containing off-axis plies. A simple analysis is presented that predicts unnotched laminate secant modulus loss due to fatigue. The analysis is based upon the elastic modulus and Poisson's ratio of the fiber and matrix, fiber volume fraction, fiber orientations, and the cyclic-hardened yield stress of the matrix material. Excellent agreement was achieved between model predictions and experimental results. With this model, designers can project the material stiffness loss for design load or strain levels and assess the feasibility of its use in stiffness-critical parts.

Johnson, W. S.↗

A fast, reliable algorithm for computing frequency responses of state space models

Computation of frequency responses for large order systems described by time invariant state space systems often provides a bottleneck in control system analysis. It is shown that banding the A-matrix in the state space model can effectively reduce the computation time for such systems while maintaining reliability in the results produced.

Wette, Matt↗

Micromechanical Modeling of High-Strain Thin-Ply Composites

This paper presents a micromechanical analysis for the elastic and viscoelastic behavior of high-strain thin-ply composites. The modeling approach is based on unit cell homogenization. The geometry of the internal weave architecture and ply configuration is characterized via micrographic analysis and explicitly modeled in the unit cell. The composites are modeled as Kirchhoff plates and the homogenization analysis computes the effective relaxation ABD matrix represented by Prony series using the elastic and viscoelastic properties of the constituent fiber and matrix. The formulation of the micromechanical model and numerical implementation d. Composite laminates with 3-ply and 4-ply configurations are studied.

Hamillage, Milinda Yapa↗

Quantum chaos on edge

Recently, the physics of many-body quantum chaotic systems close to their ground states has come under intensified scrutiny. Such studies are motivated by the emergence of model systems exhibiting chaotic fluctuations throughout the entire spectrum [the Sachdev-Ye-Kitaev (SYK) model being a renowned representative] as well as by the physics of holographic principles, which likewise unfold close to ground states. Interpreting the edge of the spectrum as a quantum critical point, here we combine a wide range of analytical and numerical methods to the identification and comprehensive description of two different universality classes: the near edge physics of “sparse” and the near edge of “dense” chaotic systems. The distinction lies in the ratio between the number of a system's random parameters and its Hilbert space dimension, which is exponentially small or algebraically small in the sparse and dense case, respectively. Notable representatives of the two classes are generic chaotic many-body models (sparse) and invariant random matrix ensembles or chaotic gravitational systems (dense). While the two families share identical spectral correlations at energy scales comparable to the level spacing, the density of states and its fluctuations near the edge are different. Considering the SYK model as a representative of the sparse class, we apply a combination of field theory and exact diagonalization to a detailed discussion of its edge spectrum. Conversely, Jackiw-Teitelboim gravity is our reference model for the dense class, where an analysis of the gravitational path integral and random matrix theory reveal universal differences to the sparse class, whose implications for the construction of holographic principles we discuss. Published by the American Physical Society 2024

Altland, Alexander (ORCID:0000000229914805)↗