Underhand cut-and-fill mining has been widely used in underground mining operations,especially when the rock mass or orebody is of poor quality or prone to rockburst due to high stress.In such cases,mining workers sho...Underhand cut-and-fill mining has been widely used in underground mining operations,especially when the rock mass or orebody is of poor quality or prone to rockburst due to high stress.In such cases,mining workers should carry out all production activities under the cemented backfill roof or sill mat instead of a highly fractured and unstable rock roof or a strong rock roof with a high potential of rockburst.Therefore,the stability and required strength of the sill mat are critical issues for mining engineers.In 1991,Mitchell considered that sill mat could fail by caving,sliding,rotation,and flexure.Mitchell also proposed an analytical solution to determine the minimum required strength of the sill mat for each type of failure based on two stiff or immobile rock walls.However,recent publications using numerical modeling and field measurements indicate that the compressive stresses in the sill mat induced by rock wall closure due to a stope excavation beneath the sill mat can be significant.It is thus highly necessary to investigate the required strength of the sill mat by considering rock wall closure.In this study,the crushing failure of sill mat due to rock wall closure generated by underground excavation and a new failure mode called"crushing and caving”is revealed by numerical modeling.An analytical solution corresponding to each failure mode is then developed to estimate the minimum required cohesion(cmin)of the sill mat.A criterion is also proposed to determine if the sill mat fails by crushing or crushing-and-caving failure.The proposed analytical solution does not involve any correction coefficients.The validity of the proposed analytical solution is demonstrated by numerical modeling.The proposed analytical solution can thus be employed to predict the cmin of sill mat subjected to wall closure generated by underlying stope excavation.展开更多
In the underhand cut-and-fill mining method,a sill mat(i.e.an artificial horizontal pillar)constructed by cemented backfill is essential to prevent mine workers from being directly exposed under problematic rock roofs...In the underhand cut-and-fill mining method,a sill mat(i.e.an artificial horizontal pillar)constructed by cemented backfill is essential to prevent mine workers from being directly exposed under problematic rock roofs.A critical issue is to determine the minimum required strength of the sill mat to ensure a safe and cost-effective design.Until now,Mitchell’s analytical solution is the only available option,considering two stiff and immobile rock walls.Unavoidable rock wall closure associated with stope excavation below the sill mat was neglected.This,along with other undefined parameters,explains why Mitchell’s solution is rarely used in sill mat design.A new analytical solution for determining the minimum required strength of the sill mat accounting for wall closure is necessary.In this study,a closed-form analytical solution for estimating rock wall closure generated by stope excavation below a sill mat is developed by using Salamon’s and Flamant’s models.The proposed analytical solution does not contain any coefficients of correction or calibration.Despite several assumptions(or somewhat of oversimplifications)necessary to render a simple analytical solution possible,good agreements are obtained between the rock wall closures predicted by applying the proposed analytical solution and those obtained numerically with FLAC3D for many cases with arbitrarily chosen geometrical and material parameters.The proposed analytical solution is therefore validated and can be used to evaluate the rock wall closure generated by stope excavation below a sill mat.展开更多
A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equ...A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equation is integrated along the wall-normal direction to link the tangential component of the effective body force for the IB method to the wall shear stress predicted by the wall model;(ii) a set of Lagrangian points near the wall are introduced to compute the normal component of the effective body force for the IB method by reconstructing the normal component of the velocity. This novel method will be a classical direct-forcing IB method if the grid is fine enough to resolve the flow near the wall. The method is used to simulate the flows around the DARPA SUBOFF model. The results obtained are well comparable to the measured experimental data and wall-resolved LES results.展开更多
Detailed investigations on the microstructure and the mechanical properties of the wing membrane of the dragonfly are carried out. It is found that in the direction of the thickness the membrane was divided into three...Detailed investigations on the microstructure and the mechanical properties of the wing membrane of the dragonfly are carried out. It is found that in the direction of the thickness the membrane was divided into three layers rather than a single entity as traditionally considered, and on the surfaces the membrane displays a random distribution rough microstructure that is composed of numerous nanometer scale columns coated by the cuticle wax secreted. The characteristics of the surface structure are measured and described. The mechanical properties of the membranes taken separately from the wings of live and dead dragonflies are investigated by the nanoindentation technique. The Young＇s moduli obtained here are approximately two times greater than the previous result, and the reasons that yield the difference are discussed.展开更多
Flapping-powered propulsion is used by many animals to locomote through air or water. Here we review recent experimental and numerical studies on self-propelled mechanical systems powered by a flapping motion. These s...Flapping-powered propulsion is used by many animals to locomote through air or water. Here we review recent experimental and numerical studies on self-propelled mechanical systems powered by a flapping motion. These studies improve our understanding of the mutual interaction between actively flapping bodies and surrounding fluids. The results obtained in these works provide not only new insights into biolocomotion but also useful information for the biomimetic design of artificial flyers and swimmers.展开更多
An improved electromechanical model of the RF MEMS(radio frequency microelec- tromechanical systems)switches is introduced,in which the effects of intrinsic residual stress from fabrication processes,axial stress due ...An improved electromechanical model of the RF MEMS(radio frequency microelec- tromechanical systems)switches is introduced,in which the effects of intrinsic residual stress from fabrication processes,axial stress due to stretching of beam,and fringing field are taken into account. Four dimensionless numbers are derived from the governing equation of the developed model.A semi- analytical method is developed to calculate the behavior of the RF MEMS switches.Subsequently the influence of the material and geometry parameters on the behavior of the structure is analyzed and compared,and the corresponding analysis with the dimensionless numbers is conducted too.The quantitative relationship between the presented parameters and the critical pull-in voltage is obtained, and the relative importance of those parameters is given.展开更多
A general formulation is developed for the contact behavior of a finite circular plate with a tensionless elastic foundation. The gap distance between the plate and elastic foundation is incorporated as an important p...A general formulation is developed for the contact behavior of a finite circular plate with a tensionless elastic foundation. The gap distance between the plate and elastic foundation is incorporated as an important parameter. Unlike the previous models with zero gap distance and large/infinite plate radius, which assumes the lift-off/separation of a flexural plate from its supporting elastic foundation, this study shows that lift-off may not occur. The results show how the contact area varies with the plate radius, boundary conditions and gap distance. When the plate radius becomes large enough and the gap distance is reduced to zero, the converged contact radius close to the previous ones is obtained.展开更多
Molecular dynamics （MD） simulations are performed to study adhesion and peeling of a short fragment of single strand DNA （ssDNA） molecule from a graphite surface. The critical peel-off force is found to depend on ...Molecular dynamics （MD） simulations are performed to study adhesion and peeling of a short fragment of single strand DNA （ssDNA） molecule from a graphite surface. The critical peel-off force is found to depend on both the peeling angle and the elasticity of ssDNA. For the short ssDNA strand under investigation, we show that the simulation results can be explained by a continuum model of an adhesive elastic band on substrate. The analysis suggests that it is often the peak value, rather than the mean value, of adhesion energy which determines the peeling of a nanoscale material.展开更多
Because of the interaction between film and substrate,the film buckling stress can vary significantly,depending on the delamination geometry,the film and substrate mechanical properties.The Mexican hat effect indicate...Because of the interaction between film and substrate,the film buckling stress can vary significantly,depending on the delamination geometry,the film and substrate mechanical properties.The Mexican hat effect indicates such interaction.An analytical method is presented,and related dimensional analysis shows that a single dimensionless parameter can effectively evaluate the effect.展开更多
In an indentation test,the effective Young's modulus of a film/substrate bilayer heterostructure varies with the indentation depth,a phenomenon known as the substrate effect.In previous studies investigating this,...In an indentation test,the effective Young's modulus of a film/substrate bilayer heterostructure varies with the indentation depth,a phenomenon known as the substrate effect.In previous studies investigating this,only the Young's modulus of the film was unknown.Once the effective Young's modulus of a film/substrate structure is determined at a given contact depth,the Young's modulus of the film can be uniquely determined,i.e.,there is a one-to-one relation between the Young's modulus of the film and the film/substrate effective Young's modulus.However,at times it is extremely challenging or even impossible to measure the film thickness.Furthermore,the precise definition of the layer/film thickness for a two-dimensional material can be problematic.In the current study,therefore,the thickness of the film and its Young's modulus are treated as two unknowns that must be determined.Unlike the case with one unknown,there are infinite combinations of film thickness and Young's modulus which can yield the same effective Young's modulus for the film/substrate.An inverse problem is formulated and solved to extract the Young's modulus and thickness of the film from the indentation depth-load curve.The accuracy and robustness of the inverse problem-solving method are also demonstrated.展开更多
Structural relaxation through isothermal annealing at temperature below glass transition is conducted on Zr46.75 Ti8.25 CU7.5Ni11Be27.5 （Vitreloy-4） bulk metallic glass. Defect concentration is correlated with the a...Structural relaxation through isothermal annealing at temperature below glass transition is conducted on Zr46.75 Ti8.25 CU7.5Ni11Be27.5 （Vitreloy-4） bulk metallic glass. Defect concentration is correlated with the annealing time t according to differential scanning calorimetry thermalgrams. The effects of structural relaxation on mechanical properties and deformation behaviour are investigated by using instrumented nanoindentation. It is found that the as-cast alloy exhibits pronounced serration flow during the loading process of nanoindentation, and the size and number of the serrations decrease with the annealing time. The change of the deformation behaviour with structural relaxation is explained using a free volume model.展开更多
In 1805, Thomas Young was the first to propose an equation(Young's equation) to predict the value of the equilibrium contact angle of a liquid on a solid. On the basis of our predecessors, we further clarify that ...In 1805, Thomas Young was the first to propose an equation(Young's equation) to predict the value of the equilibrium contact angle of a liquid on a solid. On the basis of our predecessors, we further clarify that the contact angle in Young's equation refers to the super-nano contact angle. Whether the equation is applicable to nanoscale systems remains an open question. Zhu et al. [College Phys. 4 7(1985)] obtained the most simple and convenient approximate formula, known as the Zhu–Qian approximate formula of Young's equation. Here, using molecular dynamics simulation, we test its applicability for nanodrops. Molecular dynamics simulations are performed on argon liquid cylinders placed on a solid surface under a temperature of 90 K, using Lennard–Jones potentials for the interaction between liquid molecules and between a liquid molecule and a solid molecule with the variable coefficient of strength a. Eight values of a between 0.650 and 0.825 are used. By comparison of the super-nano contact angles obtained from molecular dynamics simulation and the Zhu–Qian approximate formula of Young's equation, we find that it is qualitatively applicable for nanoscale systems.展开更多
Tendril-bearing climbing plants must recur to the tendril helices with chiral perversion or dual chirality for climbing and to obtain sun exposure. Despite researchers' prolonged fascination with climbing tendrils...Tendril-bearing climbing plants must recur to the tendril helices with chiral perversion or dual chirality for climbing and to obtain sun exposure. Despite researchers' prolonged fascination with climbing tendrils since Darwin's time and even earlier, why the soft and slender tendrils can bear heavy loads such as the self-weight of a plant or additional load caused by rain remains elusive. In this paper, we take towel gourd tendrils as an example and investigate the macroscopic and microscopic mechanical behaviors of tendrils through experiments and simulations. Our study indicates that the tendril filament exhibits rubber-like hyperelastic behaviors and can particularly endure large elongation, which is mainly attributed to the superelasticity of the cellulose fibril helix contained in the cell wall. Combination of the tendril helical structure with dual chirality or chiral perversion at a macroscale and a cellulose filament helix at a subcellular level creates superior elasticity for biological species relying on support and climbing. This study provides deep insight into the structure-property relationship of climbing tendrils, and the relationship is useful for the bioinspired design of composite systems with superior elasticity.展开更多
A scale-similarity model of a two-point two-time Lagrangian velocity correlation(LVC) was originally developed for the relative dispersion of tracer particles in isotropic turbulent flows(HE, G. W., JIN, G. D., and ZH...A scale-similarity model of a two-point two-time Lagrangian velocity correlation(LVC) was originally developed for the relative dispersion of tracer particles in isotropic turbulent flows(HE, G. W., JIN, G. D., and ZHAO, X. Scale-similarity model for Lagrangian velocity correlations in isotropic and stationary turbulence. Physical Review E, 80, 066313(2009)). The model can be expressed as a two-point Eulerian space correlation and the dispersion velocity V. The dispersion velocity denotes the rate at which one moving particle departs from another fixed particle. This paper numerically validates the robustness of the scale-similarity model at high Taylor micro-scale Reynolds numbers up to 373, which are much higher than the original values(R_λ = 66, 102). The effect of the Reynolds number on the dispersion velocity in the scale-similarity model is carefully investigated. The results show that the scale-similarity model is more accurate at higher Reynolds numbers because the two-point Lagrangian velocity correlations with different initial spatial separations collapse into a universal form compared with a combination of the initial separation and the temporal separation via the dispersion velocity.Moreover, the dispersion velocity V normalized by the Kolmogorov velocity V_η ≡ η/τ_η in which η and τ_η are the Kolmogorov space and time scales, respectively, scales with the Reynolds number R_λ as V/V_η ∝ R_λ^(1.39) obtained from the numerical data.展开更多
In view of the continued disputes on the fundamental question of whether the surface tension of a vapour bubble in liquid argon increases, or decreases, or remains unchanged with the increase of curvature radius, a cy...In view of the continued disputes on the fundamental question of whether the surface tension of a vapour bubble in liquid argon increases, or decreases, or remains unchanged with the increase of curvature radius, a cylindrical vapour bubble of argon is studied by molecular dynamics simulation in this paper instead of spherical vapour bubble so as to reduce the statistical error. So far, the surface tension of the cylindrical vapour bubble has not been studied by molecular dynamics simulation in the literature. Our results show that the surface tension decreases with radius increasing. By fitting the Tolman equation with our data, the Tolman length σ = -0.6225 sigma is given under cut-off radius 2.5σ, where σ = 0.3405 nm is the diameter of an argon atom. The Tolman length of Ar being negative is affirmed and the Tolman length of Ar being approximately zero given in the literature is negated, and it is pointed out that this error is attributed to the application of the inapplicable empirical equation of state and the neglect of the difference between surface tension and an equimolar surface.展开更多
The effects of the diameters of single-walled carbon nanotubes （SWCNTs） （7.83A to 27.40A） and temperature （20 K-45 K） on the equilibrium structure of an argon cluster are systematically studied by molecular dyna...The effects of the diameters of single-walled carbon nanotubes （SWCNTs） （7.83A to 27.40A） and temperature （20 K-45 K） on the equilibrium structure of an argon cluster are systematically studied by molecular dynamics simulation with consideration of the SWCNTs to be fixed. Since the diameters of SWCNTs with different chiralities increase when temperature is fixed at 20 K, the equilibrium structures of the argon cluster transform from monoatomic chains to helical and then to multishell coaxial cylinders. Chirality has almost no noticeable influence on these cylindrosymmetric structures. The effects of temperature and a non-equilibrium sudden heating process on the structures of argon clusters in SWCNTs are also studied by molecular dynamics simulation.展开更多
The quasistatic nanoindentation process of a spherical indenter in a single crystal copper is investigated with the molecular statistical thermodynamics(MST)method based on the embedded atom method(EAM)potential.The i...The quasistatic nanoindentation process of a spherical indenter in a single crystal copper is investigated with the molecular statistical thermodynamics(MST)method based on the embedded atom method(EAM)potential.The indentation modulus obtained in the MST simulation is 129.9 GPa,which agrees well with the theoretical prediction(129GPa).In the elastic regime,the obtained maximum displacement of the indenter is two times the contact depth and the contact area is qualitatively proportional to the contact depth,which agrees well with Hertzian elastic theory of contact.The MST simulation can reproduce the nucleation of dislocation as well.Moreover,the efficiency of the MST method is about 8 times higher than that of traditional MD simulations.展开更多
A critical Blot number, which determines both the sensitivity of spherical ceramics to quenching and the durations of the temperature-wave propagation and the thermal stresses in the ceramics subjected to thermal shoc...A critical Blot number, which determines both the sensitivity of spherical ceramics to quenching and the durations of the temperature-wave propagation and the thermal stresses in the ceramics subjected to thermal shock, is theoretically obtained. The results prove that once the Blot number of a ceramic sphere is greater than the critical number, its thermal shock failure will be such a rapid process that the failure only occurs in the initial regime of heat conduction, whereas the thermal shock failure of the ceramic sphere is uncertain in the course of heat conduction. The presented results provide a guide to the selection of the ceramics applied in the thermostructural engineering with thermal shock.展开更多
We investigate the plastic deformation and constitutive behaviour of bulk metallic glasses （BMGs）. A dimensionless Deborah number DeiD = tr/ti is proposed to characterize the rate effect in BMGs, where tr is the str...We investigate the plastic deformation and constitutive behaviour of bulk metallic glasses （BMGs）. A dimensionless Deborah number DeiD = tr/ti is proposed to characterize the rate effect in BMGs, where tr is the structural relaxing characteristic time of BMGs under shear load, ti is the macroscopic imposed characteristic time of applied stress or the characteristic time of macroscopic deformation. The results demonstrate that the modified free volume model can characterize the strain rate effect in BMGs effectively.展开更多
For the determination of surface tension of liquid droplets by molecular dynamics simulations, the most timeconsuming part is the calculation of pressure tensor in the transition layer, which makes it difficult to enh...For the determination of surface tension of liquid droplets by molecular dynamics simulations, the most timeconsuming part is the calculation of pressure tensor in the transition layer, which makes it difficult to enhance the precision of the computation. A new method for the calculation of surface tension of liquid droplets to reduce the calculation quantity of pressure tensor in transition layer to the minimum is proposed in this paper. Two thousand particles are taken as example to show how to carry out our scheme.展开更多
基金financial support from the Young Scientist Project of the National Key Research and Development Program of China(Grant No.2021YFC2900600)Beijing Nova Program(Grant No.20220484057)+1 种基金The authors acknowledge the financial support from the Natural Sciences and Engineering Research Council of Canada(Grant No.RGPIN-2018-06902)industrial partners of the Research Institute on Mines and the Environment(RIME UQAT-Polytechnique:https://irme.ca/en/).
文摘Underhand cut-and-fill mining has been widely used in underground mining operations,especially when the rock mass or orebody is of poor quality or prone to rockburst due to high stress.In such cases,mining workers should carry out all production activities under the cemented backfill roof or sill mat instead of a highly fractured and unstable rock roof or a strong rock roof with a high potential of rockburst.Therefore,the stability and required strength of the sill mat are critical issues for mining engineers.In 1991,Mitchell considered that sill mat could fail by caving,sliding,rotation,and flexure.Mitchell also proposed an analytical solution to determine the minimum required strength of the sill mat for each type of failure based on two stiff or immobile rock walls.However,recent publications using numerical modeling and field measurements indicate that the compressive stresses in the sill mat induced by rock wall closure due to a stope excavation beneath the sill mat can be significant.It is thus highly necessary to investigate the required strength of the sill mat by considering rock wall closure.In this study,the crushing failure of sill mat due to rock wall closure generated by underground excavation and a new failure mode called"crushing and caving”is revealed by numerical modeling.An analytical solution corresponding to each failure mode is then developed to estimate the minimum required cohesion(cmin)of the sill mat.A criterion is also proposed to determine if the sill mat fails by crushing or crushing-and-caving failure.The proposed analytical solution does not involve any correction coefficients.The validity of the proposed analytical solution is demonstrated by numerical modeling.The proposed analytical solution can thus be employed to predict the cmin of sill mat subjected to wall closure generated by underlying stope excavation.
基金financial support from the Young Scientist Project of the National Key Research and Development Program of China(Grant No.2021YFC2900600)the Beijing Nova Program(Grant No.20220484057)+1 种基金The authors acknowledge the financial support from the Natural Sciences and Engineering Research Council of Canada(Grant No.RGPIN-2018-06902)industrial partners of the Research Institute on Mines and the Environment(RIME UQAT-Polytechnique:https://irme.ca/en/).
文摘In the underhand cut-and-fill mining method,a sill mat(i.e.an artificial horizontal pillar)constructed by cemented backfill is essential to prevent mine workers from being directly exposed under problematic rock roofs.A critical issue is to determine the minimum required strength of the sill mat to ensure a safe and cost-effective design.Until now,Mitchell’s analytical solution is the only available option,considering two stiff and immobile rock walls.Unavoidable rock wall closure associated with stope excavation below the sill mat was neglected.This,along with other undefined parameters,explains why Mitchell’s solution is rarely used in sill mat design.A new analytical solution for determining the minimum required strength of the sill mat accounting for wall closure is necessary.In this study,a closed-form analytical solution for estimating rock wall closure generated by stope excavation below a sill mat is developed by using Salamon’s and Flamant’s models.The proposed analytical solution does not contain any coefficients of correction or calibration.Despite several assumptions(or somewhat of oversimplifications)necessary to render a simple analytical solution possible,good agreements are obtained between the rock wall closures predicted by applying the proposed analytical solution and those obtained numerically with FLAC3D for many cases with arbitrarily chosen geometrical and material parameters.The proposed analytical solution is therefore validated and can be used to evaluate the rock wall closure generated by stope excavation below a sill mat.
基金Project supported by the National Natural Science Foundation of China(Nos.91752118,11672305,11232011,and 11572331)the Strategic Priority Research Program(No.XDB22040104)the Key Research Program of Frontier Sciences of the Chinese Academy of Sciences(No.QYZDJ-SSWSYS002)
文摘A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equation is integrated along the wall-normal direction to link the tangential component of the effective body force for the IB method to the wall shear stress predicted by the wall model;(ii) a set of Lagrangian points near the wall are introduced to compute the normal component of the effective body force for the IB method by reconstructing the normal component of the velocity. This novel method will be a classical direct-forcing IB method if the grid is fine enough to resolve the flow near the wall. The method is used to simulate the flows around the DARPA SUBOFF model. The results obtained are well comparable to the measured experimental data and wall-resolved LES results.
基金The project supported by the National Natural Science Foundation of China(10372102 and 10672164)
文摘Detailed investigations on the microstructure and the mechanical properties of the wing membrane of the dragonfly are carried out. It is found that in the direction of the thickness the membrane was divided into three layers rather than a single entity as traditionally considered, and on the surfaces the membrane displays a random distribution rough microstructure that is composed of numerous nanometer scale columns coated by the cuticle wax secreted. The characteristics of the surface structure are measured and described. The mechanical properties of the membranes taken separately from the wings of live and dead dragonflies are investigated by the nanoindentation technique. The Young＇s moduli obtained here are approximately two times greater than the previous result, and the reasons that yield the difference are discussed.
基金supported by the Chinese Academy of Sciences (Grants KJCX-SW-L08, KJCX3-SYW-S01)the National Natural Science Foundation of China (Grants 11021262, 11023001, 11232011, 11372331)
文摘Flapping-powered propulsion is used by many animals to locomote through air or water. Here we review recent experimental and numerical studies on self-propelled mechanical systems powered by a flapping motion. These studies improve our understanding of the mutual interaction between actively flapping bodies and surrounding fluids. The results obtained in these works provide not only new insights into biolocomotion but also useful information for the biomimetic design of artificial flyers and swimmers.
基金The project supported by the National Natural Science Foundation of China,the Chinese Academy of Sciences,the RGC/NSFC Joint Research Scheme (N-HKUST 601/01)the Joint Laboratory of Microsystems
文摘An improved electromechanical model of the RF MEMS(radio frequency microelec- tromechanical systems)switches is introduced,in which the effects of intrinsic residual stress from fabrication processes,axial stress due to stretching of beam,and fringing field are taken into account. Four dimensionless numbers are derived from the governing equation of the developed model.A semi- analytical method is developed to calculate the behavior of the RF MEMS switches.Subsequently the influence of the material and geometry parameters on the behavior of the structure is analyzed and compared,and the corresponding analysis with the dimensionless numbers is conducted too.The quantitative relationship between the presented parameters and the critical pull-in voltage is obtained, and the relative importance of those parameters is given.
基金supported by the National Natural Science Foundation of China (11021262 and 11023001)Chinese Academyof Sciences (KJCX2-EW-L03)
文摘A general formulation is developed for the contact behavior of a finite circular plate with a tensionless elastic foundation. The gap distance between the plate and elastic foundation is incorporated as an important parameter. Unlike the previous models with zero gap distance and large/infinite plate radius, which assumes the lift-off/separation of a flexural plate from its supporting elastic foundation, this study shows that lift-off may not occur. The results show how the contact area varies with the plate radius, boundary conditions and gap distance. When the plate radius becomes large enough and the gap distance is reduced to zero, the converged contact radius close to the previous ones is obtained.
基金The project supported by the Distinguished Young Scholar Fund of NSFC(10225209)key project from the Chinese Academy of Sciences(KJCX2-SW-L2)
文摘Molecular dynamics （MD） simulations are performed to study adhesion and peeling of a short fragment of single strand DNA （ssDNA） molecule from a graphite surface. The critical peel-off force is found to depend on both the peeling angle and the elasticity of ssDNA. For the short ssDNA strand under investigation, we show that the simulation results can be explained by a continuum model of an adhesive elastic band on substrate. The analysis suggests that it is often the peak value, rather than the mean value, of adhesion energy which determines the peeling of a nanoscale material.
基金supported by the National Natural Science Foundation of China(11023001 and 11372321)
文摘Because of the interaction between film and substrate,the film buckling stress can vary significantly,depending on the delamination geometry,the film and substrate mechanical properties.The Mexican hat effect indicates such interaction.An analytical method is presented,and related dimensional analysis shows that a single dimensionless parameter can effectively evaluate the effect.
基金This research was supported by the National Natural Science Foundation of China(11772335,21622304,61674045,and 21203038)by the Ministry of Science and Technology(MOST)of China(2016YFA0200700).Z.H.Cheng was supported by the Distinguished Technical Talents Project and the Youth Innovation Promotion Association of Chinese Academy of Sciences.
文摘In an indentation test,the effective Young's modulus of a film/substrate bilayer heterostructure varies with the indentation depth,a phenomenon known as the substrate effect.In previous studies investigating this,only the Young's modulus of the film was unknown.Once the effective Young's modulus of a film/substrate structure is determined at a given contact depth,the Young's modulus of the film can be uniquely determined,i.e.,there is a one-to-one relation between the Young's modulus of the film and the film/substrate effective Young's modulus.However,at times it is extremely challenging or even impossible to measure the film thickness.Furthermore,the precise definition of the layer/film thickness for a two-dimensional material can be problematic.In the current study,therefore,the thickness of the film and its Young's modulus are treated as two unknowns that must be determined.Unlike the case with one unknown,there are infinite combinations of film thickness and Young's modulus which can yield the same effective Young's modulus for the film/substrate.An inverse problem is formulated and solved to extract the Young's modulus and thickness of the film from the indentation depth-load curve.The accuracy and robustness of the inverse problem-solving method are also demonstrated.
基金Support by the National Natural Science Foundation of China under Grant Nos 50571109, 10572142 and 10432050, the Knowledge Innovation Project of Chinese Academy of Sciences, and the Natural Science Fund for 0utstanding Youth in Hunan Province under Grant No 02JJYB010.
文摘Structural relaxation through isothermal annealing at temperature below glass transition is conducted on Zr46.75 Ti8.25 CU7.5Ni11Be27.5 （Vitreloy-4） bulk metallic glass. Defect concentration is correlated with the annealing time t according to differential scanning calorimetry thermalgrams. The effects of structural relaxation on mechanical properties and deformation behaviour are investigated by using instrumented nanoindentation. It is found that the as-cast alloy exhibits pronounced serration flow during the loading process of nanoindentation, and the size and number of the serrations decrease with the annealing time. The change of the deformation behaviour with structural relaxation is explained using a free volume model.
基金Project supported by the National Natural Science Foundation of China(Grant No.11072242)the Key Scientific Studies Program of Hebei Province Higher Education Institute,China(Grant No.ZD2018301)Cangzhou National Science Foundation,China(Grant No.177000001)
文摘In 1805, Thomas Young was the first to propose an equation(Young's equation) to predict the value of the equilibrium contact angle of a liquid on a solid. On the basis of our predecessors, we further clarify that the contact angle in Young's equation refers to the super-nano contact angle. Whether the equation is applicable to nanoscale systems remains an open question. Zhu et al. [College Phys. 4 7(1985)] obtained the most simple and convenient approximate formula, known as the Zhu–Qian approximate formula of Young's equation. Here, using molecular dynamics simulation, we test its applicability for nanodrops. Molecular dynamics simulations are performed on argon liquid cylinders placed on a solid surface under a temperature of 90 K, using Lennard–Jones potentials for the interaction between liquid molecules and between a liquid molecule and a solid molecule with the variable coefficient of strength a. Eight values of a between 0.650 and 0.825 are used. By comparison of the super-nano contact angles obtained from molecular dynamics simulation and the Zhu–Qian approximate formula of Young's equation, we find that it is qualitatively applicable for nanoscale systems.
基金the National Natural Science Foundation of China (Grants 11872273, 11472191, 11602163, and 11672297)the Major Program of the National Science Foundation of China (Grant 11890683)+1 种基金the Opening Fund of State Key Laboratory of Nonlinear Mechanicsthe Australian Endeavour Research Fellowship.
文摘Tendril-bearing climbing plants must recur to the tendril helices with chiral perversion or dual chirality for climbing and to obtain sun exposure. Despite researchers' prolonged fascination with climbing tendrils since Darwin's time and even earlier, why the soft and slender tendrils can bear heavy loads such as the self-weight of a plant or additional load caused by rain remains elusive. In this paper, we take towel gourd tendrils as an example and investigate the macroscopic and microscopic mechanical behaviors of tendrils through experiments and simulations. Our study indicates that the tendril filament exhibits rubber-like hyperelastic behaviors and can particularly endure large elongation, which is mainly attributed to the superelasticity of the cellulose fibril helix contained in the cell wall. Combination of the tendril helical structure with dual chirality or chiral perversion at a macroscale and a cellulose filament helix at a subcellular level creates superior elasticity for biological species relying on support and climbing. This study provides deep insight into the structure-property relationship of climbing tendrils, and the relationship is useful for the bioinspired design of composite systems with superior elasticity.
基金Project supported by the Science Challenge Program(No.TZ2016001)the National Natural Science Foundation of China(Nos.11472277,11572331,11232011,and 11772337)+1 种基金the Strategic Priority Research Program,Chinese Academy of Sciences(No.XDB22040104)the Key Research Program of Frontier Sciences,Chinese Academy of Sciences(No.QYZDJ-SSW-SYS002)
文摘A scale-similarity model of a two-point two-time Lagrangian velocity correlation(LVC) was originally developed for the relative dispersion of tracer particles in isotropic turbulent flows(HE, G. W., JIN, G. D., and ZHAO, X. Scale-similarity model for Lagrangian velocity correlations in isotropic and stationary turbulence. Physical Review E, 80, 066313(2009)). The model can be expressed as a two-point Eulerian space correlation and the dispersion velocity V. The dispersion velocity denotes the rate at which one moving particle departs from another fixed particle. This paper numerically validates the robustness of the scale-similarity model at high Taylor micro-scale Reynolds numbers up to 373, which are much higher than the original values(R_λ = 66, 102). The effect of the Reynolds number on the dispersion velocity in the scale-similarity model is carefully investigated. The results show that the scale-similarity model is more accurate at higher Reynolds numbers because the two-point Lagrangian velocity correlations with different initial spatial separations collapse into a universal form compared with a combination of the initial separation and the temporal separation via the dispersion velocity.Moreover, the dispersion velocity V normalized by the Kolmogorov velocity V_η ≡ η/τ_η in which η and τ_η are the Kolmogorov space and time scales, respectively, scales with the Reynolds number R_λ as V/V_η ∝ R_λ^(1.39) obtained from the numerical data.
基金Project supported by the National Natural Science Foundation of China (Grant No. 11072242)
文摘In view of the continued disputes on the fundamental question of whether the surface tension of a vapour bubble in liquid argon increases, or decreases, or remains unchanged with the increase of curvature radius, a cylindrical vapour bubble of argon is studied by molecular dynamics simulation in this paper instead of spherical vapour bubble so as to reduce the statistical error. So far, the surface tension of the cylindrical vapour bubble has not been studied by molecular dynamics simulation in the literature. Our results show that the surface tension decreases with radius increasing. By fitting the Tolman equation with our data, the Tolman length σ = -0.6225 sigma is given under cut-off radius 2.5σ, where σ = 0.3405 nm is the diameter of an argon atom. The Tolman length of Ar being negative is affirmed and the Tolman length of Ar being approximately zero given in the literature is negated, and it is pointed out that this error is attributed to the application of the inapplicable empirical equation of state and the neglect of the difference between surface tension and an equimolar surface.
基金Project supported by the National Natural Science Foundation of China(Grant No.11072242)
文摘The effects of the diameters of single-walled carbon nanotubes （SWCNTs） （7.83A to 27.40A） and temperature （20 K-45 K） on the equilibrium structure of an argon cluster are systematically studied by molecular dynamics simulation with consideration of the SWCNTs to be fixed. Since the diameters of SWCNTs with different chiralities increase when temperature is fixed at 20 K, the equilibrium structures of the argon cluster transform from monoatomic chains to helical and then to multishell coaxial cylinders. Chirality has almost no noticeable influence on these cylindrosymmetric structures. The effects of temperature and a non-equilibrium sudden heating process on the structures of argon clusters in SWCNTs are also studied by molecular dynamics simulation.
基金Supported by the National Natural Science Foundation of China under Grant Nos 10932011,10772181,11021262,10732090 and 10772012the National Basic Research Program of China under Grant No 2007CB814803.
文摘The quasistatic nanoindentation process of a spherical indenter in a single crystal copper is investigated with the molecular statistical thermodynamics(MST)method based on the embedded atom method(EAM)potential.The indentation modulus obtained in the MST simulation is 129.9 GPa,which agrees well with the theoretical prediction(129GPa).In the elastic regime,the obtained maximum displacement of the indenter is two times the contact depth and the contact area is qualitatively proportional to the contact depth,which agrees well with Hertzian elastic theory of contact.The MST simulation can reproduce the nucleation of dislocation as well.Moreover,the efficiency of the MST method is about 8 times higher than that of traditional MD simulations.
基金Supported by the National Natural Science Foundations of China under Grant No 90716004.
文摘A critical Blot number, which determines both the sensitivity of spherical ceramics to quenching and the durations of the temperature-wave propagation and the thermal stresses in the ceramics subjected to thermal shock, is theoretically obtained. The results prove that once the Blot number of a ceramic sphere is greater than the critical number, its thermal shock failure will be such a rapid process that the failure only occurs in the initial regime of heat conduction, whereas the thermal shock failure of the ceramic sphere is uncertain in the course of heat conduction. The presented results provide a guide to the selection of the ceramics applied in the thermostructural engineering with thermal shock.
基金Supported by the National Natural Science Foundation of China under Grant Nos 10725211 and 10721202, and the Knowledge Innovation Project of Chinese Academy of Sciences under Grant No KJCX-SW-L08, and the Doctorial Start-up Fund of Hunan University of Science and Technology （E50840）.
文摘We investigate the plastic deformation and constitutive behaviour of bulk metallic glasses （BMGs）. A dimensionless Deborah number DeiD = tr/ti is proposed to characterize the rate effect in BMGs, where tr is the structural relaxing characteristic time of BMGs under shear load, ti is the macroscopic imposed characteristic time of applied stress or the characteristic time of macroscopic deformation. The results demonstrate that the modified free volume model can characterize the strain rate effect in BMGs effectively.
基金Project supported by the National Natural Science Foundation of China(Grant No.10772189)the Knowledge Innovation Program of Chinese Academy of Sciences
文摘For the determination of surface tension of liquid droplets by molecular dynamics simulations, the most timeconsuming part is the calculation of pressure tensor in the transition layer, which makes it difficult to enhance the precision of the computation. A new method for the calculation of surface tension of liquid droplets to reduce the calculation quantity of pressure tensor in transition layer to the minimum is proposed in this paper. Two thousand particles are taken as example to show how to carry out our scheme.
