An official website of the United States government
Scale-free phase-field approach and corresponding finite element method simulations for multivariant martensitic phase transformation from cubic Si I to tetragonal Si II in a polycrystalline aggregate are presented. Important features of the model are large and very anisotropic transformation strain tensor 𝜺𝑡 = {0.1753; 0.1753; −0.447} and stress-tensor dependent athermal dissipative threshold for transformation, which produce essential challenges for computations. 3D polycrystals with stochastically oriented grains are subjected to uniaxial strain- and stress-controlled loadings under periodic boundary conditions and zero averaged lateral strains. Coupled evolution of discrete martensitic microstructure, volume fractions of martensitic variants and Si II, stress and transformation strain tensors, and texture are presented and analyzed. Macroscopic variables effectively representing multivariant transformational behavior are introduced. Macroscopic stress–strain and transformational behavior for 55 and 910 grains are close. Large transformation strains and grain boundaries lead to huge internal stresses of tens GPa, which affect microstructure evolution and macroscopic behavior. In contrast to a single crystal, the local mechanical instabilities due to phase transformation and negative local tangent modulus are stabilized at the macroscale by arresting/slowing the growth of Si II regions by the grain boundaries. This leads to increasing stress during transformation. The developed methodology can be used for studying similar phase transformations with large transformation strains and for further development by including plastic strain and strain-induced transformations.
more »
« less
Fifth-degree elastic energy for predictive continuum stress–strain relations and elastic instabilities under large strain and complex loading in silicon
Abstract Materials under complex loading develop large strains and often phase transformation via an elastic instability, as observed in both simple and complex systems. Here, we represent a material (exemplified for Si I) under large Lagrangian strains within a continuum description by a 5th-order elastic energy found by minimizing error relative to density functional theory (DFT) results. The Cauchy stress—Lagrangian strain curves for arbitrary complex loadings are in excellent correspondence with DFT results, including the elastic instability driving the Si I → II phase transformation (PT) and the shear instabilities. PT conditions for Si I → II under action of cubic axial stresses are linear in Cauchy stresses in agreement with DFT predictions. Such continuum elastic energy permits study of elastic instabilities and orientational dependence leading to different PTs, slip, twinning, or fracture, providing a fundamental basis for continuum physics simulations of crystal behavior under extreme loading.
more »
« less
- PAR ID:
- 10180665
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- npj Computational Materials
- Volume:
- 6
- Issue:
- 1
- ISSN:
- 2057-3960
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Recent works have shown that in contrast to classical linear elastic fracture mechanics, endowing crack fronts in a brittle Green-elastic solid with Steigmann-Ogden surface elasticity yields a model that predicts bounded stresses and strains at the crack tips for plane-strain problems. However, singularities persist for anti-plane shear (mode-III fracture) under far field loading, even when Steigmann-Ogden surface elasticity is incorporated. This work is motivated by obtaining a model of brittle fracture capable of predicting bounded stresses and strains for all modes of loading. We formulate an exact general theory of a three-dimensional solid containing a boundary surface with strain-gradient surface elasticity. For planar reference surfaces parameterized by flat coordinates, the form of surface elasticity reduces to that introduced by Hilgers and Pipkin, and when the surface energy is independent of the surface covariant derivative of the stretching, the theory reduces to that of Steigmann and Ogden. We discuss material symmetry using Murdoch and Cohen’s extension of Noll’s theory. We present a model small-strain surface energy that incorporates resistance to geodesic distortion, satisfies strong ellipticity, and requires the same material constants found in the Steigmann-Ogden theory. Finally, we derive and apply the linearized theory to mode-III fracture in an infinite plate under far-field loading. We prove that there always exists a unique classical solution to the governing integro-differential equation, and in contrast to using Steigmann-Ogden surface elasticity, our model is consistent with the linearization assumption in predicting finite stresses and strains at the crack tips.more » « less
-
Pressure-induced phase transformations (PTs) in Si, the most important electronic material, have been broadly studied. However, strain-induced PTs in Si were never studied in situ. Here, we revealed in situ various important plastic strain-induced PT phenomena. A correlation between the particle size's direct and inverse Hall-Petch effect on yield strength and pressure for strain-induced PT is found. For 100 nm particles, strain-induced PT Si-I³Si-II initiates at 0.3 GPa versus 16.2 GPa under hydrostatic conditions; Si-I³Si-III PT starts at 0.6 GPa and does not occur under hydrostatic pressure. Pressure in small Si-II and Si-III regions is ~5-7 GPa higher than in Si-I. Retaining Si-II and single-phase Si-III at ambient pressure and obtaining reverse Si-II³Si-I PT demonstrates the possibilities of manipulating different synthetic paths. The obtained results corroborate the elaborated dislocation pileup-based mechanism and have numerous applications for developing economic defect-induced synthesis of nanostructured materials, surface treatment (polishing, turning, etc.), and friction.more » « less
-
Pressure-induced phase transformations (PTs) between numerous phases of Si, the most important electronic material, have been studied for decades. This is not the case for plastic strain-induced PTs. Here, we revealed in-situ various unexpected plastic strain-induced PT phenomena. Thus, for 100 nm Si, strain-induced PT Si-I to Si-II (and Si-I to Si-III) initiates at 0.4 GPa (0.6 GPa) versus 16.2 GPa (∞, since it does not occur) under hydrostatic conditions; for 30 nm Si, it is 6.1 GPa versus ∞. The predicted theoretical correlation between the direct and inverse Hall-Petch effect of the grain size on the yield strength and the minimum pressure for strain-induced PT is confirmed for the appearance of Si-II. Retaining Si-II at ambient pressure and obtaining reverse Si-II to Si-I PT are achieved, demonstrating the possibilities of manipulating different synthetic paths.more » « less
-
Abstract 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.more » « less
