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 1,135 records · Page 63

A variational phase-field framework for thermal softening and dynamic ductile fracture

A variational phase field model for dynamic ductile fracture is presented. Here, the model is designed for elasto-viscoplastic materials subjected to rapid deformations in which the effects of heat generation and material softening are dominant. The variational framework allows for the consistent inclusion of plastic dissipation in the heat equation as well as thermal softening. It employs a coalescence function to degrade fracture energy during regimes of high plastic flow. A variationally consistent form of the Johnson–Cook model is developed for use with the framework. Results from various benchmark problems in dynamic ductile fracture are presented to demonstrate capabilities. In particular, the ability of the model to regularize shear band formation and subsequent damage evolution in two- and three-dimensional problems is demonstrated. Importantly, these phenomena are naturally captured through the underlying physics without the need for phenomenological criteria such as stability thresholds for the onset of shear band formation.

42 ENGINEERING↗

Validation of material models for puncture of 7075-T651 aluminum plate

Plate puncture simulations are challenging computational tasks that require advanced material models including high strain rate and thermal-mechanical effects on both deformation and failure, plus finite element techniques capable of representing large deformations and material failure. The focus of this work is on the material issues, which require large sets of experiments, flexible material models and challenging calibration procedures. In this study, we consider the puncture of 12.7 mm thick, 7075-T651 aluminum alloy plates by a cylindrical punch with a hemispherical nose and diameter of 12.7 mm. The plasticity and ductile failure models were isotropic with calibration data obtained from uniaxial tension tests at different temperatures and strain rates plus quasi-static notched tension tests and shear-dominated tests described here. Sixteen puncture experiments were conducted to identify the threshold penetration energy, mode of puncture and punch acceleration during impact, The punch was mounted on a 139 kg mass and dropped on the plates with different impact speeds. Since the mass was the same in all tests, the quantity of interest was the impact speed. The axis and velocity of the punch were perpendicular to the plate surface. The mean threshold punch speed was 3.05 m/s, and the mode of failure was plugging by thermal-mechanical shear banding accompanied by scabbing fragments. Application of the material models in simulations of the tests yielded accurate estimates of the threshold puncture speed and of the mode of failure. Time histories of the punch acceleration compared well between simulation and test. Remarkably, the success of the simulations occurred in spite of even the smallest element used being larger than the width of the shear bands.

36 MATERIALS SCIENCE↗

Utilizing integrated neutron diffraction and elastoplastic self-consistent crystal plasticity model to quantitatively assess the strengthening mechanism in Al–12.5Ce and Al–12.5Ce–0.4Mg alloys

An integrated in-situ neutron diffraction and elastic plastic self-consistent crystal plasticity (EPSC-CP) modeling scheme is performed on a binary Al–12Ce alloy and a ternary Al–12Ce–0.4Mg alloys. Using this scheme, the constitutive parameters, i.e. elastic constants and slip system parameters of individual phases can be calibrated which can be used in microstructure-based CP models to predict materials performance. From this study, it is shown that the elastic constants of Al 11 Ce 3 intermetallics calculated from density function theory calculation in the literature are rather accurate. When applied to the EPSC-CP model, the lattice strains of both the binary and ternary alloys are correctly predicted as compared with experiments, and large lattice strain differences between Al (100) plane and Al 11 Ce 3 (010) plane are demonstrated. The slip system parameters calibrated by the scheme shows that the addition of 0.4 wt% Mg in the alloy has little influence on the critical resolved shear stress of initial dislocation glide in the Al matrix which caused plastic yield in the material. This can be explained by the very dilute Mg solute content in the Al solid solution, causing large spacing of Al–Mg lattice misfit sites and little impact on resistance of dislocation glide at initial yield. The 0.4 wt% Mg addition, on the other hand, has a large influence on the hardening term in the slip system parameters, indicating those Al–Mg misfit sites do help dislocation accumulation during the deformation. The impact of dilute Mg addition on the Al slip system parameters is also reflected in the flow behavior of the ternary alloy: little impact on the yield stress, but a large impact on working hardening and tensile strength of the materials which is consistent with the literature.

36 MATERIALS SCIENCE↗

Grain boundary slip – twin transmission in titanium

Here, using a combination of mechanical testing, scanning electron microscopy, and a unified crystal plasticity framework for discrete intragranular shear localization, we investigate intense, localized slip bands on prismatic planes and {10$\bar{1}$2} $\langle$$\bar{1}$011$\rangle$ tensile twins, and their transmission across the grain boundaries (GBs) in commercially pure titanium. The analyses show that the orientation and curvature of the GB influence the local stress fields at the GBs, and consequently, the slip/twin transmission across the boundary. In addition to host grain properties, neighboring grain properties, such as active slip systems and instances of heterogeneity like slip bands and twins, heavily affect the deformation mechanisms in each grain. Finally, the applicability of geometric factors to predict the transmission in the experimentally observed co-located pairs is discussed. The local stress field calculated by the discrete slip and twin band model is shown to be capable of determining whether a transmission has occurred in an observed co-located pair, and also the direction of the transmission.

36 MATERIALS SCIENCE↗

Atomic-scale homogeneous plastic flow beyond near-theoretical yield stress in a metallic glass

Abstract The onset of yielding and the related atomic-scale plastic flow behavior of bulk metallic glasses at room temperature have not been fully understood due to the difficulty in performing the atomic-scale plastic deformation experiments needed to gain direct insight into the underlying fundamental deformation mechanisms. Here we overcome these limitations by combining a unique sample preparation method with atomic force microscopy-based indentation, which allows study of the yield stress, onset of yielding, and atomic-scale plastic flow of a platinum-based bulk metallic glass in volumes containing as little as approximately 1000 atoms. Yield stresses markedly higher than in conventional nanoindentation testing were observed, surpassing predictions from current models that relate yield stress to tested volumes; subsequent flow was then established to be homogeneous without exhibiting collective shear localization or loading rate dependence. Overall, variations in glass properties due to fluctuations of free volume are found to be much smaller than previously suggested.

Yu, Jiaxin↗

Investigation of the yielding transition in concentrated colloidal systems via rheo-XPCS

We probe the microstructural yielding dynamics of a concentrated colloidal system by performing creep/recovery tests with simultaneous collection of coherent scattering data via X-ray Photon Correlation Spectroscopy (XPCS). This combination of rheology and scattering allows for time-resolved observations of the microstructural dynamics as yielding occurs, which can be linked back to the applied rheological deformation to form structure–property relations. Under sufficiently small applied creep stresses, examination of the correlation in the flow direction reveals that the scattering response recorrelates with its predeformed state, indicating nearly complete microstructural recovery, and the dynamics of the system under these conditions slows considerably. Conversely, larger creep stresses increase the speed of the dynamics under both applied creep and recovery. The data show a strong connection between the microstructural dynamics and the acquisition of unrecoverable strain. By comparing this relationship to that predicted from homogeneous, affine shearing, we find that the yielding transition in concentrated colloidal systems is highly heterogeneous on the microstructural level.

42 ENGINEERING↗

Geodetic Strain Analysis Tool

A geodetic software analysis tool enables the user to analyze 2D crustal strain from geodetic ground motion, and create models of crustal deformation using a graphical interface. Users can use any geodetic measurements of ground motion and derive the 2D crustal strain interactively. This software also provides a forward-modeling tool that calculates a geodetic velocity and strain field for a given fault model, and lets the user compare the modeled strain field with the strain field obtained from the user s data. Users may change parameters on-the-fly and obtain a real-time recalculation of the resulting strain field. Four data products are computed: maximum shear, dilatation, shear angle, and principal components. The current view and data dependencies are processed first. The remaining data products and views are then computed in a round-robin fashion to anticipate view changes. When an analysis or display parameter is changed, the affected data products and views are invalidated and progressively re-displayed as available. This software is designed to facilitate the derivation of the strain fields from the GPS and strain meter data that sample it to facilitate the understanding of the strengths and weaknesses of the strain field derivation from continuous GPS (CGPS) and other geodetic data from a variety of tectonic settings, to converge on the "best practices" strain derivation strategy for the Solid Earth Science ESDR System (SESES) project given the CGPS station distribution in the western U.S., and to provide SESES users with a scientific and educational tool to explore the strain field on their own with user-defined parameters.

Kedar, Sharon↗

Large-scale landslide simulations: Global deformation, velocities and basal friction

The cause of the apparent small friction exhibited by long runout landslides has long been speculated upon. In an attempt to provide some insight into the matter, this paper describes results obtained from a discrete particle computer simulation of landslides composed of up to 1,000,000 two-dimensional discs. While simplified, the results show many of the characteristics of field data (the volumetric effect on runout, preserved strata, etc.) and with allowances made for the two-dimensional nature of the simulation, the runouts compare well with those of actual landslides. The results challenge the current view that landslides travel as a nearly solid block riding atop a low friction basal layer. Instead, they show that the mass is completely shearing and indicate that the apparent friction coefficient is an increasing function of shear rate. The volumetric effect can then be understood. With all other conditions being equal, different size slides appear to travel with nearly the same average velocity; however, as the larger landslides are thicker, they experience smaller shear rates and correspondingly smaller frictional resistance.

Campbell, Charles S.↗

Seismic anisotropy, dominant slip systems and phase transitions in the lowermost mantle

SUMMARY The presence of seismic anisotropy at the base of the Earth's mantle is well established, but there is no consensus on the deformation mechanisms in lower mantle minerals that could explain it. Strong anisotropy in magnesium post-perovskite (pPv) has been invoked, but different studies disagree on the dominant slip systems at play. Here, we aim to further constrain this by implementing the most recent results from atomistic models and high-pressure deformation experiments, coupled with a realistic composition and a 3-D geodynamic model, to compare the resulting deformation-induced anisotropy with seismic observations of the lowermost mantle. We account for forward and reverse phase transitions from bridgmanite (Pv) to pPv. We find that pPv with either dominant (001) or (010) slip can both explain the seismically observed anisotropy in colder regions where downwellings turn to horizontal flow, but only a model with dominant (001) slip matches seismic observations at the root of hotter large-scale upwellings. Allowing for partial melt does not change these conclusions, while it significantly increases the strength of anisotropy and reduces shear and compressional velocities at the base of upwellings.

58 GEOSCIENCES↗

Fatigue Analysis of Notched Laminates: A Time-Efficient Macro-Mechanical Approach

A coupled transversely isotropic deformation and damage fatigue model is implemented within the finite element method and was utilized along with a static progressive damage model to predict the fatigue life, stiffness degradation as a function of number of cycles, and post-fatigue tension and compression response of notched, multidirectional laminates. Initially, the material parameters for the fatigue model were obtained utilizing micromechanics simulations and the provided [0], [90] and [plus or minus 45] experimental composite laminate S-N (stress-cycle) data. Within the fatigue damage model, the transverse and shear properties of the plies were degraded with an isotropic scalar damage variable. The damage in the longitudinal (fiber) ply direction was suppressed, and only the strength of the fiber was degraded as a function of fatigue cycles. A maximum strain criterion was used to capture the failure in each element, and once this criterion was satisfied, the longitudinal stiffness of the element was decreased by a factor of 10 (sup 4). The resulting, degraded properties were then used to calculate the new stress state. This procedure was repeated until final failure of the composite laminate was achieved or a specified number of cycles reached. For post-fatigue tension and compression behavior, four internal state variables were used to control the damage and failure. The predictive capability of the above-mentioned approach was assessed by performing blind predictions of the notched multidirectional IM7/977-3 composite laminates response under fatigue and post-fatigue tensile and compressive loading, followed by a recalibration phase. Although three different multidirectional laminates were analyzed in the course of this study, only detailed results (i.e., stiffness degradation and post-fatigue stress-train curves as well as damage evolution states for a single laminate ([30/60/90/minus 30/minus 60] (sub 2s)) are discussed in detail here.

Composite Fatigue analysis↗

Pinning dislocations in colloidal crystals with active particles that seek stacking faults

There is growing interest in functional, adaptive devices built from colloidal subunits of micron size or smaller. A colloidal material with dynamic mechanical properties could facilitate such microrobotic machines. Here we study via computer simulation how active interstitial particles in small quantities can be used to modify the bulk mechanical properties of a colloidal crystal. Passive interstitial particles are known to pin dislocations in metals, thereby increasing resistance to plastic deformation. We extend this tactic by employing anisotropic active interstitials that travel super-diffusively and bind strongly to stacking faults associated with partial dislocations. We find that: (1) interstitials that are effective at reducing plasticity compromise between strong binding to stacking faults and high mobility in the crystal bulk. (2) Reorientation of active interstitials in the crystal depends upon rotational transitions between high-symmetry crystal directions. (3) The addition of certain active interstitial shapes at concentrations as low as 60 per million host particles (0.006%) can create a shear threshold for dislocation migration. This work demonstrates how active materials in a dense matrix can locally sense their environment and lead to bulk property changes.

36 MATERIALS SCIENCE↗

Nanotwin stability in alloyed copper under ambient and cryo-temperature dependent deformation states

Nanometer size grains, solid solution strengthening, low temperatures, and a high density of nanotwins can provide significant strengthening to a material. Nevertheless, under loading the microstructure may evolve altering the effect of grain size, twin density, etc. on this strengthening. In some cases, these effects can be coupled as in detwinning, which is a mechanism for cryogenic grain growth. Here, in this work, a series of copper nanotwinned columnar grains are solid solution alloyed with aluminum and subject to ambient and cryogenic (–196 °C) environments during mechanical indentation, with the influence of alloying on the detwinning evolution studied. Under indentation, at a lower aluminum solute content (2 at.%), the columnar grains bent corresponding to a loss of twins. At a higher aluminum solute content (8 at.%), the grains shear and retain the nanotwins. Atomistic-scale simulations reveal an increase in the incoherent twin boundary velocity with decreasing temperature, but this boundary velocity is reduced with increasing alloy solute content. The coupling of solid solution strengthening and nanotwin stabilization is then used to explain the microstructural stability as a function of indentation at different temperature conditions.

36 MATERIALS SCIENCE↗

Cyclic axial-torsional deformation behavior of a cobalt-base superalloy

Multiaxial loading, especially at elevated temperature, can cause the inelastic response of a material to differ significantly from that predicted by simple flow rules, i.e., von Mises or Tresca. To quantify some of these differences, the cyclic high-temperature, deformation behavior of a wrought cobalt-based superalloy, Haynes 188, is investigated under combined axial and torsional loads. Haynes 188 is currently used in many aerospace gas turbine and rocket engine applications, e.g., the combustor liner for the T800 turboshaft engine for the RAH-66 Comanche helicopter and the liquid oxygen posts in the main injector of the space shuttle main engine. The deformation behavior of this material is assessed through the examination of hysteresis loops generated from a biaxial fatigue test program. A high-temperature axial, torsional, and combined axial-torsional fatigue data base has been generated on Haynes 188 at 760 C. Cyclic loading tests have been conducted on uniform gauge section tubular specimens in a servohydraulic axial-torsional test rig. Test control and data acquisition were accomplished with a minicomputer. In this paper, the cyclic hardening characteristics and typical hysteresis loops in the axial stress versus axial strain, shear stress versus engineering shear strain, axial strain versus engineering shear strain, and axial stress versus shear stress spaces are presented for cyclic, in-phase and out-of-phase, axial torsional tests. For in-phase tests three different values of the proportionality constant, lambda (ratio of engineering shear strain amplitude to axial strain amplitude), are examined, viz., 0.86, 1.73, and 3.46. In the out-of-phase tests, three different values of the phase angle, phi (between the axial and engineering shear strain waveforms), are studied, viz., 30, 60, and 90 deg with lambda = 1.73. The cyclic hardening behaviors of all the tests conducted on Haynes 188 at 760 C are evaluated using the von Mises equivalent stress-strain and the maximum shear stress-maximum engineering shear strain (Tresca) curves. Comparisons are also made between the hardening behaviors of cyclic axial, torsional, and combined in-phase and out-of-phase axial-torsional fatigue tests. These comparisons are accomplished through simple Ramberg-Osgood type stress-strain functions for cyclic, axial stress-strain and shear stress-engineering shear strain curves.

Bonacuse, Peter J.↗

Phase-Shifting Liquid Crystal Interferometers for Microgravity Fluid Physics

The initial focus of this project was to eliminate both of these problems in the Liquid Crystal Point-Diffraction Interferometer (LCPDI). Progress toward that goal will be described, along with the demonstration of a phase shifting Liquid Crystal Shearing Interferometer (LCSI) that was developed as part of this work. The latest LCPDI, other than a lens to focus the light from a test section onto a diffracting microsphere within the interferometer and a collimated laser for illumination, the pink region contained within the glass plates on the rod-mounted platform is the complete interferometer. The total width is approximately 1.5 inches with 0.25 inches on each side for bonding the electrical leads. It is 1 inch high and there are only four diffracting microspheres within the interferometer. As a result, it is very easy to align, achieving the first goal. The liquid crystal electro-optical response time is a function of layer thickness, with thinner devices switching faster due to a reduction in long-range viscoelastic forces between the LC molecules. The LCPDI has a liquid crystal layer thickness of 10 microns, which is controlled by plastic or glass microspheres embedded in epoxy 'pads' at the corners of the device. The diffracting spheres are composed of polystyrene/divinyl benzene polymer with an initial diameter of 15 microns. The spheres deform slightly when the interferometer is assembled to conform to the spacing produced by the microsphere-filled epoxy spacer pads. While the speed of this interferometer has not yet been tested, previous LCPDIs fabricated at the Laboratory for Laser Energetics switched at a rate of approximately 3.3 Hz, a factor of 10 slower than desired. We anticipate better performance when the speed of these interferometers is tested since they are approximately three times thinner. Phase shifting in these devices is a function of the AC voltage level applied to the liquid crystal. As the voltage increases, the dye in the liquid crystal tends to become more transparent, thus introducing a rather large amount of error into the phase-shifting measurement. While that error can be greatly reduced by normalization, we prefer eliminating the source of the error. To that end, we have pursued development of a 'blend' of custom dyes that will not exhibit these properties. That goal has not yet been fully achieved. Guardalben, et al, presented a similar set of interferograms in a paper partially funded by this grant. Shearing interferometers are a second class of common path interferometers. Typically they consist of a thick glass plate optimized for equal reflection from the front and back surface. While not part of the original thrust of the project, through the course of laboratory work, we demonstrated a prototype of a shearing interferometer capable of phase shifting using a commercial liquid crystal retardation plate. A schematic of this liquid crystal shearing interferometer (LCSI) and a sample set of interferograms are in the reference. This work was also supported by the U.S. Department of Energy Office of Inertial Confinement Fusion under Cooperative Agreement No. DE-FC03-92SF19460, the University of Rochester, and the New York State Energy Research and Development Authority. The support of DOE does not constitute an endorsement by DOE of the views expressed in this article. Additional information is included in the original extended abstract.

Griffin, DeVon W.↗

Simplified solution for stresses and deformations

A shortcut is presented to the classical Hertzian solution for local stress and deformation of two elastic bodies in contact, which employs simplified forms for the ellipticity and for the complete elliptic integrals of the first and second kinds as a function of the geometry. Unlike previously proposed simplified formulas whose range of applicability was limited to ellipticities equal to, or greater than one, the present method extends the range of validity to include ellipticities less than one; in which the semimajor axis in the elliptical contact lies in a direction parallel to the rolling, rather than being perpendicular, as in previous studies. An auxiliary shear stress parameter is also expressed in simplified form as a function of the geometry, permitting a shortcut calculation to be made for the location and magnitude of the maximum subsurface shear stress.

Hamrock, B. J.↗

Bridging microscale to macroscale mechanical property measurements of FeCrAl alloys by crystal plasticity modeling

FeCrAl alloys are candidates for accident tolerant fuel cladding of light water reactors. In this work, a microstructure- and temperature-dependent crystal plasticity model is employed to bridge microscale to macroscale mechanical property measurements of FeCrAl alloys. With the visco-plastic self-consistent (VPSC) polycrystal plasticity framework, a mechanism-based single crystal plasticity (MSCP) model adopts the Arrhenius type rate equation to describe the dependence of the critical resolved shear stress for dislocation slips on their temperature-dependent intrinsic frictional resistance and the microstructure-dependent irradiation hardening. The intrinsic frictional resistance associated with {110}<111> and {112}<111> slip systems were measured by in-situ micromechanical testing on unirradiated/irradiated samples at 25-500 °C. The irradiation hardening is estimated by the Bacon-Kocks-Scattergood (BKS) model with density and size of radiation-induced defects measured from microstructural characterization. Several features associated with thermo-mechanical behavior of unirradiated/irradiated polycrystalline FeCrAl alloys are captured. High density of deformation-induced dislocations and radiation-induced defects results in obvious hardening at room temperature, which is weakened at high temperature, and facilitates damage evolution during deformation. Moreover, both high temperature and radiation-induced defects, which facilitate dislocation multiplication, trigger large hardening rate. Finally, the proposed method together with application of accelerator-based ion irradiation technique is a surrogate approach to simulate neutron damage, improving the efficiency associated with evaluation of mechanical properties of FeCrAl alloys exposed to temperature, stress and radiation conditions.

36 MATERIALS SCIENCE↗

Structural deformation and mechanical response of CrS 2 , CrSe 2 and Janus CrSSe

In the framework of density functional theory (DFT), we investigate the structural deformation, and mechanical behavior of the Janus CrSSe, which has out-of-plane structural asymmetry, with conventional transition metal dichalcogenides (TMDs) CrS 2 and CrSe 2 . The Janus CrSSe could be a potential candidate for machinable optoelectronic and piezoelectric applications. Here, we predict that these compounds are chemically, mechanically, and dynamically stable with the covalent bond between the TM(Cr) and chalcogen(X=S, Se) atoms. Due to the influence of tensile strain, the Cr–X bond length of each monolayers increases and the thickness decreases. Interestingly, the in-plane stiffness, shear and layer moduli, Poisson’s ratio, ultimate bi/uni-axial stress of Janus CrSSe are in between the values of CrS 2 and CrSe 2 monolayers. Similar to TMDs, the orientation-dependent in-plane stiffness and Poisson’s ratio demonstrates the isotropic behavior in Janus CrSSe. Furthermore, it can sustain a larger value of uni/bi-axial tensile strain with the critical strain equivalent to CrX 2 monolayers. By applying higher-order strain, we have also found average elastic–plastic behavior as expected. These findings demonstrate that the Janus CrSSe monolayer is a mechanically stable and ductile compound that maintains the hybrid behavior.

36 MATERIALS SCIENCE↗

A Modular Approach to Model Oscillating Control Surfaces Using Navier Stokes Equations

The use of active controls for rotorcraft is becoming more important for modern aerospace configurations. Efforts to reduce the vibrations of helicopter blades with use of active-controls are in progress. Modeling oscillating control surfaces using the linear aerodynamics theory is well established. However, higher-fidelity methods are needed to account for nonlinear effects, such as those that occur in transonic flow. The aeroelastic responses of a wing with an oscillating control surface, computed using the transonic small perturbation (TSP) theory, have been shown to cause important transonic flow effects such as a reversal of control surface effectiveness that occurs as the shock wave crosses the hinge line. In order to account for flow complexities such as blade-vortex interactions of rotor blades higher-fidelity methods based on the Navier-Stokes equations are used. Reference 6 presents a procedure that uses the Navier-Stokes equations with moving-sheared grids and demonstrates up to 8 degrees of control-surface amplitude, using a single grid. Later, this procedure was extended to accommodate larger amplitudes, based on sliding grid zones. The sheared grid method implemented in EulerlNavier-Stokes-based aeroelastic code ENS AERO was successfully applied to active control design by industry. Recently there are several papers that present results for oscillating control surface using Reynolds Averaged Navier-Stokes (RANS) equations. References 9 and 10 report 2-D cases by filling gaps with overset grids. Reference 9 compares integrated forces with the experiment at low oscillating frequencies whereas Ref. 10 reports parametric studies but with no validation. Reference II reports results for a 3D case by modeling the gap region with a deformed grid and compares force results with the experiment only at the mid-span of flap. In Ref. II grid is deformed to match the control surface deflections at the section where the measurements are made. However, there is no indication in Ref. II that the gaps are explicitly modeled as in Ref. 6. Computations using overset grids are reported in Ref. 12 for a case by adding moving control surface to an existing blade but with no validation either with an experiment or another computation.

Navier↗