The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of c...The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of curl interface underground or “v” figure complex model, it is difficult to realize. So it is important to forward the complex geoelectricity model. This paper takes two Maxwell’s vorticity equations as departure point, makes use of the principles of Yee’s space grid model theory and the basic principle finite difference time domain method, and deduces a GPR forward system of equation of two dimensional spaces. The Mur super absorbed boundary condition is adopted to solve the super strong reflection on the interceptive boundary when there is the forward simulation. And a self-made program is used to process forward simulation to two typical geoelectricity model.展开更多
A new general optimal principle of designing explicit finite difference method was obtained. Several applied cases were put forward to explain the uses of the principle. The validity of the principal was tested by a n...A new general optimal principle of designing explicit finite difference method was obtained. Several applied cases were put forward to explain the uses of the principle. The validity of the principal was tested by a numeric example.展开更多
In order to investigate the process of incremental sheet forming (ISF) through both experimental and numerical approaches, a three-dimensional elasto-plastic finite element model (FEM) was developed to simulate the pr...In order to investigate the process of incremental sheet forming (ISF) through both experimental and numerical approaches, a three-dimensional elasto-plastic finite element model (FEM) was developed to simulate the process and the simulated results were compared with those of experiment. The results of numerical simulations, such as the strain history and distribution, the stress state and distribution, sheet thickness distribution, etc, were discussed in details, and the influences of process parameters on these results were also analyzed. The simulated results of the radial strain and the thickness distribution are in good agreement with experimental results. The simulations reveal that the deformation is localized around the tool and constantly remains close to a plane strain state. With decreasing depth step, increasing tool diameter and wall inclination angle, the axial stress reduces, leading to less thinning and more homogeneous plastic strain and thickness distribution. During ISF, the plastic strain increases stepwise under the action of the tool. Each increase in plastic strain is accompanied by hydrostatic pressure, which explains why obtainable deformation using ISF exceeds the forming limits of conventional sheet forming.展开更多
Being aimed at the inside wall wrinkling and sinking phenomenon of palladium-yttrium alloy thin wall spiral tube used for preparation of high purity hydrogen, extraction of hydrogen isotope, and purification and separ...Being aimed at the inside wall wrinkling and sinking phenomenon of palladium-yttrium alloy thin wall spiral tube used for preparation of high purity hydrogen, extraction of hydrogen isotope, and purification and separation of hydrogen in the winding process, this article analyzed the reasons for above phenomena, established a numerical simulation model of winding process of above tube, using elastic-plastic Finite Element method analyzed the max. tensile stress and max. compression stress and their locations, thereby provides a theory base for the control of working forming course of thin wall spiral tube.展开更多
In this paper, the internal fluid motion of a jet system is described by the Navier Stokes mechanics equations. For the simulation of the motion, the penalty function finite element method is used, and the velocity ve...In this paper, the internal fluid motion of a jet system is described by the Navier Stokes mechanics equations. For the simulation of the motion, the penalty function finite element method is used, and the velocity vectors and stream function curves are obtained. Using the Prandtl theory, this paper derives the free jet velocity and the jet bunch width in a half-space, the latter of which is amended by experiment. The results obtained in this paper are applied to micro-type high pressure water jet cleaner and the ejector of rocket engine.展开更多
Two central schemes of finite difference (FD) up to different accuracy orders of space sampling step Dx (Fourth order and Sixth order respectively) were used to study the 1-D nonlinear P-wave propagation in the nonlin...Two central schemes of finite difference (FD) up to different accuracy orders of space sampling step Dx (Fourth order and Sixth order respectively) were used to study the 1-D nonlinear P-wave propagation in the nonlinear solid media by the numerical method. Distinctly different from the case of numerical modeling of linear elastic wave, there may be several difficulties in the numerical treatment to the nonlinear partial differential equation, such as the steep gradients, shocks and unphysical oscillations. All of them are the great obstacles to the stability and conver-gence of numerical calculation. Fortunately, the comparative study on the modeling of nonlinear wave by the two FD schemes presented in the paper can provide us with an easy method to keep the stability and convergence in the calculation field when the product of the absolute value of nonlinear coefficient and the value of u/x are small enough, namely, the value of bu/x is much smaller than 1. Several results are founded in the numerical study of nonlinear P-wave propagation, such as the waveform aberration, the generation and growth of harmonic wave and the energy redistribution among different frequency components. All of them will be more violent when the initial amplitude A0 is larger or the nonlinearity of medium is stronger. Correspondingly, we have found that the nonlinear P-wave propagation velocity will change with different initial frequency f of source wave or the wave velocity c (equal to the P-wave velocity in the same medium without considering nonlinearity).展开更多
In order to obtain satisfactory mechanical properties for the cam used in high-power ship diesel engines, a new quenching technology was proposed by designing a two-stage quenching process with an alkaline bath as the...In order to obtain satisfactory mechanical properties for the cam used in high-power ship diesel engines, a new quenching technology was proposed by designing a two-stage quenching process with an alkaline bath as the quenching medium. To demonstrate the effectiveness of the proposed new quenching technology, both numerical analysis and experimental study were performed. The new quenching technology was analyzed using finite element method. The combined effects of the temperature, stress and microstructure fields were investigated considering nonlinear material properties. Finally, an experimental study was performed to verify the effectiveness of the proposed new quenching technology. The numerical results show that internal stress is affected by both thermal stress and transformation stress. In addition, the direction of the internal stress is changed several times due to thermal interaction and microstructure evolution during the quenching process. The experimental results show that the proposed new quenching technology significantly improves the mechanical properties and microstructures of the cam. The tensile strength, the impact resistance and the hardness value of the cam by the proposed new quenching technology are improved by 4.3%, 8.9% and 3.5% compared with those by the traditional quenching technology. Moreover, the residual stress and cam shape deformation are reduced by 40.0% and 48.9% respectively for the cam manufactured by the new quenching technology.展开更多
In order to investigate the temperature distribution, a three-dimensional finite element model (FEM) was developed to simulate the temperature regime in the channels of double-loop inductor, and the simulated result...In order to investigate the temperature distribution, a three-dimensional finite element model (FEM) was developed to simulate the temperature regime in the channels of double-loop inductor, and the simulated results were compared with experimental data from low load trials of a 400 kW inductor. The results of numerical simulations, such as the temperature and Joule heating rate, show reasonable correlation with experimental data. The results indicate that Joule heating rate and the temperature reach the maximum at the comers and the minimum at the centre of the cross-section area. The temperature difference between the inlet and outlet is in an inverse proportion to mass transport. Joule heating rate and the temperature are directly proportional to power frequency. It is concluded that mass transport and power frequency play a critical role in determining the temperature regime and Joule heating rate, the relative permeability of the magnetic core shows no significant influence on temperature regime and Joule heating rate, when the relative permeability varies from 5 000 to 10 000.展开更多
地震勘探中的波动方程正演模拟受计算和存储能力的限制,只能在有限空间进行,需要对其设置边界条件。常用的完全匹配层(Perfect Match Layer,PML)是一种应用广泛的边界条件,需要对边界条件设定一定的层数,层数太大时会降低正演速度和增...地震勘探中的波动方程正演模拟受计算和存储能力的限制,只能在有限空间进行,需要对其设置边界条件。常用的完全匹配层(Perfect Match Layer,PML)是一种应用广泛的边界条件,需要对边界条件设定一定的层数,层数太大时会降低正演速度和增大内存存储。因此,文中通过将PML边界条件和改进后的梯度黏弹边界(Improved Gradient Viscoelastic Layer,IGVL)结合提出一种PML-IGVL边界条件,同时,为了减少截断误差和数值频散,将PML-IGVL边界条件应用于紧致交错网格的黏滞声波方程中进行数值模拟。均匀介质和Marmousi模型中的波场数值模拟结果表明,相对于PML边界条件,PML-IGVL边界条件使用的边界厚度更小,在相同较低边界层数下,能更好地吸收反射波,同时能节省更多的存储空间,并提高正演效率,证明了PML-IGVL边界条件的有效性和优越性,是正演模拟中一种有效的边界吸收方法。展开更多
文摘The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of curl interface underground or “v” figure complex model, it is difficult to realize. So it is important to forward the complex geoelectricity model. This paper takes two Maxwell’s vorticity equations as departure point, makes use of the principles of Yee’s space grid model theory and the basic principle finite difference time domain method, and deduces a GPR forward system of equation of two dimensional spaces. The Mur super absorbed boundary condition is adopted to solve the super strong reflection on the interceptive boundary when there is the forward simulation. And a self-made program is used to process forward simulation to two typical geoelectricity model.
文摘A new general optimal principle of designing explicit finite difference method was obtained. Several applied cases were put forward to explain the uses of the principle. The validity of the principal was tested by a numeric example.
基金Project(50175034) supported by the National Natural Science Foundation of China
文摘In order to investigate the process of incremental sheet forming (ISF) through both experimental and numerical approaches, a three-dimensional elasto-plastic finite element model (FEM) was developed to simulate the process and the simulated results were compared with those of experiment. The results of numerical simulations, such as the strain history and distribution, the stress state and distribution, sheet thickness distribution, etc, were discussed in details, and the influences of process parameters on these results were also analyzed. The simulated results of the radial strain and the thickness distribution are in good agreement with experimental results. The simulations reveal that the deformation is localized around the tool and constantly remains close to a plane strain state. With decreasing depth step, increasing tool diameter and wall inclination angle, the axial stress reduces, leading to less thinning and more homogeneous plastic strain and thickness distribution. During ISF, the plastic strain increases stepwise under the action of the tool. Each increase in plastic strain is accompanied by hydrostatic pressure, which explains why obtainable deformation using ISF exceeds the forming limits of conventional sheet forming.
文摘Being aimed at the inside wall wrinkling and sinking phenomenon of palladium-yttrium alloy thin wall spiral tube used for preparation of high purity hydrogen, extraction of hydrogen isotope, and purification and separation of hydrogen in the winding process, this article analyzed the reasons for above phenomena, established a numerical simulation model of winding process of above tube, using elastic-plastic Finite Element method analyzed the max. tensile stress and max. compression stress and their locations, thereby provides a theory base for the control of working forming course of thin wall spiral tube.
文摘In this paper, the internal fluid motion of a jet system is described by the Navier Stokes mechanics equations. For the simulation of the motion, the penalty function finite element method is used, and the velocity vectors and stream function curves are obtained. Using the Prandtl theory, this paper derives the free jet velocity and the jet bunch width in a half-space, the latter of which is amended by experiment. The results obtained in this paper are applied to micro-type high pressure water jet cleaner and the ejector of rocket engine.
基金Project of Knowledge Innovation Program from Chinese Academy of Sciences (KZCX2-109).
文摘Two central schemes of finite difference (FD) up to different accuracy orders of space sampling step Dx (Fourth order and Sixth order respectively) were used to study the 1-D nonlinear P-wave propagation in the nonlinear solid media by the numerical method. Distinctly different from the case of numerical modeling of linear elastic wave, there may be several difficulties in the numerical treatment to the nonlinear partial differential equation, such as the steep gradients, shocks and unphysical oscillations. All of them are the great obstacles to the stability and conver-gence of numerical calculation. Fortunately, the comparative study on the modeling of nonlinear wave by the two FD schemes presented in the paper can provide us with an easy method to keep the stability and convergence in the calculation field when the product of the absolute value of nonlinear coefficient and the value of u/x are small enough, namely, the value of bu/x is much smaller than 1. Several results are founded in the numerical study of nonlinear P-wave propagation, such as the waveform aberration, the generation and growth of harmonic wave and the energy redistribution among different frequency components. All of them will be more violent when the initial amplitude A0 is larger or the nonlinearity of medium is stronger. Correspondingly, we have found that the nonlinear P-wave propagation velocity will change with different initial frequency f of source wave or the wave velocity c (equal to the P-wave velocity in the same medium without considering nonlinearity).
基金Project(50875268) supported by the National Natural Science Foundation of China Project(CSTC2008AB3057) supported by Foundation of Chongqing Science and Technology Commission, China+1 种基金 Project(108107) supported by the Key Project of Ministry of Education of China Project(50925518) supported by the National Science Fund for Distinguished Young Scholars
文摘In order to obtain satisfactory mechanical properties for the cam used in high-power ship diesel engines, a new quenching technology was proposed by designing a two-stage quenching process with an alkaline bath as the quenching medium. To demonstrate the effectiveness of the proposed new quenching technology, both numerical analysis and experimental study were performed. The new quenching technology was analyzed using finite element method. The combined effects of the temperature, stress and microstructure fields were investigated considering nonlinear material properties. Finally, an experimental study was performed to verify the effectiveness of the proposed new quenching technology. The numerical results show that internal stress is affected by both thermal stress and transformation stress. In addition, the direction of the internal stress is changed several times due to thermal interaction and microstructure evolution during the quenching process. The experimental results show that the proposed new quenching technology significantly improves the mechanical properties and microstructures of the cam. The tensile strength, the impact resistance and the hardness value of the cam by the proposed new quenching technology are improved by 4.3%, 8.9% and 3.5% compared with those by the traditional quenching technology. Moreover, the residual stress and cam shape deformation are reduced by 40.0% and 48.9% respectively for the cam manufactured by the new quenching technology.
基金Project(50876116) supported by the National Natural Science Foundation of ChinaProject(2007CK3077) supported by the Innovative Program of Hunan Science and Technology Agency, ChinaProject(1343-77225) supported by the Graduate School of Central South University, China
文摘In order to investigate the temperature distribution, a three-dimensional finite element model (FEM) was developed to simulate the temperature regime in the channels of double-loop inductor, and the simulated results were compared with experimental data from low load trials of a 400 kW inductor. The results of numerical simulations, such as the temperature and Joule heating rate, show reasonable correlation with experimental data. The results indicate that Joule heating rate and the temperature reach the maximum at the comers and the minimum at the centre of the cross-section area. The temperature difference between the inlet and outlet is in an inverse proportion to mass transport. Joule heating rate and the temperature are directly proportional to power frequency. It is concluded that mass transport and power frequency play a critical role in determining the temperature regime and Joule heating rate, the relative permeability of the magnetic core shows no significant influence on temperature regime and Joule heating rate, when the relative permeability varies from 5 000 to 10 000.