SEARCH · Engineering Papers

Results for “shear modulus”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Mechanical Properties of a Solvated Biomolecule: RGD (1FUV) Peptide

The mechanical properties of proteins/peptides play an essential role in their functionalities and implications, as well as their structure and dynamic properties. Understanding mechanical properties is pivotal to our knowledge of protein folding and the molecular basis of diverse cellular processes. Herein, we present a computational approach using ab initio quantum mechanical calculations to determine the mechanical properties—such as bulk modulus, shear modulus, Young’s modulus, and Poisson’s ratio—of a solvated Arg-Gly-Asp (RGD) peptide model. Since this peptide serves as the RGD-directed integrin recognition site and may participate in cellular adhesion, it is considered a promising small peptide for medicinal applications. This successful approach paves the way for investigating larger and more complex biomolecules.

Biochemistry & Molecular Biology↗

DOE OSTI2023

AGC-2 Specimen Post Irradiation Data Package Report

This report documents results of the post-irradiation examination material property testing of the creep, control, and piggyback specimens from the irradiation creep capsule Advanced Graphite Creep (AGC)-2 are reported. This is the second of a series of six irradiation test trains planned as part of the AGC experiment to fully characterize the neutron irradiation effects and radiation creep behavior of current nuclear graphite grades. The AGC-2 capsule was irradiated in the Idaho National Laboratory Advanced Test Reactor at a nominal temperature of 600°C and to a peak dose of 5 dpa (displacements per atom). One half of the creep specimens were subjected to mechanical stresses (an applied stress of either 13.8, 17.2, or 20.7 MPa) to induce irradiation creep. All post-irradiation testing and measurement results are reported with the exception of the irradiation mechanical strength testing, which is the last destructive testing stage of the irradiation testing program. Material property tests were conducted on specimens from 15 nuclear graphite grades using a similar loading configuration as the first AGC capsule (AGC-1) to provide easy comparison between the two capsules. However, AGC-2 contained an increased number of specimens (i.e., 487 total specimens irradiated) and replaced specimens of the minor grade 2020 with the newer grade 2114. The data reported include specimen dimensions for both stressed and unstressed specimens to establish the irradiation creep rates, mass and volume data necessary to derive density, elastic constants (Young’s modulus, shear modulus, and Poisson’s ratio) from ultrasonic time of flight velocity measurements, Young’s modulus from the fundamental frequency of vibration, electrical resistivity, and thermal diffusivity and thermal expansion data from 100–500°C. No data outliers were determined after all measurements were completed. A brief statistical analysis was performed on the irradiated data and a limited comparison between pre- and post-irradiation properties is presented. A more complete evaluation of trends in the material property changes, as well as irradiation-induced creep due to irradiation, temperature, and applied load on specimens will be discussed in later AGC-2 post-irradiation examination analysis reports.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

AGC-4 Specimen Post-Irradiation Examination Data Interim Report

This interim report documents preliminary results of the post-irradiation examination material property testing from the fourth advanced graphite creep (AGC), AGC 4, capsule specimens. This is the fourth of a series of six irradiation test trains planned as part of the AGC experiment to fully characterize the neutron irradiation effects and radiation creep behavior of current nuclear graphite grades to moderate dose levels (=7 dpa). The AGC 4 capsule was irradiated in the Idaho National Laboratory Advanced Test Reactor at a nominal temperature of 800°C and to a peak dose of 8 dpa. Half of the AGC-4 specimens were subjected to compressive stresses to induce irradiation creep. Post-irradiation testing and measurement results are reported with the exception of thermal testing, which is still in progress, and irradiation mechanical strength testing. Additionally, some specimens initially deemed too hot to be examined in the ART Graphite laboratory may still be measured. The data reported includes specimen dimensions for both stressed and unstressed specimens to establish the irradiation creep rates, mass and dimensional data necessary to derive density, elastic constants (Young?s modulus, shear modulus, and Poisson?s ratio) from ultrasonic time of flight velocity measurements, Young?s modulus from the fundamental frequency of vibration, and electrical resistivity. A more complete evaluation of trends in the material property changes, as well as irradiation-induced creep due to the irradiation environment and applied load on the specimens, will be discussed later in AGC 4 post-irradiation examination analysis reports.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

AGC-4 Specimen Post-Irradiation Examination Data Interim Report

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Effect of Cr and Al on Elastic Constants of FeCrAl Alloys Investigated by Molecular Dynamics Method

The FeCrAl alloy system is recognized as one of the candidate materials for accident-tolerant fuel (ATF) cladding in the nuclear power industry due to its high oxidation resistance under irradiation and high-temperature environments. The concentrations of Cr and Al have a significant effect on elastic properties of the FeCrAl alloy. In this work, elastic constants C 11 , C 12 , C 44 , bulk modulus and shear modulus of FeCrAl alloy were calculated with molecular dynamics methods. We explored compositions with 1–15 wt.% Cr and 1–5 wt.% Al at temperatures from 0 K to 750 K. The results show that the concentrations of Al and Cr have different effects on the elastic constants. When the concentration of Al was fixed, a decrease in bulk modulus and shear modulus with increasing Cr content was observed, consistent with previous experimental results. The dependence of elastic constants on temperature was also the same as in the experiments. Investigations into elastic properties of defect-containing alloys have shown that vacancies, voids, interstitials and Cr-rich precipitations have different effects on elastic properties of FeCrAl alloys. Investigations of elastic properties of defect-containing alloys have shown that vacancies, void, interstitials and Cr-rich precipitations have different effects on elastic properties of FeCrAl alloys. Therefore, the present results indicate that both the Cr and Al concentrations and radiation defects should be considered to develop and apply the FeCrAl alloy in ATF design.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Plate motion in sheared granular fault system

Plate motion near the fault gouge layer, and the elastic interplay between the gouge layer and the plate under stick-slip conditions, is key to understanding the dynamics of sheared granular fault systems. Here, a two-dimensional implementation of the combined finite-discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), is used to explicitly simulate a sheared granular gouge fault system. In this work we focus on investigating the influence of normal load, driving shear velocity and plate stiffness on the velocities and displacements in the direction parallel to the shear direction (x-direction) measured at locations on the upper and lower plates just adjacent to the gouge. The simulations show that during slip phases the magnitudes of the measured velocities on the upper and lower plates are proportional to the normal load and may be inversely proportional to the square root of the plate's shear modulus. Whereas, the driving shear velocity does not show distinct influence on the measured velocities. Additionally, large slip velocities are generally associated with large macroscopic friction coefficient drops. For the models subjected to smaller normal loads, larger shear velocities and with stiffer shear plates, the same magnitude of slip velocity could cause a larger drop of macroscopic friction coefficient. During stick phases, the velocities of the upper and lower plates are respectively slightly greater and slightly smaller than half of the driving shear velocity and are both in the same direction of shear. The shear strain rate of the gouge is calculated from this velocity difference between the upper and lower plate during stick phases and thus the gouge effective shear modulus can be calculated. The results show that the gouge effective shear modulus increases proportionally with normal load, while the influence of shear velocity and plate stiffness on gouge effective shear modulus is minor. The simulations address the dynamics of a laboratory-scale fault gouge system and may aid in revealing the complexities of earthquake frictional dynamics.

58 GEOSCIENCES↗

DOE OSTI2020

Softening of breathing elastic mode and trigonal elastic mode in the disordered pyrochlore magnet NaCaCo 2 F 7

Cobalt pyrochlore fluoride NaCaCo 2 F 7 is a disordered frustrated magnet that is composed of Co 2+ ions with an effective spin-12 magnetic moment and that exhibits spin freezing below $T_f ~$ 2.4 K. Here, we perform ultrasound velocity measurements on a single crystal of the cubic NaCaCo 2 F 7 . The temperature dependence of the bulk modulus (the breathing elastic mode) exhibits Curie-type softening upon cooling below ~20 K down to $T_f$, which is suppressed by the magnetic field. This Curie-type softening should be a precursor to the enhancement of the strength of exchange disorder via the spin-lattice coupling, which causes the spin freezing. In contrast to the magnetic-field-suppressed Curie-type softening in the bulk modulus, the trigonal shear modulus exhibits softening with a characteristic minimum upon cooling, which is enhanced by the magnetic field at temperatures below ~20 K. This magnetic-field-enhanced elastic anomaly in the trigonal shear modulus suggests a coupling of the lattice to the dynamical spin-cluster state. For NaCaCo 2 F 7 , the observed elastic anomalies reveal an occurrence of magnetic-field-induced crossover from an isostructural lattice instability toward the spin freezing to a trigonal lattice instability arising from the emergent dynamical spin-cluster state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Stepwise Dynamic Calibration of a Hydromechanical Simulation Using Time-Lapse Vertical Seismic Profile

This study aims to develop a methodology for calibrating subsurface stress changes through time-lapse Vertical Seismic Profiling (VSP) integration. The selected study site is the 13-10A injector well within the ongoing CO2-EOR operation of the Farnsworth Field Unit. The Time-lapse VSP dataset carries the combined effects of fluid substitution and mean effective stress changes, thereby providing a dataset amenable for the calibration of production and injection-induced stress changes. The concept is similar to calibrating a reservoir simulation model in that the process honor real field data to set up an inverse problem. The solution optimizes the independent and impactful geomechanical parameters that replicate the observed time-lapse seismic velocity changes. This stress calibration is enabled by 4D geomechanical modeling and the VSP Integration workflow. This calibration benefits from extensive geological, geophysical and geomechanical characterization through 3D seismic data, geophysical well logs, and core assessed as part of the 1D MEM conducted on the 13-10A subject well. These data are used to develop a site-specific rock physics model. The Biot Gassmann workflow combines rock physics and reservoir simulation outputs to determine the fluid substitution contribution to seismic velocity change. Additionally, modeled seismic velocity attributed to mean effective stress are determined from the geomechanical simulation outputs, and the stress-velocity relationship developed from the ultrasonic seismic velocity measurements on the extracted Morrow B core. A penalty function is then formed between the modeled seismic velocities and the observed time-lapse VSP dataset. Four independent and impactful geomechanical parameters have been determined. These are the bulk modulus and shear modulus for zero porosity and the shear and compressional seismic velocity to mean effective stress derivatives. The dataset of numerous coupled hydromechanical- geomechanical simulation realizations is built by combining variations of the four stated geomechanical parameters. A machine learning-assisted workflow comprised of an artificial neural network and a particle swarm optimizer are used to converge on the optimal geomechanical parameters. The successful execution of this workflow has affirmed the suitability of acoustic time-lapse measurements for 4D-VSP geomechanical stress calibration pending measurable stress sensitivities within the anticipated effective stress changes and the availability of suitable and reliable datasets for petroelastic modeling.

02 PETROLEUM↗

DOE OSTI2023

Effects of rejuvenation modes on the microstructures and mechanical properties of metallic glasses

This study systematically investigates the effects of three rejuvenation modes, namely biaxial loading, thermal pressure loading, and cryogenic thermal cycling, on the microstructures and mechanical properties of Cu 64 Zr 36 metallic glasses (MGs) using molecular dynamics (MD) simulations. The simulation results demonstrate that biaxial tension application during loading results in a rejuvenated glassy state with high strength. Specifically, we find that the yield stress, Young's modulus, and shear modulus of the rejuvenated glass obtained by biaxial tensile loading decreased by 2.8%, 2.1%, and 2.7%, respectively, compared to an as-cast MG reference. We also find that MGs achieve the greatest degree of rejuvenation, at the expense of strength, with the use of thermal pressure loading. Interestingly, the application of cryogenic thermal cycling results in a rejuvenated state with a lower strength drop compared to thermal pressure loading. Furthermore, analysis of the potential energy and medium-range-order (MRO) size indicates that the microstructures of MGs achieves a high-energy ordered glass state after thermal pressure loading, but a high-energy disordered glass state after biaxial tensile loading and cryogenic thermal cycling. Furthermore, these findings provide meaningful insights into the effect of rejuvenation modes on the microstructures and mechanical properties of MGs, as well as useful guidelines for designing MGs with high strength and superior plasticity for prospective applications in industry.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Predicting Elastic Properties of Materials from Electronic Charge Density Using 3D Deep Convolutional Neural Networks

Materials representation plays a key role in machine learning-based prediction of materials properties and new materials discovery. Currently both graph and three-dimensional (3D) voxel representation methods are based on the heterogeneous elements of the crystal structures. Here, we propose to use electronic charge density (ECD) as a generic unified 3D descriptor for materials property prediction with the advantage of possessing close relation with the physical and chemical properties of materials. We developed an ECD-based 3D convolutional neural networks (CNNs) for predicting the elastic properties of materials, in which CNNs can learn effective hierarchical features with multiple convolving and pooling operations. Extensive benchmark experiments over 2170 $Fm\bar3m$ face-centered-cubic materials show that our ECD-based CNNs can achieve good performance for elasticity prediction. Especially, our CNN models based on the fusion of elemental Materials-Agnostic Platform for Informatics and Exploration features and ECD descriptors achieved the best fivefold cross-validation performance. More importantly, we showed that our ECD-based CNN models can achieve significantly better extrapolation performance when evaluated over nonredundant data sets, where there are few neighbor-training samples around test samples. As an additional validation, we evaluated the predictive performance of our models on 329 materials of space group $Fm\bar3m$ by comparing to density functional theory calculated values, which shows a better prediction power of our model for bulk modulus than shear modulus. Because of the unified representation power of ECD, it is expected that our ECD-based CNN approach can also be applied to predict other physical and chemical properties of crystalline materials.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Effects of Microstructural Composition, Porosity, and Microcracks on the Elastic Moduli of Nuclear Graphites

The elastic responses of nuclear graphite not only depend on the graphitic content itself but are largely dictated by the microstructural constitution of the material. The types of raw materials combined with the manufacturing processes used to produce the graphite yield the microstructural content that primarily includes graphite filler, graphitized pitch binder, and voids/defects that typically occupy approximately 20% of the volume. Among these microstructural components, porosity and microcracking (considered to be part of voids/defects) heavily influence the overall properties of the material, including the elastic moduli. It is widely accepted that the primary effect of oxidation is to increase porosity, but the related effect on the moduli cannot be explained satisfactorily by simply noting changes to porosity. Models describing the elastic moduli of porous, polycrystalline graphite materials have been developed to interpret experimental determinations of Young's modulus and shear modulus in oxidized graphites. Beyond porosity, the moduli are heavily influenced by microcracks, and their effects can be assessed using physical property models based on orientation distribution coefficients and microcrack-modified, single-crystal moduli to represent the combined effects of crystallite orientation and microcracking. This work demonstrates how most modulus measurements in nuclear graphites can be understood using relatively simple models that describe the effects of porosity and microcracking. We also present directions to be pursued to account for microstructure-related effects that occur as a result of neutron irradiation.

Spicer, James B.↗

DOE OSTI2023

Material Property Determinations for Rock Valley Legacy Samples in Support of Source Physics Experiment

The Source Physics Experiment (SPE) is a long-term NNSA research and development effort designed to improve nonproliferation verification and monitoring capabilities. The overarching goals of the SPE program are to improve understanding of prompt signals and physical signatures that develop from underground chemical explosions and associated modeling capabilities. Our work focuses on a primary factor controlling chemical explosion induced signals and signatures: the material properties of the rocks in which the chemical explosion takes place. This document reports on material property determinations of legacy core USGS Test Well F and outcrop analogs for the subsurface stratigraphy for the third phase of SPE in the Rock Valley (RV) area of the Nevada National Security Site (NNSS). The objective of this work is to establish a baseline set of lithologic descriptions and material properties expected prior to observatory borehole drilling in support of the SPE-RVDC (Rock Valley Direct Comparison) experiment. We determine for each rock type the compressional failure envelope, elastic properties as a function of stress (bulk modulus versus mean stress, shear modulus versus shear stress, Young’s modulus versus axial stress and Poisson’s ratio versus axial stress), indirect tensile strength, and porosity. Geologic characterization, both at the core-scale and microscale, provides context for using the data in modeling efforts and to inform interpretations for the material properties testing.

36 MATERIALS SCIENCE↗

DOE OSTI2022

High-Entropy Borides under Extreme Environment of Pressures and Temperatures

The high-entropy transition metal borides containing a random distribution of five or more constituent metallic elements offer novel opportunities in designing materials that show crystalline phase stability, high strength, and thermal oxidation resistance under extreme conditions. We present a comprehensive theoretical and experimental investigation of prototypical high-entropy boride (HEB) materials such as (Hf, Mo, Nb, Ta, Ti)B2 and (Hf, Mo, Nb, Ta, Zr)B2 under extreme environments of pressures and temperatures. The theoretical tools include modeling elastic properties by special quasi-random structures that predict a bulk modulus of 288 GPa and a shear modulus of 215 GPa at ambient conditions. HEB samples were synthesized under high pressures and high temperatures and studied to 9.5 GPa and 2273 K in a large-volume pressure cell. The thermal equation of state measurement yielded a bulk modulus of 276 GPa, in excellent agreement with theory. The measured compressive yield strength by radial X-ray diffraction technique in a diamond anvil cell was 28 GPa at a pressure of 65 GPa, which is a significant fraction of the shear modulus at high pressures. The high compressive strength and phase stability of this material under high pressures and high temperatures make it an ideal candidate for application as a structural material in nuclear and aerospace fields.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Effects of Salt on Phase Behavior and Rheological Properties of Alginate–Chitosan Polyelectrolyte Complexes

Oppositely charged polyelectrolytes often form polyelectrolyte complexes (PECs) due to the association through electrostatic interactions. Obtaining PECs using natural, biocompatible polyelectrolytes is of interest in the food, pharmaceutical, and biomedical industries. In this work, PECs were prepared from two biopolymers, positively charged chitosan and negatively charged alginate. We investigate the changes in the structure and properties of PECs by adding sodium chloride (salt doping) to the system. The shear modulus of PECs can be tuned from ~10 to 10 4 Pa by changing the salt concentration. The addition of salt led to a decrease in the water content of the complex phase with increasing shear modulus. However, at a very high salt concentration, the shear modulus of the complex phase decreased but did not lead to the liquid coacervate formation, typical of synthetic polyelectrolytes. This difference in phase behavior has likely been attributed to the hydrophobicity of chitosan and long semiflexible alginate and chitosan chains that restrict the conformational changes. Large amplitude oscillatory shear experiments captured nonlinear responses of PECs. The compositions of the PECs, determined as a function of salt concentration, signify the preferential partitioning of salt into the complex phase. Small-angle X-ray scattering of the salt-doped PECs indicates that the Kuhn length and radius of the alginate–chitosan associated structure qualitatively agree with the captured phase behavior and rheological data. This study provides insights into the structure–property as a function of salt concentration of natural polymer-based PECs necessary for developing functional materials from natural polyelectrolytes.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Relation between interfacial shear and friction force in 2D materials

Understanding the interfacial properties between an atomic layer and its substrate is of key interest at both the fundamental and technological levels. From Fermi level pinning to strain engineering and superlubricity, the interaction between a single atomic layer and its substrate governs electronic, mechanical and chemical properties. Here, we measure the hardly accessible interfacial transverse shear modulus of an atomic layer on a substrate. By performing measurements on bulk graphite, and on epitaxial graphene films on SiC with different stacking orders and twisting, as well as in the presence of intercalated hydrogen, we find that the interfacial transverse shear modulus is critically controlled by the stacking order and the atomic layer–substrate interaction. Importantly, we demonstrate that this modulus is a pivotal measurable property to control and predict sliding friction in supported two-dimensional materials. The experiments demonstrate a reciprocal relationship between friction force per unit contact area and interfacial shear modulus. As a result, the same relationship emerges from simulations with simple friction models, where the atomic layer– substrate interaction controls the shear stiffness and therefore the resulting friction dissipation.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Calibration and validation of the foundation for a multiphase strength model for tin

In this work, the Common Model of Multi-phase Strength and Equation of State (CMMP) model was applied to tin. Specifically, calibrations of the strength-specific elements of the CMMP foundation were developed with a combination of experiments and theory, and then the model was validated experimentally. The first element of the foundation is a multi-phase analytic treatment of the melt temperature and the shear modulus for the solid phases. These models were parameterized for each phase based on ab initio calculations using the software VASP (Vienna Ab initio Simulations Package) based on density functional theory. The shear modulus model for the ambient phase was validated with ultrasonic sound speed measurements as a function of pressure and temperature. The second element of the foundation is a viscoplastic strength model for the phase, upon which strength for inaccessible higher-pressure phases can be scaled as necessary. The stress–strain response of tin was measured at strain rates of 10 -3 to 3 x 10 3 s -1 and temperatures ranging from 87 to 373 K. The Preston–Tonks–Wallace (PTW) strength model was fit to that data using Bayesian model calibration. For validation, six forward and two reverse Taylor impact experiments were performed at different velocities to measure large plastic deformation of tin at strain rates up to ⁠10 5 s -1 . The PTW model accurately predicted the deformed shapes of the cylinders, with modest discrepancies attributed to the inability of PTW to capture the effects of twinning and dynamic recrystallization. Some material in the simulations of higher velocity Taylor cylinders reached the melting temperature, thus testing the multiphase model because of the presence of a second phase, the liquid. In simulations using a traditional modeling approach, the abrupt reduction of strength upon melt resulted in poor predictions of the deformed shape and non-physical temperatures. With CMMP, the most deformed material points evolved gradually to a mixed solid–liquid but never a fully liquid state, never fully lost strength, stayed at the melt temperature as the latent heat of fusion was absorbed, and predicted the deformed shape well.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Full-Field Strain Measurement Integrated with Two Dimension Regression Analysis to Evaluate the Bi-Modulus Elastic Properties of Isotropic and Transversely Isotropic Materials

Background: Measuring the physical properties of shale is critical for optimizing engineering activities such as geothermal energy generation and hydraulic fracturing. Shale is a transversely isotropic material. Furthermore, this material can also include micro and macro cracks at different locations and orientations that cause it to behave differently under tensile or compressive loading. Objective: In this work, a combined experimental–numerical approach is proposed to evaluate the bi-modulus elastic properties of isotropic and transversely isotropic materials. Methods: Full-field strain measurements for a circular disk under diametral compression are integrated with a regression analysis technique to evaluate the elastic properties of bi-modulus materials subjected to tensile and compressive loads using two loading configurations on the same specimen. Digital Image Correlation (DIC) is used to measure the full-field strains. Subsequently, in the case of an isotropic material, a linear least-squares approach is utilized to process the experimentally determined strains in conjunction with analytical expressions of the stress fields (in terms of far-field loading) to determine the elastic modulus E, the shear modulus G, and the Poisson’s ratio $v$. In the case of a transversely isotopic material, such as shale, a finite element model is implemented to determine the stress fields (again in terms of far-field loading), which is followed by repeating the previous regression analysis in an iterative process to estimate the elastic parameters. Results: The results show that the proposed technique successfully provides a complete set of elastic properties as a function of both the loading condition and the principal material directions. The technique is validated by measurements on a known isotropic material and then applied to determine the properties of shale. Conclusion: In this work, the proposed approach is successfully used to calculate the bi-modulus elastic response of poly(methyl meth- acrylate) (PMMA) and shale. As expected, PMMA exhibits an isotropic response with no bi-modulus effect, however, shale exhibits both transverse isotropy and a bi-modulus effect. Therefore, this approach holds promise for investigating the elastic properties of materials like rocks and fiber-reinforced composite laminates as functions of the principal material directions and the loading conditions.

42 ENGINEERING↗

DOE OSTI2021

Ab Initio Modeling of fcc Fe-Co-Cr-Ni High Entropy Alloys with Full Composition Range

The Fe-Co-Cr-Ni quaternary system has been studied extensively in the past decade, not only because of the superior properties achieved like high tensile ductility and fracture toughness, but also it is the foundation for the development of FeCrCoNi-based HEAs. However, most of the investigations are mainly focusing on the equiatomic and semi-equiatomic compositions. The physical properties of non-equiatomic compositions are barely explored. In the current work, the ab initio approach is adopted to predict the lattice parameter, structure stabilities, elastic properties, and enthalpy of formations of fcc Fe-Cr-Co-Ni HEA single crystals with the special quasirandom structure (SQS) method. To expand the design of Fe-Co-Cr-Ni HEAs into a full composition range, the current simulations cover all the binaries, ternaries, and the whole quaternary system. The 2-D composition-property contour diagrams are created to show the physical properties in the ternary and quaternary systems. Furthermore, our investigation shows that the fcc alloys' lattice parameter increases with the addition of Fe and Cr and decreases with Ni and Co. In addition to the enthalpy of formation at o K, the shape distortion rate is proposed to be another essential parameter to evaluate crystal stability. The results show the addition of Cr would destabilize the fcc lattice and cause symmetry breaking. Finally, the result of Pugh’s ratio shows that most of the fcc alloys show ductile behavior (P > 1.75), especially for the alloys around Fe 3 Ni, FeCo, and Fe 2 NiCo (P > 2.5). The brittle behavior (P < 1.75) is located around Fe 3 Cr and Fe 2 CoCr. Meanwhile, Cr 2 NiCo and Cr 2 FeNiCo are considered promising compositions to be tried and verified experimentally because of the higher bulk modulus and moderate shear modulus.

36 MATERIALS SCIENCE↗