Engineering Papers⌕ Search

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.

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

At least 973 records · Page 54

Modeling of Damage Initiation and Progression in a SiC/SiC Woven Ceramic Matrix Composite

The goal of an ongoing project at NASA Glenn is to investigate the effects of the complex microstructure of a woven ceramic matrix composite and its variability on the effective properties and the durability of the material. Detailed analysis of these complex microstructures may provide clues for the material scientists who `design the material? or to structural analysts and designers who `design with the material? regarding damage initiation and damage propagation. A model material system, specifically a five-harness satin weave architecture CVI SiC/SiC composite composed of Sylramic-iBN fibers and a SiC matrix, has been analyzed. Specimens of the material were serially sectioned and polished to capture the detailed images of fiber tows, matrix and porosity. Open source analysis tools were used to isolate various constituents and finite elements models were then generated from simplified models of those images. Detailed finite element analyses were performed that examine how the variability in the local microstructure affected the macroscopic behavior as well as the local damage initiation and progression. Results indicate that the locations where damage initiated and propagated is linked to specific microstructural features.

Mital, Subodh K.↗

Improved Subcell Model for the Prediction of Braided Composite Response

In this work, the modeling of triaxially braided composites was explored through a semi-analytical discretization. Four unique subcells, each approximated by a "mosaic" stacking of unidirectional composite plies, were modeled through the use of layered-shell elements within the explicit finite element code LS-DYNA. Two subcell discretizations were investigated: a model explicitly capturing pure matrix regions, and a novel model which absorbed pure matrix pockets into neighboring tow plies. The in-plane stiffness properties of both models, computed using bottom-up micromechanics, correlated well to experimental data. The absorbed matrix model, however, was found to best capture out-of- plane flexural properties by comparing numerical simulations of the out-of-plane displacements from single-ply tension tests to experimental full field data. This strong correlation of out-of-plane characteristics supports the current modeling approach as a viable candidate for future work involving impact simulations.

Cater, Christopher R.↗

Influence of Finite Element Size in Residual Strength Prediction of Composite Structures

The sensitivity of failure load to the element size used in a progressive failure analysis (PFA) of carbon composite center notched laminates is evaluated. The sensitivity study employs a PFA methodology previously developed by the authors consisting of Hashin-Rotem intra-laminar fiber and matrix failure criteria and a complete stress degradation scheme for damage simulation. The approach is implemented with a user defined subroutine in the ABAQUS/Explicit finite element package. The effect of element size near the notch tips on residual strength predictions was assessed for a brittle failure mode with a parametric study that included three laminates of varying material system, thickness and stacking sequence. The study resulted in the selection of an element size of 0.09 in. X 0.09 in., which was later used for predicting crack paths and failure loads in sandwich panels and monolithic laminated panels. Comparison of predicted crack paths and failure loads for these panels agreed well with experimental observations. Additionally, the element size vs. normalized failure load relationship, determined in the parametric study, was used to evaluate strength-scaling factors for three different element sizes. The failure loads predicted with all three element sizes provided converged failure loads with respect to that corresponding with the 0.09 in. X 0.09 in. element size. Though preliminary in nature, the strength-scaling concept has the potential to greatly reduce the computational time required for PFA and can enable the analysis of large scale structural components where failure is dominated by fiber failure in tension.

Satyanarayana, Arunkumar↗

Mineralogy, petrology, and trace element geochemistry of the Johnstown meteorite - A brecciated orthopyroxenite with siderophile and REE-rich components

The compositional and petrologic characteristics of the Johnstown meteorite show it to contain uncontaminated and unbrecciated orthopyroxenite clasts of cumulative origin that (1) must have undergone subsolidus recrystalization, (2) are parental to the brecciated matrix, and (3) show no evidence of a xenolithic, meteoritic contribution to the matrix except for contamination by the projectile which crushed it on impact. The trapped liquid was not introduced in the impact process. The variability of such trace elements as the light rare earth elements, and the presence of plagioclase and olivine in only one of the thin sections studied, demonstrates the heterogeneity of coarse-grained diogenites on a millimeter scale and the difficulty of obtaining representative samples of such meteorites. The data presented indicate that this meteorite is a monominct breccia.

Floran, R. J.↗

A Model for Siderophile Element Distribution in Planetary Differentiation

Planetary differentiation begins with partial melting of small planetesimals. At low degrees of partial melting, a sulfur-rich liquid segregates by physical mechanisms including deformation-assisted porous flow. Experimental studies of the physical mechanisms by which Fe-S melts segregate from the silicate matrix of a molten H chondrite are part of a companion paper. Geochemical studies of these experimental products revealed that metallic liquids were in equilibrium with residual metal in the H chondrite matrix. This contribution explores the geochemical signatures produced by early stages of core formation. Particularly, low-degree partial melt segregation of Fe-S liquids leaves residual metal in the silicate matrix. Some achondrites appear to be residues of partial melting, e.g., ureilites, which are known to contain metal. The metal in these achondrites may show a distinct elemental signature. To quantify the effect of sulfur on siderophile element contents of residual metal we have developed a model based on recent parametrizations of equilibrium solid metal-liquid metal partitioning experiments.

Humayun, M.↗

Origins and Distribution of Chondritic Olivine Inferred from Wild 2 and Chondrite Matrix

To date, only 180 particle impact tracks from Wild 2 have been extracted from the Stardust aerogel collector and even fewer have been thoroughly characterized. In order to provide a cohesive compositional dataset that can be compared to the meteorite record, we have made both major and minor element analyses (TEM/EDXS) of olivine and low-Ca pyroxene for 39 particles harvested from 26 tracks. However, the dearth of equivalent analyses for these phases in chondrite matrix hinders their comparison to the Wild 2 samples. To properly permit comparison of chondritic olivine and pyroxene to the Wild 2 samples, we have also provided a large, comprehensive EPMA dataset (greater than10(exp 3) analyses) of analogous grains (5-30 micrometers) isolated in L/LL3.0-4, CI, CM, CR, CH, CO, CV, Acfer 094, EH3, EL6, and Kakangari matrix

Frank, D. R.↗

Improved Plate and Beam Models for Thermoviscoelastic Constitutive Modeling of Composites

The effective properties of composites are influenced by the time-dependent behavior of polymer matrices very sensitive to changes in temperature. Improved plate and beam models are required to efficiently design, and simulate composite structures when the long-term performance of large anisotropic composite structures is the matter of interest. In this work, mechanics of structure genome (MSG) is used to con-struct linear thermoviscoelastic plate and beam models that can homogenize three-dimensional heterogeneous materials made of constituents with time- and temperature-dependent behavior. The formulation derives the transient strain energy based on integral formulation for thermorheologically simple materials subject to finite temperature changes with the restriction that the strain is small. The reduced time parameter is introduced to relate the time-temperature dependency of the anisotropic material by means of master curves at reference conditions. The new formulation has been implemented in SwiftCompTM, a general-purpose multiscale constitutive modeling code based on MSG. Experimental data and three-dimensional direct numerical simulations of thin-ply high-strain composites (TP-HSC) using a commercial finite element analysis (FEA) package are conducted to verify the accuracy of SwiftCompTM results. The paper also analyzes the relationship between the shift factor of the polymer matrix and the temperature dependencies of the effective beam properties.

Finite element analysis↗

The Cytoskeleton and Force Response Mechanisms

The long term aim of this project was to define the mechanisms by which cells sense and respond to the physical forces experienced at 1g and missing in microgravity. Identification and characterization of the elements of the cells force response mechanism could provide pathways and molecules to serve as targets for pharmacological intervention to mitigate the pathologic effects of microgravity. Mechanical forces experienced by the organism can be transmitted to cells through molecules that allow cells to bind to the extracellular matrix and through other types of molecules which bind cells to each other. These molecules are coupled in large complexes of proteins to structural elements such as the actin cytoskeleton that give the cell the ability to sense, resist and respond to force. Application of small forces to tissue culture cells causes local elevation of intracellular calcium through stretch activated ion channels, increased tyrosine phosphorylation and a restructuring of the actin cytoskeleton. Using collagen coated iron oxide beads and strong magnets, we can apply different levels of force to cells in culture. We have found that force application causes the cells to polymerize actin at the site of mechanical deformation and unexpectedly, to depolymerize actin across the rest of the cell. Observations of GFP- actin expressing cells demonstrate that actin accumulates at the site of deformation within the first five minutes of force application and is maintained for many tens of minutes after force is removed. Consistent with the reinforcement of the cytoskeletal structures underlying the integrin-bead interaction, force also alters the motion of bound magnetic beads. This effect is seen following the removal of the magnetic field, and is only partially ablated by actin disruption with cytochalsin B. While actin is polymerizing locally at the site of force application, force also stimulates a global reduction in actin filament content within the cells. We have examined the roles of several actin filament disassembly factors in the global reduction of cellular actin filaments. The calcium regulated actin filament severing protein gelsolin is not necessary for the increased actin turnover, as cells derived from gelsolin null and wildtype mice still show a reduction in total actin filament content. Instead, our work suggests that the actin binding protein cofilin may be important for these changes in actin dynamics. Cofilin binds to and enhances the disassembly of actin filaments. Using immunological methods, we observe transient changes in the phosphorylation state of cofilin upon force application that suggests that cofilin may mediate actin filament turnover. Early after force application, cofilin is transiently dephosphorylated, activating its actin disassembly activity. Subsequently, we find a hyper-phosphorylation of cofilin, rendering it inactive. This reduction in cofilin activity may explain the stability of the force induced actin structuttes. In testing this hypothesis, we aimed to generate cells that express the constituitively active kinase (LIM-kinase) that phosphorylates cofilin. lnidial attempts in the cell lines used for the our previous studies proved unsuccessful. While we prepare this work for pubication, we are continuing to study other cell lines and tissue sources to determine whether they show a reduction in F-actin content after force application.

Allen, Philip Goodwin↗

Verification and Validation of a Three-Dimensional Composite Impact Model with Tabulated Input

A material model which incorporates several key capabilities which have been identified by the aerospace community as lacking in the composite impact models currently available in the commercial transient dynamic finite element code LS-DYNA® has been developed. The material model utilizes experimentally based tabulated input to define the evolution of plasticity and damage as opposed to specifying discrete input parameters (such as modulus and strength. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. The capability to account for the rate and temperature dependent deformation response of composites has also been incorporated into the material model. For the damage model, a strain equivalent formulation is utilized to allow for the uncoupling of the deformation and damage analyses. In the damage model, a diagonal damage tensor is defined to account for the directionally dependent variation of damage. However, the terms in the damage matrix are semi-coupled such that the damage in a particular coordinate direction is a function of the stresses and plastic strains in all of the coordinate directions. For the failure model, a tabulated approach is utilized in which a stress or strain based invariant is defined as a function of the location of the current stress state in stress space to define the initiation of failure, which allows an arbitrarily shaped failure surface to be defined. A systematic series of validation and verification studies, at a variety of length scales ranging from single element simulations to simulations of a flat panel impact test, have been performed to fully exercise and evaluate the capabilities of the developed model.

Impact↗

Some Experiences with Nonoverlapping Schur Complement Parallel Preconditioning for CFD Calculations

In this work we consider solving matrices which arise from the discretization of advection-diffusion field equations on arbitrary triangulated domains using stabilized numerical methods. The talk will discuss several candidate matrix preconditioning algorithms based on the 2 x 2 block factorization induced by an apriori partitioning of the triangulated domain. Application of the 2 x 2 block preconditioner requires the formation and inversion of the Schur complement submatrix. We consider several strategies for simplifying this task: incomplete Schur complement factorizations, drop tolerance element filling, Schur complement probing, and localized Schur complement inversion. Numerical results will be shown comparing performance and efficiency of these approximations. The matrix preconditioner has also been embedded into a Newton algorithm for solving the nonlinear Euler and Navier-Stokes equations governing compressible flow. The remainder of the talk will show numerous examples in CFD to demonstrate the efficiency and robustness of the techniques.

Barth, Timothy J.↗

Nebular condensation of moderately volatile elements and their abundances in ordinary chondrites

Two models accounting for the concentrations of moderately volatile elements in ordinary chondrites are discussed. The first relies on an apparent correlation between ordinary chondrite/type-I carbonaceous chondrite abundance ratios and ideal solution-condensation temperatures of such elements as Sn, Cu, Li, Na, and Mn to construct a model in which volatiles are lost as gases before nebular condensation or as finely-divided solids that are incompletely agglomerated, and the condensation/agglomeration efficiency declines as a function of time. The second theory largely discounts this gas-dust diffraction effect and instead proposed condensation on matrix and volatization from chondules and coarse metal grains as mechanisms explaining the apparent existence of three distinct classes of volatile-element abundance in ordinary chondrites.

Wai, C. M.↗

Economical Implementation of a Filter Engine in an FPGA

A logic design has been conceived for a field-programmable gate array (FPGA) that would implement a complex system of multiple digital state-space filters. The main innovative aspect of this design lies in providing for reuse of parts of the FPGA hardware to perform different parts of the filter computations at different times, in such a manner as to enable the timely performance of all required computations in the face of limitations on available FPGA hardware resources. The implementation of the digital state-space filter involves matrix vector multiplications, which, in the absence of the present innovation, would ordinarily necessitate some multiplexing of vector elements and/or routing of data flows along multiple paths. The design concept calls for implementing vector registers as shift registers to simplify operand access to multipliers and accumulators, obviating both multiplexing and routing of data along multiple paths. Each vector register would be reused for different parts of a calculation. Outputs would always be drawn from the same register, and inputs would always be loaded into the same register. A simple state machine would control each filter. The output of a given filter would be passed to the next filter, accompanied by a "valid" signal, which would start the state machine of the next filter. Multiple filter modules would share a multiplication/accumulation arithmetic unit. The filter computations would be timed by use of a clock having a frequency high enough, relative to the input and output data rate, to provide enough cycles for matrix and vector arithmetic operations. This design concept could prove beneficial in numerous applications in which digital filters are used and/or vectors are multiplied by coefficient matrices. Examples of such applications include general signal processing, filtering of signals in control systems, processing of geophysical measurements, and medical imaging. For these and other applications, it could be advantageous to combine compact FPGA digital filter implementations with other application-specific logic implementations on single integrated-circuit chips. An FPGA could readily be tailored to implement a variety of filters because the filter coefficients would be loaded into memory at startup.

Kowalski, James E.↗

SPAR reference manual

The functions and operating rules of the SPAR system, which is a group of computer programs used primarily to perform stress, buckling, and vibrational analyses of linear finite element systems, were given. The following subject areas were discussed: basic information, structure definition, format system matrix processors, utility programs, static solutions, stresses, sparse matrix eigensolver, dynamic response, graphics, and substructure processors.

Whetstone, W. D.↗

Design, fabrication and test of graphite/polyimide composite joints and attachments for advanced aerospace vehicles

Bonded and bolted designs are presented for each of four major attachment types. Prepreg processing problems are discussed and quality control data are given for lots 2W4604, 2W4632 and 2W4643. Preliminary design allowables test results for tension tests and compression tests of laminates are included. The final small specimen test matrix is defined and the configuration of symmetric step-lap joint specimens are shown. Finite element modeling studies of a double lap joint were performed to evaluate the number of elements required through the adhesive thickness to assess effects of various joint parameters on stress distributions. Results of finite element analyses assessing the effect of an adhesive fillet on the stress distribution in a double lap joint are examined.

Skoumal, D. E.↗

Appendix E: Research papers. Use of remote sensing in landscape stratification for environmental impact assessment

A refinement to the matrix approach to environmental impact assessment is to use landscape units in place of separate environmental elements in the analysis. Landscape units can be delineated by integrating remotely sensed data and available single-factor data. A remote sensing approach to landscape stratification is described and the conditions under which it is superior to other approaches that require single-factor maps are indicated. Flowcharts show the steps necessary to develop classification criteria, delineate units and a map legend, and use the landscape units in impact assessment. Application of the approach to assessing impacts of a transmission line in Montana is presented to illustrate the method.

Stanturf, J. A.↗

Ply-level failure analysis of a graphite/epoxy laminate under bearing-bypass loading

A combined experimental and analytical study was conducted to investigate and predict the failure modes of a graphite/epoxy laminate subjected to combined bearing and bypass loading. Tests were conducted in a test machine that allowed the bearing-bypass load ratio to be controlled while a single-fastener coupon was loaded to failure in either tension or compression. Onset and ultimate failure modes and strengths were determined for each test case. The damage-onset modes were studied in detail by sectioning and micrographing the damaged specimens. A two-dimensional, finite-element analysis was conducted to determine lamina strains around the bolt hole. Damage onset consisted of matrix cracks, delamination, and fiber failures. Stiffness loss appeared to be caused by fiber failures rather than by matrix cracking and delamination. An unusual offset-compression mode was observed for compressive bearing-bypass laoding in which the specimen failed across its width along a line offset from the hole. The computed lamina strains in the fiber direction were used in a combined analytical and experimental approach to predict bearing-bypass diagrams for damage onset from a few simple tests.

Naik, R. A.↗

Ply-level failure analysis of a graphite/epoxy laminate under bearing-bypass loading

A combined experimental and analytical study was conducted to investigate and predict the failure modes of a graphite/epoxy laminate subjected to combined bearing and bypass loading. Tests were conducted in a test machine that allowed the bearing-bypass load ratio to be controlled while a single-fastener coupon was loaded to failure in either tension or compression. Onset and ultimate failure modes and strengths were determined for each test case. The damage-onset modes were studied in detail by sectioning and micrographing the damaged specimens. A two-dimensional, finite-element analysis was conducted to determine lamina strains around the bolt hole. Damage onset consisted of matrix cracks, delamination, and fiber failures. Stiffness loss appeared to be caused by fiber failures rather than by matrix cracking and delamination. An unusual offset-compression mode was observed for compressive bearing-bypass loading in which the specimen failed across its width along a line offset from the hole. The computed lamina strains in the fiber direction were used in a combined analytical and experimental approach to predict bearing-bypass diagrams for damage onset from a few simple tests.

Naik, R. A.↗

A finite element method for the thermochemical decomposition of polymeric materials. I - Theory

The governing differential equations are developed to model the thermomechanical behavior of chemically decomposing, polymeric materials. These equations account for thermal and gaseous diffusion through a poroelastic, transversely isotropic solid. The Bubnov-Galerkin finite element method is applied to the governing equations to cast the coupled set into a single matrix equation. A method for solving these equations simultaneously at each time step is discussed.

Sullivan, R. M.↗