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 1,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Analysis of matrix cracking and local delamination in (O/theta/-theta)sub s graphite epoxy laminates under tension load

Several 3D finite element analyses of (O/theta/-theta)sub s graphite epoxy laminates, where theta = 15, 20, 25, 30, and 45 deg, subjected to axial tension load were performed. The interlaminar stresses in the theta/-theta interface were calculated with and without a matrix crack in the central -theta plies. The interlaminar normal stress changes from a small compressive stress when no matrix crack is present to a high tensile stress at the intersection of the matrix crack and free edge. The analysis of local delamination from the -theta matrix crack indicates a high strain energy release rate and a localized mode I component near the free edge, within one ply distance from the matrix crack. In order to examine the stress state causing the matrix cracking the maximum principal normal stress in a plane perpendicular to the fiber direction in the -theta ply was calculated in an uncracked laminate. The corresponding shear stress parallel to the fiber was also calculated. The principal normal stress at the laminate edge increases through the ply thickness and reached a very high tensile value at the theta/-theta interface indicating that the crack in the -theta ply may initiate at the theta/-theta interface. Crack profiles on the laminate edge in the -theta ply were constructed from the principal stress directions. The cracks were found to be more curved for layups with smaller theta angles, which is consistent with experimental observations in the literature.

Salpekar, S. A.↗
DOE OSTI2020

Microstructures and mechanical properties of a modified 9Cr ferritic-martensitic steel in the as-built condition after additive manufacturing

A newly developed nano-structured high-Mn 9Cr ferritic-martensitic (FM) steel designed for additive manufacturing (Additive-manufactured Nanostructured Alloy, ANA) was fabricated via the direct energy deposition method. The as-printed ANA2 showed tensile properties (strength and elongation) and Charpy impact toughness (upper-shelf energy and ductile-brittle transition temperature) comparable to oxide-dispersion-strengthened (ODS) Eurofer and PM2000. Compared with wrought P91, the as-built ANA2 showed significantly higher yield strength but inferior Charpy impact toughness and ductility. Microstructural characterization revealed uniformly distributed sub-micron pores, an extremely high dislocation density, and a high density of ultrafine nano-structured precipitates in the matrix. The microstructures indicated the dislocations as the dominant strengthening element and the ultrafine precipitates as the primary contribution to estimated irradiation sink density. This work demonstrates the ability to exploit the unique consolidation properties of additive manufacturing to fabricate steels which marry the benefits of both ODS alloys and conventionally processed FM steels.

36 MATERIALS SCIENCE↗
DOE OSTI2021

FEA-aided investigation of the effective thermal conductivity in a medium with embedded spheres

For multiple applications in nuclear energy, the ability to accurately represent material behavior with a simplified model is important to facilitate practical engineering-scale simulations. Here, we focus on the homogenized thermal response of a medium containing spherical inclusions, similar to a fuel form (compact or pebble) containing (TRISO) particles. A review on effective thermal conductivity (ETC) modeling is performed considering a random distribution of mono-sized spherical inclusions in a continuous matrix, with a primary focus on the analytical models. Finite element simulations are performed to evaluate each analytical model in varying material conditions. The model predictions are compared with the expected results obtained from the finite element predictions in addition to the Wiener and Hashin-Shtrikman bounds. Lastly, we included a practical discussion on the homogenization applied to a TRISO fuel pebble.

42 ENGINEERING↗
DOE OSTI2020

Prediction of room-temperature multiferroicity in strained MoCr 2 S 6 monolayer

The contrasting d-orbital occupation required for ferroelectricity vs ferromagnetism makes it difficult for their coexistence in two-dimensional materials, especially at high temperature. To resolve this intrinsic contradiction, we propose a layered MoCr 2 S 6 multiferroics by alloying magnetic Cr element into the ferroelectric 1T phase of the MoS 2 matrix. First-principles calculations disclose that a spontaneous symmetry breaking, depending on the Mo atom displacement, leads to a robust ferroelectricity, which coexists with a ferromagnetic order originated from two neighboring Cr atoms. The effect can be further enhanced by tensile strain to bring about a room-temperature multiferroicity. Our results shed new light on the fundamental understanding of multiferroics and display promising applications in spintronics and multistate data storage.

2D materials↗
DOE OSTI2024

Thermal quasiparticle theory

The widely used thermal Hartree–Fock (HF) theory is generalized to include the effect of electron correlation while maintaining its quasi-independent-particle framework. An electron-correlated internal energy (or grand potential) is postulated in consultation with the second-order finite-temperature many-body perturbation theory (MBPT), which then dictates the corresponding thermal orbital (quasiparticle) energies in such a way that all fundamental thermodynamic relations are obeyed. The associated density matrix is of a one-electron type, whose diagonal elements take the form of the Fermi–Dirac distribution functions, when the grand potential is minimized. The formulas for the entropy and chemical potential are unchanged from those of Fermi–Dirac or thermal HF theory. The theory thus stipulates a finite-temperature extension of the second-order Dyson self-energy of one-particle many-body Green’s function theory and can be viewed as a second-order, diagonal, frequency-independent, thermal inverse Dyson equation. At low temperatures, the theory approaches finite-temperature MBPT of the same order, but it may outperform the latter at intermediate temperatures by including additional electron-correlation effects through orbital energies. Here, a physical meaning of these thermal orbital energies is proposed (encompassing that of thermal HF orbital energies, which has been elusive) as a finite-temperature version of Janak’s theorem.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Fermi-liquid theory and divergences of the two-particle irreducible vertex in the periodic Anderson lattice

Here we analyze the divergences of the irreducible vertex function in dynamical mean field theory, which may indicate either a nonphysical breakdown of the perturbation theory or a response to some physical phenomenon. To investigate this question, we construct a quasiparticle vertex from the diverging irreducible vertex functions. This vertex describes the scattering between quasiparticles and quasiholes in a Fermi liquid. We show that the quasparticle vertex does not diverge in the charge channel, wherein the irreducible vertex does diverge; and we show that the quasiparticle vertex does diverge in the spin channel, wherein the irreducible vertex does not diverge. This divergence occurs at the Mott transition, wherein the Fermi-liquid theory breaks down. Both the half filled Hubbard and Anderson lattices are investigated. In general, our results support that the divergences of the irreducible vertex function do not indicate a nonphysical failure of the perturbation theory. Instead, the divergences are the mathematical consequence of inverting a matrix (the local charge susceptibility) which accumulates increasingly negative diagonal elements as the Hubbard interaction suppresses charge fluctuations. Indeed, we find that the first divergences of the irreducible vertex in both Hubbard and Anderson lattices occurs near the maximum magnitude of the (negative) vertex-connected part of the charge susceptibility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1978

Steady, Oscillatory and Unsteady, Subsonic and Supersonic Aerodynamics (SOUSSA) for complex aircraft configurations

The Green's function method and the computer program SOUSSA (Steady Oscillatory and Unsteady Subsonic and Supersonic Aerodynamics) are reviewed. The Green's function method is applied to the fully unsteady potential equation yielding an integro-differential-delay equation. This equation is approximated by a set of differential-delay equations in time using the finite element method. The Laplace transform is used to yield a matrix relating the velocity potential to the normal wash. The matrix of the generalized aerodynamic forces is obtained by premultiplying and postmultiplying the matrices relating generalized forces to the potential and the normal wash by the generalized coordinates. The program SOUSSA is compared with existing numerical results. Results indicate that the program is not only general, flexible, and easy to use, but also accurate and fast.

Morino, L.↗
NASA NTRS1984

Investigation of Toughened Neat Resins and Their Relations to Advanced Composite Mechanical Properties

Four epoxies (Hercules 3502, Hercules 2220-1 and 2220-3, and Ciba-Geigy Fibredux 914 epoxy) were evaluated to judge their merits within the ACEE program and also to provide material property input for a micromechanics finite element program. A casting procedure was developed for the four matrix systems mentioned and mechanical properties were studied at three test temperatures and two moisture contents. A scanning electron microscope study was performed to catalogue fracture surfaces for the unreinforced (neat) epoxies. In addition, tension testing, shear testing, a limited amount of fracture toughness testing, and thermal expansion coefficient and moisture expansion coefficient testing were carried out. Mechanical material properties were generated for input into a micromechanics program and to characterize these four matrix systems.

R. S. Zimmerman↗
NASA NTRS1984

Dynamic thermal episodes in the protosolar nebula: Development of models from observations on CAI's

Evaluation of earlier observations indicated that layered rims on coarse-grained Allende CAI's were possibly the result of partial melting by ablation/drag-heating and reaction of CAI exteriors with a gas or gases of non-solar composition. Bunch and Chang reported the common occurrence of thin, fine-grained, matrix-like bands that at least partially surround rims of CAI's. Although material in these bands in general appears to be similar to matrix, SEM observations show them to be dissimilar in volatile element content, mineral composition, and grain morphology. Moreover, they appear to be related in time of formation with rim development and Na-metasomatism of CAI's. Observations indicate a short-lived but intense heating episode followed by radid cooling as the mechanism responsible for these CAI features.

Bunch, T. E.↗
NASA NTRS1984

A finite element method for nonlinear forced vibrations of beams

Techniques for defining a finite element model (FEM) for analysis of nonlinear vibrations in beam structures subjected to harmonic excitation are presented. The resulting model covers longitudinal deformation and inertial effects. The nonlinear oscillations of a beam element under forced excitation are modeled by a harmonic force matrix based on first order approximations of the Jacobian elliptic forcing function. Harmonic force and nonlinear stiffness matrices are derived and the nonlinear forced responses of beams are calculated under various boundary conditions. The results of FEM computations for simply-supported and clamped beams show that midplane stretching caused by large deflections increases the nonlinearity. Axially-restrained beams experience only hardening nonlinearity, while axially-free beams have reduced nonlinearity in deformation and inertia and an increase in linearity due to large deflection.

Mei, C.↗
NASA NTRS1986

Free vibration characteristics of multiple load path blades by the transfer matrix method

The determination of free vibrational characteristics is basic to any dynamic design, and these characteristics can form the basis for aeroelastic stability analyses. Conventional helicopter blades are typically idealized as single-load-path blades, and the transfer matrix method is well suited to analyze such blades. Several current helicopter dynamic programs employ transfer matrices to analyze the rotor blades. In this paper, however, the transfer matrix method is extended to treat multiple-load-path blades, without resorting to an equivalent single-load-path approximation. With such an extension, these current rotor dynamic programs which employ the transfer matrix method can be modified with relative ease to account for the multiple load paths. Unlike the conventional blades, the multiple-load-path blades require the introduction of the axial degree-of-freedom into the solution process to account for the differential axial displacements of the different load paths. The transfer matrix formulation is validated through comparison with the finite-element solutions.

Murthy, V. R.↗
NASA NTRS1989

Remote sensing of earth terrain

A mathematically rigorous and fully polarimetric radar clutter model used to evaluate the radar backscatter from various types of terrain clutter such as forested areas, vegetation canopies, snow covered terrains, or ice fields is presented. With this model, the radar backscattering coefficients for the multichannel polarimetric radar returns can be calculated, in addition to the complex cross correlation coefficients between elements of the polarimetric measurement vector. The complete polarization covariance matrix can be computed and the scattering properties of the clutter environment characterized over a broad range of incident angle and frequencies.

Kong, J. A.↗
NASA NTRS1992

Fiber shape effects on metal matrix composite behavior

The effects of different fiber shapes on the behavior of a SiC/Ti-15 metal matrix composite is computationally simulated. A three-dimensional finite element model consisting of a group of nine unidirectional fibers is used in the analysis. The model is employed to represent five different fiber shapes: a circle, an ellipse, a kidney, and two different cross shapes. The distribution of microstresses and the composite material properties, such as moduli, coefficients of thermal expansion, and Poisson's ratios, are obtained from the finite element analysis for the various fiber shapes. Comparisons of these results are used to determine the sensitivity of the composite behavior to the different fiber shapes and assess their potential benefits. No clear benefits result from different fiber shapes though there are some increases/decreases in isolated properties.

Brown, H. C.↗
NASA NTRS1993

The effect of nonsymmetric pressure stiffness on the dynamic characteristics of Solid Rocket Motors

This paper discusses the effect of pressure on the dynamics of pre-stiffened structures such as the Advanced Solid Rocket Motor (ASRM). Previous work in which the stiffness terms resulting from constant pressure were derived has been extended to enable modeling of nonconstant pressure applied over nonenclosed volumes. These conditions will result in nonsymmetric terms in the global stiffness matrix which will not cancel out. Three new pressure stiffness elements incorporating these nonsymmetric terms have been implemented as dummy elements in COSMIC NASTRAN and have been tested on various simple examples as well as an existing ASRM NASTRAN finite element model. The results indicate that for all load cases of practical interest to the ASRM program, the nonsymmetric terms have very little effect on the dynamic characteristics. In addition, the pressure stiffness elements developed in the previous work which assumed constant pressure gave virtually the same results as the new elements even for problems in which the pressures are not constant. The original elements appear to work well as long as the pressure gradient across any individual element is no larger than about 0.75 psi/inch. The new elements are therefore most useful for determining the conditions under which the original pressure stiffness elements can be used.

Christensen, Eric R.↗
NASA NTRS1990

Micromechanical thermal analysis of interphase region in a titanium aluminide MMC

The high reactivity between the fiber and matrix in silicon carbide/titanium aluminide MMCs leads to the formation of brittle reaction products at the fiber/matrix (F/M) interface. Also, the high thermal expansion coefficient mismatch between the fiber and matrix leads to high tensile residual stresses at the F/M interface, and this can lead to premature cracking during cooldown. One solution to these problems is the use of a metallic fiber coating like Ta which acts as an F/M diffusion barrier and reacts with the matrix to form a beta stabilized compliant layer. A finite element micromechanics analysis was performed to study the effects of Ta and beta interphase layers on the thermal residual stresses during consolidation. A 5-micron-thick beta layer reduced cool-down stresses by 8 percent compared to a 2 percent reduction computed for a Ta layer of the same thickness. Plastic yielding in the Ta was not effective in reducing cool-down stresses. Compliant alpha-2 particles next to the stiffer gamma particles reduced stresses in the gamma particles by less than 2 percent. A simple closed form analysis was developed to calculate thermal residual stresses in a fiber/interphase/matrix system. A 2-micron-thick Ag interphase layer was found to reduce residual stresses by about 11 percent.

Naik, R. A.↗
NASA NTRS2002

Genetic Optimization of a Tensegrity Structure

Marshall Space Flight Center (MSFC) is charged with developing advanced technologies for space telescopes. The next generation of space optics will be very large and lightweight. Tensegrity structures are built of compressive members (bars), and tensile members (strings). For most materials, the tensile strength of a longitudinal member is larger than its buckling strength; therefore a large stiffness to mass ratio can be achieved by increasing the use of tensile members. Tensegrities are the epitome of lightweight structures, since they take advantage of the larger tensile strength of materials. The compressive members of tensegrity structures are disjoint allowing compact storage of the structure. The structure has the potential to eliminate the requirement for assembly by man in space; it can be deployed by adjustments in its cable tension. A tensegrity structure can be more reliably modeled since none of the individual members experience bending moments. (Members that experience deformation in more than one dimension are much harder to model.) Structures that can be more precisely modeled can be more precisely controlled. Furthermore, an astoundingly wide variety of natural systems, including carbon atoms, water molecules, proteins, viruses, cells, tissues and even human and other living creatures are tensegrity structures. Through the process of evolution, nature continually improves the design of living creatures for the environment they live in. Since tensegrities are nature's structure of choice, it is conceivable that they have other benefits we are unaware of. A. Keane and S. Brown designed a satellite boom truss system with an enhanced vibration performance. They started with a standard truss system, then used a genetic algorithm to alter the design, while optimizing the vibration performance. An improvement of over 20,000% in frequency-averaged energy levels was obtained using this approach. In this report an introduction to tensegrity structures is given, along with a description of how to generate the nodal coordinates and connectivity of a multiple stage cylindrical tensegrity structure. A description of how finite elements can be used to develop a stiffness and mass matrix so that the modes of vibration can be determined from the eigenvalue problem is shown. A brief description of a micro genetic algorithm is then presented.

Jaime R. Taylor↗
NASA NTRS2004

Electron-Hydrogen Elastic Scattering

Scattering by single-electron systems is always of interest because the wave function of the target is known exactly. Various approximations have been employed to take into account distortion produced in the target. Among them are the method of polarized orbitals and the close coupling approximation. Recently, e-H and e-He+ S-wave scattering in the elastic region has been studied using the Feshbach projection operator formalism. In this approach, the usual Hartree-Fock and exchange potentials are augmented by an optical potential and the resulting phase shifts have rigorous lower bounds. Now this method is being applied to the e-H P-wave scattering in the elastic region. The number of terms in the Hylleraas-type wave function for the 1,3 P phase shifts is 84 and the resulting phase shifts (preliminary) are given. The results have been given up to five digits because to that accuracy they are rigorous lower bounds. They are in general agreement with the variational (VAR) results of Armstead, and those obtained from the intermediate energy R-matrix method (RM) of Scholz et al., and the finite element method (FEM) of Botero and Shertzer. The later two methods do not provide any bounds on phase shifts.

Bhatia, A. K.↗