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 325 records · Page 18

Transcriptional profile of a myotube starvation model of atrophy

Skeletal muscle wasting is a pervasive phenomenon that can result from a wide range of pathological conditions as well as from habitual muscular inactivity. The present work describes a cell-culture condition that induces significant atrophy in skeletal muscle C2C12 myotubes. The failure to replenish differentiation media in mature myotubes leads to rapid atrophy (53% in diameter), which is referred to here as starvation. Affymetrix microarrays were used to develop a transcriptional profile of control (fed) vs. atrophied (nonfed) myotubes. Myotube starvation was characterized by an upregulation of genes involved in translational inhibition, amino acid biosynthesis and transport, and cell cycle arrest/apoptosis, among others. Downregulated genes included several structural and regulatory elements of the extracellular matrix as well as several elements of Wnt/frizzled and TGF-beta signaling pathways. Interestingly, the characteristic transcriptional upregulation of the ubiquitin-proteasome system, calpains, and cathepsins known to occur in multiple in vivo models of atrophy were not seen during myotube starvation. With the exception of the downregulation of extracellular matrix genes, serine protease inhibitor genes, and the upregulation of the translation initiation factor PHAS-I, this model of atrophy in cell culture has a transcriptional profile quite distinct from any study published to date with atrophy in whole muscle. These data show that, although the gross morphology of atrophied muscle fibers may be similar in whole muscle vs. myotube culture, the processes by which this phenotype is achieved differ markedly.

NASA Discipline Musculoskeletal↗

FEM R-Matrix Calculations of Electron-H2 Collisions

Electron-H2 Collisions have been studied using the finite-element R-matrix method. Our approach uses a mixed finite-element and Gaussian basis to describe the continuum electron. In the inner region, the multi-centered nature of the Gaussian basis provides an efficient representation in the regions of space near the nuclei, whereas the piecemeal and energy-independent nature of the FEM basis is particular well suited for R-matrix calculations. Fixed nuclei calculations have been carried out at 1.4 a(sub 0), the equilibrium internuclear distance of the ground state. Up to six target states are included: the X(sup 1)E+(sub g), b(sup 3)E+(sub u), a(sup 3)E+(sub g), B(sup 1)E+(sub u), c(sup 3)II(sub u), and C(sup 1)II(sub u) states, and configuration-interaction functions are used to describe the target states. The results will be compared with available theoretical and experimental data.

Huo, Winifred M.↗

Vibration analysis of the Long Duration Exposure Facility (LDEF) using SPAR

The structural modeling of the Long Duration Exposure Facility (LDEF) utilizing the SPAR system of computer programs for vibration analysis is discussed. The technical areas of interest were: (1) development of the LDEF finite element model; (2) derivation of tray effective panel stiffness matrix using finite element tray models; (3) assessment of attachment conditions and end fitting flexibility by comparing SPAR with test static displacements; (4) SPAR grouping; and (5) derivation of the LDEF frequencies and mode shapes and comparing them with tests. Special detailed finite element modeling was required to obtain good agreement between analytical and test vibration modes. An orthotropic panel in the overall model was developed. Orthotropic stiffness for this panel were obtained from finely detailed statically loaded SPAR models which included stiffness and allowed for partial relative sliding of the tray clamping attachments. Sensitivity to LDEF joint boundary conditions was determined, and static test data proved valuable in assessing modeling of local end fittings.

Edighoffer, H.↗

Implementation of a Matrix Crack Spacing Parameter in a Continuum Damage Mechanics Finite Element Model

Continuum Damage Mechanics (CDM) based progressive damage and failure analysis (PDFA) methods have demonstrated success in a variety of finite element analysis (FEA) implementations. However, the technical maturity of CDM codes has not yet been proven for the full design space of composite materials in aerospace applications. CDM-based approaches represent the presence of damage by changing the local material stiffness definitions and without updating the original mesh or element integration schemes. Without discretely representing cracks and their paths through the mesh, damage in models with CDM-based materials is often distributed in a region of partially damaged elements ahead of stress concentrations. Having a series of discrete matrix cracks represented by a softened region may affect predictions of damage propagation and, thus, structural failure. This issue can be mitigated by restricting matrix damage development to discrete, fiber-aligned rows of elements; hence CDM-based matrix cracks can be implemented to be more representative of discrete matrix cracks. This paper evaluates the effect of restricting CDM matrix crack development to discrete, fiber-aligned rows where the spacing of these rows is controlled by a user-defined crack spacing parameter. Initially, the effect of incrementally increasing matrix crack spacing in a unidirectional center notch coupon is evaluated. Then, the lessons learned from the center notch specimen are applied to open-hole compression finite element models. Results are compared to test data, and the limitations, successes, and potential of the matrix crack spacing approach are discussed.

Hyder, Imran↗

A constitutive material model for nonlinear finite element structural analysis using an iterative matrix approach

A unified numerical method for the integration of stiff time dependent constitutive equations is presented. The solution process is directly applied to a constitutive model proposed by Bodner. The theory confronts time dependent inelastic behavior coupled with both isotropic hardening and directional hardening behaviors. Predicted stress-strain responses from this model are compared to experimental data from cyclic tests on uniaxial specimens. An algorithm is developed for the efficient integration of the Bodner flow equation. A comparison is made with the Euler integration method. An analysis of computational time is presented for the three algorithms.

Koenig, Herbert A.↗

Matrix methods and automation in structural engineering

The partial differential equations of motion of shell structures subject to arbitrary loads may be extremely difficult or even impossible to solve. Since the digital computer is now an available tool for the structural engineer, current research has been directed toward methods which involve matrix formulation of large systems of algebraic equations and matrix solutions for discrete elements rather than the solution of the partial differential equations. These methods require rapid and accurate computer solutions. Recognizing the accuracy problems inherent in working with large matrices, a comprehensive survey of available computer programs was performed for large matrix inversion and for eigenvalue and eigenvector solutions under a Research Grant from the National Aeronautics and Space Administration and the techniques are developed in this paper for using these programs most efficiently for structural applications. The contents of the paper include: automation of matrix compilation, methods of very large matrix inversion and solution of simultaneous equation, and techniques for finding eigenvalues and vectors. In addition, a finite element stiffness matrix approach developed at the Denver Research Institute for both plates and shells subject to arbitrary dynamic loads is described as it was instrumented with complete automation on the digital computer.

Anita S West↗

Direct finite element solution on an optical laboratory matrix-vector processor

The first optical laboratory system results employing a direct LU decomposition solution of a system of linear algebraic equations are presented for a finite element problem solution. This also represents the first laboratory demonstration of the use of sign-magnitude negative number representation as well as new bit partitioning techniques to increase the accuracy of an optical encoded processor beyond the number of bit channels available.

Casasent, David↗

Micro-area Auger analysis of a SiC/Ti fibre composite

Micro-area Auger electron spectroscopy with a spatial resolution of less than 50 nm has been used to study the concentration of elements across the reaction zone of a W-reinforced SiC fiber in a titanium matrix. Although the elemental concentrations obtained by this technique are affected by the reaction zone morphology to a greater extent than in the case of X-ray microprobe analysis, the proposed technique has the advantage of a much higher spatial resolution and avoids the problems of bulk averaging that characterize the X-ray technique.

Zironi, E. P.↗

Electronic neural network for dynamic resource allocation

A VLSI implementable neural network architecture for dynamic assignment is presented. The resource allocation problems involve assigning members of one set (e.g. resources) to those of another (e.g. consumers) such that the global 'cost' of the associations is minimized. The network consists of a matrix of sigmoidal processing elements (neurons), where the rows of the matrix represent resources and columns represent consumers. Unlike previous neural implementations, however, association costs are applied directly to the neurons, reducing connectivity of the network to VLSI-compatible 0 (number of neurons). Each row (and column) has an additional neuron associated with it to independently oversee activations of all the neurons in each row (and each column), providing a programmable 'k-winner-take-all' function. This function simultaneously enforces blocking (excitatory/inhibitory) constraints during convergence to control the number of active elements in each row and column within desired boundary conditions. Simulations show that the network, when implemented in fully parallel VLSI hardware, offers optimal (or near-optimal) solutions within only a fraction of a millisecond, for problems up to 128 resources and 128 consumers, orders of magnitude faster than conventional computing or heuristic search methods.

Thakoor, A. P.↗

Statistical significance test for transition matrices of atmospheric Markov chains

Low-frequency variability of large-scale atmospheric dynamics can be represented schematically by a Markov chain of multiple flow regimes. This Markov chain contains useful information for the long-range forecaster, provided that the statistical significance of the associated transition matrix can be reliably tested. Monte Carlo simulation yields a very reliable significance test for the elements of this matrix. The results of this test agree with previously used empirical formulae when each cluster of maps identified as a distinct flow regime is sufficiently large and when they all contain a comparable number of maps. Monte Carlo simulation provides a more reliable way to test the statistical significance of transitions to and from small clusters. It can determine the most likely transitions, as well as the most unlikely ones, with a prescribed level of statistical significance.

Vautard, Robert↗

Accuracy and speed in computing the Chebyshev collocation derivative

We studied several algorithms for computing the Chebyshev spectral derivative and compare their roundoff error. For a large number of collocation points, the elements of the Chebyshev differentiation matrix, if constructed in the usual way, are not computed accurately. A subtle cause is is found to account for the poor accuracy when computing the derivative by the matrix-vector multiplication method. Methods for accurately computing the elements of the matrix are presented, and we find that if the entities of the matrix are computed accurately, the roundoff error of the matrix-vector multiplication is as small as that of the transform-recursion algorithm. Results of CPU time usage are shown for several different algorithms for computing the derivative by the Chebyshev collocation method for a wide variety of two-dimensional grid sizes on both an IBM and a Cray 2 computer. We found that which algorithm is fastest on a particular machine depends not only on the grid size, but also on small details of the computer hardware as well. For most practical grid sizes used in computation, the even-odd decomposition algorithm is found to be faster than the transform-recursion method.

Don, Wai-Sun↗

Parallel conjugate gradient algorithms for manipulator dynamic simulation

Parallel conjugate gradient algorithms for the computation of multibody dynamics are developed for the specialized case of a robot manipulator. For an n-dimensional positive-definite linear system, the Classical Conjugate Gradient (CCG) algorithms are guaranteed to converge in n iterations, each with a computation cost of O(n); this leads to a total computational cost of O(n sq) on a serial processor. A conjugate gradient algorithms is presented that provide greater efficiency using a preconditioner, which reduces the number of iterations required, and by exploiting parallelism, which reduces the cost of each iteration. Two Preconditioned Conjugate Gradient (PCG) algorithms are proposed which respectively use a diagonal and a tridiagonal matrix, composed of the diagonal and tridiagonal elements of the mass matrix, as preconditioners. Parallel algorithms are developed to compute the preconditioners and their inversions in O(log sub 2 n) steps using n processors. A parallel algorithm is also presented which, on the same architecture, achieves the computational time of O(log sub 2 n) for each iteration. Simulation results for a seven degree-of-freedom manipulator are presented. Variants of the proposed algorithms are also developed which can be efficiently implemented on the Robot Mathematics Processor (RMP).

Fijany, Amir↗

Aircraft electric field measurements: Calibration and ambient field retrieval

An aircraft locally distorts the ambient thundercloud electric field. In order to determine the field in the absence of the aircraft, an aircraft calibration is required. In this work a matrix inversion method is introduced for calibrating an aircraft equipped with four or more electric field sensors and a high-voltage corona point that is capable of charging the aircraft. An analytic, closed form solution for the estimate of a (3 x 3) aircraft calibration matrix is derived, and an absolute calibration experiment is used to improve the relative magnitudes of the elements of this matrix. To demonstrate the calibration procedure, we analyze actual calibration date derived from a Lear jet 28/29 that was equipped with five shutter-type field mill sensors (each with sensitivities of better than 1 V/m) located on the top, bottom, port, starboard, and aft positions. As a test of the calibration method, we analyze computer-simulated calibration data (derived from known aircraft and ambient fields) and explicitly determine the errors involved in deriving the variety of calibration matrices. We extend our formalism to arrive at an analytic solution for the ambient field, and again carry all errors explicitly.

Koshak, William J.↗

Strain Rate Dependent Deformation and Strength Modeling of a Polymer Matrix Composite Utilizing a Micromechanics Approach

Potential gas turbine applications will expose polymer matrix composites to very high strain rate loading conditions, requiring an ability to understand and predict the material behavior under extreme conditions. Specifically, analytical methods designed for these applications must have the capability of properly capturing the strain rate sensitivities and nonlinearities that are present in the material response. The Ramaswamy-Stouffer constitutive equations, originally developed to analyze the viscoplastic deformation of metals, have been modified to simulate the nonlinear deformation response of ductile, crystalline polymers. The constitutive model is characterized and correlated for two representative ductile polymers. Fiberite 977-2 and PEEK, and the computed results correlate well with experimental values. The polymer constitutive equations are implemented in a mechanics of materials based composite micromechanics model to predict the nonlinear, rate dependent deformation response of a composite ply. Uniform stress and uniform strain assumptions are applied to compute the effective stresses of a composite unit cell from the applied strains. The micromechanics equations are successfully verified for two polymer matrix composites. IM7/977-2 and AS4/PEEK. The ultimate strength of a composite ply is predicted with the Hashin failure criteria that were implemented in the composite micromechanics model. The failure stresses of the two composite material systems are accurately predicted for a variety of fiber orientations and strain rates. The composite deformation model is implemented in LS-DYNA, a commercially available transient dynamic explicit finite element code. The matrix constitutive equations are converted into an incremental form, and the model is implemented into LS-DYNA through the use of a user defined material subroutine. The deformation response of a bulk polymer and a polymer matrix composite are predicted by finite element analyses. The results compare reasonably well to experimental values, with some discrepancies. The discrepancies are at least partially caused by the method used to integrate the rate equations in the polymer constitutive model.

Goldberg, Robert K.↗

Lunar meteorite Yamato-86032 - Mineralogical, petrological, and geochemical studies

Yamato-86032 is a shock-lithified anorthositic fragmental breccia. It consists mainly of highly feldspathic meta-breccias and meta-meltrocks and possibly contains a small contribution from mare lithologies, but there is no indication of a KREEP component. In many respects Y-86032 is similar to the previously described lunar meteorites Y-82192/3, but there are some notable differences. We have analyzed about 40 major and trace elements in bulk matrix, impact melt, and clast samples from two chips of Y-86032. The abundances of most lithophile and incompatible elements are lower in Y-86032 than in Y-82192 (which contains very low abundances compared to normal lunar highland rocks). The REE abundances are comparable to those of Y-82192. The elements Sc, Cr, Mn, Fe and Co have significantly lower abundances than in Y-82192, and the siderophile element pattern is also different. Since cosmic ray exposure data indicate pairing of Y-86032 with Y-82192/3, the source region for these meteorites on the moon must have been fairly heterogeneous.

Koeberl, Christian↗

Elastic-plastic finite element analyses of an unidirectional, 9 vol percent tungsten fiber reinforced copper matrix composite

Micromechanical modeling via elastic-plastic finite element analyses were performed to investigate the effects that the residual stresses and the degree of matrix work hardening (i.e., cold-worked, annealed) have upon the behavior of a 9 vol percent, unidirectional W/Cu composite, undergoing tensile loading. The inclusion of the residual stress-containing state as well as the simulated matrix material conditions proved to be significant since the Cu matrix material exhibited plastic deformation, which affected the subsequent tensile response of the composite system. The stresses generated during cooldown to room temperature from the manufacturing temperature were more of a factor on the annealed-matrix composite, since they induced the softened matrix to plastically flow. This event limited the total load-carrying capacity of this matrix-dominated, ductile-ductile type material system. Plastic deformation of the hardened-matrix composite during the thermal cooldown stage was not considerable, therefore, the composite was able to sustain a higher stress before showing any appreciable matrix plasticity. The predicted room temperature, stress-strain response, and deformation stages under both material conditions represented upper and lower bounds characteristic of the composite's tensile behavior. The initial deformation stage for the hardened material condition showed negligible matrix plastic deformation while for the annealed state, its initial deformation stage showed extensive matrix plasticity. Both material conditions exhibited a final deformation stage where the fiber and matrix were straining plastically. The predicted stress-strain results were compared to the experimental, room temperature, tensile stress-strain curve generated from this particular composite system. The analyses indicated that the actual thermal-mechanical state of the composite's Cu matrix, represented by the experimental data, followed the annealed material condition.

Sanfeliz, Jose G.↗