Search NASA⌕ Search

SEARCH · Search NASA

Results for “MATRIX STRESS CALCULATION”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Thermo-viscoelastic analysis of composite materials, volume 1

Advanced composite materials, especially graphite/epoxy, are being applied to aircraft structures in order to improve performance and save weight. An important consideration in composite design is the residual strength of a structure containing holes, delaminations, or interlaminar damage when subjected to compressive loads. Recent studies have revealed the importance of viscoelastic effects in polymer-based composites. The viscoelastic effect is particularly significant at elevated temperature/moisture conditions since the matrix material is strongly affected by the environment. The solution of viscoelastic problems in composites was limited to special cases which can be solved by classical lamination theory. A finite element procedure is presented for calculating time-dependent stresses and strains in composite structures with general configurations and complicated boundary conditions. Using this procedure the in-plane and interlaminar stress distributions and histories in notched and unnotched composites were obtained for mechanical and thermal loads. Both two-dimensional and three-dimensional viscoelastic problems are analyzed. The effects of layup orientation and load spectrum on creep response and stress relaxation were also studied.

Lin, K. Y.↗
DOE OSTI2020

Real-space formulation of the stress tensor for O (N) density functional theory: Application to high temperature calculations

We present an accurate and efficient real-space formulation of the Hellmann–Feynman stress tensor for O (N) Kohn–Sham density functional theory (DFT). While applicable at any temperature, the formulation is most efficient at high temperature where the Fermi–Dirac distribution becomes smoother and the density matrix becomes correspondingly more localized. We first rewrite the orbital-dependent stress tensor for real-space DFT in terms of the density matrix, thereby making it amenable to O (N) methods. We then describe its evaluation within the O (N) infinite-cell Clenshaw–Curtis Spectral Quadrature (SQ) method, a technique that is applicable to metallic and insulating systems, is highly parallelizable, becomes increasingly efficient with increasing temperature, and provides results corresponding to the infinite crystal without the need of Brillouin zone integration. We demonstrate systematic convergence of the resulting formulation with respect to SQ parameters to exact diagonalization results and show convergence with respect to mesh size to the established plane wave results. Here, we employ the new formulation to compute the viscosity of hydrogen at 10 6 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2004

Modeling the Nonlinear, Strain Rate Dependent Deformation of Shuttle Leading Edge Materials with Hydrostatic Stress Effects Included

An analysis method based on a deformation (as opposed to damage) approach has been developed to model the strain rate dependent, nonlinear deformation of woven ceramic matrix composites, such as the Reinforced Carbon Carbon (RCC) material used on the leading edges of the Space Shuttle. In the developed model, the differences in the tension and compression deformation behaviors have also been accounted for. State variable viscoplastic equations originally developed for metals have been modified to analyze the ceramic matrix composites. To account for the tension/compression asymmetry in the material, the effective stress and effective inelastic strain definitions have been modified. The equations have also been modified to account for the fact that in an orthotropic composite the in-plane shear response is independent of the stiffness in the normal directions. The developed equations have been implemented into LS-DYNA through the use of user defined subroutines (UMATs). Several sample qualitative calculations have been conducted, which demonstrate the ability of the model to qualitatively capture the features of the deformation response present in woven ceramic matrix composites.

Goldberg, Robert K.↗
NASA NTRS1977

Comparison of linear inversion methods by examination of the duality between iterative and inverse matrix methods

Linear numerical inversion methods applied to atmospheric remote sounding generally can be categorized in two ways: (1) iterative, and (2) inverse matrix methods. However, these two categories are not unrelated; a duality exists between them. In other words, given an iterative scheme, a corresponding inverse matrix method exists, and conversely. This duality concept is developed for the more familiar linear methods. The iterative duals are compared with the classical linear iterative approaches and their differences analyzed. The importance of the initial profile in all methods is stressed. Calculations using simulated data are made to compare accuracies and to examine the dependence of the solution on the initial profile.

Fleming, H. E.↗
NASA NTRS1989

Composite blade structural analyzer (COBSTRAN) user's manual

The installation and use of a computer code, COBSTRAN (COmposite Blade STRuctrual ANalyzer), developed for the design and analysis of composite turbofan and turboprop blades and also for composite wind turbine blades was described. This code combines composite mechanics and laminate theory with an internal data base of fiber and matrix properties. Inputs to the code are constituent fiber and matrix material properties, factors reflecting the fabrication process, composite geometry and blade geometry. COBSTRAN performs the micromechanics, macromechanics and laminate analyses of these fiber composites. COBSTRAN generates a NASTRAN model with equivalent anisotropic homogeneous material properties. Stress output from NASTRAN is used to calculate individual ply stresses, strains, interply stresses, thru-the-thickness stresses and failure margins. Curved panel structures may be modeled providing the curvature of a cross-section is defined by a single value function. COBSTRAN is written in FORTRAN 77.

Aiello, Robert A.↗
NASA NTRS1993

Analysis of delamination in cross-ply laminates initiating from impact induced matrix cracking

Two-dimensional finite element analyses of (02/90(8)/02) glass/epoxy and graphite/epoxy composite laminates were performed to investigate some of the characteristics of damage development due to an impact load. A cross section through the thickness of the laminate with fixed ends, and carrying a transverse load in the center, was analyzed. Inclined matrix cracks, such as those produced by a low-velocity impact, were modeled in the 90 deg ply group. The introduction of the matrix cracks caused large interlaminar tensile and shear stresses in the vicinity of both crack tips in the 0/90 and 90/0 interfaces, indicating that matrix cracking may give rise to delamination. The ratio of Mode I to total strain energy release rate, G(I)/G(total), at the beginning of delamination, calculated at the two (top and bottom) matrix crack tips was 60 and 28 percent, respectively, in the glass/epoxy laminate. The corresponding ratio was 97 and 77 percent in the graphite/epoxy laminate. Thus, a significant Mode I component of strain energy release rate may be present at the delamination initiation due to an impact load. The value of strain energy release rate at either crack tip increased due to an increase in the delamination length at the other crack tip and may give rise to an unstable delamination growth under constant load.

Salpekar, S. A.↗
NASA NTRS1990

Integrated modeling and analysis of a space-truss article

MSC/NASTRAN is being used in the Controls-Structures Interaction (CSI) program at NASA Langley Research Center as a key analytical tool for structural analysis as well as the basis for control law development, closed-loop performance evaluation, and system safety checks. Guest investigators from academia and industry are performing dynamics and control experiments on a flight-like deployable space truss called Mini-Mast to determine the effectiveness of various active-vibration control laws. MSC/NASTRAN was used to calculate natural frequencies and mode shapes below 100 Hz to describe the dynamics of the 20-meter-long lightweight Mini-Mast structure. Gravitational effects contribute significantly to structural stiffness and are accounted for through a two-phase solution in which the differential stiffness matrix is calculated and then used in the eigensolution. Reduced modal models are extracted for control law design and evaluation of closed-loop system performance. Predicted actuator forces from controls simulations are then applied to the extended model to predict member loads and stresses. These pre-test analyses reduce risks associated with the structural integrity of the test article, which is a major concern in closed-loop control experiments due to potential instabilities.

Alan E Stockwell↗
NASA NTRS2005

NASALife-Component Fatigue and Creep Life Prediction Program and Illustrative Examples

NASALife is a life prediction program for propulsion system components made of ceramic matrix composites (CMC) under cyclic thermo-mechanical loading and creep rupture conditions. Although, the primary focus was for CMC components the underlying methodologies are equally applicable to other material systems as well. The program references data for low cycle fatigue (LCF), creep rupture, and static material properties as part of the life prediction process. Multiaxial stresses are accommodated by Von Mises based methods and a Walker model is used to address mean stress effects. Varying loads are reduced by the Rainflow counting method. Lastly, damage due to cyclic loading (Miner s rule) and creep are combined to determine the total damage per mission and the number of missions the component can survive before failure are calculated. Illustration of code usage is provided through example problem of a CMC turbine stator vane made of melt-infiltrated, silicon carbide fiber-reinforced, silicon carbide matrix composite (MI SiC/SiC)

Murthy, Pappu L. N.↗
NASA NTRS2000

Equivalent Skin Analysis of Wing Structures Using Neural Networks

An efficient method of modeling trapezoidal built-up wing structures is developed by coupling. in an indirect way, an Equivalent Plate Analysis (EPA) with Neural Networks (NN). Being assumed to behave like a Mindlin-plate, the wing is solved using the Ritz method with Legendre polynomials employed as the trial functions. This analysis method can be made more efficient by avoiding most of the computational effort spent on calculating contributions to the stiffness and mass matrices from each spar and rib. This is accomplished by replacing the wing inner-structure with an "equivalent" material that combines to the skin and whose properties are simulated by neural networks. The constitutive matrix, which relates the stress vector to the strain vector, and the density of the equivalent material are obtained by enforcing mass and stiffness matrix equities with rec,ard to the EPA in a least-square sense. Neural networks for the material properties are trained in terms of the design variables of the wing structure. Examples show that the present method, which can be called an Equivalent Skin Analysis (ESA) of the wing structure, is more efficient than the EPA and still fairly good results can be obtained. The present ESA is very promising to be used at the early stages of wing structure design.

Liu, Youhua↗
NASA NTRS1986

Ultrasonic determination of the elastic constants of the stiffness matrix for unidirectional fiberglass epoxy composites

The elastic constants of a fiberglass epoxy unidirectional composite are determined by measuring the phase velocities of longitudinal and shear stress waves via the through transmission ultrasonic technique. The waves introduced into the composite specimens were generated by piezoceramic transducers. Geometric lengths and the times required to travel those lengths were used to calculate the phase velocities. The model of the transversely isotropic medium was adopted to relate the velocities and elastic constants.

Marques, E. R. C.↗
DOE OSTI2018

Structural Phase Transformation in Strained Monolayer MoWSe 2 Alloy

Two-dimensional (2D) materials exhibit different mechanical properties from their bulk counterparts owing to their monolayer atomic thickness. In this paper, we have examined the mechanical behavior of 2D molybdenum tungsten diselenide (MoWSe 2 ) precipitation alloy grown using chemical vapor deposition and composed of numerous nanoscopic MoSe 2 and WSe 2 regions. Applying a bending strain blue-shifted the MoSe 2 and WSe 2 A 1g Raman modes with the stress concentrated near the precipitate interfaces predominantly affecting the WSe 2 modes. In situ local Raman measurements suggested that the crack propagated primarily thorough MoSe 2 -rich regions in the monolayer alloy. Molecular dynamics (MD) simulations were performed to study crack propagation in an MoSe 2 monolayer containing nanoscopic WSe 2 regions akin to the experiment. Raman spectra calculated from MD trajectories of crack propagation confirmed the emergence of intermediate peaks in the strained monolayer alloy, mirroring experimental results. The simulations revealed that the stress buildup around the crack tip caused an irreversible structural transformation from the 2H to 1T phase both in the MoSe 2 matrix and WSe 2 patches. This was corroborated by high-angle annular dark-field images. Finally, crack branching and subsequent healing of a crack branch were also observed in WSe 2 , indicating the increased toughness and crack propagation resistance of the alloyed 2D MoWSe 2 over the unalloyed counterparts.

36 MATERIALS SCIENCE↗
NASA NTRS1988

A hybrid-stress finite element for linear anisotropic elasticity

Standard assumed displacement finite elements with anisotropic material properties perform poorly in complex stress fields such as combined bending and shear and combined bending and torsion. A set of three dimensional hybrid-stress brick elements were developed with fully anisotropic material properties. Both eight-node and twenty-node bricks were developed based on the symmetry group theory of Punch and Atluri. An eight-node brick was also developed using complete polynomials and stress basis functions and reducing the order of the resulting stress parameter matrix by applying equilibrium constraints and stress compatibility constraints. Here the stress compatibility constraints must be formulated assuming anisotropic material properties. The performance of these elements was examined in numerical examples covering a broad range of stress distributions. The stress predictions show significant improvement over the assumed displacement elements but the calculation time is increased.

Fly, Gerald W.↗
NASA NTRS2018

Influence of Fiber Content, Thermal Cycles and Creep on Electrical Resistivity Measurements of SiCf/SiC Minicomposites

Fiber-reinforced Ceramic Matrix Composites (CMCs) are currently being used in high-temperature applications such as jet engine components. In such extreme environment applications, materials are exposed to many mechanical and thermal loading cycles during the life of the component. This exposure can degrade the component, and thus, components need to be monitored and inspected during service. The electrical resistivity nondestructive evaluation technique has been shown to be very sensitive to damage in CMCs, and although its performance is highly dependent on temperature, this temperature-dependence is still not well understood. Therefore, this current study aims to establish the influence of thermal cycles, creep, CVI-SiC matrix cracks and fiber volume fraction on the electrical resistance of Hi Nicalon and Hi-Nicalon Type S SiC fiber-reinforced minicomposites with BN interphases. Minicomposites with 97% CVI-SiC matrix volume fraction were used to characterize the electrical resistivity of the CVI-SiC matrix. Next, models were constructed to calculate fiber resistivity and onset and evolution of damage during tensile testing at room-temperature. The thermal cycles were achieved by heating as-received and precracked specimens to 1200°C and then cooling them to room temperature, with or without applied stress on fibers.

creep↗
NASA NTRS1982

Mechanics aspects of NDE by sound and ultrasound

Nondestructive evaluation (NDE) is considered as a means to detect the energy release mechanism of defects and the interaction of microstructures within materials with sound waves and/or ultrasonic waves. Ultrasonic inspection involves the frequency range 20 kHz-1 GHz with amplitudes depending on the sensitivity of the test instrumentation. Pulse echo systems are most frequently used in NDE. Information is extracted from the signals through measurements of the signal velocity, attenuation, the acoustic emission when stress is applied, and calculation of the acoustoelastic coefficients. Fracture properties, tensile and shear strengths, the interlaminar shear strength, the cohesive strength, yield and impact strengths, the hardness, and the residual stress can be assayed by ultrasonic methods. Finally, attention is given to analytical treatment of the derived data, with mention given to transition matrix, integral equation, and eigenstrain approaches.

Fu, L. S.↗
DOE OSTI2020

GPU acceleration of all-electron electronic structure theory using localized numeric atom-centered basis functions

We present an implementation of all-electron density-functional theory for massively parallel GPU-based platforms, using localized atom-centered basis functions and real-space integration grids. Special attention is paid to domain decomposition of the problem on non-uniform grids, which enables compute- and memory-parallel execution across thousands of nodes for real-space operations, e.g. the update of the electron density, the integration of the real-space Hamiltonian matrix, and calculation of Pulay forces. To assess the performance of our GPU implementation, we performed benchmarks on three different architectures using a 103-material test set. We find that operations which rely on dense serial linear algebra show dramatic speedups from GPU acceleration: in particular, SCF iterations including force and stress calculations exhibit speedups ranging from 4.5 to 6.6. For the architectures and problem types investigated here, this translates to an expected overall speedup between 3–4 for the entire calculation (including non-GPU accelerated parts), for problems featuring several tens to hundreds of atoms. Additional calculations for a 375-atom Bi2Se3 bilayer show that the present GPU strategy scales for large-scale distributed-parallel simulations.

42 ENGINEERING↗
DOE OSTI2020

Design of a V–Ti–Ni alloy with superelastic nano-precipitates

Stress-induced martensitic transformations enable metastable alloys to exhibit enhanced strain hardening capacity, leading to improved formability and toughness. As is well-known from transformation-induced plasticity (TRIP) steels, however, the resulting martensite can limit ductility and fatigue life due to its intrinsic brittleness. In this work, we explore an alloy design strategy that utilizes stress-induced martensitic transformations but does not retain the martensite phase. This strategy is based on the introduction of superelastic nano-precipitates, which exhibit reverse transformation after initial stress-induced forward transformation. To this end, utilizing ab-initio simulations and thermodynamic calculations we designed and produced a V 45 Ti 30 Ni 25 (at%) alloy. In this alloy, TiNi is present as nano-precipitates uniformly distributed within a ductile V-rich base-centered cubic (bcc) β matrix, as well as being present as a larger matrix phase. We characterized the microstructure of the produced alloy using various scanning electron microscopy (SEM) and transmission electron microscopy (TEM) methods. Furthermore, the bulk mechanical properties of the alloy are demonstrated through tensile tests, and the reversible transformation in each of the TiNi morphologies were confirmed by in-situ TEM micro-pillar compression experiments, in-situ high-energy diffraction synchrotron cyclic tensile tests, indentation experiments, and differential scanning calorimetry experiments. The observed transformation pathways and variables impacting phase stability are critically discussed.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1972

On fracture phenomena in advanced fiber composite materials.

The extension of linear elastic fracture mechanics (LEFM) from metallic alloys to advanced fiber composite laminates is considered. LEFM is shown to be valid for both isotropic and anisotropic homogeneous continua; the applicability of LEFM to advanced fiber composites is thus dependent on the validity of a homogeneous model of such materials. An experimental program to determine the validity of such a model for graphite/epoxy laminates is reviewed. Such laminates are found to have an apparent fracture toughness, from which it is inferred that a homogeneous material model is valid for the particular specimen geometry and composite laminates considered. Strain energy release rates are calculated from the experimentally determined fracture toughness of the various laminates. These strain energy release rates are found to lie in one of two groups, depending upon whether crack extension required fiber failure or matrix failure. The latter case is further investigated. It is concluded that matrix failure is governed by the tensile stress normal to the crack path.

Konish, H. J., Jr.↗
NASA NTRS1988

On the viscoplastic response of a composite wafer

In the present treatment of a viscoplastic composite wafer formed from a viscoplastic matrix that is reinforced by elastic or viscoplastic fibers attached to its plane surfaces, one constitutive equation is established by considering the viscoplastic behavior of the matrix as determined by an integral-type constitutive law whose intrinsic time-measure is pertinent to endochronic viscoplasticity. Attention is given to asymptotic cases where fiber stiffnesses and the hydrostatic modulus of the wafer are much larger than the latter's shear modulus. An explicit calculation is used when the stress field is uniaxial.

Valanis, K. C.↗