In this paper, we combine the pseudo arc-length numerical method with the mathematical model of multiphase compressible flow for simulating the shock wave interaction with a deformable particle. Firstly, an arc-length...In this paper, we combine the pseudo arc-length numerical method with the mathematical model of multiphase compressible flow for simulating the shock wave interaction with a deformable particle. Firstly, an arc-length parameter is introduced to weaken the discontinuous singularity of governing equations, and an efficient pseudo arc-length numerical method of multiphase compressible flow is proposed. Then the accuracy and adaptive moving mesh property of this algorithm are tested. Finally, the multiphase pseudo arc-length numerical method is applied to the problem of interaction between shock wave and the deformable particle. Through the flow flied change and data analysis of key points, it can be found the complex wave structures are presented after the interactions between the planar incident shock wave and the metal particle, and all these wave interactions lead to the movement and deformation of metal particle, and then the deformed particle will affect the transmitted shock wave back. According to the discussion, the deformation of particle and shock wave propagation in the particle are determined by the shock wave impedance of each medium and shock speed, so the interaction between shock wave and the deformable particle can be studied on the basis of physical properties of explosive mediums.展开更多
A particle mapping transportation algorithm was proposed on the basis of the particle-in-cell method.The particles with rectangular influence domains were employed in the transportation algorithm to reduce the numeric...A particle mapping transportation algorithm was proposed on the basis of the particle-in-cell method.The particles with rectangular influence domains were employed in the transportation algorithm to reduce the numerical fluctuations.Based on the error analysis in the process of particle motion computation,a prediction-correction algorithm was introduced to improve the computational accuracy.Furthermore,the performance of the particle mapping transportation method was evaluated by using the rotation,the slotted disk and the shear advection tests,and the results were compared with other interface reconstruction methods.Finally,the hemispherical projectile penetration into a steel target was numerically simulated.The results showed that the proposed method produced less numerical fluctuations and exhibited clear material interfaces,which indicated that it is accurate and effective.展开更多
In this study, a complex multi-body structure was proposed, and the mechanism for the dynamic response of the structure under explosive driving was investigated by using the Lagrange equations of the second kind. An i...In this study, a complex multi-body structure was proposed, and the mechanism for the dynamic response of the structure under explosive driving was investigated by using the Lagrange equations of the second kind. An initial value subject to explosion loading was analyzed to develop the theoretical model of the dynamic response, and the centroid trajectory of three different structural shapes was solved. To verify the accuracy of the theoretical model, numerical simulation via finite element analysis within LS-DYNA and a dynamic experiment were conducted, and the consistent dynamic response process of the multi-body structure was obtained. In addition, the dynamic response time of the multi-body structure under different explosion loading conditions was calculated by the theoretical model, numerical simulation, and experimental investigation. It was found that the increased opening charge mass reduces the dynamic response time.展开更多
A three-dimensional(3D) embedded Eulerian-Lagrangian method is proposed to simulate the 3D fluid-structure interaction(FSI) problems subjected to explosion and impact loading. This method achieves a high-quality calcu...A three-dimensional(3D) embedded Eulerian-Lagrangian method is proposed to simulate the 3D fluid-structure interaction(FSI) problems subjected to explosion and impact loading. This method achieves a high-quality calculation of fluid and structure deformation by adding Lagrangian particles to Eulerian grids. The overall computational domain is solved by the Eulerian method, and the Lagrangian particles with specified volume and influence domains are used to track structural deformations. The bidirectional mapping of physical quantities is achieved using the weighted average of the influence domain, which are based on the topological relationship between Eulerian grids and Lagrangian particles. Then, the data dependence solution and parallel algorithm realization are presented for the large-scale numerical calculations of explosion and impact problems. Additionally, the corresponding parallel program is developed based on the message passing interface(MPI) standard, and the parallel efficiency of parallel hydrocode are tested. The numerical results of typical explosion and impact problems are compared with corresponding experimental data to verify the effectiveness of the method. These comparisons show that the embedded EulerianLagrangian method successfully combine the advantages of both the Eulerian and Lagrangian methods to efficiently calculate the processes of large deformation and dynamic damage to the materials. The results presented in this work provide a useful reference point for further research on explosion and impact problems.展开更多
In this study,we mainly focus on the structural morphology and inter-atomic bonding state of tribofilms resulting from a highly-hydrogenated amorphous carbon(a-C:H) film in order to ascertain the underlying mechanisms...In this study,we mainly focus on the structural morphology and inter-atomic bonding state of tribofilms resulting from a highly-hydrogenated amorphous carbon(a-C:H) film in order to ascertain the underlying mechanisms for its superlubric behavior(i.e.,less than 0.01 friction coefficient).Specifically,we achieved superlubricity(i.e.,friction coefficients of down to 0.003) with this film in dry nitrogen and argon atmospheres especially when the tribo-pair is made of an a-C:H coated Si disk sliding against an a-C:H coated steel ball,while the a-C:H coated disk against uncoated ball does not provide superlubricity.We also found that the state of superlubricity is more stable in argon than in nitrogen and the formation of a smooth and uniformly-thick carbonaceous tribofilm appears to be one of the key factors for the realization of such superlubricity.Besides,the interfacial morphology of sliding test pairs and the atomic-scale bond structure of the carbon-based tribofilms also play an important role in the observed superlubric behavior of a-C:H films.Using Raman spectroscopy and high resolution transmission electron microscopy,we have compared the structural differences of the tribofilms produced on bare and a-C:H coated steel balls.For the a-C:H coated ball as mating material which provided superlow friction in argon,structural morphology of the tribofilm was similar or comparable to that of the original a-C:H coating;while for the bare steel ball,the sp^2-bonded C fraction in the tribofilm increased and a fingerprint-like nanocrystalline structure was detected by high resolution transmission electron microscopy(HRTEM).We also calculated the shear stresses for different tribofilms,and established a relationship between the magnitude of the shear stresses and the extent of sp^3-sp^2 phase transformation.展开更多
This paper presented a numerical approach to solving the problem of a flat-ended punch in contact with a half-space matrix embedded with multiple three dimensional arbitrary-shaped inhomogeneities.Based on the semi-an...This paper presented a numerical approach to solving the problem of a flat-ended punch in contact with a half-space matrix embedded with multiple three dimensional arbitrary-shaped inhomogeneities.Based on the semi-analytical method(SAM)and the equivalent inclusion method,numerical procedures were developed and the effects of inclusion shape and distribution were analyzed.Fast Fourier transform technique was implemented to accelerate the calculation of surface deformation and subsurface stress.Interactions of inter-inclusions and inclusion-matrix were taken into account.Numerical results showed the presence of inhomogeneities(i.e.,microstructures in solids)indeed had a great effect on local contact pressure and a strong disturbance to the subsurface stress field in the vicinity of inclusions.The effects were dependent on the shape and distribution of inclusions and inter-inclusion interactions.The physical significance of this study is to provide an insight into the relation between the material microstructure and its response to the external load,and the solution approach and procedures may find useful applications,for example,the analysis of fatigue and crack propagation for composite materials,prediction of stress field in solids containing material defects,and study of the mechanism of chemical-mechanical polish(CMP)for inhomogeneous materials,etc.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.11390363,11325209 and 11221202)
文摘In this paper, we combine the pseudo arc-length numerical method with the mathematical model of multiphase compressible flow for simulating the shock wave interaction with a deformable particle. Firstly, an arc-length parameter is introduced to weaken the discontinuous singularity of governing equations, and an efficient pseudo arc-length numerical method of multiphase compressible flow is proposed. Then the accuracy and adaptive moving mesh property of this algorithm are tested. Finally, the multiphase pseudo arc-length numerical method is applied to the problem of interaction between shock wave and the deformable particle. Through the flow flied change and data analysis of key points, it can be found the complex wave structures are presented after the interactions between the planar incident shock wave and the metal particle, and all these wave interactions lead to the movement and deformation of metal particle, and then the deformed particle will affect the transmitted shock wave back. According to the discussion, the deformation of particle and shock wave propagation in the particle are determined by the shock wave impedance of each medium and shock speed, so the interaction between shock wave and the deformable particle can be studied on the basis of physical properties of explosive mediums.
基金supported by the National Basic Research Program of China(Grant No.2010CB832706)the National Natural Science Foundation of China(Grant Nos.11032002 and 11172041)the Project of State Key Laboratory of Explosion Science and Technology(Grant No. ZDKT10-03c)
文摘A particle mapping transportation algorithm was proposed on the basis of the particle-in-cell method.The particles with rectangular influence domains were employed in the transportation algorithm to reduce the numerical fluctuations.Based on the error analysis in the process of particle motion computation,a prediction-correction algorithm was introduced to improve the computational accuracy.Furthermore,the performance of the particle mapping transportation method was evaluated by using the rotation,the slotted disk and the shear advection tests,and the results were compared with other interface reconstruction methods.Finally,the hemispherical projectile penetration into a steel target was numerically simulated.The results showed that the proposed method produced less numerical fluctuations and exhibited clear material interfaces,which indicated that it is accurate and effective.
基金supported by the National Natural Science Foundation of China(Grant Nos.11372046,11521062)
文摘In this study, a complex multi-body structure was proposed, and the mechanism for the dynamic response of the structure under explosive driving was investigated by using the Lagrange equations of the second kind. An initial value subject to explosion loading was analyzed to develop the theoretical model of the dynamic response, and the centroid trajectory of three different structural shapes was solved. To verify the accuracy of the theoretical model, numerical simulation via finite element analysis within LS-DYNA and a dynamic experiment were conducted, and the consistent dynamic response process of the multi-body structure was obtained. In addition, the dynamic response time of the multi-body structure under different explosion loading conditions was calculated by the theoretical model, numerical simulation, and experimental investigation. It was found that the increased opening charge mass reduces the dynamic response time.
基金supported by the National Natural Science Foundation of China (Grant Nos. 11822203 and 11532012)。
文摘A three-dimensional(3D) embedded Eulerian-Lagrangian method is proposed to simulate the 3D fluid-structure interaction(FSI) problems subjected to explosion and impact loading. This method achieves a high-quality calculation of fluid and structure deformation by adding Lagrangian particles to Eulerian grids. The overall computational domain is solved by the Eulerian method, and the Lagrangian particles with specified volume and influence domains are used to track structural deformations. The bidirectional mapping of physical quantities is achieved using the weighted average of the influence domain, which are based on the topological relationship between Eulerian grids and Lagrangian particles. Then, the data dependence solution and parallel algorithm realization are presented for the large-scale numerical calculations of explosion and impact problems. Additionally, the corresponding parallel program is developed based on the message passing interface(MPI) standard, and the parallel efficiency of parallel hydrocode are tested. The numerical results of typical explosion and impact problems are compared with corresponding experimental data to verify the effectiveness of the method. These comparisons show that the embedded EulerianLagrangian method successfully combine the advantages of both the Eulerian and Lagrangian methods to efficiently calculate the processes of large deformation and dynamic damage to the materials. The results presented in this work provide a useful reference point for further research on explosion and impact problems.
基金supported by the National Basic Research Program of China (Grant No.2011CB013404)National Natural Science Foundation of China(Grant Nos.51321092,51527901 and 51375010)
文摘In this study,we mainly focus on the structural morphology and inter-atomic bonding state of tribofilms resulting from a highly-hydrogenated amorphous carbon(a-C:H) film in order to ascertain the underlying mechanisms for its superlubric behavior(i.e.,less than 0.01 friction coefficient).Specifically,we achieved superlubricity(i.e.,friction coefficients of down to 0.003) with this film in dry nitrogen and argon atmospheres especially when the tribo-pair is made of an a-C:H coated Si disk sliding against an a-C:H coated steel ball,while the a-C:H coated disk against uncoated ball does not provide superlubricity.We also found that the state of superlubricity is more stable in argon than in nitrogen and the formation of a smooth and uniformly-thick carbonaceous tribofilm appears to be one of the key factors for the realization of such superlubricity.Besides,the interfacial morphology of sliding test pairs and the atomic-scale bond structure of the carbon-based tribofilms also play an important role in the observed superlubric behavior of a-C:H films.Using Raman spectroscopy and high resolution transmission electron microscopy,we have compared the structural differences of the tribofilms produced on bare and a-C:H coated steel balls.For the a-C:H coated ball as mating material which provided superlow friction in argon,structural morphology of the tribofilm was similar or comparable to that of the original a-C:H coating;while for the bare steel ball,the sp^2-bonded C fraction in the tribofilm increased and a fingerprint-like nanocrystalline structure was detected by high resolution transmission electron microscopy(HRTEM).We also calculated the shear stresses for different tribofilms,and established a relationship between the magnitude of the shear stresses and the extent of sp^3-sp^2 phase transformation.
基金supported by the National Basic Research Program of China(Grant Nos.2009CB724200,2011CB013404 and 2011CB706602)
文摘This paper presented a numerical approach to solving the problem of a flat-ended punch in contact with a half-space matrix embedded with multiple three dimensional arbitrary-shaped inhomogeneities.Based on the semi-analytical method(SAM)and the equivalent inclusion method,numerical procedures were developed and the effects of inclusion shape and distribution were analyzed.Fast Fourier transform technique was implemented to accelerate the calculation of surface deformation and subsurface stress.Interactions of inter-inclusions and inclusion-matrix were taken into account.Numerical results showed the presence of inhomogeneities(i.e.,microstructures in solids)indeed had a great effect on local contact pressure and a strong disturbance to the subsurface stress field in the vicinity of inclusions.The effects were dependent on the shape and distribution of inclusions and inter-inclusion interactions.The physical significance of this study is to provide an insight into the relation between the material microstructure and its response to the external load,and the solution approach and procedures may find useful applications,for example,the analysis of fatigue and crack propagation for composite materials,prediction of stress field in solids containing material defects,and study of the mechanism of chemical-mechanical polish(CMP)for inhomogeneous materials,etc.