A seismic-induced landslide is a common geological catastrophe that occurs in nature.The Wangjiayan landslide,which was triggered by the Wenchuan earthquake,is a typical case in point.The Wanjiayan landslide caused ma...A seismic-induced landslide is a common geological catastrophe that occurs in nature.The Wangjiayan landslide,which was triggered by the Wenchuan earthquake,is a typical case in point.The Wanjiayan landslide caused many casualties and resulted in enormous property loss.This study constructs a simple surficial failure model based on the upper bound approach of three-dimensional(3D)limit analysis to evaluate the slope stability of the Wangjiayan case,while a traditional two-dimensional(2D)analysis is also conducted as a reference for comparison with the results of the 3D analysis.A quasi-static calculation is used to study the effect of the earthquake in terms of horizontal ground acceleration,while a parametric study is conducted to evaluate the critical cohesion of slopes.Rather than employing a 3D analysis,using the 2D analysis yields an underestimation regarding the safety factor.In the Wangjiayan landslide,the difference in the factors of safety between the 3D and 2D analyses can reach 20%.The sliding surface morphology as determined by the 3D method is similar to actual morphology,and the parameters of both are also compared to analyze the reliability of the proposed 3D method.展开更多
Urban tunnels are often built in multiples,yet the influence of surface structures on their support demands is frequently overlooked,highlighting the need for refined design under realistic loading.This study addresse...Urban tunnels are often built in multiples,yet the influence of surface structures on their support demands is frequently overlooked,highlighting the need for refined design under realistic loading.This study addresses this gap by examining the response of a uniformly loaded strip footing above twin horseshoe tunnels with structural linings.A non-dimensional Footing Stability Improvement Factor(I_(f)),defined as the ratio of ultimate bearing capacity with and without lining to that of a footing on intact rock—is evaluated using adaptive finite element limit analysis.The effects of normalized vertical depth of tunnel,normalised horizontal spacing of tunnel,and lining-thickness factor are examined,providing key insights into foundation-tunnel interaction and support optimization.This study finds that,for unlined twin tunnels,footing-tunnel interaction becomes negligible when embedment exceeds four times of footing widths and tunnel spacing surpasses eleven times of footing widths,beyond which additional support is unnecessary.For lined tunnels,the tunnel-lining thickness required to maintain footing stability equivalent to that on intact rock decreases rapidly with embedment depth.At a tunnel spacing equal to onehalf of the footing width,the critical lining-thickness ratio is about six percent of the footing width at a tunnel depth of half the footing width,reducing to approximately five,four,and three percent at depths of one,one-and-a-half,and two footing widths,and asymptotically approaching nearly one percent for embedment depths exceeding two-and-a-half footing widths.Additionally,rock mass parameters exhibit minimal impact(±2%)on I_(f).The results show that failure mechanisms in footing-tunnel systems are primarily governed by tunnels-positional parmeters,and lining thickness.The study presents a mechanistic framework that delineates the critical depth-spacing-thickness regimes controlling twin-tunnel behavior,providing engineers with clear design boundaries and performance-based lining guidelines to improve the stability of both tunnels and overlying foundations.展开更多
The three-dimensional spectral analysis method was applied to airglow data from September 2023 to August 2024 derivedfrom an OH airglow imager located at the Hejing station (42.79°N, 83.73°E) to study the pr...The three-dimensional spectral analysis method was applied to airglow data from September 2023 to August 2024 derivedfrom an OH airglow imager located at the Hejing station (42.79°N, 83.73°E) to study the propagation characteristics of gravity waves(GWs) over Northwest China. We found that obvious seasonal variations occur in the propagation of GWs. In spring, GWs mainlypropagate in the northeast direction. In summer and autumn, GWs mainly propagate in the north direction. However, GWs mainlypropagate in the south direction in winter. The direction of GW propagation in the zonal direction is controlled by the wind-filteringeffect, whereas the north–south meridional direction is mainly determined by the location of the wave source. We found that the averageenergy spectrum exhibits a 10%–20% higher intensity in summer and winter compared with spring and autumn. For the first time, wereport the seasonal variation characteristics of GWs over the inland areas of Northwest China, which is of great significance forunderstanding the regional distribution characteristics of GWs.展开更多
A method of three-dimensional loaded slope stability for anisotropic and nonhomogeneous slopes was presented based on the upper-bound theorem of the limit analysis approach. The approach can be considered as a modific...A method of three-dimensional loaded slope stability for anisotropic and nonhomogeneous slopes was presented based on the upper-bound theorem of the limit analysis approach. The approach can be considered as a modification and extension of the solutions. The influences of friction angle, anisotropy factor, nonhomogeneous factor, slope angle, ratio of width to depth, and load on the slope crest were investigated. The results show that solutions are suitable to deal with the purely cohesive soils and frictional/cohesive soils, isotropic and anisotropic, homogeneous and nonhomogeneous, loaded and unloaded cases.展开更多
A generalized nonlinear Baker failure criterion is employed with the upper bound limit analysis to study the surrounding rock stability of underground cavities. A three-dimensional(3D) failure mode is established by e...A generalized nonlinear Baker failure criterion is employed with the upper bound limit analysis to study the surrounding rock stability of underground cavities. A three-dimensional(3D) failure mode is established by extending the two-dimensional(2D) failure mode, which offers an upper bound expression of the surrounding rock pressure. This method is validated with a series of examples before the influence of four parameters of scale parameter, curvature parameter, shift parameter and lateral pressure coefficient, on the surrounding rock pressure is analyzed. According to these results, failure ranges of the underground cavities are determined. The following conclusions are reached:(1) the proposed approach is more accurate to predict surrounding rock pressure than the Mohr-Coulomb failure criterion;(2) the surrounding rock with large scale parameter, curvature parameter, shift parameter, and lateral pressure coefficient can lead to a more stable underground cavity;(3) the failure range in 3D mode can be predicted according to the upper bound solutions.展开更多
In the framework of upper bound theorem of limit analysis, the progressive collapse of shallow rectangular tunnels with double-layer rock mass has been theoretically analyzed based on the three-dimensional (3D) veloci...In the framework of upper bound theorem of limit analysis, the progressive collapse of shallow rectangular tunnels with double-layer rock mass has been theoretically analyzed based on the three-dimensional (3D) velocity discontinuity surfaces. According to the virtual work principle, the difference theorem and the variation method, the collapse surface of double-layer rock mass is determined based on the Hoek-Brown failure criterion. The formula can be degenerated to a single-layer rock collapsing problem when the rock mass is homogeneous. To estimate the validity of the result, the numerical simulation software PLAXIS 3D is used to simulate the collapse of shallow tunnels with double-layer rock mass, and the comparative analysis shows that numerical results are in good agreement with upper-bound solutions. According to the results of parametric analysis, the potential range of collapse of a double-layer rock mass above a shallow cavity decreases with a decrease in A1/A2,σci1/σci2 and σtm1/σtm2 and an increase in B1/B2,γ1/γ2. The range will decrease with a decrease in support pressure q and increase with a decrease in surface overload σs. Therefore, reinforced supporting is beneficial to improve the stability of the cavity during actual construction.展开更多
Based on the lower bound theorem of limit analysis, a solution procedure for limit analysis of three_dimensional elastoplastic structures was established using conventional boundary element method (BEM). The elastic s...Based on the lower bound theorem of limit analysis, a solution procedure for limit analysis of three_dimensional elastoplastic structures was established using conventional boundary element method (BEM). The elastic stress field for lower bound limit analysis was computed directly by three_dimensional boundary element method (3_D BEM). The self_equilibrium stress field was constructed by the linear combination of several self_equilibrium “basis vectors” which can be computed by elastic_plastic incremental iteration of 3_D BEM analysis. The lower bound limit analysis problem was finally reduced to a series of nonlinear programming sub_problems with relatively few optimal variables. The complex method was used to solve the nonlinear programming sub_problems. The numerical results show that the present solution procedure has good accuracy and high efficiency.展开更多
The theory of limit analysis is presented for a three-dimensional stability problem of excavation. In frictional soil, the failure surface has the shape of logarithm helicoid, with the outline profile defined by a log...The theory of limit analysis is presented for a three-dimensional stability problem of excavation. In frictional soil, the failure surface has the shape of logarithm helicoid, with the outline profile defined by a log- spiral curve. The internal dissipation rate of energy caused by soil cohesion and gravity power of the failure soil is obtained through theoretical derivation. By solving the energy balance equation, the stability factor for the excavation is obtained. Influence of the ratio of width to height, the slope angle, and the top angle on the stability is examined. Numerical results of the proposed algorithm are presented in the form of non dimensional graph. Examples illustrate the practical use of the results.展开更多
DAYU3D is a modern three-dimensional(3D)computer code for thermal-hydraulic design and accident analysis in hightemperature gas-cooled reactors(HTGRs),developed by the Institute of Nuclear and New Energy Technology(IN...DAYU3D is a modern three-dimensional(3D)computer code for thermal-hydraulic design and accident analysis in hightemperature gas-cooled reactors(HTGRs),developed by the Institute of Nuclear and New Energy Technology(INET)at Tsinghua University.Compared to the traditional codes like TINTE,the DAYU3D code has advantages due to its refined framework,improved models,and more efficient algorithms.It is able to simulate the continuous movement of control rods and is more rigorous in treating radiation heat transfer and the break mass flow.Advanced computational methods significantly improve the computational efficiency of DAYU3D,achieving a time reduction of over 60%compared to TINTE.Extensive verification and validation with more than 100 cases demonstrate that DAYU3D is promising for HTGR 3D thermal-hydraulic design and accident analyses.展开更多
Based on the upper bound of limit analysis, the plane-strain analysis of the slopes reinforced with a row of piles to the 3D case was extended. A 3D rotational failure mechanism was adopted to yield the upper bound of...Based on the upper bound of limit analysis, the plane-strain analysis of the slopes reinforced with a row of piles to the 3D case was extended. A 3D rotational failure mechanism was adopted to yield the upper bound of the factor of safety. Parametric studies were carried out to explore the end effects of the slope failures and the effects of the pile location and diameter on the safety of the reinforced slopes. The results demonstrate that the end effects nearly have no effects on the most suitable location of the installed piles but have significant influence on the safety of the slopes. For a slope constrained to a narrow width, the slope becomes more stable owing to the contribution of the end effects. When the slope is reinforced with a row of piles in small space between piles, the effects of group piles are significant for evaluating the safety of slopes. The presented method is more appropriate for assessing the stability of slopes reinforced with piles and can be also utilized in the design of plies stabilizing the unstable slopes.展开更多
The increasing occurrence of sinkholes caused by water main bursts has attracted significant research attention in recent years.This study addresses the gap in evaluating soil blowout stability resulting from water ma...The increasing occurrence of sinkholes caused by water main bursts has attracted significant research attention in recent years.This study addresses the gap in evaluating soil blowout stability resulting from water main failures by investigating the three-dimensional stability of blowouts with circular,hemispherical,and spherical openings.Advanced finite element limit analysis(FELA)combined with adaptive meshing is employed to analyze critical factors,including soil cover depth,surcharge pressure,and internal water pressure,that contribute to blowout failure.In addition,dimensionless ratios are used throughout the paper to assess the influence of these factors.Numerical findings are rigorously validated,ensuring reliability and accuracy.Practical design charts are provided to accommodate a wide range of design scenarios,offering valuable guidance for engineers.This study introduces a pioneering sinkhole simulation methodology,leading to the understanding of three-dimensional blowout stability mechanisms.展开更多
Based on nonlinear failure criterion,a three-dimensional failure mechanism of the possible collapse of deep tunnel is presented with limit analysis theory.Support pressure is taken into consideration in the virtual wo...Based on nonlinear failure criterion,a three-dimensional failure mechanism of the possible collapse of deep tunnel is presented with limit analysis theory.Support pressure is taken into consideration in the virtual work equation performed under the upper bound theorem.It is necessary to point out that the properties of surrounding rock mass plays a vital role in the shape of collapsing rock mass.The first order reliability method and Monte Carlo simulation method are then employed to analyze the stability of presented mechanism.Different rock parameters are considered random variables to value the corresponding reliability index with an increasing applied support pressure.The reliability indexes calculated by two methods are in good agreement.Sensitivity analysis was performed and the influence of coefficient variation of rock parameters was discussed.It is shown that the tensile strength plays a much more important role in reliability index than dimensionless parameter,and that small changes occurring in the coefficient of variation would make great influence of reliability index.Thus,significant attention should be paid to the properties of surrounding rock mass and the applied support pressure to maintain the stability of tunnel can be determined for a given reliability index.展开更多
A three-dimensional finite element simulation was carried out to investigate the effects of tunnel construction on nearby pile foundation.The displacement controlled model (DCM) was used to simulate the tunneling-indu...A three-dimensional finite element simulation was carried out to investigate the effects of tunnel construction on nearby pile foundation.The displacement controlled model (DCM) was used to simulate the tunneling-induced volume loss effects.The numerical model was verified based on the results of a centrifuge test and a set of parametric studies was implemented based on this model.There is good agreement between the trend of the results of the centrifuge test and the present model.The results of parametric studies show that the tunnelling-induced pile internal force and deformation depend mainly on the pile?tunnel distance,the pile length to tunnel depth ratio and the volume loss.Two different zones are separated by a 45° line projected from the tunnel springline.Within the zone of influence,the pile is subjected to tensile force and large settlement;whereas outside the zone of influence,dragload and small settlement are induced.It is also established that the impact of tunnelling on a pile group is substantially smaller as compared with a single pile in the same location with the rear pile in a group,demonstrating a positive pile group effect.展开更多
In the last decade, three dimensional discontin- uous deformation analyses （3D DDA） has attracted more and more attention of researchers and geotechnical engineers worldwide. The original DDA formulation utilizes a ...In the last decade, three dimensional discontin- uous deformation analyses （3D DDA） has attracted more and more attention of researchers and geotechnical engineers worldwide. The original DDA formulation utilizes a linear displacement function to describe the block movement and deformation, which would cause block expansion under rigid body rotation and thus limit its capability to model block de- formation. In this paper, 3D DDA is coupled with tetrahe- dron finite elements to tackle these two problems. Tetrahe- dron is the simplest in the 3D domain and makes it easy to implement automatic discretization, even for complex topol- ogy shape. Furthermore, element faces will remain planar and element edges will remain straight after deformation for tetrahedron finite elements and polyhedral contact detection schemes can be used directly. The matrices of equilibrium equations for this coupled method are given in detail and an effective contact searching algorithm is suggested. Valida- tion is conducted by comparing the results of the proposed coupled method with that of physical model tests using one of the most common failure modes, i.e., wedge failure. Most of the failure modes predicted by the coupled method agree with the physical model results except for 4 cases out of the total 65 cases. Finally, a complex rockslide example demon- strates the robustness and versatility of the coupled method.展开更多
Simulated photo-degradation of fluorescent dissolved organic matter(FDOM) in Lake Baihua(BH) and Lake Hongfeng(HF) was investigated with three-dimensional excitationemission matrix(3 DEEM) fluorescence combined with t...Simulated photo-degradation of fluorescent dissolved organic matter(FDOM) in Lake Baihua(BH) and Lake Hongfeng(HF) was investigated with three-dimensional excitationemission matrix(3 DEEM) fluorescence combined with the fluorescence regional integration(FRI),parallel factor(PARAFAC) analysis,and multi-order kinetic models.In the FRI analysis,fulvic-like and humic-like materials were the main constituents for both BH-FDOM and HF-FDOM.Four individual components were identified by use of PARAFAC analysis as humic-like components(C1),fulvic-like components(C2),protein-like components(C3) and unidentified components(C4).The maximum 3 DEEM fluorescence intensity of PARAFAC components C1-C3 decreased by about 60%,70% and 90%,respectively after photo-degradation.The multi-order kinetic model was acceptable to represent the photo-degradation of FDOM with correlation coefficient(Radj2)(0.963-0.998).The photo-degradation rate constants(kn) showed differences of three orders of magnitude,from 1.09 × 10-6 to 4.02 × 10-4 min-1,and half-life of multi-order model(T1/2n)ranged from 5.26 to 64.01 min.The decreased values of fluorescence index(FI) and biogenic index(BI),the fact that of percent fluorescence response parameter of Region I(PⅠ,n) showed the greatest change ratio,followed by percent fluorescence response parameter of Region II(PⅡ,n,while the largest decrease ratio was found for C3 components,and the lowest T1/2n was observed for C3,indicated preferential degradation of protein-like materials/components derived from biological sources during photodegradation.This research on the degradation of FDOM by 3 DEEM/FRI-PARAFAC would be beneficial to understanding the photo-degradation of FD OM in natural environments and accurately predicting the environmental behaviors of contaminants in the presence of FDOM.展开更多
A three-dimensional analysis of a simply-supported functionally graded rectangular plate with an arbitrary distribution of material properties is made using a simple and effective method based on the Haar wavelet. Wit...A three-dimensional analysis of a simply-supported functionally graded rectangular plate with an arbitrary distribution of material properties is made using a simple and effective method based on the Haar wavelet. With good features in treating singularities, Haar series solution converges rapidly for arbitrary distributions, especially for the case where the material properties change rapidly in some regions. Through numerical examples the influences of the ratio of material constants on the top and bottom surfaces and different material gradient distributions on the structural response of the plate to mechanical stimuli are studied.展开更多
A three-dimensional analysis model based on the finite element method （FEM） is developed, which can derive the evolution and distribution characteristics of heat flux deposited on the divertor plate from the surface...A three-dimensional analysis model based on the finite element method （FEM） is developed, which can derive the evolution and distribution characteristics of heat flux deposited on the divertor plate from the surface temperature measured by infrared thermography diagnostics. The numerical simulations of surface heating due to localized power bursts and the power deposition calculations demonstrate that this analysis can provide accurate results and useful information about localized hot spots compared with the normal one- and two-dimensional calculations. In this paper, the details of this three- dimensional analysis are presented, and some results in ohmic heating and electron cyclotron resonant heating （ECRH） discharge on HL-2A are given.展开更多
Using bathymetry and altimetric gravity anomalies, a 1°×9 1° lithospheric effective elastic thickness（Te） model over the Louisville Ridge and its adjacent regions is calculated using the moving window...Using bathymetry and altimetric gravity anomalies, a 1°×9 1° lithospheric effective elastic thickness（Te） model over the Louisville Ridge and its adjacent regions is calculated using the moving window admittance technique. For comparison, three bathymetry models are used： general bathymetric charts of the oceans, SIO V15.1,and BAT_VGG. The results show that BAT_VGG is more suitable for calculating T e than the other two models. T e along the Louisville Ridge was re-evaluated. The southeast of the ridge has a medium Te of 10–20 km, while Te increases dramatically seaward of the Tonga-Kermadec trench as a result of the collision of the Pacific and IndoAustralian plates.展开更多
This paper proposes a based on 3D-VLE （three-dimensional nonlinear viscoelastic theory） three-parameters viscoelastic model for studying the time-dependent behaviour of concrete filled steel tube （CFT） columns. Th...This paper proposes a based on 3D-VLE （three-dimensional nonlinear viscoelastic theory） three-parameters viscoelastic model for studying the time-dependent behaviour of concrete filled steel tube （CFT） columns. The method of 3D-VLE was developed to analyze the effects of concrete creep behavior on CFT structures. After the evaluation of the parameters in the proposed creep model, experimental measurements of two prestressed reinforced concrete beams were used to investigate the creep phenomenon of three CFT columns under long-term axial and eccentric load was investigated. The experimentally obtained time-dependent creep behaviour accorded well with the cu~＇es obtained from the proposed method. Many factors （such as ratio of long-term load to strength, slenderness ratio, steel ratio, and eccentricity ratio） were considered to obtain the regularity of influence of concrete creep on CFT structures. The analytical results can be consulted in the engineering practice and design.展开更多
Based on the working of Lighthill and Hunt et al., in the present paper the author has established the topological rules adapting to analysing the skin-friction lines and the section streamlines in cascades. These rul...Based on the working of Lighthill and Hunt et al., in the present paper the author has established the topological rules adapting to analysing the skin-friction lines and the section streamlines in cascades. These rules are (1) for a rotor cascade without shroud band, the total number of nodal points equals that the saddle points on the skin-friction line vector fields in eachpitch range; (2) for an annular or straight cascade with no-clearances at blade ends, the total number of saddle points is two more than that of nodal points on the skin-friction line fields in a pitch; (3) the total number of saddles in the secondary flow fields on cross-sections in cascade is one less than that of nodes; (4) in the section streamline vector fields on a meridian surface penetrating a flow passage, and on leading and trailing edge sections, the total number of nodes is equal to that of saddles; (5) on the streamline vector fields of a blade-to-blade surface, the total number of nodes is one less than that of saddles.展开更多
基金National Natural Science Foundation of China under Grant Nos.52378335 and 52322808.
文摘A seismic-induced landslide is a common geological catastrophe that occurs in nature.The Wangjiayan landslide,which was triggered by the Wenchuan earthquake,is a typical case in point.The Wanjiayan landslide caused many casualties and resulted in enormous property loss.This study constructs a simple surficial failure model based on the upper bound approach of three-dimensional(3D)limit analysis to evaluate the slope stability of the Wangjiayan case,while a traditional two-dimensional(2D)analysis is also conducted as a reference for comparison with the results of the 3D analysis.A quasi-static calculation is used to study the effect of the earthquake in terms of horizontal ground acceleration,while a parametric study is conducted to evaluate the critical cohesion of slopes.Rather than employing a 3D analysis,using the 2D analysis yields an underestimation regarding the safety factor.In the Wangjiayan landslide,the difference in the factors of safety between the 3D and 2D analyses can reach 20%.The sliding surface morphology as determined by the 3D method is similar to actual morphology,and the parameters of both are also compared to analyze the reliability of the proposed 3D method.
文摘Urban tunnels are often built in multiples,yet the influence of surface structures on their support demands is frequently overlooked,highlighting the need for refined design under realistic loading.This study addresses this gap by examining the response of a uniformly loaded strip footing above twin horseshoe tunnels with structural linings.A non-dimensional Footing Stability Improvement Factor(I_(f)),defined as the ratio of ultimate bearing capacity with and without lining to that of a footing on intact rock—is evaluated using adaptive finite element limit analysis.The effects of normalized vertical depth of tunnel,normalised horizontal spacing of tunnel,and lining-thickness factor are examined,providing key insights into foundation-tunnel interaction and support optimization.This study finds that,for unlined twin tunnels,footing-tunnel interaction becomes negligible when embedment exceeds four times of footing widths and tunnel spacing surpasses eleven times of footing widths,beyond which additional support is unnecessary.For lined tunnels,the tunnel-lining thickness required to maintain footing stability equivalent to that on intact rock decreases rapidly with embedment depth.At a tunnel spacing equal to onehalf of the footing width,the critical lining-thickness ratio is about six percent of the footing width at a tunnel depth of half the footing width,reducing to approximately five,four,and three percent at depths of one,one-and-a-half,and two footing widths,and asymptotically approaching nearly one percent for embedment depths exceeding two-and-a-half footing widths.Additionally,rock mass parameters exhibit minimal impact(±2%)on I_(f).The results show that failure mechanisms in footing-tunnel systems are primarily governed by tunnels-positional parmeters,and lining thickness.The study presents a mechanistic framework that delineates the critical depth-spacing-thickness regimes controlling twin-tunnel behavior,providing engineers with clear design boundaries and performance-based lining guidelines to improve the stability of both tunnels and overlying foundations.
基金supported by the National Science Foundation of China(Grant Nos.42374205 and 41974179)the Specialized Research Fund of the National Space Science Center,Chinese Academy of Sciences(Grant No.E4PD3010)supported by the Specialized Research Fund for State Key Laboratories.
文摘The three-dimensional spectral analysis method was applied to airglow data from September 2023 to August 2024 derivedfrom an OH airglow imager located at the Hejing station (42.79°N, 83.73°E) to study the propagation characteristics of gravity waves(GWs) over Northwest China. We found that obvious seasonal variations occur in the propagation of GWs. In spring, GWs mainlypropagate in the northeast direction. In summer and autumn, GWs mainly propagate in the north direction. However, GWs mainlypropagate in the south direction in winter. The direction of GW propagation in the zonal direction is controlled by the wind-filteringeffect, whereas the north–south meridional direction is mainly determined by the location of the wave source. We found that the averageenergy spectrum exhibits a 10%–20% higher intensity in summer and winter compared with spring and autumn. For the first time, wereport the seasonal variation characteristics of GWs over the inland areas of Northwest China, which is of great significance forunderstanding the regional distribution characteristics of GWs.
基金Projects(41002095,41172251,41272317) supported by the National Natural Science Foundation of China
文摘A method of three-dimensional loaded slope stability for anisotropic and nonhomogeneous slopes was presented based on the upper-bound theorem of the limit analysis approach. The approach can be considered as a modification and extension of the solutions. The influences of friction angle, anisotropy factor, nonhomogeneous factor, slope angle, ratio of width to depth, and load on the slope crest were investigated. The results show that solutions are suitable to deal with the purely cohesive soils and frictional/cohesive soils, isotropic and anisotropic, homogeneous and nonhomogeneous, loaded and unloaded cases.
基金Projects(51679117,11772358,51774322,51474249,51404179,51274249)supported by the National Natural Science Foundation of China。
文摘A generalized nonlinear Baker failure criterion is employed with the upper bound limit analysis to study the surrounding rock stability of underground cavities. A three-dimensional(3D) failure mode is established by extending the two-dimensional(2D) failure mode, which offers an upper bound expression of the surrounding rock pressure. This method is validated with a series of examples before the influence of four parameters of scale parameter, curvature parameter, shift parameter and lateral pressure coefficient, on the surrounding rock pressure is analyzed. According to these results, failure ranges of the underground cavities are determined. The following conclusions are reached:(1) the proposed approach is more accurate to predict surrounding rock pressure than the Mohr-Coulomb failure criterion;(2) the surrounding rock with large scale parameter, curvature parameter, shift parameter, and lateral pressure coefficient can lead to a more stable underground cavity;(3) the failure range in 3D mode can be predicted according to the upper bound solutions.
基金Projects(51478477,51878074)supported by the National Natural Science Foundation of ChinaProject(2017-123-033)supported by the Guizhou Provincial Department of Transportation Foundation,ChinaProjects(2018zzts663,2018zzts656)supported by the Fundamental Research Funds for the Central Universities,China
文摘In the framework of upper bound theorem of limit analysis, the progressive collapse of shallow rectangular tunnels with double-layer rock mass has been theoretically analyzed based on the three-dimensional (3D) velocity discontinuity surfaces. According to the virtual work principle, the difference theorem and the variation method, the collapse surface of double-layer rock mass is determined based on the Hoek-Brown failure criterion. The formula can be degenerated to a single-layer rock collapsing problem when the rock mass is homogeneous. To estimate the validity of the result, the numerical simulation software PLAXIS 3D is used to simulate the collapse of shallow tunnels with double-layer rock mass, and the comparative analysis shows that numerical results are in good agreement with upper-bound solutions. According to the results of parametric analysis, the potential range of collapse of a double-layer rock mass above a shallow cavity decreases with a decrease in A1/A2,σci1/σci2 and σtm1/σtm2 and an increase in B1/B2,γ1/γ2. The range will decrease with a decrease in support pressure q and increase with a decrease in surface overload σs. Therefore, reinforced supporting is beneficial to improve the stability of the cavity during actual construction.
文摘Based on the lower bound theorem of limit analysis, a solution procedure for limit analysis of three_dimensional elastoplastic structures was established using conventional boundary element method (BEM). The elastic stress field for lower bound limit analysis was computed directly by three_dimensional boundary element method (3_D BEM). The self_equilibrium stress field was constructed by the linear combination of several self_equilibrium “basis vectors” which can be computed by elastic_plastic incremental iteration of 3_D BEM analysis. The lower bound limit analysis problem was finally reduced to a series of nonlinear programming sub_problems with relatively few optimal variables. The complex method was used to solve the nonlinear programming sub_problems. The numerical results show that the present solution procedure has good accuracy and high efficiency.
基金the National Natural Science Foundation of China(Nos.41002095,41172251 and 41272317)
文摘The theory of limit analysis is presented for a three-dimensional stability problem of excavation. In frictional soil, the failure surface has the shape of logarithm helicoid, with the outline profile defined by a log- spiral curve. The internal dissipation rate of energy caused by soil cohesion and gravity power of the failure soil is obtained through theoretical derivation. By solving the energy balance equation, the stability factor for the excavation is obtained. Influence of the ratio of width to height, the slope angle, and the top angle on the stability is examined. Numerical results of the proposed algorithm are presented in the form of non dimensional graph. Examples illustrate the practical use of the results.
文摘DAYU3D is a modern three-dimensional(3D)computer code for thermal-hydraulic design and accident analysis in hightemperature gas-cooled reactors(HTGRs),developed by the Institute of Nuclear and New Energy Technology(INET)at Tsinghua University.Compared to the traditional codes like TINTE,the DAYU3D code has advantages due to its refined framework,improved models,and more efficient algorithms.It is able to simulate the continuous movement of control rods and is more rigorous in treating radiation heat transfer and the break mass flow.Advanced computational methods significantly improve the computational efficiency of DAYU3D,achieving a time reduction of over 60%compared to TINTE.Extensive verification and validation with more than 100 cases demonstrate that DAYU3D is promising for HTGR 3D thermal-hydraulic design and accident analyses.
基金Projects(51278382,51479050)supported by the National Natural Science Foundation of ChinaProject(2015CB057901)supported by the National Key Basic Research Program of China+3 种基金Project(201501035-03)supported by the Public Service Sector R&D Project of Ministry of Water Resource of ChinaProject(2014B06814)supported by the Fundamental Research Funds for the Central Universities,ChinaProject(B13024)supported by the"111"ProjectProject(YK913004)supported by the Open Foundation of Key Laboratory of Failure Mechanism and Safety Control Techniques of Earthrock Dam of the Ministry of Water Resources,China
文摘Based on the upper bound of limit analysis, the plane-strain analysis of the slopes reinforced with a row of piles to the 3D case was extended. A 3D rotational failure mechanism was adopted to yield the upper bound of the factor of safety. Parametric studies were carried out to explore the end effects of the slope failures and the effects of the pile location and diameter on the safety of the reinforced slopes. The results demonstrate that the end effects nearly have no effects on the most suitable location of the installed piles but have significant influence on the safety of the slopes. For a slope constrained to a narrow width, the slope becomes more stable owing to the contribution of the end effects. When the slope is reinforced with a row of piles in small space between piles, the effects of group piles are significant for evaluating the safety of slopes. The presented method is more appropriate for assessing the stability of slopes reinforced with piles and can be also utilized in the design of plies stabilizing the unstable slopes.
文摘The increasing occurrence of sinkholes caused by water main bursts has attracted significant research attention in recent years.This study addresses the gap in evaluating soil blowout stability resulting from water main failures by investigating the three-dimensional stability of blowouts with circular,hemispherical,and spherical openings.Advanced finite element limit analysis(FELA)combined with adaptive meshing is employed to analyze critical factors,including soil cover depth,surcharge pressure,and internal water pressure,that contribute to blowout failure.In addition,dimensionless ratios are used throughout the paper to assess the influence of these factors.Numerical findings are rigorously validated,ensuring reliability and accuracy.Practical design charts are provided to accommodate a wide range of design scenarios,offering valuable guidance for engineers.This study introduces a pioneering sinkhole simulation methodology,leading to the understanding of three-dimensional blowout stability mechanisms.
基金Project (2013CB036004) supported by National Basic Research Program of China
文摘Based on nonlinear failure criterion,a three-dimensional failure mechanism of the possible collapse of deep tunnel is presented with limit analysis theory.Support pressure is taken into consideration in the virtual work equation performed under the upper bound theorem.It is necessary to point out that the properties of surrounding rock mass plays a vital role in the shape of collapsing rock mass.The first order reliability method and Monte Carlo simulation method are then employed to analyze the stability of presented mechanism.Different rock parameters are considered random variables to value the corresponding reliability index with an increasing applied support pressure.The reliability indexes calculated by two methods are in good agreement.Sensitivity analysis was performed and the influence of coefficient variation of rock parameters was discussed.It is shown that the tensile strength plays a much more important role in reliability index than dimensionless parameter,and that small changes occurring in the coefficient of variation would make great influence of reliability index.Thus,significant attention should be paid to the properties of surrounding rock mass and the applied support pressure to maintain the stability of tunnel can be determined for a given reliability index.
文摘A three-dimensional finite element simulation was carried out to investigate the effects of tunnel construction on nearby pile foundation.The displacement controlled model (DCM) was used to simulate the tunneling-induced volume loss effects.The numerical model was verified based on the results of a centrifuge test and a set of parametric studies was implemented based on this model.There is good agreement between the trend of the results of the centrifuge test and the present model.The results of parametric studies show that the tunnelling-induced pile internal force and deformation depend mainly on the pile?tunnel distance,the pile length to tunnel depth ratio and the volume loss.Two different zones are separated by a 45° line projected from the tunnel springline.Within the zone of influence,the pile is subjected to tensile force and large settlement;whereas outside the zone of influence,dragload and small settlement are induced.It is also established that the impact of tunnelling on a pile group is substantially smaller as compared with a single pile in the same location with the rear pile in a group,demonstrating a positive pile group effect.
基金supported by the Key Project of Chinese National Programs for Fundamental Research and Development(2010CB731502)the National Natural Science Foundation of China(50978745)
文摘In the last decade, three dimensional discontin- uous deformation analyses （3D DDA） has attracted more and more attention of researchers and geotechnical engineers worldwide. The original DDA formulation utilizes a linear displacement function to describe the block movement and deformation, which would cause block expansion under rigid body rotation and thus limit its capability to model block de- formation. In this paper, 3D DDA is coupled with tetrahe- dron finite elements to tackle these two problems. Tetrahe- dron is the simplest in the 3D domain and makes it easy to implement automatic discretization, even for complex topol- ogy shape. Furthermore, element faces will remain planar and element edges will remain straight after deformation for tetrahedron finite elements and polyhedral contact detection schemes can be used directly. The matrices of equilibrium equations for this coupled method are given in detail and an effective contact searching algorithm is suggested. Valida- tion is conducted by comparing the results of the proposed coupled method with that of physical model tests using one of the most common failure modes, i.e., wedge failure. Most of the failure modes predicted by the coupled method agree with the physical model results except for 4 cases out of the total 65 cases. Finally, a complex rockslide example demon- strates the robustness and versatility of the coupled method.
基金financially supported by the National Natural Science Foundation of China(No.41573130)BNU Interdisciplinary Research Foundation for First-Year Doctoral Candidates(No.BNUXKJC1802)
文摘Simulated photo-degradation of fluorescent dissolved organic matter(FDOM) in Lake Baihua(BH) and Lake Hongfeng(HF) was investigated with three-dimensional excitationemission matrix(3 DEEM) fluorescence combined with the fluorescence regional integration(FRI),parallel factor(PARAFAC) analysis,and multi-order kinetic models.In the FRI analysis,fulvic-like and humic-like materials were the main constituents for both BH-FDOM and HF-FDOM.Four individual components were identified by use of PARAFAC analysis as humic-like components(C1),fulvic-like components(C2),protein-like components(C3) and unidentified components(C4).The maximum 3 DEEM fluorescence intensity of PARAFAC components C1-C3 decreased by about 60%,70% and 90%,respectively after photo-degradation.The multi-order kinetic model was acceptable to represent the photo-degradation of FDOM with correlation coefficient(Radj2)(0.963-0.998).The photo-degradation rate constants(kn) showed differences of three orders of magnitude,from 1.09 × 10-6 to 4.02 × 10-4 min-1,and half-life of multi-order model(T1/2n)ranged from 5.26 to 64.01 min.The decreased values of fluorescence index(FI) and biogenic index(BI),the fact that of percent fluorescence response parameter of Region I(PⅠ,n) showed the greatest change ratio,followed by percent fluorescence response parameter of Region II(PⅡ,n,while the largest decrease ratio was found for C3 components,and the lowest T1/2n was observed for C3,indicated preferential degradation of protein-like materials/components derived from biological sources during photodegradation.This research on the degradation of FDOM by 3 DEEM/FRI-PARAFAC would be beneficial to understanding the photo-degradation of FD OM in natural environments and accurately predicting the environmental behaviors of contaminants in the presence of FDOM.
基金Project supported by the National Natural Sciences Foundation of China(No.10432030).
文摘A three-dimensional analysis of a simply-supported functionally graded rectangular plate with an arbitrary distribution of material properties is made using a simple and effective method based on the Haar wavelet. With good features in treating singularities, Haar series solution converges rapidly for arbitrary distributions, especially for the case where the material properties change rapidly in some regions. Through numerical examples the influences of the ratio of material constants on the top and bottom surfaces and different material gradient distributions on the structural response of the plate to mechanical stimuli are studied.
基金Project supported by the National Natural Science Foundation of China (Grant No. 10805016)the National Magnetic Confinement Fusion Science Program,China (Grant No. 2009GB104008).
文摘A three-dimensional analysis model based on the finite element method （FEM） is developed, which can derive the evolution and distribution characteristics of heat flux deposited on the divertor plate from the surface temperature measured by infrared thermography diagnostics. The numerical simulations of surface heating due to localized power bursts and the power deposition calculations demonstrate that this analysis can provide accurate results and useful information about localized hot spots compared with the normal one- and two-dimensional calculations. In this paper, the details of this three- dimensional analysis are presented, and some results in ohmic heating and electron cyclotron resonant heating （ECRH） discharge on HL-2A are given.
基金supported financially by the Key Foundation of the Institute of Seismology,China Earthquake Administration (No. IS201506205)the National Natural Science Foundation of China (Nos. 41504017, 41204019, 41304003)
文摘Using bathymetry and altimetric gravity anomalies, a 1°×9 1° lithospheric effective elastic thickness（Te） model over the Louisville Ridge and its adjacent regions is calculated using the moving window admittance technique. For comparison, three bathymetry models are used： general bathymetric charts of the oceans, SIO V15.1,and BAT_VGG. The results show that BAT_VGG is more suitable for calculating T e than the other two models. T e along the Louisville Ridge was re-evaluated. The southeast of the ridge has a medium Te of 10–20 km, while Te increases dramatically seaward of the Tonga-Kermadec trench as a result of the collision of the Pacific and IndoAustralian plates.
文摘This paper proposes a based on 3D-VLE （three-dimensional nonlinear viscoelastic theory） three-parameters viscoelastic model for studying the time-dependent behaviour of concrete filled steel tube （CFT） columns. The method of 3D-VLE was developed to analyze the effects of concrete creep behavior on CFT structures. After the evaluation of the parameters in the proposed creep model, experimental measurements of two prestressed reinforced concrete beams were used to investigate the creep phenomenon of three CFT columns under long-term axial and eccentric load was investigated. The experimentally obtained time-dependent creep behaviour accorded well with the cu~＇es obtained from the proposed method. Many factors （such as ratio of long-term load to strength, slenderness ratio, steel ratio, and eccentricity ratio） were considered to obtain the regularity of influence of concrete creep on CFT structures. The analytical results can be consulted in the engineering practice and design.
文摘Based on the working of Lighthill and Hunt et al., in the present paper the author has established the topological rules adapting to analysing the skin-friction lines and the section streamlines in cascades. These rules are (1) for a rotor cascade without shroud band, the total number of nodal points equals that the saddle points on the skin-friction line vector fields in eachpitch range; (2) for an annular or straight cascade with no-clearances at blade ends, the total number of saddle points is two more than that of nodal points on the skin-friction line fields in a pitch; (3) the total number of saddles in the secondary flow fields on cross-sections in cascade is one less than that of nodes; (4) in the section streamline vector fields on a meridian surface penetrating a flow passage, and on leading and trailing edge sections, the total number of nodes is equal to that of saddles; (5) on the streamline vector fields of a blade-to-blade surface, the total number of nodes is one less than that of saddles.
