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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2020

A multi-physics constitutive model to predict hydrolytic aging in quasi-static behaviour of thin cross-linked polymers

The effect of hydrolytic aging on mechanical quasi-static responses of rubber-like materials, in particular, the idealized Mullins effect and permanent set have been modeled. The effect of hydrolytic damage on the mechanical integrity of the polymer matrix is modeled as the direct competition of two micro-structural phenomena (i) chain scission and (ii) reduction of cross-links. Both phenomena and their correlation were modeled and thus, the strain energy of the polymer matrix is written with respect to three independent mechanisms; i) the shrinking original matrix that has not been attacked by water, ii) conversion of the first network to a new network due to the reduction of the cross-links, and iii) energy loss from network degradation due to water attacks to polymer active agents. The proposed model satisfies the Clausius-Duhem inequality and is thus physically feasible. The model is validated with respect to sets of our experimental data and other sets available in the literature. The proposed model is based on the assumption of homogeneous diffusion and mainly relevant for thin samples. In view of its accuracy, interpret-ability, and deep insight it provides into the nature of damage accumulation, the model is a good choice for further implementation in FE applications.

42 ENGINEERING↗

NASA NTRS2010

Microscopic Scale Simulation of the Ablation of Fibrous Materials

This slide presentation reviews the ablation by oxidation of carbon-fiber preforms impregnated in carbonized phenolic-formaldehyde matrix is modeled at microscopic scale. Direct numerical simulations show that the matrix ablates in volume leaving the carbon fibers exposed. This is due to the fact that the reactivity of carbonized phenolics is higher than the reactivity of carbon fibers. After the matrix is depleted, the fibers ablate showing progressive reduction of their diameter. The overall material recession occurs when the fibers are consumed. Two materials with the same carbon-fiber preform, density and chemical composition, but with different matrix distributions are studied. These studies show that at moderate temperatures (< 1000K) the microstructure of the material influences its recession rate; a fact that is not captured by current models that are based on chemical composition only. Surprisingly, the response of these impregnated-fiber materials is weakly dependent on the microstructure at very high temperatures (e.g., Stardust peak heating conditions: 3360K).

Lachaud, Jean↗

NASA NTRS1990

A creep model for metallic composites based on matrix testing: Application to Kanthal composites

An anisotropic creep model is formulated for metallic composites with strong fibers and low to moderate fiber volume percent (less than 40 percent). The idealization admits no creep in the local fiber direction and assumes equal creep strength in longitudinal and transverse shear. Identification of the matrix behavior with that of the isotropic limit of the theory permits characterization of the composite through uniaxial creep tests on the matrix material. Constant and step-wise creep tests are required as a data base. The model provides an upper bound on the transverse creep strength of a composite having strong fibers embedded in a particular matrix material. Comparison of the measured transverse strength with the upper bound gives an assessment of the integrity of the composite. Application is made to a Kanthal composite, a model high-temperature composite system. Predictions are made of the creep response of fiber reinforced Kanthal tubes under interior pressure.

Binienda, W. K.↗

DOE OSTI2022

Inverse Calculation of Burden Distribution Matrix Using B-spline Model Based PDF control in Blast Furnace Burden Charging Process

The inverse calculation of burden distribution matrix (BDM) is one of the most important challenges in the blast furnace operation in iron-making processes. In general, blast furnace consumes 65% of the total energy for the whole steel-making. Focusing on this practical challenge, this article proposes a new burden distribution spatial model in calculating burden charging process, and develops a B-spline approximation-based probability density function (PDF) control algorithm to assign the expected thickness distribution of burden layer and, thus, develops a new method for the required inverse calculation of BDM. First, a novel method for the thickness distribution of burden layer is given using B-spline model to produce an expected distribution shape subjected to a desired tracking within a specific spatial constraint. Then, according to the coexistence of continuous and bounded discrete variables in BDM, a novel hybrid optimization control method by combining integer programming and PDF tracking is further established for the effective inverse calculation of BDM. Finally, the proposed PDF-based iterative inverse calculation of BDM using B-spline models are tested using various data from industrial examples. Furthermore, the simulation results show that the proposed method is well suited to solve the BDM inverse calculation problem in practice.

42 ENGINEERING↗

NASA NTRS2005

Life Modeling and Design Analysis for Ceramic Matrix Composite Materials

The primary research efforts focused on characterizing and modeling static failure, environmental durability, and creep-rupture behavior of two classes of ceramic matrix composites (CMC), silicon carbide fibers in a silicon carbide matrix (SiC/SiC) and carbon fibers in a silicon carbide matrix (C/SiC). An engineering life prediction model (Probabilistic Residual Strength model) has been developed specifically for CMCs. The model uses residual strength as the damage metric for evaluating remaining life and is posed probabilistically in order to account for the stochastic nature of the material s response. In support of the modeling effort, extensive testing of C/SiC in partial pressures of oxygen has been performed. This includes creep testing, tensile testing, half life and residual tensile strength testing. C/SiC is proposed for airframe and propulsion applications in advanced reusable launch vehicles. Figures 1 and 2 illustrate the models predictive capabilities as well as the manner in which experimental tests are being selected in such a manner as to ensure sufficient data is available to aid in model validation.

Source record↗

NASA NTRS1978

Residual-flexibility corrections for transient modal rotordynamic models

A modal residual-flexibility approach due to Schwendler and MacNeal (1962) is adapted to account for the 'static' contribution of higher-frequency modes without requiring their integration. It is assumed that each rotor mode acts as a lightly damped second-order system. The additional accuracy provided by residual-flexibility corrections becomes progressively more important as a modal model's actual boundary conditions are forced to deviate from the boundary conditions used to define the rotor's original structural model (stiffness matrix) and its associated eigendata input to the transient modal model. An analysis of the high-pressure-oxygen turbopump of the Space Shuttle main engine shows that the residual-flexibility approach ensures a substantial improvement in accuracy for a relatively moderate increase in computer-time requirements.

Childs, D. W.↗

NASA NTRS2000

A Melting Layer Model for Passive/Active Microwave Remote Sensing Applications: Model Formulation and Comparison with Observations - Part 1

In this study, a 1-D steady-state microphysical model which describes the vertical distribution of melting precipitation particles is developed. The model is driven by the ice-phase precipitation distributions just above the freezing level at applicable gridpoints of "parent" 3-D cloud-resolving model (CRM) simulations. It extends these simulations by providing the number density and meltwater fraction of each particle in finely separated size categories through the melting layer. The depth of the modeled melting layer is primarily determined by the initial material density of the ice-phase precipitation. The radiative properties of melting precipitation at microwave frequencies are calculated based upon different methods for describing the dielectric properties of mixed phase particles. Particle absorption and scattering efficiencies at the Tropical Rainfall Measuring Mission Microwave Imager frequencies (10.65 to 85.5 GHz) are enhanced greatly for relatively small (approx. 0.1) meltwater fractions. The relatively large number of partially-melted particles just below the freezing level in stratiform regions leads to significant microwave absorption, well-exceeding the absorption by rain at the base of the melting layer. Calculated precipitation backscatter efficiencies at the Precipitation Radar frequency (13.8 GHz) increase in proportion to the particle meltwater fraction, leading to a "bright-band" of enhanced radar reflectivities in agreement with previous studies. The radiative properties of the melting layer are determined by the choice of dielectric models and the initial water contents and material densities of the "seeding" ice-phase precipitation particles. Simulated melting layer profiles based upon snow described by the Fabry-Szyrmer core-shell dielectric model and graupel described by the Maxwell-Garnett water matrix dielectric model lead to reasonable agreement with radar-derived melting layer optical depth distributions. Moreover, control profiles that do not contain mixed-phase precipitation particles yield optical depths that are systematically lower than those observed. Therefore, the use of the melting layer model to extend 3-D CRM simulations appears justified, at least until more realistic spectral methods for describing melting precipitation in high-resolution, 3-D CRM's are implemented.

Olson, William S.↗

NASA NTRS1996

Neutron diffraction measurements and modeling of residual strains in metal matrix composites

Neutron diffraction measurements at room temperature are used to characterize the residual strains in tungsten fiber-reinforced copper matrix, tungsten fiber-reinforced Kanthal matrix, and diamond particulate-reinforced copper matrix composites. Results of finite element modeling are compared with the neutron diffraction data. In tungsten/Kanthal composites, the fibers are in compression, the matrix is in tension, and the thermal residual strains are a strong function of the volume fraction of fibers. In copper matrix composites, the matrix is in tension and the stresses are independent of the volume fraction of tungsten fibers or diamond particles and the assumed stress free temperature because of the low yield strength of the matrix phase.

Saigal, A.↗

NASA NTRS1993

Convection equation modeling: A non-iterative direct matrix solution algorithm for use with SINDA

The determination of the boundary conditions for a component-level analysis, applying discrete finite element and finite difference modeling techniques often requires an analysis of complex coupled phenomenon that cannot be described algebraically. For example, an analysis of the temperature field of a coldplate surface with an integral fluid loop requires a solution to the parabolic heat equation and also requires the boundary conditions that describe the local fluid temperature. However, the local fluid temperature is described by a convection equation that can only be solved with the knowledge of the locally-coupled coldplate temperatures. Generally speaking, it is not computationally efficient, and sometimes, not even possible to perform a direct, coupled phenomenon analysis of the component-level and boundary condition models within a single analysis code. An alternative is to perform a disjoint analysis, but transmit the necessary information between models during the simulation to provide an indirect coupling. For this approach to be effective, the component-level model retains full detail while the boundary condition model is simplified to provide a fast, first-order prediction of the phenomenon in question. Specifically for the present study, the coldplate structure is analyzed with a discrete, numerical model (SINDA) while the fluid loop convection equation is analyzed with a discrete, analytical model (direct matrix solution). This indirect coupling allows a satisfactory prediction of the boundary condition, while not subjugating the overall computational efficiency of the component-level analysis. In the present study a discussion of the complete analysis of the derivation and direct matrix solution algorithm of the convection equation is presented. Discretization is analyzed and discussed to extend of solution accuracy, stability and computation speed. Case studies considering a pulsed and harmonic inlet disturbance to the fluid loop are analyzed to assist in the discussion of numerical dissipation and accuracy. In addition, the issues of code melding or integration with standard class solvers such as SINDA are discussed to advise the user of the potential problems to be encountered.

Schrage, Dean S.↗

NASA NTRS1990

Local stresses in metal matrix composites subjected to thermal and mechanical loading

An elasticity solution has been used to analyze matrix stresses near the fiber/matrix interface in continuous fiber-reinforced metal-matrix composites, modeling the micromechanics in question in terms of a cylindrical fiber and cylindrical matrix sheath which is embedded in an orthotropic medium representing the composite. The model's predictions for lamina thermal and mechanical properties are applied to a laminate analysis determining ply-level stresses due to thermomechanical loading. A comparison is made between these results, which assume cylindrical symmetry, and the predictions yielded by a FEM model in which the fibers are arranged in a square array.

Highsmith, Alton L.↗

NASA NTRS1996

Computer Simulation of Spatial Arrangement and Connectivity of Particles in Three-Dimensional Microstructure: Application to Model Electrical Conductivity of Polymer Matrix Composite

Computer simulation is a powerful tool for analyzing the geometry of three-dimensional microstructure. A computer simulation model is developed to represent the three-dimensional microstructure of a two-phase particulate composite where particles may be in contact with one another but do not overlap significantly. The model is used to quantify the "connectedness" of the particulate phase of a polymer matrix composite containing hollow carbon particles in a dielectric polymer resin matrix. The simulations are utilized to estimate the morphological percolation volume fraction for electrical conduction, and the effective volume fraction of the particles that actually take part in the electrical conduction. The calculated values of the effective volume fraction are used as an input for a self-consistent physical model for electrical conductivity. The predicted values of electrical conductivity are in very good agreement with the corresponding experimental data on a series of specimens having different particulate volume fraction.

Louis, P.↗

NASA NTRS2010

Microscopic Scale Simulation of the Ablation of Fibrous Materials

Ablation by oxidation of carbon-fiber preforms impregnated in carbonized phenolic matrix is modeled at microscopic scale. Direct numerical simulations show that the carbonized phenolic matrix ablates in volume leaving the carbon fibers exposed. This is due to the fact that the reactivity of carbonized phenolic is higher than the reactivity of carbon fibers. After the matrix is depleted, the fibers ablate showing progressive reduction of their diameter. The overall material recession occurs when the fibers are consumed. Two materials with the same carbon-fiber preform, density and chemical composition, but with different matrix distributions are studied. These studies show that at moderate temperatures (< 1000 K) the microstructure of the material influences its recession rate; a fact that is not captured by current models that are based on chemical composition only. Surprisingly, the response of these impregnated-fiber materials is weakly dependent on the microstructure at very high temperatures (e.g., Stardust peak heating conditions: 3360K).

Lachaud, Jean Romain↗

NASA NTRS1985

Models of the upper atmosphere

Several computer models of the thermosphere are discussed. J70MM, the Jacchia 1970 model with matrix and matrix mean output, and J703X, the Jacchia 1970, 1971, and 1977 models, are of primary interest. The subroutines in these programs were studied in detail, and several mistakes were found and corrected. It is proposed that a weighted average of the three-hour geomagnetic index be used in the models rather than a single index value. The densities (for a given date and time) generated by the J703X program (all three models), using both a single value of the index and weighted averages for three different lag times, are tabulated and discussed. Most of the equations used in the models are the results of empirical curve fits. An attempt was made to generate a theoretical prediction of the thermospheric temperature profile based on a simple physical model of atmospheric heat conduction in a spherically-symmetric shell. The exponential heating theory and the temperature dependence of k are discussed.

Davis, T. M.↗

NASA NTRS2006

Modeling the Stress Strain Behavior of Woven Ceramic Matrix Composites

Woven SiC fiber reinforced SiC matrix composites represent one of the most mature composite systems to date. Future components fabricated out of these woven ceramic matrix composites are expected to vary in shape, curvature, architecture, and thickness. The design of future components using woven ceramic matrix composites necessitates a modeling approach that can account for these variations which are physically controlled by local constituent contents and architecture. Research over the years supported primarily by NASA Glenn Research Center has led to the development of simple mechanistic-based models that can describe the entire stress-strain curve for composite systems fabricated with chemical vapor infiltrated matrices and melt-infiltrated matrices for a wide range of constituent content and architecture. Several examples will be presented that demonstrate the approach to modeling which incorporates a thorough understanding of the stress-dependent matrix cracking properties of the composite system.

Morscher, Gregory N.↗

NASA NTRS1991

Analysis of structural response data using discrete modal filters

The application of reciprocal modal vectors to the analysis of structural response data is described. Reciprocal modal vectors are constructed using an existing experimental modal model and an existing frequency response matrix of a structure, and can be assembled into a matrix that effectively transforms the data from the physical space to a modal space within a particular frequency range. In other words, the weighting matrix necessary for modal vector orthogonality (typically the mass matrix) is contained within the reciprocal model matrix. The underlying goal of this work is mostly directed toward observing the modal state responses in the presence of unknown, possibly closed loop forcing functions, thus having an impact on both operating data analysis techniques and independent modal space control techniques. This study investigates the behavior of reciprocol modal vectors as modal filters with respect to certain calculation parameters and their performance with perturbed system frequency response data.

Freudinger, Lawrence C.↗

NASA NTRS1993

Identification of a second-order mechanical system model from a state space realization

Many system identification algorithms produce models in modern state variable form. When a mechanical system is identified one knows that there must be a model in second order form. In structural dynamics one is interested in either a modal model of the system, or in a model including a mass matrix, stiffness matrix and a damping matrix. In this paper, algorithms are developed that convert a modern state space realization into the above two representations of interest in mechanical systems. The algorithm can be used with any identification methods that produce a modern state variable representation such as the Eigensystem Realization Algorithm (ERA), its modified version by data correlation (ERA/DC), and combinations of these with Observer/Kalman Filter Identification (OKID). An algorithm is developed that allows one to identify the damping matrix in a model representation, so that one determines how near to modal the damping is in the system, and can therefore understand the degree to which the modes are coupled by the damping in the system. A second algorithm produced the mass, damping, and stiffness matrices, given the input and output matrices. It is assumed that the number of sensors (or the number of actuators) is greater than or equal to the number of modes in the system, and it is proved that this is a necessary condition to be able to uniquely identify these matrices. Experience with examples suggests that a singular value truncation involved in these algorithms can help one determine the true system order. This truncation has the benefit of knowledge of the correct form for a mechanical system model, which is not present in the modern control identification algorithm.

Tseng, Dong-Huei↗

NASA NTRS2003

Failure Criteria for FRP Laminates in Plane Stress

A new set of six failure criteria for fiber reinforced polymer laminates is described. Derived from Dvorak's fracture mechanics analyses of cracked plies and from Puck's action plane concept, the physically-based criteria, denoted LaRC03, predict matrix and fiber failure accurately without requiring curve-fitting parameters. For matrix failure under transverse compression, the fracture plane is calculated by maximizing the Mohr-Coulomb effective stresses. A criterion for fiber kinking is obtained by calculating the fiber misalignment under load, and applying the matrix failure criterion in the coordinate frame of the misalignment. Fracture mechanics models of matrix cracks are used to develop a criterion for matrix in tension and to calculate the associated in-situ strengths. The LaRC03 criteria are applied to a few examples to predict failure load envelopes and to predict the failure mode for each region of the envelope. The analysis results are compared to the predictions using other available failure criteria and with experimental results. Predictions obtained with LaRC03 correlate well with the experimental results.

Davila, Carlos G.↗

DOE OSTI2024

Comparing models for a unitary black hole S matrix

This paper compares features, challenges, and puzzles of different models for a unitary black hole S -matrix, focusing on both recent nonisometric models, as well as “nonviolent unitarization,” which is based on new quantum interactions of a black hole. As a foundation for comparison, the description of real-time Hawking evolution is first overviewed, including leading effects of gravitational dressing and backreaction. Connection is then made to qubit models for evolution, and some technology is outlined to facilitate their description. Important features of both nonisometric models and nonviolent unitarization are investigated in qubit models, which illustrate essential differences between the respective approaches. The nonisometric models present puzzles for understanding evolution of internal outgoing excitations, which can be excited by interactions such as particle decay. Qubit models for nonviolent unitarization are further developed, and nicely illustrate aspects of that approach. Some remaining questions in generalizing to more complete models for evolution are discussed. Published by the American Physical Society 2024

Astronomy & Astrophysics↗