The mechanical properties of Mg–Al–Ca alloys are significantly affected by their Laves phases,including the Al_(2)Ca phase.Laves phases are generally considered to be brittle and have a detrimental effect on the duc...The mechanical properties of Mg–Al–Ca alloys are significantly affected by their Laves phases,including the Al_(2)Ca phase.Laves phases are generally considered to be brittle and have a detrimental effect on the ductility of Mg.Recently,the Al_(2)Ca phase was shown to undergo plastic deformation in a dilute Mg-Al-Ca alloy to increase the ductility and work hardening of the alloy.In the present study,we investigated the extent to which the deformation of Al_(2)Ca is driven by dislocations in the Mg matrix by simulating the interactions between the basal edge dislocations and Al_(2)Ca particles.In particular,the effects of the interparticle spacing,particle orientation,and particle size were considered.Shearing of small particles and dislocation cross-slips near large particles were observed.Both events contribute to strengthening,and accommodate to plasticity.The shear resistance of the dislocation to bypass the particles increased as the particle size increased.The critical resolved shear stress(CRSS)for activating dislocations and stacking faults was easier to reach for small Al_(2)Ca particles owing to the higher local shear stress,which is consistent with the experimental observations.Overall,this work elucidates the driving force for Al_(2)Ca particles in Mg–Al–Ca alloys to undergo plastic deformation.展开更多
Grain boundary(GB)segregation substantially influences the mechanical properties and performance of magnesium(Mg).Atomic-scale modeling,typically using ab-initio or semi-empirical approaches,has mainly focused on GB s...Grain boundary(GB)segregation substantially influences the mechanical properties and performance of magnesium(Mg).Atomic-scale modeling,typically using ab-initio or semi-empirical approaches,has mainly focused on GB segregation at highly symmetric GBs in Mg alloys,often failing to capture the diversity of local atomic environments and segregation energies,resulting in inaccurate structure-property predictions.This study employs atomistic simulations and machine learning models to systematically investigate the segregation behavior of common solute elements in polycrystalline Mg at both 0 K and finite temperatures.The machine learning models accurately predict segregation thermodynamics by incorporating energetic and structural descriptors.We found that segregation energy and vibrational free energy follow skew-normal distributions,with hydrostatic stress,an indicator of excess free volume,emerging as an important factor influencing segregation tendency.The local atomic environment's flexibility,quantified by flexibility volume,is also crucial in predicting GB segregation.Comparing the grain boundary solute concentrations calculated via the Langmuir-Mc Lean isotherm with experimental data,we identified a pronounced segregation tendency for Nd,highlighting its potential for GB engineering in Mg alloys.This work demonstrates the powerful synergy of atomistic simulations and machine learning,paving the way for designing advanced lightweight Mg alloys with tailored properties.展开更多
Spurious forces are a significant challenge for multi-scale methods,e.g.,the coupled atomistic/discrete dislocation(CADD)method.The assumption of isotropic matter in the continuum domain is a critical factor leading t...Spurious forces are a significant challenge for multi-scale methods,e.g.,the coupled atomistic/discrete dislocation(CADD)method.The assumption of isotropic matter in the continuum domain is a critical factor leading to such forces.This study aims to minimize spurious forces,ensuring that atomic dislocations experience more precise forces from the continuum domain.The authors have already implemented this idea using a simplified and unrealistic slipping system.To create a comprehensive and realistic model,this paper considers all possible slip systems in the face center cubic(FCC)lattice structure,and derives the required relationships for the displacement fields.An anisotropic version of the three-dimensional CADD(CADD3D)method is presented,which generates the anisotropic displacement fields for the partial dislocations in all the twelve slip systems of the FCC lattice structure.These displacement fields are tested for the most probable slip systems of aluminum,nickel,and copper with different anisotropic levels.Implementing these anisotropic displacement fields significantly reduces the spurious forces on the slip systems of FCC materials.This improvement is particularly pronounced at greater distances from the interface and in more anisotropic materials.Furthermore,the anisotropic CADD3D method enhances the spurious stress difference between the slip systems,particularly for materials with higher anisotropy.展开更多
Porous materials are widely used in the field of protection because of their excellent energy absorption characteristics.In this work,a series of polyurethane microscopic models are established and the effect of poros...Porous materials are widely used in the field of protection because of their excellent energy absorption characteristics.In this work,a series of polyurethane microscopic models are established and the effect of porosity on the shock waves is studied with classical molecular dynamics simulations.Firstly,shock Hugoniot relations for different porosities are obtained,which compare well with the experimental data.The pores collapse and form local stress wave,which results in the complex multi-wave structure of the shock wave.The microstructure analysis shows that the local stress increases and the local velocity decreases gradually during the process of pore collapse to complete compaction.Finally,it leads to stress relaxation and velocity homogenization.The shock stress peaks can be fitted with two exponential functions,and the amplitude of attenuation coefficient decreases with the increase of density.Besides,the pore collapse under shock or non-shock are discussed by the entropy increase rate of the system.The energy is dissipated mainly through the multiple interactions of the waves under shock.The energy is dissipated mainly by the friction between atoms under non-shock.展开更多
MicroMagnetic.jl is an open-source Julia package for micromagnetic and atomistic simulations.Using the features of the Julia programming language,MicroMagnetic.jl supports CPU and various GPU platforms,including NVIDI...MicroMagnetic.jl is an open-source Julia package for micromagnetic and atomistic simulations.Using the features of the Julia programming language,MicroMagnetic.jl supports CPU and various GPU platforms,including NVIDIA,AMD,Intel,and Apple GPUs.Moreover,MicroMagnetic.jl supports Monte Carlo simulations for atomistic models and implements the nudged-elastic-band method for energy barrier computations.With built-in support for double and single precision modes and a design allowing easy extensibility to add new features,MicroMagnetic.jl provides a versatile toolset for researchers in micromagnetics and atomistic simulations.展开更多
The tension and compression of face-centered-cubic high-entropy alloy(HEA) nanowires are significantly asymmetric, but the tension–compression asymmetry in nanoscale body-centered-cubic(BCC) HEAs is still unclear. In...The tension and compression of face-centered-cubic high-entropy alloy(HEA) nanowires are significantly asymmetric, but the tension–compression asymmetry in nanoscale body-centered-cubic(BCC) HEAs is still unclear. In this study,the tension–compression asymmetry of the BCC Al Cr Fe Co Ni HEA nanowire is investigated using molecular dynamics simulations. The results show a significant asymmetry in both the yield and flow stresses, with BCC HEA nanowire stronger under compression than under tension. The strength asymmetry originates from the completely different deformation mechanisms in tension and compression. In compression, atomic amorphization dominates plastic deformation and contributes to the strengthening, while in tension, deformation twinning prevails and weakens the HEA nanowire.The tension–compression asymmetry exhibits a clear trend of increasing with the increasing nanowire cross-sectional edge length and decreasing temperature. In particular, the compressive strengths along the [001] and [111] crystallographic orientations are stronger than the tensile counterparts, while the [110] crystallographic orientation shows the exactly opposite trend. The dependences of tension–compression asymmetry on the cross-sectional edge length, crystallographic orientation,and temperature are explained in terms of the deformation behavior of HEA nanowire as well as its variations caused by the change in these influential factors. These findings may deepen our understanding of the tension–compression asymmetry of the BCC HEA nanowires.展开更多
Molecular dynamics （MD） simulations of monocrystalline copper （100） surface during nanomachining process were performed based on a new 3D simulation model. The material removal mechanism and system temperature dis...Molecular dynamics （MD） simulations of monocrystalline copper （100） surface during nanomachining process were performed based on a new 3D simulation model. The material removal mechanism and system temperature distribution were discussed. The simulation results indicate that the system temperature distribution presents a roughly concentric shape, a steep temperature gradient is observed in diamond cutting tool, and the highest temperature is located in chip. Centrosymmetry parameter method was used to monitor defect structures. Dislocations and vacancies are the two principal types of defect structures. Residual defect structures impose a major change on the workpiece physical properties and machined surface quality. The defect structures in workpiece are temperature dependent. As the temperature increases, the dislocations are mainly mediated from the workpiece surface, while the others are dissociated into point defects. The relatively high cutting speed used in nanomachining results in less defect structures, beneficial to obtain highly machined surface quality.展开更多
Deformation twinning, i.e., twin nucleation and twin growth （or twin boundary migration, TBM） activated by impinged basal slip at a symmetrical tilt grain boundary in HCP Mg, was examined with molecular dynamics （M...Deformation twinning, i.e., twin nucleation and twin growth （or twin boundary migration, TBM） activated by impinged basal slip at a symmetrical tilt grain boundary in HCP Mg, was examined with molecular dynamics （MD） simulations. The results show that the {1^-1^-21}-type twinning acts as the most preferential mode of twinning. Once such twins are formed, they are almost ready to grow. The TBM of such twins is led by pure atomic shuffling events. A secondary mode of twinning can also occur in our simulations. The {112^-2} twinning is observed at 10 K as the secondary twin. This secondary mode of twinning shows different energy barriers for nucleation as well as for growth compared with the {1^-1^-21}-type twining. In particular, TBMs in this case is triggered intrinsically by pyramidal slip at its twin boundary.展开更多
Coupled atomistic/dislocation/continuum simulation of interfacial fracture is performed in this paper.The model consists of a nanoscopic core made by atomistic assembly and a surrounding elastic continuum with discret...Coupled atomistic/dislocation/continuum simulation of interfacial fracture is performed in this paper.The model consists of a nanoscopic core made by atomistic assembly and a surrounding elastic continuum with discrete dislocations. Atomistic dislocations nucleate from the crack tip and move to the continuum layer where they glide according to the dislocation dynamics curve.An atoms/continuum overlapping belt is devised to facilitate the transition between the two scales.The continuum constraint on the atomic assembly is imposed through the mechanics at- mosphere along the overlapping belt.Transmissions of mechanics parameters such as displacements,stresses,masses and momenta across the belt are realized.The present model allows us to explore interfacial fracture processes under different mode mixity.The effect of atomistic zigzag interface on the fracture process is revealed:it hinders dislocation emission from the crack tip,especially under high mode mixity.展开更多
The phenomenon of interfacial fracture, as manifested by atom- istic cleavage, debonding and dislocation emission, provides a challenge for combined atomistic-continuum analysis. As a precursor for fully coupled atomi...The phenomenon of interfacial fracture, as manifested by atom- istic cleavage, debonding and dislocation emission, provides a challenge for combined atomistic-continuum analysis. As a precursor for fully coupled atomistic-continuum simulation of interfacial fracture, we focus here on the atomistic behavior within a nanoscopic core surrounding the crack tip. The inter-atomic potential under Em- bedded Atom Method is recapitulated to form an essential framework of atomistic simulation. The calculations are performed for a side-cracked disc configuration un- der a remote K field loading. It is revealed that a critical loading rate defines the brittle-to-ductile transition of homogeneous materials. We further observe that the near tip mode mixity dictates the nanoscopic profile near an interfacial crack tip. A zigzag interface structure is simulated which plays a significant role in the dislocation emission from an interfacial crack tip, as will be explored in the second part of this investigation.展开更多
The phase field crystal method and Continuum Modeling are applied to study the cooperative dislocation motion of the grain boundary(GB)migration,the manner of the nucleation of the grain and of the grain growth in two...The phase field crystal method and Continuum Modeling are applied to study the cooperative dislocation motion of the grain boundary(GB)migration,the manner of the nucleation of the grain and of the grain growth in two dimensions(2 D)under the deviatoric deformation at high temperature.Three types of the nucleation modes of new finding are observed by the phase field crystal simulation:The first mode of the nucleation is generated by the GB splitting into two sub-GBs;the second mode is of the reaction of the sub-GB dislocations,such as,the generation and annihilation of a pair of partial Frank sessile dislocation in 2 D.The process can be considered as the nucleation of dynamic recrystallization;the third mode is caused by two oncoming rows of the dislocations of these sub-GBs,crossing and passing each other to form new gap which is the nucleation place of the new deformed grain.The research is shown that due to the nucleation of different modes the mechanism of the grain growth by means of the sub-GB migration is different,and therefore,the grain growth rates are also different.Under the deviatoric deformation of the applied biaxial strain,the grain growth is faster than that of the grain growth without external applied stress.It is observed that the cooperative dislocation motion of the GB migration under the deviatoric deformation accompanies with local plastic flow and the state of the stress of the system changes sharply.When the system is in the process of recrystallized grain growth,the system energy is in an unstable state due to the release of the strain energy to cause that the reverse movement of the plastic flow occurs.The area growth of the deformed grain is approximately proportional to the strain square and also to the time square.The rule of the time square of the deformed grain growth can also be deduced by establishing the continuum dynamic equation of the biaxial strain-driven migration of the GB.The copper metal is taken as an example of the calculation,and the obtained result is a good agreement with that of the experiment.展开更多
How to correctly extract Cauchy stress from the atomistic simulations is a crucial issue in studying the mechanical behaviours of atomic systems, but is still in controversy. In this paper, three typical atomistic sim...How to correctly extract Cauchy stress from the atomistic simulations is a crucial issue in studying the mechanical behaviours of atomic systems, but is still in controversy. In this paper, three typical atomistic simulation examples are used to validate various existing stress definitions. It is found that the classical virial stress fails in predicting the stresses in these examples, because the velocity depends on the choice of the local average volume or the reference frame velocity and other factors. In contrast, the Lagrangian cross-section stress and Lagrangian virial stress are validated by these examples, and the instantaneous Lagrangian atomic stress definition is also proposed for dynamical problems.展开更多
It is well known that precipitation hardening in magnesium(Mg)alloys is far less effective than in aluminum alloys.Thus,it is important to understand the surface and interfacial structure and energetics between precip...It is well known that precipitation hardening in magnesium(Mg)alloys is far less effective than in aluminum alloys.Thus,it is important to understand the surface and interfacial structure and energetics between precipitates and matrix.In upscale modeling of magnesium alloys,these energy data are of great significance.In this work,we calculated the surface and interfacial energies of Mg_(17)Al_(12)-Mg system by carefully selecting the surface or interface termination,using atomistic simulations.The results show that,the higher fraction of Mg atoms on the surface,the lower the surface energy of Mg_(17)Al_(12).The interfacial energy of Mg/Mg_(17)Al_(12)was calculated in which the Burgers orientation relationship(OR)was satisfied.It was found that the(011)P|(0002)Mg interface has the lowest interfacial energy(248 mJ/m 2).Because the Burgers OR breaks when{10¯12}twin occurs,which reorients the matrix,the interfacial energy for Mg_(17)Al_(12)and a{10¯12}twin was also calculated.The results show that after twinning,the lowest interfacial energy increases by 244 mJ/m^(2),and the interface becomes highly incoherent due to the change in orientation relationship between Mg_(17)Al_(12)and the matrix.展开更多
The deformation behavior in magnesium single crystal under c-axis tension is investigated in a temperature range between 250 K and 570 K by molecular dynamics simulations. At a low temperature, twinning and shear band...The deformation behavior in magnesium single crystal under c-axis tension is investigated in a temperature range between 250 K and 570 K by molecular dynamics simulations. At a low temperature, twinning and shear bands are found to be the main deformation mechanisms. In particular, the {102} tension twins with the reorientation angle of about 90 °are observed in the simulations. The mechanisms of {102} twinning are illustrated by the simulated motion of atoms. Moreover, grain nucleation and growth are found to be accompanied with the {102} twinning. At temperatures above 450 K, the twin frequency decreases with increasing temperature. The {102} extension twin almost disappears at the temperature of 570 K. The non-basal slip plays an important role on the tensile deformation in magnesium single crystal at high temperatures.展开更多
Two types of Lennard-Jones potential are widely used in modeling adhesive contacts. However, the relationships between the parameters of the two types of Lennard-Jones potential are not well defined. This paper employ...Two types of Lennard-Jones potential are widely used in modeling adhesive contacts. However, the relationships between the parameters of the two types of Lennard-Jones potential are not well defined. This paper employs a self- consistent method to derive the Lennard-Jones surface force law from the interatomic Lennard-Jones potential with emphasis on the relationships between the parameters. The ei^ect of using correct parameters in the adhesion models is demonstrated in single sphere-flat contact via continuum models and an atomistic model. Furthermore, the adhesion hysteresis behaviour is investigated, and the S-shaped force-distance relation is revealed by the atomistic model. It shows that the adhesion hysteresis loop is generated by the jump-to-contact and jump-off-contact, which are illustrated by the S-shaped force-distance curve.展开更多
Free transverse vibration of monolayer graphene, boron nitride (BN), and silicon carbide (SiC) sheets is investigated by using molecular dynamics finite element method. Eigenfrequencies and eigenmodes of these three s...Free transverse vibration of monolayer graphene, boron nitride (BN), and silicon carbide (SiC) sheets is investigated by using molecular dynamics finite element method. Eigenfrequencies and eigenmodes of these three sheets in rectangular shape are studied with different aspect ratios with respect to various boundary conditions. It is found that aspect ratios and boundary conditions affect in a similar way on natural frequencies of graphene, BN, and SiC sheets. Natural frequencies in all modes decrease with an increase of the sheet’s size. Graphene exhibits the highest natural frequencies, and SiC sheet possesses the lowest ones. Missing atoms have minor effects on natural frequencies in this study.展开更多
We investigate the mechanical and microstructural changes of the densified silica glass under uniaxial loading-unloading via atomistic simulations with a modified BKS potential. The stress–strain relationship is foun...We investigate the mechanical and microstructural changes of the densified silica glass under uniaxial loading-unloading via atomistic simulations with a modified BKS potential. The stress–strain relationship is found to include three respective stages: elastic, plastic and hardening regions. The bulk modulus increases with the initial densification and will undergo a rapid increase after complete densification. The yield pressure varies from 5 to 12 GPa for different densified samples. In addition, the Si–O–Si bond angle reduces during elastic deformation under compression, and 5-fold Si will increase linearly in the plastic deformation. In the hardening region, the peak splitting and the new peak are both found on the Si–Si and O–O pair radial distribution functions, where the 6-fold Si is increased. Instead, the lateral displacement of the atoms always varies linearly with strain, without evident periodic characteristic. As is expected, the samples are permanently densified after release from the plastic region, and the maximum density of recovered samples is about 2.64 g/cm^3, which contains 15 % 5-fold Si, and the Si–O–Si bond angle is less than the ordinary silica glass. All these findings are of great significance for understanding the deformation process of densified silica glass.展开更多
A large number of Embedded Atom Method (EAM) potentials have been developed for the Ni/Al system. These potentials are compared to a common data base. It is found that there is significant difference in quality in the...A large number of Embedded Atom Method (EAM) potentials have been developed for the Ni/Al system. These potentials are compared to a common data base. It is found that there is significant difference in quality in these potentials. One of the potentials has also been extended to represent the properties of hydrogen in Ni/Al intermetallics. This potential describes the solution and migration behavior of hydrogen in both Ni and Al.A number of calculations using the Ni/Al/H potential have been performed. It is found that hydrogen strongly prefers sites in Mi3AI that are surrounded by 6 Ni atoms. Calculations of trapping of hydrogen to a number of grain boundaries in Ni3Al have been performed as a function of hydrogen chemical potential at room temperature. The failure of these bicrystals under tensile stress has been examined and will be compared to the failure of pure Mi3AI boundaries.In order to investigate the potential embrittlement of γ/γ'alloys, trapping of hydrogen to a spherical Hi3Al precipitate in Ni as a function of chemical potential at room temperature has been calculated. It appears that the boundary is not a strong trap for hydrogen, hence embrittlement in these alloys is not primarily due to interactions of hydrogen with the γ/γ'interface. (Supported by the U. S. DOE through Contract DE-AC04-94AL85000.)展开更多
A series of material parameters are derived from atomistic simulations and implemented into a phase field（PF） model to simulate void evolution in body-centered cubic（bcc） iron subjected to different irradiation do...A series of material parameters are derived from atomistic simulations and implemented into a phase field（PF） model to simulate void evolution in body-centered cubic（bcc） iron subjected to different irradiation doses at different temperatures.The simulation results show good agreement with experimental observations — the porosity as a function of temperature varies in a bell-shaped manner and the void density monotonically decreases with increasing temperatures; both porosity and void density increase with increasing irradiation dose at the same temperature. Analysis reveals that the evolution of void number and size is determined by the interplay among the production, diffusion and recombination of vacancy and interstitial.展开更多
As one of the fundamental outcomes of dislocation self-interaction,dislocation dipoles have an important influence on the plastic deformation of materials,especially on fatigue and creep.In this work,superdislocation ...As one of the fundamental outcomes of dislocation self-interaction,dislocation dipoles have an important influence on the plastic deformation of materials,especially on fatigue and creep.In this work,superdislocation dipoles inγ-TiAl andα_(2)-Ti_(3)Al were systematically investigated by atomistic simulations,with a variety of dipole heights,orientations and annealing tempe ratures.The results indicate that non-screw super-dipoles transform into locally stable dipolar or reconstructed cores at low temperature,while into isolated or interconnected point defect clusters and stacking fault tetrahedra at high temperature via short-range diffu sion.Non-screw super-dipoles inγ-TiAl andα_(2)-Ti_(3)Al exhibit similar features as fcc and hcp metals,respectively.Generally,over long-term annealing where diffusion is significant,60°superdipoles inγ-TiAl are stable,whereas the stability of super-dipoles inα2-Ti3 Al increases with dipole height and orientation angle.The influence on mechanical properties can be well evaluated by integrating these results into mesoscale or constitutive models.展开更多
基金funded by the National Natural Science Foundation of China(nos.51631006 and 51825101)。
文摘The mechanical properties of Mg–Al–Ca alloys are significantly affected by their Laves phases,including the Al_(2)Ca phase.Laves phases are generally considered to be brittle and have a detrimental effect on the ductility of Mg.Recently,the Al_(2)Ca phase was shown to undergo plastic deformation in a dilute Mg-Al-Ca alloy to increase the ductility and work hardening of the alloy.In the present study,we investigated the extent to which the deformation of Al_(2)Ca is driven by dislocations in the Mg matrix by simulating the interactions between the basal edge dislocations and Al_(2)Ca particles.In particular,the effects of the interparticle spacing,particle orientation,and particle size were considered.Shearing of small particles and dislocation cross-slips near large particles were observed.Both events contribute to strengthening,and accommodate to plasticity.The shear resistance of the dislocation to bypass the particles increased as the particle size increased.The critical resolved shear stress(CRSS)for activating dislocations and stacking faults was easier to reach for small Al_(2)Ca particles owing to the higher local shear stress,which is consistent with the experimental observations.Overall,this work elucidates the driving force for Al_(2)Ca particles in Mg–Al–Ca alloys to undergo plastic deformation.
基金Z.X.and T.A.S.acknowledge the financial support by the German Research Foundation(DFG)(Grant Nr.505716422)T.A.S.are grateful for the financial support from the DFG(Grant Nr.AL1343/7-1,AL1343/8-1 and Yi 103/3-1)+4 种基金Z.X.,S.K.K.and U.K.acknowledge financial support by the DFG through the projects A05,A07 and C02 of the SFB1394 StructuralChemical Atomic Complexity-From Defect Phase Diagrams to Material Properties,project ID 409476157Additionally,Z.X.and S.K.K.are grateful for funding from the European Research Council(ERC)under the European Union’s Horizon 2020 research and innovation program(grant agreement No.852096 FunBlocks)J.G.acknowledges funding from the French National Research Agency(ANR),Grant ANR-21-CE08-0001(ATOUUM)and ANR-22-CE92-0058-01(SILA)The authors gratefully acknowledge the computing time provided to them at the NHR Center NHR4CES at RWTH Aachen University(project number p0020431 and p0020267)。
文摘Grain boundary(GB)segregation substantially influences the mechanical properties and performance of magnesium(Mg).Atomic-scale modeling,typically using ab-initio or semi-empirical approaches,has mainly focused on GB segregation at highly symmetric GBs in Mg alloys,often failing to capture the diversity of local atomic environments and segregation energies,resulting in inaccurate structure-property predictions.This study employs atomistic simulations and machine learning models to systematically investigate the segregation behavior of common solute elements in polycrystalline Mg at both 0 K and finite temperatures.The machine learning models accurately predict segregation thermodynamics by incorporating energetic and structural descriptors.We found that segregation energy and vibrational free energy follow skew-normal distributions,with hydrostatic stress,an indicator of excess free volume,emerging as an important factor influencing segregation tendency.The local atomic environment's flexibility,quantified by flexibility volume,is also crucial in predicting GB segregation.Comparing the grain boundary solute concentrations calculated via the Langmuir-Mc Lean isotherm with experimental data,we identified a pronounced segregation tendency for Nd,highlighting its potential for GB engineering in Mg alloys.This work demonstrates the powerful synergy of atomistic simulations and machine learning,paving the way for designing advanced lightweight Mg alloys with tailored properties.
文摘Spurious forces are a significant challenge for multi-scale methods,e.g.,the coupled atomistic/discrete dislocation(CADD)method.The assumption of isotropic matter in the continuum domain is a critical factor leading to such forces.This study aims to minimize spurious forces,ensuring that atomic dislocations experience more precise forces from the continuum domain.The authors have already implemented this idea using a simplified and unrealistic slipping system.To create a comprehensive and realistic model,this paper considers all possible slip systems in the face center cubic(FCC)lattice structure,and derives the required relationships for the displacement fields.An anisotropic version of the three-dimensional CADD(CADD3D)method is presented,which generates the anisotropic displacement fields for the partial dislocations in all the twelve slip systems of the FCC lattice structure.These displacement fields are tested for the most probable slip systems of aluminum,nickel,and copper with different anisotropic levels.Implementing these anisotropic displacement fields significantly reduces the spurious forces on the slip systems of FCC materials.This improvement is particularly pronounced at greater distances from the interface and in more anisotropic materials.Furthermore,the anisotropic CADD3D method enhances the spurious stress difference between the slip systems,particularly for materials with higher anisotropy.
基金financial support from National Natural Science Foundation of China(Grant No.12172325)。
文摘Porous materials are widely used in the field of protection because of their excellent energy absorption characteristics.In this work,a series of polyurethane microscopic models are established and the effect of porosity on the shock waves is studied with classical molecular dynamics simulations.Firstly,shock Hugoniot relations for different porosities are obtained,which compare well with the experimental data.The pores collapse and form local stress wave,which results in the complex multi-wave structure of the shock wave.The microstructure analysis shows that the local stress increases and the local velocity decreases gradually during the process of pore collapse to complete compaction.Finally,it leads to stress relaxation and velocity homogenization.The shock stress peaks can be fitted with two exponential functions,and the amplitude of attenuation coefficient decreases with the increase of density.Besides,the pore collapse under shock or non-shock are discussed by the entropy increase rate of the system.The energy is dissipated mainly through the multiple interactions of the waves under shock.The energy is dissipated mainly by the friction between atoms under non-shock.
基金supported by the National Key R&D Program of China(Grant No.2022YFA1403603)the Strategic Priority Research Program of Chinese Academy of Sciences(Grant No.XDB33030100)+2 种基金the National Natural Science Fund for Distinguished Young Scholar(Grant No.52325105)the National Natural Science Foundation of China(Grant Nos.12374098,11974021,and 12241406)the CAS Project for Young Scientists in Basic Research(Grant No.YSBR-084).
文摘MicroMagnetic.jl is an open-source Julia package for micromagnetic and atomistic simulations.Using the features of the Julia programming language,MicroMagnetic.jl supports CPU and various GPU platforms,including NVIDIA,AMD,Intel,and Apple GPUs.Moreover,MicroMagnetic.jl supports Monte Carlo simulations for atomistic models and implements the nudged-elastic-band method for energy barrier computations.With built-in support for double and single precision modes and a design allowing easy extensibility to add new features,MicroMagnetic.jl provides a versatile toolset for researchers in micromagnetics and atomistic simulations.
基金Project supported by the National Natural Science Foundation of China (Grant No.12272118)the National Key Research and Development Program of China (Grant No.2022YFE03030003)。
文摘The tension and compression of face-centered-cubic high-entropy alloy(HEA) nanowires are significantly asymmetric, but the tension–compression asymmetry in nanoscale body-centered-cubic(BCC) HEAs is still unclear. In this study,the tension–compression asymmetry of the BCC Al Cr Fe Co Ni HEA nanowire is investigated using molecular dynamics simulations. The results show a significant asymmetry in both the yield and flow stresses, with BCC HEA nanowire stronger under compression than under tension. The strength asymmetry originates from the completely different deformation mechanisms in tension and compression. In compression, atomic amorphization dominates plastic deformation and contributes to the strengthening, while in tension, deformation twinning prevails and weakens the HEA nanowire.The tension–compression asymmetry exhibits a clear trend of increasing with the increasing nanowire cross-sectional edge length and decreasing temperature. In particular, the compressive strengths along the [001] and [111] crystallographic orientations are stronger than the tensile counterparts, while the [110] crystallographic orientation shows the exactly opposite trend. The dependences of tension–compression asymmetry on the cross-sectional edge length, crystallographic orientation,and temperature are explained in terms of the deformation behavior of HEA nanowire as well as its variations caused by the change in these influential factors. These findings may deepen our understanding of the tension–compression asymmetry of the BCC HEA nanowires.
基金Project (50925521) supported by the National Natural Science Fund for Distinguished Young Scholars of China
文摘Molecular dynamics （MD） simulations of monocrystalline copper （100） surface during nanomachining process were performed based on a new 3D simulation model. The material removal mechanism and system temperature distribution were discussed. The simulation results indicate that the system temperature distribution presents a roughly concentric shape, a steep temperature gradient is observed in diamond cutting tool, and the highest temperature is located in chip. Centrosymmetry parameter method was used to monitor defect structures. Dislocations and vacancies are the two principal types of defect structures. Residual defect structures impose a major change on the workpiece physical properties and machined surface quality. The defect structures in workpiece are temperature dependent. As the temperature increases, the dislocations are mainly mediated from the workpiece surface, while the others are dissociated into point defects. The relatively high cutting speed used in nanomachining results in less defect structures, beneficial to obtain highly machined surface quality.
基金Project(2012CB932202)supported by the National Basic Research Program of ChinaProjects(50890174,50971088)supported by the National Natural Science Foundation of China
文摘Deformation twinning, i.e., twin nucleation and twin growth （or twin boundary migration, TBM） activated by impinged basal slip at a symmetrical tilt grain boundary in HCP Mg, was examined with molecular dynamics （MD） simulations. The results show that the {1^-1^-21}-type twinning acts as the most preferential mode of twinning. Once such twins are formed, they are almost ready to grow. The TBM of such twins is led by pure atomic shuffling events. A secondary mode of twinning can also occur in our simulations. The {112^-2} twinning is observed at 10 K as the secondary twin. This secondary mode of twinning shows different energy barriers for nucleation as well as for growth compared with the {1^-1^-21}-type twining. In particular, TBMs in this case is triggered intrinsically by pyramidal slip at its twin boundary.
基金The project supported by the National Natural Science Foundation of China
文摘Coupled atomistic/dislocation/continuum simulation of interfacial fracture is performed in this paper.The model consists of a nanoscopic core made by atomistic assembly and a surrounding elastic continuum with discrete dislocations. Atomistic dislocations nucleate from the crack tip and move to the continuum layer where they glide according to the dislocation dynamics curve.An atoms/continuum overlapping belt is devised to facilitate the transition between the two scales.The continuum constraint on the atomic assembly is imposed through the mechanics at- mosphere along the overlapping belt.Transmissions of mechanics parameters such as displacements,stresses,masses and momenta across the belt are realized.The present model allows us to explore interfacial fracture processes under different mode mixity.The effect of atomistic zigzag interface on the fracture process is revealed:it hinders dislocation emission from the crack tip,especially under high mode mixity.
基金The project supported by the National Natural Science Foundation of China
文摘The phenomenon of interfacial fracture, as manifested by atom- istic cleavage, debonding and dislocation emission, provides a challenge for combined atomistic-continuum analysis. As a precursor for fully coupled atomistic-continuum simulation of interfacial fracture, we focus here on the atomistic behavior within a nanoscopic core surrounding the crack tip. The inter-atomic potential under Em- bedded Atom Method is recapitulated to form an essential framework of atomistic simulation. The calculations are performed for a side-cracked disc configuration un- der a remote K field loading. It is revealed that a critical loading rate defines the brittle-to-ductile transition of homogeneous materials. We further observe that the near tip mode mixity dictates the nanoscopic profile near an interfacial crack tip. A zigzag interface structure is simulated which plays a significant role in the dislocation emission from an interfacial crack tip, as will be explored in the second part of this investigation.
基金supported by National Nature Science Foundation of China(Nos.51161003 and 51561031)Nature Science Foundation of Guangxi Province(No.2018GXNSFAA138150)。
文摘The phase field crystal method and Continuum Modeling are applied to study the cooperative dislocation motion of the grain boundary(GB)migration,the manner of the nucleation of the grain and of the grain growth in two dimensions(2 D)under the deviatoric deformation at high temperature.Three types of the nucleation modes of new finding are observed by the phase field crystal simulation:The first mode of the nucleation is generated by the GB splitting into two sub-GBs;the second mode is of the reaction of the sub-GB dislocations,such as,the generation and annihilation of a pair of partial Frank sessile dislocation in 2 D.The process can be considered as the nucleation of dynamic recrystallization;the third mode is caused by two oncoming rows of the dislocations of these sub-GBs,crossing and passing each other to form new gap which is the nucleation place of the new deformed grain.The research is shown that due to the nucleation of different modes the mechanism of the grain growth by means of the sub-GB migration is different,and therefore,the grain growth rates are also different.Under the deviatoric deformation of the applied biaxial strain,the grain growth is faster than that of the grain growth without external applied stress.It is observed that the cooperative dislocation motion of the GB migration under the deviatoric deformation accompanies with local plastic flow and the state of the stress of the system changes sharply.When the system is in the process of recrystallized grain growth,the system energy is in an unstable state due to the release of the strain energy to cause that the reverse movement of the plastic flow occurs.The area growth of the deformed grain is approximately proportional to the strain square and also to the time square.The rule of the time square of the deformed grain growth can also be deduced by establishing the continuum dynamic equation of the biaxial strain-driven migration of the GB.The copper metal is taken as an example of the calculation,and the obtained result is a good agreement with that of the experiment.
基金supported by the National Natural Science Foundation of China (Grant Nos10702034,10732050 and 90816006)the National Basic Research Program of China (973 Program 2007CB936803,2010CB832701)
文摘How to correctly extract Cauchy stress from the atomistic simulations is a crucial issue in studying the mechanical behaviours of atomic systems, but is still in controversy. In this paper, three typical atomistic simulation examples are used to validate various existing stress definitions. It is found that the classical virial stress fails in predicting the stresses in these examples, because the velocity depends on the choice of the local average volume or the reference frame velocity and other factors. In contrast, the Lagrangian cross-section stress and Lagrangian virial stress are validated by these examples, and the instantaneous Lagrangian atomic stress definition is also proposed for dynamical problems.
基金Bin Li gratefully thank support from the U.S.National Science Foundation(CMMI-1635088).
文摘It is well known that precipitation hardening in magnesium(Mg)alloys is far less effective than in aluminum alloys.Thus,it is important to understand the surface and interfacial structure and energetics between precipitates and matrix.In upscale modeling of magnesium alloys,these energy data are of great significance.In this work,we calculated the surface and interfacial energies of Mg_(17)Al_(12)-Mg system by carefully selecting the surface or interface termination,using atomistic simulations.The results show that,the higher fraction of Mg atoms on the surface,the lower the surface energy of Mg_(17)Al_(12).The interfacial energy of Mg/Mg_(17)Al_(12)was calculated in which the Burgers orientation relationship(OR)was satisfied.It was found that the(011)P|(0002)Mg interface has the lowest interfacial energy(248 mJ/m 2).Because the Burgers OR breaks when{10¯12}twin occurs,which reorients the matrix,the interfacial energy for Mg_(17)Al_(12)and a{10¯12}twin was also calculated.The results show that after twinning,the lowest interfacial energy increases by 244 mJ/m^(2),and the interface becomes highly incoherent due to the change in orientation relationship between Mg_(17)Al_(12)and the matrix.
基金supported by National Natural Science Foundation of China (GrantNos.11072026 and 10632020)the Fundamental Research Funds for the Central Universities, and finalized during a sabbatical leave of D.S. at the Graduate Institute of Ferrous Technology (G.I.F.T.) of POSTECHPohang, Korea as part of an International Outgoing Fellowship (Marie Curie Actions) of the 7th Programme of the European Commission
文摘The deformation behavior in magnesium single crystal under c-axis tension is investigated in a temperature range between 250 K and 570 K by molecular dynamics simulations. At a low temperature, twinning and shear bands are found to be the main deformation mechanisms. In particular, the {102} tension twins with the reorientation angle of about 90 °are observed in the simulations. The mechanisms of {102} twinning are illustrated by the simulated motion of atoms. Moreover, grain nucleation and growth are found to be accompanied with the {102} twinning. At temperatures above 450 K, the twin frequency decreases with increasing temperature. The {102} extension twin almost disappears at the temperature of 570 K. The non-basal slip plays an important role on the tensile deformation in magnesium single crystal at high temperatures.
基金supported by the National Natural Science Foundation of China (Grant No. 10476019)the Fundamental Research Funds for the Central Universities of China (Grant No. JY10000904018)
文摘Two types of Lennard-Jones potential are widely used in modeling adhesive contacts. However, the relationships between the parameters of the two types of Lennard-Jones potential are not well defined. This paper employs a self- consistent method to derive the Lennard-Jones surface force law from the interatomic Lennard-Jones potential with emphasis on the relationships between the parameters. The ei^ect of using correct parameters in the adhesion models is demonstrated in single sphere-flat contact via continuum models and an atomistic model. Furthermore, the adhesion hysteresis behaviour is investigated, and the S-shaped force-distance relation is revealed by the atomistic model. It shows that the adhesion hysteresis loop is generated by the jump-to-contact and jump-off-contact, which are illustrated by the S-shaped force-distance curve.
文摘Free transverse vibration of monolayer graphene, boron nitride (BN), and silicon carbide (SiC) sheets is investigated by using molecular dynamics finite element method. Eigenfrequencies and eigenmodes of these three sheets in rectangular shape are studied with different aspect ratios with respect to various boundary conditions. It is found that aspect ratios and boundary conditions affect in a similar way on natural frequencies of graphene, BN, and SiC sheets. Natural frequencies in all modes decrease with an increase of the sheet’s size. Graphene exhibits the highest natural frequencies, and SiC sheet possesses the lowest ones. Missing atoms have minor effects on natural frequencies in this study.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.51727807 and 11875318)Beijing Institute of Technology Research Fund Program for Young ScholarsYue Qi Young Scholar Project in CUMTB。
文摘We investigate the mechanical and microstructural changes of the densified silica glass under uniaxial loading-unloading via atomistic simulations with a modified BKS potential. The stress–strain relationship is found to include three respective stages: elastic, plastic and hardening regions. The bulk modulus increases with the initial densification and will undergo a rapid increase after complete densification. The yield pressure varies from 5 to 12 GPa for different densified samples. In addition, the Si–O–Si bond angle reduces during elastic deformation under compression, and 5-fold Si will increase linearly in the plastic deformation. In the hardening region, the peak splitting and the new peak are both found on the Si–Si and O–O pair radial distribution functions, where the 6-fold Si is increased. Instead, the lateral displacement of the atoms always varies linearly with strain, without evident periodic characteristic. As is expected, the samples are permanently densified after release from the plastic region, and the maximum density of recovered samples is about 2.64 g/cm^3, which contains 15 % 5-fold Si, and the Si–O–Si bond angle is less than the ordinary silica glass. All these findings are of great significance for understanding the deformation process of densified silica glass.
文摘A large number of Embedded Atom Method (EAM) potentials have been developed for the Ni/Al system. These potentials are compared to a common data base. It is found that there is significant difference in quality in these potentials. One of the potentials has also been extended to represent the properties of hydrogen in Ni/Al intermetallics. This potential describes the solution and migration behavior of hydrogen in both Ni and Al.A number of calculations using the Ni/Al/H potential have been performed. It is found that hydrogen strongly prefers sites in Mi3AI that are surrounded by 6 Ni atoms. Calculations of trapping of hydrogen to a number of grain boundaries in Ni3Al have been performed as a function of hydrogen chemical potential at room temperature. The failure of these bicrystals under tensile stress has been examined and will be compared to the failure of pure Mi3AI boundaries.In order to investigate the potential embrittlement of γ/γ'alloys, trapping of hydrogen to a spherical Hi3Al precipitate in Ni as a function of chemical potential at room temperature has been calculated. It appears that the boundary is not a strong trap for hydrogen, hence embrittlement in these alloys is not primarily due to interactions of hydrogen with the γ/γ'interface. (Supported by the U. S. DOE through Contract DE-AC04-94AL85000.)
基金Project supported by the National Magnetic Confinement Fusion Energy Research Project of China(Grant No.2015GB118001)the Fundamental Research Funds for the Central Universities,China(Grant No.DUT16RC(3)052)+1 种基金the National Basic Research Program of China(Grant No.2012CB619402)the NETL Project(Grant No.DE-FE0027776)
文摘A series of material parameters are derived from atomistic simulations and implemented into a phase field（PF） model to simulate void evolution in body-centered cubic（bcc） iron subjected to different irradiation doses at different temperatures.The simulation results show good agreement with experimental observations — the porosity as a function of temperature varies in a bell-shaped manner and the void density monotonically decreases with increasing temperatures; both porosity and void density increase with increasing irradiation dose at the same temperature. Analysis reveals that the evolution of void number and size is determined by the interplay among the production, diffusion and recombination of vacancy and interstitial.
基金the National Key Research and Development Program of China(No.2016YFB0701304 and 2017YFB0306201)the Natural Science Foundation of China(Nos.51671195 and 91960202)+4 种基金the Frontier and Key Projects of the Chinese Academy of Sciences(Nos.QYZDJ-SSW-JSC031-01 and XXH13506-304)the Natural Science Foundation of Liaoning(No.20180510032)the Aeronautical Science Foundation of China(No.20160292002)the Strategic Priority Research Program of Chinese Academy of Sciences(No.XDC01000000)The Project is sponsored by the“Liaoning BaiQianWan”Talents Program。
文摘As one of the fundamental outcomes of dislocation self-interaction,dislocation dipoles have an important influence on the plastic deformation of materials,especially on fatigue and creep.In this work,superdislocation dipoles inγ-TiAl andα_(2)-Ti_(3)Al were systematically investigated by atomistic simulations,with a variety of dipole heights,orientations and annealing tempe ratures.The results indicate that non-screw super-dipoles transform into locally stable dipolar or reconstructed cores at low temperature,while into isolated or interconnected point defect clusters and stacking fault tetrahedra at high temperature via short-range diffu sion.Non-screw super-dipoles inγ-TiAl andα_(2)-Ti_(3)Al exhibit similar features as fcc and hcp metals,respectively.Generally,over long-term annealing where diffusion is significant,60°superdipoles inγ-TiAl are stable,whereas the stability of super-dipoles inα2-Ti3 Al increases with dipole height and orientation angle.The influence on mechanical properties can be well evaluated by integrating these results into mesoscale or constitutive models.
