Shaped charge liner(SCL)has been extensively applied in oil recovery and defense industries.Achieving superior penetration capability through optimizing SCL structures presents a substantial challenge due to intricate...Shaped charge liner(SCL)has been extensively applied in oil recovery and defense industries.Achieving superior penetration capability through optimizing SCL structures presents a substantial challenge due to intricate rate-dependent processes involving detonation-driven liner collapse,high-speed jet stretching,and penetration.This study introduces an innovative optimization strategy for SCL structures that employs jet penetration efficiency as the primary objective function.The strategy combines experimentally validated finite element method with machine learning(FEM-ML).We propose a novel jet penetration efficiency index derived from enhanced cutoff velocity and shape characteristics of the jet via machine learning.This index effectively evaluates the jet penetration performance.Furthermore,a multi-model fusion based on a machine learning optimization method,called XGBOOST-MFO,is put forward to optimize SCL structure over a large input space.The strategy's feasibility is demonstrated through the optimization of copper SCL implemented via the FEM-ML strategy.Finally,this strategy is extended to optimize the structure of the recently emerging CrMnFeCoNi high-entropy alloy conical liners and hemispherical copper liners.Therefore,the strategy can provide helpful guidance for the engineering design of SCL.展开更多
Metals have been mankind’s most essential materials for thousands of years.In recent years,however,innovation-driven development of major national security strategy and core areas of the national economy is highly im...Metals have been mankind’s most essential materials for thousands of years.In recent years,however,innovation-driven development of major national security strategy and core areas of the national economy is highly impeded by a shortage of advanced higher-strength-toughness metals.One of the main reasons is that metals inherently exhibit the inverted-relationship of strength-toughness.The emergence of two types of disordered metals:amorphous alloys and high entropy alloys,provides a fully-fresh strategy for strength-toughening by tailoring the topological and/or chemical disorder.In this paper,we first briefly review the history of strength-toughening of metals,and summarize the development route-map.We then introduce amorphous alloys and high entropy alloys,as well as some case studies in tailoring disorder to successfully achieve coexisting high strength and high ductility/toughness.Relevant challenges that await further research are summarized in concluding remarks.展开更多
Physical aging is an inherent property of glassy matter, but understanding its microscopic mechanism remains a challenge particularly at the particle level. In this work, we use a confocal microscope to in-situ trace ...Physical aging is an inherent property of glassy matter, but understanding its microscopic mechanism remains a challenge particularly at the particle level. In this work, we use a confocal microscope to in-situ trace the particle trajectories in a 3D colloidal glass for 73000 s, aiming at resolving the aging dynamics. By calculating the mean square displacement of particle motions, we find that the glass aging with time can be divided into three stages: β relaxation, α relaxation and free diffusion. The system's mean square displacement at each aging state is quantitatively resolved into three contributions of particle dynamics modes: vibration within the nearest-neighbor cages, hopping between cages and cooperative rearrangement. We further calculate the particle's free volume and find that the β-to-α transition is accompanied by the temporary increase of the system-averaged free volume due to pronounced hops of particles. Nevertheless, the temporal autocorrelation of the free volume spatial distribution still obeys a monotonically stretched exponential decay with an exponent of 0.76, which is related to the sub-diffusion dynamics of cooperative rearrangements and hops mixed in α relaxation. According to the resolved vibrational displacements,we calculate the vibrational density of states of this 3D glass, and the characteristic boson peak is reproduced at low frequencies.Our findings shed insight into the particle-level aging dynamics of a real glass under purely thermal activation.展开更多
The competition between physical aging and structural rejuvenation determines the physical and mechanical properties of glassy materials.Thus,the rejuvenation-aging boundary must be identified quantitatively.In this w...The competition between physical aging and structural rejuvenation determines the physical and mechanical properties of glassy materials.Thus,the rejuvenation-aging boundary must be identified quantitatively.In this work,we unravel a stress boundary to distinguish rejuvenation from aging via the thermo-mechanical creep of a typical Zr-based metallic glass.The crept glasses were rejuvenated into high-enthalpy disordered states when the applied stress exceeded a threshold that was numerically close to the steady-state flow stress;otherwise,the glasses were aged.A theoretical model for glass creep was adopted to demystify the observed stress threshold of rejuvenation.The model revealed that the thermo-mechanical creep beyond the threshold stress could activate sufficient shear transformations to create a net free volume,thus leading to structural rejuvenation.Furthermore,we derived the analytical expressions for the threshold and flow stresses.Both stresses can act as the rejuvenation-aging boundary,which is well supported by experimental creep data.The present work procures a deeper understanding of the rejuvenation mechanism of glasses and provides useful implications for abstaining from glass aging.展开更多
Glasses with rejuvenated structures usually exhibit improved room-temperature plasticity,which facilitates their applications.However,glass rejuvenation requires external energy injection to“shake up”the frozen-in d...Glasses with rejuvenated structures usually exhibit improved room-temperature plasticity,which facilitates their applications.However,glass rejuvenation requires external energy injection to“shake up”the frozen-in disordered structure.In this work,we give the answer to how much the required energy is.According to the constitutive model of amorphous plasticity,we find that the applied stress higher than the steady-state flow value can effectively induce the structural disordering in terms of the generation of free volume.Therefore,the effective energy density(EED)of structural rejuvenation is defined as the integral of this effective stress on the corresponding strain.By tailoring the applied strain,strain rate,temperature and initial free volume,different degrees of structural rejuvenation are achieved,which show a generally linear correlation with the defined EED.This work deepens the understanding of glass rejuvenation from an energy perspective.展开更多
In this study, mechanical tests were conducted oil a face-centered cubic FeCoNiCrMn high-entropy alloy, both in tension and compression, in a wide range of strain rates (10^-4-10^4 s^-1) to systematically investigat...In this study, mechanical tests were conducted oil a face-centered cubic FeCoNiCrMn high-entropy alloy, both in tension and compression, in a wide range of strain rates (10^-4-10^4 s^-1) to systematically investigate its dynamic response and underlying deformation mechanism. Materials with different grain sizes were tested to understand the effect of grain size, thus grain boundary volume, on the mechanical prop-erties. Microstructures of various samples both before and after deformation were examined using elec-tron backscatter diffraction and transmission electron microscopy. The dislocation structure as well as deformation-induced twins were analyzed and correlated with the measured mechanical properties. Plastic stability during tension of the current high-entropy alloy (HEA), in particular, at dynamic strain rates, was discussed in lights of strain-rate sensitivity and work hardening rate. It was found that, under dynamic conditions, the strength and uniform ductility increased simultaneously as a result of the mas-sive formation of deformation twins. Specifically, an ultimate tensile strength of 734 MPa and uniform elongation of-63% are obtained at 2.3×10^3 s^-1, indicating that the alloy has great potential for energy absorption upon impact loading.展开更多
The martensitic transformation(MT)lays the foundation for microstructure and performance tailoring of many engineering materials,especially steels,which are with>1.8 billion tons produced per year the most importan...The martensitic transformation(MT)lays the foundation for microstructure and performance tailoring of many engineering materials,especially steels,which are with>1.8 billion tons produced per year the most important material class.The atomic-scale migration path is a long-term challenge for MT dur-ing quenching in high-carbon(nitrogen)steels.Here,we provide direct evidence of(11^(-)2)body-centred tetragonal(BCT)twinned martensite in carbon steels by transmission electron microscopy(TEM)investi-gation,and the increase in tetragonality with the C content matches X-ray diffraction(XRD)results.The specific{11^(-)2}_(BCT)twin planes which are related to the elongated c axis provide essential structural details to revisit the migration path of the atoms in MT.Therefore,the face-centred cubic(FCC)to BCT twin to body-centred cubic(BCC)twin transition pathway and its underlying mechanisms are revealed through direct experimental observation and atomistic simulations.Our findings shed new light on the nature of the martensitic transition,thus providing new opportunities for the nanostructural control of metals and alloys.展开更多
The shear banding instability occurs as the homogenous deformation in metallic glasses (MGs) develops to a critical point, at which the discontinuity in deformation rate is incipient across nano-scale shear bands. W...The shear banding instability occurs as the homogenous deformation in metallic glasses (MGs) develops to a critical point, at which the discontinuity in deformation rate is incipient across nano-scale shear bands. When and where the shear instability takes place is an important issue for understanding the shear band origin. However, such condition and direction of shear localization concerning the unique properties of MGs is still lacking for general stress state. In this paper, a new constitutive is introduced for MGs accounting for the pressure sensitivity, dilatancy and structural evolution; the shear banding is regarded as the appearance of instability in the constitutive description of inelastic deformation. Tying the bifurcation theory to the new constitutive, the general condition of deformation localization is derived. The shear band orientation corresponding to the easiest direction of shear instability is then obtained in dependence on pressure sensitivity, dilatancy and Poisson's ratio for MGs. The range of the predicted shear band angles is consistent with the experimental observations.展开更多
In this work,we comprehensively investigate the response of amorphous and crystalline Zr-based alloys under nanosecond pulse laser ablation.The in situ multiphysics processes and ablation morphologies of the two alloy...In this work,we comprehensively investigate the response of amorphous and crystalline Zr-based alloys under nanosecond pulse laser ablation.The in situ multiphysics processes and ablation morphologies of the two alloy targets are explored and compared.The results indicate that the dynamics of laser-induced plasma and shock waves obey the idea blast wave theory and are insensitive to the topological structures of targets.Both targets experience significant superheating and culminate in explosive boiling.This ablation process leads to the formation of a hierarchical structure in the resultant ablation crater:microdents covered by widespread nanovoids.The amorphous target shows shallower microdents and smaller nanovoids than their crystalline counterparts because the former has a smaller heat-affected zone and experiences a higher degree of superheating.The hierarchical structure can adjust the surface wettability of targets from initial hydrophilic to hydrophobic,showing an increase of the contact angle approximately 119% for amorphous alloy compared with the crystal approximately 64%.This work demonstrates that amorphous alloys have a better performance against nanosecond pulse laser ablation and provides a feasible and one-step method of wettability modification for either amorphous or crystalline alloys.展开更多
Metallic glasses(MGs)constitute an emerging class of advanced structural materials due to their excellent mechanical properties.However,brittle failure at room temperature and the resultant complicated fracture behavi...Metallic glasses(MGs)constitute an emerging class of advanced structural materials due to their excellent mechanical properties.However,brittle failure at room temperature and the resultant complicated fracture behavior greatly limit their wide engineering applications.Over the past decades,the deformation and fracture in ductile or brittle mode referring to material compositions,load conditions,sample size,etc.,have been widely studied,and significant progress has been made in understanding the failure behavior of MGs.Micromechanisms of fracture have been revealed involving shear banding,cavitation and the nature of the crack tip field.The ductile-to-brittle transition and inherent governing parameters have been found.To well describe and predict the failure behavior of MGs,failure criteria for ductile and brittle MGs have been established empirically or based on atomic interactions.In this paper,we provide a detailed review of the above advances and identify outstanding issues in the failure of MGs that need to be further clarified.展开更多
Shear banding in amorphous metals originates from the activation and percolation of flow units.To uncover the self-assembly dynamics of flow units in metallic glasses,a rectangular sample with two flow units embedded ...Shear banding in amorphous metals originates from the activation and percolation of flow units.To uncover the self-assembly dynamics of flow units in metallic glasses,a rectangular sample with two flow units embedded in the matrix undergoing simple shearing was analyzed using finite element simulations.The vortex evolution behavior,including activation,growth,and collapse during the self-assembly of flow units,was revealed.It was found that the formation of a mature vortex indicates the onset of yielding,and the collapse of the vortex represents the percolation of flow units or shear localization.The effects of initial free volume distribution and the distance between flow units on vortex behavior were also studied.Increasing the initial free volume concentration within flow units or the matrix leads to a gentler vortex evolution process and better homogeneous plasticity.The shape of vortex tends to be"flatter"with the increase in flow units'spacing,and the optimal spacing was found to maximize the strength of the material.展开更多
基金supported by the NSFC Basic Science Center Program for"Multi-scale Problems in Nonlinear Mechanics" (Grant No.11988102)the NSFC (Grant Nos.U2141204,12172367)+2 种基金the Key Research Program of the Chinese Academy of Sciences (Grant No.ZDRW-CN-2021-2-3)the National Key Research and Development Program of China (Grant No.2022YFC3320504-02)the opening project of State Key Laboratory of Explosion Science and Technology (Grant No.KFJJ21-01 and No.KFJJ18-14 M)。
文摘Shaped charge liner(SCL)has been extensively applied in oil recovery and defense industries.Achieving superior penetration capability through optimizing SCL structures presents a substantial challenge due to intricate rate-dependent processes involving detonation-driven liner collapse,high-speed jet stretching,and penetration.This study introduces an innovative optimization strategy for SCL structures that employs jet penetration efficiency as the primary objective function.The strategy combines experimentally validated finite element method with machine learning(FEM-ML).We propose a novel jet penetration efficiency index derived from enhanced cutoff velocity and shape characteristics of the jet via machine learning.This index effectively evaluates the jet penetration performance.Furthermore,a multi-model fusion based on a machine learning optimization method,called XGBOOST-MFO,is put forward to optimize SCL structure over a large input space.The strategy's feasibility is demonstrated through the optimization of copper SCL implemented via the FEM-ML strategy.Finally,this strategy is extended to optimize the structure of the recently emerging CrMnFeCoNi high-entropy alloy conical liners and hemispherical copper liners.Therefore,the strategy can provide helpful guidance for the engineering design of SCL.
基金This work was supported by the National Outstanding Youth Science Fund Project of National Natural Science Foundation of China(NSFC)(No.12125206)the NSFC Basic Science Center for“Multiscale Problems in Nonlinear Mechanics”(No.11988102)the NSFC(Nos.11972345 and 11790292).
文摘Metals have been mankind’s most essential materials for thousands of years.In recent years,however,innovation-driven development of major national security strategy and core areas of the national economy is highly impeded by a shortage of advanced higher-strength-toughness metals.One of the main reasons is that metals inherently exhibit the inverted-relationship of strength-toughness.The emergence of two types of disordered metals:amorphous alloys and high entropy alloys,provides a fully-fresh strategy for strength-toughening by tailoring the topological and/or chemical disorder.In this paper,we first briefly review the history of strength-toughening of metals,and summarize the development route-map.We then introduce amorphous alloys and high entropy alloys,as well as some case studies in tailoring disorder to successfully achieve coexisting high strength and high ductility/toughness.Relevant challenges that await further research are summarized in concluding remarks.
基金supported by the National Outstanding Youth Science Fund Project (Grant No. 12125206)Basic Science Center for “Multiscale Problems in Nonlinear Mechanics”(Grant No. 11988102)+1 种基金General Project of National Natural Science Foundation of China (Grant No. 11972345)CAS Project for Young Scientists in Basic Research (Grant No. YSBR-096)。
文摘Physical aging is an inherent property of glassy matter, but understanding its microscopic mechanism remains a challenge particularly at the particle level. In this work, we use a confocal microscope to in-situ trace the particle trajectories in a 3D colloidal glass for 73000 s, aiming at resolving the aging dynamics. By calculating the mean square displacement of particle motions, we find that the glass aging with time can be divided into three stages: β relaxation, α relaxation and free diffusion. The system's mean square displacement at each aging state is quantitatively resolved into three contributions of particle dynamics modes: vibration within the nearest-neighbor cages, hopping between cages and cooperative rearrangement. We further calculate the particle's free volume and find that the β-to-α transition is accompanied by the temporary increase of the system-averaged free volume due to pronounced hops of particles. Nevertheless, the temporal autocorrelation of the free volume spatial distribution still obeys a monotonically stretched exponential decay with an exponent of 0.76, which is related to the sub-diffusion dynamics of cooperative rearrangements and hops mixed in α relaxation. According to the resolved vibrational displacements,we calculate the vibrational density of states of this 3D glass, and the characteristic boson peak is reproduced at low frequencies.Our findings shed insight into the particle-level aging dynamics of a real glass under purely thermal activation.
基金This work was supported by the National Outstanding Youth Science Fund Project(Grant No.12125206)the National Natural Science Foundation of China(Grant Nos.11972345,and 11790292)the NSFC Basic Science Center for“Multiscale Problems in Nonlinear Mechanics”(Grant No.11988102).
文摘The competition between physical aging and structural rejuvenation determines the physical and mechanical properties of glassy materials.Thus,the rejuvenation-aging boundary must be identified quantitatively.In this work,we unravel a stress boundary to distinguish rejuvenation from aging via the thermo-mechanical creep of a typical Zr-based metallic glass.The crept glasses were rejuvenated into high-enthalpy disordered states when the applied stress exceeded a threshold that was numerically close to the steady-state flow stress;otherwise,the glasses were aged.A theoretical model for glass creep was adopted to demystify the observed stress threshold of rejuvenation.The model revealed that the thermo-mechanical creep beyond the threshold stress could activate sufficient shear transformations to create a net free volume,thus leading to structural rejuvenation.Furthermore,we derived the analytical expressions for the threshold and flow stresses.Both stresses can act as the rejuvenation-aging boundary,which is well supported by experimental creep data.The present work procures a deeper understanding of the rejuvenation mechanism of glasses and provides useful implications for abstaining from glass aging.
基金supported by the General Project(No.11972345)the National Outstanding Youth Science Fund Project(No.12125206)+1 种基金the Major Project(No.11790292)of National Natural Science Foundation of China(NSFC)the NSFC Basic Science Center for“Multiscale Problems in Nonlinear Mechanics”(No.11988102).
文摘Glasses with rejuvenated structures usually exhibit improved room-temperature plasticity,which facilitates their applications.However,glass rejuvenation requires external energy injection to“shake up”the frozen-in disordered structure.In this work,we give the answer to how much the required energy is.According to the constitutive model of amorphous plasticity,we find that the applied stress higher than the steady-state flow value can effectively induce the structural disordering in terms of the generation of free volume.Therefore,the effective energy density(EED)of structural rejuvenation is defined as the integral of this effective stress on the corresponding strain.By tailoring the applied strain,strain rate,temperature and initial free volume,different degrees of structural rejuvenation are achieved,which show a generally linear correlation with the defined EED.This work deepens the understanding of glass rejuvenation from an energy perspective.
基金supported by the National Natural Science Foundation of China(51671018,51531001,51422101,51371003,and 51671021)111 Project(B07003)+5 种基金International S&T Cooperation Program of China(2015DFG52600)Program for Changjiang Scholars and Innovative Research Team in University of China(IRT_14R05)the Projects of SKL-AMM-USTB(2016Z-04,2016-09,2016Z-16)the financial support from the Top-Notch Young Talents Programthe Fundamental Research Funds for the Central Universitiesthe financial support by US-NSF under contract DMR-1408722
文摘In this study, mechanical tests were conducted oil a face-centered cubic FeCoNiCrMn high-entropy alloy, both in tension and compression, in a wide range of strain rates (10^-4-10^4 s^-1) to systematically investigate its dynamic response and underlying deformation mechanism. Materials with different grain sizes were tested to understand the effect of grain size, thus grain boundary volume, on the mechanical prop-erties. Microstructures of various samples both before and after deformation were examined using elec-tron backscatter diffraction and transmission electron microscopy. The dislocation structure as well as deformation-induced twins were analyzed and correlated with the measured mechanical properties. Plastic stability during tension of the current high-entropy alloy (HEA), in particular, at dynamic strain rates, was discussed in lights of strain-rate sensitivity and work hardening rate. It was found that, under dynamic conditions, the strength and uniform ductility increased simultaneously as a result of the mas-sive formation of deformation twins. Specifically, an ultimate tensile strength of 734 MPa and uniform elongation of-63% are obtained at 2.3×10^3 s^-1, indicating that the alloy has great potential for energy absorption upon impact loading.
基金supported by the NSFC(Nos.51931005,51901235,11790292)the NSFC Basic Science Center Program for“Multi-scale Problems in Nonlinear Mechanics”(No.11988102)+1 种基金the Key Research Program of Frontier Sciences(No.QYZDJSSW-JSC011)the Strategic Priority Research Program of the Chinese Academy of Sciences(No.XDB22040302),and the CityU grant 9360161.
文摘The martensitic transformation(MT)lays the foundation for microstructure and performance tailoring of many engineering materials,especially steels,which are with>1.8 billion tons produced per year the most important material class.The atomic-scale migration path is a long-term challenge for MT dur-ing quenching in high-carbon(nitrogen)steels.Here,we provide direct evidence of(11^(-)2)body-centred tetragonal(BCT)twinned martensite in carbon steels by transmission electron microscopy(TEM)investi-gation,and the increase in tetragonality with the C content matches X-ray diffraction(XRD)results.The specific{11^(-)2}_(BCT)twin planes which are related to the elongated c axis provide essential structural details to revisit the migration path of the atoms in MT.Therefore,the face-centred cubic(FCC)to BCT twin to body-centred cubic(BCC)twin transition pathway and its underlying mechanisms are revealed through direct experimental observation and atomistic simulations.Our findings shed new light on the nature of the martensitic transition,thus providing new opportunities for the nanostructural control of metals and alloys.
基金Financial support is from the National Key Basic Research Program of China (Grant No. 2012CB937500)the National Natural Science Foundation of China (Grants Nos. 11202221 and 11132011)the CAS/SAFEA International Partnership Program for Creative Research Teams
文摘The shear banding instability occurs as the homogenous deformation in metallic glasses (MGs) develops to a critical point, at which the discontinuity in deformation rate is incipient across nano-scale shear bands. When and where the shear instability takes place is an important issue for understanding the shear band origin. However, such condition and direction of shear localization concerning the unique properties of MGs is still lacking for general stress state. In this paper, a new constitutive is introduced for MGs accounting for the pressure sensitivity, dilatancy and structural evolution; the shear banding is regarded as the appearance of instability in the constitutive description of inelastic deformation. Tying the bifurcation theory to the new constitutive, the general condition of deformation localization is derived. The shear band orientation corresponding to the easiest direction of shear instability is then obtained in dependence on pressure sensitivity, dilatancy and Poisson's ratio for MGs. The range of the predicted shear band angles is consistent with the experimental observations.
基金supported by the National Outstanding Youth Science Fund Project(Grant No.12125206)of the National Natural Science Foundation of China(NSFC)the NSFC Basic Science Center for"Multiscale Problems in Nonlinear Mechanics"(Grant No.11988102)the NSFC(Grant Nos.11972345 and 11790292).
文摘In this work,we comprehensively investigate the response of amorphous and crystalline Zr-based alloys under nanosecond pulse laser ablation.The in situ multiphysics processes and ablation morphologies of the two alloy targets are explored and compared.The results indicate that the dynamics of laser-induced plasma and shock waves obey the idea blast wave theory and are insensitive to the topological structures of targets.Both targets experience significant superheating and culminate in explosive boiling.This ablation process leads to the formation of a hierarchical structure in the resultant ablation crater:microdents covered by widespread nanovoids.The amorphous target shows shallower microdents and smaller nanovoids than their crystalline counterparts because the former has a smaller heat-affected zone and experiences a higher degree of superheating.The hierarchical structure can adjust the surface wettability of targets from initial hydrophilic to hydrophobic,showing an increase of the contact angle approximately 119% for amorphous alloy compared with the crystal approximately 64%.This work demonstrates that amorphous alloys have a better performance against nanosecond pulse laser ablation and provides a feasible and one-step method of wettability modification for either amorphous or crystalline alloys.
基金supported by the National Natural Science Foundation Basic Science Center Program for"Multiscale Problems in Nonlinear Mechanics"(Grant No.11988102)the National Natural Science Foundation of China(Grant Nos.11972346 and 11790292)+3 种基金the National Key Research and Development Program of China(Grant No.2017YFB0702003)the Strategic Priority Research Program(Grant Nos.XDB22040302 and XDB22040303)the Key Research Program of Frontier Sciences(Grant No.QYZDJSSW-JSC011)the Science Challenge Project(Grant No.TZ2018001).
文摘Metallic glasses(MGs)constitute an emerging class of advanced structural materials due to their excellent mechanical properties.However,brittle failure at room temperature and the resultant complicated fracture behavior greatly limit their wide engineering applications.Over the past decades,the deformation and fracture in ductile or brittle mode referring to material compositions,load conditions,sample size,etc.,have been widely studied,and significant progress has been made in understanding the failure behavior of MGs.Micromechanisms of fracture have been revealed involving shear banding,cavitation and the nature of the crack tip field.The ductile-to-brittle transition and inherent governing parameters have been found.To well describe and predict the failure behavior of MGs,failure criteria for ductile and brittle MGs have been established empirically or based on atomic interactions.In this paper,we provide a detailed review of the above advances and identify outstanding issues in the failure of MGs that need to be further clarified.
基金supported by the NSFC(Nos.11972346 and 11790292)the NSFCBasic ScienceCenter Program for“Multi-scale Problems in Nonlinear Mechanics”(No.11988102).
文摘Shear banding in amorphous metals originates from the activation and percolation of flow units.To uncover the self-assembly dynamics of flow units in metallic glasses,a rectangular sample with two flow units embedded in the matrix undergoing simple shearing was analyzed using finite element simulations.The vortex evolution behavior,including activation,growth,and collapse during the self-assembly of flow units,was revealed.It was found that the formation of a mature vortex indicates the onset of yielding,and the collapse of the vortex represents the percolation of flow units or shear localization.The effects of initial free volume distribution and the distance between flow units on vortex behavior were also studied.Increasing the initial free volume concentration within flow units or the matrix leads to a gentler vortex evolution process and better homogeneous plasticity.The shape of vortex tends to be"flatter"with the increase in flow units'spacing,and the optimal spacing was found to maximize the strength of the material.