Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shear deformation”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

A mixed shear flexible finite element for the analysis of laminated plates

A mixed shear flexible finite element based on the Hencky-Mindlin type shear deformation theory of laminated plates is presented and their behavior in bending is investigated. The element consists of three displacements, two rotations, and three moments as the generalized degrees of freedom per node. The numerical convergence and accuracy characteristics of the element are investigated by comparing the finite element solutions with the exact solutions. The present study shows that reduced-order integration of the stiffness coefficients due to shear is necessary to obtain accurate results for thin plates.

Putcha, N. S.↗
NASA NTRS1987

Nonlinear beam kinematics by decomposition of the rotation tensor

A simple matrix expression is obtained for the strain components of a beam in which the displacements and rotations are large. The only restrictions are on the magnitudes of the strain and of the local rotation, a newly-identified kinematical quantity. The local rotation is defined as the change of orientation of material elements relative to the change of orientation of the beam reference triad. The vectors and tensors in the theory are resolved along orthogonal triads of base vectors centered along the undeformed and deformed beam reference axes, so Cartesian tensor notation is used. Although a curvilinear coordinate system is natural to the beam problem, the complications usually associated with its use are circumvented. Local rotations appear explicitly in the resulting strain expressions, facilitating the treatment of beams with both open and closed cross sections in applications of the theory. The theory is used to obtain the kinematical relations for coupled bending, torsion, extension, shear deformation, and warping of an initially curved and twisted beam.

Danielson, D. A.↗
NASA NTRS1987

Nonlinear Beam Kinematics by Decomposition of the Rotation Tensor

A simple matrix expression is obtained for the strain components of a beam in which the displacements and rotations are large. The only restrictions are on the magnitudes of the strain and of the local rotation, a newly-identified kinematical quantity. The local rotation is defined as the change of orientation of material elements relative to the change of orientation of the beam reference triad. The vectors and tensors in the theory are resolved along orthogonal triads of base vectors centered along the undeformed and deformed beam reference axes, so Cartesian tensor notation is used. Although a curvilinear coordinate system is natural to the beam problem, the complications usually associated with its use are circumvented. Local rotations appear explicitly in the resulting strain expressions, facilitating the treatment of beams with both open and closed cross sections in applications of the theory. The theory is used to obtain the kinematical relations for coupled bending, torsion extension, shear deformation, and warping of an initially curved and twisted beam.

Danielson, D. A.↗
NASA NTRSDate not supplied

Nonlinear Homogenization of Finitely Deformed Viscoelastic-Viscoplastic Composites Using Mechanics of Structure Genome

The objective of this paper is to develop a micromechanics approach to homogenizing finitely deformed viscoelastic-viscoplastic composites using the mechanics of structure genome. The incremental constitutive relation for glassy polymers, formulated in the spatial configuration, is implemented in the present approach.This involves (1) pulling-back the constitutive model to the material configuration and (2)choosing the deformation gradient tensor and the first Piola–Kirchhoff stress tensor as the strain and the stress measures during homogenization, respectively. An Euler–Newton predictor–corrector method is developed for homogenization. Each step involves formulating a variational statement using the mechanics of structure genome, discretizing the statement in a finite-dimensional space, and solving the problem using an Euler/multilevel Newton method. The present approach is demonstrated by homogenizing fiber- and particle-reinforced composites undergoing uniaxial, biaxial, or shear deformation, at different stain rates.

Multi-scale modeling, High Strain Composites, Visc↗
NASA NTRS1999

Accurate Determination of Tranverse Normal Stresses in Sandwich Panels Subjected to Thermomechanical Loadings

A two-stage computational procedure is presented for the accurate determination of transverse normal stresses in sandwich panels subjected to thermomechanical loadings. The procedure involves the use of a first-order shear deformation model in the first stage, and an iterational process for successive improvement of the accuracy of the displacement and stress fields in the second stage. The effectiveness of the procedure is demonstrated by means of numerical studies (if thin and moderately thick flat rectangular sandwich panels. Two sets of boundary conditions are considered; namely, one with all edges simply supported. and the other with two opposite edges simply supported and the remaining two clamped. The displacement components and the transverse shear and normal stresses obtained by the proposed computational procedure are found to be in close agreement with the Solutions of the three-dimensional (3-D) continuum models.

Noor, Ahmed K.↗
NASA NTRS2004

Modeling the Nonlinear, Strain Rate Dependent Deformation of Woven Ceramic Matrix Composites 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 with a plain weave fiber architecture. In the developed model, the differences in the tension and compression response have also been considered. State variable based 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 stiffness is independent of the stiffness in the normal directions. The developed equations have been implemented into a commercially available transient dynamic finite element code, LS-DYNA, through the use of user defined subroutines (UMATs). The tensile, compressive, and shear deformation of a representative plain weave woven ceramic matrix composite are computed and compared to experimental results. The computed values correlate well to the experimental data, demonstrating the ability of the model to accurately compute the deformation response of woven ceramic matrix composites.

Goldberg, Robert K.↗
DOE OSTI2024

Transformational faulting in Mn 2 GeO 4 from olivine to wadsleyite structure: Implications for physical mechanism of deep-focus earthquakes

High-pressure and temperature deformation experiments interfaced with acoustic emission (AE) monitoring have been conducted to study transformational faulting in Mn 2 GeO 4 olivine, which transforms to the β phase, isostructural to wadsleyite. Metastable Mn 2 GeO 4 olivine exhibits a marked embrittlement behavior at temperatures between 800 and 1100 K, emitting numerous AEs. At each temperature, brittle deformation is characterized by a two-stage process: (1) a “preparation” stage with numerous diffusedly located low-magnitude AEs and large b values (>2), and (2) a failure stage where larger-magnitude AEs form a planar distribution with b values about 1. Microstructure analysis reveals extensive kink band development in olivine grains in the recovered samples. Kink band boundaries (KBBs), with a typical thickness of ∼100 nm, are filled with a nanometric β-Mn 2 GeO 4 “gouge”. A dense array of secondary shear localizations is often present within the kink bands, suggesting significant shear deformation therein. The combined observations suggest that faulting in metastable Mn 2 GeO 4 olivine is a self-similar process, from grain-scale to the sample-scale. Both observed embrittlement behavior and the microstructure of metastable Mn 2 GeO 4 olivine are essentially identical to those in Mg 2 GeO 4 olivine we have reported previously, indicating that the physical mechanism of faulting in metastable olivine is insensitive to the specific crystallographic structure of the high-pressure phase. The low b values (about 1) observed in the faulting process in our experiments are similar to those of deep focus earthquakes in cold subduction zones. Our observed mechanism explains deep focus seismicity in cold metastable mantle wedges, provided that the self-similarity assumption holds to geological scales.

58 GEOSCIENCES↗
DOE OSTI2022

Manufacture aluminum alloy tube from powder with a single-step extrusion via ShAPE

Mechanical performance of aluminum in terms of strength, wear and corrosion resistance, especially high-temperature strength has been shown to improve with the addition of transition metal (TM) elements of Fe, Cr, and Ti. However, these feedstock materials occur as powders. As such, the traditional fabrication process is complex and expensive because it requires multiple procedures and consumes considerable energy. This study developed Shear Assisted Processing and Extrusion (ShAPE) as a single-step process that manufactures tubes directly from Al-TM powders obtained via gas atomization. Meter-long Al-TM alloy tubes are extruded from powders with different processing conditions. The ShAPE tubes have very low porosity and the average density is 2.94 kg/cm 3 , which is equal to or higher than parts fabricated by hot extrusion and sintering. The powder-to-tube fabrication process was revealed and discussed by examining the microstructural evolution. The Vickers hardness of ShAPE tubes ranges from 110 to 140 HV through wall thickness and at different extrusion speeds. The variation in hardness was correlated with the extent of refinement of intermetallics and attributed to the shear deformation per unit extrusion length. Energy cost analysis shows that ShAPE save ~60% energy compared to the conventional sintering and extrusion processes. Results indicate ShAPE is a low-cost, high-efficiency manufacturing process for producing tubes from metallic powders.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Error Analysis on Numerical Integration Algorithms in a Hypoelasticity Framework

This report determines local truncation errors for common stress integration algorithms used in explicit finite element codes with hypoelastic material models. The hypoelastic integration algorithms in question utilize an operator splitting procedure in a rotation neutralized configuration, where the stress response is determined from de- coupling the total deformation into rotational and strain dependent components. This document analyzes the error in evolving the stress given a one-step time increment Δt and compares the errors associated with both the rotational and strain components of the operator splitting method. A slight modification to a traditional algorithm is proposed and studied, where the rate of deformation is appropriately rotated from the midstep configuration at t n+1/2 to the end step configuration at t n+1 before the constitutive evaluation. The proposed modification either completely eliminates the error associated with the rotation rate or is of the same order of magnitude as the original algorithm for the three test cases considered in this report. These cases consist of an unaxial stretch with a constant true strain rate with a rigid body rotation, an uniaxial stretch with a constant engineering strain rate with a rigid body rotation, and a simple shear deformation. All three cases are compared to a closed form solution, and in almost every test case the alternative algorithm yields the most accurate one-step local truncation error.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1987

Stresses and deformations in angle-ply composite tubes

The stress and deformations in angle-ply composite tubes subjected to axisymmetric thermal loading were investigated both experimentally and analytically. For the theoretical portion a generalized plane strain elasticity analysis was developed. The analysis included mechanical and thermal loading, and temperature-dependent material properties. The elasticity analysis was also used to study the effect of including a thin metallic coating on a graphite-epoxy tube. The stresses in the coatings were found to be quite high, exceeding the yield stress of aluminum. An important finding in the analytical studies was the fact that even tubes with a balanced-symmetric lamination sequence exhibit shear deformation, or twist. For the experimental portion an apparatus was developed to measure torsional and axial response in the temperature range of 140 to 360 K. Eighteen specimens were tested, combining three material systems, eight lamination sequences, and three off-axis ply orientation angles. For the twist response, agreement between analysis and experiment was found to be good. The axial response of the tubes tested was found to be greater than predicted by a factor of three. As a result, it is recommended that the thermally induced axial deformations be investigated, both experimentally and analytically.

Rousseau, Carl Q.↗
DOE OSTI2020

Parallel three-dimensional simulations of quasi-static elastoplastic solids

Hypo-elastoplasticity is a flexible framework for modeling the mechanics of many hard materials under small elastic deformation and large plastic deformation. Under typical loading rates, most laboratory tests of these materials happen in the quasi-static limit, but there are few existing numerical methods tailor-made for this physical regime. Here, we extend to three dimensions a recent projection method for simulating quasi-static hypo-elastoplastic materials. The method is based on a mathematical correspondence to the incompressible Navier–Stokes equations, where the projection method of Chorin (1968) is an established numerical technique. We develop and utilize a three-dimensional parallel geometric multigrid solver employed to solve a linear system for the quasi-static projection. Our method is tested through simulation of three-dimensional shear band nucleation and growth, a precursor to failure in many materials. As an example system, we employ a physical model of a bulk metallic glass based on the shear transformation zone theory, but the method can be applied to any elastoplasticity model. We consider several examples of three-dimensional shear banding, and examine shear band formation in physically realistic materials with heterogeneous initial conditions under both simple shear deformation and boundary conditions inspired by friction welding.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1975

Analysis of bonded joints

A refined elastic analysis of bonded joints which accounts for transverse shear deformation and transverse normal stress was developed to obtain the stresses and displacements in the adherends and in the bond. The displacements were expanded in terms of polynomials in the thicknesswise coordinate; the coefficients of these polynomials were functions of the axial coordinate. The stress distribution was obtained in terms of these coefficients by using strain-displacement and stress-strain relations. The governing differential equations were obtained by integrating the equations of equilibrium, and were solved. The boundary conditions (interface or support) were satisfied to complete the analysis. Single-lap, flush, and double-lap joints were analyzed, along with the effects of adhesive properties, plate thicknesses, material properties, and plate taper on maximum peel and shear stresses in the bond. The results obtained by using the thin-beam analysis available in the literature were compared with the results obtained by using the refined analysis. In general, thin-beam analysis yielded reasonably accurate results, but in certain cases the errors were high. Numerical investigations showed that the maximum peel and shear stresses in the bond can be reduced by (1) using a combination of flexible and stiff bonds, (2) using stiffer lap plates, and (3) tapering the plates.

Srinivas, S.↗
DOE OSTI2022

Sand Migration Simulation during Gas Production from Gas Hydrate Reservoir at Kuparuk 7–11–12 site in the Prodhoe Bay Unit, Alaska

Uncontrolled sand production impedes continuous gas production from a hydrate reservoir as observed in the past field-scale gas production tests. Sand mobilization is strongly linked with sediment deformation and high pressure gradient. Sand production occurs when the mobilized sands reach the well. Throughout gas production from a hydrate reservoir, because of hydrate dissociation, deformation and pressure gradient change both in time and space. Therefore, it is necessary to consider the entire process to identify where sand is likely mobilized and how much mobilized sand could reach the well. Here, this study utilizes a coupled thermal, hydrological, chemical, and mechanical (thermo-hydro-chemo-mechanical) sand migration model to simulate a year long gas production from the hydrate reservoir at the Kuparuk 7–11–12 site in the Prodhoe Bay Unit, Alaska. It is found that sand would mainly come from the lower portion of the production zone where faster hydrate dissociation occurs. The relatively faster hydrate dissociation coupled with the fact that hydrate-bearing sediments and hydrate-free sediments have similar stiffness but different strengths causes complex stress transfer, which results in excessive shear deformation when the upper portion of the production zone starts to dissociate. This is evident not only near the well but also away from the well, leading to a large amount of sand mobilization. Another focus of the modeling study is evaluation of the effect of sand migration on gas production. The comparison of two extreme cases–with sandscreen and without sandscreen–suggests that migrated and settled solids around a sand control device can prevent pressure drop across the near-wellbore zone and lead to reduction of the gas production rate.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2004

Combined In-Plane and Through-the-Thickness Analysis for Failure Prediction of Bolted Composite Joints

Although two-dimensional methods provide accurate predictions of contact stresses and bolt load distribution in bolted composite joints with multiple bolts, they fail to capture the effect of thickness on the strength prediction. Typically, the plies close to the interface of laminates are expected to be the most highly loaded, due to bolt deformation, and they are usually the first to fail. This study presents an analysis method to account for the variation of stresses in the thickness direction by augmenting a two-dimensional analysis with a one-dimensional through the thickness analysis. The two-dimensional in-plane solution method based on the combined complex potential and variational formulation satisfies the equilibrium equations exactly, and satisfies the boundary conditions and constraints by minimizing the total potential. Under general loading conditions, this method addresses multiple bolt configurations without requiring symmetry conditions while accounting for the contact phenomenon and the interaction among the bolts explicitly. The through-the-thickness analysis is based on the model utilizing a beam on an elastic foundation. The bolt, represented as a short beam while accounting for bending and shear deformations, rests on springs, where the spring coefficients represent the resistance of the composite laminate to bolt deformation. The combined in-plane and through-the-thickness analysis produces the bolt/hole displacement in the thickness direction, as well as the stress state in each ply. The initial ply failure predicted by applying the average stress criterion is followed by a simple progressive failure. Application of the model is demonstrated by considering single- and double-lap joints of metal plates bolted to composite laminates.

Kradinov, V.↗
DOE OSTI2020

Medium amplitude parallel superposition (MAPS) rheology. Part 2: Experimental protocols and data analysis

An experimental protocol is developed to directly measure the new material functions revealed by medium amplitude parallel superposition (MAPS) rheology. This protocol measures the medium amplitude response of a material to a simple shear deformation composed of three sine waves at different frequencies, revealing a rich dataset consisting of up to 19 measurements of the third-order complex modulus at distinct three-frequency coordinates. We discuss how the choice of input frequencies influences the features of the MAPS domain studied by the experiment. A polynomial interpolation method for reducing the bias of measured values from spectral leakage and reducing variance due to noise is discussed, including a derivation of the optimal range of amplitudes for the input signal. This leads to the conclusion that conducting the experiment in a stress-controlled fashion possesses a distinct advantage to the strain-controlled mode. The experimental protocol is demonstrated through measurements of the MAPS response of a model complex fluid: a surfactant solution of wormlike micelles. The resulting dataset is indeed large and feature-rich, while still acquired in a time comparable to similar medium amplitude oscillatory shear (MAOS) experiments. We demonstrate that the data represent measurements of an intrinsic material function by studying its internal consistency, compatibility with low-frequency predictions for Coleman–Noll simple fluids, and agreement with data obtained via MAOS amplitude sweeps. Finally, the data are compared to predictions from the corotational Maxwell model to demonstrate the power of MAPS rheology in determining whether a constitutive model is consistent with a material’s time-dependent response.

Lennon, Kyle R. (ORCID:0000000212515461)↗
NASA NTRS1983

Measurement of the linear viscoelastic behavior of antimisting kerosene

Measurements of dynamic viscoelastic properties in very small oscillating shear deformations was made on solutions of a jet fuel, Jet A, containing an antimisting polymeric additive, FM-9. A few measurements were also made on solutions of FM-9 in a mixed solvent of mineral oil, Tetralin, and 0-terphenyl. Two samples of FM-9 had approximate number-average molecular weights of 12,000,000 and 8,100,000 as deduced from analysis of the measurements. The ranges of variables were 2.42 to 4.03 g/1 in concentration (0.3 to 0.5% by weight), 1 to 35 in temperature, 1.3 to 9.4 cp in solvent viscosity, and 103 to 6100 Hz in frequency. Measurements in the Jet A solvent were made both with and without a modifying carrier. The results were compared with the Zimm theory and the viscoelastic behavior was found to resemble rather closely that of ordinary non-polar polymers in theta solvents. The relation of the results to the antithixotropic behavior of such solutions a high shear rates is discussed in terms of intramolecular and intermolecular interactions.

Ferry, J. D.↗
DOE OSTI2021

A theory of localized excitations in supercooled liquids

Here, a new connection between the structure and dynamics in glass-forming liquids is presented. We show how the origin of spatially localized excitations, as defined by the dynamical facilitation (DF) theory, can be understood from a structure-based framework. This framework is constructed by associating excitation events in the DF theory to hopping events between energy minima in the potential energy landscape (PEL). By reducing the PEL to an equal energy well picture and applying a harmonic approximation, we develop a field theory to describe elastic fluctuations about inherent states, which are energy minimizing configurations of the PEL. We model an excitation as a shear transformation zone (STZ) inducing a localized pure shear deformation onto an inherent state. We connect STZs to T1 transition events that break the elastic bonds holding the local structure of an inherent state. A formula for the excitation energy barrier, denoted as J σ , is obtained as a function of inherent-state elastic moduli and the radial distribution function. The energy barrier from the current theory is compared to the one predicted by the DF theory where good agreement is found in various two-dimensional continuous poly-disperse atomistic models of glass formers. These results strengthen the role of structure and elasticity in driving glassy dynamics through the creation and relaxation of localized excitations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

High speed manufacturing of aluminum alloy 7075 tubing by Shear Assisted Processing and Extrusion (ShAPE)

Shear assisted processing and extrusion (ShAPE) was used to extrude aluminum alloy 7075 tubing at speeds up to 12.2 m/min. This is in comparison to 2.0 m/min which is generally the limit for conventionally extruded 7075. The increased speed is primarily attributed to more extensive shear deformation, compared to conventional extrusion, which results in a high density of low angle grain boundaries that facilitate continued deformation and delay the onset of surface tearing. Mechanical testing after heat treating to the T6 condition gave an ultimate tensile strength of 565.3 ± 4.6 MPa, yield strength of 495.7 ± 8.7 MPa, and elongation of 16.8 ± 0.8%. Strength values exceed the ASTM minimum standard and are on par with ASM typical values, while elongation was substantially improved compared to 7 and 11% for the ASTM and ASM values respectively. Finally, it was observed that low temperature extrusion at 341 °C and 40 rpm gave superior material properties in the T6 condition compared to high temperature extrusion at 441 °C and 120 rpm due to variances in nanoscale second phase size and distribution.

36 MATERIALS SCIENCE↗