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.展开更多
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.展开更多
Multi-layer slopes are widely found in clay residue receiving fields.A generalized horizontal slice method(GHSM)for assessing the stability of multi-layer slopes that considers the energy dissipation between adjacent ...Multi-layer slopes are widely found in clay residue receiving fields.A generalized horizontal slice method(GHSM)for assessing the stability of multi-layer slopes that considers the energy dissipation between adjacent horizontal slices is presented.In view of the upper-bound limit analysis theory,the energy equation is derived and the ultimate failure mode is generated by comparing the sliding surface passing through the slope toe(mode A)with that below(mode B).In addition,the influence of the number of slices on the stability coefficients in the GHSM is studied and the stable value is obtained.Compared to the original method(Chen’s method),the GHSM can acquire more precise results,which takes into account the energy dissipation in the inner sliding soil mass.Moreover,the GHSM,limit equilibrium method(LEM)and numerical simulation method(NSM)are applied to analyze the stability of a multi-layer slope with different slope angles and the results of the safety factor and failure mode are very close in each case.The ultimate failure modes are shown to be mode B when the slope angle is not more than 28°.It illustrates that the determination of the ultimate sliding surface requires comparison of multiple failure modes,not only mode A.展开更多
Extended reach wells （ERWs）, especially horizontal extended reach well with a high HD （horizontal displacement） to TVD （true vertical depth） ratio, represent a frontier technology and challenge the drilling limi...Extended reach wells （ERWs）, especially horizontal extended reach well with a high HD （horizontal displacement） to TVD （true vertical depth） ratio, represent a frontier technology and challenge the drilling limitations. Oil and gas reservoir in beaches or lakes and offshore can be effectively exploited by using extended reach drilling （ERD） technology. This paper focuses on the difficult technological problems encountered during exploiting the Liuhua 11-1 oil field in the South China Sea, China. Emphasis is on investigating the key subjects including prediction and control of open hole limit extension in offshore ERD, prediction of casing wear and its prevention and torque reduction, φ244.5mm casing running with floating collars to control drag force, and steerable drilling modes. The basic concept of limit extension in ERD is presented and the prediction method for open hole limit extension is given in this paper. A set of advanced drilling mechanics and control technology has been established and its practical results are verified by field cases. All those efforts may be significant for further investigating and practicing ERD limit theory and control technology in the future.展开更多
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 method of analyzing the stability of twin shallow tunnels was presented using both limit analysis with nonlinear failure criterion and reliability theory.In the condition of nonlinear failure criterion,the critical ...A method of analyzing the stability of twin shallow tunnels was presented using both limit analysis with nonlinear failure criterion and reliability theory.In the condition of nonlinear failure criterion,the critical clear distancesof twin shallow tunnels were obtained by analyzing the change of surrounding pressure.A reliability model was established based on limit state equation,and the failure probability was solved by virtue of Monte Carlo method.Safety factor and corresponding clear distance of different safety levels were obtained by introducing a target reliability index.The scope of clear distance for different safety levels is described,which can be used as a supplement and improvement to the design codes of tunnels.展开更多
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.展开更多
According to the lower-bound theorem of limit analysis the Rigid Finite Element Meth-od(RFEM)is applied to structural limit analysis and the linear programmings for limit analysis are deducedin this paper.Moreover,the...According to the lower-bound theorem of limit analysis the Rigid Finite Element Meth-od(RFEM)is applied to structural limit analysis and the linear programmings for limit analysis are deducedin this paper.Moreover,the Thermo-Parameter Method(TPM)and Parametric Variational principles(PVP)are used to reduce the computational effort while maintaining the accuracy of solutions.A better solution isalso obtained in this paper.展开更多
The natural element method （NEM） is a newly- developed numerical method based on Voronoi diagram and Delaunay triangulation of scattered points, which adopts natural neighbour interpolation to construct trial functi...The natural element method （NEM） is a newly- developed numerical method based on Voronoi diagram and Delaunay triangulation of scattered points, which adopts natural neighbour interpolation to construct trial functions in the framework of Galerkin method. Owing to its distinctive advantages, the NEM is used widely in many problems of computational mechanics. Utilizing the NEM, this paper deals with numerical limit analysis of structures made up of perfectly rigid-plastic material. According to kinematic the- orem of plastic limit analysis, a mathematical programming natural element formulation is established for determining the upper bound multiplier of plane problems, and a direct iteration algorithm is proposed accordingly to solve it. In this algorithm, the plastic incompressibility condition is handled by two different treatments, and the nonlinearity and nons- moothness of the goal function are overcome by distinguishing the rigid zones from the plastic zones at each iteration. The procedure implementation of iterative process is quite simple and effective because each iteration is equivalent to solving an associated elastic problem. The obtained limit load multiplier is proved to monotonically converge to the upper bound of true solution. Several benchmark examples are investigated to validate the significant performance of the NEM in the application field of limit analysis.展开更多
To simplify the stability analysis of frozen soil slope, a pseudo-coupled numerical approach is developed. In this approach, the coupled heat transfer and water flow in frozen soils are simulated first, and based on t...To simplify the stability analysis of frozen soil slope, a pseudo-coupled numerical approach is developed. In this approach, the coupled heat transfer and water flow in frozen soils are simulated first, and based on the computed thermal-hydro field, the stability of frozen soil slope is evaluated. Although the shear strength for frozen soil is very complicated and is usually represented by a nonlinear MC failure criterion, a simple linear MC yield criterion is utilized. In this method, the internal friction angle is expressed as a function of volumetric ice content and the cohesion is fitted as a simple bilinear expression of Tand volumetric water content. To assess slope stability, the limit analysis is employed in conjunction with the recently developed a-section search algorithm. A frozen soil slope example is used to examine the proposed pseudo-coupled numerical approach, and numerical studies validate its effectiveness. Based on numerical results, it is seen that slope stability may be remarkably influenced by warming air （or grotmd surface） temperature. With increasing ground surface temperature, slope stability indicated by FOS may reduce to 1.0, implying that wanning air temperature could be a trigger of frozen soil slope failure.展开更多
The load-bearing capacity of ductile composite structures comprised of periodic composites is studied by a combined micro/macromechanicai approach. Firstly, on the microscopic level, a representative volume element （...The load-bearing capacity of ductile composite structures comprised of periodic composites is studied by a combined micro/macromechanicai approach. Firstly, on the microscopic level, a representative volume element （RVE） is selected to reflect the microstructures of the composite materials and the constituents are assumed to be elastic perfectly-plastic. Based on the homogenization theory and the static limit theorem, an optimization formulation to directly calculate the macroscopic strength domain of the RVE is obtained. The finite element modeling of the static limit analysis is formulated as a nonlinear mathematical programming and solved by the sequential quadratic programming method, where the temperature parameter method is used to construct the self-stress field. Secondly, Hill＇s yield criterion is adopted to connect the micromechanicai and macromechanical analyses. And the limit loads of composite structures are worked out on the macroscopic scale. Finally, some examples and comparisons are shown.展开更多
Considering the fact that in some complex cases,plate anchors are buried in multi-layered geotechnical materials,the ultimate dynamic analysis was performed to investigate the uplift capacity and failure mechanism of ...Considering the fact that in some complex cases,plate anchors are buried in multi-layered geotechnical materials,the ultimate dynamic analysis was performed to investigate the uplift capacity and failure mechanism of shallow strips and circular plate anchors in multi-layered soils.The nonlinear strength criterion and non-associated flow rule of geotechnical materials were introduced to investigate the influence of nonuniformity on the pullout performance and failure mechanism of shallow plate anchors.The expressions of the detaching curves or surfaces were obtained to reflect the failure mechanism,which can be used to figure out the ultimate uplift capacity and failure range.The results are generally in agreement with the numerical simulations and previous research.The effects of various parameters on the ultimate uplift capacity and failure mechanism of plate anchors in multi-layered soils were investigated,and it is found that the ultimate uplift capacity and failure range of shallow anchors increase with the increase of initial cohesion and dilatancy coefficient,but decrease with the unit weight,axial tensile stress and nonlinear coefficient.展开更多
Seismic permanent displacement of the soil walls plays an important role in design of these structures. Due to the increase in growth of urban areas and the limitations in use of flat grounds, many structures are buil...Seismic permanent displacement of the soil walls plays an important role in design of these structures. Due to the increase in growth of urban areas and the limitations in use of flat grounds, many structures are built near slopes and retaining walls. During earthquakes, these structures can apply an additional surcharge on the wall. The intensity and location of the surcharge is of considerable importance on the seismic displacements of the soil wall. In this study, by using the limit analysis and upper bound theorem, seismic permanent displacement of the soil wall under surcharge has been analyzed. Thus, a formulation is presented for calculating the yield acceleration and seismic displacement for different surcharge conditions. The effect of seismic acceleration, surcharge intensity, its location and soil properties is investigated. A parameter called the ＂displacement coefficient＂ is proposed, and is a potential modification for Newmark’s sliding-block method.展开更多
The plastic load-bearing capacity of ductile composites such as metal matrix composites is studied with an insight into the microstructures. The macroscopic strength of a composite is obtained by combining the homogen...The plastic load-bearing capacity of ductile composites such as metal matrix composites is studied with an insight into the microstructures. The macroscopic strength of a composite is obtained by combining the homogenization theory with static limit analysis, where the temperature parameter method is used to construct the self-equilibrium stress field. An interface failure model is proposed to account for the effects of the interface on the failure of composites. The static limit analysis with the finite-element method is then formulated as a constrained nonlinear programming problem, which is solved by the Sequential Quadratic Programming （SQP） method. Finally, the macroscopic transverse strength of perforated materials, the macroscopic transverse and off-axis strength of fiber-reinforced composites are obtained through numerical calculation. The computational results are in good agreement with the experimental data.展开更多
The ultimate strength of reinforced concrete(RC) rectangular members subjected to combined bending,shear and torsion is obtained from the limit analysis proposed in the present paper. Based on a warped failure surface...The ultimate strength of reinforced concrete(RC) rectangular members subjected to combined bending,shear and torsion is obtained from the limit analysis proposed in the present paper. Based on a warped failure surface determined by external loads, and a reasonable assumed stress distribution balancing external loads but not violating the yield condition, the bending-shear-torsion interaction can be derived from equilibrium conditions.According to the definition of lower-bound theorem in limit analysis, the calculated ultimate loads will be carried safely by the structure. The present method is a simple approach to obtain carrying capacities for RC elements under complex external loads. After comparing with the test results, a good agreement has been observed. The present method can be extended to explain the failure mechanism of RC members subjected to axial loads, and it is possible to develop a simple unified theory of RC members for engineering.展开更多
The solution of a slope stability problem can be approached by its least upper-bound and maximum lower-bound with high accuracy. The limit equilibrium methods that employ vertical slices imply a lower bound of the fac...The solution of a slope stability problem can be approached by its least upper-bound and maximum lower-bound with high accuracy. The limit equilibrium methods that employ vertical slices imply a lower bound of the factor of safety. It has been successfully extended to the area of active earth pressure analysis that accounts for different input of locations of earth pressure applications. Those methods that employ slices with inclined interfaces give an upper-bound approach to the stability analysis. It enjoys a sound mechanical background and is able to provide accurate solutions of soil plasticity. It has been successfully extended to the area of bearing capacity analysis in which various empirical coefficients are no longer necessary. The 3D upper- and lower-bound methods under this framework have been made possible and show great potential for solving various engineering problems.展开更多
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.展开更多
The analytical solutions for predicting the exact shape of collapse mechanisms in shallow tunnels with arbitrary excavation profiles were obtained by virtue of the upper bound theorem of limit analysis and variation p...The analytical solutions for predicting the exact shape of collapse mechanisms in shallow tunnels with arbitrary excavation profiles were obtained by virtue of the upper bound theorem of limit analysis and variation principle according to Hoek-Brown failure criterion. The seepage force was included in the upper bound limit analysis, and it was computed from the gradient of excess pore pressure distribution. The seepage was regarded as a work rate of external force. The numerical results of roof collapse in square and circular tunnels with different rock parameters were derived and discussed, which proves to be valid in comparison with the previous work. The influences of different parameters on the shape of collapsing blocks were also discussed.展开更多
Stabilizing pile is a kind of earth shoring structure frequently used in slope engineering. When the piles have cantilever segments above the ground,laggings are usually installed to avoid collapse of soil between pil...Stabilizing pile is a kind of earth shoring structure frequently used in slope engineering. When the piles have cantilever segments above the ground,laggings are usually installed to avoid collapse of soil between piles. Evaluating the earth pressure acting on laggings is of great importance in design process.Since laggings are usually less stiff than piles,the lateral pressure on lagging is much closer to active earth pressure. In order to estimate the lateral earth pressure on lagging more accurately,first,a model test of cantilever stabilizing pile and lagging systems was carried out. Then,basing the experimental results a three-dimensional sliding wedge model was established. Last,the calculation process of the total active force on lagging is presented based on the kinematic approach of limit analysis. A comparison is made between the total active force on lagging calculated by the formula presented in this study and the force on a same-size rigid retaining wall obtained from Rankine's theory. It is found that the proposed method fits well with the experimental results.Parametric studies show that the total active force on lagging increases with the growth of the lagging height and the lagging clear span; while decreases asthe soil internal friction angle and soil cohesion increase.展开更多
文摘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.
基金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.
基金support provided by the National Key R&D Program of China(No.2017YFC1501304)the National Natural Science Foundation of China(Nos.42090054,41922055 and 41931295)the Fundamental Research Funds for the Central Universities,China University of Geosciences(Wuhan)(No.CUGGC09).
文摘Multi-layer slopes are widely found in clay residue receiving fields.A generalized horizontal slice method(GHSM)for assessing the stability of multi-layer slopes that considers the energy dissipation between adjacent horizontal slices is presented.In view of the upper-bound limit analysis theory,the energy equation is derived and the ultimate failure mode is generated by comparing the sliding surface passing through the slope toe(mode A)with that below(mode B).In addition,the influence of the number of slices on the stability coefficients in the GHSM is studied and the stable value is obtained.Compared to the original method(Chen’s method),the GHSM can acquire more precise results,which takes into account the energy dissipation in the inner sliding soil mass.Moreover,the GHSM,limit equilibrium method(LEM)and numerical simulation method(NSM)are applied to analyze the stability of a multi-layer slope with different slope angles and the results of the safety factor and failure mode are very close in each case.The ultimate failure modes are shown to be mode B when the slope angle is not more than 28°.It illustrates that the determination of the ultimate sliding surface requires comparison of multiple failure modes,not only mode A.
基金support from the project of CNOOC China Limited-Shenzhen (Grant No. Z2007SLSZ-034)the foundation project of the State Key Laboratory of Petroleum Resource and Prospecting (Grant No. PRPDX2008-08) is gratefully acknowledged
文摘Extended reach wells （ERWs）, especially horizontal extended reach well with a high HD （horizontal displacement） to TVD （true vertical depth） ratio, represent a frontier technology and challenge the drilling limitations. Oil and gas reservoir in beaches or lakes and offshore can be effectively exploited by using extended reach drilling （ERD） technology. This paper focuses on the difficult technological problems encountered during exploiting the Liuhua 11-1 oil field in the South China Sea, China. Emphasis is on investigating the key subjects including prediction and control of open hole limit extension in offshore ERD, prediction of casing wear and its prevention and torque reduction, φ244.5mm casing running with floating collars to control drag force, and steerable drilling modes. The basic concept of limit extension in ERD is presented and the prediction method for open hole limit extension is given in this paper. A set of advanced drilling mechanics and control technology has been established and its practical results are verified by field cases. All those efforts may be significant for further investigating and practicing ERD limit theory and control technology in the future.
基金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.
基金Project(51378514)supported by the National Natural Science Foundation of China
文摘A method of analyzing the stability of twin shallow tunnels was presented using both limit analysis with nonlinear failure criterion and reliability theory.In the condition of nonlinear failure criterion,the critical clear distancesof twin shallow tunnels were obtained by analyzing the change of surrounding pressure.A reliability model was established based on limit state equation,and the failure probability was solved by virtue of Monte Carlo method.Safety factor and corresponding clear distance of different safety levels were obtained by introducing a target reliability index.The scope of clear distance for different safety levels is described,which can be used as a supplement and improvement to the design codes of tunnels.
基金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.
基金The project supported by National Natural Science Foundation of China
文摘According to the lower-bound theorem of limit analysis the Rigid Finite Element Meth-od(RFEM)is applied to structural limit analysis and the linear programmings for limit analysis are deducedin this paper.Moreover,the Thermo-Parameter Method(TPM)and Parametric Variational principles(PVP)are used to reduce the computational effort while maintaining the accuracy of solutions.A better solution isalso obtained in this paper.
基金supported by the National Foundation for Excellent Doctoral Thesis of China (200025)the Program for New Century Excellent Talents in University (NCET-04-0075)the National Natural Science Foundation of China (19902007)
文摘The natural element method （NEM） is a newly- developed numerical method based on Voronoi diagram and Delaunay triangulation of scattered points, which adopts natural neighbour interpolation to construct trial functions in the framework of Galerkin method. Owing to its distinctive advantages, the NEM is used widely in many problems of computational mechanics. Utilizing the NEM, this paper deals with numerical limit analysis of structures made up of perfectly rigid-plastic material. According to kinematic the- orem of plastic limit analysis, a mathematical programming natural element formulation is established for determining the upper bound multiplier of plane problems, and a direct iteration algorithm is proposed accordingly to solve it. In this algorithm, the plastic incompressibility condition is handled by two different treatments, and the nonlinearity and nons- moothness of the goal function are overcome by distinguishing the rigid zones from the plastic zones at each iteration. The procedure implementation of iterative process is quite simple and effective because each iteration is equivalent to solving an associated elastic problem. The obtained limit load multiplier is proved to monotonically converge to the upper bound of true solution. Several benchmark examples are investigated to validate the significant performance of the NEM in the application field of limit analysis.
基金supported in part by the Scientific Research Foundation for the 973 Program of China (No. 2012CB026104)Research Fund of Young Teachers for the Doctoral Program of Higher Education of China (No. 20110009120020)the Fundamental Research Funds of the Central Universities (No. 2013JBM059)
文摘To simplify the stability analysis of frozen soil slope, a pseudo-coupled numerical approach is developed. In this approach, the coupled heat transfer and water flow in frozen soils are simulated first, and based on the computed thermal-hydro field, the stability of frozen soil slope is evaluated. Although the shear strength for frozen soil is very complicated and is usually represented by a nonlinear MC failure criterion, a simple linear MC yield criterion is utilized. In this method, the internal friction angle is expressed as a function of volumetric ice content and the cohesion is fitted as a simple bilinear expression of Tand volumetric water content. To assess slope stability, the limit analysis is employed in conjunction with the recently developed a-section search algorithm. A frozen soil slope example is used to examine the proposed pseudo-coupled numerical approach, and numerical studies validate its effectiveness. Based on numerical results, it is seen that slope stability may be remarkably influenced by warming air （or grotmd surface） temperature. With increasing ground surface temperature, slope stability indicated by FOS may reduce to 1.0, implying that wanning air temperature could be a trigger of frozen soil slope failure.
基金supported by the National Natural Science Foundation of China (No.50809003)the National Foundation for Excellent Doctorial Dissertation of China (200025).
文摘The load-bearing capacity of ductile composite structures comprised of periodic composites is studied by a combined micro/macromechanicai approach. Firstly, on the microscopic level, a representative volume element （RVE） is selected to reflect the microstructures of the composite materials and the constituents are assumed to be elastic perfectly-plastic. Based on the homogenization theory and the static limit theorem, an optimization formulation to directly calculate the macroscopic strength domain of the RVE is obtained. The finite element modeling of the static limit analysis is formulated as a nonlinear mathematical programming and solved by the sequential quadratic programming method, where the temperature parameter method is used to construct the self-stress field. Secondly, Hill＇s yield criterion is adopted to connect the micromechanicai and macromechanical analyses. And the limit loads of composite structures are worked out on the macroscopic scale. Finally, some examples and comparisons are shown.
基金Project(51874202) supported by the National Natural Science Foundation of ChinaProject(2017JQ0003) supported by the Sichuan Youth Fund,China。
文摘Considering the fact that in some complex cases,plate anchors are buried in multi-layered geotechnical materials,the ultimate dynamic analysis was performed to investigate the uplift capacity and failure mechanism of shallow strips and circular plate anchors in multi-layered soils.The nonlinear strength criterion and non-associated flow rule of geotechnical materials were introduced to investigate the influence of nonuniformity on the pullout performance and failure mechanism of shallow plate anchors.The expressions of the detaching curves or surfaces were obtained to reflect the failure mechanism,which can be used to figure out the ultimate uplift capacity and failure range.The results are generally in agreement with the numerical simulations and previous research.The effects of various parameters on the ultimate uplift capacity and failure mechanism of plate anchors in multi-layered soils were investigated,and it is found that the ultimate uplift capacity and failure range of shallow anchors increase with the increase of initial cohesion and dilatancy coefficient,but decrease with the unit weight,axial tensile stress and nonlinear coefficient.
文摘Seismic permanent displacement of the soil walls plays an important role in design of these structures. Due to the increase in growth of urban areas and the limitations in use of flat grounds, many structures are built near slopes and retaining walls. During earthquakes, these structures can apply an additional surcharge on the wall. The intensity and location of the surcharge is of considerable importance on the seismic displacements of the soil wall. In this study, by using the limit analysis and upper bound theorem, seismic permanent displacement of the soil wall under surcharge has been analyzed. Thus, a formulation is presented for calculating the yield acceleration and seismic displacement for different surcharge conditions. The effect of seismic acceleration, surcharge intensity, its location and soil properties is investigated. A parameter called the ＂displacement coefficient＂ is proposed, and is a potential modification for Newmark’s sliding-block method.
基金Project supported by the Key Grant Project of Chinese Ministry of Education (No.0306)the National Foundationfor Excellent Doctoral Dissertation of China (No.200025).
文摘The plastic load-bearing capacity of ductile composites such as metal matrix composites is studied with an insight into the microstructures. The macroscopic strength of a composite is obtained by combining the homogenization theory with static limit analysis, where the temperature parameter method is used to construct the self-equilibrium stress field. An interface failure model is proposed to account for the effects of the interface on the failure of composites. The static limit analysis with the finite-element method is then formulated as a constrained nonlinear programming problem, which is solved by the Sequential Quadratic Programming （SQP） method. Finally, the macroscopic transverse strength of perforated materials, the macroscopic transverse and off-axis strength of fiber-reinforced composites are obtained through numerical calculation. The computational results are in good agreement with the experimental data.
基金the National Natural Science Foundation of China(No.51178265)
文摘The ultimate strength of reinforced concrete(RC) rectangular members subjected to combined bending,shear and torsion is obtained from the limit analysis proposed in the present paper. Based on a warped failure surface determined by external loads, and a reasonable assumed stress distribution balancing external loads but not violating the yield condition, the bending-shear-torsion interaction can be derived from equilibrium conditions.According to the definition of lower-bound theorem in limit analysis, the calculated ultimate loads will be carried safely by the structure. The present method is a simple approach to obtain carrying capacities for RC elements under complex external loads. After comparing with the test results, a good agreement has been observed. The present method can be extended to explain the failure mechanism of RC members subjected to axial loads, and it is possible to develop a simple unified theory of RC members for engineering.
基金Project (Nos. 50539100,50679035 and 50509027) supported by the National Natural ScienceFoundation of China
文摘The solution of a slope stability problem can be approached by its least upper-bound and maximum lower-bound with high accuracy. The limit equilibrium methods that employ vertical slices imply a lower bound of the factor of safety. It has been successfully extended to the area of active earth pressure analysis that accounts for different input of locations of earth pressure applications. Those methods that employ slices with inclined interfaces give an upper-bound approach to the stability analysis. It enjoys a sound mechanical background and is able to provide accurate solutions of soil plasticity. It has been successfully extended to the area of bearing capacity analysis in which various empirical coefficients are no longer necessary. The 3D upper- and lower-bound methods under this framework have been made possible and show great potential for solving various engineering problems.
文摘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.
基金Project(2013CB036004)supported by the National Basic Research Program of ChinaProjects(51178468,51378510)supported by National Natural Science Foundation of China
文摘The analytical solutions for predicting the exact shape of collapse mechanisms in shallow tunnels with arbitrary excavation profiles were obtained by virtue of the upper bound theorem of limit analysis and variation principle according to Hoek-Brown failure criterion. The seepage force was included in the upper bound limit analysis, and it was computed from the gradient of excess pore pressure distribution. The seepage was regarded as a work rate of external force. The numerical results of roof collapse in square and circular tunnels with different rock parameters were derived and discussed, which proves to be valid in comparison with the previous work. The influences of different parameters on the shape of collapsing blocks were also discussed.
基金financially supported by the National Key Technology Research and Development Program of the Ministry of Science and Technology of China under Grant No. 2012BAJ22B06
文摘Stabilizing pile is a kind of earth shoring structure frequently used in slope engineering. When the piles have cantilever segments above the ground,laggings are usually installed to avoid collapse of soil between piles. Evaluating the earth pressure acting on laggings is of great importance in design process.Since laggings are usually less stiff than piles,the lateral pressure on lagging is much closer to active earth pressure. In order to estimate the lateral earth pressure on lagging more accurately,first,a model test of cantilever stabilizing pile and lagging systems was carried out. Then,basing the experimental results a three-dimensional sliding wedge model was established. Last,the calculation process of the total active force on lagging is presented based on the kinematic approach of limit analysis. A comparison is made between the total active force on lagging calculated by the formula presented in this study and the force on a same-size rigid retaining wall obtained from Rankine's theory. It is found that the proposed method fits well with the experimental results.Parametric studies show that the total active force on lagging increases with the growth of the lagging height and the lagging clear span; while decreases asthe soil internal friction angle and soil cohesion increase.
