Search NASA⌕ Search

SEARCH · Search NASA

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 469 records · Page 26

Humans to Mars, but How Many? Trade Space Analyses for Mars Mission Crew Size

Missions to Mars will differ from previous human spaceflight missions in that the crew of astronauts will be required to operate in an Earth-independent manner due to the long communication delays. Without a systematic, repeatable process to determine the number and composition of crew necessary to successfully accomplish these missions, NASA increases the risk that crew sizes may be too small to meet primary mission objectives under nominal conditions and, more consequentially, that crewmembers may not have the expertise needed to successfully respond to unforeseen failures without the real-time expertise of the Mission Control Central (MCC) team NASA currently relies upon. The NASA Engineering and Safety Center (NESC) developed a methodology for assessing the trade space of factors that affect the number of crew for Mars missions based on human performance and expertise modeling. This methodology includes the consideration of results from three human performance models developed using the Improved Performance Research and Integration Tool (IMPRINT) modeling platform as well as a custom-built model on expertise trained within the crew based on NASA’s crew qualification and responsibilities matrix (CQRM). We present example model results and discuss implications on trade space analyses for Mars mission crew size.

Mars Mission↗

Power flow for spacecraft power systems

A method for constructing the generalized system-level admittance matrix for use with a Newton-Raphson power flow is presented. The network modeling technique presented does not use the standard pi-equivalent models, which assume a lossless return path, for the transmission line and transformer. If the return path cannot be assumed lossless, then the standard algorithms for constructing the system admittance matrix cannot be used. The method presented here uses concepts from linear graph theory to combine network modules to form the system-level admittance matrix. The modeling technique is presented, and the resulting matrix is used with a standard Newton-Raphson power flow to calculate all system voltages and current (power) flows.

Halpin, S. M.↗

An improved compression molding technology for continuous fiber reinforced composite laminate. Part 1: AS-4/LaRC-TPI 1500 (HFG) Prepreg system

Poor processability of fiber reinforced high performance polyimide thermoplastic resin composites is a well recognized issue which, in many cases, prohibits the fabrication of composite parts with satisfactorily consolidated quality. Without modifying the resin matrix chemistry, improved compression modeling procedures were proposed and investigated with the AS-4/LaRC-TPI 1500 High Flow Grade (HFG) prepreg system. Composite panels with excellent C-scans can be consistently molded by this method under 700 F and a consolidation pressure as low as 100 psi. A mechanism for the consolidation of the composite under this improved molding technique is discussed. This mechanism reveals that a certain degree of matrix shear and tow filament slippage and nesting between plies occur during consolidation, which leads to a reduction of the consolidating pressure necessary to offset the otherwise intimate inter fiber-fiber contact and consequently achieves a better consolidation quality. Outstanding short beam shear strength and flexural strength were obtained from the molded panels. A prolonged consolidation step under low pressure, i.e., 100 psi at 700 F for 75 minutes, was found to significantly enhance the composite mechanical properties.

Hou, Tan-Hung↗

Linearized T-Matrix and Mie Scattering Computations

We present a new linearization of T-Matrix and Mie computations for light scattering by non-spherical and spherical particles, respectively. In addition to the usual extinction and scattering cross-sections and the scattering matrix outputs, the linearized models will generate analytical derivatives of these optical properties with respect to the real and imaginary parts of the particle refractive index, and (for non-spherical scatterers) with respect to the ''shape'' parameter (the spheroid aspect ratio, cylinder diameter/height ratio, Chebyshev particle deformation factor). These derivatives are based on the essential linearity of Maxwell's theory. Analytical derivatives are also available for polydisperse particle size distribution parameters such as the mode radius. The T-matrix formulation is based on the NASA Goddard Institute for Space Studies FORTRAN 77 code developed in the 1990s. The linearized scattering codes presented here are in FORTRAN 90 and will be made publicly available.

Analytic linearization↗

Investigation of a Macromechanical Approach to Analyzing Triaxially-Braided Polymer Composites

A macro level finite element-based model has been developed to simulate the mechanical and impact response of triaxially-braided polymer matrix composites. In the analytical model, the triaxial braid architecture is simulated by using four parallel shell elements, each of which is modeled as a laminated composite. The commercial transient dynamic finite element code LS-DYNA is used to conduct the simulations, and a continuum damage mechanics model internal to LS-DYNA is used as the material constitutive model. The material stiffness and strength values required for the constitutive model are determined based on coupon level tests on the braided composite. Simulations of quasi-static coupon tests of a representative braided composite are conducted. Varying the strength values that are input to the material model is found to have a significant influence on the effective material response predicted by the finite element analysis, sometimes in ways that at first glance appear non-intuitive. A parametric study involving the input strength parameters provides guidance on how the analysis model can be improved.

Goldberg, Robert K.↗

Modification of a Macromechanical Finite-Element Based Model for Impact Analysis of Triaxially-Braided Composites

A macro level finite element-based model has been developed to simulate the mechanical and impact response of triaxially-braided polymer matrix composites. In the analytical model, the triaxial braid architecture is simulated by using four parallel shell elements, each of which is modeled as a laminated composite. For the current analytical approach, each shell element is considered to be a smeared homogeneous material. The commercial transient dynamic finite element code LS-DYNA is used to conduct the simulations, and a continuum damage mechanics model internal to LS-DYNA is used as the material constitutive model. The constitutive model requires stiffness and strength properties of an equivalent unidirectional composite. Simplified micromechanics methods are used to determine the equivalent stiffness properties, and results from coupon level tests on the braided composite are utilized to back out the required strength properties. Simulations of quasi-static coupon tests of several representative braided composites are conducted to demonstrate the correlation of the model. Impact simulations of a represented braided composites are conducted to demonstrate the capability of the model to predict the penetration velocity and damage patterns obtained experimentally.

Goldberg, Robert K.↗

Stochastic process approximation for recursive estimation with guaranteed bound on the error covariance

An approach, is proposed for the design of approximate, fixed order, discrete time realizations of stochastic processes from the output covariance over a finite time interval, was proposed. No restrictive assumptions are imposed on the process; it can be nonstationary and lead to a high dimension realization. Classes of fixed order models are defined, having the joint covariance matrix of the combined vector of the outputs in the interval of definition greater or equal than the process covariance; (the difference matrix is nonnegative definite). The design is achieved by minimizing, in one of those classes, a measure of the approximation between the model and the process evaluated by the trace of the difference of the respective covariance matrices. Models belonging to these classes have the notable property that, under the same measurement system and estimator structure, the output estimation error covariance matrix computed on the model is an upper bound of the corresponding covariance on the real process. An application of the approach is illustrated by the modeling of random meteorological wind profiles from the statistical analysis of historical data.

Menga, G.↗

Composites boost 21st-century aircraft engines

Research and development in light-weight, high-temperature composite materials for ultrahigh-bypass engines to be used in high-speed civil transport/rotocraft is presented. It is noted that the expected benefits to be attained by this R&D include weight reduction, lowered fuel consumption, and lower direct operating costs. A major effort underway in this area is the Advanced High Temperature Engine Materials Technology Program (HITEMP) of NASA, which focuses on providing revolutionary high-temperature composite materials: to 425 C (800 F) for polymer-matrix composites (PMCs), to 1250 C (2280 F) for metal-matrix/intermetallic-matrix composites (MMCs/IMCs), and to as high as 1650 C (3000 F) for ceramic-matrix composites (CMCs). Analytical modeling is being used to investigate the structural behavior of these advanced materials in six distinct areas: micromechanics, deformation and damage, fatigue, fracture, trade-off studies, and load definition. It is concluded that the development of advanced materials such as high-temperature composites is highly dependent on the availability of high-temperature fibers. The wide range of fiber characteristics needed will require the development of more than one fiber. In general, a candidate fiber should have low density, high strength, high stiffness, a CTE matching the matrix, chemical compatibility with the matrix, environmental stability and appropriate fiber diameter.

Stephens, Joseph R.↗

The upper bounds of reduced axial and shear moduli in cross-ply laminates with matrix cracks

The present study proposes a mathematical model utilizing the internal state variable concept for predicting the upper bounds of the reduced axial and shear stiffnesses in cross-ply laminates with matrix cracks. The displacement components at the matrix crack surfaces are explicitly expressed in terms of the observable axial and shear strains and the undamaged material properties. The reduced axial and shear stiffnesses are predicted for glass/epoxy and graphite/epoxy laminates. Comparison of the model with other theoretical and experimental studies is also presented to confirm direct applicability of the model to angle-ply laminates with matrix cracks subjected to general in-plane loading.

Lee, Jong-Won↗

Applying transfer matrix method to the estimation of the modal characteristics of the NASA Mini-Mass Truss

It is beneficial to use a distributed parameter model for large space structures because the approach minimizes the number of model parameters. Holzer's transfer matrix method provides a useful means to simplify and standardize the procedure for solving the system of partial differential equations. Any large space structures can be broken down into sub-structures with simple elastic and dynamical properties. For each single element, such as beam, tether, or rigid body, we can derive the corresponding transfer matrix. Combining these elements' matrices enables the solution of the global system equations. The characteristics equation can then be formed by satisfying the appropriate boundary conditions. Then natural frequencies and mode shapes can be determined by searching the roots of the characteristic equation at frequencies within the range of interest. This paper applies this methodology, and the maximum likelihood estimation method, to refine the modal characteristics of the NASA Mini-Mast Truss by successively matching the theoretical response to the test data of the truss. The method is being applied to more complex configurations.

Shen, Ji-Yao↗

Development of a Subcell Based Modeling Approach for Modeling the Architecturally Dependent Impact Response of Triaxially Braided Polymer Matrix Composites

Understanding the high velocity impact response of polymer matrix composites with complex architectures is critical to many aerospace applications, including engine fan blade containment systems where the structure must be able to completely contain fan blades in the event of a blade-out. Despite the benefits offered by these materials, the complex nature of textile composites presents a significant challenge for the prediction of deformation and damage under both quasi-static and impact loading conditions. The relatively large mesoscale repeating unit cell (in comparison to the size of structural components) causes the material to behave like a structure rather than a homogeneous material. Impact experiments conducted at NASA Glenn Research Center have shown the damage patterns to be a function of the underlying material architecture. Traditional computational techniques that involve modeling these materials using smeared homogeneous, orthotropic material properties at the macroscale result in simulated damage patterns that are a function of the structural geometry, but not the material architecture. In order to preserve heterogeneity at the highest length scale in a robust yet computationally efficient manner, and capture the architecturally dependent damage patterns, a previously-developed subcell modeling approach where the braided composite unit cell is approximated as a series of four adjacent laminated composites is utilized. This work discusses the implementation of the subcell methodology into the commercial transient dynamic finite element code LS-DYNA (Livermore Software Technology Corp.). Verification and validation studies are also presented, including simulation of the tensile response of straight-sided and notched quasi-static coupons composed of a T700/PR520 triaxially braided [0deg/60deg/-60deg] composite. Based on the results of the verification and validation studies, advantages and limitations of the methodology as well as plans for future work are discussed.

finite element method↗

Residual Stresses Induced due to Curing of the Bulk Matrix in a Simplified Three-Dimensional (3D) Woven Repeating Unit Cell

It has been observed that 3D woven materials are highly susceptible to processes induced defects including tow misalignment, intratow and intertow cracking and/or voids. Adequate process modeling of 3D woven Polymer Matrix Composites (PMCs) is necessary to predict and estimate the effects of the manufacturing process on these defects. A preliminary step towards this is achieved with the commercial Finite Element Analysis (FEA) tool Abaqus complemented with written user subroutines to account for the effect of shrinkage and thermoelastic properties evolution as a function of the degree of cure of the bulk matrix phase in a 3D woven composite repeating unit cell (RUC). As a pathfinder, the process modeling framework is demonstrated on a simplified 3D woven RUC of an AS4/RIMR 135 system. Curing simulations are performed on the 3D woven RUC to observe the processing induced residual stresses. The results of the analysis show that high stress concentrations, present within the 3D intertow matrix bordering the binder tow, coincide with X-Ray computed tomography (X-Ray-CT) scan data depicting cracking in the same location

3D Woven↗

Modeling of Time-Dependent Strength Degradation of SiC/SiC Ceramic Matrix Composites via MAC/GMC Computer Code

Silicon carbide fiber reinforced silicon carbide (SiC/SiC) ceramic matrix composites (CMCs) display time-dependent strength degradation at intermediate temperatures (600 to 900 degrees Centigrade). This is generally believed to be an oxidation induced phenomenon. The understanding of the effect of temperature with environment (oxidation) is key towards development of SiC/SiC CMCs with a reliable load carrying capacity. Various theories have been proposed to explain the strength degradation. One suggests that the boron nitride (BN) coating deposited on the fibers oxidizes causing fusion of fibers. Another theory proposes that the SiC fibers are oxidized forming a silica scale leading to premature fiber failure. A more recent theory suggests that SiC fiber strength is intrinsically time-dependent due to slow crack growth of flaws in the fibers. An empirical model, termed as a "fiber classic model," which is based on a standard slow crack growth type power-law, has been implemented within NASA's micromechanics-based MAC/GMC computer code as a user routine. Model parameters for this "classic model" were calibrated from stress-rupture data of Hi-Nicalon (TM) monofilaments using the maximum likelihood estimation (MLE) technique. This new capability in the MAC/GMC computer code was then used to predict the stress-rupture behavior of Hi-Nicalon (TM) tows as well as 2-D SiC/SiC composites reinforced with Hi-Nicalon (TM) fibers. Results demonstrate that the MAC/GMC with this new capability successfully predicts the time-to-failure vs. applied stress within the intermediate temperature range at various scales as well as laminated composites in an oxidizing environment.

Environmental degradation↗

Micromechanical Modeling of the Thermal Expansion of Graphite/copper Composites with Nonuniform Microstructure

Two micromechanical models were developed to investigate the thermal expansion of graphite/copper (Gr/Cu) composites. The models incorporate the effects of temperature-dependent material properties, matrix inelasticity, initial residual stresses due to processing history, and nonuniform fiber distribution. The first model is based on the multiple concentric cylinder geometry, with each cylinder treated as a two-phase composite with a characteristic fiber volume fractions. By altering the fiber volume fraction of the individual cylinders, unidirectional composites with radially nonuniform fiber distributions can be investigated using this model. The second model is based on the inelastic lamination theory. By varying the fiber content in the individual laminae, composites with nonuniform fiber distribution in the thickness direction can be investigated. In both models, the properties of the individual regions (cylinders or laminae) are calculated using the method of cells micromechanical model. Classical incremental plasticity theory is used to model the inelastic response of the copper matrix at the microlevel. The models were used to characterize the effects of nonuniform fiber distribution on the thermal expansion of Gr/Cu. These effects were compared to the effects of matrix plasticity, choice of stress-free temperature, and slight fiber misalignment. It was found that the radially nonuniform fiber distribution has little effect on the thermal expansion of Gr/Cu but could become significant for composites with large fiber-matrix transverse CTE and Young's modulus mismatch. The effect of nonuniform fiber distribution in the through-thickness direction of a laminate was more significant, but only approached that of the stress-free temperature for the most extreme cases that include large amounts of bending. Subsequent comparison with experimental thermal expansion data indicated the need for more accurate characterization of the graphite fiber thermomechanical properties. Correlation with cyclic data revealed the presence of a mechanism not considered in the developed models. The predicted response did, however, exhibit ratcheting behavior that has been observed experimentally in Gr/Cu. Finally, simulation of the actual fiber distribution of particular specimens had little effect on the predicted thermal expansion.

Bednarcyk, Brett A.↗

Probability Bounds Analysis Applied to Multi-Purpose Crew Vehicle Nonlinearity

The Multi-Purpose Crew Vehicle (MPCV) Program Orion vehicle finite element model (FEM) was updated based on a modal test performed by Lockheed Martin. Due to nonlinearity observed in the test results, linear low force level (LL) and high force level (HL) FEMs were developed for use during various Space Launch System (SLS) flight regimes depending on expected forcing levels. Uncertainty models were derived for the combined MPCV and MPCV Stage Adaptor LL and HL Hurty/Craig-Bampton (HCB) components based on the MPCV structural test article Configuration 4 modal test-analysis correlation results. Subsequently, system-level uncertainty quantification analyses were performed using both models for various SLS flight configurations to determine the impact of the nonlinearity on important system metrics. The system metrics included both transfer functions associated with attitude control and dynamic loads associated with aerodynamic buffeting during ascent. In each case, an independent Monte Carlo (MC) analysis was performed, and no attempt was made to combine the results. The Hybrid Parametric Variation (HPV) method was used to develop the LL and HL MPCV HCB uncertainty models. The HPV method provides both parametric and non-parametric components of uncertainty. The non-parametric uncertainty accounts for the difference in model-form between the linearized analytical model and the corresponding linearized component test results in the form of mode shapes and frequencies at that force level. This linear model-form uncertainty is implemented in the HPV method using random matrix theory. However, the HPV uncertainty models developed for the linear LL and HL MPCV components do not account for the nonlinearity in the MPCV. With respect to the linearized models, this nonlinearity is also an uncertainty in model form, but in this case, it must be treated independently as an epistemic uncertainty. It represents a lack of knowledge, in contrast to an aleatory uncertainty due to the randomness of a variable. In the case of an epistemic variable, the true value is unknown, only the interval within which it lies is known. Epistemic uncertainty can be reduced with increased knowledge, while in general, aleatory uncertainty cannot. This work combines the epistemic uncertainty due to the MPCV nonlinearity with the parametric and non-parametric uncertainty within the HPV method using a second order propagation approach. The LL and HL test data is augmented with surrogate test data derived from a nonlinear MPCV representation. The impact of the MPCV nonlinearity on system response statistics is determined using a series of cumulative distribution functions in the form of a horsetail plot, or p-box. This results in an interval of probabilities for a specific response value, or an interval of response values at a specific probability.

Daniel C Kammer↗

Statistical scaling relationships and size effects in the strength and creep rupture of fibrous composites

In this presentation we discuss a new theoretical model and supporting experimental results for the strength and lifetime in creep rupture of unidirectional, carbon fiber/epoxy matrix composites at ambient conditions. First we review the 'standard' Weibull/power-law methodology that has been standard practice. Then we discuss features of a recent model which build on the statistical aspects of fiber strength, micromechanical aspects of stress transfer around fiber breaks, and time-dependent creep of the matrix. The model is applied to 'microcomposites' consisting of seven fibers in a matrix for which strength and creep-rupture data are available. The model yields Weibull distributions in an envelope format for both strength and lifetime. The respective shape, scale and power-law parameters depend on such parameters as the Weibull shape parameter for fiber strength, the exponent for matrix creep, the effective load transfer length (which grows in time due to matrix creep) and the critical cluster size for failed fibers. The experimental results are consistent with the theory, though time-dependent debonding appears to be part of the failure process.

Phoenix, S. Leigh↗

Elastic and Piezoelectric Properties of Boron Nitride Nanotube Composites: Finite Element Model - Part II

This paper is the second part of a two-part series where the first part presents a molecular dynamics model of a single Boron Nitride Nanotube (BNNT) and this paper scales up to multiple BNNTs in a polymer matrix. This paper presents finite element (FE) models to investigate the effective elastic and piezoelectric properties of (BNNT) nanocomposites. The nanocomposites studied in this paper are thin films of polymer matrix with aligned co-planar BNNTs. The FE modelling approach provides a computationally efficient way to gain an understanding of the material properties. We examine several FE models to identify the most suitable models and investigate the effective properties with respect to the BNNT volume fraction and the number of nanotube walls. The FE models are constructed to represent aligned and randomly distributed BNNTs in a matrix of resin using 2D and 3D hollow and 3D filled cylinders. The homogenisation approach is employed to determine the overall elastic and piezoelectric constants for a range of volume fractions. These models are compared with an analytical model based on Mori-Tanaka formulation suitable for finite length cylindrical inclusions. The model applies to primarily single-wall BNNTs but is also extended to multi-wall BNNTs, for which preliminary results will be presented. Results from the Part 1 of this series can help to establish a constitutive relationship for input into the finite element model to enable the modeling of multiple BNNTs in a polymer matrix.

Kim, H. Alicia↗

Ground-state properties of the periodic Anderson model

The ground-state energy, hybridization matrix element, local moment, and spin-density correlations of a one-dimensional, finite-chain, periodic, symmetric Anderson model are obtained by numerical simulations and compared with perturbation theory and strong-coupling results. It is found that the local f-electron spins are compensated by correlation with other f-electrons as well as band electrons leading to a nonmagnetic ground state.

Blankenbecler, R.↗