Numerical simulation is the most powerful computational and analysis tool for a large variety of engineering and physical problems.For a complex problem relating to multi-field,multi-process and multi-scale,different ...Numerical simulation is the most powerful computational and analysis tool for a large variety of engineering and physical problems.For a complex problem relating to multi-field,multi-process and multi-scale,different computing tools have to be developed so as to solve particular fields at different scales and for different processes.Therefore,the integration of different types of software is inevitable.However,it is difficult to perform the transfer of the meshes and simulated results among software packages because of the lack of shared data formats or encrypted data formats.An image processing based method for three-dimensional model reconstruction for numerical simulation was proposed,which presents a solution to the integration problem by a series of slice or projection images obtained by the post-processing modules of the numerical simulation software.By means of mapping image pixels to meshes of either finite difference or finite element models,the geometry contour can be extracted to export the stereolithography model.The values of results,represented by color,can be deduced and assigned to the meshes.All the models with data can be directly or indirectly integrated into other software as a continued or new numerical simulation.The three-dimensional reconstruction method has been validated in numerical simulation of castings and case studies were provided in this study.展开更多
Objective With the cranial aneurysm’s 3D-computational hemodynamics numerical simulation,we can get the cranial aneurysms and its ’s mother arteries ’blood flow rate,path and wall shearing stress and so on. Taking ...Objective With the cranial aneurysm’s 3D-computational hemodynamics numerical simulation,we can get the cranial aneurysms and its ’s mother arteries ’blood flow rate,path and wall shearing stress and so on. Taking use of the contrast analysis of the hemodynamics result and the findings in operation,we try to find hemodynamic character and mechanism of rupture展开更多
Based on the 2D horizontal plane numerical model,a quasi-3D numerical model is establishedfor coastal regions of shallow water.The characteristics of this model are that the velocity profiles can be ob-tained at the s...Based on the 2D horizontal plane numerical model,a quasi-3D numerical model is establishedfor coastal regions of shallow water.The characteristics of this model are that the velocity profiles can be ob-tained at the same time when the equations of the value of difference between the horizontal current velocityand its depth-averaged velocity in the vertical direction are solved and the results obtained are consistent withthe results of the 2D model.The circulating flow in the rectangular area induced by wind is simulated and ap-plied to the tidal flow field of the radial sandbanks in the South Yellow Sea.The computational results fromthis quasi-3D model are in good agreement with analytical results and observed data.The solution of the finitedifference equations has been found to be stable,and the model is simple,effective and practical.展开更多
Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of...Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.展开更多
A 3- D free surface flow in open channels based on the Reynolds equations with the k-εturbulence closure model is presented in this paper. Insted of the "rigid lid" approximation, the solutionof the free su...A 3- D free surface flow in open channels based on the Reynolds equations with the k-εturbulence closure model is presented in this paper. Insted of the "rigid lid" approximation, the solutionof the free surface equation is implemented in the velocity-pressure iterative procedure on the basis ofthe conventional SIMPLE method. This model was used to compute the flow in rectangular channelswith trenches dredged across the bottom. The velocity, eddy viscosity coefficient, turbulent shear stress,turbulent kinetic energy and elevation of the free surface can be obtained. The computed results are ingood agreement with previous experimental data.展开更多
In this study,the compressive,split tensile,and flexural strengths of concrete with nano-CaCO_(3) only were compared with those of concrete with nano-CaCO_(3) and basalt fibers through field experiments,and the underl...In this study,the compressive,split tensile,and flexural strengths of concrete with nano-CaCO_(3) only were compared with those of concrete with nano-CaCO_(3) and basalt fibers through field experiments,and the underlying mechanisms were analyzed by the Scanning Electron Microscope (SEM) techniques.On the mesoscale,a damage model of concrete was established based on the continuum progressive damage theory,which was used to investigate the influence of different lengths and contents of fibers on the mechanical properties of concrete.Then,the experimental and numerical simulation results were compared and analyzed to verify the feasibility of model.The results show that nano-CaCO_(3) can enhance the compressive strength of the concrete,with an optimal content of 2.0%.Adding basalt fibers into the nano-CaCO_(3) reinforced concrete may further enhance the compressive,split tensile,and flexural strengths of the concrete;however,the higher content of basalt fiber can not lead to higher performance of concrete.The optimal length and content of fiber are 6 mm and 0.20%,respectively.The SEM result shows that the aggregation of basalt fibers is detrimental to the mechanical properties of concrete.The numerical simulation results are in good agreement with the experimental results.展开更多
A 3-D finite-element numerical simulation model of temperature field for CIESC casting solidification process was developed with the aid of ANSYS software and a series ofcorresponding experiments were made. The result...A 3-D finite-element numerical simulation model of temperature field for CIESC casting solidification process was developed with the aid of ANSYS software and a series ofcorresponding experiments were made. The results showed that the good agreement was obtained between the numerical simulation and the experiments. Based on the numerical simulation results, the characteristics of tem- perature distribution in the castings during CIESC solidification process were analyzed and summarized. According to the G/R1/2 method and numerical simulation results, there is no any shrinkage defect in the CIESE casting and structure of casting is fine and compact.展开更多
The weak structure plane is an important factor affecting the stability of rock slope, and detecting the spatial structure of the weak structural plane is beneficial to analyze the stability of the slope and estimate ...The weak structure plane is an important factor affecting the stability of rock slope, and detecting the spatial structure of the weak structural plane is beneficial to analyze the stability of the slope and estimate the quantity of the landslide. Based on 3D electrical resistivity tomography, a model of rock slope with weak structure plane is established, and the characteristics of three-dimensional resistivity imaging of weak structure plane under different ground water conditions are simulated. The results show that the weak structural plane has a better reflection in 3D electrical resistivity tomography;the distribution of weak structural plane of 3D resistivity imaging can be roughly determined under different ground water conditions;the three-dimensional electrical resistivity tomography is feasible in the detection of weak structural plane of rock slope.展开更多
POM97, an oceanic model, has been used for the first time to the numerical study on the tidal waves of the as Regions around Taiwan. In this paper, we have got the result that the semidiurnal tidal waves of these area...POM97, an oceanic model, has been used for the first time to the numerical study on the tidal waves of the as Regions around Taiwan. In this paper, we have got the result that the semidiurnal tidal waves of these area mainly are the co-operating tides which come from the south of 23’N of the western Pacific. Those semidiurnal tidal waves affecting the Taiwan Straits come respectively from the south and the north entrance of the channel, and the north tidal wave is stronger than the south one. The strongest tidal field is the area from the Meizhou Bay to the Xinhua Bay along the coast of Fujian Province, where the biggest amplitude of the M2 partial tide can reach 240 cm. The strongest tidal cur- rent fields lie in the Penghu watercourse, where the maximum velocity of the M2 partial tide can arrive at 196 m/s. In the horizontal structure of the tidal currets, we have found that there is a stream dot in the north of the channel, besides, there still exist four new ones. As for the vertical structure, it mainly is biassed to the right at the surface, and to the left near the bottom layer.展开更多
—The numerical simulation is based on the authors'high-order models with a dissipative termfor nonlinear and dispersive wave in water of varying depth.Corresponding finite-difference equations andgeneral conditions ...—The numerical simulation is based on the authors'high-order models with a dissipative termfor nonlinear and dispersive wave in water of varying depth.Corresponding finite-difference equations andgeneral conditions for open and fixed natural boundaries with an arbitrary reflection coefficient and phaseshift are also given in this paper.The systematical tests of numerical simulation show that the theoreticalmodels,the finite-difference algorithms and the boundary conditions can give good calculation results forthe wave propagating in shallow and deep water with an arbitrary slope varying from gentle to steep.展开更多
Design of forming dies and whole process of simulation of cold rolling involutes spline can be realized by using of CAD software of PRO-E and CAE software of DEFORM-3D. Software DEFORM-3D provides an automatic and opt...Design of forming dies and whole process of simulation of cold rolling involutes spline can be realized by using of CAD software of PRO-E and CAE software of DEFORM-3D. Software DEFORM-3D provides an automatic and optimized remeshing function, especially for the large deformation. In order to use this function sufficiently, simulation of cold rolling involutes spline can be implemented indirectly. The relationship between die and workpiece, forming force and characteristic of deformation in the forming process of cold rolling involutes spline are analyzed and researched. Meanwhile, reliable proofs for the design of dies and deforming equipment are provided.展开更多
Based on specific geology and mining conditions of certain coal working face in China, a series of numerical models under different unconsolidated layers thickness were respectively established by employing FLAC3D. Th...Based on specific geology and mining conditions of certain coal working face in China, a series of numerical models under different unconsolidated layers thickness were respectively established by employing FLAC3D. The relationship between the unconsolidated layers thickness and surface movement laws was studied. Maximum surface subsidence, Maximum horizontal displacement and surface subsidence degree were obtained. Contours of surface subsidence/horizontal displacement and curves were drawn. Some laws of surface subsidence/horizontal displacement were analyzed. The role of the unconsolidated layers in surface subsidence was revealed. It is significant to predict surface subsidence of thick unconsolidated layers for coal mine and take effective measures to control surface subsidence.展开更多
The numerical simulation of internal features,such as inclusions and voids,is important to analyze their impact on the performance of composite materials.However,the complex geometries of internal features and the ind...The numerical simulation of internal features,such as inclusions and voids,is important to analyze their impact on the performance of composite materials.However,the complex geometries of internal features and the induced continuous-discontinuous(C-D)deformation fields are challenges to their numerical simulation.In this study,a 3D approach using a simple mesh to simulate irregular internal geometries is developed for the first time.With the help of a developed voxel crack model,image models that are efficient when recording complex geometries are directly imported into the simulation.Surface reconstructions,which are usually labor-intensive,are excluded from this approach.Moreover,using image models as the geometric input,image processing techniques are applied to detect material interfaces and develop contact pairs.Then,the C-D deformations of the complex internal features are directly calculated based on the numerical manifold method.The accuracy and convergence of the developed3D approach are examined based on multiple benchmarks.Successful 3D C-D simulation of sandstones with naturally formed complex microfeatures demonstrates the capability of the developed approach.展开更多
Flat subduction can significantly influence the distribution of volcanism,stress state,and surface topography of the overriding plate.However,the mechanisms for inducing flat subduction remain controversial.Previous t...Flat subduction can significantly influence the distribution of volcanism,stress state,and surface topography of the overriding plate.However,the mechanisms for inducing flat subduction remain controversial.Previous two-dimensional(2-D)numerical models and laboratory analogue models suggested that a buoyant impactor(aseismic ridge,oceanic plateau,or the like)may induce flat subduction.However,three-dimensional(3-D)systematic studies on the relationship between flat subduction and buoyant blocks are still lacking.Here,we use a 3-D numerical model to investigate the influence of the aseismic ridge,especially its width(which is difficult to consider in 2-D numerical models),on the formation of flat subduction.Our model results suggest that the aseismic ridge needs to be wide and thick enough to induce flat subduction,a condition that is difficult to satisfy on the Earth.We also find that the subduction of an aseismic ridge parallel to the trench or a double aseismic ridge normal to the trench has a similar effect on super-wide aseismic ridge subduction in terms of causing flat subduction,which can explain the flat subduction observed beneath regions such as Chile and Peru.展开更多
Forward modeling is the basis of inversion imaging and quantitative interpretation for DC resistivity exploration.Currently,a numerical model of the DC resistivity method must be finely divided to obtain a highly accu...Forward modeling is the basis of inversion imaging and quantitative interpretation for DC resistivity exploration.Currently,a numerical model of the DC resistivity method must be finely divided to obtain a highly accurate solution under complex conditions,resulting in a long calculation time and large storage.Therefore,we propose a 3D numerical simulation method in a mixed space-wavenumber domain to overcome this challenge.The partial differential equation about abnormal potential is transformed into many independent ordinary differential equations with different wavenumbers using a 2D Fourier transform along the x axis and y axis direction.In this way,a large-scale 3D numerical simulation problem is decomposed into several 1D numerical simulation problems,which significantly reduces the computational and storage requirements.In addition,these ordinary 1D differential equations with different wavenumbers are independent of each other and high parallelelism of the algorithm.They are solved using a finite-element algorithm combined with a chasing method,and the obtained solution is modified using a contraction operator.In this method,the vertical direction is reserved as the spatial domain,then grid size can be determined flexibly based on the underground current density distribution,which considers the solution accuracy and calculation efficiency.In addition,for the first time,we use the contraction operator in the integral equation method to iterate the algorithm.The algorithm takes advantage of the high efficiency of the standard Fourier transform and chasing method,as well as the fast convergence of the contraction operator.We verified the accuracy of the algorithm and the convergence of the contraction operator.Compared with a volume integral method and goal-oriented adaptive finite-element method,the proposed algorithm has lower memory requirements and high computational efficiency,making it suitable for calculating a model with large-scale nodes.Moreover,different examples are used to verify the high adaptability and parallelism of the proposed algorithm.The findings show that the 3D numerical simulation method of DC resistivity method in a mixed space-wavenumber domain is highly efficient,precise,and parallel.展开更多
Basing on the analysis of the traits of the roll forging process, a system-model of computer simulation has been established. Three-dimensional rigid-plastic FEM has been used for the simulation of the deformation pro...Basing on the analysis of the traits of the roll forging process, a system-model of computer simulation has been established. Three-dimensional rigid-plastic FEM has been used for the simulation of the deformation process in the oval and round pass rolling, including the entering, rolling, and separating stages. The analysis was conducted using the Deform-3D ver.5.0 code. The important information concerned with the deformation area characteristic, material flow, and velocity field has been presented. Otherwise,the location of the neutral plane in the deformation area was shown clearly.展开更多
A modified piano key weir with a rounded nose and a parapet wall (MPKW) can improve the discharge capacity significantly compared to a standard piano key weir. However, the optimum of the inlet/outlet width ratio (Wi/...A modified piano key weir with a rounded nose and a parapet wall (MPKW) can improve the discharge capacity significantly compared to a standard piano key weir. However, the optimum of the inlet/outlet width ratio (Wi/Wo) on the discharge efficiency of MPKW is still not investigated numerically. The present work utilized the numerical modeling to investigate and analyze the effects of the inlet/outlet key width ratios on the hydraulic characteristics and discharge capacity of the MPKW. To validate the numerical model with the experimental data, the results indicate that the average relative error is 2.96%, which confirms that the numerical model is fairly well to predictthe specifications of flow over on the MPKW. Numerical simulation results indicated that the discharge capacity of the MPKW can be improved up to 8.5% by optimizing the Wi/Wo ratio ranging from 1.53 to 1.67 even if the other parameters of the MPKW keep unchanged. A big Wi/Wo ratio generally leads to an increase in discharge capacity at low heads and a little effect on the discharge efficiency at high heads. The discharge efficiency of the inlet and outlet crests increases up to 9.6% for high heads, while discharge efficiency of the lateral crest decreases up to 23.5% compared with the reference model. The findings of the study revealed that the intrinsic influencing mechanism of the Wi/Wo ratio on the discharge performance of MPKWs.展开更多
In meandering rivers, the flow pattern is highly complex, with specific characteristics at bends that are not observed along straight paths. A numerical model can be effectively used to predict such flow fields. Since...In meandering rivers, the flow pattern is highly complex, with specific characteristics at bends that are not observed along straight paths. A numerical model can be effectively used to predict such flow fields. Since river bends are not uniform-some are divergent and others convergent-in this study, after the SSIIM 3-D model was calibrated using the result of measurements along a uniform 180° bend with a width of 0.6 m, a similar but convergent 180° bend, 0.6 m to 0.45 m wide, was simulated using the SSIIM 3-D numerical model. Flow characteristics of the convergent 180° bend, including lengthwise and vertical velocity profiles, primary and secondary flows, lengthwise and widthwise slopes of the water surface, and the helical flow strength, were compared with those of the uniform 180° bend. The verification results of the model show that the numerical model can effectively simulate the flow field in the uniform bend. In addition, this research indicates that, in a convergent channel, the maximum velocity path at a plane near the water surface crosses the channel's centerline at about a 30° to 40° cross-section, while in the uniform bend, this occurs at about the 50° cross-section. The varying range of the water surface elevation is wider in the convergent channel than in the uniform one, and the strength of the helical flow is generally greater in the uniform channel than in the convergent one. Also, unlike the uniform bend, the convergent bend exhibits no rotational cell against the main direction of secondary flow rotation at the 135° cross-section.展开更多
基金funded by National Key R&D Program of China(No.2021YFB3401200)the National Natural Science Foundation of China(No.51875308)the Beijing Nature Sciences Fund-Haidian Originality Cooperation Project(L212002).
文摘Numerical simulation is the most powerful computational and analysis tool for a large variety of engineering and physical problems.For a complex problem relating to multi-field,multi-process and multi-scale,different computing tools have to be developed so as to solve particular fields at different scales and for different processes.Therefore,the integration of different types of software is inevitable.However,it is difficult to perform the transfer of the meshes and simulated results among software packages because of the lack of shared data formats or encrypted data formats.An image processing based method for three-dimensional model reconstruction for numerical simulation was proposed,which presents a solution to the integration problem by a series of slice or projection images obtained by the post-processing modules of the numerical simulation software.By means of mapping image pixels to meshes of either finite difference or finite element models,the geometry contour can be extracted to export the stereolithography model.The values of results,represented by color,can be deduced and assigned to the meshes.All the models with data can be directly or indirectly integrated into other software as a continued or new numerical simulation.The three-dimensional reconstruction method has been validated in numerical simulation of castings and case studies were provided in this study.
文摘Objective With the cranial aneurysm’s 3D-computational hemodynamics numerical simulation,we can get the cranial aneurysms and its ’s mother arteries ’blood flow rate,path and wall shearing stress and so on. Taking use of the contrast analysis of the hemodynamics result and the findings in operation,we try to find hemodynamic character and mechanism of rupture
基金National Natural Science Foundation of China(Grant No.49236120)
文摘Based on the 2D horizontal plane numerical model,a quasi-3D numerical model is establishedfor coastal regions of shallow water.The characteristics of this model are that the velocity profiles can be ob-tained at the same time when the equations of the value of difference between the horizontal current velocityand its depth-averaged velocity in the vertical direction are solved and the results obtained are consistent withthe results of the 2D model.The circulating flow in the rectangular area induced by wind is simulated and ap-plied to the tidal flow field of the radial sandbanks in the South Yellow Sea.The computational results fromthis quasi-3D model are in good agreement with analytical results and observed data.The solution of the finitedifference equations has been found to be stable,and the model is simple,effective and practical.
基金Project(51375498)supported by the National Natural Science Foundation of China
文摘Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.
文摘A 3- D free surface flow in open channels based on the Reynolds equations with the k-εturbulence closure model is presented in this paper. Insted of the "rigid lid" approximation, the solutionof the free surface equation is implemented in the velocity-pressure iterative procedure on the basis ofthe conventional SIMPLE method. This model was used to compute the flow in rectangular channelswith trenches dredged across the bottom. The velocity, eddy viscosity coefficient, turbulent shear stress,turbulent kinetic energy and elevation of the free surface can be obtained. The computed results are ingood agreement with previous experimental data.
基金Funded by the National Natural Science Foundation of China (No. 51969026)the Natural Science Foundation of Qinghai Province in China (No. 2018-ZJ-750)。
文摘In this study,the compressive,split tensile,and flexural strengths of concrete with nano-CaCO_(3) only were compared with those of concrete with nano-CaCO_(3) and basalt fibers through field experiments,and the underlying mechanisms were analyzed by the Scanning Electron Microscope (SEM) techniques.On the mesoscale,a damage model of concrete was established based on the continuum progressive damage theory,which was used to investigate the influence of different lengths and contents of fibers on the mechanical properties of concrete.Then,the experimental and numerical simulation results were compared and analyzed to verify the feasibility of model.The results show that nano-CaCO_(3) can enhance the compressive strength of the concrete,with an optimal content of 2.0%.Adding basalt fibers into the nano-CaCO_(3) reinforced concrete may further enhance the compressive,split tensile,and flexural strengths of the concrete;however,the higher content of basalt fiber can not lead to higher performance of concrete.The optimal length and content of fiber are 6 mm and 0.20%,respectively.The SEM result shows that the aggregation of basalt fibers is detrimental to the mechanical properties of concrete.The numerical simulation results are in good agreement with the experimental results.
文摘A 3-D finite-element numerical simulation model of temperature field for CIESC casting solidification process was developed with the aid of ANSYS software and a series ofcorresponding experiments were made. The results showed that the good agreement was obtained between the numerical simulation and the experiments. Based on the numerical simulation results, the characteristics of tem- perature distribution in the castings during CIESC solidification process were analyzed and summarized. According to the G/R1/2 method and numerical simulation results, there is no any shrinkage defect in the CIESE casting and structure of casting is fine and compact.
文摘The weak structure plane is an important factor affecting the stability of rock slope, and detecting the spatial structure of the weak structural plane is beneficial to analyze the stability of the slope and estimate the quantity of the landslide. Based on 3D electrical resistivity tomography, a model of rock slope with weak structure plane is established, and the characteristics of three-dimensional resistivity imaging of weak structure plane under different ground water conditions are simulated. The results show that the weak structural plane has a better reflection in 3D electrical resistivity tomography;the distribution of weak structural plane of 3D resistivity imaging can be roughly determined under different ground water conditions;the three-dimensional electrical resistivity tomography is feasible in the detection of weak structural plane of rock slope.
文摘POM97, an oceanic model, has been used for the first time to the numerical study on the tidal waves of the as Regions around Taiwan. In this paper, we have got the result that the semidiurnal tidal waves of these area mainly are the co-operating tides which come from the south of 23’N of the western Pacific. Those semidiurnal tidal waves affecting the Taiwan Straits come respectively from the south and the north entrance of the channel, and the north tidal wave is stronger than the south one. The strongest tidal field is the area from the Meizhou Bay to the Xinhua Bay along the coast of Fujian Province, where the biggest amplitude of the M2 partial tide can reach 240 cm. The strongest tidal cur- rent fields lie in the Penghu watercourse, where the maximum velocity of the M2 partial tide can arrive at 196 m/s. In the horizontal structure of the tidal currets, we have found that there is a stream dot in the north of the channel, besides, there still exist four new ones. As for the vertical structure, it mainly is biassed to the right at the surface, and to the left near the bottom layer.
文摘—The numerical simulation is based on the authors'high-order models with a dissipative termfor nonlinear and dispersive wave in water of varying depth.Corresponding finite-difference equations andgeneral conditions for open and fixed natural boundaries with an arbitrary reflection coefficient and phaseshift are also given in this paper.The systematical tests of numerical simulation show that the theoreticalmodels,the finite-difference algorithms and the boundary conditions can give good calculation results forthe wave propagating in shallow and deep water with an arbitrary slope varying from gentle to steep.
文摘Design of forming dies and whole process of simulation of cold rolling involutes spline can be realized by using of CAD software of PRO-E and CAE software of DEFORM-3D. Software DEFORM-3D provides an automatic and optimized remeshing function, especially for the large deformation. In order to use this function sufficiently, simulation of cold rolling involutes spline can be implemented indirectly. The relationship between die and workpiece, forming force and characteristic of deformation in the forming process of cold rolling involutes spline are analyzed and researched. Meanwhile, reliable proofs for the design of dies and deforming equipment are provided.
基金Project(2007BAK28B03)supported by the National Eleventh-Five Year Research Program of ChinaProject(2010YD05)supported by the Fundamental Research Funds for the Central UniversitiesProject(200911036)supported by the Ministry of Land and Resources Research Special
文摘Based on specific geology and mining conditions of certain coal working face in China, a series of numerical models under different unconsolidated layers thickness were respectively established by employing FLAC3D. The relationship between the unconsolidated layers thickness and surface movement laws was studied. Maximum surface subsidence, Maximum horizontal displacement and surface subsidence degree were obtained. Contours of surface subsidence/horizontal displacement and curves were drawn. Some laws of surface subsidence/horizontal displacement were analyzed. The role of the unconsolidated layers in surface subsidence was revealed. It is significant to predict surface subsidence of thick unconsolidated layers for coal mine and take effective measures to control surface subsidence.
基金supported by the National Natural Science Foundation of China(Grant Nos.41807277,42172306,and U1965204)the Natural Science Foundation of Hebei Province(Grant No.D2019202440)。
文摘The numerical simulation of internal features,such as inclusions and voids,is important to analyze their impact on the performance of composite materials.However,the complex geometries of internal features and the induced continuous-discontinuous(C-D)deformation fields are challenges to their numerical simulation.In this study,a 3D approach using a simple mesh to simulate irregular internal geometries is developed for the first time.With the help of a developed voxel crack model,image models that are efficient when recording complex geometries are directly imported into the simulation.Surface reconstructions,which are usually labor-intensive,are excluded from this approach.Moreover,using image models as the geometric input,image processing techniques are applied to detect material interfaces and develop contact pairs.Then,the C-D deformations of the complex internal features are directly calculated based on the numerical manifold method.The accuracy and convergence of the developed3D approach are examined based on multiple benchmarks.Successful 3D C-D simulation of sandstones with naturally formed complex microfeatures demonstrates the capability of the developed approach.
基金supported by the Strategic Priority Research Program of the Chinese Academy of Sciences(Grant No.XDB 41000000)the National Natural Science Foundation of China(Grant No.41820104004)the Fundamental Research Funds for the Central Universities(Grant No.WK2080000144).
文摘Flat subduction can significantly influence the distribution of volcanism,stress state,and surface topography of the overriding plate.However,the mechanisms for inducing flat subduction remain controversial.Previous two-dimensional(2-D)numerical models and laboratory analogue models suggested that a buoyant impactor(aseismic ridge,oceanic plateau,or the like)may induce flat subduction.However,three-dimensional(3-D)systematic studies on the relationship between flat subduction and buoyant blocks are still lacking.Here,we use a 3-D numerical model to investigate the influence of the aseismic ridge,especially its width(which is difficult to consider in 2-D numerical models),on the formation of flat subduction.Our model results suggest that the aseismic ridge needs to be wide and thick enough to induce flat subduction,a condition that is difficult to satisfy on the Earth.We also find that the subduction of an aseismic ridge parallel to the trench or a double aseismic ridge normal to the trench has a similar effect on super-wide aseismic ridge subduction in terms of causing flat subduction,which can explain the flat subduction observed beneath regions such as Chile and Peru.
文摘Forward modeling is the basis of inversion imaging and quantitative interpretation for DC resistivity exploration.Currently,a numerical model of the DC resistivity method must be finely divided to obtain a highly accurate solution under complex conditions,resulting in a long calculation time and large storage.Therefore,we propose a 3D numerical simulation method in a mixed space-wavenumber domain to overcome this challenge.The partial differential equation about abnormal potential is transformed into many independent ordinary differential equations with different wavenumbers using a 2D Fourier transform along the x axis and y axis direction.In this way,a large-scale 3D numerical simulation problem is decomposed into several 1D numerical simulation problems,which significantly reduces the computational and storage requirements.In addition,these ordinary 1D differential equations with different wavenumbers are independent of each other and high parallelelism of the algorithm.They are solved using a finite-element algorithm combined with a chasing method,and the obtained solution is modified using a contraction operator.In this method,the vertical direction is reserved as the spatial domain,then grid size can be determined flexibly based on the underground current density distribution,which considers the solution accuracy and calculation efficiency.In addition,for the first time,we use the contraction operator in the integral equation method to iterate the algorithm.The algorithm takes advantage of the high efficiency of the standard Fourier transform and chasing method,as well as the fast convergence of the contraction operator.We verified the accuracy of the algorithm and the convergence of the contraction operator.Compared with a volume integral method and goal-oriented adaptive finite-element method,the proposed algorithm has lower memory requirements and high computational efficiency,making it suitable for calculating a model with large-scale nodes.Moreover,different examples are used to verify the high adaptability and parallelism of the proposed algorithm.The findings show that the 3D numerical simulation method of DC resistivity method in a mixed space-wavenumber domain is highly efficient,precise,and parallel.
基金supported by the National Natural Science Foundation of China(No.50675014).
文摘Basing on the analysis of the traits of the roll forging process, a system-model of computer simulation has been established. Three-dimensional rigid-plastic FEM has been used for the simulation of the deformation process in the oval and round pass rolling, including the entering, rolling, and separating stages. The analysis was conducted using the Deform-3D ver.5.0 code. The important information concerned with the deformation area characteristic, material flow, and velocity field has been presented. Otherwise,the location of the neutral plane in the deformation area was shown clearly.
文摘A modified piano key weir with a rounded nose and a parapet wall (MPKW) can improve the discharge capacity significantly compared to a standard piano key weir. However, the optimum of the inlet/outlet width ratio (Wi/Wo) on the discharge efficiency of MPKW is still not investigated numerically. The present work utilized the numerical modeling to investigate and analyze the effects of the inlet/outlet key width ratios on the hydraulic characteristics and discharge capacity of the MPKW. To validate the numerical model with the experimental data, the results indicate that the average relative error is 2.96%, which confirms that the numerical model is fairly well to predictthe specifications of flow over on the MPKW. Numerical simulation results indicated that the discharge capacity of the MPKW can be improved up to 8.5% by optimizing the Wi/Wo ratio ranging from 1.53 to 1.67 even if the other parameters of the MPKW keep unchanged. A big Wi/Wo ratio generally leads to an increase in discharge capacity at low heads and a little effect on the discharge efficiency at high heads. The discharge efficiency of the inlet and outlet crests increases up to 9.6% for high heads, while discharge efficiency of the lateral crest decreases up to 23.5% compared with the reference model. The findings of the study revealed that the intrinsic influencing mechanism of the Wi/Wo ratio on the discharge performance of MPKWs.
文摘In meandering rivers, the flow pattern is highly complex, with specific characteristics at bends that are not observed along straight paths. A numerical model can be effectively used to predict such flow fields. Since river bends are not uniform-some are divergent and others convergent-in this study, after the SSIIM 3-D model was calibrated using the result of measurements along a uniform 180° bend with a width of 0.6 m, a similar but convergent 180° bend, 0.6 m to 0.45 m wide, was simulated using the SSIIM 3-D numerical model. Flow characteristics of the convergent 180° bend, including lengthwise and vertical velocity profiles, primary and secondary flows, lengthwise and widthwise slopes of the water surface, and the helical flow strength, were compared with those of the uniform 180° bend. The verification results of the model show that the numerical model can effectively simulate the flow field in the uniform bend. In addition, this research indicates that, in a convergent channel, the maximum velocity path at a plane near the water surface crosses the channel's centerline at about a 30° to 40° cross-section, while in the uniform bend, this occurs at about the 50° cross-section. The varying range of the water surface elevation is wider in the convergent channel than in the uniform one, and the strength of the helical flow is generally greater in the uniform channel than in the convergent one. Also, unlike the uniform bend, the convergent bend exhibits no rotational cell against the main direction of secondary flow rotation at the 135° cross-section.