This paper systematically investigates the numerical simulation model construction and methods for hot dry rock geothermal resource development.It highlights the influence law and characterization differences of multi...This paper systematically investigates the numerical simulation model construction and methods for hot dry rock geothermal resource development.It highlights the influence law and characterization differences of multi-physics field coupling mechanism across various stages of development and utilization.The technical features and applicable scenarios of typical numerical simulation methods,as well as the application potential and advantages of emerging technologies such as intelligent algorithms in numerical simulation for hot dry rock geothermal development,are comprehensively reviewed.In addition,the functional characteristics and engineering application cases of mainstream geothermal numerical simulation software in China and abroad are summarized.On this basis,the core challenges for existing techniques are identified,and future development directions are proposed.At present,numerical simulation for hot dry rock geothermal resource development still faces several challenges,including insufficient accuracy in characterizing complex reservoir structures,incomplete representation of multi-physics field coupling mechanisms,limited cross-scale simulation capability,inadequate adaptability of software to diverse scenarios,and insufficient support from field monitoring and fundamental data.In the future,numerical simulation technologies for hot dry rock geothermal resource development should advance theoretical and technical research in full-chain integrated modeling,refined characterization of multi-physics field coupling,deep integration of intelligent algorithms with numerical simulation,and establishment of an independent and controllable software ecosystem,thereby providing theoretical and technical support for the sustainable and efficient development of hot dry rock geothermal resources in China.展开更多
The uplift resistance of the soil overlying shield tunnels significantly impacts their anti-floating stability.However,research on uplift resistance concerning special-shaped shield tunnels is limited.This study combi...The uplift resistance of the soil overlying shield tunnels significantly impacts their anti-floating stability.However,research on uplift resistance concerning special-shaped shield tunnels is limited.This study combines numerical simulation with machine learning techniques to explore this issue.It presents a summary of special-shaped tunnel geometries and introduces a shape coefficient.Through the finite element software,Plaxis3D,the study simulates six key parameters—shape coefficient,burial depth ratio,tunnel’s longest horizontal length,internal friction angle,cohesion,and soil submerged bulk density—that impact uplift resistance across different conditions.Employing XGBoost and ANN methods,the feature importance of each parameter was analyzed based on the numerical simulation results.The findings demonstrate that a tunnel shape more closely resembling a circle leads to reduced uplift resistance in the overlying soil,whereas other parameters exhibit the contrary effects.Furthermore,the study reveals a diminishing trend in the feature importance of buried depth ratio,internal friction angle,tunnel longest horizontal length,cohesion,soil submerged bulk density,and shape coefficient in influencing uplift resistance.展开更多
Machine learning-assisted methods for rapid and accurate prediction of temperature field,mushy zone,and grain size were proposed for the heating−cooling combined mold(HCCM)horizontal continuous casting of C70250 alloy...Machine learning-assisted methods for rapid and accurate prediction of temperature field,mushy zone,and grain size were proposed for the heating−cooling combined mold(HCCM)horizontal continuous casting of C70250 alloy plates.First,finite element simulations of casting processes were carried out with various parameters to build a dataset.Subsequently,different machine learning algorithms were employed to achieve high precision in predicting temperature fields,mushy zone locations,mushy zone inclination angle,and billet grain size.Finally,the process parameters were quickly optimized using a strategy consisting of random generation,prediction,and screening,allowing the mushy zone to be controlled to the desired target.The optimized parameters are 1234℃for heating mold temperature,47 mm/min for casting speed,and 10 L/min for cooling water flow rate.The optimized mushy zone is located in the middle of the second heat insulation section and has an inclination angle of roughly 7°.展开更多
Energy shortage has become one of themost concerning issues in the world today,and improving energy utilization efficiency is a key area of research for experts and scholars worldwide.Small-diameter heat exchangers of...Energy shortage has become one of themost concerning issues in the world today,and improving energy utilization efficiency is a key area of research for experts and scholars worldwide.Small-diameter heat exchangers offer advantages such as reduced material usage,lower refrigerant charge,and compact structure.However,they also face challenges,including increased refrigerant pressure drop and smaller heat transfer area inside the tubes.This paper combines the advantages and disadvantages of both small and large-diameter tubes and proposes a combined-diameter heat exchanger,consisting of large and small diameters,for use in the indoor units of split-type air conditioners.There are relatively few studies in this area.In this paper,A theoretical and numerical computation method is employed to establish a theoretical-numerical calculation model,and its reliability is verified through experiments.Using this model,the optimal combined diameters and flow path design for a combined-diameter heat exchanger using R32 as the working fluid are derived.The results show that the heat transfer performance of all combined diameter configurations improves by 2.79%to 8.26%compared to the baseline design,with the coefficient of performance(COP)increasing from 4.15 to 4.27~4.5.These designs can save copper material,but at the cost of an increase in pressure drop by 66.86%to 131.84%.The scheme IIIH,using R32,is the optimal combined-diameter and flow path configuration that balances both heat transfer performance and economic cost.展开更多
By simultaneously introducing a finite-difference-based numerical loss term and a clustering-reconstruction mechanism,we propose an enhanced physics-informed neural network named the informed reconstruction-oriented n...By simultaneously introducing a finite-difference-based numerical loss term and a clustering-reconstruction mechanism,we propose an enhanced physics-informed neural network named the informed reconstruction-oriented numerical network(IRON-Net)and subsequently apply it to the Manakov equations-a well-known two-component nonlinear physical model.Numerical experiments are conducted on a dataset containing eight analytical solu-tions with noise.The results indicate that,compared to conventional PINNs and other mainstream algorithms,IRON-Net demonstrates significant advantages in training accuracy,convergence rate,and robustness,achiev-ing a stepwise improvement in the neural network’s ability to enforce physical constraints.Additional ablation experiments further confirm the necessity of the consistency constraint within IRON-Net.This study provides an effective approach for modeling and parameter identification in complex nonlinear optical systems as well as other nonlinear physical scenarios.展开更多
The complex wavefields in ocean-bottom node(OBN)four-component(4C)seismic data,particularly interface waves like Scholte waves,pose significant challenges for processing and inversion.A key issue is how to achieve hig...The complex wavefields in ocean-bottom node(OBN)four-component(4C)seismic data,particularly interface waves like Scholte waves,pose significant challenges for processing and inversion.A key issue is how to achieve high-fidelity and accurate simulation of seismic wave propagation phenomena at the ocean bottom,thereby providing a theoretical foundation for fully leveraging ocean-bottom 4C data in imaging and inversion.To address this,we employ a three-dimensional acoustic-elastic coupling equation(AECE)for high-fidelity numerical simulation of OBN 4C data.This method inherently satisfies the fluid-solid interface conditions without manual enforcement,enabling physically consistent generation of all wave types,including Scholte waves.Numerical examples and field data from the East China Sea demonstrate that the AECE effectively replicates the waveform components and characteristics of field 4C data.This establishes the AECE as a robust foundation for OBN 4C data processing,elastic imaging,and full-waveform inversion.展开更多
The mechanical properties of chondrocytes are closely related to the onset and progression of osteoarthritis;however,current research on the mechanical behavior of chondrocytes undergoing large deformations is insuffi...The mechanical properties of chondrocytes are closely related to the onset and progression of osteoarthritis;however,current research on the mechanical behavior of chondrocytes undergoing large deformations is insufficient.In this work,via micropipette aspiration(MPA)and atomic force microscopy(AFM)experiments on chondrocytes,finite element simulations combined with numerical optimization were conducted to obtain the mechanical parameters of three viscohyperelastic models(neo-Hookean(NH),Mooney-Rivlin(MR),and Arruda-Boyce(AB)).The results showed that for the elastic responses of chondrocytes,all three models can capture the mechanical behaviors of cells with good accuracy for both the MPA and AFM experiments,among which the AB model had the best fit.In terms of the viscoelastic behavior of chondrocytes,the single-term Prony series of the three models can describe the creep response of the MPA experiment well,whereas for the pressure relaxation behavior of the AFM experiment,the fitting degree of the single-term Prony series of the three models was low.However,the prediction ability can be significantly improved by using the two-term Prony series,for both the MPA and the AFM experiments,the AB model still yielded the best prediction of viscoelastic responses.Thus,compared with the NH and MR models,the AB model is more suitable for characterizing the elastic and viscoelastic mechanical responses of chondrocytes undergoing large deformations.This study provides an alternative methodology for investigating the large deformation mechanical properties of chondrocytes,which may help to further study and reveal the mechanotransduction mechanisms of chondrocytes.展开更多
This study presents a numerical investigation of impinging-jet atomization across various Weber numbers(We) under high backpressure conditions.Using the volume-of-fluid method,adaptive mesh refinement,and the iso Adve...This study presents a numerical investigation of impinging-jet atomization across various Weber numbers(We) under high backpressure conditions.Using the volume-of-fluid method,adaptive mesh refinement,and the iso Advector interface reconstruction technique,atomization characteristics are simulated and analyzed for different values of We.The results indicate that the geometry induces turbulent jets,which drive turbulent atomization through the shear interactions at the gas–liquid interface.The key observed phenomena include the interaction of impact waves with liquid sheet perforation and the breakup of web of ligaments,both of which are prominent under high backpressure conditions.A novel method,based on the threshold velocity of spray droplet groups,is employed to quantitatively measure the spreading angle,showing that the angle increases with We in both front and side views.Additionally,the Sauter mean diameter of droplets follows power-law scaling with exponents of-1/3 in the upstream region and-1/2 in the downstream region,while the droplet size distribution conforms to a log-normal profile.This research provides valuable insights into interface evolution and droplet characteristics during impingingjet atomization under high backpressure,offering essential guidance for optimizing industrial atomization processes.展开更多
The numerical manifold method,extensively utilized in numerical computations,faces significant challenges in generating complex manifold elements,particularly for three-dimensional applications.To overcome this challe...The numerical manifold method,extensively utilized in numerical computations,faces significant challenges in generating complex manifold elements,particularly for three-dimensional applications.To overcome this challenge,the meshfree numerical manifold method is developed by integrating the moving least-squares method into the numerical manifold method,effectively bypassing the need for meshing complex geometric objects.However,the implementation of the moving least-squares method introduces computational efficiency issues.To mitigate these,parallel computing methods have been incorporated,resulting in a tenfold increase in the speed of assembling the stiffness matrix with central processing unit parallelism,and a twentyfold increase with graphics processing unit parallelism.The static mechanical system equations for the meshfree numerical manifold method are derived using the Galerkin method.The method’s effectiveness and accuracy are then validated through a series of numerical experiments.The experiments demonstrated that the meshfree numerical manifold method achieves a high precision with minimal nodes and integration points.Additionally,positioning nodes outside the domain significantly improves computational accuracy at the boundaries.展开更多
This study investigated the seismic behavior of columns under different design parameters. A finite element model was established based on a low-cyclic loading experiment of five steel pipe-aeolian sand recycled concr...This study investigated the seismic behavior of columns under different design parameters. A finite element model was established based on a low-cyclic loading experiment of five steel pipe-aeolian sand recycled concrete columns. Their seismic performance was analyzed by adjusting the axial compression ratio, the slenderness ratio, the diameter thickness ratio, and steel pipe strength. The study focused on the hysteretic curve, the skeleton curve, and the displacement ductility coefficient to understand the impact of these parameters. Sensitivity analysis was performed to evaluate the influence of each parameter on peak load and ductility. Test results revealed that stiffness initially increased and then decreased with a displacement rate of aeolian sand. Notably, when the sand replacement rate reached 30%, the stiffness degradation rate was the slowest. The finite element model showed that the P-δ effect of a recycled concrete column varied with the axial compression ratio, and increasing the slenderness ratio decreased ultimate bearing capacity and the displacement ductility coefficient. Conversely, increasing the steel pipe thickness ratio reduced the ultimate bearing capacity and peak displacement but enhanced the displacement ductility coefficient, thereby improving seismic performance. High-strength steel enhanced load capacity but reduced ductility. Sensitivity analysis indicated that the axial compression ratio and the thickness ratio significantly affected peak load, while the thickness ratio had a greater impact on ductility.展开更多
The swimming mechanism of octopuses has been a hot topic in the field of biology in recent years.Existing research primarily focuses on the influence of morphological parameters on octopus swimming ability;the impact ...The swimming mechanism of octopuses has been a hot topic in the field of biology in recent years.Existing research primarily focuses on the influence of morphological parameters on octopus swimming ability;the impact of kinematic parameters remains relatively unexplored.This study focuses on the swimming posture and movement characteristics of octopuses,formulating a synchronized flexible undulation equation for their eight arms.Based on this,computational fluid dynamics combined with dynamic mesh technology is employed to numerically simulate their propulsion mechanism and swimming characteristics.Furthermore,by varying the hover coast time ratio and duty cycle,the influence of different kinematic parameters on propulsion efficiency is explored(close phase after hover coast time ratio is defined as DSC;open phase after hover coast time ratio is defined as DSO,and the duty cycle is defined as DC).The results indicate that the periodic opening and closing motion of the octopus leads to the alternating generation and dissipation of reverse vortices in the wake field,with the jet effect between them being the primary source of propulsion.When DSC=0.30,the forward speed of the octopus after 1 s is only 6%lower than that without hover coast behavior.When DSO=0.30,the peak thrust coefficient of the octopus reaches 0.82,which is 15%higher than that without hover coast behavior,indicating the highest burst acceleration.On the other hand,when DC=0.67,under asymmetric periodic swimming with fast-closing-slow-opening,the average thrust coefficient reaches 0.26,an increase of 23.8%compared to uniform opening-closing periodic swimming,while the average lateral force coefficient decreases by 66.7%.In this scenario,the octopus can achieve better forward propulsion and maintain stable movement.The findings of this study provide valuable insights for research on the intermittent swimming behavior of octopuses and underwater vehicles.展开更多
Collaborative technology for the remote,large-scale deployment of drones using dispersal systems holds significant potential in applications such as post-disaster rescue,which must balance low overload with high thrus...Collaborative technology for the remote,large-scale deployment of drones using dispersal systems holds significant potential in applications such as post-disaster rescue,which must balance low overload with high thrust,in addition to precisely controlling the separation attitude.To address these issues,this paper introduces a multi-gasbag propulsion system with a high aspect ratio that coordinates multiple gasbags to generate sufficient thrust.By adjusting the inlet size of the gasbag,the separation behavior of the release unit can be accurately controlled.A multidimensional two-phase flow model is established,accompanied by both combustion and flow experiments and a double-gasbag propulsion experiment.The results demonstrate that the proposed mathematical model is accurate,effectively captures the pressure fluctuations and spatiotemporal distribution of flow field parameters,and determines the separation attitude of the release unit.For the cases studied in this paper,the pressure at the gasbag inlet(z=650 mm)is the dominant factor during the gasbag propulsion response,causing the release unit to rotate counterclockwise when the gasbag inlet sizes are identical.Increasing the inlet size at z=50 mm compensates for the adverse effects of uneven axial pressure distribution,thereby achieving a neutral separation for the release unit.When the radii r1 and r2 vary between 2 and 12 mm,the angular velocity and attitude angle of the release unit are found to range from-15.50 to 15.20 rad/s and from-0.109 to 0.106 rad,respectively.展开更多
Bubble dynamics near complex boundaries is critical for engineering applications like underwater explosions and cavitation control.This study investigates the collapsing behavior of near-wall bubbles adjacent to three...Bubble dynamics near complex boundaries is critical for engineering applications like underwater explosions and cavitation control.This study investigates the collapsing behavior of near-wall bubbles adjacent to three boundary conditions(planar,elliptical convex,and elliptical concave surfaces)using a compressible multi-component flow model.The finite volume method combined with fifth-order Weighted Essentially Non-Oscillation(WENO)reconstruction and the Harten-Lax-van Leer Contact(HLLC)Riemann solver is employed for spatial discretization,while the third-order Total Variation Diminishing(TVD)Runge-Kutta scheme handles temporal discretization.Results show that elliptical convex and concave surfaces exhibit opposite regulatory effects:the convex surface accelerates bubble collapse,reduces oscillation periods,and increases the water jet pressure peak,whereas the concave surface delays collapse,prolongs periods,and decreases pressure peaks.With increasing stand-off distance ratio,bubble oscillation periods decrease,and minimum equivalent radii also reduce for all boundaries.This work provides insights into complex boundary-induced bubble dynamics,supporting the optimization of cavitation-resistant structures and underwater explosion protection.展开更多
Long-hole raise blasting(LHRB)is a highly efficientexcavation method that is used extensively in underground mining and civil engineering.However,deviations in drilling are usually not accounted for in LHRB,which may ...Long-hole raise blasting(LHRB)is a highly efficientexcavation method that is used extensively in underground mining and civil engineering.However,deviations in drilling are usually not accounted for in LHRB,which may adversely affect the efficiencyand progress of raise excavation.In this paper,the effect of drilling deviation on the optimization of LHRB is investigated.The actual trajectories of the blastholes were measured,and the deviation rates between different diameters were compared.The effects of the drilling deviation on the burn-cut blasting mode(BCBM)and spherical cartridge blasting mode(SCBM)of LHRB were theoretically analyzed.Numerical models with vertical holes that consider the actual hole location at different positions of the raise were subsequently developed to simulate the raise blasting damage.The results indicated that the BCBM relied substantially on the drilling accuracy to provide free surface and compensation space for further blasting,whereas the requirements of drilling deviation for the SCBM were less strict.With increasing hole deviation in the BCBM,the height of the failure area of the raise blasting increased.Optimization designs that combine the BCBM and SCBM were proposed to strike a balance between the efficiencyand reliability of LHRB.A 40 m high slot raise in a large-diameter long-hole(LDL)stope was successfully formed by multimode LHRB.The fieldtest results reveal that the optimization of LHRB is feasible in practical engineering.展开更多
With increase in the number of operations involving relief wells,radial wells,U-shaped wells,and other complex well structures,challenges such as collision prevention,obstacle bypassing,and adjacent-well connectivity ...With increase in the number of operations involving relief wells,radial wells,U-shaped wells,and other complex well structures,challenges such as collision prevention,obstacle bypassing,and adjacent-well connectivity achievement during drilling have become inevitable.These challenges necessitate a technology that can accurately detect adjacent wells in real time during drilling operations.As current borehole acoustic reflection imaging technology heavily relies on cable-based logging,it cannot perform real-time detection of adjacent wells during drilling,thereby limiting the drilling efficiency.This study proposes a new adjacent-well-acoustic-detection-while-drilling method that integrates wireline borehole acoustic reflection imaging with drilling technology,along with an adjacent-well imaging method based on compressed sensing(CS).Together,these methods enable high-resolution,real-time detection of the adjacent target wells during drilling,ensuring safe and efficient underground drilling operations.The finite-difference method was used to simulate three-dimensional numerical models under drilling conditions for two scenarios—with and without target wells adjacent to the drilling well.Experimental validation was conducted in a water tank using an adjacent-well-acoustic-detection-while-drilling tool.The simulated target well was imaged using the CS method,and the imaging results were compared with those obtained from numerical and physical simulations,thereby validating the feasibility of the proposed acoustic detection and imaging methods.The results demonstrate that as the radial distance from the target well increases,the PP echo exhibits delayed arrival times and approaches a plane wave while exhibiting amplitude attenuation.Conversely,a linear increase in the target well diameter advances the PP echo arrival time and enhances its amplitude proportionally.When the target and drilling wells are approximately parallel with a small intersection angle,PP echoes yield better detection results than SS echoes;when the wells are coplanar with a large intersection angle,SS echoes provide better detection results.The receiver element aligned with the target well's azimuth detects all echo modes with the earliest arrival times and highest amplitudes.The adjacent-well imaging method based on CS offers very high spatial resolution,with target wells appearing as local amplitude maxima.This feature enables the precise determination of their azimuth and inclination relative to the drilling wells.The findings offer a solid physical and methodological foundation for real-time detection of adjacent wells during drilling operations and demonstrate enormous theoretical and engineering application potential.展开更多
With the evolution of geophysical surveys from traditional two-dimensional(2 D)to three-dimensional(3 D)models,the resulting large data volumes pose significant challenges to inversion,particularly when resolving larg...With the evolution of geophysical surveys from traditional two-dimensional(2 D)to three-dimensional(3 D)models,the resulting large data volumes pose significant challenges to inversion,particularly when resolving large-scale 3 D structures.A direct solver for solving an ill-conditioned linear system resulting from the finite-difference approximation of a boundary value problem requires more memory and time than iterative solvers.To overcome this limitation,an efficient iterative solver for 3 D finite-difference approach is introduced to calculate the 3 D gravitational potential and the associated gravitational field.Firstly,the boundary value problem associated with 3 D gravitational potential is discretized using central finite-difference technique based on right rectangular prismatic grids.The resulting large unsymmetric sparse systems are then solved using the generalized minimal residual algorithm(GMRES)iterative solver in combination with incomplete LU factorization.Secondly,to obtain high-accuracy partial derivatives of gravitational potential,a high-degree Lagrange interpolation scheme is employed.Finally,three density models are applied to test the accuracy,reliability,and flexibility of our 3 D finite-difference algorithm.All computational results demonstrate that our method provides an accurate approximation of the gravitational field and is applicable to 3 D forward modeling.展开更多
Railways in desert suffer from windblown sand hazards,in which the problems at bridge-road transitions are usually more serious.However,the evolution of wind-sand flow fields and the associated deposition-erosion patt...Railways in desert suffer from windblown sand hazards,in which the problems at bridge-road transitions are usually more serious.However,the evolution of wind-sand flow fields and the associated deposition-erosion patterns around bridge-road transitions under different wind directions remain limited.In this study,numerical simulations are employed to study the influence of wind angle on the flow field structure,wind profiles,horizontal wind speed distribution,and sand deposition-erosion behavior around the bridge-road transition.The results indicate that as the wind angle increases,the leeward vortex recirculation zone expands,the peak of the vertical wind speed profile increases,and the horizontal wind speed decreases.With respect to deposition-erosion patterns,the sand deposition length on the windward sides of both the bridge pier and the bridge-road transition decreases with increasing wind angle.In contrast,the deposition length on the leeward side increases.Specifically,when the wind angle increases from 60°to 90°,the windward deposition length at the bridge-road transition decreases from 1.99H to 1.49H(H,the clear height of bridge),while the leeward deposition expands from 0.70H to 13.28H.In addition,the bridge-road transition section serves as a channel for sand transport,promoting the migration of sand toward the subgrade.A smaller wind angle enhances the lateral guidance effect on the wind-sand flow,resulting in a greater sand deposition length on the leeward side of the subgrade.Therefore,the prevention and control of sand hazards in the bridge-road transition section should not be neglected.It is necessary to regularly remove accumulated sand to prevent its spread from the transition section to the subgrade,which could threaten traffic safety.展开更多
Large-diameter shield tunneling in composite strata presents significant challenges due to the heterogeneous mechanical properties,posing substantial risks to excavation safety.Existing research,however,provides limit...Large-diameter shield tunneling in composite strata presents significant challenges due to the heterogeneous mechanical properties,posing substantial risks to excavation safety.Existing research,however,provides limited insights into how specific layer configurations influence soil strength,failure mechanisms,and the soil arching effect.This study aims to bridge these knowledge gaps by systematically investigating the mechanical behaviors of sand–clay and sand–pebble composites,as well as their excavation-induced responses.Methodologically,we integrate laboratory triaxial tests on reconstituted composite specimens with 3-dimensional(3D)finite element simulations of the tunneling process.Key contributions of this work include:1)the identification of a critical clay layer thickness(20–40 mm)that governs the transition from interfacial slippage to bulging failure in sand–clay composites;2)the quantification of a progressive strength reduction(up to 14.52%)in sand–pebble composites as the sand layer shifts downward;and 3)a novel comparative analysis using a soil stress coefficient method,which reveals a more extensive destruction zone and a more fragile soil arching effect in sand–pebble strata,thereby indicating an elevated risk of collapse.These findings provide crucial insights into the performance of large-diameter shield tunnels in complex geological conditions and offer guidance for safer design and operational control.展开更多
Accurately predicting ski-jump flood discharge atomization is crucial for designing effective disaster-mitigation measures,particularly because low ambient pressure increases the risk of atomized protection in high-al...Accurately predicting ski-jump flood discharge atomization is crucial for designing effective disaster-mitigation measures,particularly because low ambient pressure increases the risk of atomized protection in high-altitude regions.However,owing to the complex effects of low ambient pressure on strongly coupled atomized field sources,it is difficult to fully describe the comprehensive behaviour of such sources theoretically,which limits the further development of random splashing numerical models.In this paper,a refined random splashing numerical model characterized by low ambient pressure is developed based on experimental results and applied to high-altitude earth‒rockfill dam projects.Compared with the reference ambient pressure condition(P0=101.457 kPa),which corresponds to the same flood discharge flow,a decrease in ambient pressure by 0.1 P0leads to a maximum change rate not exceeding 10 m for the characteristic boundary of the 10 mm/h atomized rain intensity line at the QX Hydropower Station.This observation also applies to both the 40 mm/h and 10 mm/h atomized rain intensity lines at the RM Hydropower Station.For the two groups of flip bucket types designed for the RM Hydropower station,the loads associated with atomized protection are predominantly concentrated on the left bank.The maximum height of the 10 mm/h atomized rain intensity line affected by atomized rain ranges from 0.83 to 0.85 times the maximum dam height of 315 m.The distance between the farthest downstream boundary and the Spillway No.3 outlet is between 666.80 and 692.40 metres.Since the flip bucket shape variations only slightly affect the atomization zone extent,further optimization is needed.The study can provide valuable methodological and decision-making support for safeguarding against existing and potential impacts within areas affected by flood discharge atomization from high-altitude hydropower stations.展开更多
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.展开更多
基金Supported by the National Natural Science Foundation of China(52192620)Group Project of Natural Science Innovation Research of China(52421002)Joint Fund Project of the National Natural Science Foundation of China(U24B20325).
摘要This paper systematically investigates the numerical simulation model construction and methods for hot dry rock geothermal resource development.It highlights the influence law and characterization differences of multi-physics field coupling mechanism across various stages of development and utilization.The technical features and applicable scenarios of typical numerical simulation methods,as well as the application potential and advantages of emerging technologies such as intelligent algorithms in numerical simulation for hot dry rock geothermal development,are comprehensively reviewed.In addition,the functional characteristics and engineering application cases of mainstream geothermal numerical simulation software in China and abroad are summarized.On this basis,the core challenges for existing techniques are identified,and future development directions are proposed.At present,numerical simulation for hot dry rock geothermal resource development still faces several challenges,including insufficient accuracy in characterizing complex reservoir structures,incomplete representation of multi-physics field coupling mechanisms,limited cross-scale simulation capability,inadequate adaptability of software to diverse scenarios,and insufficient support from field monitoring and fundamental data.In the future,numerical simulation technologies for hot dry rock geothermal resource development should advance theoretical and technical research in full-chain integrated modeling,refined characterization of multi-physics field coupling,deep integration of intelligent algorithms with numerical simulation,and establishment of an independent and controllable software ecosystem,thereby providing theoretical and technical support for the sustainable and efficient development of hot dry rock geothermal resources in China.
基金Guangzhou Metro Scientific Research Project(No.JT204-100111-23001)Chongqing Municipal Special Project for Technological Innovation and Application Development(No.CSTB2022TIAD-KPX0101)Science and Technology Research and Development Program of China State Railway Group Co.,Ltd.(No.N2023G045)。
摘要The uplift resistance of the soil overlying shield tunnels significantly impacts their anti-floating stability.However,research on uplift resistance concerning special-shaped shield tunnels is limited.This study combines numerical simulation with machine learning techniques to explore this issue.It presents a summary of special-shaped tunnel geometries and introduces a shape coefficient.Through the finite element software,Plaxis3D,the study simulates six key parameters—shape coefficient,burial depth ratio,tunnel’s longest horizontal length,internal friction angle,cohesion,and soil submerged bulk density—that impact uplift resistance across different conditions.Employing XGBoost and ANN methods,the feature importance of each parameter was analyzed based on the numerical simulation results.The findings demonstrate that a tunnel shape more closely resembling a circle leads to reduced uplift resistance in the overlying soil,whereas other parameters exhibit the contrary effects.Furthermore,the study reveals a diminishing trend in the feature importance of buried depth ratio,internal friction angle,tunnel longest horizontal length,cohesion,soil submerged bulk density,and shape coefficient in influencing uplift resistance.
基金financially supported by the National Key Research and Development Program of China (No. 2023YFB3812601)the National Natural Science Foundation of China (No. 51925401)the Young Elite Scientists Sponsorship Program by CAST, China (No. 2022QNRC001)。
摘要Machine learning-assisted methods for rapid and accurate prediction of temperature field,mushy zone,and grain size were proposed for the heating−cooling combined mold(HCCM)horizontal continuous casting of C70250 alloy plates.First,finite element simulations of casting processes were carried out with various parameters to build a dataset.Subsequently,different machine learning algorithms were employed to achieve high precision in predicting temperature fields,mushy zone locations,mushy zone inclination angle,and billet grain size.Finally,the process parameters were quickly optimized using a strategy consisting of random generation,prediction,and screening,allowing the mushy zone to be controlled to the desired target.The optimized parameters are 1234℃for heating mold temperature,47 mm/min for casting speed,and 10 L/min for cooling water flow rate.The optimized mushy zone is located in the middle of the second heat insulation section and has an inclination angle of roughly 7°.
基金supported by Supported by the Scientific Research Foundation for High-Level Talents of Zhoukou Normal University(ZKNUC2024018).
摘要Energy shortage has become one of themost concerning issues in the world today,and improving energy utilization efficiency is a key area of research for experts and scholars worldwide.Small-diameter heat exchangers offer advantages such as reduced material usage,lower refrigerant charge,and compact structure.However,they also face challenges,including increased refrigerant pressure drop and smaller heat transfer area inside the tubes.This paper combines the advantages and disadvantages of both small and large-diameter tubes and proposes a combined-diameter heat exchanger,consisting of large and small diameters,for use in the indoor units of split-type air conditioners.There are relatively few studies in this area.In this paper,A theoretical and numerical computation method is employed to establish a theoretical-numerical calculation model,and its reliability is verified through experiments.Using this model,the optimal combined diameters and flow path design for a combined-diameter heat exchanger using R32 as the working fluid are derived.The results show that the heat transfer performance of all combined diameter configurations improves by 2.79%to 8.26%compared to the baseline design,with the coefficient of performance(COP)increasing from 4.15 to 4.27~4.5.These designs can save copper material,but at the cost of an increase in pressure drop by 66.86%to 131.84%.The scheme IIIH,using R32,is the optimal combined-diameter and flow path configuration that balances both heat transfer performance and economic cost.
基金supported by the Hubei Provincial Natural Science Foundation(Grant No.2023AFB873)the National Natural Science Foun-dation of China(Grant Nos.12505006,11975172,122611-31495,and 12381240286).
摘要By simultaneously introducing a finite-difference-based numerical loss term and a clustering-reconstruction mechanism,we propose an enhanced physics-informed neural network named the informed reconstruction-oriented numerical network(IRON-Net)and subsequently apply it to the Manakov equations-a well-known two-component nonlinear physical model.Numerical experiments are conducted on a dataset containing eight analytical solu-tions with noise.The results indicate that,compared to conventional PINNs and other mainstream algorithms,IRON-Net demonstrates significant advantages in training accuracy,convergence rate,and robustness,achiev-ing a stepwise improvement in the neural network’s ability to enforce physical constraints.Additional ablation experiments further confirm the necessity of the consistency constraint within IRON-Net.This study provides an effective approach for modeling and parameter identification in complex nonlinear optical systems as well as other nonlinear physical scenarios.
基金supported by National Natural Science Foundation of China(42074149)Scientific Research Foundation for High-level Talents of Anhui University of Science and Technology(2025yjrc0001)。
摘要The complex wavefields in ocean-bottom node(OBN)four-component(4C)seismic data,particularly interface waves like Scholte waves,pose significant challenges for processing and inversion.A key issue is how to achieve high-fidelity and accurate simulation of seismic wave propagation phenomena at the ocean bottom,thereby providing a theoretical foundation for fully leveraging ocean-bottom 4C data in imaging and inversion.To address this,we employ a three-dimensional acoustic-elastic coupling equation(AECE)for high-fidelity numerical simulation of OBN 4C data.This method inherently satisfies the fluid-solid interface conditions without manual enforcement,enabling physically consistent generation of all wave types,including Scholte waves.Numerical examples and field data from the East China Sea demonstrate that the AECE effectively replicates the waveform components and characteristics of field 4C data.This establishes the AECE as a robust foundation for OBN 4C data processing,elastic imaging,and full-waveform inversion.
基金supported by the Natural Science Foundation of Shanxi Province(Grant No.202203021221076)Shanxi Huajin Orthopaedic Public Foundation(Grant No.2021066).
摘要The mechanical properties of chondrocytes are closely related to the onset and progression of osteoarthritis;however,current research on the mechanical behavior of chondrocytes undergoing large deformations is insufficient.In this work,via micropipette aspiration(MPA)and atomic force microscopy(AFM)experiments on chondrocytes,finite element simulations combined with numerical optimization were conducted to obtain the mechanical parameters of three viscohyperelastic models(neo-Hookean(NH),Mooney-Rivlin(MR),and Arruda-Boyce(AB)).The results showed that for the elastic responses of chondrocytes,all three models can capture the mechanical behaviors of cells with good accuracy for both the MPA and AFM experiments,among which the AB model had the best fit.In terms of the viscoelastic behavior of chondrocytes,the single-term Prony series of the three models can describe the creep response of the MPA experiment well,whereas for the pressure relaxation behavior of the AFM experiment,the fitting degree of the single-term Prony series of the three models was low.However,the prediction ability can be significantly improved by using the two-term Prony series,for both the MPA and the AFM experiments,the AB model still yielded the best prediction of viscoelastic responses.Thus,compared with the NH and MR models,the AB model is more suitable for characterizing the elastic and viscoelastic mechanical responses of chondrocytes undergoing large deformations.This study provides an alternative methodology for investigating the large deformation mechanical properties of chondrocytes,which may help to further study and reveal the mechanotransduction mechanisms of chondrocytes.
基金partly supported by the National Natural Science Foundation of China(Nos.U23B6009 and 12272050)。
摘要This study presents a numerical investigation of impinging-jet atomization across various Weber numbers(We) under high backpressure conditions.Using the volume-of-fluid method,adaptive mesh refinement,and the iso Advector interface reconstruction technique,atomization characteristics are simulated and analyzed for different values of We.The results indicate that the geometry induces turbulent jets,which drive turbulent atomization through the shear interactions at the gas–liquid interface.The key observed phenomena include the interaction of impact waves with liquid sheet perforation and the breakup of web of ligaments,both of which are prominent under high backpressure conditions.A novel method,based on the threshold velocity of spray droplet groups,is employed to quantitatively measure the spreading angle,showing that the angle increases with We in both front and side views.Additionally,the Sauter mean diameter of droplets follows power-law scaling with exponents of-1/3 in the upstream region and-1/2 in the downstream region,while the droplet size distribution conforms to a log-normal profile.This research provides valuable insights into interface evolution and droplet characteristics during impingingjet atomization under high backpressure,offering essential guidance for optimizing industrial atomization processes.
基金supported by the National Natural Science Foundation of China(Grant Nos.42272338 and 41902275)China Railway Tunnel Group Co.,Ltd.(Grant No.CZ02-08)+4 种基金Sichuan Transportation Science and Technology Program(Grant No.2018-ZL-02)Department of Transportation of Zhejiang Province(Grant No.202213)China Railway First Survey and Design Institute Group Co.,Ltd.(Grant No.2022KY53ZD(CYH)-10)Chongqing Institute of Geology and Mineral Resources(Grant No.TICG-K2024001)Special Project for Performance Incentive and Guidance of Scientific Research Institutions in Chongqing(Grant No.CSTB2023JXJL-YFX0006).
摘要The numerical manifold method,extensively utilized in numerical computations,faces significant challenges in generating complex manifold elements,particularly for three-dimensional applications.To overcome this challenge,the meshfree numerical manifold method is developed by integrating the moving least-squares method into the numerical manifold method,effectively bypassing the need for meshing complex geometric objects.However,the implementation of the moving least-squares method introduces computational efficiency issues.To mitigate these,parallel computing methods have been incorporated,resulting in a tenfold increase in the speed of assembling the stiffness matrix with central processing unit parallelism,and a twentyfold increase with graphics processing unit parallelism.The static mechanical system equations for the meshfree numerical manifold method are derived using the Galerkin method.The method’s effectiveness and accuracy are then validated through a series of numerical experiments.The experiments demonstrated that the meshfree numerical manifold method achieves a high precision with minimal nodes and integration points.Additionally,positioning nodes outside the domain significantly improves computational accuracy at the boundaries.
基金National Natural Science Foundation of China under Grant No. 51868059Inner Mongolia Natural Science Foundation under Grant No. 2021MS05022Basic Scientific Research Expenses Program of Universities directly under Inner Mongolia Autonomous Region under Grant No. JY20220005。
摘要This study investigated the seismic behavior of columns under different design parameters. A finite element model was established based on a low-cyclic loading experiment of five steel pipe-aeolian sand recycled concrete columns. Their seismic performance was analyzed by adjusting the axial compression ratio, the slenderness ratio, the diameter thickness ratio, and steel pipe strength. The study focused on the hysteretic curve, the skeleton curve, and the displacement ductility coefficient to understand the impact of these parameters. Sensitivity analysis was performed to evaluate the influence of each parameter on peak load and ductility. Test results revealed that stiffness initially increased and then decreased with a displacement rate of aeolian sand. Notably, when the sand replacement rate reached 30%, the stiffness degradation rate was the slowest. The finite element model showed that the P-δ effect of a recycled concrete column varied with the axial compression ratio, and increasing the slenderness ratio decreased ultimate bearing capacity and the displacement ductility coefficient. Conversely, increasing the steel pipe thickness ratio reduced the ultimate bearing capacity and peak displacement but enhanced the displacement ductility coefficient, thereby improving seismic performance. High-strength steel enhanced load capacity but reduced ductility. Sensitivity analysis indicated that the axial compression ratio and the thickness ratio significantly affected peak load, while the thickness ratio had a greater impact on ductility.
基金supported by the National Key R&D Program of China(Grant No.2024YFD2400200).
摘要The swimming mechanism of octopuses has been a hot topic in the field of biology in recent years.Existing research primarily focuses on the influence of morphological parameters on octopus swimming ability;the impact of kinematic parameters remains relatively unexplored.This study focuses on the swimming posture and movement characteristics of octopuses,formulating a synchronized flexible undulation equation for their eight arms.Based on this,computational fluid dynamics combined with dynamic mesh technology is employed to numerically simulate their propulsion mechanism and swimming characteristics.Furthermore,by varying the hover coast time ratio and duty cycle,the influence of different kinematic parameters on propulsion efficiency is explored(close phase after hover coast time ratio is defined as DSC;open phase after hover coast time ratio is defined as DSO,and the duty cycle is defined as DC).The results indicate that the periodic opening and closing motion of the octopus leads to the alternating generation and dissipation of reverse vortices in the wake field,with the jet effect between them being the primary source of propulsion.When DSC=0.30,the forward speed of the octopus after 1 s is only 6%lower than that without hover coast behavior.When DSO=0.30,the peak thrust coefficient of the octopus reaches 0.82,which is 15%higher than that without hover coast behavior,indicating the highest burst acceleration.On the other hand,when DC=0.67,under asymmetric periodic swimming with fast-closing-slow-opening,the average thrust coefficient reaches 0.26,an increase of 23.8%compared to uniform opening-closing periodic swimming,while the average lateral force coefficient decreases by 66.7%.In this scenario,the octopus can achieve better forward propulsion and maintain stable movement.The findings of this study provide valuable insights for research on the intermittent swimming behavior of octopuses and underwater vehicles.
基金supported by the National Natural Science Foundation of China(Grant No.52406186)。
摘要Collaborative technology for the remote,large-scale deployment of drones using dispersal systems holds significant potential in applications such as post-disaster rescue,which must balance low overload with high thrust,in addition to precisely controlling the separation attitude.To address these issues,this paper introduces a multi-gasbag propulsion system with a high aspect ratio that coordinates multiple gasbags to generate sufficient thrust.By adjusting the inlet size of the gasbag,the separation behavior of the release unit can be accurately controlled.A multidimensional two-phase flow model is established,accompanied by both combustion and flow experiments and a double-gasbag propulsion experiment.The results demonstrate that the proposed mathematical model is accurate,effectively captures the pressure fluctuations and spatiotemporal distribution of flow field parameters,and determines the separation attitude of the release unit.For the cases studied in this paper,the pressure at the gasbag inlet(z=650 mm)is the dominant factor during the gasbag propulsion response,causing the release unit to rotate counterclockwise when the gasbag inlet sizes are identical.Increasing the inlet size at z=50 mm compensates for the adverse effects of uneven axial pressure distribution,thereby achieving a neutral separation for the release unit.When the radii r1 and r2 vary between 2 and 12 mm,the angular velocity and attitude angle of the release unit are found to range from-15.50 to 15.20 rad/s and from-0.109 to 0.106 rad,respectively.
摘要Bubble dynamics near complex boundaries is critical for engineering applications like underwater explosions and cavitation control.This study investigates the collapsing behavior of near-wall bubbles adjacent to three boundary conditions(planar,elliptical convex,and elliptical concave surfaces)using a compressible multi-component flow model.The finite volume method combined with fifth-order Weighted Essentially Non-Oscillation(WENO)reconstruction and the Harten-Lax-van Leer Contact(HLLC)Riemann solver is employed for spatial discretization,while the third-order Total Variation Diminishing(TVD)Runge-Kutta scheme handles temporal discretization.Results show that elliptical convex and concave surfaces exhibit opposite regulatory effects:the convex surface accelerates bubble collapse,reduces oscillation periods,and increases the water jet pressure peak,whereas the concave surface delays collapse,prolongs periods,and decreases pressure peaks.With increasing stand-off distance ratio,bubble oscillation periods decrease,and minimum equivalent radii also reduce for all boundaries.This work provides insights into complex boundary-induced bubble dynamics,supporting the optimization of cavitation-resistant structures and underwater explosion protection.
基金the National Natural Science Foundation of China(Grant No.52374152),the Guangxi Key R&D Plan(Grant No.2022AB31023)the China Postdoctoral Science Foundation(Grant No.2024M752145).
摘要Long-hole raise blasting(LHRB)is a highly efficientexcavation method that is used extensively in underground mining and civil engineering.However,deviations in drilling are usually not accounted for in LHRB,which may adversely affect the efficiencyand progress of raise excavation.In this paper,the effect of drilling deviation on the optimization of LHRB is investigated.The actual trajectories of the blastholes were measured,and the deviation rates between different diameters were compared.The effects of the drilling deviation on the burn-cut blasting mode(BCBM)and spherical cartridge blasting mode(SCBM)of LHRB were theoretically analyzed.Numerical models with vertical holes that consider the actual hole location at different positions of the raise were subsequently developed to simulate the raise blasting damage.The results indicated that the BCBM relied substantially on the drilling accuracy to provide free surface and compensation space for further blasting,whereas the requirements of drilling deviation for the SCBM were less strict.With increasing hole deviation in the BCBM,the height of the failure area of the raise blasting increased.Optimization designs that combine the BCBM and SCBM were proposed to strike a balance between the efficiencyand reliability of LHRB.A 40 m high slot raise in a large-diameter long-hole(LDL)stope was successfully formed by multimode LHRB.The fieldtest results reveal that the optimization of LHRB is feasible in practical engineering.
基金supported in part by the National Natural Science Foundation of China under Grant Nos.12334019,12274465,12504558 and U25B20244in part by the China Postdoctoral Science Foundation under Grant No.2025M770469+1 种基金in part by the Postdoctoral Fellowship Program of CPSF under Grant No.GZC20251952in part by the Science Foundation of China University of Petroleum,Beijing under Grant No.2462025XKBH014.
摘要With increase in the number of operations involving relief wells,radial wells,U-shaped wells,and other complex well structures,challenges such as collision prevention,obstacle bypassing,and adjacent-well connectivity achievement during drilling have become inevitable.These challenges necessitate a technology that can accurately detect adjacent wells in real time during drilling operations.As current borehole acoustic reflection imaging technology heavily relies on cable-based logging,it cannot perform real-time detection of adjacent wells during drilling,thereby limiting the drilling efficiency.This study proposes a new adjacent-well-acoustic-detection-while-drilling method that integrates wireline borehole acoustic reflection imaging with drilling technology,along with an adjacent-well imaging method based on compressed sensing(CS).Together,these methods enable high-resolution,real-time detection of the adjacent target wells during drilling,ensuring safe and efficient underground drilling operations.The finite-difference method was used to simulate three-dimensional numerical models under drilling conditions for two scenarios—with and without target wells adjacent to the drilling well.Experimental validation was conducted in a water tank using an adjacent-well-acoustic-detection-while-drilling tool.The simulated target well was imaged using the CS method,and the imaging results were compared with those obtained from numerical and physical simulations,thereby validating the feasibility of the proposed acoustic detection and imaging methods.The results demonstrate that as the radial distance from the target well increases,the PP echo exhibits delayed arrival times and approaches a plane wave while exhibiting amplitude attenuation.Conversely,a linear increase in the target well diameter advances the PP echo arrival time and enhances its amplitude proportionally.When the target and drilling wells are approximately parallel with a small intersection angle,PP echoes yield better detection results than SS echoes;when the wells are coplanar with a large intersection angle,SS echoes provide better detection results.The receiver element aligned with the target well's azimuth detects all echo modes with the earliest arrival times and highest amplitudes.The adjacent-well imaging method based on CS offers very high spatial resolution,with target wells appearing as local amplitude maxima.This feature enables the precise determination of their azimuth and inclination relative to the drilling wells.The findings offer a solid physical and methodological foundation for real-time detection of adjacent wells during drilling operations and demonstrate enormous theoretical and engineering application potential.
基金Project(2025ZD1009704)supported by the National Science and Technology Major Project of ChinaProjects(2023JJ30659,2022JJ30706)supported by Hunan Provincial Natural Science Foundation,China。
摘要With the evolution of geophysical surveys from traditional two-dimensional(2 D)to three-dimensional(3 D)models,the resulting large data volumes pose significant challenges to inversion,particularly when resolving large-scale 3 D structures.A direct solver for solving an ill-conditioned linear system resulting from the finite-difference approximation of a boundary value problem requires more memory and time than iterative solvers.To overcome this limitation,an efficient iterative solver for 3 D finite-difference approach is introduced to calculate the 3 D gravitational potential and the associated gravitational field.Firstly,the boundary value problem associated with 3 D gravitational potential is discretized using central finite-difference technique based on right rectangular prismatic grids.The resulting large unsymmetric sparse systems are then solved using the generalized minimal residual algorithm(GMRES)iterative solver in combination with incomplete LU factorization.Secondly,to obtain high-accuracy partial derivatives of gravitational potential,a high-degree Lagrange interpolation scheme is employed.Finally,three density models are applied to test the accuracy,reliability,and flexibility of our 3 D finite-difference algorithm.All computational results demonstrate that our method provides an accurate approximation of the gravitational field and is applicable to 3 D forward modeling.
基金supported by the Major Science and Technology Projects of Inner Mongolia Autonomous Region‘Open Bidding for Selecting the Best Candidates’(2024JBGS0009-01)the National Natural Science Foundation of China(42461011)the Key Research and Development Program of Gansu Province(25YFGA040).
摘要Railways in desert suffer from windblown sand hazards,in which the problems at bridge-road transitions are usually more serious.However,the evolution of wind-sand flow fields and the associated deposition-erosion patterns around bridge-road transitions under different wind directions remain limited.In this study,numerical simulations are employed to study the influence of wind angle on the flow field structure,wind profiles,horizontal wind speed distribution,and sand deposition-erosion behavior around the bridge-road transition.The results indicate that as the wind angle increases,the leeward vortex recirculation zone expands,the peak of the vertical wind speed profile increases,and the horizontal wind speed decreases.With respect to deposition-erosion patterns,the sand deposition length on the windward sides of both the bridge pier and the bridge-road transition decreases with increasing wind angle.In contrast,the deposition length on the leeward side increases.Specifically,when the wind angle increases from 60°to 90°,the windward deposition length at the bridge-road transition decreases from 1.99H to 1.49H(H,the clear height of bridge),while the leeward deposition expands from 0.70H to 13.28H.In addition,the bridge-road transition section serves as a channel for sand transport,promoting the migration of sand toward the subgrade.A smaller wind angle enhances the lateral guidance effect on the wind-sand flow,resulting in a greater sand deposition length on the leeward side of the subgrade.Therefore,the prevention and control of sand hazards in the bridge-road transition section should not be neglected.It is necessary to regularly remove accumulated sand to prevent its spread from the transition section to the subgrade,which could threaten traffic safety.
基金supported by the Shanghai Science and Technology Innovation Program under Grant(No.19DZ1201004).
摘要Large-diameter shield tunneling in composite strata presents significant challenges due to the heterogeneous mechanical properties,posing substantial risks to excavation safety.Existing research,however,provides limited insights into how specific layer configurations influence soil strength,failure mechanisms,and the soil arching effect.This study aims to bridge these knowledge gaps by systematically investigating the mechanical behaviors of sand–clay and sand–pebble composites,as well as their excavation-induced responses.Methodologically,we integrate laboratory triaxial tests on reconstituted composite specimens with 3-dimensional(3D)finite element simulations of the tunneling process.Key contributions of this work include:1)the identification of a critical clay layer thickness(20–40 mm)that governs the transition from interfacial slippage to bulging failure in sand–clay composites;2)the quantification of a progressive strength reduction(up to 14.52%)in sand–pebble composites as the sand layer shifts downward;and 3)a novel comparative analysis using a soil stress coefficient method,which reveals a more extensive destruction zone and a more fragile soil arching effect in sand–pebble strata,thereby indicating an elevated risk of collapse.These findings provide crucial insights into the performance of large-diameter shield tunnels in complex geological conditions and offer guidance for safer design and operational control.
基金funded by the State Key Laboratory of Hydraulic Engineering Intelligent Construction and Operation,Tianjin University(Grants No.HESS-2216)Nanxun scholars program of ZJWEU(Grant No.RC2024011034)+9 种基金Joint Funds of the Zhejiang Provincial Natural Science Foundation of China(Grant No.LZJWY22E090003)Project funded by China Postdoctoral Science Foundation(Grant No.2023M732603)Science and Technology Projects of Xizang Autonomous Region,China(Grant No.XZ202501ZY0109)National Natural Science Foundation of China(Grant No.U1765202)Natural Science Foundation of Tianjin(Grant No.22JCYBJC01180)Major Science and Technology Program of Zhejiang Province(Grant No.2021C03019)Zhejiang Key Laboratory of River-Lake Water Network Health Restoration(Grant No.HHSWKF202501)Central Guidance Funds for Science and Technology Local Development Project(Grant No.2025ZY01091)Scientific Research Fund Key Project of Zhejiang Institute of Hydraulics&Estuary(Zhejiang Institute of Marine Planning&Design)(Grant No.ZIHE25Z002)University-Level Key Course of Zhejiang University of Water Resources and Electric Power(Grant No.ZDKC202319).
摘要Accurately predicting ski-jump flood discharge atomization is crucial for designing effective disaster-mitigation measures,particularly because low ambient pressure increases the risk of atomized protection in high-altitude regions.However,owing to the complex effects of low ambient pressure on strongly coupled atomized field sources,it is difficult to fully describe the comprehensive behaviour of such sources theoretically,which limits the further development of random splashing numerical models.In this paper,a refined random splashing numerical model characterized by low ambient pressure is developed based on experimental results and applied to high-altitude earth‒rockfill dam projects.Compared with the reference ambient pressure condition(P0=101.457 kPa),which corresponds to the same flood discharge flow,a decrease in ambient pressure by 0.1 P0leads to a maximum change rate not exceeding 10 m for the characteristic boundary of the 10 mm/h atomized rain intensity line at the QX Hydropower Station.This observation also applies to both the 40 mm/h and 10 mm/h atomized rain intensity lines at the RM Hydropower Station.For the two groups of flip bucket types designed for the RM Hydropower station,the loads associated with atomized protection are predominantly concentrated on the left bank.The maximum height of the 10 mm/h atomized rain intensity line affected by atomized rain ranges from 0.83 to 0.85 times the maximum dam height of 315 m.The distance between the farthest downstream boundary and the Spillway No.3 outlet is between 666.80 and 692.40 metres.Since the flip bucket shape variations only slightly affect the atomization zone extent,further optimization is needed.The study can provide valuable methodological and decision-making support for safeguarding against existing and potential impacts within areas affected by flood discharge atomization from high-altitude hydropower stations.
基金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.