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 721 records · Page 40

Additional Trials & Tribulations in Synthesizing a Sulfide Standard

Well-characterized sulfide reference materials that can serve as matrix-matched calibrants for in-situ trace element analyses via laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) remain elusive. Here, we describe the creation an in-house sulfide standard at NASA JSC that will be used to measure siderophile and chalcophile trace elements in low pressure experimental products, specifically pentlandite ([FeNi]9S8) and pyrrhotite ([FeNi]1-xS), as well as natural sulfides in terrestrial and meteoritic mineral assemblages. Recent methods in creating a homogenous sulfide standard include pressed pellets, synthesized chips, or fused glass [e.g. 1-3]. Other studies, however, have had success in forming homogenous FeS standards via synthesization methods [e.g. 4-6]. We present a variation of the method described in [4], in which we create a pyrrhotite standard doped with a variety of trace elements (Zn, As, Se, Mo, Ru, Rh, Pd, Sn, Sb, Te, W, Os, Ir, Pt, Au) at ~10-40 ppm and Cu at ~200 ppm. Sulfide compositions were constructed using Fe and Ni metal and elemental sulfur powders. Trace elements were added to elemental sulfur from Atomic Absorption Spectroscopy (AAS) elemental standards as either nitrate or chloride solutions to prevent oxidation of the metal powders. The dried sulfur mixture was mixed with Fe and Ni powders and mechanically mixed before two aliquots were placed in separate SiO2 tubes. Each tube was held under vacuum for ~30 minutes, sealed under vacuum, and then heated at 800C for 48 hours. Like [4], the synthesis products were composed of pourous sulfide crystals. Major element analyses of both experimental aliquots, collected using an Electron Probe Microanalyser (EPMA), yield indistinguishable major element compositions (uncertainties in 2SE & 2RSE[%]), with an average of 57.90  0.09 (0.16 %), 4.95  0.03 (0.55 %), and 39.70  0.14 (0.34 %), for Fe, Ni, and S, respectively. Trace element data were measured using a Photon Machines 193nm laser ablation system coupled to a Thermo-Scientific Element-XR ICP-MS. Spot sizes were limited to 50 μm due to the porous nature of the sulfide target material. Trace element abundances, normalized to Fe as an internal standard, were also found to be homogenous between the two aliquots, with weighted mean 2RSE (%) values of <3.0 for all trace elements. Synthesized products were re-powdered and absolute concentrations measured via solution ICP-MS. Although the sulfide appears to be homogenous, sintering experiments will be performed to more closely match the standard density to natural sulfides and minimize differences in ablation behavior. Improved density also allows for higher sensitivity (i.e. more compact target material) and larger spot sizes or traverses, as void space is eliminated. Additional major and trace element analyses on the products of the sintering experiments will be undertaken.

Jacob B Setera↗

Boundary formulations for sensitivity analysis without matrix derivatives

A new hybrid approach to continuum structural shape sensitivity analysis employing boundary element analysis (BEA) is presented. The approach uses iterative reanalysis to obviate the need to factor perturbed matrices in the determination of surface displacement and traction sensitivities via a univariate perturbation/finite difference (UPFD) step. The UPFD approach makes it possible to immediately reuse existing subroutines for computation of BEA matrix coefficients in the design sensitivity analysis process. The reanalysis technique computes economical response of univariately perturbed models without factoring perturbed matrices. The approach provides substantial computational economy without the burden of a large-scale reprogramming effort.

Kane, J. H.↗

Evaluation of thermal and mechanical loading effects on the structural behavior of a SIC/titanium composite

Composite specimens of titanium-15-3 matrix reinforced with continuous SCS-6 silicon carbide fibers were tested under a variety of thermal and mechanical loadings. A combined experimental/finite element approach was used to estimate the effective in situ modulus of the matrix material and to evaluate changes in modulus due to the applied loads. Several fiber orientations were tested. Results indicate that the effect of thermal loads on composite stiffness varies with fiber orientation. Applications of this method to test specimens damaged by uniaxial tension, thermal cycling, and isothermal fatigue loadings are used to illustrate that by monitoring overall structural behavior, changes in stiffness caused by thermomechanical loading can be detected.

Grady, J. E.↗

Exploiting symmetries in the modeling and analysis of tires

A computational procedure is presented for reducing the size of the analysis models of tires having unsymmetric material, geometry and/or loading. The two key elements of the procedure when applied to anisotropic tires are: (1) decomposition of the stiffness matrix into the sum of an orthotropic and nonorthotropic parts; and (2) successive application of the finite-element method and the classical Rayleigh-Ritz technique. The finite-element method is first used to generate few global approximation vectors (or modes). Then the amplitudes of these modes are computed by using the Rayleigh-Ritz technique. The proposed technique has high potential for handling practical tire problems with anisotropic materials, unsymmetric imperfections and asymmetric loading. It is also particularly useful for use with three-dimensional finite-element models of tires.

Noor, Ahmed K.↗

A three-dimensional nonlinear Timoshenko beam based on the core-congruential formulation

A three-dimensional, geometrically nonlinear two-node Timoshenkoo beam element based on the total Larangrian description is derived. The element behavior is assumed to be linear elastic, but no restrictions are placed on magnitude of finite rotations. The resulting element has twelve degrees of freedom: six translational components and six rotational-vector components. The formulation uses the Green-Lagrange strains and second Piola-Kirchhoff stresses as energy-conjugate variables and accounts for the bending-stretching and bending-torsional coupling effects without special provisions. The core-congruential formulation (CCF) is used to derived the discrete equations in a staged manner. Core equations involving the internal force vector and tangent stiffness matrix are developed at the particle level. A sequence of matrix transformations carries these equations to beam cross-sections and finally to the element nodal degrees of freedom. The choice of finite rotation measure is made in the next-to-last transformation stage, and the choice of over-the-element interpolation in the last one. The tangent stiffness matrix is found to retain symmetry if the rotational vector is chosen to measure finite rotations. An extensive set of numerical examples is presented to test and validate the present element.

Crivelli, Luis A.↗

Basaltic volcanism on the eucrite parent body - Petrology and chemistry of the polymict eucrite ALHA80102

The polymict eucrite meteorite ALHA80102 is an unequilibrated breccia of basaltic and gabbroic clasts in a fragmental matrix. Clasts include basalts of many textural types, cumulate gabbro, black 'glass', and ferroan troctolite (plagioclase, silica, Fe-rich olivine, ilmenite, mesostasis). Ferroan troctolite has not been previously reported from eucrites or howardites; it is interpreted as the end-product of fractional crystallization of eucritic magmas. Bulk and trace element compositions (by electron microprobe and INAA) of clasts and matrix from ALHA80102 are similar to those of other eucrites; the meteorite contains clasts similar to Juvinas and to Stannern. A clast of cumulate eucrite gabbro is enriched in the light rare earths (La/Lu = 2XCI). This clast is interpreted as an unrepresentative sample of metamorphically equilibrated gabbro; LREE-enriched magmas need not be invoked. ALHA80102 is similar to other polymict eucrites from the Allan Hills and may be paired with ALHA76005, ALHA77302, and ALHA78040.

Treiman, A. H.↗

Parallel triangularization of substructured finite element problems

Much of the computational effort of the finite element process involves the solution of a system of linear equations. The coefficient matrix of this system, known as the global stiffness matrix, is symmetric, positive definite, and generally sparse. An important technique for reducing the time required to solve this system is substructuring or matrix partitioning. Substructuring is based on the idea of dividing a structure into pieces, each of which can then be analyzed relatively indepenently. As a result of this division, each point in the finite element discretization is either interior to a substructure or on a boundary between substructures. Contributions to the global stiffness matrix from connections between boundary points from the K(bb) matrix are reported. The triangularization of a general K(bb) matrix on a parallel machine is specifically discussed.

Leuze, M. R.↗

A kinematically driven anisotropic viscoelastic constitutive model applied to tires

Aircraft tires are composite structures manufactured with viscoelastic materials such as carbon black filled rubber and nylon cords. When loaded they experience large deflections and moderately large strains. Detailed structural models of tires require the use of either nonlinear shell or nonlinear three dimensional solid finite elements. Computational predictions of the dynamic response of tires must consider the composite viscoelastic material behavior in a realistic fashion. We describe a modification to a nonlinear anisotropic shell finite element so it can be used to model viscoelastic stresses during general deformations. The model is developed by introducing internal variables of the type used to model elastic strain energy. The internal variables are strains, curvatures, and transverse shear angles which are in a one-to-one correspondence with the generalized coordinates used to model the elastic strain energy for nonlinear response. A difference-relaxation equation is used to relate changes in the observable strain field to changes in the internal strain field. The internal stress state is introduced into the equilibrium equations by converting it to nodal loads associated with the element's displacement degrees of freedom. In this form the tangent matrix in the Newton-Raphson solution algorithm is not modified from its form for the nonlinear statics problem. Only the gradient vector is modified and the modification is not computationally costly. The existing finite element model for the Space Shuttle nose gear tire is used to provide examples of the algorithm. In the first example, the tire's rim is displaced at a constant rate up to a fixed value. In the second example, the tire's rim is enforced to follow a saw tooth load and unload curve to generate hysteresis loops.

Johnson, Arthur R.↗

Stratospheric aerosols - Undissolved granules and physical state

The physical state of stratospheric aerosol particles was studied along with the nature of included undissolved granules. It was found that: (1) undissolved granules are present in only a third of the aerosol particles and are not always associated with a slurry matrix, (2) those undissolved particles usually contain only sulfur and sodium or undetectable light elements, (3) all stratospheric aerosols have a similar appearance: a volatile slurry mixture of crystalline-like material in a liquid matrix, and (4) variations in the relative amounts of liquid and crystalline materials at different times cause variations in particle spreading.

Farlow, N. H.↗

Stress Stiffening Of STARDYNE Plate Elements

Stress-stiffening program, written in FORTRAN language, calculates geometric stiffness matrix, resulting from stressed condition. Creates file appended to STARDYNE bulk data deck containing entries added to STARDYNE stiffness matrix. Format of this file compatible with STARDYNE.

Philpot, Dennis C.↗

On the use of the consistent mass matrix in flexible multibody systems

Two different methods are used in finite element structural analysis to define element mass distribution: consistent and lumped mass approximations. The lumped mass method is currently favored in multibody dynamics software packages because of its computational efficiency. This paper, however, presents a multibody dynamics formulation based on the consistent mass approach which should prove to be as computationally efficient as the lumped mass approach. It is shown that when component flexibility is modelled using modal coordinates, the equations of motion, as formulated by these two methods, are very similar in appearance, thereby facilitating their comparison. It is also shown, through simulation of a simple multibody system, that a system formulated with the lumped mass approximation will often require more finite elements and smaller integration step sizes than the corresponding consistent mass system.

Anthony, Tobin C.↗

Advanced reinforcement systems for intermetallic applications

A 2-D axisymmetric model was employed to determine the magnitude of the radial, axial, and hoop stresses caused by the thermal expansion difference between fiber and matrix and which result from the fabrication temperature cycle. Finite element analysis was conducted for single fiber model systems based on SCS-6/Ti3Al+Nb and Al2O3/NiAl. The stress distribution due to the imposition of a graded intermediate layer for each system was determined and included variables of layer thickness and gradation in interlayer chemistry in order to vary the expansion gradient between fiber and matrix. Thermal cycling tests were conducted on sputter coated SCS-6 fibers selectively coated with Ti3Al+Nb, with and without an intermediate layer. Cracking of the Ti3Al+Nb layers was prevented by an interlayer based on Ti-TiN-Ti. The interlayer thickness appeared critical to its efficiency. Similarly, for the case of Al2O3/NiAl, an intermediate layer consisting of a Ni bond coat on the sapphire fiber followed by a graded Al2O3-NiAl layer did not crack when given a thermal excursion to 1100 C and then cooled to room temperature. Acoustic emission tests on single fiber specimens were unsuccessful in detecting load drops associated with the successive fracture of the fiber. For the SCS-6/Ti3Al system this was the result of several factors which included the matrix/fiber ratio and poor bonding of the matrix and fiber. In the case of the Al2O3/NiAl system brittle failure of the NiAl matrix precluded fiber breakdown during tensile loading.

Merrick, Howard F.↗

Scanning and Transmission Electron Microscopy of High Temperature Materials

Software and hardware updates to further extend the capability of the electron microscope were carried out. A range of materials such as intermetallics, metal-matrix composites, ceramic-matrix composites, ceramics and intermetallic compounds, based on refractory elements were examined under this research. Crystal structure, size, shape and volume fraction distribution of various phases which constitute the microstructures were examined. Deformed materials were studied to understand the effect of interfacial microstructure on the deformation and fracture behavior of these materials. Specimens tested for a range of mechanical property requirements, such as stress rupture, creep, low cycle fatigue, high cycle fatigue, thermomechanical fatigue, etc. were examined. Microstructural and microchemical stability of these materials exposed to simulated operating environments were investigated. The EOIM Shuttle post-flight samples were also examined to understand the influence of low gravity processing on microstructure. In addition, fractographic analyses of Nb-Zr-W, titanium aluminide, molybdenum silicide and silicon carbide samples were carried out. Extensive characterization of sapphire fibers in the fiber-reinforced composites made by powder cloth processing was made. Finally, pressure infiltration casting of metal-matrix composites was carried out.

Source record↗

Hamilton's principle - Finite-element methods and flexible body dynamics

A variational formulation is given for the equations of motion for an unconstrained elastic body, and Hamilton's principle is used to derive the equations of motion and deformation of the body. Finite-element approximations are developed for these dynamical equations with respect to a body axis system satisfying the mean axis conditions. The free-body influence matrix for the body then is developed in terms of the finite-element model parameters.

Cavin, R. K., III↗

Evaluation of thermal and mechanical loading effects on the structural behavior of a SiC/titanium composite

Composite specimens of titanium-15-3 matrix reinforced with continuous SCS-6 silicon carbide fibers were tested under a variety of thermal and mechanical loadings. A combined experimental/finite element approach was used to estimate the effective in-situ modulus of the matrix material and to evaluate changes in modulus due to the applied loads. Several fiber orientations were tested. Results indicate that the effect of the thermal loads on composite stiffness varies with fiber orientation. Applications of this method to test specimens damaged by uniaxial tension, thermal cycling, and isothermal fatigue loadings are used to illustrate that by monitoring overall structural behavior, changes in stiffness caused by thermomechanical loading can be detected.

Grady, Joseph E.↗

Finite element mesh refinement criteria for stress analysis

This paper discusses procedures for finite-element mesh selection and refinement. The objective is to improve accuracy. The procedures are based on (1) the minimization of the stiffness matrix race (optimizing node location); (2) the use of h-version refinement (rezoning, element size reduction, and increasing the number of elements); and (3) the use of p-version refinement (increasing the order of polynomial approximation of the elements). A step-by-step procedure of mesh selection, improvement, and refinement is presented. The criteria for 'goodness' of a mesh are based on strain energy, displacement, and stress values at selected critical points of a structure. An analysis of an aircraft lug problem is presented as an example.

Kittur, Madan G.↗

Error analysis for semi-analytic displacement derivatives with respect to shape and sizing variables

Sensitivity analysis is fundamental to the solution of structural optimization problems. Consequently, much research has focused on the efficient computation of static displacement derivatives. As originally developed, these methods relied on analytical representations for the derivatives of the structural stiffness matrix (K) with respect to the design variables (b sub i). To extend these methods for use with complex finite element formulations and facilitate their implementation into structural optimization programs using the general finite element method analysis codes, the semi-analytic method was developed. In this method the matrix the derivative of K/the derivative b sub i is approximated by finite difference. Although it is well known that the accuracy of the semi-analytic method is dependent on the finite difference parameter, recent work has suggested that more fundamental inaccuracies exist in the method when used for shape optimization. Another study has argued qualitatively that these errors are related to nonuniform errors in the stiffness matrix derivatives. The accuracy of the semi-analytic method is investigated. A general framework was developed for the error analysis and then it is shown analytically that the errors in the method are entirely accounted for by errors in delta K/delta b sub i. Furthermore, it is demonstrated that acceptable accuracy in the derivatives can be obtained through careful selection of the finite difference parameter.

Fenyes, Peter A.↗