High rock temperature is a great challenge frequently encountered during subsurface resource recovery and deep underground space utilization,and it is still unclear how the granitic rock responds to realtime high temp...High rock temperature is a great challenge frequently encountered during subsurface resource recovery and deep underground space utilization,and it is still unclear how the granitic rock responds to realtime high temperature upon shear loading.To better understand the shear fracture behavior and underlying processes of intact granite exposed to thermal-mechanical coupling loading,direct shear tests were conducted utilizing a newly built testing apparatus at varied normal stresses and high temperatures.Influencesof different temperatures and different heating methods(real-time heating and thermal treatment)on the shear mechanical behavior were compared and discussed.Results indicate that shear stress fluctuationswith some small stress drops occur as shear stress is approaching the peak strength under real-time heating,accompanied by more and earlier AE signal uprushes.This suggests that greater cracking events occur earlier during real-time heating than after thermal treatment,resulting in a lower peak shear strength.Furthermore,the peak shear strength,post-peak stress drop,and cohesion rise from room temperature(RT)to 200℃(the peak strength increases by 8%,5.8%,and 9.9%under normal stress of 5 MPa,15 MPa,and 20 MPa,correspondingly),and subsequently decline from 200℃to 400℃.Temperature has a limited impact on shear stiffness from RT to 200℃,but significantlyreduces it from 200℃to 400℃,with drops of 15%,7.9%,and 10%under normal stress of 5 MPa,15 MPa,and 20 MPa,respectively.Moreover,the shear strength and stiffness under real-time heating are lower than those for the thermally treated specimens.The strengthening of intact granite below 200℃upon shear is associated with loss of water and a more compacted structure,while the weakening effect of temperature on shear strength from 200℃to 400℃is due to the new thermal cracks and less brittle and stiff of minerals.展开更多
High geo-stress and high temperature in deep rock engineering increase the possibility of engineering and geological disasters in discontinuous rocks.However,the influence of thermomechanical coupling on the shear beh...High geo-stress and high temperature in deep rock engineering increase the possibility of engineering and geological disasters in discontinuous rocks.However,the influence of thermomechanical coupling on the shear behavior and damage evolution of prefractured granite remains immature.In this context,true triaxial laboratory tests and discrete element method simulations under different confining pressures(σ3=3 MPa,σ2=4 MPa,andσ3=80 MPa,σ2=100 MPa)and temperatures(25℃-500℃)were carried out on rough granite fractures with two different orientations.Results indicate that high temperature and high confining pressure increase the peak strength of the prefractured specimen,leading to more microcracks in the host rock and more gouges between the surfaces.Thermal strengthening at low temperatures(<300℃)and residual stick-slip only occur under a greater confining pressure for prefractured specimens.High confining pressure suppresses generation of the thermal microcracks in the heating stage.Cracks first initiate in the asperities on the fracture surfaces,and then propagate into the rock matrix during the mechanical loading stage.In addition,prefractured granite with a larger fracture angle is characterized by smaller peak and residual strength,faster residual slip,fewer new cracks on the specimen surface,and a more pronounced thermal strengthening effect on peak strength.The slip tendency analysis indicates that a higher maximum principal stress(s1)and a large fracture angle(45°-75°)generally result in a higher potential for fracture slip or activation.This study will contribute to a better understanding of the fracture shear mechanism under true triaxial thermomechanical coupling conditions and provides new insights into the stability evaluation of deep dynamic geological hazards.展开更多
The dynamic characteristics of a single liquid-filled pipe have been broadly studied in the previous literature.The parallel liquid-filled pipe(PLFP)system is also widely used in engineering,and its structure is more ...The dynamic characteristics of a single liquid-filled pipe have been broadly studied in the previous literature.The parallel liquid-filled pipe(PLFP)system is also widely used in engineering,and its structure is more complex than that of a single pipe.However,there are few reports about the dynamic characteristics of the PLFPs.Therefore,this paper proposes improved frequency modeling and solution for the PLFPs,involving the logical alignment principle and coupled matrix processing.The established model incorporates both the fluid-structure interaction(FSI)and the structural coupling of the PLFPs.The validity of the established model is verified by modal experiments.The effects of some unique parameters on the dynamic characteristics of the PLFPs are discussed.This work provides a feasible method for solving the FSI of multiple pipes in parallel and potential theoretical guidance for the dynamic analysis of the PLFPs in engineering.展开更多
The vibration of a Francis turbine is analyzed with the additional quality matrix method based on fluid-structure coupling （FSC）. Firstly, the vibration frequency and mode of blade and runner in air and water are ca...The vibration of a Francis turbine is analyzed with the additional quality matrix method based on fluid-structure coupling （FSC）. Firstly, the vibration frequency and mode of blade and runner in air and water are calculated. Secondly, the influences to runner frequency domain by large flow, small flow and design flow working conditions are compared. Finally the influences to runner modes by centrifugal forces under three rotating speeds of 400 r/rain, 500 r/min and 600 r/rain are compared. The centrifugal force and small flow working condition have greatly influence on the vibration of small runner. With the increase of centrifugal force, the vibration frequency of the runner is sharply increased. Some order frequencies are even close to the runner natural frequency in the air. Because the low frequency vibration will severely damage the stability of the turbine, low frequency vibration of units should be avoided as soon as possible.展开更多
Based on global initiatives such as the clean energy transition and the development of renewable energy,the pumped storage power station has become a new and significant way of energy storage and regulation,and its co...Based on global initiatives such as the clean energy transition and the development of renewable energy,the pumped storage power station has become a new and significant way of energy storage and regulation,and its construction environment is more complex than that of a traditional reservoir.In particular,the stability of the rock strata in the underground reservoirs is affected by the seepage pressure and rock stress,which presents some challenges in achieving engineering safety and stability.Using the advantages of the numerical simulation method in dealing deal with nonlinear problems in engineering stability,in this study,the stability of the underground reservoir of the Shidangshan(SDS)pumped storage power station was numerically calculated and quantitatively analyzed based on fluid-structure coupling theory,providing an important reference for the safe operation and management of the underground reservoir.First,using the COMSOL software,a suitablemechanicalmodel was created in accordance with the geological structure and project characteristics of the underground reservoir.Next,the characteristics of the stress field,displacement field,and seepage field after excavation of the underground reservoir were simulated in light of the seepage effect of groundwater on the nearby rock of the underground reservoir.Finally,based on the construction specifications and Molar-Coulomb criterion,a thorough evaluation of the stability of the underground reservoir was performed through simulation of the filling and discharge conditions and anti-seepage strengthening measures.The findings demonstrate that the numerical simulation results have a certain level of reliability and are in accordance with the stress measured in the project area.The underground reservoir excavation resulted in a maximum displacement value of the rock mass around the caverns of 3.56 mm in a typical section,and the safety coefficient of the parts,as determined using the Molar-Coulomb criterion,was higher than 1,indicating that the project as a whole is in a stable state.展开更多
Based on the theory of Housner, the transverse seismic response of beam aqueduct considering fluid-structure coupling is established. With the variation of aqueduct cross-section ratio of depth to width, the aqueduct ...Based on the theory of Housner, the transverse seismic response of beam aqueduct considering fluid-structure coupling is established. With the variation of aqueduct cross-section ratio of depth to width, the aqueduct transverse seismic response change. The transverse seismic response of a large-scale aqueduct in several work condition are calculated. It shows that the transverse seismic response is greatly influenced by the water mass in the aqueduct, but the shaking water play a TLD role. ff the whole water is appended aqueduct body, it will magnify seismic inertia action. When aqueduct cross-section is selected, the influence of ratio of depth and width to pier seismic response should be considered in order to reduce seismic action.展开更多
The development and rapid usage of numerical codes for fluid-structure interaction(FSI) problems are of great relevance to researchers in many engineering fields such as civil engineering and ocean engineering. This m...The development and rapid usage of numerical codes for fluid-structure interaction(FSI) problems are of great relevance to researchers in many engineering fields such as civil engineering and ocean engineering. This multidisciplinary field known as FSI has been expanded to engineering fields such as offshore structures, tall slender structures and other flexible structures applications. The motivation of this paper is to investigate the numerical model of two-way coupling FSI partitioned flexible plate structure under fluid flow. The adopted partitioned method and approach utilized the advantage of the existing numerical algorithms in solving the two-way coupling fluid and structural interactions. The flexible plate was subjected to a fluid flow which causes large deformation on the fluid domain from the oscillation of the flexible plate. Both fluid and flexible plate are subjected to the interaction of load transfer within two physics by using the strong and weak coupling methods of MFS and Load Transfer Physics Environment, respectively. The oscillation and deformation results have been validated which demonstrate the reliability of both strong and weak method in resolving the two-way coupling problem in contribution of knowledge to the feasibility field study of ocean engineering and civil engineering.展开更多
The multi-physics simulation of coupled fluid-structure interaction problems, with disjoint fluid and solid domains, requires one to choose a method for enforcing the fluid-structure coupling at the interface between ...The multi-physics simulation of coupled fluid-structure interaction problems, with disjoint fluid and solid domains, requires one to choose a method for enforcing the fluid-structure coupling at the interface between solid and fluid. While it is common knowledge that the choice of coupling technique can be very problem dependent, there exists no satisfactory coupling comparison methodology that allows for conclusions to be drawn with respect to the comparison of computational cost and solution accuracy for a given scenario. In this work, we develop a computational framework where all aspects of the computation can be held constant, save for the method in which the coupled nature of the fluid-structure equations is enforced. To enable a fair comparison of coupling methods, all simulations presented in this work are implemented within a single numerical framework within the deal.ii [1] finite element library. We have chosen the two-dimensional benchmark test problem of Turek and Hron [2] as an example to examine the relative accuracy of the coupling methods studied;however, the comparison technique is equally applicable to more complex problems. We show that for the specific case considered herein the monolithic approach outperforms partitioned and quasi-direct methods;however, this result is problem dependent and we discuss computational and modeling aspects which may affect other comparison studies.展开更多
The upper crustal anisotropy of Yunnan area, SE margin of Tibetan Plateau, is investigated by measuring the shear wave splitting of local earthquakes. The mean value of the measured delay times is 0.054 s and far less...The upper crustal anisotropy of Yunnan area, SE margin of Tibetan Plateau, is investigated by measuring the shear wave splitting of local earthquakes. The mean value of the measured delay times is 0.054 s and far less than that from Pms splitting analysis, indicating that the crustal anisotropy is contributed mostly from mid-lower crust. The fast polarization directions are mostly sub-parallel to the maximum horizontal compression directions while the stations near fault zones show fault-parallel fast polarization directions, suggesting both stress and geological structure contribute to the upper crust anisotropy.Comparing fast polarization directions from shear wave splitting of local earthquakes and Pms, large angle differences are shown at most stations, implying different anisotropy properties between upper and mid-lower crust. However, in southwestern Yunnan, the fast polarization directions of Pms and Swave splitting are nearly parallel, and the stress and surface strain rate directions show strong correlation, which may indicate that the surface and deep crust deformations can be explained by the same mechanism and the surface deformation can represent the deformation of the whole crust. Therefore,the high correlation between surface strain and mantle deformation in this area suggests the mechanical coupling between crust and mantle in southwestern Yunnan. In the rest region of Yunnan, the crustmantle coupling mechanisms are supported by the lack of significant crustal anisotropy with Ne S fast polarization directions from Pms splitting. Therefore, we conclude that the crust and upper mantle are coupled in Yunnan, SE margin of Tibetan Plateau.展开更多
Cemented paste backfill(CPB)is a key technology for green mining in metal mines,in which tailings thickening comprises the primary link of CPB technology.However,difficult flocculation and substandard concentrations o...Cemented paste backfill(CPB)is a key technology for green mining in metal mines,in which tailings thickening comprises the primary link of CPB technology.However,difficult flocculation and substandard concentrations of thickened tailings often occur.The rheological properties and concentration evolution in the thickened tailings remain unclear.Moreover,traditional indoor thickening experiments have yet to quantitatively characterize their rheological properties.An experiment of flocculation condition optimization based on the Box-Behnken design(BBD)was performed in the study,and the two response values were investigated:concentration and the mean weighted chord length(MWCL)of flocs.Thus,optimal flocculation conditions were obtained.In addition,the rheological properties and concentration evolution of different flocculant dosages and ultrafine tailing contents under shear,compression,and compression-shear coupling experimental conditions were tested and compared.The results show that the shear yield stress under compression and compression-shear coupling increases with the growth of compressive yield stress,while the shear yield stress increases slightly under shear.The order of shear yield stress from low to high under different thickening conditions is shear,compression,and compression-shear coupling.Under compression and compression-shear coupling,the concentration first rapidly increases with the growth of compressive yield stress and then slowly increases,while concentration increases slightly under shear.The order of concentration from low to high under different thickening conditions is shear,compression,and compression-shear coupling.Finally,the evolution mechanism of the flocs and drainage channels during the thickening of the thickened tailings under different experimental conditions was revealed.展开更多
The shear failure of intact rock under thermo-mechanical(TM)coupling conditions is common,such as in enhanced geothermal mining and deep mine construction.Under the effect of a continuous engineering disturbance,shear...The shear failure of intact rock under thermo-mechanical(TM)coupling conditions is common,such as in enhanced geothermal mining and deep mine construction.Under the effect of a continuous engineering disturbance,shear-formed fractures are prone to secondary instability,posing a severe threat to deep engineering.Although numerous studies regarding three-dimensional(3D)morphologies of fracture surfaces have been conducted,the understanding of shear-formed fractures under TM coupling conditions is limited.In this study,direct shear tests of intact granite under various TM coupling conditions were conducted,followed by 3D laser scanning tests of shear-formed fractures.Test results demonstrated that the peak shear strength of intact granite is positively correlated with the normal stress,whereas it is negatively correlated with the temperature.The internal friction angle and cohesion of intact granite significantly decrease with an increase in the temperature.The anisotropy,roughness value,and height of the asperities on the fracture surfaces are reduced as the normal stress increases,whereas their variation trends are the opposite as the temperature increases.The macroscopic failure mode of intact granite under TM coupling conditions is dominated by mixed tensileeshear and shear failures.As the normal stress increases,intragranular fractures are developed ranging from a local to a global distribution,and the macroscopic failure mode of intact granite changes from mixed tensileeshear to shear failure.Finally,3D morphological characteristics of the asperities on the shear-formed fracture surfaces were analyzed,and a quadrangular pyramid conceptual model representing these asperities was proposed and sufficiently verified.展开更多
Underwater cylindrical shell structures have been found a wide of application in many engineering fields, such as the element of marine, oil platforms, etc. The coupled vibration analysis is a hot issue for these unde...Underwater cylindrical shell structures have been found a wide of application in many engineering fields, such as the element of marine, oil platforms, etc. The coupled vibration analysis is a hot issue for these underwater structures. The vibration characteristics of underwater structures are influenced not only by hydrodynamic pressure but also by hydrostatic pressure corresponding to different water depths. In this study, an acoustic finite element method was used to evaluate the underwater structures. Taken the hydrostatic pressure into account in terms of initial stress stiffness, an acoustical fluid-structure coupled analysis of underwater cylindrical shells has been made to study the effect of hydrodynamic pressures on natural frequency and sound radiation. By comparing with the frequencies obtained by the acoustic finite element method and by the added mass method based on the Bessel function, the validity of present analysis was checked. Finally, test samples of the sound radiation of stiffened cylindrical shells were acquired by a harmonic acoustic analysis. The results showed that hydrostatic pressure plays an important role in determining a large submerged body motion, and the characteristics of sound radiation change with water depth. Furthermore, the analysis methods and the results are of significant reference value for studies of other complicated submarine structures.展开更多
We propose a suite of strategies for the parallel solution of fully implicit monolithic fluid-structure interaction(FSI).The solver is based on a modeling approach that uses the velocity and pressure as the primitive ...We propose a suite of strategies for the parallel solution of fully implicit monolithic fluid-structure interaction(FSI).The solver is based on a modeling approach that uses the velocity and pressure as the primitive variables,which offers a bridge between computational fluid dynamics(CFD)and computational structural dynamics.The spatiotemporal discretization leverages the variational multiscale formulation and the generalized-αmethod as a means of providing a robust discrete scheme.In particular,the time integration scheme does not suffer from the overshoot phenomenon and optimally dissipates high-frequency spurious modes in both subproblems of FSI.Based on the chosen fully implicit scheme,we systematically develop a combined suite of nonlinear and linear solver strategies.Invoking a block factorization of the Jacobian matrix,the Newton-Raphson procedure is reduced to solving two smaller linear systems in the multi-corrector stage.The first is of the elliptic type,indicating that the algebraic multigrid method serves as a well-suited option.The second exhibits a two-by-two block structure that is analogous to the system arising in CFD.Inspired by prior studies,the additive Schwarz domain decomposition method and the block-factorization-based preconditioners are invoked to address the linear problem.Since the number of unknowns matches in both subdomains,it is straightforward to balance loads when parallelizing the algorithm for distributed-memory architectures.We use two representative FSI benchmarks to demonstrate the robustness,efficiency,and scalability of the overall FSI solver framework.In particular,it is found that the developed FSI solver is comparable to the CFD solver in several aspects,including fixed-size and isogranular scalability as well as robustness.展开更多
Rotating machinery in the aviation industry is increasingly embracing high speeds and miniaturization,and foil dynamic pressure gas bearing has great application value due to its self-lubrication and self-adaptive def...Rotating machinery in the aviation industry is increasingly embracing high speeds and miniaturization,and foil dynamic pressure gas bearing has great application value due to its self-lubrication and self-adaptive deformation characteristics.This study explores the interaction mechanism between micro-scale variable-sectional shearing flow with hyper-rotation speeds and a three-layer elastic foil assembly through bidirectional aero-elastic coupling in a Multi-layer Thrust Gas Foil Bearing(MTGFB).The bearing capacity of the MTGFB varies non-linearly with the decrease of gas film clearance,while the collaborative deformation of the three-layer elastic foil assembly can deal with different load conditions.As the load capacity increases,the enhanced dynamic pressure effect causes the top foil to evolve from a single arch to multiple arches.The hydrodynamic effects in the gas film evolve to form multiple segmented wedges with different pitch ratios,while the peak pressure of the gas film always occurs near the vaults of the top foil.As the rotational speed frequency approaches the natural frequency,the resonance of the gas film and elastic foil assembly system occurs,and a phase delay occurs between the pressure pulsation and the vibration of foils.The load capacity of the MTGFB also depends on the elastic moduli of the elastic foil assembly.Increasing the elastic modulus decreases the deformation amplitude of the top foil,whereas it increases those of the backboard and middle foil,increasing the load capacity.展开更多
Replaceable flexural and shear fuse-type coupling beams are used in hybrid coupled shear wall(HCSW)systems,enabling concrete buildings to be promptly recovered after severe earthquakes.This study aimed to analytically...Replaceable flexural and shear fuse-type coupling beams are used in hybrid coupled shear wall(HCSW)systems,enabling concrete buildings to be promptly recovered after severe earthquakes.This study aimed to analytically evaluate the seismic behavior of flexural and shear fuse beams situated in short-,medium-and high-rise RC buildings that have HCSWs.Three building groups hypothetically located in a high seismic hazard zone were studied.A series of 2D nonlinear time history analyses was accomplished in OpenSees,using the ground motion records scaled at the design basis earthquake level.It was found that the effectiveness of fuses in HCSWs depends on various factors such as size and scale of the building,allowable rotation value,inter-story drift ratio,residual drift quantity,energy dissipation value of the fuses,etc.The results show that shear fuses better meet the requirements of rotations and drifts.In contrast,flexural fuses dissipate more energy,but their sectional stiffness should increase to meet other requirements.It was concluded that adoption of proper fuses depends on the overall scale of the building and on how associated factors are considered.展开更多
The development of thermal stress in the exhaust manifold of a gasoline engine is considered.The problem is addresses in the frame of a combined approach wherefluid and structure are coupled using the GT-POWER and STA...The development of thermal stress in the exhaust manifold of a gasoline engine is considered.The problem is addresses in the frame of a combined approach wherefluid and structure are coupled using the GT-POWER and STAR-CCM+software.First,the external characteristic curve of the engine is compared with a one-dimen-sional simulation model,then the parameters of the model are modified until the curve matches the available experimental values.GT-POWER is then used to transfer the inlet boundary data under transient conditions to STAR-CCM+in real-time.The temperature profiles of the inner and outer walls of the exhaust manifold are obtained in this way,together with the thermal stress and thermal deformation of the exhaust manifold itself.Using this information,the original model is improved through the addition of connections.Moreover,the local branch pipes are optimized,leading to significant improvements in terms of thermal stress and thermal deforma-tion of the exhaust manifold(a 7%reduction in the maximum thermal stress).展开更多
This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the...This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the rigid structure is taken as ＂fictitious＂ fluid with zero strain rate. Both fluid and structure are described by velocity and pressure. The whole domain, including fluid region and structure region, is modeled by the incompressible Navier-Stokes equations which are discretized with fixed Eulerian mesh. However, to keep the structure＇ s rigid body shape and behavior, a rigid body constraint is enforced on the ＂fictitious＂ fluid domain by use of the Distributed Lagrange Multipher/Fictitious Domain （DLM/ FD） method which is originally introduced to solve particulate flow problems by Glowinski et al. For the verification of the model presented herein, a 2D numerical wave tank is established to simulate small amplitude wave propagations, and then numerical results are compared with analytical solutions. Finally, a 2D example of fluid-structure interaction under wave dynamic forces provides convincing evidences for the method excellent solution quality and fidelity.展开更多
Fluid-structure interaction （FSI） is essentially a dynamic phenomenon and always exists in fluid-filled pipe system. The four-equation model, which has been proved to be effective to describe and predict the phenome...Fluid-structure interaction （FSI） is essentially a dynamic phenomenon and always exists in fluid-filled pipe system. The four-equation model, which has been proved to be effective to describe and predict the phenomenon of FSI due to friction coupling and Poisson coupling being taken into account, is utilized to describe the FSI of fluid-filled pipe system. Terse compatibility equations are educed by the method of characteristics （MOC） to describe the fluid-filled pipe system. To shorten computing time needed to get the solutions under the condition of keeping accuracy requirement, two steps are adopted, firstly the time step Δt and divided number of the straight pipe are optimized, sec-ondly the mesh spacing Δz close to boundary is subdivided in several submeshes automatically ac-cording to the speed gradient of fluid. The mathematical model and arithmetic are validated by com-parisons between simulation solutions of two straight pipe systems and experiment known from lit-erature.展开更多
A novel algorithm is proposed for the simulation of fluid-structure interaction problems.In particular,much attention is paid to natural phenomena such as debris flow.The fluid part(debris flow fluid)is simulated in t...A novel algorithm is proposed for the simulation of fluid-structure interaction problems.In particular,much attention is paid to natural phenomena such as debris flow.The fluid part(debris flow fluid)is simulated in the framework of the smoothed particle hydrodynamics(SPH)approach,while the solid part(downstream obstacles)is treated using the finite element method(FEM).Fluid-structure coupling is implemented through dynamic boundary conditions.In particular,the software“TensorFlow”and an algorithm based on Python are combined to conduct the required calculations.The simulation results show that the dynamics of viscous and non-viscous debris flows can be extremely different when there are obstacles in the downstream direction.The implemented SPH-FEM coupling method can simulate the fluid-structure coupling problem with a reasonable approximation.展开更多
To analyze wind-induced response characteristics of a wind turbine tower more accurately, the blade-tower coupling effect was investigated. The mean wind velocity of the rotating blades and tower was simulated accordi...To analyze wind-induced response characteristics of a wind turbine tower more accurately, the blade-tower coupling effect was investigated. The mean wind velocity of the rotating blades and tower was simulated according to wind shear effects, and the fluctuating wind velocity time series of the wind turbine were simulated by a harmony superposition method. A dynamic finite element method (FEM) was used to calculate the wind-induced response of the blades and tower. Wind-induced responses of the tower were calculated in two cases (one included the blade-tower coupling effect, and the other only added the mass of blades and the hub at the top of the tower), and then the maximal displacements at the top of the tower of the tow cases were compared with each other. As a result of the influence of the blade-tower coupling effect and the total base shear of the blades, the maximal displacement of the first case increased nearly by 300% compared to the second case. To obtain more precise analysis, the blade-tower coupling effect and the total base shear of the blades should be considered simultaneously in the design of wind turbine towers.展开更多
基金support from the Taishan Scholars Program,Key Research Program of Frontier Sciences,Chinese Academy of Sciences(CAS),Grant No.ZDBS-LY-DQC022Open Research Fund of State Key Laboratory of Geomechanics and Geotechnical Engineering Safety,Grant No.SKLGME023003.
文摘High rock temperature is a great challenge frequently encountered during subsurface resource recovery and deep underground space utilization,and it is still unclear how the granitic rock responds to realtime high temperature upon shear loading.To better understand the shear fracture behavior and underlying processes of intact granite exposed to thermal-mechanical coupling loading,direct shear tests were conducted utilizing a newly built testing apparatus at varied normal stresses and high temperatures.Influencesof different temperatures and different heating methods(real-time heating and thermal treatment)on the shear mechanical behavior were compared and discussed.Results indicate that shear stress fluctuationswith some small stress drops occur as shear stress is approaching the peak strength under real-time heating,accompanied by more and earlier AE signal uprushes.This suggests that greater cracking events occur earlier during real-time heating than after thermal treatment,resulting in a lower peak shear strength.Furthermore,the peak shear strength,post-peak stress drop,and cohesion rise from room temperature(RT)to 200℃(the peak strength increases by 8%,5.8%,and 9.9%under normal stress of 5 MPa,15 MPa,and 20 MPa,correspondingly),and subsequently decline from 200℃to 400℃.Temperature has a limited impact on shear stiffness from RT to 200℃,but significantlyreduces it from 200℃to 400℃,with drops of 15%,7.9%,and 10%under normal stress of 5 MPa,15 MPa,and 20 MPa,respectively.Moreover,the shear strength and stiffness under real-time heating are lower than those for the thermally treated specimens.The strengthening of intact granite below 200℃upon shear is associated with loss of water and a more compacted structure,while the weakening effect of temperature on shear strength from 200℃to 400℃is due to the new thermal cracks and less brittle and stiff of minerals.
基金support from the National Key Research and Development Program of China(Grant No.2022YFE0137200)supported by the Taishan Scholars Program and Open Research Fund of State Key Laboratory of Geomechanics and Geotechnical Engineering Safety(Grant No.SKLGME023003).
文摘High geo-stress and high temperature in deep rock engineering increase the possibility of engineering and geological disasters in discontinuous rocks.However,the influence of thermomechanical coupling on the shear behavior and damage evolution of prefractured granite remains immature.In this context,true triaxial laboratory tests and discrete element method simulations under different confining pressures(σ3=3 MPa,σ2=4 MPa,andσ3=80 MPa,σ2=100 MPa)and temperatures(25℃-500℃)were carried out on rough granite fractures with two different orientations.Results indicate that high temperature and high confining pressure increase the peak strength of the prefractured specimen,leading to more microcracks in the host rock and more gouges between the surfaces.Thermal strengthening at low temperatures(<300℃)and residual stick-slip only occur under a greater confining pressure for prefractured specimens.High confining pressure suppresses generation of the thermal microcracks in the heating stage.Cracks first initiate in the asperities on the fracture surfaces,and then propagate into the rock matrix during the mechanical loading stage.In addition,prefractured granite with a larger fracture angle is characterized by smaller peak and residual strength,faster residual slip,fewer new cracks on the specimen surface,and a more pronounced thermal strengthening effect on peak strength.The slip tendency analysis indicates that a higher maximum principal stress(s1)and a large fracture angle(45°-75°)generally result in a higher potential for fracture slip or activation.This study will contribute to a better understanding of the fracture shear mechanism under true triaxial thermomechanical coupling conditions and provides new insights into the stability evaluation of deep dynamic geological hazards.
基金Project supported by the National Natural Science Foundation of China(No.11972112)the Fundamental Research Funds for the Central Universities of China(Nos.N2103024 and N2103002)the Major Projects of Aero-Engines and Gasturbines(No.J2019-I-0008-0008)。
文摘The dynamic characteristics of a single liquid-filled pipe have been broadly studied in the previous literature.The parallel liquid-filled pipe(PLFP)system is also widely used in engineering,and its structure is more complex than that of a single pipe.However,there are few reports about the dynamic characteristics of the PLFPs.Therefore,this paper proposes improved frequency modeling and solution for the PLFPs,involving the logical alignment principle and coupled matrix processing.The established model incorporates both the fluid-structure interaction(FSI)and the structural coupling of the PLFPs.The validity of the established model is verified by modal experiments.The effects of some unique parameters on the dynamic characteristics of the PLFPs are discussed.This work provides a feasible method for solving the FSI of multiple pipes in parallel and potential theoretical guidance for the dynamic analysis of the PLFPs in engineering.
基金Outstanding Youth Science Fund Subsidization of Sichuan Province, China (No. 05204033).
文摘The vibration of a Francis turbine is analyzed with the additional quality matrix method based on fluid-structure coupling （FSC）. Firstly, the vibration frequency and mode of blade and runner in air and water are calculated. Secondly, the influences to runner frequency domain by large flow, small flow and design flow working conditions are compared. Finally the influences to runner modes by centrifugal forces under three rotating speeds of 400 r/rain, 500 r/min and 600 r/rain are compared. The centrifugal force and small flow working condition have greatly influence on the vibration of small runner. With the increase of centrifugal force, the vibration frequency of the runner is sharply increased. Some order frequencies are even close to the runner natural frequency in the air. Because the low frequency vibration will severely damage the stability of the turbine, low frequency vibration of units should be avoided as soon as possible.
基金funded by the BeijingNatural Science Foundation of China(8222003)National Natural Science Foundation of China(41807180).
文摘Based on global initiatives such as the clean energy transition and the development of renewable energy,the pumped storage power station has become a new and significant way of energy storage and regulation,and its construction environment is more complex than that of a traditional reservoir.In particular,the stability of the rock strata in the underground reservoirs is affected by the seepage pressure and rock stress,which presents some challenges in achieving engineering safety and stability.Using the advantages of the numerical simulation method in dealing deal with nonlinear problems in engineering stability,in this study,the stability of the underground reservoir of the Shidangshan(SDS)pumped storage power station was numerically calculated and quantitatively analyzed based on fluid-structure coupling theory,providing an important reference for the safe operation and management of the underground reservoir.First,using the COMSOL software,a suitablemechanicalmodel was created in accordance with the geological structure and project characteristics of the underground reservoir.Next,the characteristics of the stress field,displacement field,and seepage field after excavation of the underground reservoir were simulated in light of the seepage effect of groundwater on the nearby rock of the underground reservoir.Finally,based on the construction specifications and Molar-Coulomb criterion,a thorough evaluation of the stability of the underground reservoir was performed through simulation of the filling and discharge conditions and anti-seepage strengthening measures.The findings demonstrate that the numerical simulation results have a certain level of reliability and are in accordance with the stress measured in the project area.The underground reservoir excavation resulted in a maximum displacement value of the rock mass around the caverns of 3.56 mm in a typical section,and the safety coefficient of the parts,as determined using the Molar-Coulomb criterion,was higher than 1,indicating that the project as a whole is in a stable state.
基金National Natural Science Foundation of China (50279024)Lanzhou Jiaotong University Qinglan Talented Person Fund Project (QL-05-13A).
文摘Based on the theory of Housner, the transverse seismic response of beam aqueduct considering fluid-structure coupling is established. With the variation of aqueduct cross-section ratio of depth to width, the aqueduct transverse seismic response change. The transverse seismic response of a large-scale aqueduct in several work condition are calculated. It shows that the transverse seismic response is greatly influenced by the water mass in the aqueduct, but the shaking water play a TLD role. ff the whole water is appended aqueduct body, it will magnify seismic inertia action. When aqueduct cross-section is selected, the influence of ratio of depth and width to pier seismic response should be considered in order to reduce seismic action.
文摘The development and rapid usage of numerical codes for fluid-structure interaction(FSI) problems are of great relevance to researchers in many engineering fields such as civil engineering and ocean engineering. This multidisciplinary field known as FSI has been expanded to engineering fields such as offshore structures, tall slender structures and other flexible structures applications. The motivation of this paper is to investigate the numerical model of two-way coupling FSI partitioned flexible plate structure under fluid flow. The adopted partitioned method and approach utilized the advantage of the existing numerical algorithms in solving the two-way coupling fluid and structural interactions. The flexible plate was subjected to a fluid flow which causes large deformation on the fluid domain from the oscillation of the flexible plate. Both fluid and flexible plate are subjected to the interaction of load transfer within two physics by using the strong and weak coupling methods of MFS and Load Transfer Physics Environment, respectively. The oscillation and deformation results have been validated which demonstrate the reliability of both strong and weak method in resolving the two-way coupling problem in contribution of knowledge to the feasibility field study of ocean engineering and civil engineering.
文摘The multi-physics simulation of coupled fluid-structure interaction problems, with disjoint fluid and solid domains, requires one to choose a method for enforcing the fluid-structure coupling at the interface between solid and fluid. While it is common knowledge that the choice of coupling technique can be very problem dependent, there exists no satisfactory coupling comparison methodology that allows for conclusions to be drawn with respect to the comparison of computational cost and solution accuracy for a given scenario. In this work, we develop a computational framework where all aspects of the computation can be held constant, save for the method in which the coupled nature of the fluid-structure equations is enforced. To enable a fair comparison of coupling methods, all simulations presented in this work are implemented within a single numerical framework within the deal.ii [1] finite element library. We have chosen the two-dimensional benchmark test problem of Turek and Hron [2] as an example to examine the relative accuracy of the coupling methods studied;however, the comparison technique is equally applicable to more complex problems. We show that for the specific case considered herein the monolithic approach outperforms partitioned and quasi-direct methods;however, this result is problem dependent and we discuss computational and modeling aspects which may affect other comparison studies.
基金supported by the National 973 Project of China (No.2013CB733303)the open fund of Key Laboratory of Geospace Environment and Geodesy,Ministry of Education (No.15-02-07)
文摘The upper crustal anisotropy of Yunnan area, SE margin of Tibetan Plateau, is investigated by measuring the shear wave splitting of local earthquakes. The mean value of the measured delay times is 0.054 s and far less than that from Pms splitting analysis, indicating that the crustal anisotropy is contributed mostly from mid-lower crust. The fast polarization directions are mostly sub-parallel to the maximum horizontal compression directions while the stations near fault zones show fault-parallel fast polarization directions, suggesting both stress and geological structure contribute to the upper crust anisotropy.Comparing fast polarization directions from shear wave splitting of local earthquakes and Pms, large angle differences are shown at most stations, implying different anisotropy properties between upper and mid-lower crust. However, in southwestern Yunnan, the fast polarization directions of Pms and Swave splitting are nearly parallel, and the stress and surface strain rate directions show strong correlation, which may indicate that the surface and deep crust deformations can be explained by the same mechanism and the surface deformation can represent the deformation of the whole crust. Therefore,the high correlation between surface strain and mantle deformation in this area suggests the mechanical coupling between crust and mantle in southwestern Yunnan. In the rest region of Yunnan, the crustmantle coupling mechanisms are supported by the lack of significant crustal anisotropy with Ne S fast polarization directions from Pms splitting. Therefore, we conclude that the crust and upper mantle are coupled in Yunnan, SE margin of Tibetan Plateau.
基金financially supported by the National Natural Science Foundation of China(Nos.52130404 and 52304121)the Fundamental Research Funds for the Central Universities(No.FRF-TP-22-112A1)+4 种基金the Guangdong Basic and Applied Basic Research Foundation(No.2021A 1515110161)the ANID(Chile)through Fondecyt project 1210610the Centro de Modelamiento Matemático(BASAL funds for Centers of Excellence FB210005)the CRHIAM project ANID/FONDAP/15130015 and ANID/FONDAP/1523A0001the Anillo project ANID/ACT210030。
文摘Cemented paste backfill(CPB)is a key technology for green mining in metal mines,in which tailings thickening comprises the primary link of CPB technology.However,difficult flocculation and substandard concentrations of thickened tailings often occur.The rheological properties and concentration evolution in the thickened tailings remain unclear.Moreover,traditional indoor thickening experiments have yet to quantitatively characterize their rheological properties.An experiment of flocculation condition optimization based on the Box-Behnken design(BBD)was performed in the study,and the two response values were investigated:concentration and the mean weighted chord length(MWCL)of flocs.Thus,optimal flocculation conditions were obtained.In addition,the rheological properties and concentration evolution of different flocculant dosages and ultrafine tailing contents under shear,compression,and compression-shear coupling experimental conditions were tested and compared.The results show that the shear yield stress under compression and compression-shear coupling increases with the growth of compressive yield stress,while the shear yield stress increases slightly under shear.The order of shear yield stress from low to high under different thickening conditions is shear,compression,and compression-shear coupling.Under compression and compression-shear coupling,the concentration first rapidly increases with the growth of compressive yield stress and then slowly increases,while concentration increases slightly under shear.The order of concentration from low to high under different thickening conditions is shear,compression,and compression-shear coupling.Finally,the evolution mechanism of the flocs and drainage channels during the thickening of the thickened tailings under different experimental conditions was revealed.
基金supported by the National Natural Science Foundation of China(Grant No.51974173)the Natural Science Foundation of Shandong Province,China(Grant No.ZR2020QD122).
文摘The shear failure of intact rock under thermo-mechanical(TM)coupling conditions is common,such as in enhanced geothermal mining and deep mine construction.Under the effect of a continuous engineering disturbance,shear-formed fractures are prone to secondary instability,posing a severe threat to deep engineering.Although numerous studies regarding three-dimensional(3D)morphologies of fracture surfaces have been conducted,the understanding of shear-formed fractures under TM coupling conditions is limited.In this study,direct shear tests of intact granite under various TM coupling conditions were conducted,followed by 3D laser scanning tests of shear-formed fractures.Test results demonstrated that the peak shear strength of intact granite is positively correlated with the normal stress,whereas it is negatively correlated with the temperature.The internal friction angle and cohesion of intact granite significantly decrease with an increase in the temperature.The anisotropy,roughness value,and height of the asperities on the fracture surfaces are reduced as the normal stress increases,whereas their variation trends are the opposite as the temperature increases.The macroscopic failure mode of intact granite under TM coupling conditions is dominated by mixed tensileeshear and shear failures.As the normal stress increases,intragranular fractures are developed ranging from a local to a global distribution,and the macroscopic failure mode of intact granite changes from mixed tensileeshear to shear failure.Finally,3D morphological characteristics of the asperities on the shear-formed fracture surfaces were analyzed,and a quadrangular pyramid conceptual model representing these asperities was proposed and sufficiently verified.
基金China National 111 Project Under Grant No. B07019.
文摘Underwater cylindrical shell structures have been found a wide of application in many engineering fields, such as the element of marine, oil platforms, etc. The coupled vibration analysis is a hot issue for these underwater structures. The vibration characteristics of underwater structures are influenced not only by hydrodynamic pressure but also by hydrostatic pressure corresponding to different water depths. In this study, an acoustic finite element method was used to evaluate the underwater structures. Taken the hydrostatic pressure into account in terms of initial stress stiffness, an acoustical fluid-structure coupled analysis of underwater cylindrical shells has been made to study the effect of hydrodynamic pressures on natural frequency and sound radiation. By comparing with the frequencies obtained by the acoustic finite element method and by the added mass method based on the Bessel function, the validity of present analysis was checked. Finally, test samples of the sound radiation of stiffened cylindrical shells were acquired by a harmonic acoustic analysis. The results showed that hydrostatic pressure plays an important role in determining a large submerged body motion, and the characteristics of sound radiation change with water depth. Furthermore, the analysis methods and the results are of significant reference value for studies of other complicated submarine structures.
基金This work was supported by the National Natural Science Foundation of China(Grant No.12172160)Shenzhen Science and Technology Program(Grant No.JCYJ20220818100600002)+1 种基金South-ern University of Science and Technology(Grant No.Y01326127)the Department of Science and Technology of Guangdong Province(Grant Nos.2020B1212030001 and 2021QN020642).
文摘We propose a suite of strategies for the parallel solution of fully implicit monolithic fluid-structure interaction(FSI).The solver is based on a modeling approach that uses the velocity and pressure as the primitive variables,which offers a bridge between computational fluid dynamics(CFD)and computational structural dynamics.The spatiotemporal discretization leverages the variational multiscale formulation and the generalized-αmethod as a means of providing a robust discrete scheme.In particular,the time integration scheme does not suffer from the overshoot phenomenon and optimally dissipates high-frequency spurious modes in both subproblems of FSI.Based on the chosen fully implicit scheme,we systematically develop a combined suite of nonlinear and linear solver strategies.Invoking a block factorization of the Jacobian matrix,the Newton-Raphson procedure is reduced to solving two smaller linear systems in the multi-corrector stage.The first is of the elliptic type,indicating that the algebraic multigrid method serves as a well-suited option.The second exhibits a two-by-two block structure that is analogous to the system arising in CFD.Inspired by prior studies,the additive Schwarz domain decomposition method and the block-factorization-based preconditioners are invoked to address the linear problem.Since the number of unknowns matches in both subdomains,it is straightforward to balance loads when parallelizing the algorithm for distributed-memory architectures.We use two representative FSI benchmarks to demonstrate the robustness,efficiency,and scalability of the overall FSI solver framework.In particular,it is found that the developed FSI solver is comparable to the CFD solver in several aspects,including fixed-size and isogranular scalability as well as robustness.
基金the financial support from the National Natural Science Foundation of China(No.52206091)the Aeronautical Science Foundation of China(No.201928052008)the Natural Science Foundation of Jiangsu Province,China(No.BK20210303)。
文摘Rotating machinery in the aviation industry is increasingly embracing high speeds and miniaturization,and foil dynamic pressure gas bearing has great application value due to its self-lubrication and self-adaptive deformation characteristics.This study explores the interaction mechanism between micro-scale variable-sectional shearing flow with hyper-rotation speeds and a three-layer elastic foil assembly through bidirectional aero-elastic coupling in a Multi-layer Thrust Gas Foil Bearing(MTGFB).The bearing capacity of the MTGFB varies non-linearly with the decrease of gas film clearance,while the collaborative deformation of the three-layer elastic foil assembly can deal with different load conditions.As the load capacity increases,the enhanced dynamic pressure effect causes the top foil to evolve from a single arch to multiple arches.The hydrodynamic effects in the gas film evolve to form multiple segmented wedges with different pitch ratios,while the peak pressure of the gas film always occurs near the vaults of the top foil.As the rotational speed frequency approaches the natural frequency,the resonance of the gas film and elastic foil assembly system occurs,and a phase delay occurs between the pressure pulsation and the vibration of foils.The load capacity of the MTGFB also depends on the elastic moduli of the elastic foil assembly.Increasing the elastic modulus decreases the deformation amplitude of the top foil,whereas it increases those of the backboard and middle foil,increasing the load capacity.
文摘Replaceable flexural and shear fuse-type coupling beams are used in hybrid coupled shear wall(HCSW)systems,enabling concrete buildings to be promptly recovered after severe earthquakes.This study aimed to analytically evaluate the seismic behavior of flexural and shear fuse beams situated in short-,medium-and high-rise RC buildings that have HCSWs.Three building groups hypothetically located in a high seismic hazard zone were studied.A series of 2D nonlinear time history analyses was accomplished in OpenSees,using the ground motion records scaled at the design basis earthquake level.It was found that the effectiveness of fuses in HCSWs depends on various factors such as size and scale of the building,allowable rotation value,inter-story drift ratio,residual drift quantity,energy dissipation value of the fuses,etc.The results show that shear fuses better meet the requirements of rotations and drifts.In contrast,flexural fuses dissipate more energy,but their sectional stiffness should increase to meet other requirements.It was concluded that adoption of proper fuses depends on the overall scale of the building and on how associated factors are considered.
基金supported by the Basic Ability Improvement Project for Young and Middle-Aged Teachers in Guangxi Universities,Project No.2021KY0792.
文摘The development of thermal stress in the exhaust manifold of a gasoline engine is considered.The problem is addresses in the frame of a combined approach wherefluid and structure are coupled using the GT-POWER and STAR-CCM+software.First,the external characteristic curve of the engine is compared with a one-dimen-sional simulation model,then the parameters of the model are modified until the curve matches the available experimental values.GT-POWER is then used to transfer the inlet boundary data under transient conditions to STAR-CCM+in real-time.The temperature profiles of the inner and outer walls of the exhaust manifold are obtained in this way,together with the thermal stress and thermal deformation of the exhaust manifold itself.Using this information,the original model is improved through the addition of connections.Moreover,the local branch pipes are optimized,leading to significant improvements in terms of thermal stress and thermal deforma-tion of the exhaust manifold(a 7%reduction in the maximum thermal stress).
基金This study is supported by the National Natural Science Foundation of China (Grant No50579046) the Science Foundation of Tianjin Municipal Commission of Science and Technology (Grant No043114711)
文摘This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the rigid structure is taken as ＂fictitious＂ fluid with zero strain rate. Both fluid and structure are described by velocity and pressure. The whole domain, including fluid region and structure region, is modeled by the incompressible Navier-Stokes equations which are discretized with fixed Eulerian mesh. However, to keep the structure＇ s rigid body shape and behavior, a rigid body constraint is enforced on the ＂fictitious＂ fluid domain by use of the Distributed Lagrange Multipher/Fictitious Domain （DLM/ FD） method which is originally introduced to solve particulate flow problems by Glowinski et al. For the verification of the model presented herein, a 2D numerical wave tank is established to simulate small amplitude wave propagations, and then numerical results are compared with analytical solutions. Finally, a 2D example of fluid-structure interaction under wave dynamic forces provides convincing evidences for the method excellent solution quality and fidelity.
文摘Fluid-structure interaction （FSI） is essentially a dynamic phenomenon and always exists in fluid-filled pipe system. The four-equation model, which has been proved to be effective to describe and predict the phenomenon of FSI due to friction coupling and Poisson coupling being taken into account, is utilized to describe the FSI of fluid-filled pipe system. Terse compatibility equations are educed by the method of characteristics （MOC） to describe the fluid-filled pipe system. To shorten computing time needed to get the solutions under the condition of keeping accuracy requirement, two steps are adopted, firstly the time step Δt and divided number of the straight pipe are optimized, sec-ondly the mesh spacing Δz close to boundary is subdivided in several submeshes automatically ac-cording to the speed gradient of fluid. The mathematical model and arithmetic are validated by com-parisons between simulation solutions of two straight pipe systems and experiment known from lit-erature.
文摘A novel algorithm is proposed for the simulation of fluid-structure interaction problems.In particular,much attention is paid to natural phenomena such as debris flow.The fluid part(debris flow fluid)is simulated in the framework of the smoothed particle hydrodynamics(SPH)approach,while the solid part(downstream obstacles)is treated using the finite element method(FEM).Fluid-structure coupling is implemented through dynamic boundary conditions.In particular,the software“TensorFlow”and an algorithm based on Python are combined to conduct the required calculations.The simulation results show that the dynamics of viscous and non-viscous debris flows can be extremely different when there are obstacles in the downstream direction.The implemented SPH-FEM coupling method can simulate the fluid-structure coupling problem with a reasonable approximation.
基金supported by the National Natural Science Foundation of China (No. 50708015)the Program for New Century Excellent Talents in University (No. NCET-06-0270), China
文摘To analyze wind-induced response characteristics of a wind turbine tower more accurately, the blade-tower coupling effect was investigated. The mean wind velocity of the rotating blades and tower was simulated according to wind shear effects, and the fluctuating wind velocity time series of the wind turbine were simulated by a harmony superposition method. A dynamic finite element method (FEM) was used to calculate the wind-induced response of the blades and tower. Wind-induced responses of the tower were calculated in two cases (one included the blade-tower coupling effect, and the other only added the mass of blades and the hub at the top of the tower), and then the maximal displacements at the top of the tower of the tow cases were compared with each other. As a result of the influence of the blade-tower coupling effect and the total base shear of the blades, the maximal displacement of the first case increased nearly by 300% compared to the second case. To obtain more precise analysis, the blade-tower coupling effect and the total base shear of the blades should be considered simultaneously in the design of wind turbine towers.
