SEARCH · Engineering Papers

Results for “Matrix elements”

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 577 records · Page 32

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Stability analysis of bridged cracks in brittle matrix composites

The bridging of matrix cracks by fibers is an important toughening mechanism in fiber reinforced brittle matrix composites. This paper presents the results of a nonlinear finite element analysis of the Mode-I propagation of a bridged matrix crack in a finite size specimen. The composite is modeled as an orthotropic continuum and the bridging due to the fibers is modeled as a distribution of tractions which resist crack opening. A critical stress intensity factor criterion is employed for matrix crack propagation while a critical crack opening condition is used for fiber failure. The structural response of the specimen (load-deflection curves) as well as the stress intensity factor of the propagating crack are calculated for various constituent properties and specimen configurations for both tensile and bending loading. By controlling the length of the bridged crack results are obtained which highlight the transition from stable to unstable behavior of the propagating crack.

Ballarini, Roberto↗

NASA NTRS2014

Theoretical Development of an Orthotropic Elasto-Plastic Generalized Composite Material Model

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites is becoming critical as these materials are gaining increased usage in the aerospace and automotive industries. While there are several composite material models currently available within LS-DYNA (Registered), there are several features that have been identified that could improve the predictive capability of a composite model. To address these needs, a combined plasticity and damage model suitable for use with both solid and shell elements is being developed and is being implemented into LS-DYNA as MAT_213. A key feature of the improved material model is the use of tabulated stress-strain data in a variety of coordinate directions to fully define the stress-strain response of the material. To date, the model development efforts have focused on creating the plasticity portion of the model. The Tsai-Wu composite failure model has been generalized and extended to a strain-hardening based orthotropic material model with a non-associative flow rule. The coefficients of the yield function, and the stresses to be used in both the yield function and the flow rule, are computed based on the input stress-strain curves using the effective plastic strain as the tracking variable. The coefficients in the flow rule are computed based on the obtained stress-strain data. The developed material model is suitable for implementation within LS-DYNA for use in analyzing the nonlinear response of polymer composites.

Plasticity↗

NASA NTRS2024

Bulk Major and Trace Elemental Composition of an Aggregate Sample From Asteroid Bennu

On September 24, 2023, NASA’s OSIRIS-REx mission returned a capsule to Earth carrying material from asteroid Bennu. This event was the first time a U.S. mission delivered pristine samples from an asteroid and is the largest asteroid sample return to date. As these samples represent some of the oldest, most primitive, and pristine materials available to us, and which originate from a known and well-studied asteroid, they allow us a rare opportunity to gain a better understanding of the formation and evolution of our solar system. Key to understanding the material returned from asteroid Bennu is establishing its bulk chemical composition. Previous studies have shown that each chondrite group has a distinct elemental composition. For the carbonaceous chondrites specifically, each group exhibits a distinct pattern of moderately and highly volatile elemental depletions, relative to CI chondrites. CI chondrites are considered the most primitive chondrite group and broadly represent the solar photosphere composition. The two most striking features of these depletion patterns is that the moderately volatile element depletions increase with decreasing 50% condensation temperature and then plateau out at abundances that roughly correlate with matrix abundance. Due to these distinctive patterns, bulk elemental composition has become an important classification tool for establishing the different chondrite groups and the connections between them. As such, determining the bulk chemical composition of the Bennu aggregates will help to test two mission hypothesis: “Bennu’s bulk elemental composition reflects that of its main parent asteroid and is similar to the composition of the Sun, with depletions in moderately to highly volatile elements” and “Bennu’s dominant lithologies are comparable in bulk mineralogy, petrology, and composition to the most aqueously altered carbonaceous chondrites” . Furthermore as the carbonaceous and non-carbonaceous chondrites are thought to have likely formed in the inner and outer protoplanetary disk, respectively (e.g., [8]), determining the bulk chemical compositions of the Bennu aggregates will also help test the major mission hypothesis that “Bennu's parent body formed beyond the snow line by accretion of material in the pro-toplanetary disk” . Regarding chondrite formation, observed elemental patterns have been successfully used to model how the mixing of volatile-rich and volatile-poor chondritic components can produce the observed carbonaceous chondrites groups. They have also been significant in investigating how volatilization processes influenced chondrite formation. As such, establishing bulk chemical compositions of the pristine Bennu samples is vital for understanding the asteroid, and solar system formation. To begin this processes, we analyzed the bulk major and trace elemental compositions of Bennu aggregates.

P Koefoed↗

NASA NTRS1994

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↗

NASA NTRS1974

Matrix isolation technique for the study of some factors affecting the partitioning of trace elements

The factors that affect the preferred positions of cations in ionic solid solutions were investigated utilizing vibrational spectroscopy. Solid solutions of the sulfate and chromate ions codoped with La(+3) and Ca(+2) in a KBr host lattice were examined as a function of the polyvalent cation concentration. The cation-anion pairing process was found to be random for Ca(+2), whereas the formation of La(+3)-SO4(-2) ion pairs with a C2 sub v bonding geometry is highly preferential to any type of La(+3)-CrO4(-2) ion pair formation. The relative populations of ion pair site configurations are discussed in terms of an energy-entropy competition model which can be applied to the partition of trace elements during magmatic processes.

Grzybowski, J. M.↗

NASA NTRS1982

Micromechanical predictions of crack initiation, propagation and crack growth resistance in boron/aluminum composites

An elastoplastic, axisymmetric finite element model has been used to predict the initiation and propagation of a crack in a composite model consisting of a single broken boron fiber embedded in an annular sheath of aluminum matrix. The accuracy of the axisymmetric finite element model for crack problems has been established by solving the classical problem of a penny-shaped crack in a thick cylindrical rod under axial tension. Also, the stress intensity factors predicted by the present numerical model are compared with continuum results. A constant displacement boundary condition applied during an increment of crack growth permits a substantial amount of stable crack growth in the matrix material. The concept of Crack Growth Resistance Curves (KR-curves) has been used to determine the point of crack instability

Mahishi, J. M.↗

NASA NTRS1978

A high order accurate finite element algorithm for high Reynolds number flow prediction

A Galerkin-weighted residuals formulation is employed to establish an implicit finite element solution algorithm for generally nonlinear initial-boundary value problems. Solution accuracy, and convergence rate with discretization refinement, are quantized in several error norms, by a systematic study of numerical solutions to several nonlinear parabolic and a hyperbolic partial differential equation characteristic of the equations governing fluid flows. Solutions are generated using selective linear, quadratic and cubic basis functions. Richardson extrapolation is employed to generate a higher-order accurate solution to facilitate isolation of truncation error in all norms. Extension of the mathematical theory underlying accuracy and convergence concepts for linear elliptic equations is predicted for equations characteristic of laminar and turbulent fluid flows at nonmodest Reynolds number. The nondiagonal initial-value matrix structure introduced by the finite element theory is determined intrinsic to improved solution accuracy and convergence. A factored Jacobian iteration algorithm is derived and evaluated to yield a consequential reduction in both computer storage and execution CPU requirements while retaining solution accuracy.

Baker, A. J.↗

NASA NTRS1982

Calculation of exact vibration modes for plane grillages by the dynamic stiffness method

A dynamic stiffness method is developed for the calculation of the exact modal parameters for plane grillages which consist of straight and uniform beams with coincident elastic and inertial axes. Elementary bending-torsion beam theory is utilized, and bending translation is restricted to one direction. The exact bending-torsion dynamic stiffness matrix is obtained for a straight and uniform beam element with coincident elastic and inertial axes. The element stiffness matrices are assembled using the standard procedure of the static stiffness method to form the dynamic stiffness matrix of the complete grillage. The exact natural frequencies, mode shapes, and generalized masses of the grillage are then calculated by solving a nonlinear eigenvalue problem based on the dynamic stiffness matrix. The exact modal solutions for an example grillage are calculated and compared with the approximate solutions obtained by using the finite element method.

Hallauer, W. L., Jr.↗

NASA NTRS1982

Aspects of the history of 66095 based on trace elements in clasts and whole rock

Halogens, P, U and Na are reported in anorthositic and basaltic clasts and matrix from rusty rock 66095. Large fractions of Cl and Br associated with the separated phases from 66095 are soluble in H2O. Up to two orders of magnitude variation in concentrations of these elements in the breccia components and varying H2O-soluble Cl/Br ratios indicate different sources of volatiles. An approximately constant ratio of the H2O- to 0.1 M HNO3-soluble Br in the various components suggests no appreciable alteration in the original distributions of this element in the breccia forming processes. Up to 50% or more of the phosphorus and of the non-H2O-soluble Cl was dissolved from most of the breccia components by 0.1 M HNO3. Clast and matrix residues from the leaching steps contain, in most cases, the Cl/P2O5 ratio found in 66095 whole rock and in a number of other Apollo 16 samples. Evidence that phosphates are the major P-phases in the breccia is based on the 0.1 M acid solubility of Cl and P in the matrix sample and on elemental concentrations which are consistent with those of KREEP.

Jovanovic, S.↗

NASA NTRS2015

Verification and Validation of a Three-Dimensional Generalized Composite Material Model

A general purpose orthotropic elasto-plastic computational constitutive material model has been developed to improve predictions of the response of composites subjected to high velocity impact. The three-dimensional orthotropic elasto-plastic composite material model is being implemented initially for solid elements in LS-DYNA as MAT213. In order to accurately represent the response of a composite, experimental stress-strain curves are utilized as input, allowing for a more general material model that can be used on a variety of composite applications. The theoretical details are discussed in a companion paper. This paper documents the implementation, verification and qualitative validation of the material model using the T800-F3900 fiber/resin composite material

Plasticity↗

NASA NTRS1993

Proposed framework for thermomechanical life modeling of metal matrix composites

The framework of a mechanics of materials model is proposed for thermomechanical fatigue (TMF) life prediction of unidirectional, continuous-fiber metal matrix composites (MMC's). Axially loaded MMC test samples are analyzed as structural components whose fatigue lives are governed by local stress-strain conditions resulting from combined interactions of the matrix, interfacial layer, and fiber constituents. The metallic matrix is identified as the vehicle for tracking fatigue crack initiation and propagation. The proposed framework has three major elements. First, TMF flow and failure characteristics of in situ matrix material are approximated from tests of unreinforced matrix material, and matrix TMF life prediction equations are numerically calibrated. The macrocrack initiation fatigue life of the matrix material is divided into microcrack initiation and microcrack propagation phases. Second, the influencing factors created by the presence of fibers and interfaces are analyzed, characterized, and documented in equation form. Some of the influences act on the microcrack initiation portion of the matrix fatigue life, others on the microcrack propagation life, while some affect both. Influencing factors include coefficient of thermal expansion mismatch strains, residual (mean) stresses, multiaxial stress states, off-axis fibers, internal stress concentrations, multiple initiation sites, nonuniform fiber spacing, fiber debonding, interfacial layers and cracking, fractured fibers, fiber deflections of crack fronts, fiber bridging of matrix cracks, and internal oxidation along internal interfaces. Equations exist for some, but not all, of the currently identified influencing factors. The third element is the inclusion of overriding influences such as maximum tensile strain limits of brittle fibers that could cause local fractures and ensuing catastrophic failure of surrounding matrix material. Some experimental data exist for assessing the plausibility of the proposed framework.

Halford, Gary R.↗

NASA NTRS2013

The Effect of Scale Dependent Discretization on the Progressive Failure of Composite Materials Using Multiscale Analyses

A multiscale modeling methodology, which incorporates a statistical distribution of fiber strengths into coupled micromechanics/ finite element analyses, is applied to unidirectional polymer matrix composites (PMCs) to analyze the effect of mesh discretization both at the micro- and macroscales on the predicted ultimate tensile (UTS) strength and failure behavior. The NASA code FEAMAC and the ABAQUS finite element solver were used to analyze the progressive failure of a PMC tensile specimen that initiates at the repeating unit cell (RUC) level. Three different finite element mesh densities were employed and each coupled with an appropriate RUC. Multiple simulations were performed in order to assess the effect of a statistical distribution of fiber strengths on the bulk composite failure and predicted strength. The coupled effects of both the micro- and macroscale discretizations were found to have a noticeable effect on the predicted UTS and computational efficiency of the simulations.

Ricks, Trenton M.↗

NASA NTRS1984

Algorithms for solving large sparse systems of simultaneous linear equations on vector processors

Very efficient algorithms for solving large sparse systems of simultaneous linear equations have been developed for serial processing computers. These involve a reordering of matrix rows and columns in order to obtain a near triangular pattern of nonzero elements. Then an LU factorization is developed to represent the matrix inverse in terms of a sequence of elementary Gaussian eliminations, or pivots. In this paper it is shown how these algorithms are adapted for efficient implementation on vector processors. Results obtained on the CYBER 200 Model 205 are presented for a series of large test problems which show the comparative advantages of the triangularization and vector processing algorithms.

David, R. E.↗

NASA NTRS2016

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↗

NASA NTRS1989

Vibration transmission through rolling element bearings in geared rotor system, part 1

A mathematical model is proposed to examine the vibration transmission through rolling element bearings in geared rotor systems. Current bearing models, based on either ideal boundary conditions for the shaft or purely translational stiffness element description, cannot explain how the vibratory motion may be transmitted from the rotating shaft to the casing. This study clarifies this issue qualitatively and quantitatively by developing a comprehensive bearing stiffness matrix of dimension 6 model for the precision rolling element bearings from basic principles. The proposed bearing formulation is extended to analyze the overall geared rotor system dynamics including casing and mounts. The bearing stiffness matrix is included in discrete system models using lumped parameter and/or dynamic finite element techniques. Eigensolution and forced harmonic response due to rotating mass unbalance or kinematic transmission error excitation for a number of examples are computed.

Singh, Rajendra↗

NASA NTRS1991

Computational simulation of microfracture in high temperature metal matrix composites

A computational simulation procedure, using three-dimensional finite element method and global strain energy release rates, is described to predict the microfracture process and identify/quantify the hierarchy of respective fracture modes in metal matrix composites. The procedure is used to predict the microfracture in unidirectional composites under longitudinal loads and accounts for the interphase strength in the in situ state. A novel procedure to computationally simulate the fiber pushthrough process is also described. In this simulation, the interface material is replaced by an anisotropic material with greatly reduced shear modulus in order to simulate the fiber pushthrough process using a linear analysis. Such a procedure is easily implemented and is computationally very effective.

Mital, Subodh K.↗

NASA NTRS1991

Experimental and analytical analysis of stress-strain behavior in a (90/0 deg)2s, SiC/Ti-15-3 laminate

The nonlinear stress strain behavior of 90 degree/0 degree sub 2s, SiC/Ti-15-3 composite laminate was numerically investigated with a finite element, unit cell approach. Tensile stress-strain curves from room temperature experiments depicted three distinct regions of deformation, and these regions were predicted by finite element analysis. The first region of behavior, which was linear elastic, occurred at low applied stresses. As applied stresses increased, fiber/matrix debonding in the 90 degree plies caused a break in the stress-strain curve and initiated a second linear region. In this second region, matrix plasticity in the 90 degree plies developed. The third region, which was typified by nonlinear, stress-strain behavior occr red at high stresses. In this region, the onset of matrix plasticity in the 0 degree plies stiffened the laminate in the direction transverse to the applied load. Metallographic sections confirmed the existence of matrix plasticity in specific areas of the structure. Finite element analysis also predicted these locations of matrix slip.

Lerch, Bradley A.↗

NASA NTRS1989

A three-node C deg element for analysis of laminated composite sandwich shells

A three-node flat shell element with C deg rotation fields has been developed for analysis of arbitrary composite shells. The element may consist of any number of orthotropic layers, each layer having different material properties and angular orientation. The formulation includes coupling between bending and extension, which is essential for analysis of unsymmetric laminates. Shearing deflections are included, since laminated and sandwich construction frequently results in shear stiffness much smaller than bending stiffness. Formulation of the element is straightforward, and calculation of its stiffness matrix is simple and fast. Convergence of solutions with mesh refinement is uniform for both thin and thick shells and is insensitive to element shape, although not as rapid as some other elements that lack one or more capabilities of the newly developed element. An experimental verification of the shall element is reported in the appendix.

Martin, C. Wayne↗