Search NASA⌕ Search

SEARCH · Search NASA

Results for “large material deformation and rotation”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Upcycling of Mixed Aluminum Alloy Shredder Scrap using Shear Processing

Conservation of critical materials is an increasing area of focus in the Unites States. In 2023, aluminum was added to the US Department of Energy Final Critical Materials List which has spurred public and private research into sustainable management of these resources. Additionally, efficiency in manufacturing and conservation of natural resources are growing concerns with targets to lower global carbon emissions, as primary aluminum alloy production is energy intensive, requiring 14 MWh of electricity plus 0.4 tonnes of CO2 per tonne of Al. Due to these factors, is essential that more sustainable manufacturing methods for aluminum alloys are developed going forward. To this end, much research is ongoing on the topic of more efficient utilization and recovery. However, most of this research still requires primary aluminum in the production process. Here, it will be attempted to bypass the use of primary aluminum and produce useful material recycled from 100% post-consumer scrap. Even considering recent developments in recycling of Al scrap, there is still a large amount of post-consumer scrap that is underutilized due to high impurity content, and that amount will increase significantly as more and more aluminum alloys are utilized in vehicles. This “scrap wave” is expected to cover 80% of the demand for automotive aluminum alloys by 2050 . A challenge to be addressed before the coming scrap wave can be fully utilized is that the tolerance of manufacturing techniques to impurities or off-spec alloy compositions must be increased. Particularly, in 5000-and 6000-series alloys (the most common wrought alloys in durable products), excess iron, copper, and silicon create brittle intermetallics during casting that remain in the extruded microstructure which limit the formability, ductility, and corrosion resistance of the alloy. Concerningly, many of the highest-volume post-consumer aluminum scrap streams such as automotive shredder scrap contain a mix of alloys including both wrought and cast alloys. Their compositions can vary widely depending on geography and the time of year. Because they are mixed, they often contain high content of multiple alloying elements such as Si and Cu in higher concentrations than are found in typical wrought alloys. They may also be contaminated with non-Al alloys from fasteners that get mixed in and often have high content of unwanted elements such as Fe. As a method for utilizing these underused scrap streams that are high in tramp elements, an emerging extrusion technology is being developed at the Pacific Northwest National Laboratory (PNNL) that aims to upcycle 100% post-consumer aluminum scrap directly into extruded components without the addition of primary aluminum . This new technology, called Shear Assisted Processing and Extrusion (ShAPE), is enabling a shift away from today’s recycling paradigm by reaching deeper into lower-value scrap streams, using shredder scrap as the extrusion billet material. Sometimes referred to as Twitch or Tweak, these scrap streams result from shredding and sorting of automobiles, building materials, appliances, and consumer goods. ShAPE combines the linear axis of conventional extrusion with a rotating extrusion die. This rotating die applies large strain to the material during extrusion, which breaks up large impurity-containing intermetallic particles, reducing their deleterious effects. This has been demonstrated for 6063 machining scrap spiked with excess Fe, and Twitch scrap high in Fe, Si and Cu where strength and ductility were retained for both feedstock compositions. Additionally, the extreme plastic deformation during ShAPE enables extrusion of billets with high Si that are too brittle for processing by conventional extrusion. By using 100% post-consumer shedder scrap as feedstock, ShAPE has the potential to slash embodied energy and carbon in extruded components by >80% compared to conventional extrusion of primary aluminum alloys.

Milligan, Brian K.↗
DOE OSTI2020

Creep Deformation and Dynamic Grain Growth in an Interstitial-Free Steel

Dynamic grain growth is demonstrated to be much faster than static grain growth in a body-centered-cubic, interstitial-free steel sheet material at 850°C. Dynamic grain growth occurs during concurrent plastic deformation at elevated temperature, whereas static grain growth occurs during static annealing. Grain growth during steady-state plastic flow in tension at 850°C to a true strain of 0.2 at a true-strain rate of 10 -4 s -1 doubled grain size, while static annealing for the same time produced no increase in grain size. This is described as dynamic normal grain growth (DNGG) because no abnormally large grains were observed. The recrystallized microstructure of the steel demonstrated a log-normal distribution of grain sizes. DNGG produced bimodal grain size distributions that deviate from the theoretical expectation of a simple shift to larger sizes during normal growth. The bimodal distributions contained a remnant of small grains that were not consumed during grain growth. DNGG produced a crystallographic texture that is unique from both the recrystallized material and that produced by lattice rotation alone. DNGG strengthened the {111}⟨110⟩ and {111}⟨112⟩ components of the strong γ-fiber component in the original recrystallization texture. Lattice rotation from tensile deformation, by contrast, strengthened the α-fiber components that intersect the original γ-fiber.

36 MATERIALS SCIENCE↗
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.↗
DOE OSTI2026

Controlling Deformation in Al/Ti: How Interface Roughness and Orientation Drive Bimetal Mechanics

The microstructural characteristics and morphology of interfaces in metals can be crucial in governing the initiation of plasticity and early deformation mechanisms under extreme conditions. During high-strain-rate deformation, these interfaces significantly affect dislocation nucleation, twinning, and other mechanisms that directly impact material strength and failure. Despite their importance, a substantial knowledge gap remains between the observed macroscopic material behavior and the underlying role of bimetal interfaces in plasticity initiation. Here, to address this gap, large-scale molecular dynamics simulations are performed on Al/Ti bimetal structures to examine the effect of interface characteristics on the onset of plasticity under uniaxial compression. Specifically, this study investigates how interface roughness (flat vs waveform interfaces) modifies the initiation of plastic events in the microstructure. Atomistic simulations indicate that interface roughness (a microscopic behavior) reduces the stress required for dislocation nucleation, thereby reducing the peak stress relative to a flat interface. For the square interface, the step height strongly influences plasticity initiation by setting the separation of locally flat regions. Varying the interface rotation relative to the loading direction (macroscopic behavior) reveals that the peak stress for both flat and waveform interfaces initially decreases and then increases with rotation, accompanied by a shift in dominant mechanisms—from Al twinning (flat) to interface sliding and Ti-dominated phase transformation and twinning at lower angles. The change in the stress-strain slope during initial compression reflects how rotation alters the resolved shear stress and activates different slip systems. Overall, the study provides valuable insights into the role of bimetal interfaces in controlling plasticity initiation and early deformation pathways for Al/Ti under extreme loading conditions.

36 MATERIALS SCIENCE↗
NASA NTRS2013

Jet Engine Bird Ingestion Simulations: Comparison of Rotating to Non-Rotating Fan Blades

Bird strike events in commercial airliners are a fairly common occurrence. According to data collected by the US Department of Agriculture, over 80,000 bird strikes were reported in the period 1990-2007 in the US alone [1]. As a result, bird ingestion is an important factor in aero engine design and FAA certification. When it comes to bird impacts on engine fan blades, the FAA requires full-scale bird ingestion tests on an engine running at full speed to pass certification requirements. These rotating tests are complex and very expensive. To reduce development costs associated with new materials for fan blades, it is desirable to develop more cost effective testing procedures than full-scale rotating engine tests for material evaluation. An impact test on a non-rotating single blade that captures most of the salient physics of the rotating test would go a long way towards enabling large numbers of evaluative material screening tests. NASA Glenn Research Center has been working to identify a static blade test procedure that would be effective at reproducing similar results as seen in rotating tests. The current effort compares analytical simulations of a bird strike on various nonrotating blades to a bird strike simulation on a rotating blade as a baseline case. Several different concepts for simulating the rotating loads on a non-rotating blade were analyzed with little success in duplicating the deformation results seen in the rotating case. The rotating blade behaves as if it were stiffer than the non-rotating blade resulting in less plastic deformation from a given bird impact. The key factor limiting the success of the non-rotating blade simulations is thought to be the effect of gyroscopics. Prior to this effort, it was anticipated the difficulty would be in matching the pre-stress in the blade due to centrifugal forces Additional work is needed to verify this assertion, and to determine if a static test procedure can simulate the gyroscopic effects in a suitable manner. This paper describes the various non-rotating concepts analyzed, and demonstrates the effect believed to be gyroscopic in nature on the results.

Howard, Samuel A.↗
DOE OSTI2024

Hydrogen effects on the deformation and slip localization in a single crystal austenitic stainless steel

Hydrogen is known to embrittle austenitic stainless steels, which are widely used in high-pressure hydrogen storage and delivery systems, but the mechanisms that lead to such material degradation are still being elucidated. The current work investigates the deformation behavior of single crystal austenitic stainless steel 316L through combined uniaxial tensile testing, characterization and atomistic simulations. Thermally precharged hydrogen is shown to increase the critical resolved shear stress (CRSS) without previously reported deviations from Schmid’s law. Molecular dynamics simulations further expose the statistical nature of the hydrogen and vacancy contributions to the CRSS in the presence of alloying. Slip distribution quantification over large in-plane distances (> 1 mm), achieved via atomic force microscopy (AFM), highlights the role of hydrogen increasing the degree of slip localization in both single and multiple slip configurations. The most active slip bands accumulate significantly more deformation in hydrogen precharged specimens, with potential implications for damage nucleation. For $\langle$110$\rangle$ tensile loading, slip localization further enhances the activity of secondary slip, increases the density of geometrically necessary dislocations and leads to a distinct lattice rotation behavior compared to hydrogen-free specimens, as evidenced by electron backscatter diffraction (EBSD) maps. Finally, the results of this study provide a more comprehensive picture of the deformation aspect of hydrogen embrittlement in austenitic stainless steels.

36 MATERIALS SCIENCE↗
NASA NTRS2013

Jet Engine Bird Ingestion Simulations: Comparison of Rotating to Non-Rotating Fan Blades

Bird strike events in commercial airliners are a fairly common occurrence. According to data collected by the US Department of Agriculture, over 80,000 bird strikes were reported in the period 1990 to 2007 in the US alone (Ref. 1). As a result, bird ingestion is an important factor in aero engine design and FAA certification. When it comes to bird impacts on engine fan blades, the FAA requires full-scale bird ingestion tests on an engine running at full speed to pass certification requirements. These rotating tests are complex and very expensive. To reduce development costs associated with new materials for fan blades, it is desirable to develop more cost effective testing procedures than full-scale rotating engine tests for material evaluation. An impact test on a nonrotating single blade that captures most of the salient physics of the rotating test would go a long way towards enabling large numbers of evaluative material screening tests. NASA Glenn Research Center has been working to identify a static blade test procedure that would be effective at reproducing similar results as seen in rotating tests. The current effort compares analytical simulations of a bird strike on various non-rotating blades to a bird strike simulation on a rotating blade as a baseline case. Several different concepts for simulating the rotating loads on a non-rotating blade were analyzed with little success in duplicating the deformation results seen in the rotating case. The rotating blade behaves as if it were stiffer than the non-rotating blade resulting in less plastic deformation from a given bird impact. The key factor limiting the success of the non-rotating blade simulations is thought to be the effect of gyroscopics. Prior to this effort, it was anticipated the difficulty would be in matching the prestress in the blade due to centrifugal forces Additional work is needed to verify this assertion, and to determine if a static test procedure can simulate the gyroscopic effects in a suitable manner. This paper describes the various non-rotating concepts analyzed, and demonstrates the effect believed to be gyroscopic in nature on the results

Howard, Samuel A.↗
NASA NTRS1981

A hybrid state vector approach to aeroelastic analysis with application to composite lifting surfaces

A computational technique has been developed for performing preliminary design aeroelastic analyses of large aspect ratio lifting surfaces. This technique, applicable to both fixed and rotating wing configurations, is based upon a formulation of the structural equilibrium equations in terms of a hybrid state vector containing generalized force and displacement variables. An integrating matrix is employed to solve these equations for divergence and flutter eigenvalues and steady aeroelastic deformation. Results are presented for simple examples which verify the technique and demonstrate how it can be applied to analyze lifting surfaces, including those constructed from composite materials.

Lehman, L. L.↗
DOE OSTI2022

Nontrivial nanostructure, stress relaxation mechanisms, and crystallography for pressure-induced Si-I → Si-II phase transformation

Crystallographic theory based on energy minimization suggests austenite-twinned martensite interfaces with specific orientation, which are confirmed experimentally for various materials. Pressure-induced phase transformation (PT) from semiconducting Si-I to metallic Si-II, due to very large and anisotropic transformation strain, may challenge this theory. Here, unexpected nanostructure evolution during Si-I → Si-II PT is revealed by combining molecular dynamics (MD), crystallographic theory, generalized for strained crystals, and in situ real-time Laue X-ray diffraction (XRD). Twinned Si-II, consisting of two martensitic variants, and unexpected nanobands, consisting of alternating strongly deformed and rotated residual Si-I and third variant of Si-II, form {111} interface with Si-I and produce almost self-accommodated nanostructure despite the large transformation volumetric strain of -0.237. The interfacial bands arrest the {111} interfaces, leading to repeating nucleation-growth-arrest process and to growth by propagating {110} interface, which (as well as {111} interface) do not appear in traditional crystallographic theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1990

Stress field rotation or block rotation: An example from the Lake Mead fault system

The Coulomb criterion, as applied by Anderson (1951), has been widely used as the basis for inferring paleostresses from in situ fault slip data, assuming that faults are optimally oriented relative to the tectonic stress direction. Consequently if stress direction is fixed during deformation so must be the faults. Freund (1974) has shown that faults, when arranged in sets, must generally rotate as they slip. Nur et al., (1986) showed how sufficiently large rotations require the development of new sets of faults which are more favorably oriented to the principal direction of stress. This leads to the appearance of multiple fault sets in which older faults are offset by younger ones, both having the same sense of slip. Consequently correct paleostress analysis must include the possible effect of fault and material rotation, in addition to stress field rotation. The combined effects of stress field rotation and material rotation were investigated in the Lake Meade Fault System (LMFS) especially in the Hoover Dam area. Fault inversion results imply an apparent 60 degrees clockwise (CW) rotation of the stress field since mid-Miocene time. In contrast structural data from the rest of the Great Basin suggest only a 30 degrees CW stress field rotation. By incorporating paleomagnetic and seismic evidence, the 30 degrees discrepancy can be neatly resolved. Based on paleomagnetic declination anomalies, it is inferred that slip on NW trending right lateral faults caused a local 30 degrees counter-clockwise (CCW) rotation of blocks and faults in the Lake Mead area. Consequently the inferred 60 degrees CW rotation of the stress field in the LMFS consists of an actual 30 degrees CW rotation of the stress field (as for the entire Great Basin) plus a local 30 degrees CCW material rotation of the LMFS fault blocks.

Ron, Hagai↗
NASA NTRS1990

Multiple strike slip faults sets: A case study from the Dead Sea transform

In many strike slip tectonic settings, large rotations of crust blocks about vertical axes have been inferred from paleomagnetic data. These blocks are bounded by sets of parallel faults which presumably accommodate the relative motion between the blocks as regional deformation progress. A mechanical model by Nur et al., (1986) suggests that rotations greater than phi sub c equals 25 to 45 degrees must be accommodated by more than one set of faults, with angle phi sub c between their direction; consequently the sum of the angles between sets must be roughly equal to the total tectonic material rotation. To test this model, the authors investigated the fault geometry and field relation of fault sets in the Mt. Hermon area in northern Israel, where paleomagnetic declination implies data 69 degrees plus or minus 13 degrees counter-clockwise block rotation. The statistical and field relation analysis of over 315 faults shows that the faulting is predominantly right lateral strike slip consisting of three distinct sets. The oldest set strikes 253 degrees, the second oldest set strikes 293 degrees and the youngest strikes 339 degrees. This last direction is consistent also with the current north-south direction of the maximum principle stress axis. The angle phi sub c between the first and second sets is 39 degrees and between the second and third sets 46 degrees, in good agreement with the phi sub c angle predicted from mechanical considerations. The sum of the two angles is 85 degrees, in good agreement with the 69 degrees plus or minus 13 degrees CCW paleomagnetically derived rotation. The results suggest specifically that the sequential development of multiple intersecting fault sets is responsible for the faulting in the Mt. Hermon area; and generally that the model of block rotation with multiple faults provides very good simple rules for analyzing very complex fault patterns.

Ron, Hagai↗
DOE OSTI2024

Anisotropic neutron star crust, solar system mountains, and gravitational waves

"Mountains" or non-axisymmetric deformations of rotating neutron stars (NS) efficiently radiate gravitational waves (GW). We consider analogies between NS mountains and surface features of solar system bodies. Both NS and moons such as Europa or Enceladus have thin crusts over deep oceans while Mercury has a thin crust over a large metallic core. Thin sheets may wrinkle in universal ways. Europa has linear features, Enceladus has "Tiger" stripes, and Mercury has lobate scarps. NS may have analogous features. The innermost inner core of the Earth is anisotropic with a shear modulus that depends on direction. If NS crust material is also anisotropic this will produce an ellipticity, when the crust is stressed, that grows with spin frequency. This yields a braking index (log derivative of spin down rate assuming only GW spin down) very different from n=5 and could explain the maximum spin observed for neutron stars and a possible minimum ellipticity of millisecond pulsars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Large rotations of the grain-scale stress tensor during yielding set the stage for failure

The stress-state within individual grains in a polycrystal determine the fate of the aggregate including mechanical failure. By tracking the evolution of the stress tensor throughout the elastoplastic transition, large rotations of the stress-state, which have long been theorized to occur, are observed experimentally for the first time. These stress rotations (~15°) are more than an order of magnitude larger than the concomitant crystallographic lattice reorientations (~0.9°) well-known to occur during metal plasticity. Furthermore, these rotations are accompanied by a decrease in stress triaxiality within certain grains, promote strain softening, and set the stage for failure at an early stage of deformation. Finally, these results provide a completely new perspective through which to contemplate the question of “hot-spots” responsible for failure of high-performance structural materials.

36 MATERIALS SCIENCE↗
NASA NTRS2016

A Continuum Damage Mechanics Model to Predict Kink-Band Propagation Using Deformation Gradient Tensor Decomposition

A new model is proposed that represents the kinematics of kink-band formation and propagation within the framework of a mesoscale continuum damage mechanics (CDM) model. The model uses the recently proposed deformation gradient decomposition approach to represent a kink band as a displacement jump via a cohesive interface that is embedded in an elastic bulk material. The model is capable of representing the combination of matrix failure in the frame of a misaligned fiber and instability due to shear nonlinearity. In contrast to conventional linear or bilinear strain softening laws used in most mesoscale CDM models for longitudinal compression, the constitutive response of the proposed model includes features predicted by detailed micromechanical models. These features include: 1) the rotational kinematics of the kink band, 2) an instability when the peak load is reached, and 3) a nonzero plateau stress under large strains.

Bergan, Andrew C.↗
DOE OSTI2024

Mechanical properties and deformation mechanisms of single crystal Mg micropillars subjected to high-strain-rate C-axis compression

Here, the mechanical properties and deformation mechanisms of single crystal magnesium under c-axis quasi-static and high-strain rate compressions are investigated through in situ scanning electron microscope (SEM) experiments and post-mortem transmission electron microscope (TEM) characterization. The findings revealed that ductility and high rates of hardening are preserved for pillars as large as 15 μm. Furthermore, rate effects result in a mild increase in flow stress with plastic deformations controlled primarily by the slip of type dislocations. Importantly and in contrast to other literature reports, plastic deformation occurs in the absence of twining. As the strain increases and plastic deformation exceeds about 4%, crystal rotation activates basal slip, <$\mathrm{a}$> type dislocations, resulting in a more rate independent flow stress. TEM observation on micropillars compressed at a strain rate of 250/s, revealed the activation of {${11}$$\bar{2}$$\bar{2}$} < $\bar{1}$$\bar{1}23$ > slip systems and high mobility of screw dislocations as major contributors to plastic strains in excess of 10% without fracture. These findings are relevant to the design of lightweight materials used in transportation systems, e.g., selection of material grain size. Moreover, the experimental data here reported provides the materials science community with a unique opportunity to validate discrete dislocation dynamics (DDD) formulations employed in multiscale design of materials.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Sub-ply level scaling approach investigated for graphite-epoxy composite beam columns

Scale model graphite-epoxy composite specimens were fabricated using the 'sub-ply level' approach and tested as beam-columns under an eccentric axial load to determine the effect of specimen size on flexural response and failure. In the current research project, although the fiber diameters are not scaled, the thickness of the pre-preg material itself has been scaled by adjusting the number of fibers through the thickness of a single ply. Three different grades of graphite-epoxy composite material (AS4/3502) were obtained from Hercules, Inc., in which the number of fibers through the thickness of a single ply was reduced (Grade 190 with 12 to 16 fibers, Grade 95 with 6 to 8 fibers, and Grade 48 with 3 to 4 fibers). Thus, using the sub-ply level approach, a baseline eight ply quasi-isotropic laminate could be fabricated using either the Grade 48 or Grade 95 material and the corresponding full-scale laminate would be constructed from Grade 95 or standard Grade 190 material, respectively. Note that in the sub-ply level approach, the number of ply interfaces is constant for the baseline and full-scale laminates. This is not true for the ply level and sublaminate level scaled specimens. The three grades of graphite-epoxy composite material were used to fabricate scale model beam-column specimens with in-plane dimensions of 0.5*n x 5.75*n, where n=1,2,4 corresponsing to 1/4, 1/2, and full-scale factors. Angle ply, cross ply, and quasi-isotropic laminate stacking sequences were chosen for the investigation and the test matrices for each laminate type are given. Specimens in each laminate family with the same in-plane dimensions but different thicknesses were tested to isolate the influence of the thickness dimension on the flexural response and failure. Also, specific lay-ups were chosen with blocked plies and dispersed plies for each laminate type. Specimens were subjected to an eccentric axial load until failure. The load offset was introduced through a set of hinges which were attached to the platens of a standard load test machine. Three sets of geometrically scaled hinges were used to ensure that scaled loading conditions were applied. This loading condition was chosen because it promotes large flexural deformations and specimens fail at the center of the beam, away from the grip supports. Five channels of data including applied vertical load, end shortening displacement, strain from gages applied back-to-back at the midspan of the beam, and rotation of the hinge from a bubble inclinometer were recorded for each specimen. The beam-column test configuration was used previously to study size effects in ply level scaled composite specimens of the same material system, sizes, and stacking sequences. Thus, a direct comparison between the two scaling approaches is possible. Ply level scaled beam-columns with angle ply, cross ply, and quasi-isotropic lay-ups exhibited no size dependencies in the flexural response, but significant size effects in strength. The reduction in strength with increasing specimen size was not predicted successfully by analysis techniques. It is anticipated that results from this investigation will lead to a better understanding of the strength scale effect in composite structures.

Jackson, Karen E.↗
DOE OSTI2021

Domain patterns and super-elasticity of freestanding BiFeO 3 membranes via phase-field simulations

Super-elasticity of functional ferroelectric oxides offers promises for integrating ferroelectric films into flexible electronics. However, super-elastic deformation is a complex phenomenon related to possibly multiple concurrent mechanisms. Fundamentally understanding how multiple mechanisms contribute to the super-elasticity of ferroelectric oxides is crucial to realizing their potential flexible electronic applications. In this study, we employ phase-field simulations to model the dynamics of ferroelectric domain patterns of freestanding BiFeO 3 membranes to understand the origin of their super-elasticity under substantial bending deformation (5% strain). It is demonstrated that both a reversible Rhombohedral-Tetragonal (R-T) phase transition and a nearly reversible domain evolution of BiFeO 3 membranes contribute to accommodating the large deformation and thus their super-elasticity. The dynamics of domain evolution also reveal the formation of an exotic ferroelectric vortex and polarization rotation before the phase transition. We constructed a diagram of phases and domain patterns as a function of the membrane thickness and bending angle, which allows one to readily predict the emergence of T phase and ferroelectric vortex in bent BFO membranes. These results not only provide fundamental understanding of mesoscale super-elastic mechanisms but also reveal exotic domain states of ferroelectric membranes.

36 MATERIALS SCIENCE↗