The utilization of prefabricated light modular radiant heating system has demonstrated significant increases in heat transfer efficiency and energy conservation capabilities.Within prefabricated building construction,...The utilization of prefabricated light modular radiant heating system has demonstrated significant increases in heat transfer efficiency and energy conservation capabilities.Within prefabricated building construction,this new heating method presents an opportunity for the development of comprehensive facilities.The parameters for evaluating the effectiveness of such a system are the upper surface layer’s heat flux and temperature.In this paper,thermal resistance analysis calculation based on a simplified model for this unique radiant heating system analysis is presented with the heat transfer mechanism’s evaluation.The results obtained from thermal resistance analysis calculation and numerical simulation indicate that the thermal resistance analysis method is highly accurate with temperature discrepancies ranging from 0.44℃ to−0.44℃ and a heat flux discrepancy of less than 7.54%,which can meet the requirements of practical engineering applications,suggesting a foundation for the prefabricated radiant heating system.展开更多
This study considers an MHD Jeffery-Hamel nanofluid flow with distinct nanoparticles such as copper,Al_(2)O_(3)and SiO_(2)between two rigid non-parallel plane walls with the fuzzy extension of the generalized dual par...This study considers an MHD Jeffery-Hamel nanofluid flow with distinct nanoparticles such as copper,Al_(2)O_(3)and SiO_(2)between two rigid non-parallel plane walls with the fuzzy extension of the generalized dual parametric homotopy algorithm.The nanofluids have been formulated to enhance the thermophysical characteristics of fluids,including thermal diffusivity,conductivity,convective heat transfer coefficients and viscosity.Due to the presence of distinct nanofluids,a change in the value of volume fraction occurs that influences the velocity profiles of the flow.The short value of nanoparticles volume fraction is considered an uncertain parameter and represented in a triangular fuzzy number range among[0.0,0.1,0.2].A novel generalized dual parametric homotopy algorithm with fuzzy extension is used here to study the fuzzy velocities at various channel positions.Finally,the effectiveness of the proposed approach has been demonstrated through a comparison with the available results in the crisp case.展开更多
Understanding the steady mechanism of biomass smoldering plays a great role in the utilization of smoldering technology.In this study numerical analysis of steady smoldering of biomass rods was performed.A two-dimensi...Understanding the steady mechanism of biomass smoldering plays a great role in the utilization of smoldering technology.In this study numerical analysis of steady smoldering of biomass rods was performed.A two-dimensional(2D)steady model taking into account both char oxidation and pyrolysis was developed on the basis of a calculated propagation velocity according to empirical correlation.The model was validated against the smoldering experiment of biomass rods under natural conditions,and the maximum error was smaller than 31%.Parameter sensitivity analysis found that propagation velocity decreases significantly while oxidation area and pyrolysis zone increase significantly with the increasing diameter of rod fuel.展开更多
The center cutter of a hard rock tunnel boring machine(TBM) is installed on the cutterhead at a small radius and thus bears complex side force.Given this fact,the formation mechanism and change law of the side force s...The center cutter of a hard rock tunnel boring machine(TBM) is installed on the cutterhead at a small radius and thus bears complex side force.Given this fact,the formation mechanism and change law of the side force suffered by the center cutter were studied.Based on the rock shear failure criterion in combination with the lateral rolling width,a model for predicting the average side force was set up.Besides,a numerical analysis model of the rock fragmentation of the center cutter was established,and the instantaneous load changing features were investigated.Results shows that the inner side of the center cutter can form lateral rolling annulus in rock during the rotary cutting process.The smaller the installation radius is,the greater the cutter side force will be.In a working condition,the side force of the innermost center cutter is 11.66 k N,while it decreases sharply when installation radius increases.Variation tends to be gentle when installation radius is larger than 500 mm,and the side force of the outermost center cutter is reduced to 0.74 k N.展开更多
Wellbore stability analysis is a growing concern in oil industries. There are many parameters affecting the stability of a wellbore including geomechanical properties (e.g., elastic modulus, uni-axial compressive stre...Wellbore stability analysis is a growing concern in oil industries. There are many parameters affecting the stability of a wellbore including geomechanical properties (e.g., elastic modulus, uni-axial compressive strength (UCS) and cohesion) and acting forces (e.g., field stresses and mud pressure). Accurate determination of these parameters is time-consuming, expensive and sometimes even impossible. This work offers a systematic sensitivity analysis to quantify the amount of each parameter’s effect on the stability of a wellbore. Maximum wellbore wall displacement is used as a stability factor to study the stability of a wellbore. A 3D finite difference method with Mohr model is used for the numerical modeling. The numerical model is verified against an analytical solution. A dimensionless sensitivity factor is developed in order to compare the results of various parameters in the sensitivity analysis. The results show a different order of importance of parameters based on rock strength. The most sensitive properties for a weak rock are the maximum horizontal stress, internal friction angle and formation pressure, respectively, while for a strong rock, the most sensitive parameters are the maximum horizontal stress, mud pressure and pore pressure, respectively. The amount of error in wellbore stability analysis inflicted by the error in estimation of each parameter was also derived.展开更多
A coupled thermo-hydro-mechanical-migratory model of dual-porosity medium for saturated-unsaturated ubiquitous-joint rockmass was established,in which the stress field and the temperature field were single,but the see...A coupled thermo-hydro-mechanical-migratory model of dual-porosity medium for saturated-unsaturated ubiquitous-joint rockmass was established,in which the stress field and the temperature field were single,but the seepage field and the concentration field were double,and the influences of sets,spaces,angles,continuity ratios,stiffnesses of fractures on the constitutive relationship of the medium were considered.Also,the relative two-dimensional program of finite element method was developed.Taking a hypothetical nuclear waste repository as a calculation example,the case in which the rockmass was unsaturated dual-porosity medium and radioactive nuclide leak was simulated numerically,and the temperatures,negative pore pressures,saturations,flow velocities,nuclide concentrations and principal stresses in the rockmass were investigated.The results show that the negative pore pressures and nuclide concentrations in the porosity and fracture present different changes and distributions.Even though the saturation degree in porosity is only about 1/10 that in fracture,the flow velocity of underground water in fracture is about three times that in porosity because the permeability coefficient of fracture is almost four orders higher than that of porosity.The value of nuclide concentration in fracture is close to that in porosity.展开更多
There were differences between real boundary and blast hole controlling boundary of irregular mined-out area in underground metal mines. There were errors in numerical analysis of stability for goaf, if it was analyze...There were differences between real boundary and blast hole controlling boundary of irregular mined-out area in underground metal mines. There were errors in numerical analysis of stability for goaf, if it was analyzed as regular 3D mined-out area and the influence of coupling stress-seepage-disturbance was not considered adequately. Taking a lead zinc mine as the background, the model was built by the coupling of Surpac and Midas-Gts based on the goaf model precisely measured by CMS.According to seepage stress fundamental equations based on the equivalent continuum mechanical and the theory about equivalent load of dynamic disturbance in deep-hole blasting, the stability of mined-out area under multi-field coupling of stress-seepage-dynamic disturbance was numerically analyzed. The results show that it is more consistent between the numerical analysis model based on the real model of irregular 3D shape goaf and the real situation, which could faithfully reappear the change rule of stress–strain about the surrounding rock under synthetic action of blasting dynamic loading and the seepage pressure. The mined-out area multi-field coupling formed by blasting excavation is stable. Based on combination of the advantages of the CMS,Surpac and Midas-Gts, and fully consideration of the effects of multi-field coupling, the accurate and effective way could be provided for numerical analysis of stability for mined-out area.展开更多
In order to study the safety factor and instability process of cohesive soil slope, the discrete element method(DEM) was applied. DEM software PFC2 D was used to simulate the triaxial test to study the influence of th...In order to study the safety factor and instability process of cohesive soil slope, the discrete element method(DEM) was applied. DEM software PFC2 D was used to simulate the triaxial test to study the influence of the particle micro parameters on the macroscopic characteristics of cohesive soil and calibrate the micro parameters of DEM model on this basis. Embankment slope stability analysis was carried out by strength reduction and gravity increase method, it is shown that the safety factor obtained by strength reduction method is more conservative, and the arc-shaped feature of the sliding surface under the gravity increase method is more obvious. Throughout the progressive failure process, the failure trends, maximum displacements, and velocity changes obtained by the two methods were consistent. When slope was destroyed, the upper part was cracked, the middle part was sheared, and the lower part was destroyed by extrusion. The conclusions of this paper can be applied to the safety factor calculation of cohesive soil slopes and the analysis of the instability process.展开更多
The distributed shear stress and the displacement across shear band, the evolution of plastic zones, and the load-carrying capacity of rock specimen were investigated in plane strain direct shear test according to Fas...The distributed shear stress and the displacement across shear band, the evolution of plastic zones, and the load-carrying capacity of rock specimen were investigated in plane strain direct shear test according to Fast Lagrangian Analysis of Continua （FLAC）. And then the shear displacement distribution in normal direction of system composed of localized shear band and elastic rock was analyzed based on gradient-dependent plasticity. The adopted failure criterion was a composite of Mohr-Coulomb criterion, that is, the relation between tension cut-off and postpeak constitutive of rock was linear strain-softening. Numerical results show that shear stress field approximately undergoes three different stages. At first, shear stress is only concentrated in the middle of top and base of specimen. Next, shear stress in the middle of specimen tends to increase, owing to superposition of shear stresses. Interestingly, two peaks of shear stress appear far from the loading ends of specimen, and the peaks approach with the increase in timestep until elements at the center of specimen yield. Finally, relatively lower shear stress level is reached in large part of specimen except in the regions near the two ends. As flow stress decreases, the analytical shear displacement distribution in shear band based on gradient-dependent plasticity becomes steeps outside the band, it is linear and its slope tends to decrease. These theoretical results qualitatively agree with that of the present numerical predicted results. Main advantage of the analytical solution over the numerical results according to FLAC is that it is continuous, smooth and non-linear （except at elastic stage）.展开更多
A novel spiral non-circular bevel gear that could be applied to variable-speed driving in intersecting axes was proposed by combining the design principles of non-circular bevel gears and the manufacturing principles ...A novel spiral non-circular bevel gear that could be applied to variable-speed driving in intersecting axes was proposed by combining the design principles of non-circular bevel gears and the manufacturing principles of face-milling spiral bevel gears.Unlike straight non-circular bevel gears,spiral non-circular bevel gears have numerous advantages,such as a high contact ratio,high intensity,good dynamic performance,and an adjustable contact region.In addition,while manufacturing straight non-circular bevel gears is difficult,spiral non-circular bevel gears can be efficiently and precisely fabricated with a 6-axis bevel gear cutting machine.First,the generating principles of spiral non-circular bevel gears were introduced.Next,a mathematical model,including a generating tooth profile,tooth spiral,pressure angle,and generated tooth profile for this gear type was established.Then the precision of the model was verified by a tooth contact analysis using FEA,and the contact patterns and stress distributions of the spiral non-circular bevel gears were investigated.展开更多
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.展开更多
Aiming at accuracy control of the thermal crown of work rolls in cold rolling,new parameters such as regulation domain and control-efficiency factors were proposed and a numerical analysis model of the thermal crown o...Aiming at accuracy control of the thermal crown of work rolls in cold rolling,new parameters such as regulation domain and control-efficiency factors were proposed and a numerical analysis model of the thermal crown of work rolls was established using finite difference method to study roll's thermal deformation.Based on simulation results,the influences of control-efficiency factors on thermal crown are presented and the thermal crown of work rolls is analyzed after taking sub-cooling of sprinkling beam into consideration.It has been found that the control-efficiency factor of any position on the roll's surface is linear function of the temperature and the control ability of water temperature is stronger than other control parameters.In addition,the verification of the model has been carried out based on the producing technology data in some factories and the numerical simulation results coincide well with the experimental data.Therefore,this work has important value for on-line control of roll's crown in cold rolling.展开更多
Dimensional analysis and numerical simulations were carried out to research prediction method of breakthrough time of horizontal wells in bottom water reservoir. Four dimensionless independent variables and dimensionl...Dimensional analysis and numerical simulations were carried out to research prediction method of breakthrough time of horizontal wells in bottom water reservoir. Four dimensionless independent variables and dimensionless time were derived from 10 influencing factors of the problem by using dimensional analysis. Simulations of horizontal well in reservoir with bottom water were run to find the prediction correlation. A general and concise functional relationship for predicting breakthrough time was established based on simulation results and theoretical analysis. The breakthrough time of one conceptual model predicted by the correlation is very close to the result by Eclipse with less than 2% error. The practical breakthrough time of one well in Helder oilfield is 10 d, and the predicted results by the method is 11.2 d, which is more accurate than the analytical result. Case study indicates that the method could predict breakthrough time of horizontal well under different reservoir conditions accurately. For its university and ease of use, the method is suitable for quick prediction of breakthrough time.展开更多
A footing may get an eccentric load caused by earthquake or wind, thus the bearing capacity of footing subjected to eccentric load become a fundamental geotechnical problem. The conventional limit equilibrium method u...A footing may get an eccentric load caused by earthquake or wind, thus the bearing capacity of footing subjected to eccentric load become a fundamental geotechnical problem. The conventional limit equilibrium method used for this problem usually evaluates the material properties only by its final strength. But the classical finite element method(FEM) does not necessarily provide a clear collapse mechanism associated with the yield condition of elements. To overcome these defects, a numerical procedure is proposed to create an explicit collapse mode combining a modified smeared shear band approach with a modified initial stress method. To understand the practical performance of sand foundation and verify the performance of the proposed procedure applied to the practical problems, the computing results were compared with the laboratory model tests results and some conventional solutions. Furthermore, because the proposed numerical procedure employs a simple elasto-plastic model which requires a small number of soil parameters, it may be applied directly to practical design works.展开更多
The parameter sensitivities affecting the flutter speed of the NREL （National Renewable Energy Laboratory） 5-MW baseline HAWT （horizontal axis wind turbine） blades are analyzed. An aeroelastic model, which compris...The parameter sensitivities affecting the flutter speed of the NREL （National Renewable Energy Laboratory） 5-MW baseline HAWT （horizontal axis wind turbine） blades are analyzed. An aeroelastic model, which comprises an aerodynamic part to calculate the aerodynamic loads and a structural part to determine the structural dynamic responses, is established to describe the classical flutter of the blades. For the aerodynamic part, Theodorsen unsteady aerodynamics model is used. For the structural part, Lagrange’s equation is employed. The flutter speed is determined by introducing “V–g” method to the aeroelastic model, which converts the issue of classical flutter speed determination into an eigenvalue problem. Furthermore, the time domain aeroelastic response of the wind turbine blade section is obtained with employing Runge-Kutta method. The results show that four cases （i.e., reducing the blade torsional stiffness, moving the center of gravity or the elastic axis towards the trailing edge of the section, and placing the turbine in high air density area） will decrease the flutter speed. Therefore, the judicious selection of the four parameters （the torsional stiffness, the chordwise position of the center of gravity, the elastic axis position and air density） can increase the relative inflow speed at the blade section associated with the onset of flutter.展开更多
The finite element analysis （FEA） software Ansys was employed to study the stress state of the dies of both plane and non-plane parting face structures with uniform interference and the die of plane parting face str...The finite element analysis （FEA） software Ansys was employed to study the stress state of the dies of both plane and non-plane parting face structures with uniform interference and the die of plane parting face structure with non-uniform interference. Considering the symmetry of the die, a half gear tooth model of the two-ring assembled die with 2.5 GPa inner pressure was constructed. Four paths were defined to investigate the stress state at the bottom comer of the die where stress concentration was serious. FEA results show that the change of parting face from non-plane to plane can greatly reduce the stress at the teeth tips of the die so that the tip fracture is avoided. The interference structure of the die is the most important influencing factor for the stress concentration at the bottom comer. When non-uniform interference is adopted the first principal stress at the comer on the defined paths of the die is much lower than that with uniform interference. The bottom hole radius is another important influencing factor for the comer stress concentration. The first principal stress at the comer of the plane parting face die with non-uniform interference is reduced from 2.3 to 1.9 GPa when the hole radius increases from 12.5 to 16.0 mm. The optimization of the die structure increases the life of the die from 100 to 6 000 hits.展开更多
基金Project(NB-2020-JG-07)supported by the Research and Engineering Application of Key Technologies for New Building Industrialization Project of China Northwest Architectural Design and Research Institute Co.,Ltd.Project(2023-CXTD-29)supported by the Key Scientific and Technological Innovation Team of Shaanxi Province,ChinaProject supported by the K.C.Wong Education Foundation。
文摘The utilization of prefabricated light modular radiant heating system has demonstrated significant increases in heat transfer efficiency and energy conservation capabilities.Within prefabricated building construction,this new heating method presents an opportunity for the development of comprehensive facilities.The parameters for evaluating the effectiveness of such a system are the upper surface layer’s heat flux and temperature.In this paper,thermal resistance analysis calculation based on a simplified model for this unique radiant heating system analysis is presented with the heat transfer mechanism’s evaluation.The results obtained from thermal resistance analysis calculation and numerical simulation indicate that the thermal resistance analysis method is highly accurate with temperature discrepancies ranging from 0.44℃ to−0.44℃ and a heat flux discrepancy of less than 7.54%,which can meet the requirements of practical engineering applications,suggesting a foundation for the prefabricated radiant heating system.
文摘This study considers an MHD Jeffery-Hamel nanofluid flow with distinct nanoparticles such as copper,Al_(2)O_(3)and SiO_(2)between two rigid non-parallel plane walls with the fuzzy extension of the generalized dual parametric homotopy algorithm.The nanofluids have been formulated to enhance the thermophysical characteristics of fluids,including thermal diffusivity,conductivity,convective heat transfer coefficients and viscosity.Due to the presence of distinct nanofluids,a change in the value of volume fraction occurs that influences the velocity profiles of the flow.The short value of nanoparticles volume fraction is considered an uncertain parameter and represented in a triangular fuzzy number range among[0.0,0.1,0.2].A novel generalized dual parametric homotopy algorithm with fuzzy extension is used here to study the fuzzy velocities at various channel positions.Finally,the effectiveness of the proposed approach has been demonstrated through a comparison with the available results in the crisp case.
文摘Understanding the steady mechanism of biomass smoldering plays a great role in the utilization of smoldering technology.In this study numerical analysis of steady smoldering of biomass rods was performed.A two-dimensional(2D)steady model taking into account both char oxidation and pyrolysis was developed on the basis of a calculated propagation velocity according to empirical correlation.The model was validated against the smoldering experiment of biomass rods under natural conditions,and the maximum error was smaller than 31%.Parameter sensitivity analysis found that propagation velocity decreases significantly while oxidation area and pyrolysis zone increase significantly with the increasing diameter of rod fuel.
基金Project (2013CB035401) supported by the National Basic Research Program of ChinaProject (51475478) supported by the National Natural Science Foundation of ChinaProject (2012AA041803) supported by the High-Tech Research and Development Program of China
文摘The center cutter of a hard rock tunnel boring machine(TBM) is installed on the cutterhead at a small radius and thus bears complex side force.Given this fact,the formation mechanism and change law of the side force suffered by the center cutter were studied.Based on the rock shear failure criterion in combination with the lateral rolling width,a model for predicting the average side force was set up.Besides,a numerical analysis model of the rock fragmentation of the center cutter was established,and the instantaneous load changing features were investigated.Results shows that the inner side of the center cutter can form lateral rolling annulus in rock during the rotary cutting process.The smaller the installation radius is,the greater the cutter side force will be.In a working condition,the side force of the innermost center cutter is 11.66 k N,while it decreases sharply when installation radius increases.Variation tends to be gentle when installation radius is larger than 500 mm,and the side force of the outermost center cutter is reduced to 0.74 k N.
文摘Wellbore stability analysis is a growing concern in oil industries. There are many parameters affecting the stability of a wellbore including geomechanical properties (e.g., elastic modulus, uni-axial compressive strength (UCS) and cohesion) and acting forces (e.g., field stresses and mud pressure). Accurate determination of these parameters is time-consuming, expensive and sometimes even impossible. This work offers a systematic sensitivity analysis to quantify the amount of each parameter’s effect on the stability of a wellbore. Maximum wellbore wall displacement is used as a stability factor to study the stability of a wellbore. A 3D finite difference method with Mohr model is used for the numerical modeling. The numerical model is verified against an analytical solution. A dimensionless sensitivity factor is developed in order to compare the results of various parameters in the sensitivity analysis. The results show a different order of importance of parameters based on rock strength. The most sensitive properties for a weak rock are the maximum horizontal stress, internal friction angle and formation pressure, respectively, while for a strong rock, the most sensitive parameters are the maximum horizontal stress, mud pressure and pore pressure, respectively. The amount of error in wellbore stability analysis inflicted by the error in estimation of each parameter was also derived.
基金Project(2010CB732101) supported by the National Basic Research Program of ChinaProject(51079145) supported by the National Natural Science Foundation of China
文摘A coupled thermo-hydro-mechanical-migratory model of dual-porosity medium for saturated-unsaturated ubiquitous-joint rockmass was established,in which the stress field and the temperature field were single,but the seepage field and the concentration field were double,and the influences of sets,spaces,angles,continuity ratios,stiffnesses of fractures on the constitutive relationship of the medium were considered.Also,the relative two-dimensional program of finite element method was developed.Taking a hypothetical nuclear waste repository as a calculation example,the case in which the rockmass was unsaturated dual-porosity medium and radioactive nuclide leak was simulated numerically,and the temperatures,negative pore pressures,saturations,flow velocities,nuclide concentrations and principal stresses in the rockmass were investigated.The results show that the negative pore pressures and nuclide concentrations in the porosity and fracture present different changes and distributions.Even though the saturation degree in porosity is only about 1/10 that in fracture,the flow velocity of underground water in fracture is about three times that in porosity because the permeability coefficient of fracture is almost four orders higher than that of porosity.The value of nuclide concentration in fracture is close to that in porosity.
基金Project(2012BAK09B02-05)supported by the National"Twelfth Five"Science and Technology Support Program,ChinaProject(51274250)supported by the National Natural Science Foundation of China+2 种基金Project(2013zzts057)supported by the Fundamental Research Funds for the Central Universities,ChinaProject(11KF02)supported by the Research Fund of the State Key Laboratory of Coal Resources and Mine safety,CUMT,ChinaProject(2012M511417)supported by China Postdoctoral Science Foundation
文摘There were differences between real boundary and blast hole controlling boundary of irregular mined-out area in underground metal mines. There were errors in numerical analysis of stability for goaf, if it was analyzed as regular 3D mined-out area and the influence of coupling stress-seepage-disturbance was not considered adequately. Taking a lead zinc mine as the background, the model was built by the coupling of Surpac and Midas-Gts based on the goaf model precisely measured by CMS.According to seepage stress fundamental equations based on the equivalent continuum mechanical and the theory about equivalent load of dynamic disturbance in deep-hole blasting, the stability of mined-out area under multi-field coupling of stress-seepage-dynamic disturbance was numerically analyzed. The results show that it is more consistent between the numerical analysis model based on the real model of irregular 3D shape goaf and the real situation, which could faithfully reappear the change rule of stress–strain about the surrounding rock under synthetic action of blasting dynamic loading and the seepage pressure. The mined-out area multi-field coupling formed by blasting excavation is stable. Based on combination of the advantages of the CMS,Surpac and Midas-Gts, and fully consideration of the effects of multi-field coupling, the accurate and effective way could be provided for numerical analysis of stability for mined-out area.
基金Project(51808116) supported by the National Natural Science Foundation of ChinaProject(BK20180404) supported by the Natural Science Foundation of Jiangsu Province, China+1 种基金Project(KFJ170106) supported by the Changsha University of Science & Technology via Open Fund of National Engineering Laboratory of Highway Maintenance Technology,ChinaProject(242020R40133) supported by Fundamental Research Funds for the Central Universities, China。
文摘In order to study the safety factor and instability process of cohesive soil slope, the discrete element method(DEM) was applied. DEM software PFC2 D was used to simulate the triaxial test to study the influence of the particle micro parameters on the macroscopic characteristics of cohesive soil and calibrate the micro parameters of DEM model on this basis. Embankment slope stability analysis was carried out by strength reduction and gravity increase method, it is shown that the safety factor obtained by strength reduction method is more conservative, and the arc-shaped feature of the sliding surface under the gravity increase method is more obvious. Throughout the progressive failure process, the failure trends, maximum displacements, and velocity changes obtained by the two methods were consistent. When slope was destroyed, the upper part was cracked, the middle part was sheared, and the lower part was destroyed by extrusion. The conclusions of this paper can be applied to the safety factor calculation of cohesive soil slopes and the analysis of the instability process.
文摘The distributed shear stress and the displacement across shear band, the evolution of plastic zones, and the load-carrying capacity of rock specimen were investigated in plane strain direct shear test according to Fast Lagrangian Analysis of Continua （FLAC）. And then the shear displacement distribution in normal direction of system composed of localized shear band and elastic rock was analyzed based on gradient-dependent plasticity. The adopted failure criterion was a composite of Mohr-Coulomb criterion, that is, the relation between tension cut-off and postpeak constitutive of rock was linear strain-softening. Numerical results show that shear stress field approximately undergoes three different stages. At first, shear stress is only concentrated in the middle of top and base of specimen. Next, shear stress in the middle of specimen tends to increase, owing to superposition of shear stresses. Interestingly, two peaks of shear stress appear far from the loading ends of specimen, and the peaks approach with the increase in timestep until elements at the center of specimen yield. Finally, relatively lower shear stress level is reached in large part of specimen except in the regions near the two ends. As flow stress decreases, the analytical shear displacement distribution in shear band based on gradient-dependent plasticity becomes steeps outside the band, it is linear and its slope tends to decrease. These theoretical results qualitatively agree with that of the present numerical predicted results. Main advantage of the analytical solution over the numerical results according to FLAC is that it is continuous, smooth and non-linear （except at elastic stage）.
基金Project(52175361)supported by the National Natural Science Foundation of ChinaProject(2019 CFA 041)supported by the Natural Science Foundation of Hubei Province,ChinaProject(WUT:202407002)supported by the Fundamental Research Funds for the Central Universities,China。
文摘A novel spiral non-circular bevel gear that could be applied to variable-speed driving in intersecting axes was proposed by combining the design principles of non-circular bevel gears and the manufacturing principles of face-milling spiral bevel gears.Unlike straight non-circular bevel gears,spiral non-circular bevel gears have numerous advantages,such as a high contact ratio,high intensity,good dynamic performance,and an adjustable contact region.In addition,while manufacturing straight non-circular bevel gears is difficult,spiral non-circular bevel gears can be efficiently and precisely fabricated with a 6-axis bevel gear cutting machine.First,the generating principles of spiral non-circular bevel gears were introduced.Next,a mathematical model,including a generating tooth profile,tooth spiral,pressure angle,and generated tooth profile for this gear type was established.Then the precision of the model was verified by a tooth contact analysis using FEA,and the contact patterns and stress distributions of the spiral non-circular bevel gears were investigated.
基金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.
基金Project(2007BAF02B12)supported by the National Science Technology Support Program of ChinaProjects(E2011203090,E2012203028)supported by the Natural Science Foundation of Hebei Province,China
文摘Aiming at accuracy control of the thermal crown of work rolls in cold rolling,new parameters such as regulation domain and control-efficiency factors were proposed and a numerical analysis model of the thermal crown of work rolls was established using finite difference method to study roll's thermal deformation.Based on simulation results,the influences of control-efficiency factors on thermal crown are presented and the thermal crown of work rolls is analyzed after taking sub-cooling of sprinkling beam into consideration.It has been found that the control-efficiency factor of any position on the roll's surface is linear function of the temperature and the control ability of water temperature is stronger than other control parameters.In addition,the verification of the model has been carried out based on the producing technology data in some factories and the numerical simulation results coincide well with the experimental data.Therefore,this work has important value for on-line control of roll's crown in cold rolling.
基金Project(2011ZX05009-004)supported by the National Science and Technology Major Projects of China
文摘Dimensional analysis and numerical simulations were carried out to research prediction method of breakthrough time of horizontal wells in bottom water reservoir. Four dimensionless independent variables and dimensionless time were derived from 10 influencing factors of the problem by using dimensional analysis. Simulations of horizontal well in reservoir with bottom water were run to find the prediction correlation. A general and concise functional relationship for predicting breakthrough time was established based on simulation results and theoretical analysis. The breakthrough time of one conceptual model predicted by the correlation is very close to the result by Eclipse with less than 2% error. The practical breakthrough time of one well in Helder oilfield is 10 d, and the predicted results by the method is 11.2 d, which is more accurate than the analytical result. Case study indicates that the method could predict breakthrough time of horizontal well under different reservoir conditions accurately. For its university and ease of use, the method is suitable for quick prediction of breakthrough time.
基金Projects(cstc2012jjA0510,cstc2013jcyjA30014)supported by Chongqing Natural Science Foundation in ChinaProject(CDJZR12200011)supported by the Fundamental Research Funds for the Central Universities in China+1 种基金Project(KJTD201305)supported by the Innovation Team Building Programs of Chongqing Universities in ChinaProject supported by the Scientific Research Foundation for the Returned Oversea Chinese Scholars
文摘A footing may get an eccentric load caused by earthquake or wind, thus the bearing capacity of footing subjected to eccentric load become a fundamental geotechnical problem. The conventional limit equilibrium method used for this problem usually evaluates the material properties only by its final strength. But the classical finite element method(FEM) does not necessarily provide a clear collapse mechanism associated with the yield condition of elements. To overcome these defects, a numerical procedure is proposed to create an explicit collapse mode combining a modified smeared shear band approach with a modified initial stress method. To understand the practical performance of sand foundation and verify the performance of the proposed procedure applied to the practical problems, the computing results were compared with the laboratory model tests results and some conventional solutions. Furthermore, because the proposed numerical procedure employs a simple elasto-plastic model which requires a small number of soil parameters, it may be applied directly to practical design works.
基金Project(2015B37714)supported by the Fundamental Research Funds for the Central Universities of ChinaProject(51605005)supported by the National Natural Science Foundation of China+1 种基金Project(ZK16-03-03)supported by the Open Foundation of Jiangsu Wind Technology Center,ChinaProject([2013]56)supported by the First Group of 2011 Plan of Jiangsu Province,China
文摘The parameter sensitivities affecting the flutter speed of the NREL （National Renewable Energy Laboratory） 5-MW baseline HAWT （horizontal axis wind turbine） blades are analyzed. An aeroelastic model, which comprises an aerodynamic part to calculate the aerodynamic loads and a structural part to determine the structural dynamic responses, is established to describe the classical flutter of the blades. For the aerodynamic part, Theodorsen unsteady aerodynamics model is used. For the structural part, Lagrange’s equation is employed. The flutter speed is determined by introducing “V–g” method to the aeroelastic model, which converts the issue of classical flutter speed determination into an eigenvalue problem. Furthermore, the time domain aeroelastic response of the wind turbine blade section is obtained with employing Runge-Kutta method. The results show that four cases （i.e., reducing the blade torsional stiffness, moving the center of gravity or the elastic axis towards the trailing edge of the section, and placing the turbine in high air density area） will decrease the flutter speed. Therefore, the judicious selection of the four parameters （the torsional stiffness, the chordwise position of the center of gravity, the elastic axis position and air density） can increase the relative inflow speed at the blade section associated with the onset of flutter.
基金Project(2006BAF04B06) supported by the National Key Technology R & D Program of ChinaProject(2005AA101B19) supported by the Key Technology R & D Program of Hubei Province, China
文摘The finite element analysis （FEA） software Ansys was employed to study the stress state of the dies of both plane and non-plane parting face structures with uniform interference and the die of plane parting face structure with non-uniform interference. Considering the symmetry of the die, a half gear tooth model of the two-ring assembled die with 2.5 GPa inner pressure was constructed. Four paths were defined to investigate the stress state at the bottom comer of the die where stress concentration was serious. FEA results show that the change of parting face from non-plane to plane can greatly reduce the stress at the teeth tips of the die so that the tip fracture is avoided. The interference structure of the die is the most important influencing factor for the stress concentration at the bottom comer. When non-uniform interference is adopted the first principal stress at the comer on the defined paths of the die is much lower than that with uniform interference. The bottom hole radius is another important influencing factor for the comer stress concentration. The first principal stress at the comer of the plane parting face die with non-uniform interference is reduced from 2.3 to 1.9 GPa when the hole radius increases from 12.5 to 16.0 mm. The optimization of the die structure increases the life of the die from 100 to 6 000 hits.
