Landslides triggered by seismic activity have led to substantial human and economic losses.Nevertheless,the fundamental physical mechanisms underlying the vibration and rupture of rock slopes during earthquakes remain...Landslides triggered by seismic activity have led to substantial human and economic losses.Nevertheless,the fundamental physical mechanisms underlying the vibration and rupture of rock slopes during earthquakes remain poorly understood.In this study,finite element method-based numerical simulations were conducted based on the rock slope at Dagangshan Hydropower Station in Sichuan province,China.Firstly,systematic analysis in both the time and frequency domains were performed to examine the seismic dynamic characteristics of the slope.Subsequently,the transfer function method and the multiple stepwise linear regression method were employed to clarify the underlying mechanism and determine critical factors influencing the slope instability during earthquakes.Time-domain analysis reveals that rock slope dynamic response exhibits notable elevation,surface,and local amplification effects.Specifically,the Peak Ground Acceleration(PGA)amplification coefficient(MPGA)is significantly higher at elevated locations,near the slope surface and in areas with protrusions.Moreover,the existence of fracture zones and anti-shear galleries minimally influences the dynamic responses but considerably affect the rupture.Specifically,fracture zones exacerbate rupture,while anti-shear galleries mitigate it.Frequency-domain analysis indicates that the dynamic responses of the slope are closely correlated with the degree of slope rupture.As earthquake magnitude increases,the rupture degree of the slope intensifies,and the dominant frequency of the response within the slope decreases,e.g.,its value shifts from 3.63 to 2.75 Hz at measurement point 9near the slope surface.The transfer function of rock slope,calculated under the excitation of wide flat spectrum white noise can reflect the interrelationships between the inherent properties and the rupture degree.Notably,the peak of the transfer function undergoes inversion as the degree of rupture increases.Furthermore,through multiple stepwise linear regression analysis,four key factors influencing the surface dynamic response of the slope were identified:rock strength,slope angle,elevation,and seismic dominant frequency.These findings provide valuable insights into the underlying mechanisms of rock slope dynamic responses triggered by earthquakes,offering essential guidance for understanding and mitigating seismic impacts on rock slopes.展开更多
Experimental studies were conducted on two high-strength steel plate-frame structures with different truss spacings under various impact velocities to investigate the dynamic mechanical properties of hull plate-frame ...Experimental studies were conducted on two high-strength steel plate-frame structures with different truss spacings under various impact velocities to investigate the dynamic mechanical properties of hull plate-frame structures under drop weight impact.The results showed that decreasing the main beam spacing can effectively increase the structural stiffness,reduce the maximum deformation,and increase the damage range.Furthermore,to simulate the impact tests accurately,static and dynamic tensile tests at different strain rates were carried out,and the Cowper-Symonds model parameters were fitted via experimental data.The material properties obtained from the tensile tests were used as inputs for numerical simulations with the numerical results coincide with the experimental results.A systematic analysis and discussion were conducted on the effects of truss spacing and truss width on the dynamic response of the reinforced plates,and an optimal range for the ratio of truss spacing to truss width was proposed.In addition,a mesh size sensitivity analysis for ship hull plate frame collision simulations was performed.The applicability of the EPS,MMC,and RTCL failure criteria in the simulation of plate-frame structures was investigated via finite element simulations of falling weight impact tests.The research findings provide a reference for ship hull structure design and resilience assessment.展开更多
The turning performance of a ship is an important aspect of its maneuverability,and accurately predicting the hydrodynamic forces during ship turning motion is of great significance for the safe maneuvering design of ...The turning performance of a ship is an important aspect of its maneuverability,and accurately predicting the hydrodynamic forces during ship turning motion is of great significance for the safe maneuvering design of ships.This paper investigated the hydrodynamic performance of a KRISO container ship in steady turning using experimental and numerical approaches.The rotating arm tests were carried out in rotating arm basin of Zhejiang University,while the numerical simulations were conducted in commercial computational fluid dynamics software.Hydrodynamic forces and moments,hull surface wave height,wave patterns,and vorticity are studied under different velocities,radii,and drift angles.The results show that the increase in velocity has a significant impact on the forces and moments of the hull.The changes in longitudinal and transverse forces reflect the complex fluid dynamic interactions between the hull and water.Under conditions of small radius and large drift angle,the hull experiences greater forces and moments,indicating that stability and maneuverability will be more challenged during sudden turns.This study can provide experimental data and numerical simulation references for the research of ship turning maneuvers.展开更多
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.展开更多
0 INTRODUCTION In recent years,modern railways have been actively under construction in the complex mountainous area of Southwest China.However,rockfall poses a significant threat to both construction and operation ph...0 INTRODUCTION In recent years,modern railways have been actively under construction in the complex mountainous area of Southwest China.However,rockfall poses a significant threat to both construction and operation phases of railway projects(Yan et al.,2023;Chen et al.,2022;Fanos and Pradhan,2018).展开更多
The flexural strength of glass is a critical design parameter for applications encountering impact loadings.However,the micro defects,specimen geometry,loading rate,and load transformation from a quasi-dynamic to quas...The flexural strength of glass is a critical design parameter for applications encountering impact loadings.However,the micro defects,specimen geometry,loading rate,and load transformation from a quasi-dynamic to quasi-impulsive state may influence the measurement accuracy.Due to the stochastic and amorphous nature of the material,an accurate determination of the flexural strength remains a challenge.In this two-fold study,a coupled experimental-numerical strategy was devised to evaluate the dynamic flexural strength.In the first phase,three-point bending experiments were conducted on a novel“Electromagnetic Split Hopkinson Pressure Bar(ESHPB)”.The incident stress signal and fracture time were recorded from experimental data,while the flexural strength was indirectly computed from a numerical algorithm.A quantitative comparison of the flexural strength with those in existing literature established the accuracy of the proposed methodology.Results of the study indicate that the specimen response became independent of the support conditions under impulsive loading.That being said,the specimen behaved like it had an infinite span length,and the measured flexural strength remained the same whether the specimen was supported or not.Besides,the specimen also maintained contact at the interfaces of the incident bar and fixture supports for the entire loading duration.In the second part of this study,the computed flexural strength was used to calibrate the existing JH-2 model.Numerical prediction of the damage propagation corroborated with that obtained from reprography images,though qualitatively.This work presents a precise and robust methodology to determine the dynamic flexural strength of brittle ceramics like Aluminosilicate glass over traditional experimental procedures to facilitate its adoption.展开更多
This study aims to elucidate the dynamic evolution mechanism of the fracturing fracture system during the exploration and development of complex oil and gas reservoirs.By integrating methods of rock mechanical testing...This study aims to elucidate the dynamic evolution mechanism of the fracturing fracture system during the exploration and development of complex oil and gas reservoirs.By integrating methods of rock mechanical testing,logging calculation,and seismic inversion technology,we obtained the current insitu stress characteristics of a single well and rock mechanical parameters.Simultaneously,significant controlling factors of rock mechanical properties were analyzed.Subsequently,by coupling hydraulic fracturing physical experiments with finite element numerical simulation,three different fracturing models were configured:single-cluster,double-cluster,and triple-cluster perforations.Combined with acoustic emission technology,the fracture initiation mode and evolution characteristics during the loading process were determined.The results indicate the following findings:(1)The extension direction and length of the fracture are significantly controlled by the direction of the maximum horizontal principal stress.(2)Areas with poor cementation and compactness exhibit complex fracture morphology,prone to generating network fractures.(3)The interlayer development of fracturing fractures is controlled by the strata occurrence.(4)Increasing the displacement of fracturing fluid enlarges the fracturing fracture length and height.This research provides theoretical support and effective guidance for hydraulic fracturing design in tight oil and gas reservoirs.展开更多
Differential interferometric synthetic aperture radar （DInSAR） technology is a new method to monitor the dynamic surface subsidence. It can monitor the large scope of dynamic deformation process of surface subsidenc...Differential interferometric synthetic aperture radar （DInSAR） technology is a new method to monitor the dynamic surface subsidence. It can monitor the large scope of dynamic deformation process of surface subsidence basin and better reflect the surface subsidence form in different stages. But under the influence of factors such as noise and other factors, the tilt and horizontal deformation curves regularity calculated by DInSAR data are poorer and the actual deviation is larger. The tilt and horizontal deformations are the important indices for the safety of surface objects protection. Numerical simulation method was used to study the dynamic deformation of LW32 of West Cliff colliery in Australia based on the DInSAR monitoring data. The result indicates that the subsidence curves of two methods fit well and the correlation coefficient is more than 95%. The other deformations calculated by numerical simulation results are close to the theory form. Therefore, considering the influence, the surface and its subsidiary structures and buildings due to mining, the numerical simulation method based on the DInSAR data can reveal the distribution rules of the surface dynamic deformation values and supply the shortcomings of DInSAR technology. The research shows that the method has good applicability and can provide reference for similar situation.展开更多
A mathematical model of principal elements of the aircraft hydraulic system is presented based on the heat transfer theory. The dynamic heat transfer process of the hydraulic oil and the pump shells within an aircraft...A mathematical model of principal elements of the aircraft hydraulic system is presented based on the heat transfer theory. The dynamic heat transfer process of the hydraulic oil and the pump shells within an aircraft hydraulic system are analyzed by the difference method. A kind of means for the prediction to variational trends of the aircraft hydraulic system temperature is provided during operation. The numerical prediction and simulation under the operational conditions are presented for ground trial running and the decelerated operation in flight. Computational results show that there is a good coincidence between the experimental data and the numerical predictions.展开更多
Numerical modelling is an effective technique to improve the understanding of outburst initiation mechanisms and to take appropriate measures to address their threats.Based on the existing two-way sequential coupling ...Numerical modelling is an effective technique to improve the understanding of outburst initiation mechanisms and to take appropriate measures to address their threats.Based on the existing two-way sequential coupling method,two typical types of outbursts,i.e.the gas pocket outburst and the dynamic fracturing outburst,have been successfully simulated using field data from a coalfield in central China.The geological structure commonly observed in the coalfield,known as the‘bedding shear zone’,contributes to the gas pocket outbursts in the region.The model for this type of outburst simulates mininginduced stress and gas pressure distributions during the outburst initiation stage and the subsequent development stage.Both coal ejection and gas release are observed in the model,and the simulation results are consistent with mine site observations,i.e.the amount of ejected coal,outburst cavity profile,and gas release rate change prior to an outburst.The second type of outburst is attributed to gas accumulation and elevated gas pressure due to the gassy floor seam and the heterogeneity in the floor strata,which is explained by the dynamic fracturing theory.While the dynamic coal ejection phenomenon is not captured in the simulation,the abrupt release of retained gas from a floor coal seam is successfully replicated.Both outburst models reveal that abnormal gas emission trends can be used as indicators of an upcoming outburst.The results of this study are expected to provide new insights into the outburst initiation mechanisms and outburst prevention measures.展开更多
The dynamic performance of high-speed trains is significantly influenced by sudden changes in aerodynamic loads(ADLs)when exiting a tunnel in a windy environment.Focusing on a double-track tunnel under construction in...The dynamic performance of high-speed trains is significantly influenced by sudden changes in aerodynamic loads(ADLs)when exiting a tunnel in a windy environment.Focusing on a double-track tunnel under construction in a mountain railway,we established an aerodynamic model involving a train exiting the tunnel,and verified it in the Fluent environment.Overset mesh technology was adopted to characterize the train’s movement.The flow field involving the train,tunnel,and crosswinds was simulated using the Reynolds-averaged turbulence model.Then,we built a comprehensive train-track coupled dynamic model considering the influences of ADLs,to investigate the vehicles’dynamic responses.The aerodynamics and dynamic behaviors of the train when affected by crosswinds with different velocities and directions are analyzed and discussed.The results show that the near-wall side crosswind leads to sharper variations in ADLs than the far-wall side crosswind.The leading vehicle suffers from more severe ADLs than other vehicles,which worsens the wheel-rail interaction and causes low-frequency vibration of the car body.When the crosswind velocity exceeds 20 m/s,significant wheel-rail impacts occur,and the running safety of the train worsens rapidly.展开更多
The fragile and intricate geological environment of the Qinghai-Tibet Plateau gives rise to numerous precarious rocks along the riverbanks,posing significant risks for the upcoming construction of hydropower stations....The fragile and intricate geological environment of the Qinghai-Tibet Plateau gives rise to numerous precarious rocks along the riverbanks,posing significant risks for the upcoming construction of hydropower stations.In order to identify potential rockfalls that could endanger the Zixia hydropower project,a comprehensive analysis employing various methods was conducted to investigate the kinematic characteristics and dynamic fragmentation of such precarious rocks.Initially,UAV oblique photography and field survey were used to create a digital elevation model with a resolution of 0.25 m and map the spatial distribution of precarious rocks.Subsequently,the development characteristics of joints within rock masses were analyzed through an adit investigation.Following these preliminary steps,a transportation simulation utilizing RocPro3D,considering stochastic initiation orientation,was employed to predict the trajectories of 18 precarious rocks.As a result,two hazardous rocks that pose a direct threat to the cofferdam were identified.Finally,considering the influence of internal structure planes,a discrete element method was applied for accurately simulating the kinematic characteristics and dynamic fragmentation of these hazardous rocks.The findings underscore several key observations:(1)Slopeparallel structure planes within these hazardous rocks play a pivotal role in both the progressive failure during initiation and dynamic fragmentation during transportation;(2)Hazardous rocksⅢ-1 andⅣ-1 would pose a direct threat to the cofferdam.Notably,block b4 from hazardous rockⅢ-1,could potentially impact the cofferdam with an energy of 4598.65 kJ and an impact force of 3007.5 kN;and(3)Continuous collisions encountered during transportation facilitate the disintegration of rock masses along structure planes and generate substantial high-velocity fragments.Finally,to cope with the impact risk of collapsing blocks,a reinforced retaining wall as the mitigation measure is recommended.展开更多
This study analyzes Brazilian stromatolites in Lagoa Salgada,serving as analogs for pre-salt rocks in the Santos and Campos basins.Despite their excellent petrophysical properties,such as high porosity and permeabilit...This study analyzes Brazilian stromatolites in Lagoa Salgada,serving as analogs for pre-salt rocks in the Santos and Campos basins.Despite their excellent petrophysical properties,such as high porosity and permeability,these reservoirs present challenges in fluid flow modeling and simulation.The research investigates various factors influencing the development of carbonate reservoirs,including diagenetic processes employing several techniques,such as microcomputed tomography(micro-CT)and digital rock physics(DRP),to study petrophysical and geological characteristics.Additionally,through numerical simulations,the properties of fluid flow in different microfacies of stromatolites are estimated,with particular emphasis on understanding and highlighting changes in the direction of fluid flow in the three characterized microfacies.These findings offer crucial insights into optimizing oil and gas exploration and production techniques in carbonate reservoirs,providing a comprehensive understanding of the dynamics of fluid transport in porous media,especially in terms of directional changes within stromatolites.展开更多
In this study,four types of spiral fins with varying parameters were mounted on an upstream cylinder,and the effects of spiral fins on the vibration response of heat exchange tubes and the vortex structure in cross fl...In this study,four types of spiral fins with varying parameters were mounted on an upstream cylinder,and the effects of spiral fins on the vibration response of heat exchange tubes and the vortex structure in cross flow were studied through experiments and numerical simulations.The results indicate a strong dependency of the cylinder's vibration response on the fin parameters.The results indicate that the vibration response and wake structure of the cylinder are significantly influenced by the parameters of the fins.The introduction of a finned cylinder affects both its own vibration amplitude and frequency,as well as the downstream cylinder.The amplitudes of finned cylinders Ⅰ and Ⅲ are reduced by 57.8% and 59.9%,respectively,compared to the bare cylinder.This reduction helps to restrain vibration and diminishes the amplitudes of the downstream cylinder.Although finned cylinder Ⅱ slightly decreases its own vibration,it increases the amplitude of the downstream cylinder by 13.7%.The mean drag coefficient and the root mean square of the lift coefficient of the finned cylinder are higher than those of the bare cylinder when the finned cylinder is positioned upstream.Smaller pitch and larger equivalent diameter will lead to increased drag,resulting in enhanced vortex shedding in the wake,which amplifies the vibrations of the cylinder in that wake.The downstream of finned cylinder Ⅱ has the widest wake and higher vortex strength,and the dynamic load and vibration of the downstream cylinder are increased.The vortex intensity decays faster in the wake of finned cylinder Ⅲ,and the vibration of the downstream cylinder is weaker.展开更多
In this study,a new linear friction welding(LFW)process,embedded LFW process,was put forward,which was mainly applied to combination manufacturing of long or overlong loadcarrying titanium alloy structural components ...In this study,a new linear friction welding(LFW)process,embedded LFW process,was put forward,which was mainly applied to combination manufacturing of long or overlong loadcarrying titanium alloy structural components in aircraft.The interfacial plastic flow behavior and bonding mechanism of this process were investigated by a developed coupling EulerianLagrangian numerical model using software ABAQUS and a novel thermo-physical simulation method with designed embedded hot compression specimen.In addition,the formation mechanism and control method of welding defects caused by uneven plastic flow were discussed.The results reveal that the plastic flow along oscillating direction of this process is even and sufficient.In the direction perpendicular to oscillation,thermo-plastic metals mainly flow downward along welding interface under coupling of shear stress and interfacial pressure,resulting in the interfacial plastic zone shown as an inverted“V”shape.The upward plastic flow in this direction is relatively weak,and only a small amount of flash is extruded from top of joint.Moreover,the wedge block and welding components at top of joint are always in un-steady friction stage,leading to nonuniform temperature field distribution and un-welded defects.According to the results of numerical simulation,high oscillating frequency combined with low pressure and small amplitude is considered as appropriate parameter selection scheme to improve the upward interfacial plastic flow at top of joint and suppress the un-welded defects.The results of thermo-physical simulation illustrate that continuous dynamic recrystallization(CDRX)induces the bonding of interface,accompanying by intense dislocation movement and creation of many low-angle grain boundaries.In the interfacial bonding area,grain orientation is random with relatively low texture density(5.0 mud)owing to CDRX.展开更多
Semisubmersible naval ships are versatile military crafts that combine the advantageous features of high-speed planing crafts and submarines.At-surface,these ships are designed to provide sufficient speed and maneuver...Semisubmersible naval ships are versatile military crafts that combine the advantageous features of high-speed planing crafts and submarines.At-surface,these ships are designed to provide sufficient speed and maneuverability.Additionally,they can perform shallow dives,offering low visual and acoustic detectability.Therefore,the hydrodynamic design of a semisubmersible naval ship should address at-surface and submerged conditions.In this study,Numerical analyses were performed using a semisubmersible hull form to analyze its hydrodynamic features,including resistance,powering,and maneuvering.The simulations were conducted with Star CCM+version 2302,a commercial package program that solves URANS equations using the SST k-ωturbulence model.The flow analysis was divided into two parts:at-surface simulations and shallowly submerged simulations.At-surface simulations cover the resistance,powering,trim,and sinkage at transition and planing regimes,with corresponding Froude numbers ranging from 0.42 to 1.69.Shallowly submerged simulations were performed at seven different submergence depths,ranging from D/LOA=0.0635 to D/LOA=0.635,and at two different speeds with Froude numbers of 0.21 and 0.33.The behaviors of the hydrodynamic forces and pitching moment for different operation depths were comprehensively analyzed.The results of the numerical analyses provide valuable insights into the hydrodynamic performance of semisubmersible naval ships,highlighting the critical factors influencing their resistance,powering,and maneuvering capabilities in both at-surface and submerged conditions.展开更多
Rock masses are often exposed to dynamic loads such as earthquakes and mechanical disturbances in practical engineering scenarios.The existence of underground caverns and weak geological structures like columnar joint...Rock masses are often exposed to dynamic loads such as earthquakes and mechanical disturbances in practical engineering scenarios.The existence of underground caverns and weak geological structures like columnar jointed rock masses(CJRMs)and interlayer shear weakness zones(ISWZs)with inferior mechanical properties,significantly undermines the overall structural stability.To tackle the dynamic loading issues in the process of constructing subterranean caverns,a programmable modeling approach was utilized to reconstruct a large-scale underground cavern model incorporating ISWZs and columnar joints(CJs).By conducting dynamic simulations with varying load orientations,the analyses focused on the failure patterns,deformation characteristics,and acoustic emission activity within the caverns.Results revealed that the failure modes of the underground caverns under dynamic loading were predominantly tensile failures.Under X-direction loading,the failed elements were mainly distributed parallel to the CJs,while under Y-direction loading,they were distributed parallel to the transverse weak structural planes.Furthermore,the dynamic stability of the overall structure varied with the number of caverns.The dual-cavern model demonstrated the highest stability under X-direction loading,while the single-cavern model was the least stable.Under Y-direction loading,the cavern stability increased with the number of caverns.Importantly,different weak structures affected the dynamic response of caverns in different ways;the CJRMs were the primary contributors to structural failure,while ISWZs could mitigate the rock mass failure induced by CJs.The findings could offer valuable insights for the dynamic stability analysis of caverns containing CJRMs and ISWZs.展开更多
This numerical study of the Bohai Sea wintertime circulation by means of a two dimensional barotropic model with resolution of 1/24° in longitude and latitude showed that the Bohai Sea wintertime circulation is d...This numerical study of the Bohai Sea wintertime circulation by means of a two dimensional barotropic model with resolution of 1/24° in longitude and latitude showed that the Bohai Sea wintertime circulation is dominated by local monsoon winds. The major current components include the Bohai Warm Current, the North Shandong Coastal Current, and the Liaodong Gyre. The Bohai Warm Current originates from the Yellow Sea Warm Current at the northern part of Bohai Strait, meanders westwards and finally enters the northern part of Bohai Bay. The North Shandong Coastal Current flows along the southwest shore of Bohai Bay and Laizhou Bay and exits from the Bohai Sea through the south Bohai Strait. The anticyclonic Liaodong Gyre is located in the north of Liaodong Bay. A pair of eddies and the small scale Jinzhou Gyre are found between the Bohai Warm Current and the Liaodong Gyre. The computed volume transport for both the Bohai Warm Current and North Shandong Coastal Current is about 0.03 Sv (1 Sv=1×10 6 m 3/s). The numerical experiments showed that the combined effect of local monsoon winds and bottom topography dominate the formation of the circulation pattern. The Coriolis force and the wind stress curl are of certain importance. The beta effect, the momentum advection and the open boundary condition have little influence on the circulation pattern.展开更多
To evaluate the effects of possible ground explosion on a shallow-buried metro tunnel, this paper attempts to analyze the dynamic responses of the operating metro tunnel in soft soil, using a widely applied explicit d...To evaluate the effects of possible ground explosion on a shallow-buried metro tunnel, this paper attempts to analyze the dynamic responses of the operating metro tunnel in soft soil, using a widely applied explicit dynamic nonlinear finite element software ANSYS/LS-DYNA. The blast induced wave propagation in the soil and the tunnel, and the von Mises effective stress and acceleration of the tunnel lining were presented, and the safety of the tunnel lining was evaluated based on the failure criterion. Besides, the parametric study of the soil was also carried out. The numerical results indicate that the upper part of the tunnel lining cross-section with directions ranging from 0° to 22.5° and horizontal distances 0 to 7 m away from the explosive center are the vulnerable areas, and the metro tunnel might be safe when tunnel depth is more than 7 m and TNT charge on the ground is no more than 500 kg, and the selection of soil parameters should be paid more attentions to conduct a more precise analysis.展开更多
In this paper, a full-scale 3-D finite element model of the Jundushan cable-stayed aqueduct bridge is established with ANSYS Code. The shell, fluid, tension-only spar and beam elements are used for modeling the aquedu...In this paper, a full-scale 3-D finite element model of the Jundushan cable-stayed aqueduct bridge is established with ANSYS Code. The shell, fluid, tension-only spar and beam elements are used for modeling the aqueduct deck, filled water, cables and support towers, respectively. A multi-element cable formulation is introduced to simulate the cable vibration. The dry （without water） and wet （with water） modes of the aqueduct bridge are both extracted and investigated in detail. The dry modes of the aqueduct bridge are basically similar to those of highway cable-stayed bridges. A dry mode may correspond to two types of wet modes, which are called the in-phase （with lower frequency） and out-of-phase （with higher frequency） modes. When the water-structure system vibrates in the in-phase/out-of-phase modes, the aqueduct deck moves and water sloshes in the same/opposite phase-angle, and the sloshing water may take different surface-wave modes. The wet modes of the system reflect the properties of interaction among the deck, towers, cables and water. The in-phase wet frequency generally decreases as the water depth increases, and the out-of-phase wet frequency may increase or decrease as the water depth increases.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.52274075,42122052,52379098)。
文摘Landslides triggered by seismic activity have led to substantial human and economic losses.Nevertheless,the fundamental physical mechanisms underlying the vibration and rupture of rock slopes during earthquakes remain poorly understood.In this study,finite element method-based numerical simulations were conducted based on the rock slope at Dagangshan Hydropower Station in Sichuan province,China.Firstly,systematic analysis in both the time and frequency domains were performed to examine the seismic dynamic characteristics of the slope.Subsequently,the transfer function method and the multiple stepwise linear regression method were employed to clarify the underlying mechanism and determine critical factors influencing the slope instability during earthquakes.Time-domain analysis reveals that rock slope dynamic response exhibits notable elevation,surface,and local amplification effects.Specifically,the Peak Ground Acceleration(PGA)amplification coefficient(MPGA)is significantly higher at elevated locations,near the slope surface and in areas with protrusions.Moreover,the existence of fracture zones and anti-shear galleries minimally influences the dynamic responses but considerably affect the rupture.Specifically,fracture zones exacerbate rupture,while anti-shear galleries mitigate it.Frequency-domain analysis indicates that the dynamic responses of the slope are closely correlated with the degree of slope rupture.As earthquake magnitude increases,the rupture degree of the slope intensifies,and the dominant frequency of the response within the slope decreases,e.g.,its value shifts from 3.63 to 2.75 Hz at measurement point 9near the slope surface.The transfer function of rock slope,calculated under the excitation of wide flat spectrum white noise can reflect the interrelationships between the inherent properties and the rupture degree.Notably,the peak of the transfer function undergoes inversion as the degree of rupture increases.Furthermore,through multiple stepwise linear regression analysis,four key factors influencing the surface dynamic response of the slope were identified:rock strength,slope angle,elevation,and seismic dominant frequency.These findings provide valuable insights into the underlying mechanisms of rock slope dynamic responses triggered by earthquakes,offering essential guidance for understanding and mitigating seismic impacts on rock slopes.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.52171311 and 5227127).
文摘Experimental studies were conducted on two high-strength steel plate-frame structures with different truss spacings under various impact velocities to investigate the dynamic mechanical properties of hull plate-frame structures under drop weight impact.The results showed that decreasing the main beam spacing can effectively increase the structural stiffness,reduce the maximum deformation,and increase the damage range.Furthermore,to simulate the impact tests accurately,static and dynamic tensile tests at different strain rates were carried out,and the Cowper-Symonds model parameters were fitted via experimental data.The material properties obtained from the tensile tests were used as inputs for numerical simulations with the numerical results coincide with the experimental results.A systematic analysis and discussion were conducted on the effects of truss spacing and truss width on the dynamic response of the reinforced plates,and an optimal range for the ratio of truss spacing to truss width was proposed.In addition,a mesh size sensitivity analysis for ship hull plate frame collision simulations was performed.The applicability of the EPS,MMC,and RTCL failure criteria in the simulation of plate-frame structures was investigated via finite element simulations of falling weight impact tests.The research findings provide a reference for ship hull structure design and resilience assessment.
基金supported by the China Scholarship Council(Grant No.202306320084).
文摘The turning performance of a ship is an important aspect of its maneuverability,and accurately predicting the hydrodynamic forces during ship turning motion is of great significance for the safe maneuvering design of ships.This paper investigated the hydrodynamic performance of a KRISO container ship in steady turning using experimental and numerical approaches.The rotating arm tests were carried out in rotating arm basin of Zhejiang University,while the numerical simulations were conducted in commercial computational fluid dynamics software.Hydrodynamic forces and moments,hull surface wave height,wave patterns,and vorticity are studied under different velocities,radii,and drift angles.The results show that the increase in velocity has a significant impact on the forces and moments of the hull.The changes in longitudinal and transverse forces reflect the complex fluid dynamic interactions between the hull and water.Under conditions of small radius and large drift angle,the hull experiences greater forces and moments,indicating that stability and maneuverability will be more challenged during sudden turns.This study can provide experimental data and numerical simulation references for the research of ship turning maneuvers.
基金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 Open Research Fund of Key Laboratory of Geological Hazards on Three Gorges Reservoir Area(China Three Gorges University),Ministry of Education(No.2022KDZ03)the Science and Technology Projects of Yunnan Provincial Science and Technology Department(No.202401AT070328)+1 种基金the Young talents project of“Xingdian Talent Support Program”in Yunnan Province(No.YNWR-QNBJ-2020-019)the Fund Project of China Academy of Railway Sciences Co.,Ltd.(No.2021YJ178)。
文摘0 INTRODUCTION In recent years,modern railways have been actively under construction in the complex mountainous area of Southwest China.However,rockfall poses a significant threat to both construction and operation phases of railway projects(Yan et al.,2023;Chen et al.,2022;Fanos and Pradhan,2018).
基金funded by the 111 Project,China(No.BP0719007)the National Key R&D Program of China(No.2017YFB1103500)the National Natural Science Foundation of China(Nos.11772268,12002178 and 12025205).
文摘The flexural strength of glass is a critical design parameter for applications encountering impact loadings.However,the micro defects,specimen geometry,loading rate,and load transformation from a quasi-dynamic to quasi-impulsive state may influence the measurement accuracy.Due to the stochastic and amorphous nature of the material,an accurate determination of the flexural strength remains a challenge.In this two-fold study,a coupled experimental-numerical strategy was devised to evaluate the dynamic flexural strength.In the first phase,three-point bending experiments were conducted on a novel“Electromagnetic Split Hopkinson Pressure Bar(ESHPB)”.The incident stress signal and fracture time were recorded from experimental data,while the flexural strength was indirectly computed from a numerical algorithm.A quantitative comparison of the flexural strength with those in existing literature established the accuracy of the proposed methodology.Results of the study indicate that the specimen response became independent of the support conditions under impulsive loading.That being said,the specimen behaved like it had an infinite span length,and the measured flexural strength remained the same whether the specimen was supported or not.Besides,the specimen also maintained contact at the interfaces of the incident bar and fixture supports for the entire loading duration.In the second part of this study,the computed flexural strength was used to calibrate the existing JH-2 model.Numerical prediction of the damage propagation corroborated with that obtained from reprography images,though qualitatively.This work presents a precise and robust methodology to determine the dynamic flexural strength of brittle ceramics like Aluminosilicate glass over traditional experimental procedures to facilitate its adoption.
基金supported by the Major Scientific and Technological Projects of CNPC under grant ZD2019-183-006the National Science and Technology Major Project of China(2016ZX05014002-006)the National Natural Science Foundation of China(42072234)。
文摘This study aims to elucidate the dynamic evolution mechanism of the fracturing fracture system during the exploration and development of complex oil and gas reservoirs.By integrating methods of rock mechanical testing,logging calculation,and seismic inversion technology,we obtained the current insitu stress characteristics of a single well and rock mechanical parameters.Simultaneously,significant controlling factors of rock mechanical properties were analyzed.Subsequently,by coupling hydraulic fracturing physical experiments with finite element numerical simulation,three different fracturing models were configured:single-cluster,double-cluster,and triple-cluster perforations.Combined with acoustic emission technology,the fracture initiation mode and evolution characteristics during the loading process were determined.The results indicate the following findings:(1)The extension direction and length of the fracture are significantly controlled by the direction of the maximum horizontal principal stress.(2)Areas with poor cementation and compactness exhibit complex fracture morphology,prone to generating network fractures.(3)The interlayer development of fracturing fractures is controlled by the strata occurrence.(4)Increasing the displacement of fracturing fluid enlarges the fracturing fracture length and height.This research provides theoretical support and effective guidance for hydraulic fracturing design in tight oil and gas reservoirs.
基金Project (20110023110014) supported by Specialized Research Fund for the Doctoral Program of Higher Education of ChinaProject (2010QD01) supported by Fundamental Research Funds for the Central Universities,China
文摘Differential interferometric synthetic aperture radar （DInSAR） technology is a new method to monitor the dynamic surface subsidence. It can monitor the large scope of dynamic deformation process of surface subsidence basin and better reflect the surface subsidence form in different stages. But under the influence of factors such as noise and other factors, the tilt and horizontal deformation curves regularity calculated by DInSAR data are poorer and the actual deviation is larger. The tilt and horizontal deformations are the important indices for the safety of surface objects protection. Numerical simulation method was used to study the dynamic deformation of LW32 of West Cliff colliery in Australia based on the DInSAR monitoring data. The result indicates that the subsidence curves of two methods fit well and the correlation coefficient is more than 95%. The other deformations calculated by numerical simulation results are close to the theory form. Therefore, considering the influence, the surface and its subsidiary structures and buildings due to mining, the numerical simulation method based on the DInSAR data can reveal the distribution rules of the surface dynamic deformation values and supply the shortcomings of DInSAR technology. The research shows that the method has good applicability and can provide reference for similar situation.
文摘A mathematical model of principal elements of the aircraft hydraulic system is presented based on the heat transfer theory. The dynamic heat transfer process of the hydraulic oil and the pump shells within an aircraft hydraulic system are analyzed by the difference method. A kind of means for the prediction to variational trends of the aircraft hydraulic system temperature is provided during operation. The numerical prediction and simulation under the operational conditions are presented for ground trial running and the decelerated operation in flight. Computational results show that there is a good coincidence between the experimental data and the numerical predictions.
基金the National Natural Science Foundation of China(52304105)National Natural Science Foundation of China-National major scientific research instrument development project(52227901)Jiangsu Province International Collaboration Program-Key national industrial technology research and development cooperation projects(BZ2023050).
文摘Numerical modelling is an effective technique to improve the understanding of outburst initiation mechanisms and to take appropriate measures to address their threats.Based on the existing two-way sequential coupling method,two typical types of outbursts,i.e.the gas pocket outburst and the dynamic fracturing outburst,have been successfully simulated using field data from a coalfield in central China.The geological structure commonly observed in the coalfield,known as the‘bedding shear zone’,contributes to the gas pocket outbursts in the region.The model for this type of outburst simulates mininginduced stress and gas pressure distributions during the outburst initiation stage and the subsequent development stage.Both coal ejection and gas release are observed in the model,and the simulation results are consistent with mine site observations,i.e.the amount of ejected coal,outburst cavity profile,and gas release rate change prior to an outburst.The second type of outburst is attributed to gas accumulation and elevated gas pressure due to the gassy floor seam and the heterogeneity in the floor strata,which is explained by the dynamic fracturing theory.While the dynamic coal ejection phenomenon is not captured in the simulation,the abrupt release of retained gas from a floor coal seam is successfully replicated.Both outburst models reveal that abnormal gas emission trends can be used as indicators of an upcoming outburst.The results of this study are expected to provide new insights into the outburst initiation mechanisms and outburst prevention measures.
基金National Natural Science Foundation of China(No.52388102)New Cornerstone Science Foundation through the Xplorer Prize.
文摘The dynamic performance of high-speed trains is significantly influenced by sudden changes in aerodynamic loads(ADLs)when exiting a tunnel in a windy environment.Focusing on a double-track tunnel under construction in a mountain railway,we established an aerodynamic model involving a train exiting the tunnel,and verified it in the Fluent environment.Overset mesh technology was adopted to characterize the train’s movement.The flow field involving the train,tunnel,and crosswinds was simulated using the Reynolds-averaged turbulence model.Then,we built a comprehensive train-track coupled dynamic model considering the influences of ADLs,to investigate the vehicles’dynamic responses.The aerodynamics and dynamic behaviors of the train when affected by crosswinds with different velocities and directions are analyzed and discussed.The results show that the near-wall side crosswind leads to sharper variations in ADLs than the far-wall side crosswind.The leading vehicle suffers from more severe ADLs than other vehicles,which worsens the wheel-rail interaction and causes low-frequency vibration of the car body.When the crosswind velocity exceeds 20 m/s,significant wheel-rail impacts occur,and the running safety of the train worsens rapidly.
基金gratefully acknowledge support from the National Key Research and Development Program of China(Grant No.2022YFC3005704)the National Natural Science Foundation of China(Grant No.42277143)the Sichuan Province natural Science Foundation project(Grant No.2024NSFSC0100).
文摘The fragile and intricate geological environment of the Qinghai-Tibet Plateau gives rise to numerous precarious rocks along the riverbanks,posing significant risks for the upcoming construction of hydropower stations.In order to identify potential rockfalls that could endanger the Zixia hydropower project,a comprehensive analysis employing various methods was conducted to investigate the kinematic characteristics and dynamic fragmentation of such precarious rocks.Initially,UAV oblique photography and field survey were used to create a digital elevation model with a resolution of 0.25 m and map the spatial distribution of precarious rocks.Subsequently,the development characteristics of joints within rock masses were analyzed through an adit investigation.Following these preliminary steps,a transportation simulation utilizing RocPro3D,considering stochastic initiation orientation,was employed to predict the trajectories of 18 precarious rocks.As a result,two hazardous rocks that pose a direct threat to the cofferdam were identified.Finally,considering the influence of internal structure planes,a discrete element method was applied for accurately simulating the kinematic characteristics and dynamic fragmentation of these hazardous rocks.The findings underscore several key observations:(1)Slopeparallel structure planes within these hazardous rocks play a pivotal role in both the progressive failure during initiation and dynamic fragmentation during transportation;(2)Hazardous rocksⅢ-1 andⅣ-1 would pose a direct threat to the cofferdam.Notably,block b4 from hazardous rockⅢ-1,could potentially impact the cofferdam with an energy of 4598.65 kJ and an impact force of 3007.5 kN;and(3)Continuous collisions encountered during transportation facilitate the disintegration of rock masses along structure planes and generate substantial high-velocity fragments.Finally,to cope with the impact risk of collapsing blocks,a reinforced retaining wall as the mitigation measure is recommended.
基金the Coordenacao de Aperfeicoamento de Pessoal de Nivel Superior of the Ministry of Education(CAPES/MEC/BRAZIL)-Finance Code 001,Equinor Energy of Brazil Company(Project No.4600025270)the Brazilian National Agency of Petroleum,Gas,and Biofuels(ANP),and Petrobras(Project No.46000579151)+1 种基金INCT/Petroleum Geophysics for financial supportthe National Council for Scientific and Technological Development (CNPq) for their Research Grants of Productivity in Technological Development and Innovation-DT II (313522/2019-7 and 313746/ 2019-2)
文摘This study analyzes Brazilian stromatolites in Lagoa Salgada,serving as analogs for pre-salt rocks in the Santos and Campos basins.Despite their excellent petrophysical properties,such as high porosity and permeability,these reservoirs present challenges in fluid flow modeling and simulation.The research investigates various factors influencing the development of carbonate reservoirs,including diagenetic processes employing several techniques,such as microcomputed tomography(micro-CT)and digital rock physics(DRP),to study petrophysical and geological characteristics.Additionally,through numerical simulations,the properties of fluid flow in different microfacies of stromatolites are estimated,with particular emphasis on understanding and highlighting changes in the direction of fluid flow in the three characterized microfacies.These findings offer crucial insights into optimizing oil and gas exploration and production techniques in carbonate reservoirs,providing a comprehensive understanding of the dynamics of fluid transport in porous media,especially in terms of directional changes within stromatolites.
基金financially supported by the National Natural Science Foundation of China(22478286)。
文摘In this study,four types of spiral fins with varying parameters were mounted on an upstream cylinder,and the effects of spiral fins on the vibration response of heat exchange tubes and the vortex structure in cross flow were studied through experiments and numerical simulations.The results indicate a strong dependency of the cylinder's vibration response on the fin parameters.The results indicate that the vibration response and wake structure of the cylinder are significantly influenced by the parameters of the fins.The introduction of a finned cylinder affects both its own vibration amplitude and frequency,as well as the downstream cylinder.The amplitudes of finned cylinders Ⅰ and Ⅲ are reduced by 57.8% and 59.9%,respectively,compared to the bare cylinder.This reduction helps to restrain vibration and diminishes the amplitudes of the downstream cylinder.Although finned cylinder Ⅱ slightly decreases its own vibration,it increases the amplitude of the downstream cylinder by 13.7%.The mean drag coefficient and the root mean square of the lift coefficient of the finned cylinder are higher than those of the bare cylinder when the finned cylinder is positioned upstream.Smaller pitch and larger equivalent diameter will lead to increased drag,resulting in enhanced vortex shedding in the wake,which amplifies the vibrations of the cylinder in that wake.The downstream of finned cylinder Ⅱ has the widest wake and higher vortex strength,and the dynamic load and vibration of the downstream cylinder are increased.The vortex intensity decays faster in the wake of finned cylinder Ⅲ,and the vibration of the downstream cylinder is weaker.
基金co-supported by the National Natural Science Foundation of China(Nos.52105411,52105400and 52305420)the China Postdoctoral Science Foundation(No.2023M742830)Innovation Foundation for Doctor Dissertation of Northwestern Polytechnical University(No.CX2023008).
文摘In this study,a new linear friction welding(LFW)process,embedded LFW process,was put forward,which was mainly applied to combination manufacturing of long or overlong loadcarrying titanium alloy structural components in aircraft.The interfacial plastic flow behavior and bonding mechanism of this process were investigated by a developed coupling EulerianLagrangian numerical model using software ABAQUS and a novel thermo-physical simulation method with designed embedded hot compression specimen.In addition,the formation mechanism and control method of welding defects caused by uneven plastic flow were discussed.The results reveal that the plastic flow along oscillating direction of this process is even and sufficient.In the direction perpendicular to oscillation,thermo-plastic metals mainly flow downward along welding interface under coupling of shear stress and interfacial pressure,resulting in the interfacial plastic zone shown as an inverted“V”shape.The upward plastic flow in this direction is relatively weak,and only a small amount of flash is extruded from top of joint.Moreover,the wedge block and welding components at top of joint are always in un-steady friction stage,leading to nonuniform temperature field distribution and un-welded defects.According to the results of numerical simulation,high oscillating frequency combined with low pressure and small amplitude is considered as appropriate parameter selection scheme to improve the upward interfacial plastic flow at top of joint and suppress the un-welded defects.The results of thermo-physical simulation illustrate that continuous dynamic recrystallization(CDRX)induces the bonding of interface,accompanying by intense dislocation movement and creation of many low-angle grain boundaries.In the interfacial bonding area,grain orientation is random with relatively low texture density(5.0 mud)owing to CDRX.
文摘Semisubmersible naval ships are versatile military crafts that combine the advantageous features of high-speed planing crafts and submarines.At-surface,these ships are designed to provide sufficient speed and maneuverability.Additionally,they can perform shallow dives,offering low visual and acoustic detectability.Therefore,the hydrodynamic design of a semisubmersible naval ship should address at-surface and submerged conditions.In this study,Numerical analyses were performed using a semisubmersible hull form to analyze its hydrodynamic features,including resistance,powering,and maneuvering.The simulations were conducted with Star CCM+version 2302,a commercial package program that solves URANS equations using the SST k-ωturbulence model.The flow analysis was divided into two parts:at-surface simulations and shallowly submerged simulations.At-surface simulations cover the resistance,powering,trim,and sinkage at transition and planing regimes,with corresponding Froude numbers ranging from 0.42 to 1.69.Shallowly submerged simulations were performed at seven different submergence depths,ranging from D/LOA=0.0635 to D/LOA=0.635,and at two different speeds with Froude numbers of 0.21 and 0.33.The behaviors of the hydrodynamic forces and pitching moment for different operation depths were comprehensively analyzed.The results of the numerical analyses provide valuable insights into the hydrodynamic performance of semisubmersible naval ships,highlighting the critical factors influencing their resistance,powering,and maneuvering capabilities in both at-surface and submerged conditions.
基金funded by the National Natural Science Foundation of China(Grant Nos.42077251,41807269,and U1865203).
文摘Rock masses are often exposed to dynamic loads such as earthquakes and mechanical disturbances in practical engineering scenarios.The existence of underground caverns and weak geological structures like columnar jointed rock masses(CJRMs)and interlayer shear weakness zones(ISWZs)with inferior mechanical properties,significantly undermines the overall structural stability.To tackle the dynamic loading issues in the process of constructing subterranean caverns,a programmable modeling approach was utilized to reconstruct a large-scale underground cavern model incorporating ISWZs and columnar joints(CJs).By conducting dynamic simulations with varying load orientations,the analyses focused on the failure patterns,deformation characteristics,and acoustic emission activity within the caverns.Results revealed that the failure modes of the underground caverns under dynamic loading were predominantly tensile failures.Under X-direction loading,the failed elements were mainly distributed parallel to the CJs,while under Y-direction loading,they were distributed parallel to the transverse weak structural planes.Furthermore,the dynamic stability of the overall structure varied with the number of caverns.The dual-cavern model demonstrated the highest stability under X-direction loading,while the single-cavern model was the least stable.Under Y-direction loading,the cavern stability increased with the number of caverns.Importantly,different weak structures affected the dynamic response of caverns in different ways;the CJRMs were the primary contributors to structural failure,while ISWZs could mitigate the rock mass failure induced by CJs.The findings could offer valuable insights for the dynamic stability analysis of caverns containing CJRMs and ISWZs.
文摘This numerical study of the Bohai Sea wintertime circulation by means of a two dimensional barotropic model with resolution of 1/24° in longitude and latitude showed that the Bohai Sea wintertime circulation is dominated by local monsoon winds. The major current components include the Bohai Warm Current, the North Shandong Coastal Current, and the Liaodong Gyre. The Bohai Warm Current originates from the Yellow Sea Warm Current at the northern part of Bohai Strait, meanders westwards and finally enters the northern part of Bohai Bay. The North Shandong Coastal Current flows along the southwest shore of Bohai Bay and Laizhou Bay and exits from the Bohai Sea through the south Bohai Strait. The anticyclonic Liaodong Gyre is located in the north of Liaodong Bay. A pair of eddies and the small scale Jinzhou Gyre are found between the Bohai Warm Current and the Liaodong Gyre. The computed volume transport for both the Bohai Warm Current and North Shandong Coastal Current is about 0.03 Sv (1 Sv=1×10 6 m 3/s). The numerical experiments showed that the combined effect of local monsoon winds and bottom topography dominate the formation of the circulation pattern. The Coriolis force and the wind stress curl are of certain importance. The beta effect, the momentum advection and the open boundary condition have little influence on the circulation pattern.
基金Supported by the National Natural Science Foundation of China (40874074, 50950110347)the National High Technology Research and Development Program (863 Program) of China (2006AA11ZAA8)Shanghai Science and Technology Development Funds (07ZR14117)
文摘To evaluate the effects of possible ground explosion on a shallow-buried metro tunnel, this paper attempts to analyze the dynamic responses of the operating metro tunnel in soft soil, using a widely applied explicit dynamic nonlinear finite element software ANSYS/LS-DYNA. The blast induced wave propagation in the soil and the tunnel, and the von Mises effective stress and acceleration of the tunnel lining were presented, and the safety of the tunnel lining was evaluated based on the failure criterion. Besides, the parametric study of the soil was also carried out. The numerical results indicate that the upper part of the tunnel lining cross-section with directions ranging from 0° to 22.5° and horizontal distances 0 to 7 m away from the explosive center are the vulnerable areas, and the metro tunnel might be safe when tunnel depth is more than 7 m and TNT charge on the ground is no more than 500 kg, and the selection of soil parameters should be paid more attentions to conduct a more precise analysis.
基金National Natural Science Foundation of China Under Grant No.50678121Open Research Fund Program of State key Laboratory of Hydro-science and Engineering
文摘In this paper, a full-scale 3-D finite element model of the Jundushan cable-stayed aqueduct bridge is established with ANSYS Code. The shell, fluid, tension-only spar and beam elements are used for modeling the aqueduct deck, filled water, cables and support towers, respectively. A multi-element cable formulation is introduced to simulate the cable vibration. The dry （without water） and wet （with water） modes of the aqueduct bridge are both extracted and investigated in detail. The dry modes of the aqueduct bridge are basically similar to those of highway cable-stayed bridges. A dry mode may correspond to two types of wet modes, which are called the in-phase （with lower frequency） and out-of-phase （with higher frequency） modes. When the water-structure system vibrates in the in-phase/out-of-phase modes, the aqueduct deck moves and water sloshes in the same/opposite phase-angle, and the sloshing water may take different surface-wave modes. The wet modes of the system reflect the properties of interaction among the deck, towers, cables and water. The in-phase wet frequency generally decreases as the water depth increases, and the out-of-phase wet frequency may increase or decrease as the water depth increases.
