Drilling and blasting,characterized by their efficiency,ubiquity,and cost-effectiveness,have emerged as predominant techniques in rock excavation;however,they are accompanied by enormous destructive power.Accurately c...Drilling and blasting,characterized by their efficiency,ubiquity,and cost-effectiveness,have emerged as predominant techniques in rock excavation;however,they are accompanied by enormous destructive power.Accurately controlling the blasting energy and achieving the directional fracture of a rock mass have become common problems in the field.A two-dimensional blasting(2D blasting)technique was proposed that utilizes the characteristic that the tensile strength of a rock mass is significantly lower than its compressive strength.After blasting,only a 2D crack surface is generated along the predetermined direction,eliminating the damage to the reserved rock mass caused by conventional blasting.However,the interior of a natural rock mass is a"black box",and the process of crack propagation is difficult to capture,resulting in an unclear 2D blasting mechanism.To this end,a single-hole polymethyl methacrylate(PMMA)test piece was used to conduct a 2D blasting experiment with the help of a high-speed camera to capture the dynamic crack propagation process and the digital image correlation(DIC)method to analyze the evolution law of surface strain on the test piece.On this basis,a three-dimensional(3D)finite element model was established based on the progressive failure theory to simulate the stress,strain,damage,and displacement evolution process of the model under 2D blasting.The simulation results were consistent with the experimental results.The research results reveal the 2D blasting mechanism and provide theoretical support for the application of 2D blasting technology in the field of rock excavation.展开更多
Deeply buried mountain tunnels are often exposed to the risk of rock bursts,which always cause serious damage to the supporting structures and threaten the safety of the engineers.Due to the limited data available,a s...Deeply buried mountain tunnels are often exposed to the risk of rock bursts,which always cause serious damage to the supporting structures and threaten the safety of the engineers.Due to the limited data available,a suitable approach to predict the rockburst tendency at the preliminary stage becomes very important.In this study,an integrated methodology combining 3D initial stress inversion and rockburst tendency prediction was developed and subsequently applied to a case study of the Sangzhuling Tunnel on the Sichuan–Tibet Railway.The numerical modelling involved inverting the initial stress field using a multiple linear regression method.The tunnel excavation was simulated separately by FDM and DEM,based on a stress boundary condition from the inverted stress field.The comparative analysis demonstrates that the rockburst ratio calculated using DEM(76.70%)exhibits a slight increase compared to FDM(75.38%),and the rockburst location is consistent with the actual situation.This suggests that DEM is more suitable for simulating the stress redistribution during excavation in a jointed rock mass.The numerical simulation combined with the deviatoric stress approach effectively predicts rockburst tendency,meeting the engineering requirements.Despite its limitations,numerical simulation remains a reliable method for predicting rock bursts.展开更多
The multi-scale modeling combined with the cohesive zone model(CZM)and the molecular dynamics(MD)method were preformed to simulate the crack propagation in NiTi shape memory alloys(SMAs).The metallographic microscope ...The multi-scale modeling combined with the cohesive zone model(CZM)and the molecular dynamics(MD)method were preformed to simulate the crack propagation in NiTi shape memory alloys(SMAs).The metallographic microscope and image processing technology were employed to achieve a quantitative grain size distribution of NiTi alloys so as to provide experimental data for molecular dynamics modeling at the atomic scale.Considering the size effect of molecular dynamics model on material properties,a reasonable modeling size was provided by taking into account three characteristic dimensions from the perspective of macro,meso,and micro scales according to the Buckinghamπtheorem.Then,the corresponding MD simulation on deformation and fracture behavior was investigated to derive a parameterized traction-separation(T-S)law,and then it was embedded into cohesive elements of finite element software.Thus,the crack propagation behavior in NiTi alloys was reproduced by the finite element method(FEM).The experimental results show that the predicted initiation fracture toughness is in good agreement with experimental data.In addition,it is found that the dynamics initiation fracture toughness increases with decreasing grain size and increasing loading velocity.展开更多
Flash boiling atomization(FBA)is a promising approach for enhancing spray atomization,which can generate a fine and more evenly distributed spray by increasing the fuel injection temperature or reducing the ambient pr...Flash boiling atomization(FBA)is a promising approach for enhancing spray atomization,which can generate a fine and more evenly distributed spray by increasing the fuel injection temperature or reducing the ambient pressure.However,when the outlet speed of the nozzle exceeds 400 m/s,investigating high-speed flash boiling atomization(HFBA)becomes quite challenging.This difficulty arises fromthe involvement ofmany complex physical processes and the requirement for a very fine mesh in numerical simulations.In this study,an HFBA model for gasoline direct injection(GDI)is established.This model incorporates primary and secondary atomization,as well as vaporization and boilingmodels,to describe the development process of the flash boiling spray.Compared to lowspeed FBA,these physical processes significantly impact HFBA.In this model,the Eulerian description is utilized for modeling the gas,and the Lagrangian description is applied to model the droplets,which effectively captures the movement of the droplets and avoids excessive mesh in the Eulerian coordinates.Under various conditions,numerical solutions of the Sauter mean diameter(SMD)for GDI show good agreement with experimental data,validating the proposed model’s performance.Simulations based on this HFBA model investigate the influences of fuel injection temperature and ambient pressure on the atomization process.Numerical analyses of the velocity field,temperature field,vapor mass fraction distribution,particle size distribution,and spray penetration length under different superheat degrees reveal that high injection temperature or low ambient pressure significantly affects the formation of small and dispersed droplet distribution.This effect is conducive to the refinement of spray particles and enhances atomization.展开更多
Small-scale measurements of the radon exhalation rate using the flow-through and closed-loop methods were conducted on the surface of a uranium tailing pond to better understand the differences between the two methods...Small-scale measurements of the radon exhalation rate using the flow-through and closed-loop methods were conducted on the surface of a uranium tailing pond to better understand the differences between the two methods.An abnormal radon exhalation behavior was observed,leading to computational fluid dynamics(CFD)-based simulations in which dynamic radon migration in a porous medium and accumulation chamber was considered.Based on the in-situ experimental and numerical simulation results,variations in the radon exhalation rate subject to permeability,flow rate,and insertion depth were quantified and analyzed.The in-situ radon exhalation rates measured using the flow-through method were higher than those measured using the closed-loop method,which could be explained by the negative pressure difference between the inside and outside of the chamber during the measurements.The consistency of the variations in the radon exhalation rate between the experiments and simulations suggests the reliability of CFD-based techniques in obtaining the dynamic evolution of transient radon exhalation rates for diffusion and convection at the porous medium-air interface.The synergistic effects of the three factors(insertion depth,flow rate,and permeability)on the negative pressure difference and measured exhalation rate were quantified,and multivariate regression models were established,with positive correlations in most cases;the exhalation rate decreased with increasing insertion depth at a permeability of 1×10^(−11) m^(2).CFD-based simulations can provide theoretical guidance for improving the flow-through method and thus achieve accurate measurements.展开更多
Accurate measurements of the radon exhalation rate help identify and evaluate radon risk regions in the environment.Among these measurement methods,the closed-loop method is frequently used.However,traditional experim...Accurate measurements of the radon exhalation rate help identify and evaluate radon risk regions in the environment.Among these measurement methods,the closed-loop method is frequently used.However,traditional experiments are insufficient or cannot analyze the radon migration and exhalation patterns at the gas–solid interface in the accumulation chamber.The CFD-based technique was applied to predict the radon concentration distribution in a limited space,allowing radon accumulation and exhalation inside the chamber intuitively and visually.In this study,three radon exhalation rates were defined,and two structural ventilation tubes were designed for the chamber.The consistency of the simulated results with the variation in the radon exhalation rate in a previous experiment or analytical solution was verified.The effects of the vent tube structure and flow rate on the radon uniformity in the chamber;permeability,insertion depth,and flow rate on the radon exhalation rate and the effective diffusion coefficient on back-diffusion were investigated.Based on the results,increasing the inser-tion depth from 1 to 5 cm decreased the effective decay constant by 19.55%,whereas the curve-fitted radon exhalation rate decreased(lower than the initial value)as the deviation from the initial value increased by approximately 7%.Increasing the effective diffusion coefficient from 2.77×10^(-7) to 7.77×10^(-6) m^(2) s^(-1) made the deviation expand from 2.14 to 15.96%.The conclusion is that an increased insertion depth helps reduce leakage in the chamber,subject to notable back-diffusion,and that the closed-loop method is reasonably used for porous media with a low effective diffusion coefficient in view of the back-diffusion effect.The CFD-based simulation is expected to provide guidance for the optimization of the radon exhalation rate measurement method and,thus,the accurate measurement of the radon exhalation rate.展开更多
As an important lightning protection device in substations,lightning rods are susceptible to vibration and potential structural damage under wind loads.In order to understand their vibration mechanism,it is necessary ...As an important lightning protection device in substations,lightning rods are susceptible to vibration and potential structural damage under wind loads.In order to understand their vibration mechanism,it is necessary to conduct flow analysis.In this study,numerical simulations of the flow field around a 330 kV cylindrical lightning rod with different diameters were performed using the SST k-ωmodel.The flow patterns in different segments of the lightning rod at the same reference wind speed(wind speed at a height of 10 m)and the flow patterns in the same segment at different reference wind speeds were investigated.The variations of lift coefficient,drag coefficient,and vorticity distribution were obtained.The results showed that vortex shedding phenomena occurred in all segments of the lightning rod,and the strength of vortex shedding increased with decreasing diameter.The vorticity magnitude and the root mean square magnitudes of the lift coefficient and drag coefficient also increased accordingly.The time history curves of the lift coefficient and drag coefficient on the surface of the lightning rod exhibited sinusoidal patterns with a single dominant frequency.For the same segment,as the wind speed increased in a certain range,the root mean square values of the lift coefficient and drag coefficient decreased,while their dominant frequencies increased.Moreover,there was a proportional relationship between the dominant frequencies of the lift coefficient and drag coefficient.The findings of this study can provide valuable insights for the refined design of lightning rods with similar structures.展开更多
It is one concern of the researchers how magnesium(Mg)alloys solidify under different conditions and how their microstructure evolves during solidification,and what are the relationship between the macroscopic propert...It is one concern of the researchers how magnesium(Mg)alloys solidify under different conditions and how their microstructure evolves during solidification,and what are the relationship between the macroscopic properties and various microstructures.Such issues are difficult to be revealed through experiments only,especially for the newly developed Mg alloys,for which there is a lack of more systematic and mature system.However,multi-scale modeling and simulation can promote and deepen our understanding of the microstructure and its deformation mechanism.In this paper,we review and summarize the recent research progress of numerical simulation of Mg alloys in forming and microstructure,namely casting,extrusion,rolling,and welding,using crystal plasticity finite element(CPFEM)and molecular dynamics(DM)methods.Besides,the methods and innovations of modeling are also summarized.Lastly,the paper discusses the development prospects and challenges of the numerical simulation in the field of Mg alloys.展开更多
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.展开更多
This study investigates the flexural performance of ultra-high performance concrete(UHPC)in reinforced concrete T-beams,focusing on the effects of interfacial treatments.Three concrete T-beam specimens were fabricated...This study investigates the flexural performance of ultra-high performance concrete(UHPC)in reinforced concrete T-beams,focusing on the effects of interfacial treatments.Three concrete T-beam specimens were fabricated and tested:a control beam(RC-T),a UHPC-reinforced beam with a chiseled interface(UN-C-50F),and a UHPC-reinforced beam featuring both a chiseled interface and anchored steel rebars(UN-CS-50F).The test results indicated that both chiseling and the incorporation of anchored rebars effectively created a synergistic combination between the concrete T-beam and the UHPC reinforcement layer,with the UN-CS-50F exhibiting the highest flexural resistance.The cracking load and ultimate load of UN-CS-50F were 221.5%and 40.8%,respectively,higher than those of the RC-T.Finite element(FE)models were developed to provide further insights into the behavior of the UHPCreinforced T-beams,showing a maximumdeviation of just 8%when validated against experimental data.A parametric analysis varied the height,thickness,andmaterial strength of the UHPC reinforcement layer based on the validated FE model,revealing that increasing the UHPC layer thickness from 30 to 50 mm improved the ultimate resistance by 20%while reducing the UHPC reinforcement height from 440 to 300 mm led to a 10%decrease in bending resistance.The interfacial anchoring rebars significantly reduced crack propagation and enhanced stress redistribution,highlighting the importance of strengthening interfacial bonds and optimizing geometric parameters ofUHPCfor improved T-beam performance.These findings offer valuable insights for the design and retrofitting of UHPC-reinforced bridge girders.展开更多
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.展开更多
The detonation of fuel-rich explosives yields combustible products that persistently burn upon mixing with ambient oxygen,releasing additional energy through a phenomenon known as the afterburning effect.This process ...The detonation of fuel-rich explosives yields combustible products that persistently burn upon mixing with ambient oxygen,releasing additional energy through a phenomenon known as the afterburning effect.This process greatly influences the evolution of confined blast loading and the subsequent structural response,which is crucial in confined blast scenarios.Given the complex nature of the reaction process,accurate analysis of the afterburning effect remains challenging.Previous studies have either overlooked the mechanisms of detonation product combustion or failed to provide experimental validation.This study introduces a three-dimensional model to effectively characterize the combustion of detonation products.The model integrates chemical reaction source terms into the governing equations to consider the combustion processes.Numerical simulations and experimental tests were conducted to analyze the combustion and energy release from the detonation products of fuel-rich explosives in confined spaces.Approximately 50%of the energy was released during the combustion of detonation products in a confined TNT explosion.Although the combustion of these products was much slower than the detonation process,it aligned with the dynamic response of the structure,which enhanced the explosive yield.Excluding afterburning from the analysis reduced the center-point deformation of the structure by 30%.Following the inclusion of afterburning,the simulated quasistatic pressure increased by approximately 45%.Subsequent comparisons highlighted the merits of the proposed approach over conventional methods.This approach eliminates the reliance on empirical parameters,such as the amount and rate of energy release during afterburning,thereby laying the foundation for understanding load evolution in more complex environments,such as ships,buildings,and underground tunnels.展开更多
This paper investigates the impact of the model top and damping layer on the numerical simulation of tropical cyclones(TCs)and reveals the significant role of stratospheric gravity waves(SGWs).TCs can generate SGWs,wh...This paper investigates the impact of the model top and damping layer on the numerical simulation of tropical cyclones(TCs)and reveals the significant role of stratospheric gravity waves(SGWs).TCs can generate SGWs,which propagate upward and outward into the stratosphere.These SGWs can reach the damping layer,which is a consequence of the numerical scheme employed,where they can affect the tangential circulation through the dragging and forcing processes.In models with a higher top boundary,this tangential circulation develops far from the TC and has minimal direct impact on TC intensity.By comparison,in models with a lower top(e.g.,20 km),the damping layer is located just above the top of the TC.The SGW dragging in the damping layer and the consequent tangential force can thus induce ascent outside the eyewall,promote latent heat release,tilt the eyewall,and enlarge the inner-core radius.This process will reduce inner-core vorticity advection within the boundary layer,and eventually inhibits the intensification of the TC.This suggests that when the thickness of the damping layer is 5 km,the TC numerical model top height should be at least higher than 20 km to generate more accurate simulations.展开更多
A three-dimensional coupled sea ice-ice shelf-ocean numerical model is developed for the Prydz Bay,Antarctica,using the Regional Ocean Modeling System with a grid resolution of approximately 2 km.The influence of the ...A three-dimensional coupled sea ice-ice shelf-ocean numerical model is developed for the Prydz Bay,Antarctica,using the Regional Ocean Modeling System with a grid resolution of approximately 2 km.The influence of the grounding giant iceberg D15 on the distribution of sea ice and polynyas in the Prydz Bay is analyzed through two numerical experiments.Iceberg D15,grounded off the western edge of the West Ice Shelf(WIS),obstructs the southwestward transport of sea ice along the east coast of Prydz Bay,causing sea ice to accumulate to the east of the iceberg and form multi-year fast ice.Grounding of Iceberg D15 also decreases sea ice coverage off its south edge and creates ice-free openings in spring near Davis Station and Zhongshan Station,facilitating the accessibility of vessels to the research stations.These simulated sea ice patterns closely match current satellite observations.When Iceberg D15 is removed,the previously blocked sea ice north of the iceberg,which moved westward,shifts southwesterly along the coastline,leading to a reduction in sea ice thickness during winter and spring,as well as lower sea ice concentrations in spring across large areas north and west of the iceberg.In contrast,the sea ice thickness increases considerably southwest of the WIS,extending to the front of the Amery Ice Shelf during seasons covered by sea ice.The increase in sea ice concentration can also extend to as far as 75°E in spring.Without Iceberg D15,which previously contributed to the ice barrier of Barrier Polynya(BP),the shape of BP changes,the area of BP and Davis Polynya(DP)decreases,and the polynya off the northwest edge of the WIS near 83°E expands.These polynya patterns are much similar to the satellite remote sensing observations before Iceberg D15 was grounded.From April to October,the total area of BP and DP decreases by 2.83×10^(4)km^(2)(60%)and 2.20×10^(3)km^(2)(20%),respectively,while the total sea ice production decreases by 4.11×10^(10)m^(3)(66%)and 1.52×10^(10)m^(3)(52%)compared to the experiment with iceberg.These results indicate the substantial effects of grounding giant icebergs on the spatio-temporal distribution of sea ice,the area of polynyas,and sea ice production.High-resolution Antarctic coastal numerical models,typically with grid scales of kilometers,are sufficient to represent large icebergs,and adding the grounding giant icebergs is crucial for producing realistic simulations of sea ice and polynyas.展开更多
It has been experimentally observed that,in the perforation of metal plates by a flat-nosed projectile,there exists a plateau phenomenon where the ballistic limit increases slightly with increasing plate thickness,whi...It has been experimentally observed that,in the perforation of metal plates by a flat-nosed projectile,there exists a plateau phenomenon where the ballistic limit increases slightly with increasing plate thickness,which is related to a change in the mode of failure.No theoretical model has so far explained this phenomenon satisfactorily.This paper presents a combined numerical and theoretical study on the perforation of 2024-T351 aluminum plates struck by flat-nosed projectiles.First,numerical simulations are performed to investigate the failure mechanisms/deformation modes of the aluminum plates.Then,a theoretical model is proposed based on the numerical results and the experimental observations within a unified framework.The model takes into account the main energy absorbing mechanisms and the corresponding energies absorbed are determined analytically.In particular,a dimensionless equation is suggested to describe the relationship between global deformations and impact velocity.It transpires that the model predictions are in good agreement with the test data and the numerical results for the perforation of 2024-T351 aluminum plates struck by rigid flat-nosed projectiles in terms of residual velocity,ballistic limit,relationship between global deformations and impact velocity,and transition of failure modes.It also transpires that the present model can predict the“plateau”phenomenon,which shows a slight increase in ballistic limit as plate thickness increases.Furthermore,the energy absorption mechanisms are discussed on the basis of the theoretical analysis.展开更多
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.展开更多
Plastometric experiments,supplemented with numerical simulations using the finite element method(FEM),can be advantageously used to characterize the deformation behavior of metallic materials.The accuracy of such simu...Plastometric experiments,supplemented with numerical simulations using the finite element method(FEM),can be advantageously used to characterize the deformation behavior of metallic materials.The accuracy of such simulations predicting deformation behaviors of materials is,however,primarily affected by the applied rheology law.The presented study focuses on the characterization of the deformation behavior of AISI 1045 type carbon steel,widely used e.g.,in automotive and power engineering,under extreme conditions(i.e.,high temperatures,strain rates).The study consists of two main parts:experimentally analyzing the flow stress development of the steel under different thermomechanical conditions via uniaxial hot compression tests and establishing the rheology law via numerical simulations implementing the experimentally acquired flow stress curves.The numerical simulations then not only serve to establish the rheology law but also to verify the reliability of the selected experimental process.The results of the numerical simulations showed that the established rheology law characterizes the behavior of the investigated steel with sufficient accuracy also at high temperatures and/or strain rates,and can,therefore,be used for practical purposes.Last but not least,supplementary microstructure analyses performed for the samples subjected to the highest deformation temperature provided a deeper insight into the effects of the applied(extreme)thermomechanical conditions on the behavior of the investigated steel.展开更多
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%.展开更多
With the increasing development of deepburied engineering projects,rockburst disasters have become a frequent concern.Studies have indicated that tunnel diameter is a critical factor influencing the occurrence of rock...With the increasing development of deepburied engineering projects,rockburst disasters have become a frequent concern.Studies have indicated that tunnel diameter is a critical factor influencing the occurrence of rockbursts.To investigate the influence of tunnel diameter on the deformation and failure characteristics of surrounding rock,large-sized rocklike gypsum specimens were tested using a selfdeveloped true triaxial rockburst loading system containing circular tunnels with three different diameters(D=0.07 m,0.11 m,and 0.15 m).Acoustic emission monitoring,together with a miniature intelligent camera,was employed to analyze the entire process,focusing on macroscopic failure patterns,fragment characteristics,and underlying failure mechanisms.In addition,theoretical analyses were carried out and combined with numerical simulations to investigate the differences in energy evolution associated with rockburst physical models.The results indicate that:(1)The rockburst process with different tunnel diameters consistently evolved through three distinct stages—initial particle ejection,crack propagation accompanied by flake spalling,and,finally,fragment ejection leading to the formation of a‘V'-shaped notch.(2)Increasing tunnel diameter reduces rockburst failure load while increasing surrounding rock damage extent,total mass and average size of ejected fragments.Additionally,shear failure proportion decreases with tensile failure becoming increasingly dominant.(3)Larger tunnel diameters reduce the attenuation rate of elastic strain energy,thereby expanding the zone of elastic strain energy accumulation and disturbance and creating conditions for larger volume rockburst.(4)Larger tunnel diameters result in a smaller principal stress ratio at equivalent distances in the surrounding rock,indicating a higher likelihood of tensile failure.(5)Numerical analyses further reveal that larger tunnel diameters reduce the maximum elastic strain energy density around the tunnel,lowering the energy released per unit volume of rockburst fragments and their ejection velocities.However,both the total failure volume and overall energy release from rockburst increase.Model experiments with different tunnel diameters are of great significance for optimizing engineering design and parameter selection,as well as guiding tunnel construction under complex geological conditions.展开更多
The pressure wave generated by an urban-rail vehicle when passing through a tunnel affects the comfort of passengersand may even cause damage to the train and related tunnel structures.Therefore,controlling the trains...The pressure wave generated by an urban-rail vehicle when passing through a tunnel affects the comfort of passengersand may even cause damage to the train and related tunnel structures.Therefore,controlling the trainspeed in the tunnel is extremely important.In this study,this problem is investigated numerically in the frameworkof the standard k-εtwo-equation turbulence model.In particular,an eight-car urban rail train passingthrough a tunnel at different speeds(140,160,180 and 200 km/h)is considered.The results show that the maximumaerodynamic drag of the head and tail cars is most affected by the running speed.The pressure at selectedmeasuring points on the windward side of the head car is very high,and the negative pressure at the side windowof the driver’s cab of the tail car is also very large.From the head car to the tail car,the pressure at the same heightgradually decreases.The positive pressure peak at the head car and the negative pressure peak at the tail car aregreatly affected by the speed.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.52404155 and 52304111)State Key Laboratory for Geomechanics and Deep Underground Engineering,China University of Mining&Technology,Beijing(Grant No.XD2024006).
文摘Drilling and blasting,characterized by their efficiency,ubiquity,and cost-effectiveness,have emerged as predominant techniques in rock excavation;however,they are accompanied by enormous destructive power.Accurately controlling the blasting energy and achieving the directional fracture of a rock mass have become common problems in the field.A two-dimensional blasting(2D blasting)technique was proposed that utilizes the characteristic that the tensile strength of a rock mass is significantly lower than its compressive strength.After blasting,only a 2D crack surface is generated along the predetermined direction,eliminating the damage to the reserved rock mass caused by conventional blasting.However,the interior of a natural rock mass is a"black box",and the process of crack propagation is difficult to capture,resulting in an unclear 2D blasting mechanism.To this end,a single-hole polymethyl methacrylate(PMMA)test piece was used to conduct a 2D blasting experiment with the help of a high-speed camera to capture the dynamic crack propagation process and the digital image correlation(DIC)method to analyze the evolution law of surface strain on the test piece.On this basis,a three-dimensional(3D)finite element model was established based on the progressive failure theory to simulate the stress,strain,damage,and displacement evolution process of the model under 2D blasting.The simulation results were consistent with the experimental results.The research results reveal the 2D blasting mechanism and provide theoretical support for the application of 2D blasting technology in the field of rock excavation.
基金financially supported by the State Key Laboratory of Geohazard Prevention and Geoenvironment Protection(Chengdu University of Technology)(Grant No.SKLGP2020Z007)。
文摘Deeply buried mountain tunnels are often exposed to the risk of rock bursts,which always cause serious damage to the supporting structures and threaten the safety of the engineers.Due to the limited data available,a suitable approach to predict the rockburst tendency at the preliminary stage becomes very important.In this study,an integrated methodology combining 3D initial stress inversion and rockburst tendency prediction was developed and subsequently applied to a case study of the Sangzhuling Tunnel on the Sichuan–Tibet Railway.The numerical modelling involved inverting the initial stress field using a multiple linear regression method.The tunnel excavation was simulated separately by FDM and DEM,based on a stress boundary condition from the inverted stress field.The comparative analysis demonstrates that the rockburst ratio calculated using DEM(76.70%)exhibits a slight increase compared to FDM(75.38%),and the rockburst location is consistent with the actual situation.This suggests that DEM is more suitable for simulating the stress redistribution during excavation in a jointed rock mass.The numerical simulation combined with the deviatoric stress approach effectively predicts rockburst tendency,meeting the engineering requirements.Despite its limitations,numerical simulation remains a reliable method for predicting rock bursts.
基金Funded by the National Natural Science Foundation of China Academy of Engineering Physics and Jointly Setup"NSAF"Joint Fund(No.U1430119)。
文摘The multi-scale modeling combined with the cohesive zone model(CZM)and the molecular dynamics(MD)method were preformed to simulate the crack propagation in NiTi shape memory alloys(SMAs).The metallographic microscope and image processing technology were employed to achieve a quantitative grain size distribution of NiTi alloys so as to provide experimental data for molecular dynamics modeling at the atomic scale.Considering the size effect of molecular dynamics model on material properties,a reasonable modeling size was provided by taking into account three characteristic dimensions from the perspective of macro,meso,and micro scales according to the Buckinghamπtheorem.Then,the corresponding MD simulation on deformation and fracture behavior was investigated to derive a parameterized traction-separation(T-S)law,and then it was embedded into cohesive elements of finite element software.Thus,the crack propagation behavior in NiTi alloys was reproduced by the finite element method(FEM).The experimental results show that the predicted initiation fracture toughness is in good agreement with experimental data.In addition,it is found that the dynamics initiation fracture toughness increases with decreasing grain size and increasing loading velocity.
基金supported by the National Natural Science Foundation of China(Project Nos.12272270,11972261).
文摘Flash boiling atomization(FBA)is a promising approach for enhancing spray atomization,which can generate a fine and more evenly distributed spray by increasing the fuel injection temperature or reducing the ambient pressure.However,when the outlet speed of the nozzle exceeds 400 m/s,investigating high-speed flash boiling atomization(HFBA)becomes quite challenging.This difficulty arises fromthe involvement ofmany complex physical processes and the requirement for a very fine mesh in numerical simulations.In this study,an HFBA model for gasoline direct injection(GDI)is established.This model incorporates primary and secondary atomization,as well as vaporization and boilingmodels,to describe the development process of the flash boiling spray.Compared to lowspeed FBA,these physical processes significantly impact HFBA.In this model,the Eulerian description is utilized for modeling the gas,and the Lagrangian description is applied to model the droplets,which effectively captures the movement of the droplets and avoids excessive mesh in the Eulerian coordinates.Under various conditions,numerical solutions of the Sauter mean diameter(SMD)for GDI show good agreement with experimental data,validating the proposed model’s performance.Simulations based on this HFBA model investigate the influences of fuel injection temperature and ambient pressure on the atomization process.Numerical analyses of the velocity field,temperature field,vapor mass fraction distribution,particle size distribution,and spray penetration length under different superheat degrees reveal that high injection temperature or low ambient pressure significantly affects the formation of small and dispersed droplet distribution.This effect is conducive to the refinement of spray particles and enhances atomization.
基金National Natural Science Foundation of China(No.11575080)Hunan Provincial Natural Science Foundation of China(No.2022JJ30482)Hunan Provincial Innovation Foundation for Postgraduate(No.QL20220206).
文摘Small-scale measurements of the radon exhalation rate using the flow-through and closed-loop methods were conducted on the surface of a uranium tailing pond to better understand the differences between the two methods.An abnormal radon exhalation behavior was observed,leading to computational fluid dynamics(CFD)-based simulations in which dynamic radon migration in a porous medium and accumulation chamber was considered.Based on the in-situ experimental and numerical simulation results,variations in the radon exhalation rate subject to permeability,flow rate,and insertion depth were quantified and analyzed.The in-situ radon exhalation rates measured using the flow-through method were higher than those measured using the closed-loop method,which could be explained by the negative pressure difference between the inside and outside of the chamber during the measurements.The consistency of the variations in the radon exhalation rate between the experiments and simulations suggests the reliability of CFD-based techniques in obtaining the dynamic evolution of transient radon exhalation rates for diffusion and convection at the porous medium-air interface.The synergistic effects of the three factors(insertion depth,flow rate,and permeability)on the negative pressure difference and measured exhalation rate were quantified,and multivariate regression models were established,with positive correlations in most cases;the exhalation rate decreased with increasing insertion depth at a permeability of 1×10^(−11) m^(2).CFD-based simulations can provide theoretical guidance for improving the flow-through method and thus achieve accurate measurements.
基金This work was supported by the National Natural Science Foundation of China(No.11575080)the National Natural Science Foundation of Hunan Province,China(No.2022JJ30482)the Hunan Provincial Innovation Foundation for Postgraduates(No.QL20220206).
文摘Accurate measurements of the radon exhalation rate help identify and evaluate radon risk regions in the environment.Among these measurement methods,the closed-loop method is frequently used.However,traditional experiments are insufficient or cannot analyze the radon migration and exhalation patterns at the gas–solid interface in the accumulation chamber.The CFD-based technique was applied to predict the radon concentration distribution in a limited space,allowing radon accumulation and exhalation inside the chamber intuitively and visually.In this study,three radon exhalation rates were defined,and two structural ventilation tubes were designed for the chamber.The consistency of the simulated results with the variation in the radon exhalation rate in a previous experiment or analytical solution was verified.The effects of the vent tube structure and flow rate on the radon uniformity in the chamber;permeability,insertion depth,and flow rate on the radon exhalation rate and the effective diffusion coefficient on back-diffusion were investigated.Based on the results,increasing the inser-tion depth from 1 to 5 cm decreased the effective decay constant by 19.55%,whereas the curve-fitted radon exhalation rate decreased(lower than the initial value)as the deviation from the initial value increased by approximately 7%.Increasing the effective diffusion coefficient from 2.77×10^(-7) to 7.77×10^(-6) m^(2) s^(-1) made the deviation expand from 2.14 to 15.96%.The conclusion is that an increased insertion depth helps reduce leakage in the chamber,subject to notable back-diffusion,and that the closed-loop method is reasonably used for porous media with a low effective diffusion coefficient in view of the back-diffusion effect.The CFD-based simulation is expected to provide guidance for the optimization of the radon exhalation rate measurement method and,thus,the accurate measurement of the radon exhalation rate.
基金supported by State Grid Ningxia Electric Power Co.,Ltd.under Grant 5229CG220006Natural Science Foundation of Ningxia Province under Grant 2022AAC03629.
文摘As an important lightning protection device in substations,lightning rods are susceptible to vibration and potential structural damage under wind loads.In order to understand their vibration mechanism,it is necessary to conduct flow analysis.In this study,numerical simulations of the flow field around a 330 kV cylindrical lightning rod with different diameters were performed using the SST k-ωmodel.The flow patterns in different segments of the lightning rod at the same reference wind speed(wind speed at a height of 10 m)and the flow patterns in the same segment at different reference wind speeds were investigated.The variations of lift coefficient,drag coefficient,and vorticity distribution were obtained.The results showed that vortex shedding phenomena occurred in all segments of the lightning rod,and the strength of vortex shedding increased with decreasing diameter.The vorticity magnitude and the root mean square magnitudes of the lift coefficient and drag coefficient also increased accordingly.The time history curves of the lift coefficient and drag coefficient on the surface of the lightning rod exhibited sinusoidal patterns with a single dominant frequency.For the same segment,as the wind speed increased in a certain range,the root mean square values of the lift coefficient and drag coefficient decreased,while their dominant frequencies increased.Moreover,there was a proportional relationship between the dominant frequencies of the lift coefficient and drag coefficient.The findings of this study can provide valuable insights for the refined design of lightning rods with similar structures.
基金supported by the National Natural Science Foundation of China(No.52271091)Natural Science Foundation Project of Ningxia Province(No.2023AAC03324)the National Key Research and Development Program of China(No.2021YFB3701100).
文摘It is one concern of the researchers how magnesium(Mg)alloys solidify under different conditions and how their microstructure evolves during solidification,and what are the relationship between the macroscopic properties and various microstructures.Such issues are difficult to be revealed through experiments only,especially for the newly developed Mg alloys,for which there is a lack of more systematic and mature system.However,multi-scale modeling and simulation can promote and deepen our understanding of the microstructure and its deformation mechanism.In this paper,we review and summarize the recent research progress of numerical simulation of Mg alloys in forming and microstructure,namely casting,extrusion,rolling,and welding,using crystal plasticity finite element(CPFEM)and molecular dynamics(DM)methods.Besides,the methods and innovations of modeling are also summarized.Lastly,the paper discusses the development prospects and challenges of the numerical simulation in the field of Mg alloys.
基金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.
基金The National Natural Science Foundation of China(Grant#52278161)the Science and Technology Project of Guangzhou(Grant#2024A04J9888)the Guangdong Basic and Applied Basic Research Foundation(Grant#2023A1515010535).
文摘This study investigates the flexural performance of ultra-high performance concrete(UHPC)in reinforced concrete T-beams,focusing on the effects of interfacial treatments.Three concrete T-beam specimens were fabricated and tested:a control beam(RC-T),a UHPC-reinforced beam with a chiseled interface(UN-C-50F),and a UHPC-reinforced beam featuring both a chiseled interface and anchored steel rebars(UN-CS-50F).The test results indicated that both chiseling and the incorporation of anchored rebars effectively created a synergistic combination between the concrete T-beam and the UHPC reinforcement layer,with the UN-CS-50F exhibiting the highest flexural resistance.The cracking load and ultimate load of UN-CS-50F were 221.5%and 40.8%,respectively,higher than those of the RC-T.Finite element(FE)models were developed to provide further insights into the behavior of the UHPCreinforced T-beams,showing a maximumdeviation of just 8%when validated against experimental data.A parametric analysis varied the height,thickness,andmaterial strength of the UHPC reinforcement layer based on the validated FE model,revealing that increasing the UHPC layer thickness from 30 to 50 mm improved the ultimate resistance by 20%while reducing the UHPC reinforcement height from 440 to 300 mm led to a 10%decrease in bending resistance.The interfacial anchoring rebars significantly reduced crack propagation and enhanced stress redistribution,highlighting the importance of strengthening interfacial bonds and optimizing geometric parameters ofUHPCfor improved T-beam performance.These findings offer valuable insights for the design and retrofitting of UHPC-reinforced bridge girders.
基金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.
基金supported by the National Natural Science Foundation of China(Grant Nos.52171318 and 12202329)Joint Foundation of the Ministry of Education(Grant No.8091B022105)。
文摘The detonation of fuel-rich explosives yields combustible products that persistently burn upon mixing with ambient oxygen,releasing additional energy through a phenomenon known as the afterburning effect.This process greatly influences the evolution of confined blast loading and the subsequent structural response,which is crucial in confined blast scenarios.Given the complex nature of the reaction process,accurate analysis of the afterburning effect remains challenging.Previous studies have either overlooked the mechanisms of detonation product combustion or failed to provide experimental validation.This study introduces a three-dimensional model to effectively characterize the combustion of detonation products.The model integrates chemical reaction source terms into the governing equations to consider the combustion processes.Numerical simulations and experimental tests were conducted to analyze the combustion and energy release from the detonation products of fuel-rich explosives in confined spaces.Approximately 50%of the energy was released during the combustion of detonation products in a confined TNT explosion.Although the combustion of these products was much slower than the detonation process,it aligned with the dynamic response of the structure,which enhanced the explosive yield.Excluding afterburning from the analysis reduced the center-point deformation of the structure by 30%.Following the inclusion of afterburning,the simulated quasistatic pressure increased by approximately 45%.Subsequent comparisons highlighted the merits of the proposed approach over conventional methods.This approach eliminates the reliance on empirical parameters,such as the amount and rate of energy release during afterburning,thereby laying the foundation for understanding load evolution in more complex environments,such as ships,buildings,and underground tunnels.
基金supported by the National Natural Science Foundation of China(Grant Nos.42475016,42192555 and 42305085)the China Postdoctoral Science Foundation(Grant No.2023M741615)the 2023 Graduate Research Innovation Project of Hunan Province(Grant No.CX20230011)。
文摘This paper investigates the impact of the model top and damping layer on the numerical simulation of tropical cyclones(TCs)and reveals the significant role of stratospheric gravity waves(SGWs).TCs can generate SGWs,which propagate upward and outward into the stratosphere.These SGWs can reach the damping layer,which is a consequence of the numerical scheme employed,where they can affect the tangential circulation through the dragging and forcing processes.In models with a higher top boundary,this tangential circulation develops far from the TC and has minimal direct impact on TC intensity.By comparison,in models with a lower top(e.g.,20 km),the damping layer is located just above the top of the TC.The SGW dragging in the damping layer and the consequent tangential force can thus induce ascent outside the eyewall,promote latent heat release,tilt the eyewall,and enlarge the inner-core radius.This process will reduce inner-core vorticity advection within the boundary layer,and eventually inhibits the intensification of the TC.This suggests that when the thickness of the damping layer is 5 km,the TC numerical model top height should be at least higher than 20 km to generate more accurate simulations.
基金The National Natural Science Foundation of China under contract Nos 41976217 and 42306249the National Key Research and Development Program of China under contract No.2018YFA0605701.
文摘A three-dimensional coupled sea ice-ice shelf-ocean numerical model is developed for the Prydz Bay,Antarctica,using the Regional Ocean Modeling System with a grid resolution of approximately 2 km.The influence of the grounding giant iceberg D15 on the distribution of sea ice and polynyas in the Prydz Bay is analyzed through two numerical experiments.Iceberg D15,grounded off the western edge of the West Ice Shelf(WIS),obstructs the southwestward transport of sea ice along the east coast of Prydz Bay,causing sea ice to accumulate to the east of the iceberg and form multi-year fast ice.Grounding of Iceberg D15 also decreases sea ice coverage off its south edge and creates ice-free openings in spring near Davis Station and Zhongshan Station,facilitating the accessibility of vessels to the research stations.These simulated sea ice patterns closely match current satellite observations.When Iceberg D15 is removed,the previously blocked sea ice north of the iceberg,which moved westward,shifts southwesterly along the coastline,leading to a reduction in sea ice thickness during winter and spring,as well as lower sea ice concentrations in spring across large areas north and west of the iceberg.In contrast,the sea ice thickness increases considerably southwest of the WIS,extending to the front of the Amery Ice Shelf during seasons covered by sea ice.The increase in sea ice concentration can also extend to as far as 75°E in spring.Without Iceberg D15,which previously contributed to the ice barrier of Barrier Polynya(BP),the shape of BP changes,the area of BP and Davis Polynya(DP)decreases,and the polynya off the northwest edge of the WIS near 83°E expands.These polynya patterns are much similar to the satellite remote sensing observations before Iceberg D15 was grounded.From April to October,the total area of BP and DP decreases by 2.83×10^(4)km^(2)(60%)and 2.20×10^(3)km^(2)(20%),respectively,while the total sea ice production decreases by 4.11×10^(10)m^(3)(66%)and 1.52×10^(10)m^(3)(52%)compared to the experiment with iceberg.These results indicate the substantial effects of grounding giant icebergs on the spatio-temporal distribution of sea ice,the area of polynyas,and sea ice production.High-resolution Antarctic coastal numerical models,typically with grid scales of kilometers,are sufficient to represent large icebergs,and adding the grounding giant icebergs is crucial for producing realistic simulations of sea ice and polynyas.
文摘It has been experimentally observed that,in the perforation of metal plates by a flat-nosed projectile,there exists a plateau phenomenon where the ballistic limit increases slightly with increasing plate thickness,which is related to a change in the mode of failure.No theoretical model has so far explained this phenomenon satisfactorily.This paper presents a combined numerical and theoretical study on the perforation of 2024-T351 aluminum plates struck by flat-nosed projectiles.First,numerical simulations are performed to investigate the failure mechanisms/deformation modes of the aluminum plates.Then,a theoretical model is proposed based on the numerical results and the experimental observations within a unified framework.The model takes into account the main energy absorbing mechanisms and the corresponding energies absorbed are determined analytically.In particular,a dimensionless equation is suggested to describe the relationship between global deformations and impact velocity.It transpires that the model predictions are in good agreement with the test data and the numerical results for the perforation of 2024-T351 aluminum plates struck by rigid flat-nosed projectiles in terms of residual velocity,ballistic limit,relationship between global deformations and impact velocity,and transition of failure modes.It also transpires that the present model can predict the“plateau”phenomenon,which shows a slight increase in ballistic limit as plate thickness increases.Furthermore,the energy absorption mechanisms are discussed on the basis of the theoretical analysis.
基金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.
文摘Plastometric experiments,supplemented with numerical simulations using the finite element method(FEM),can be advantageously used to characterize the deformation behavior of metallic materials.The accuracy of such simulations predicting deformation behaviors of materials is,however,primarily affected by the applied rheology law.The presented study focuses on the characterization of the deformation behavior of AISI 1045 type carbon steel,widely used e.g.,in automotive and power engineering,under extreme conditions(i.e.,high temperatures,strain rates).The study consists of two main parts:experimentally analyzing the flow stress development of the steel under different thermomechanical conditions via uniaxial hot compression tests and establishing the rheology law via numerical simulations implementing the experimentally acquired flow stress curves.The numerical simulations then not only serve to establish the rheology law but also to verify the reliability of the selected experimental process.The results of the numerical simulations showed that the established rheology law characterizes the behavior of the investigated steel with sufficient accuracy also at high temperatures and/or strain rates,and can,therefore,be used for practical purposes.Last but not least,supplementary microstructure analyses performed for the samples subjected to the highest deformation temperature provided a deeper insight into the effects of the applied(extreme)thermomechanical conditions on the behavior of the investigated steel.
文摘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%.
基金funded by the National Natural Science Foundation of China(Nos.42077228,52174085)。
文摘With the increasing development of deepburied engineering projects,rockburst disasters have become a frequent concern.Studies have indicated that tunnel diameter is a critical factor influencing the occurrence of rockbursts.To investigate the influence of tunnel diameter on the deformation and failure characteristics of surrounding rock,large-sized rocklike gypsum specimens were tested using a selfdeveloped true triaxial rockburst loading system containing circular tunnels with three different diameters(D=0.07 m,0.11 m,and 0.15 m).Acoustic emission monitoring,together with a miniature intelligent camera,was employed to analyze the entire process,focusing on macroscopic failure patterns,fragment characteristics,and underlying failure mechanisms.In addition,theoretical analyses were carried out and combined with numerical simulations to investigate the differences in energy evolution associated with rockburst physical models.The results indicate that:(1)The rockburst process with different tunnel diameters consistently evolved through three distinct stages—initial particle ejection,crack propagation accompanied by flake spalling,and,finally,fragment ejection leading to the formation of a‘V'-shaped notch.(2)Increasing tunnel diameter reduces rockburst failure load while increasing surrounding rock damage extent,total mass and average size of ejected fragments.Additionally,shear failure proportion decreases with tensile failure becoming increasingly dominant.(3)Larger tunnel diameters reduce the attenuation rate of elastic strain energy,thereby expanding the zone of elastic strain energy accumulation and disturbance and creating conditions for larger volume rockburst.(4)Larger tunnel diameters result in a smaller principal stress ratio at equivalent distances in the surrounding rock,indicating a higher likelihood of tensile failure.(5)Numerical analyses further reveal that larger tunnel diameters reduce the maximum elastic strain energy density around the tunnel,lowering the energy released per unit volume of rockburst fragments and their ejection velocities.However,both the total failure volume and overall energy release from rockburst increase.Model experiments with different tunnel diameters are of great significance for optimizing engineering design and parameter selection,as well as guiding tunnel construction under complex geological conditions.
基金supported by the Beijing Postdoctoral Research Foundation(No.2023-ZZ-133)Scientific Research Foundation of Beijing Infrastructure Investment Co.,Ltd.(No.2023-ZB-03)Fundamental Research Funds for the Central Universities(No.2682023ZTPY036).
文摘The pressure wave generated by an urban-rail vehicle when passing through a tunnel affects the comfort of passengersand may even cause damage to the train and related tunnel structures.Therefore,controlling the trainspeed in the tunnel is extremely important.In this study,this problem is investigated numerically in the frameworkof the standard k-εtwo-equation turbulence model.In particular,an eight-car urban rail train passingthrough a tunnel at different speeds(140,160,180 and 200 km/h)is considered.The results show that the maximumaerodynamic drag of the head and tail cars is most affected by the running speed.The pressure at selectedmeasuring points on the windward side of the head car is very high,and the negative pressure at the side windowof the driver’s cab of the tail car is also very large.From the head car to the tail car,the pressure at the same heightgradually decreases.The positive pressure peak at the head car and the negative pressure peak at the tail car aregreatly affected by the speed.
