In order to determine the slip plane of slope directly by the calculation results of strength reduction method, and analyze the influential factors of slope stability, a numerical model was established in plane strain...In order to determine the slip plane of slope directly by the calculation results of strength reduction method, and analyze the influential factors of slope stability, a numerical model was established in plane strain mode by FLAC3D for homogeneous soil slope, whose parameters were reduced until the slope reached the critical state. Then FISH program was used to get the location data of slip plane from displacement contour lines. Furthermore, the method to determine multiple slip planes was also proposed by setting different heights of elastic areas. The influential factors for the stability were analyzed, including cohesion, internal friction angle, and tensile strength. The calculation results show that with the increase of cohesion, failure mode of slope changes from shallow slipping to the deep slipping, while inclination of slip plane becomes slower and slipping volume becomes larger; with the increase of friction angle, failure mode of slope changes from deep slipping to shallow slipping, while slip plane becomes steeper and upper border of slip plane comes closer to the vertex of slope; the safety factor increases little and slip plane goes far away from vertex of slope with the increase of tensile strength.展开更多
Based on the research on rock burst phenomenon induced by the breakage of thick and hard roof around roadways and working faces in coal mines, a criterion of rock burst induced by roof breakage （RBRB） was proposed a...Based on the research on rock burst phenomenon induced by the breakage of thick and hard roof around roadways and working faces in coal mines, a criterion of rock burst induced by roof breakage （RBRB） was proposed and the model was built. Through the model, a method calculating the varied stresses induced by roof breakage in support objects and coal body was proposed and a unified formula was derived for the calculation of stress increment on support objects and coal body under different breaking forms of roof. Whilst the formula for calculating dynamic load was derived by introducing dynamic index Kd. The formula was verified in Huating Mine by stress measurement. According to the formula for stress increment calculating, the sensitivities of dynamic load parameters were further studied. The results show that the thickness and breaking depth of roof, width of support objeet are the sensitive factors. Based on the discussion of the model, six associated effective methods for rock burst prevention are obtained.展开更多
Based on the fact that a static problem has an equivalent wave speed of infinity and a dynamic problem has a wave speed of finite value, an effective loading algorithm associated with the explicit dynamic relaxation m...Based on the fact that a static problem has an equivalent wave speed of infinity and a dynamic problem has a wave speed of finite value, an effective loading algorithm associated with the explicit dynamic relaxation method was presented to produce meaningful numerical solutions for static problems. The central part of the explicit dynamic relaxation method is to turn a time-independent static problem into an artificial time-dependent dynamic problem. The related numerical testing results demonstrate that: (1) the proposed effective loading algorithm is capable of enabling an applied load in a static problem to be propagated throughout the whole system within a given loading increment, so that the time-independent solution of the static problem can be obtained; (2) the proposed effective loading algorithm can be straightforwardly applied to the particle simulation method for solving a wide range of static problems.展开更多
For calculating the thermal storage time for an annular tube with phase change material （PCM）, a novel method is proposed. The method is suitable for either low-temperature PCM or high-temperature PCM whose initial ...For calculating the thermal storage time for an annular tube with phase change material （PCM）, a novel method is proposed. The method is suitable for either low-temperature PCM or high-temperature PCM whose initial temperature is near the melting point. The deviation fit is smaller than 8% when the time is below 2x104 s. Comparison between the predictions and the reported experimental data of thermal storage time at same conditions is investigated and good agreements have been got. Based on this method, the performance of the thermal storage unit and the role of natural convection are also investigated. Results show a linear relation between the maximum amount of stored heat and thermal storage time, and their ratio increases with the height of the thermal storage unit. As the thickness of the cavity increases, natural convection plays an increasingly important role in promoting the melting behavior of paraffin. When the thickness of the cavity is small, natural convection restrains the melting behavior of paraffin.展开更多
In order to predict the local scour hole and its evaluation around a cylindrical bridge pier, the computational fluid dynamics （CFD） and theories of sediment movement and transport were employed to carry out numeric...In order to predict the local scour hole and its evaluation around a cylindrical bridge pier, the computational fluid dynamics （CFD） and theories of sediment movement and transport were employed to carry out numerical simulations. In the numerical method, the time-averaged Reynolds Navier-Stokes equations and the standard k-e model were first used to simulate the three-dimensional flow field around a bridge pier fixed on river bed. The transient shear stress on river bed was treated as a crucial hydrodynamic mechanism when handling sediment incipience and transport. Then, river-bed volumetric sediment transport was calculated, followed by the modification of the river bed altitude and configuration. Boundary adaptive mesh technique was employed to modify the grid system with changed river-bed boundary. The evolution of local scour around a cylindrical bridge pier was presented. The numerical results represent the flow pattern and mechanism during the pier scouring, with a good prediction of the maximum scour hole depth compared with test results.展开更多
A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forwar...A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.展开更多
A lightweight aggregate concrete-filled steel tube(LACFST) spatial truss beam was tested under bending load. The performance was studied by the analysis of the beam deflection and strains in its chords and webs. Accor...A lightweight aggregate concrete-filled steel tube(LACFST) spatial truss beam was tested under bending load. The performance was studied by the analysis of the beam deflection and strains in its chords and webs. According to the test results, several assumptions were made to deduce the bearing capacity calculation method based on the force balance of the whole section. An optimal dimension relationship for the truss beam chords was proposed and verified by finite element analysis. Results show that the LACFST spatial truss beam failed after excessive deflection. The strain distribution agreed with Bernoulli-Euler theoretical prediction. The truss beam flexural bearing capacity calculation results matched test evidence with only a 3% difference between the two. Finite element analyses with different chord dimensions show that the ultimate bearing capacity increases as the chord dimensions increase when the chords have a diameter smaller than optimal one; otherwise, it remains almost unchanged as the chord dimensions increase.展开更多
Ab initio calculations are used to understand the fundamental mechanism of the solid solution softening/hardening of the Mo-binary system.The results reveal that the Mo-Ti,Mo-Ta,Mo-Nb,and Mo-W interactions are primari...Ab initio calculations are used to understand the fundamental mechanism of the solid solution softening/hardening of the Mo-binary system.The results reveal that the Mo-Ti,Mo-Ta,Mo-Nb,and Mo-W interactions are primarily attractive with negative heats of formation,while the interactions of Mo-Re,and Mo-Zr would be mainly repulsive with positive heats of formation.It is also shown that the addition of Re and Zr would cause the solid solution softening of Mo by the decrease of the unstable stacking fault energy and the increase of ductility.On the contrary,the elements of W,Ta,Ti,and Nb could bring about the solid-solution hardening of Mo through the impediment of the slip of the dislocation and the decrease of ductility.Electronic structures indicate that the weaker/stronger chemical bonding due to the alloying elements should fundamentally induce the solid solution softening/hardening of Mo.The results are discussed and compared with available evidence in literatures,which could deepen the fundamental understanding of the solid solution softening/hardening of the binary metallic system.展开更多
A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditiona...A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditional numerical method of the same equations corroborates well the reliability and rate of FEFDM.Moreover,a flow rate estimate method was developed for the project whose injection rate has not been clearly determined.A wellhead pressure regime determined by this method was successfully applied to the trial injection operations in Shihezi formation of Shenhua CCS Project,which is a good practice verification of FEFDM.At last,this method was used to evaluate the effect of friction and acceleration terms on the flow equation on the wellhead pressure.The result shows that for deep wellbore,the friction term can be omitted when flow rate is low and in a wide range of velocity the acceleration term can always be deleted.It is also shown that with flow rate increasing,the friction term can no longer be neglected.展开更多
The bulk electronic structure of kaolinite (001) plane was studied with quantum mechanical calculations. The CASTEP parameterization of ultrasoft pseudopotentials without core corrections was used to optimize the stru...The bulk electronic structure of kaolinite (001) plane was studied with quantum mechanical calculations. The CASTEP parameterization of ultrasoft pseudopotentials without core corrections was used to optimize the structure of kaolinite bulk and slab models. The results show that Fermi energy of kaolinite (001) plane is 3.05 eV, and the band gap is 4.52 eV. The partial density of states (PDOS) of kaolinite (001) plane indicates that Al-O and Si-O bonds on the mineral surface are highly polar. The oxygen atoms of hydroxyl groups in surface layer are capable of forming hydrogen bond with the head group of cationic collectors. The properties of dodecylamine (DDA) cation were also calculated by density function theory (DFT) method at B3LYP/6-31G (d) level for illuminating the flotation processes of kaolinite. Besides the electrostatic attraction, the mechanism between kaolinite and DDA is found to be hydrogen bonds under acidic condition.展开更多
The machining principle and realization method for the continuous generative grinding face gear by a worm wheel are introduced. Based on a five-axis linked CNC grinding machine, a new method is presented to deprive th...The machining principle and realization method for the continuous generative grinding face gear by a worm wheel are introduced. Based on a five-axis linked CNC grinding machine, a new method is presented to deprive the equation of face gear error tooth surface by assuming the tool surface as the error surface, where actual tool installation position error is introduced into the equation of virtual shaper cutter. Surface equations and 3-D models for the face gear and the worm wheel involving four kinds of tool installation errors are established. When compared, the face gear tooth surface machined in VERICUT software for simulation based on this new method and the one obtained based on real process(grinding face gear by using a theoretical worm wheel with actual position errors) are found to be coincident, which proves the validity and feasibility of this new method. By using mesh planning for the rotating projection plane of the face gear work tooth surface, the deviation values of the tooth surface and the difference surface are acquired, and the influence of four kinds of errors on the face gear tooth surface is analyzed. Accordingly, this work provides a theoretical reference for assembly craft of worm wheel, improvement of face gear machining accuracy and modification of error tooth surface.展开更多
Analysis of the aerodynamic performance of high-speed trains in special cuts would provide references for the critical overturning velocity and complement the operation safety management under strong winds.This work w...Analysis of the aerodynamic performance of high-speed trains in special cuts would provide references for the critical overturning velocity and complement the operation safety management under strong winds.This work was conducted to investigate the flow structure around trains under different cut depths,slope angles using computational fluid dynamics(CFD).The high-speed train was considered with bogies and inter-carriage gaps.And the accuracy of the numerical method was validated by combining with the experimental data of wind tunnel tests.Then,the variations of aerodynamic forces and surface pressure distribution of the train were mainly analyzed.The results show that the surroundings of cuts along the railway line have a great effect on the crosswind stability of trains.With the slope angle and depth of the cut increasing,the coefficients of aerodynamic forces tend to reduce.An angle of 75°is chosen as the optimum one for the follow-up research.Under different depth conditions,the reasonable cut depth for high-speed trains to run safely is 3 m lower than that of the conventional cut whose slope ratio is 1:1.5.Furthermore,the windward slope angle is more important than the leeward one for the train aerodynamic performance.Due to the shield of appropriate cuts,the train body is in a minor positive pressure environment.Thus,designing a suitable cut can contribute to improving the operation safety of high-speed trains.展开更多
The hybrid genetic algorithm is utilized to facilitate model parameter estimation.The tri-dimensional compression tests of soil are performed to supply experimental data for identifying nonlinear constitutive model of...The hybrid genetic algorithm is utilized to facilitate model parameter estimation.The tri-dimensional compression tests of soil are performed to supply experimental data for identifying nonlinear constitutive model of soil.In order to save computing time during parameter inversion,a new procedure to compute the calculated strains is presented by multi-linear simplification approach instead of finite element method(FEM).The real-coded hybrid genetic algorithm is developed by combining normal genetic algorithm with gradient-based optimization algorithm.The numerical and experimental results for conditioned soil are compared.The forecast strains based on identified nonlinear constitutive model of soil agree well with observed ones.The effectiveness and accuracy of proposed parameter estimation approach are validated.展开更多
High quality mesh plays an important role for finite element methods in science computation and numerical simulation.Whether the mesh quality is good or not,to some extent,it determines the calculation results of the ...High quality mesh plays an important role for finite element methods in science computation and numerical simulation.Whether the mesh quality is good or not,to some extent,it determines the calculation results of the accuracy and efficiency.Different from classic Lloyd iteration algorithm which is convergent slowly,a novel accelerated scheme was presented,which consists of two core parts:mesh points replacement and local edges Delaunay swapping.By using it,almost all the equilateral triangular meshes can be generated based on centroidal Voronoi tessellation(CVT).Numerical tests show that it is significantly effective with time consuming decreasing by 40%.Compared with other two types of regular mesh generation methods,CVT mesh demonstrates that higher geometric average quality increases over 0.99.展开更多
According to the recently developed single-trough floating machine with the world's largest volume(inflatable mechanical agitation flotation machine with volume of 320 m3) in China, the gas-fluid two-phase flow in...According to the recently developed single-trough floating machine with the world's largest volume(inflatable mechanical agitation flotation machine with volume of 320 m3) in China, the gas-fluid two-phase flow in flotation cell was simulated using computational fluid dynamics method. It is shown that hexahedral mesh scheme is more suitable for the complex structure of the flotation cell than tetrahedral mesh scheme, and a mesh quality ranging from 0.7 to 1.0 is obtained. Comparative studies of the standard k-ε, k-ω and realizable k-ε turbulence models were carried out. It is indicated that the standard k-ε turbulence model could give a result relatively close to the practice and the liquid phase flow field is well characterized. In addition, two obvious recirculation zones are formed in the mixing zones, and the pressure on the rotor and stator is well characterized. Furthermore, the simulation results using improved standard k-ε turbulence model show that surface tension coefficient of 0.072, drag model of Grace and coefficient of 4, and lift coefficient of 0.001 can be achieved. The research results suggest that gas-fluid two-phase flow in large flotation cell can be well simulated using computational fluid dynamics method.展开更多
An empirical expression of cohesion （C） and friction angle （Ф） for layered rock was suggested. This expression was compared with a test result made by the former researchers. The constitutive relationship of a tr...An empirical expression of cohesion （C） and friction angle （Ф） for layered rock was suggested. This expression was compared with a test result made by the former researchers. The constitutive relationship of a transversely isotropic medium and Mohr-Coulomb criterion in which C and Ф vary with directions were employed, and a relative 3D elasto-plastic FEM code was developed, in which the important thing was to adopt a search-trial method to find the orientation angle （p） of shear failure plane （or weakest shear plane） with respect to the major principal stress as well as the corresponding C and Ф Taking an underground opening as the calculation object, the numerical analyses were carried out by using the FEM code for two cases of transversely isotropic rock and isotropic rock, respectively, and the computation results were compared. The results show that when the rock is a transversely isotropic one, the distributions of displacements, plastic zones and stress contours in the surrounding rock will be non-axisymmetric along the tunnel＇s vertical axis, which is very different from that of isotropic rock. The stability of the tunnel in transversely isotropic rock is relatively low.展开更多
An algorithm was proposed to fast recognize three types of underwater micro-terrain, i.e. the level, the gradient and the uneven. With pendulum single beam bathymeter, the hard level concrete floor, the random uneven ...An algorithm was proposed to fast recognize three types of underwater micro-terrain, i.e. the level, the gradient and the uneven. With pendulum single beam bathymeter, the hard level concrete floor, the random uneven floor and the gradient wood panel (8-) were ultrasonically detected 20 times, respectively. The results show that the algorithm is right from fact that the first clustering values of the uneven are all less than the threshold value of 60.0% that is obtained by the level and gradient samples. The algorithm based on the dynamic clustering theory can effectively eliminate the influences of the exceptional elevation values produced by the disturbances resulted from the grazing angle, the characteristic of bottom material and environmental noises, and its real-time capability is good. Thus, the algorithm provides a foundation for the next restructuring of the micro-terrain.展开更多
基金Project(20060533071) supported by the Doctoral Program Foundation of Higher Education of ChinaProject (20060400264) supported by China Postdoctoral Science Foundation+1 种基金Project (50774093) supported by the National Natural Science Foundation of ChinaProject (1343-74236000014) supported by Graduate Student Innovation Foundation of Hunan Province, China
文摘In order to determine the slip plane of slope directly by the calculation results of strength reduction method, and analyze the influential factors of slope stability, a numerical model was established in plane strain mode by FLAC3D for homogeneous soil slope, whose parameters were reduced until the slope reached the critical state. Then FISH program was used to get the location data of slip plane from displacement contour lines. Furthermore, the method to determine multiple slip planes was also proposed by setting different heights of elastic areas. The influential factors for the stability were analyzed, including cohesion, internal friction angle, and tensile strength. The calculation results show that with the increase of cohesion, failure mode of slope changes from shallow slipping to the deep slipping, while inclination of slip plane becomes slower and slipping volume becomes larger; with the increase of friction angle, failure mode of slope changes from deep slipping to shallow slipping, while slip plane becomes steeper and upper border of slip plane comes closer to the vertex of slope; the safety factor increases little and slip plane goes far away from vertex of slope with the increase of tensile strength.
基金Project(2010CB226805) supported by the National Basic Research Program of ChinaProjects(30370412,30670558) supported by the National Natural Science Foundation of ChinaProject(SKLCRSM10X05) supported by the Self-research Program of the Key Laboratory of Coal Resources and Safe Mining,China
文摘Based on the research on rock burst phenomenon induced by the breakage of thick and hard roof around roadways and working faces in coal mines, a criterion of rock burst induced by roof breakage （RBRB） was proposed and the model was built. Through the model, a method calculating the varied stresses induced by roof breakage in support objects and coal body was proposed and a unified formula was derived for the calculation of stress increment on support objects and coal body under different breaking forms of roof. Whilst the formula for calculating dynamic load was derived by introducing dynamic index Kd. The formula was verified in Huating Mine by stress measurement. According to the formula for stress increment calculating, the sensitivities of dynamic load parameters were further studied. The results show that the thickness and breaking depth of roof, width of support objeet are the sensitive factors. Based on the discussion of the model, six associated effective methods for rock burst prevention are obtained.
基金Projects(10872219 10672190) supported by the National Natural Science Foundation of China
文摘Based on the fact that a static problem has an equivalent wave speed of infinity and a dynamic problem has a wave speed of finite value, an effective loading algorithm associated with the explicit dynamic relaxation method was presented to produce meaningful numerical solutions for static problems. The central part of the explicit dynamic relaxation method is to turn a time-independent static problem into an artificial time-dependent dynamic problem. The related numerical testing results demonstrate that: (1) the proposed effective loading algorithm is capable of enabling an applied load in a static problem to be propagated throughout the whole system within a given loading increment, so that the time-independent solution of the static problem can be obtained; (2) the proposed effective loading algorithm can be straightforwardly applied to the particle simulation method for solving a wide range of static problems.
基金Projects(51666006,51406071,51174105,51366005)supported by the National Natural Science Foundation of ChinaProject(2014CB460605)supported by the National Basic Research Program of China
文摘For calculating the thermal storage time for an annular tube with phase change material （PCM）, a novel method is proposed. The method is suitable for either low-temperature PCM or high-temperature PCM whose initial temperature is near the melting point. The deviation fit is smaller than 8% when the time is below 2x104 s. Comparison between the predictions and the reported experimental data of thermal storage time at same conditions is investigated and good agreements have been got. Based on this method, the performance of the thermal storage unit and the role of natural convection are also investigated. Results show a linear relation between the maximum amount of stored heat and thermal storage time, and their ratio increases with the height of the thermal storage unit. As the thickness of the cavity increases, natural convection plays an increasingly important role in promoting the melting behavior of paraffin. When the thickness of the cavity is small, natural convection restrains the melting behavior of paraffin.
基金Project(50978095) supported by the National Natural Science Foundation of ChinaProject(IRT0917) supported by the Program for Changjiang Scholars and Innovative Research Team in Chinese UniversityProject supported by China Scholarship Council
文摘In order to predict the local scour hole and its evaluation around a cylindrical bridge pier, the computational fluid dynamics （CFD） and theories of sediment movement and transport were employed to carry out numerical simulations. In the numerical method, the time-averaged Reynolds Navier-Stokes equations and the standard k-e model were first used to simulate the three-dimensional flow field around a bridge pier fixed on river bed. The transient shear stress on river bed was treated as a crucial hydrodynamic mechanism when handling sediment incipience and transport. Then, river-bed volumetric sediment transport was calculated, followed by the modification of the river bed altitude and configuration. Boundary adaptive mesh technique was employed to modify the grid system with changed river-bed boundary. The evolution of local scour around a cylindrical bridge pier was presented. The numerical results represent the flow pattern and mechanism during the pier scouring, with a good prediction of the maximum scour hole depth compared with test results.
基金Project(60672042) supported by the National Natural Science Foundation of China
文摘A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.
基金Project(51208176)supported by the National Natural Science Foundation of ChinaProjects(2012M511187,2013T60493)supported by the China Postdoctoral Science FoundationProject(2015B17414)supported by the Fundamental Research Funds for the Central Universities,China
文摘A lightweight aggregate concrete-filled steel tube(LACFST) spatial truss beam was tested under bending load. The performance was studied by the analysis of the beam deflection and strains in its chords and webs. According to the test results, several assumptions were made to deduce the bearing capacity calculation method based on the force balance of the whole section. An optimal dimension relationship for the truss beam chords was proposed and verified by finite element analysis. Results show that the LACFST spatial truss beam failed after excessive deflection. The strain distribution agreed with Bernoulli-Euler theoretical prediction. The truss beam flexural bearing capacity calculation results matched test evidence with only a 3% difference between the two. Finite element analyses with different chord dimensions show that the ultimate bearing capacity increases as the chord dimensions increase when the chords have a diameter smaller than optimal one; otherwise, it remains almost unchanged as the chord dimensions increase.
基金Project(51801129)supported by the National Natural Science Foundation of ChinaProject supported by the State Key Laboratory of Powder Metallurgy,China。
文摘Ab initio calculations are used to understand the fundamental mechanism of the solid solution softening/hardening of the Mo-binary system.The results reveal that the Mo-Ti,Mo-Ta,Mo-Nb,and Mo-W interactions are primarily attractive with negative heats of formation,while the interactions of Mo-Re,and Mo-Zr would be mainly repulsive with positive heats of formation.It is also shown that the addition of Re and Zr would cause the solid solution softening of Mo by the decrease of the unstable stacking fault energy and the increase of ductility.On the contrary,the elements of W,Ta,Ti,and Nb could bring about the solid-solution hardening of Mo through the impediment of the slip of the dislocation and the decrease of ductility.Electronic structures indicate that the weaker/stronger chemical bonding due to the alloying elements should fundamentally induce the solid solution softening/hardening of Mo.The results are discussed and compared with available evidence in literatures,which could deepen the fundamental understanding of the solid solution softening/hardening of the binary metallic system.
基金Project(Z110803)supported by the State Key Laboratory of Geomechanics and Geotechnical Engineering,ChinaProject(2008AA062303)supported by the National High Technology Research and Development Program of China
文摘A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditional numerical method of the same equations corroborates well the reliability and rate of FEFDM.Moreover,a flow rate estimate method was developed for the project whose injection rate has not been clearly determined.A wellhead pressure regime determined by this method was successfully applied to the trial injection operations in Shihezi formation of Shenhua CCS Project,which is a good practice verification of FEFDM.At last,this method was used to evaluate the effect of friction and acceleration terms on the flow equation on the wellhead pressure.The result shows that for deep wellbore,the friction term can be omitted when flow rate is low and in a wide range of velocity the acceleration term can always be deleted.It is also shown that with flow rate increasing,the friction term can no longer be neglected.
基金Project(2005CB623701) supported by the Major State Basic Research and Development Program of ChinaProject(50874118) supported by the National Nature Science Foundation of ChinaProject(2007B52) supported by the Foundation for the Author of National Excellent Doctoral Dissertation of China
文摘The bulk electronic structure of kaolinite (001) plane was studied with quantum mechanical calculations. The CASTEP parameterization of ultrasoft pseudopotentials without core corrections was used to optimize the structure of kaolinite bulk and slab models. The results show that Fermi energy of kaolinite (001) plane is 3.05 eV, and the band gap is 4.52 eV. The partial density of states (PDOS) of kaolinite (001) plane indicates that Al-O and Si-O bonds on the mineral surface are highly polar. The oxygen atoms of hydroxyl groups in surface layer are capable of forming hydrogen bond with the head group of cationic collectors. The properties of dodecylamine (DDA) cation were also calculated by density function theory (DFT) method at B3LYP/6-31G (d) level for illuminating the flotation processes of kaolinite. Besides the electrostatic attraction, the mechanism between kaolinite and DDA is found to be hydrogen bonds under acidic condition.
基金Projects(51535012,U1604255)supported by the National Natural Science Foundation of ChinaProject(2016JC2001)supported by the Key Research and Development Project of Hunan Province,China
文摘The machining principle and realization method for the continuous generative grinding face gear by a worm wheel are introduced. Based on a five-axis linked CNC grinding machine, a new method is presented to deprive the equation of face gear error tooth surface by assuming the tool surface as the error surface, where actual tool installation position error is introduced into the equation of virtual shaper cutter. Surface equations and 3-D models for the face gear and the worm wheel involving four kinds of tool installation errors are established. When compared, the face gear tooth surface machined in VERICUT software for simulation based on this new method and the one obtained based on real process(grinding face gear by using a theoretical worm wheel with actual position errors) are found to be coincident, which proves the validity and feasibility of this new method. By using mesh planning for the rotating projection plane of the face gear work tooth surface, the deviation values of the tooth surface and the difference surface are acquired, and the influence of four kinds of errors on the face gear tooth surface is analyzed. Accordingly, this work provides a theoretical reference for assembly craft of worm wheel, improvement of face gear machining accuracy and modification of error tooth surface.
基金Projects(51075401,U1334205)supported by the National Natural Science Foundation of ChinaProject supported by the Scholarship Award for Excellent Innovative Doctoral Student granted by Central South University of ChinaProject(132014)supported by the Fok Ying Tong Education Foundation,China
文摘Analysis of the aerodynamic performance of high-speed trains in special cuts would provide references for the critical overturning velocity and complement the operation safety management under strong winds.This work was conducted to investigate the flow structure around trains under different cut depths,slope angles using computational fluid dynamics(CFD).The high-speed train was considered with bogies and inter-carriage gaps.And the accuracy of the numerical method was validated by combining with the experimental data of wind tunnel tests.Then,the variations of aerodynamic forces and surface pressure distribution of the train were mainly analyzed.The results show that the surroundings of cuts along the railway line have a great effect on the crosswind stability of trains.With the slope angle and depth of the cut increasing,the coefficients of aerodynamic forces tend to reduce.An angle of 75°is chosen as the optimum one for the follow-up research.Under different depth conditions,the reasonable cut depth for high-speed trains to run safely is 3 m lower than that of the conventional cut whose slope ratio is 1:1.5.Furthermore,the windward slope angle is more important than the leeward one for the train aerodynamic performance.Due to the shield of appropriate cuts,the train body is in a minor positive pressure environment.Thus,designing a suitable cut can contribute to improving the operation safety of high-speed trains.
基金Project(2007CB714006) supported by the National Basic Research Program of China Project(90815023) supported by the National Natural Science Foundation of China
文摘The hybrid genetic algorithm is utilized to facilitate model parameter estimation.The tri-dimensional compression tests of soil are performed to supply experimental data for identifying nonlinear constitutive model of soil.In order to save computing time during parameter inversion,a new procedure to compute the calculated strains is presented by multi-linear simplification approach instead of finite element method(FEM).The real-coded hybrid genetic algorithm is developed by combining normal genetic algorithm with gradient-based optimization algorithm.The numerical and experimental results for conditioned soil are compared.The forecast strains based on identified nonlinear constitutive model of soil agree well with observed ones.The effectiveness and accuracy of proposed parameter estimation approach are validated.
基金Project(11002121) supported by the National Natural Science Foundation of ChinaProject(09QDZ09) supported by Doctor Foundation of Xiangtan University, China+2 种基金Project(2009LCSSE11) supported by Hunan Key Laboratory for CSSE, ChinaProject(2011FJ3231) supported by Planned Science and Technology Project of Hunan Province,ChinaProject(12JJ3054) supported by the Provincial Natural Science Foundation of Hunan,China
文摘High quality mesh plays an important role for finite element methods in science computation and numerical simulation.Whether the mesh quality is good or not,to some extent,it determines the calculation results of the accuracy and efficiency.Different from classic Lloyd iteration algorithm which is convergent slowly,a novel accelerated scheme was presented,which consists of two core parts:mesh points replacement and local edges Delaunay swapping.By using it,almost all the equilateral triangular meshes can be generated based on centroidal Voronoi tessellation(CVT).Numerical tests show that it is significantly effective with time consuming decreasing by 40%.Compared with other two types of regular mesh generation methods,CVT mesh demonstrates that higher geometric average quality increases over 0.99.
基金Project(51074027)supported by the National Natural Science Foundation of China
文摘According to the recently developed single-trough floating machine with the world's largest volume(inflatable mechanical agitation flotation machine with volume of 320 m3) in China, the gas-fluid two-phase flow in flotation cell was simulated using computational fluid dynamics method. It is shown that hexahedral mesh scheme is more suitable for the complex structure of the flotation cell than tetrahedral mesh scheme, and a mesh quality ranging from 0.7 to 1.0 is obtained. Comparative studies of the standard k-ε, k-ω and realizable k-ε turbulence models were carried out. It is indicated that the standard k-ε turbulence model could give a result relatively close to the practice and the liquid phase flow field is well characterized. In addition, two obvious recirculation zones are formed in the mixing zones, and the pressure on the rotor and stator is well characterized. Furthermore, the simulation results using improved standard k-ε turbulence model show that surface tension coefficient of 0.072, drag model of Grace and coefficient of 4, and lift coefficient of 0.001 can be achieved. The research results suggest that gas-fluid two-phase flow in large flotation cell can be well simulated using computational fluid dynamics method.
基金Project(2010CB732101) supported by the National Basic Research Program of China Project(51079145) supported by the National Natural Science Foundation of China
文摘An empirical expression of cohesion （C） and friction angle （Ф） for layered rock was suggested. This expression was compared with a test result made by the former researchers. The constitutive relationship of a transversely isotropic medium and Mohr-Coulomb criterion in which C and Ф vary with directions were employed, and a relative 3D elasto-plastic FEM code was developed, in which the important thing was to adopt a search-trial method to find the orientation angle （p） of shear failure plane （or weakest shear plane） with respect to the major principal stress as well as the corresponding C and Ф Taking an underground opening as the calculation object, the numerical analyses were carried out by using the FEM code for two cases of transversely isotropic rock and isotropic rock, respectively, and the computation results were compared. The results show that when the rock is a transversely isotropic one, the distributions of displacements, plastic zones and stress contours in the surrounding rock will be non-axisymmetric along the tunnel＇s vertical axis, which is very different from that of isotropic rock. The stability of the tunnel in transversely isotropic rock is relatively low.
基金Project(50474052) supported by the National Natural Foundation of China
文摘An algorithm was proposed to fast recognize three types of underwater micro-terrain, i.e. the level, the gradient and the uneven. With pendulum single beam bathymeter, the hard level concrete floor, the random uneven floor and the gradient wood panel (8-) were ultrasonically detected 20 times, respectively. The results show that the algorithm is right from fact that the first clustering values of the uneven are all less than the threshold value of 60.0% that is obtained by the level and gradient samples. The algorithm based on the dynamic clustering theory can effectively eliminate the influences of the exceptional elevation values produced by the disturbances resulted from the grazing angle, the characteristic of bottom material and environmental noises, and its real-time capability is good. Thus, the algorithm provides a foundation for the next restructuring of the micro-terrain.
