By solving steady model of air flow diffusion and chemical reaction in loose coal, distribution of oxygen concentration and flow velocity magnitude were obtained. Compared the simulating results with critic value as w...By solving steady model of air flow diffusion and chemical reaction in loose coal, distribution of oxygen concentration and flow velocity magnitude were obtained. Compared the simulating results with critic value as well as duration of spontaneous combustion from large-scale spontaneous combustion experiment, 'three zones' of spontaneous combustion were partitioned and mining conditions to avoid spontaneous combustion were obtained. The above method was employed to partition 'three zones' in gob of fully mechanized top-coal caving long wall face and got fairly good result. Calculation of the above method is much smaller than simulating the whole process of coal spontaneous combustion, but the prediction precision can satisfy the demand of predicting and extinguishing spontaneous combustion in mining.展开更多
It is estimated that the aggregate consumption behavior of families by computational procedures in an abstract model contemplates the maximization of agent utility in each period （two dates）, the attendance of budge...It is estimated that the aggregate consumption behavior of families by computational procedures in an abstract model contemplates the maximization of agent utility in each period （two dates）, the attendance of budgetary restrictions and the conditions of general equilibrium （without production）. It uses the technique of selection of candidate points in a simulation process with a portfolio of efficient assets and a hypothesis for the process of determining the returns of the securities： A GARCH process. By this technique, it compares the stochastic volatility patterns between the artificial series, obtained in the simulation, and the real series of household aggregate consumption in the US and Brazil.展开更多
Considering the level distribution of soil layers, the soils surrounding pile are simulated with level finite layer elements. Supposing that the vertical deformation of the soil elements surrounding pile varies in the...Considering the level distribution of soil layers, the soils surrounding pile are simulated with level finite layer elements. Supposing that the vertical deformation of the soil elements surrounding pile varies in the form of exponent function with radial distance, and considering the nonlinear constitutive relation of stress and strain, the stiffness matrix is established. The mechanics behavior of the pile—soil interface is simulated with a nonlinear interface element. This method can truly express the behavior of the pile-soil system. The load-settlement relation Q-S curves of two big diameter prototype piles on bearing test are analyzed, and satisfying results are obtained. This method is reasonable in theory and feasible in engineering.展开更多
Wire arc additive manufacturing(WAAM)has emerged as a promising approach for fabricating large-scale components.However,conventional WAAM still faces challenges in optimizing microstructural evolution,minimizing addit...Wire arc additive manufacturing(WAAM)has emerged as a promising approach for fabricating large-scale components.However,conventional WAAM still faces challenges in optimizing microstructural evolution,minimizing additive-induced defects,and alleviating residual stress and deformation,all of which are critical for enhancing the mechanical performance of the manufactured parts.Integrating interlayer friction stir processing(FSP)into WAAM significantly enhances the quality of deposited materials.However,numerical simulation research focusing on elucidating the associated thermomechanical coupling mechanisms remains insufficient.A comprehensive numerical model was developed to simulate the thermomechanical coupling behavior in friction stir-assisted WAAM.The influence of post-deposition FSP on the coupled thermomechanical response of the WAAM process was analyzed quantitatively.Moreover,the residual stress distribution and deformation behavior under both single-layer and multilayer deposition conditions were investigated.Thermal analysis of different deposition layers in WAAM and friction stir-assisted WAAM was conducted.Results show that subsequent layer deposition induces partial remelting of the previously solidified layer,whereas FSP does not cause such remelting.Furthermore,thermal stress and deformation analysis confirm that interlayer FSP effectively mitigates residual stresses and distortion in WAAM components,thereby improving their structural integrity and mechanical properties.展开更多
Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The t...Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The traditional thermal elastic-plastic finite element method(TEP-FEM)can accurately predict welding deformation.However,its efficiency is low because of the complex nonlinear transient computation,making it difficult to meet the needs of rapid engineering evaluation.To address this challenge,this study proposes an efficient prediction method for welding deformation in marine thin plate butt welds.This method is based on the coupled temperature gradient-thermal strain method(TG-TSM)that integrates inherent strain theory with a shell element finite element model.The proposed method first extracts the distribution pattern and characteristic value of welding-induced inherent strain through TEP-FEM analysis.This strain is then converted into the equivalent thermal load applied to the shell element model for rapid computation.The proposed method-particularly,the gradual temperature gradient-thermal strain method(GTG-TSM)-achieved improved computational efficiency and consistent precision.Furthermore,the proposed method required much less computation time than the traditional TEP-FEM.Thus,this study lays the foundation for future prediction of welding deformation in more complex marine thin plates.展开更多
Bedding parallel stepped rock slopes exist widely in nature and are used in slope engineering.They are characterized by complex topography and geological structure and are vulnerable to shattering under strong earthqu...Bedding parallel stepped rock slopes exist widely in nature and are used in slope engineering.They are characterized by complex topography and geological structure and are vulnerable to shattering under strong earthquakes.However,no previous studies have assessed the mechanisms underlying seismic failure in rock slopes.In this study,large-scale shaking table tests and numerical simulations were conducted to delineate the seismic failure mechanism in terms of acceleration,displacement,and earth pressure responses combined with shattering failure phenomena.The results reveal that acceleration response mutations usually occur within weak interlayers owing to their inferior performance,and these mutations may transform into potential sliding surfaces,thereby intensifying the nonlinear seismic response characteristics.Cumulative permanent displacements at the internal corners of the berms can induce quasi-rigid displacements at the external corners,leading to greater permanent displacements at the internal corners.Therefore,the internal corners are identified as the most susceptible parts of the slope.In addition,the concept of baseline offset was utilized to explain the mechanism of earth pressure responses,and the result indicates that residual earth pressures at the internal corners play a dominant role in causing deformation or shattering damage.Four evolutionary deformation phases characterize the processes of seismic responses and shattering failure of the bedding parallel stepped rock slope,i.e.the formation of tensile cracks at the internal corners of the berm,expansion of tensile cracks and bedding surface dislocation,development of vertical tensile cracks at the rear edge,and rock mass slipping leading to slope instability.Overall,this study provides a scientific basis for the seismic design of engineering slopes and offers valuable insights for further studies on preventing seismic disasters in bedding parallel stepped rock slopes.展开更多
The graded density impactor(GDI)dynamic loading technique is crucial for acquiring the dynamic physical property parameters of materials used in weapons.The accuracy and timeliness of GDI structural design are key to ...The graded density impactor(GDI)dynamic loading technique is crucial for acquiring the dynamic physical property parameters of materials used in weapons.The accuracy and timeliness of GDI structural design are key to achieving controllable stress-strain rate loading.In this study,we have,for the first time,combined one-dimensional fluid computational software with machine learning methods.We first elucidated the mechanisms by which GDI structures control stress and strain rates.Subsequently,we constructed a machine learning model to create a structure-property response surface.The results show that altering the loading velocity and interlayer thickness has a pronounced regulatory effect on stress and strain rates.In contrast,the impedance distribution index and target thickness have less significant effects on stress regulation,although there is a matching relationship between target thickness and interlayer thickness.Compared with traditional design methods,the machine learning approach offers a10^(4)—10^(5)times increase in efficiency and the potential to achieve a global optimum,holding promise for guiding the design of GDI.展开更多
Hydrogen displays the potential to partially replace pulverized coal injection(PCI)in the blast furnace,and it can reduce CO_(2)emissions.In this paper,a three-dimensional mathematical model of hydrogen and pulverized...Hydrogen displays the potential to partially replace pulverized coal injection(PCI)in the blast furnace,and it can reduce CO_(2)emissions.In this paper,a three-dimensional mathematical model of hydrogen and pulverized coal co-injection in blast furnace tuyere was established through numerical simulation,and the effect of hydrogen injection and oxygen enrichment interaction on pulverized coal combustion and raceway smelting was investigated.The simulation results indicate that when the coal injection rate decreased from 36 to 30t/h and the hydrogen injection increased from 0 to 3600 m^(3)/h,the CO_(2)emissions decreased from 1860 to 1551 kg/t,which represents a16.6%reduction,and the pulverized coal burnout decreased from 70.1%to 63.7%.The heat released from hydrogen combustion can not only promote the volatilization of pulverized coal but also affect the combustion reaction between volatilization and oxygen,which resulted in a decrease in the temperature at the end of the raceway.Co-injection of hydrogen with PCI increased the wall temperature near the upper half part of the raceway and at the outlet of the tuyere,which required a high cooling efficiency to extend the service life of the blast furnace.The increase in oxygen level compensated for the decreased average temperature in the raceway due to hydrogen injection.The increase in the oxygen content by 3%while maintaining constant hydrogen and PCI injection rates increased the burnout and average raceway temperature by 4.2%and 43 K,respectively.The mole fraction of CO and H_(2) production increased by 0.04 and 0.02,respectively.Burnout can be improved through optimization of the particle size distribution of pulverized coal.展开更多
To further reduce the explosive thickness and to improve the bonding quality of titanium/steel composite plates,explosive welding experiments of TA1/Q235 were conducted using a low detonation velocity explosive(53#)un...To further reduce the explosive thickness and to improve the bonding quality of titanium/steel composite plates,explosive welding experiments of TA1/Q235 were conducted using a low detonation velocity explosive(53#)under the guidance of the explosive welding lower limit principle with the flyer plate thicknesses of 1,2,and 4 mm and gaps of 3,6,and 8 mm.The weldability window for titanium/steel explosive welding was calculated,and a quantitative relationship between dynamic and static process parameters was developed.Aβ-V_(p) high-speed inclined collision model was proposed,and two-dimensional numerical simulations for the explosive welding tests were performed using the smoothed particle hydrodynamics(SPH)algorithm,revealing the growth evolution mechanisms of the typical waveform morphology characteristics.Through microstructural characterization techniques,such as optical microscope,scanning electron microscope,energy dispersive spectrometer,and electron backscattered diffractometer,and mechanical property tests in terms of shear strength,bending performance,and impact toughness,the microstructure and mechanical properties of the interfaces of explosively welded TA1/Q235 composite plates were investigated.Results show that the quality of interface bonding is excellent,presenting typical waveform morphology with an average interface shear strength above 330 MPa and an average impact toughness exceeding 81 J.All samples can be bent to 180°without significant delamination or cracking defects.展开更多
A suction casting experiment was conducted on Zr_(55)Cu_(30)Al_(10)Ni_(5)(at%)amorphous alloy.Using ProCAST software,numerical simulations were performed to analyze the filling and solidification processes.The velocit...A suction casting experiment was conducted on Zr_(55)Cu_(30)Al_(10)Ni_(5)(at%)amorphous alloy.Using ProCAST software,numerical simulations were performed to analyze the filling and solidification processes.The velocity field during the filling process and the temperature field during the solidification process of the alloy melt under different process parameters were obtained.Based on the simulation results,a Zr-based amorphous alloy micro-gear was prepared via casting.The results indicate that increasing the suction casting temperature enhances the fluidity of alloy melt but induces unstable flow rate during filling,which is detrimental to complete filling.Zr-based amorphous micro-gears with a module of 0.6 mm,a tooth top diameter of 8 mm,and 10 teeth were prepared through the suction casting.X-ray diffraction and differential scanning calorimetry analyses confirm that the fabricated micro-gear exhibits characteristic amorphous structural features,demonstrating well-defined geometrical contours and satisfactory forming completeness.展开更多
To elucidate the influence of confining pressure on microcrack evolution and macroscopic failure mechanisms in granite,a multi-perspective approach was adopted.This approach combined triaxial compression tests,acousti...To elucidate the influence of confining pressure on microcrack evolution and macroscopic failure mechanisms in granite,a multi-perspective approach was adopted.This approach combined triaxial compression tests,acoustic emission(AE)monitoring,and PFC simulations.The results show that:1)Confining pressure exhibits a pronounced linear correlation with both yield strength and compressive strength.The enhancement of confining pressure significantly improves the deformability of granite and promotes a progressive shift in failure mechanism from brittle rupture to ductile deformation;2)Increasing confining pressure elevates the stress threshold for microcrack initiation and suppresses crack propagation.As a result,the proportion of shear cracks increases(based on AE analysis)from 18.71%to 61.2%,marking a transition in the dominant failure mode from tensile to shear;3)Confining pressure facilitates the development of grain boundary shear cracks(GBSCs),establishing the primary damage pathways.In addition,local stress concentrations under high confinement conditions trigger intragranular cracking.This highlights the regulatory effect of confining pressure on microcrack evolution.展开更多
Natural gas hydrate widely exists in the South China Sea as clean energy.A three-phase transition layer widely exists in low permeability Class I hydrates in the Shenhu offshore area.Therefore,taking into account the ...Natural gas hydrate widely exists in the South China Sea as clean energy.A three-phase transition layer widely exists in low permeability Class I hydrates in the Shenhu offshore area.Therefore,taking into account the low-permeability characteristics with an average permeability of 5.5 mD and moderate heterogeneity,a 3-D geological model of heterogeneous Class I hydrate reservoirs with three-phase transition layers is established by Kriging interpolation and stochastic modeling method,and a numerical simulation model is used to describe the depressurization production performance of the reservoir.With the development of depressurization,a specific range of complete decomposition zones appear both in the hydrate and transition layers.The entire decomposition zone of the whole reservoir tends to outward and upward diffusion.There is apparent methane escape in the three-phase transition layer.Due to the improvement of local permeability caused by the phase transition of hydrate dissociation,some methane accumulation occurs at the bottom of the hydrate layer,forming a local methane enrichment zone.The methane migration trends in reservoirs are mainly characterized by movement toward production wells and hydrate layers under the influence of gravity.However,due to the permeability limitation of hydrate reservoirs,many fluids have not been effectively produced and remain in the reservoir.Therefore,to improve the effective pressure drop of the reservoir,the perforation method and pressure reduction method were optimized by analyzing the influencing factors based on the gas production rate.The comparative study demonstrates that perforating through the free gas layer combined with one-time depressurization can enhance the effective depressurization and improve production performance.The gas production rate from perforating through the free gas layer can be twice as high as that from perforating through the transition layer.This study can provide theoretical support for the utilization of marine energy.展开更多
Heat treatment processes, such as annealing and quenching, are crucial in determining residual stress evolution, microstructural changes and mechanical properties of metallic materials, with residual stresses playing ...Heat treatment processes, such as annealing and quenching, are crucial in determining residual stress evolution, microstructural changes and mechanical properties of metallic materials, with residual stresses playing a greater role in the performance of components. This paper investigates the effect of heat treatment on residual stresses induced in AISI 1025, manufactured using LENS. Finite element model was developed and simulated to analyze residual stress development. AISI 1025 samples suitable for tool and die applications in Fused Deposition Modelling (FDM) filament production, were fabricated using Laser Engineered Net Shaping (LENS) process, followed by heat treatment where annealing and quenching processes were done. The material’s microstructure, residual stress and hardness of heat-treated samples under investigation, were compared against the as-built samples. The results indicated that after annealing, tensile residual stresses were reduced by 93%, resulting in a reduced crack growth rate, compared to the as-built sample, although the hardness was reduced significantly by 25%. On the other hand, high tensile residual stresses of 425 ± 14 MPa were recorded after quenching process with an improvement of hardness by 21%.展开更多
Hypersonic magnetohydrodynamic(MHD)control effectively enhances the aerothermal environment of aerospace vehicles,demonstrating considerable potential in plasma flow regulation and aerodynamic optimiza-tion.As aerospa...Hypersonic magnetohydrodynamic(MHD)control effectively enhances the aerothermal environment of aerospace vehicles,demonstrating considerable potential in plasma flow regulation and aerodynamic optimiza-tion.As aerospace vehicles progress toward mid-low-altitude hypersonic regimes,their external aerothermal conditions become increasingly severe.This study addresses the challenges of complex aerodynamic force/heat environments and the difficulties in MHD control numerical simulations for hypersonic vehicles at mid-low al-titudes.On the basis of the perfect gas model and the low magnetic Reynolds number assumption,we conduct numerical simulations of MHD control under mid-low altitudes,high-Mach-number conditions.The findings reveal the following:(1)the low magnetic Reynolds number assumption is valid and computationally accurate,as corroborated by a comparative analysis with the literature;(2)in the mid-low altitude hypersonic regime,magnetic fields significantly suppress the shock standoffdistance and reduce the surface heat flux.Both the mag-netically controlled shock wave and the thermal protection exhibit nonlinear variations with the Mach number,increasing and then decreasing as the Mach number increases.The optimal Mach number for shock wave control is 13,whereas optimal thermal protection is achieved at Mach 15.At an altitude of 40 km,the optimal magne-tohydrodynamic Mach range spans 13-17,achieving a maximum heat flux attenuation of 28.81%.Additionally,the effects of magnetic shock wave control correlate approximately exponentially with altitude within certain parameters,whereas the efficacy of thermal protection behaves linearly with altitude variations.展开更多
After a long period of water flooding development,the oilfield has entered the middle and high water cut stage.The physical properties of reservoirs are changed by water erosion,which directly impacts reservoir develo...After a long period of water flooding development,the oilfield has entered the middle and high water cut stage.The physical properties of reservoirs are changed by water erosion,which directly impacts reservoir development.Conventional numerical reservoir simulation methodologies typically employ static assumptions for model construction,presuming invariant reservoir geological parameters throughout the development process while neglecting the reservoir’s temporal evolution characteristics.Although such simplifications reduce computational complexity,they introduce substantial descriptive inaccuracies.Therefore,this paper proposes a meshless numerical simulation method for reservoirs that considers time-varying characteristics.This method avoids the meshing in traditional numerical simulation methods.From the fluid flow perspective,the reservoir’s computational domain is discretized into a series of connection units.An influence domain with a certain radius centered on the nodes is selected,and one-dimensional connection units are established between the nodes to achieve the characterization of the flow topology structure of the reservoir.In order to reflect the dynamic evolution of the reservoir’s physical properties during the water injection development process,the time-varying characteristics are incorporated into the formula of the seepage characteristic parameters in the meshless calculation.The change relationship of the permeability under different surface fluxes is considered to update the calculated connection conductivity in real time.By combining with the seepage control equation for solution,a time-varying meshless numerical simulation method is formed.The results show that compared with the numerical simulationmethod of the connection elementmethod(CEM)that only considers static parameters,this method has higher simulation accuracy and can better simulate the real migration and distribution of oil and water in the reservoir.Thismethod improves the accuracy of reservoir numerical simulation and the development effect of oilfields,providing a scientific basis for optimizing the water injection strategy,adjusting the production plan,and extending the effective production cycle of the oilfield.展开更多
To mitigate the sand burial of highways in sandy regions,a separated subgrade design was widely adopted in the embankments of high-grade highways,but the problem of sand deposition on subgrade slopes and pavements sti...To mitigate the sand burial of highways in sandy regions,a separated subgrade design was widely adopted in the embankments of high-grade highways,but the problem of sand deposition on subgrade slopes and pavements still happens frequently.Based on the theory of wind-sand two-phase flow,this paper constructed a three-dimensional model of the separated subgrade,the wind-sand flow transport law around the subgrade with varying median strip widths and concave depths was simulated by Fluent software.After comparison and analysis of seven subgrade models,the flow field distribution,wind speed horizontal variation,and erosion-deposition characteristics were investigated.The findings are as follows:(1)The width of the median strip in the separated subgrade had significant influences on the wind-sand flow.The smooth passage of wind-sand flow over the road surface was facilitated with the increase of the median strip width.However,sand deposition in the median strip happened.It can lead to secondary sand damage of downwind subgrade and increase the work load of road sand removal for subsequent maintenance.(2)The obstruction to airflow and sand accumulation was aggravated with greater concave depth of the median strip.Therefore,it is advisable to minimize the concave depth of the median strip in case of more sand damage.(3)A median strip width exceeding 12 m(possibly without guardrails)for an integral embankment without enough road land is recommended.Conversely,median strip width of over 40 m for separate subgrade with unrestricted land is suggested.(4)In the case of sand deposition in the existing separated subgrade,the median strip can be filled by sand deposition or other materials,then was covered with gravel to form a flat ground like Gobi smooth surface,which can let the wind-blown sand flow pass through the subgrade section without sand deposition.展开更多
The study aimed to address the issue of elevated residual stress levels in dissimilar girth welds of cast steel joints.To achieve this,the hybrid welding technology,which yields high welding speeds while simultaneousl...The study aimed to address the issue of elevated residual stress levels in dissimilar girth welds of cast steel joints.To achieve this,the hybrid welding technology,which yields high welding speeds while simultaneously reducing residual stresses,has been introduced.This study utilizes a numerical simulation method to investigate the temperature and residual stress field in the hybrid welding of G20Mn5 casting-Q355 low-alloy steel welded pipe.A com-parison of the findings of this study with those of other welding processes revealed the technological advantages of hybrid welding.The research outcomes show that due to geometric discontinuities and material differences,the temperature field of the joint exhibits uneven distribution characteristics,and the peak temperatures on the Q355 steel side exceeds those on the G20Mn5 steel side.An evident stress gra-dient is present in the residual stress field of the joint post-welding,with peak stress located at the weld root on the Q355 steel.Compared with arc welding,the hybrid welding leads to decreased residual stresses and deformation,with high stress outside the heat-affected zone diminishing rapidly.Furthermore,it significantly improves the welding efficiency.This study elucidates the distribution and underlying causes of thermal and residual stress fields in dissimilar girth welds.This serves as a foundation for the application of hybrid welding technology in welded cast steel joints.展开更多
To study the combustion behavior of municipal solid waste(MSW)in blast furnace,the combustion process within the raceway was simulated using computational fluid dynamics.Based on the parameters of an actual blast furn...To study the combustion behavior of municipal solid waste(MSW)in blast furnace,the combustion process within the raceway was simulated using computational fluid dynamics.Based on the parameters of an actual blast furnace,a three-dimensional model including coal lance,blowpipe,tuyere,and raceway was established.The model was then used to compare the combustion characteristics of pulverized coal and MSW in the raceway and to investigate the effects of blast temperature and particle size on the combustion characteristics of MSW in the raceway.The results showed that MSW combusted more rapidly,achieving a maximum temperature of 3839 K in the raceway,comparing to 2974 K during pulverized coal injection.However,the average temperature during MSW injection was 1790 K,which was 73 K lower than that of pulverized coal injection.The maximum velocity during MSW injection was 120 m/s,lower than 188 m/s obtained during pulverized coal injection.MSW could be completely burned out in the middle of the raceway,while the burnout of pulverized coal at this position was only 50%.The combustion effect of MSW makes no difference when the blast temperature increased from 1273 to 1673 K,due to its excellent combustion characteristic.When the MSW particle size was increased from 0.074 to 2 mm,the burnout was 75%,which was still higher than that of pulverized coal injection with a particle size of 0.074 mm.However,injecting larger-sized fuel might increase the risk of tuyere wear.To ensure stable furnace conditions and great combustion,a blast temperature of 1473 K and a MSW particle size of about 1 mm will be better.展开更多
Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 me...Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 mesh ceramic particles and quartz sand proppant clusters were elaborated using computed tomography(CT)scanning.These models were used to develop a numerical simulation framework based on the lattice Boltzmann method(LBM),enabling the investigation of gas-water flow behavior within proppant-filled fractures under varying driving forces and surface tensions.Simulation results at a closure pressure of 15 MPa have revealed that ceramic particles exhibit a simpler and more porous internal structure than quartz sand of the same size.Under identical flow conditions,ceramic proppants demonstrate higher fluid replacement efficiency.Replacement efficiency increases with higher porosity,greater driving force,and lower surface tension.Furthermore,fluid displacement is strongly influenced by pore geometry:flow is faster in straighter and wider channels,with preferential movement through larger pores forming dominant flow paths.The replacement velocity exhibits a characteristic time evolution,initially rapid,then gradually decreasing,correlating positively with the development of these dominant channels.展开更多
Deepwater drilling riser is the key equipment connecting the subsea wellhead and floating drilling platform.Due to complex marine environment,vortex-induced vibration(ViV)will be generated on riser,which will induce f...Deepwater drilling riser is the key equipment connecting the subsea wellhead and floating drilling platform.Due to complex marine environment,vortex-induced vibration(ViV)will be generated on riser,which will induce fatigue failure and even cause unpredictable drilling accidents.Therefore,it is important to study the ViV characteristics of deepwater drilling riser and reveal the main controlling factors for ensuring the safe and efficient operation of deepwater drilling engineering.In this paper,the ViV of deepwater drilling riser is numerically simulated in time domain based on the discrete vortex method(DvM).A hydrodynamic analysis model and governing equation of VIV is proposed with considering the effect of riser motion using DVM and slice method,where the governing equation is solved by Runge-Kutta method.Model validation is performed,which verified the correctness and accuracy of the mechanical model and the solution method.On this basis,the influence of the number of control points,current velocity,riser outer diameter,shear flow and top tension on the ViV characteristics of deepwater drilling risers are discussed in detail.The results show that with the increase of current velocity,the vibration amplitude of deepwater drilling riser decreases obviously,while the vibration frequency increases gradually.However,if the outer diameter of riser increases,the vibration amplitude increases,while the vibration frequency decreases gradually.The top tension also has great influence on the VIV of riser.When the top tension is 1.25 G,the VIV is suppressed to a certain extent.This study has guiding significance for optimal design and engineering control of deepwater drilling riser.展开更多
基金Supported by Natural Science Program of Shaanxi Province Education Department (05JK261)
文摘By solving steady model of air flow diffusion and chemical reaction in loose coal, distribution of oxygen concentration and flow velocity magnitude were obtained. Compared the simulating results with critic value as well as duration of spontaneous combustion from large-scale spontaneous combustion experiment, 'three zones' of spontaneous combustion were partitioned and mining conditions to avoid spontaneous combustion were obtained. The above method was employed to partition 'three zones' in gob of fully mechanized top-coal caving long wall face and got fairly good result. Calculation of the above method is much smaller than simulating the whole process of coal spontaneous combustion, but the prediction precision can satisfy the demand of predicting and extinguishing spontaneous combustion in mining.
文摘It is estimated that the aggregate consumption behavior of families by computational procedures in an abstract model contemplates the maximization of agent utility in each period （two dates）, the attendance of budgetary restrictions and the conditions of general equilibrium （without production）. It uses the technique of selection of candidate points in a simulation process with a portfolio of efficient assets and a hypothesis for the process of determining the returns of the securities： A GARCH process. By this technique, it compares the stochastic volatility patterns between the artificial series, obtained in the simulation, and the real series of household aggregate consumption in the US and Brazil.
文摘Considering the level distribution of soil layers, the soils surrounding pile are simulated with level finite layer elements. Supposing that the vertical deformation of the soil elements surrounding pile varies in the form of exponent function with radial distance, and considering the nonlinear constitutive relation of stress and strain, the stiffness matrix is established. The mechanics behavior of the pile—soil interface is simulated with a nonlinear interface element. This method can truly express the behavior of the pile-soil system. The load-settlement relation Q-S curves of two big diameter prototype piles on bearing test are analyzed, and satisfying results are obtained. This method is reasonable in theory and feasible in engineering.
基金National Key Research and Development Program of China(2022YFB4600902)Shandong Provincial Science Foundation for Outstanding Young Scholars(ZR2024YQ020)。
文摘Wire arc additive manufacturing(WAAM)has emerged as a promising approach for fabricating large-scale components.However,conventional WAAM still faces challenges in optimizing microstructural evolution,minimizing additive-induced defects,and alleviating residual stress and deformation,all of which are critical for enhancing the mechanical performance of the manufactured parts.Integrating interlayer friction stir processing(FSP)into WAAM significantly enhances the quality of deposited materials.However,numerical simulation research focusing on elucidating the associated thermomechanical coupling mechanisms remains insufficient.A comprehensive numerical model was developed to simulate the thermomechanical coupling behavior in friction stir-assisted WAAM.The influence of post-deposition FSP on the coupled thermomechanical response of the WAAM process was analyzed quantitatively.Moreover,the residual stress distribution and deformation behavior under both single-layer and multilayer deposition conditions were investigated.Thermal analysis of different deposition layers in WAAM and friction stir-assisted WAAM was conducted.Results show that subsequent layer deposition induces partial remelting of the previously solidified layer,whereas FSP does not cause such remelting.Furthermore,thermal stress and deformation analysis confirm that interlayer FSP effectively mitigates residual stresses and distortion in WAAM components,thereby improving their structural integrity and mechanical properties.
基金Supported by the National Natural Science Foundation of China under Grant No.51975138the High-Tech Ship Scientific Research Project from the Ministry of Industry and Information Technology under Grant No.CJ05N20the National Defense Basic Research Project under Grant No.JCKY2023604C006.
文摘Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The traditional thermal elastic-plastic finite element method(TEP-FEM)can accurately predict welding deformation.However,its efficiency is low because of the complex nonlinear transient computation,making it difficult to meet the needs of rapid engineering evaluation.To address this challenge,this study proposes an efficient prediction method for welding deformation in marine thin plate butt welds.This method is based on the coupled temperature gradient-thermal strain method(TG-TSM)that integrates inherent strain theory with a shell element finite element model.The proposed method first extracts the distribution pattern and characteristic value of welding-induced inherent strain through TEP-FEM analysis.This strain is then converted into the equivalent thermal load applied to the shell element model for rapid computation.The proposed method-particularly,the gradual temperature gradient-thermal strain method(GTG-TSM)-achieved improved computational efficiency and consistent precision.Furthermore,the proposed method required much less computation time than the traditional TEP-FEM.Thus,this study lays the foundation for future prediction of welding deformation in more complex marine thin plates.
基金supported by the National Natural Science Foundation of China (Grant No.52108361)the Sichuan Science and Technology Program of China (Grant No.2023YFS0436)the State Key Laboratory of Geohazard Prevention and Geoenvironment Protection Independent Research Project (Grant No.SKLGP2022Z015).
文摘Bedding parallel stepped rock slopes exist widely in nature and are used in slope engineering.They are characterized by complex topography and geological structure and are vulnerable to shattering under strong earthquakes.However,no previous studies have assessed the mechanisms underlying seismic failure in rock slopes.In this study,large-scale shaking table tests and numerical simulations were conducted to delineate the seismic failure mechanism in terms of acceleration,displacement,and earth pressure responses combined with shattering failure phenomena.The results reveal that acceleration response mutations usually occur within weak interlayers owing to their inferior performance,and these mutations may transform into potential sliding surfaces,thereby intensifying the nonlinear seismic response characteristics.Cumulative permanent displacements at the internal corners of the berms can induce quasi-rigid displacements at the external corners,leading to greater permanent displacements at the internal corners.Therefore,the internal corners are identified as the most susceptible parts of the slope.In addition,the concept of baseline offset was utilized to explain the mechanism of earth pressure responses,and the result indicates that residual earth pressures at the internal corners play a dominant role in causing deformation or shattering damage.Four evolutionary deformation phases characterize the processes of seismic responses and shattering failure of the bedding parallel stepped rock slope,i.e.the formation of tensile cracks at the internal corners of the berm,expansion of tensile cracks and bedding surface dislocation,development of vertical tensile cracks at the rear edge,and rock mass slipping leading to slope instability.Overall,this study provides a scientific basis for the seismic design of engineering slopes and offers valuable insights for further studies on preventing seismic disasters in bedding parallel stepped rock slopes.
基金supported by the Guangdong Major Project of Basic and Applied Basic Research(Grant No.2021B0301030001)the National Key Research and Development Program of China(Grant No.2021YFB3802300)the Foundation of National Key Laboratory of Shock Wave and Detonation Physics(Grant No.JCKYS2022212004)。
文摘The graded density impactor(GDI)dynamic loading technique is crucial for acquiring the dynamic physical property parameters of materials used in weapons.The accuracy and timeliness of GDI structural design are key to achieving controllable stress-strain rate loading.In this study,we have,for the first time,combined one-dimensional fluid computational software with machine learning methods.We first elucidated the mechanisms by which GDI structures control stress and strain rates.Subsequently,we constructed a machine learning model to create a structure-property response surface.The results show that altering the loading velocity and interlayer thickness has a pronounced regulatory effect on stress and strain rates.In contrast,the impedance distribution index and target thickness have less significant effects on stress regulation,although there is a matching relationship between target thickness and interlayer thickness.Compared with traditional design methods,the machine learning approach offers a10^(4)—10^(5)times increase in efficiency and the potential to achieve a global optimum,holding promise for guiding the design of GDI.
基金financially supported by the National Natural Science Foundation of China(No.51904026)the Fundamental Research Funds for the Central Universities(No.06500108)。
文摘Hydrogen displays the potential to partially replace pulverized coal injection(PCI)in the blast furnace,and it can reduce CO_(2)emissions.In this paper,a three-dimensional mathematical model of hydrogen and pulverized coal co-injection in blast furnace tuyere was established through numerical simulation,and the effect of hydrogen injection and oxygen enrichment interaction on pulverized coal combustion and raceway smelting was investigated.The simulation results indicate that when the coal injection rate decreased from 36 to 30t/h and the hydrogen injection increased from 0 to 3600 m^(3)/h,the CO_(2)emissions decreased from 1860 to 1551 kg/t,which represents a16.6%reduction,and the pulverized coal burnout decreased from 70.1%to 63.7%.The heat released from hydrogen combustion can not only promote the volatilization of pulverized coal but also affect the combustion reaction between volatilization and oxygen,which resulted in a decrease in the temperature at the end of the raceway.Co-injection of hydrogen with PCI increased the wall temperature near the upper half part of the raceway and at the outlet of the tuyere,which required a high cooling efficiency to extend the service life of the blast furnace.The increase in oxygen level compensated for the decreased average temperature in the raceway due to hydrogen injection.The increase in the oxygen content by 3%while maintaining constant hydrogen and PCI injection rates increased the burnout and average raceway temperature by 4.2%and 43 K,respectively.The mole fraction of CO and H_(2) production increased by 0.04 and 0.02,respectively.Burnout can be improved through optimization of the particle size distribution of pulverized coal.
基金Jiangsu Provincial Natural Science Foundation of China(BK20211232)2023 Major Science and Technology Projects of Nanjing City(202309011)。
文摘To further reduce the explosive thickness and to improve the bonding quality of titanium/steel composite plates,explosive welding experiments of TA1/Q235 were conducted using a low detonation velocity explosive(53#)under the guidance of the explosive welding lower limit principle with the flyer plate thicknesses of 1,2,and 4 mm and gaps of 3,6,and 8 mm.The weldability window for titanium/steel explosive welding was calculated,and a quantitative relationship between dynamic and static process parameters was developed.Aβ-V_(p) high-speed inclined collision model was proposed,and two-dimensional numerical simulations for the explosive welding tests were performed using the smoothed particle hydrodynamics(SPH)algorithm,revealing the growth evolution mechanisms of the typical waveform morphology characteristics.Through microstructural characterization techniques,such as optical microscope,scanning electron microscope,energy dispersive spectrometer,and electron backscattered diffractometer,and mechanical property tests in terms of shear strength,bending performance,and impact toughness,the microstructure and mechanical properties of the interfaces of explosively welded TA1/Q235 composite plates were investigated.Results show that the quality of interface bonding is excellent,presenting typical waveform morphology with an average interface shear strength above 330 MPa and an average impact toughness exceeding 81 J.All samples can be bent to 180°without significant delamination or cracking defects.
基金National Natural Science Foundation of China(51971103)Key Research and Development Program in Gansu Province(20YF8GA052)。
文摘A suction casting experiment was conducted on Zr_(55)Cu_(30)Al_(10)Ni_(5)(at%)amorphous alloy.Using ProCAST software,numerical simulations were performed to analyze the filling and solidification processes.The velocity field during the filling process and the temperature field during the solidification process of the alloy melt under different process parameters were obtained.Based on the simulation results,a Zr-based amorphous alloy micro-gear was prepared via casting.The results indicate that increasing the suction casting temperature enhances the fluidity of alloy melt but induces unstable flow rate during filling,which is detrimental to complete filling.Zr-based amorphous micro-gears with a module of 0.6 mm,a tooth top diameter of 8 mm,and 10 teeth were prepared through the suction casting.X-ray diffraction and differential scanning calorimetry analyses confirm that the fabricated micro-gear exhibits characteristic amorphous structural features,demonstrating well-defined geometrical contours and satisfactory forming completeness.
基金Projects(U23A2060,42177143)supported by the National Natural Science Foundation of China。
文摘To elucidate the influence of confining pressure on microcrack evolution and macroscopic failure mechanisms in granite,a multi-perspective approach was adopted.This approach combined triaxial compression tests,acoustic emission(AE)monitoring,and PFC simulations.The results show that:1)Confining pressure exhibits a pronounced linear correlation with both yield strength and compressive strength.The enhancement of confining pressure significantly improves the deformability of granite and promotes a progressive shift in failure mechanism from brittle rupture to ductile deformation;2)Increasing confining pressure elevates the stress threshold for microcrack initiation and suppresses crack propagation.As a result,the proportion of shear cracks increases(based on AE analysis)from 18.71%to 61.2%,marking a transition in the dominant failure mode from tensile to shear;3)Confining pressure facilitates the development of grain boundary shear cracks(GBSCs),establishing the primary damage pathways.In addition,local stress concentrations under high confinement conditions trigger intragranular cracking.This highlights the regulatory effect of confining pressure on microcrack evolution.
基金supported by the Sinopec Technology Research and Development Project(No.30000000-22-ZC0607-0235,No.33550000-22-ZC0607-0009)the National Natural Science Foundation of China(No.52334002).
文摘Natural gas hydrate widely exists in the South China Sea as clean energy.A three-phase transition layer widely exists in low permeability Class I hydrates in the Shenhu offshore area.Therefore,taking into account the low-permeability characteristics with an average permeability of 5.5 mD and moderate heterogeneity,a 3-D geological model of heterogeneous Class I hydrate reservoirs with three-phase transition layers is established by Kriging interpolation and stochastic modeling method,and a numerical simulation model is used to describe the depressurization production performance of the reservoir.With the development of depressurization,a specific range of complete decomposition zones appear both in the hydrate and transition layers.The entire decomposition zone of the whole reservoir tends to outward and upward diffusion.There is apparent methane escape in the three-phase transition layer.Due to the improvement of local permeability caused by the phase transition of hydrate dissociation,some methane accumulation occurs at the bottom of the hydrate layer,forming a local methane enrichment zone.The methane migration trends in reservoirs are mainly characterized by movement toward production wells and hydrate layers under the influence of gravity.However,due to the permeability limitation of hydrate reservoirs,many fluids have not been effectively produced and remain in the reservoir.Therefore,to improve the effective pressure drop of the reservoir,the perforation method and pressure reduction method were optimized by analyzing the influencing factors based on the gas production rate.The comparative study demonstrates that perforating through the free gas layer combined with one-time depressurization can enhance the effective depressurization and improve production performance.The gas production rate from perforating through the free gas layer can be twice as high as that from perforating through the transition layer.This study can provide theoretical support for the utilization of marine energy.
文摘Heat treatment processes, such as annealing and quenching, are crucial in determining residual stress evolution, microstructural changes and mechanical properties of metallic materials, with residual stresses playing a greater role in the performance of components. This paper investigates the effect of heat treatment on residual stresses induced in AISI 1025, manufactured using LENS. Finite element model was developed and simulated to analyze residual stress development. AISI 1025 samples suitable for tool and die applications in Fused Deposition Modelling (FDM) filament production, were fabricated using Laser Engineered Net Shaping (LENS) process, followed by heat treatment where annealing and quenching processes were done. The material’s microstructure, residual stress and hardness of heat-treated samples under investigation, were compared against the as-built samples. The results indicated that after annealing, tensile residual stresses were reduced by 93%, resulting in a reduced crack growth rate, compared to the as-built sample, although the hardness was reduced significantly by 25%. On the other hand, high tensile residual stresses of 425 ± 14 MPa were recorded after quenching process with an improvement of hardness by 21%.
基金the results of the research project funded by National Numerical Wind Tunnel Project of China.
文摘Hypersonic magnetohydrodynamic(MHD)control effectively enhances the aerothermal environment of aerospace vehicles,demonstrating considerable potential in plasma flow regulation and aerodynamic optimiza-tion.As aerospace vehicles progress toward mid-low-altitude hypersonic regimes,their external aerothermal conditions become increasingly severe.This study addresses the challenges of complex aerodynamic force/heat environments and the difficulties in MHD control numerical simulations for hypersonic vehicles at mid-low al-titudes.On the basis of the perfect gas model and the low magnetic Reynolds number assumption,we conduct numerical simulations of MHD control under mid-low altitudes,high-Mach-number conditions.The findings reveal the following:(1)the low magnetic Reynolds number assumption is valid and computationally accurate,as corroborated by a comparative analysis with the literature;(2)in the mid-low altitude hypersonic regime,magnetic fields significantly suppress the shock standoffdistance and reduce the surface heat flux.Both the mag-netically controlled shock wave and the thermal protection exhibit nonlinear variations with the Mach number,increasing and then decreasing as the Mach number increases.The optimal Mach number for shock wave control is 13,whereas optimal thermal protection is achieved at Mach 15.At an altitude of 40 km,the optimal magne-tohydrodynamic Mach range spans 13-17,achieving a maximum heat flux attenuation of 28.81%.Additionally,the effects of magnetic shock wave control correlate approximately exponentially with altitude within certain parameters,whereas the efficacy of thermal protection behaves linearly with altitude variations.
基金funded by the 14th Five-Year Plan Major Science and Technology Project of CNOOC project number KJGG2021-0506.
文摘After a long period of water flooding development,the oilfield has entered the middle and high water cut stage.The physical properties of reservoirs are changed by water erosion,which directly impacts reservoir development.Conventional numerical reservoir simulation methodologies typically employ static assumptions for model construction,presuming invariant reservoir geological parameters throughout the development process while neglecting the reservoir’s temporal evolution characteristics.Although such simplifications reduce computational complexity,they introduce substantial descriptive inaccuracies.Therefore,this paper proposes a meshless numerical simulation method for reservoirs that considers time-varying characteristics.This method avoids the meshing in traditional numerical simulation methods.From the fluid flow perspective,the reservoir’s computational domain is discretized into a series of connection units.An influence domain with a certain radius centered on the nodes is selected,and one-dimensional connection units are established between the nodes to achieve the characterization of the flow topology structure of the reservoir.In order to reflect the dynamic evolution of the reservoir’s physical properties during the water injection development process,the time-varying characteristics are incorporated into the formula of the seepage characteristic parameters in the meshless calculation.The change relationship of the permeability under different surface fluxes is considered to update the calculated connection conductivity in real time.By combining with the seepage control equation for solution,a time-varying meshless numerical simulation method is formed.The results show that compared with the numerical simulationmethod of the connection elementmethod(CEM)that only considers static parameters,this method has higher simulation accuracy and can better simulate the real migration and distribution of oil and water in the reservoir.Thismethod improves the accuracy of reservoir numerical simulation and the development effect of oilfields,providing a scientific basis for optimizing the water injection strategy,adjusting the production plan,and extending the effective production cycle of the oilfield.
基金supported by the Third Xinjiang Scientific Expedition and Research Program-Investigation and Risk Assessment of Drought and Aeolian Disasters in Tarim River Basin(No.2021xjkk0300)the National Natural Science Foundation of China(No.62466056)the subject of'the technical scheme and application demonstration of sand disaster prevention and control of Xinjiang expressway to engineering practice,Xinjiang Transportation Investment(Group)Co.,Ltd.(No.XJJTZKX-FWCG-202401-0043).
文摘To mitigate the sand burial of highways in sandy regions,a separated subgrade design was widely adopted in the embankments of high-grade highways,but the problem of sand deposition on subgrade slopes and pavements still happens frequently.Based on the theory of wind-sand two-phase flow,this paper constructed a three-dimensional model of the separated subgrade,the wind-sand flow transport law around the subgrade with varying median strip widths and concave depths was simulated by Fluent software.After comparison and analysis of seven subgrade models,the flow field distribution,wind speed horizontal variation,and erosion-deposition characteristics were investigated.The findings are as follows:(1)The width of the median strip in the separated subgrade had significant influences on the wind-sand flow.The smooth passage of wind-sand flow over the road surface was facilitated with the increase of the median strip width.However,sand deposition in the median strip happened.It can lead to secondary sand damage of downwind subgrade and increase the work load of road sand removal for subsequent maintenance.(2)The obstruction to airflow and sand accumulation was aggravated with greater concave depth of the median strip.Therefore,it is advisable to minimize the concave depth of the median strip in case of more sand damage.(3)A median strip width exceeding 12 m(possibly without guardrails)for an integral embankment without enough road land is recommended.Conversely,median strip width of over 40 m for separate subgrade with unrestricted land is suggested.(4)In the case of sand deposition in the existing separated subgrade,the median strip can be filled by sand deposition or other materials,then was covered with gravel to form a flat ground like Gobi smooth surface,which can let the wind-blown sand flow pass through the subgrade section without sand deposition.
基金The SEU Innovation Capability Enhancement Plan for Doctoral Students(No.CXJH_SEU 24115)Marine Economic Development Project of Guangdong Province(No.GDNRC[2022]25).
文摘The study aimed to address the issue of elevated residual stress levels in dissimilar girth welds of cast steel joints.To achieve this,the hybrid welding technology,which yields high welding speeds while simultaneously reducing residual stresses,has been introduced.This study utilizes a numerical simulation method to investigate the temperature and residual stress field in the hybrid welding of G20Mn5 casting-Q355 low-alloy steel welded pipe.A com-parison of the findings of this study with those of other welding processes revealed the technological advantages of hybrid welding.The research outcomes show that due to geometric discontinuities and material differences,the temperature field of the joint exhibits uneven distribution characteristics,and the peak temperatures on the Q355 steel side exceeds those on the G20Mn5 steel side.An evident stress gra-dient is present in the residual stress field of the joint post-welding,with peak stress located at the weld root on the Q355 steel.Compared with arc welding,the hybrid welding leads to decreased residual stresses and deformation,with high stress outside the heat-affected zone diminishing rapidly.Furthermore,it significantly improves the welding efficiency.This study elucidates the distribution and underlying causes of thermal and residual stress fields in dissimilar girth welds.This serves as a foundation for the application of hybrid welding technology in welded cast steel joints.
基金supported by the National Natural Science Foundation of China(Nos.51804024 and 52474342)the Fundamental Research Funds for the Central Universities(No.FRF-IC-20-09)the State Key Laboratory of Advanced Metallurgy of University of Science and Technology Beijing(No.41621002).
文摘To study the combustion behavior of municipal solid waste(MSW)in blast furnace,the combustion process within the raceway was simulated using computational fluid dynamics.Based on the parameters of an actual blast furnace,a three-dimensional model including coal lance,blowpipe,tuyere,and raceway was established.The model was then used to compare the combustion characteristics of pulverized coal and MSW in the raceway and to investigate the effects of blast temperature and particle size on the combustion characteristics of MSW in the raceway.The results showed that MSW combusted more rapidly,achieving a maximum temperature of 3839 K in the raceway,comparing to 2974 K during pulverized coal injection.However,the average temperature during MSW injection was 1790 K,which was 73 K lower than that of pulverized coal injection.The maximum velocity during MSW injection was 120 m/s,lower than 188 m/s obtained during pulverized coal injection.MSW could be completely burned out in the middle of the raceway,while the burnout of pulverized coal at this position was only 50%.The combustion effect of MSW makes no difference when the blast temperature increased from 1273 to 1673 K,due to its excellent combustion characteristic.When the MSW particle size was increased from 0.074 to 2 mm,the burnout was 75%,which was still higher than that of pulverized coal injection with a particle size of 0.074 mm.However,injecting larger-sized fuel might increase the risk of tuyere wear.To ensure stable furnace conditions and great combustion,a blast temperature of 1473 K and a MSW particle size of about 1 mm will be better.
文摘Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 mesh ceramic particles and quartz sand proppant clusters were elaborated using computed tomography(CT)scanning.These models were used to develop a numerical simulation framework based on the lattice Boltzmann method(LBM),enabling the investigation of gas-water flow behavior within proppant-filled fractures under varying driving forces and surface tensions.Simulation results at a closure pressure of 15 MPa have revealed that ceramic particles exhibit a simpler and more porous internal structure than quartz sand of the same size.Under identical flow conditions,ceramic proppants demonstrate higher fluid replacement efficiency.Replacement efficiency increases with higher porosity,greater driving force,and lower surface tension.Furthermore,fluid displacement is strongly influenced by pore geometry:flow is faster in straighter and wider channels,with preferential movement through larger pores forming dominant flow paths.The replacement velocity exhibits a characteristic time evolution,initially rapid,then gradually decreasing,correlating positively with the development of these dominant channels.
基金the financial support from National Key R&D Program of China(Grant number:2024YFC2815100)Natural Science Foundation of China(Grant number:52322110)Beijing Nova Program(Grant number:20230484341).
文摘Deepwater drilling riser is the key equipment connecting the subsea wellhead and floating drilling platform.Due to complex marine environment,vortex-induced vibration(ViV)will be generated on riser,which will induce fatigue failure and even cause unpredictable drilling accidents.Therefore,it is important to study the ViV characteristics of deepwater drilling riser and reveal the main controlling factors for ensuring the safe and efficient operation of deepwater drilling engineering.In this paper,the ViV of deepwater drilling riser is numerically simulated in time domain based on the discrete vortex method(DvM).A hydrodynamic analysis model and governing equation of VIV is proposed with considering the effect of riser motion using DVM and slice method,where the governing equation is solved by Runge-Kutta method.Model validation is performed,which verified the correctness and accuracy of the mechanical model and the solution method.On this basis,the influence of the number of control points,current velocity,riser outer diameter,shear flow and top tension on the ViV characteristics of deepwater drilling risers are discussed in detail.The results show that with the increase of current velocity,the vibration amplitude of deepwater drilling riser decreases obviously,while the vibration frequency increases gradually.However,if the outer diameter of riser increases,the vibration amplitude increases,while the vibration frequency decreases gradually.The top tension also has great influence on the VIV of riser.When the top tension is 1.25 G,the VIV is suppressed to a certain extent.This study has guiding significance for optimal design and engineering control of deepwater drilling riser.
