Upper bound analysis of ultimate pullout capacity of shallow 3-D circular plate anchors based on nonlinear Mohr-Coulomb failure criterion

Based on the nonlinear Mohr-Coulomb failure criterion and the associated flow rules,the three-dimensional(3-D)axisymmetric failure mechanism of shallow horizontal circular plate anchors that are subjected to the ultimate pullout capacity(UPC)is determined.A derivative function of the projection function for projecting the 3-D axisymmetric failure surface on plane is deduced using the variation theory.By using difference principle,the primitive function of failure surface satisfying boundary condition and numerical solution to its corresponding ultimate pullout capacity function are obtained.The influences of nonlinear Mohr-Coulomb parameters on UPC and failure mechanism are studied.The result shows that UPC decreases with dimensionless parameter m and uniaxial tensile strength increases but increases when depth and radius of plate anchor,surface overload,initial cohesion,geomaterial density and friction angle increase.The failure surface is similar to a symmetrical spatial funnel,and its shape is mainly determined by dimensionless parameter m;the surface damage range expands with the increase of radius and depth of the plate anchor as well as initial cohesion but decreases with the increase of dimensionless parameter m and uniaxial tensile strength as well as geomaterial density.As the dimensionless parameter m=2.0,the numerical solution of UPC based on the difference principle is proved to be feasible and effective through the comparison with the exact solution.In addition,the comparison between solutions of UPC computed by variation method and those computed by upper bound method indicate that variation method outperforms upper bound method.

Three-dimensional upper bound limit analysis of underground cavities using nonlinear Baker failure criterion

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.

Upper bound seismic rotational stability analysis of gravity retaining walls considering embedment depth

Stability analysis of gravity retaining wall was currently based on the assumption that the wall had no embedment depth. The effect of earth berm was usually neglected. The present work highlighted the importance of embedment depth when assessing the seismic stability of gravity retaining walls with the pattern of pure rotation. In the framework of upper bound theorem of limit analysis, pseudo-static method was applied into two groups of parallel rigid soil slices methods in order to account for the effect of embedment depth on evaluating the critical acceleration of wall-soil system. The present analytical solution is identical to the results obtained from using limit equilibrium method, and the two methods are based on different theory backgrounds. Parameter analysis indicates that the critical acceleration increases slowly when the ratio of the embedment depth to the total height of the wall is from 0 to 0.15 and increases drastically when the ratio exceeds 0.15.

Prediction of central bursting defects in rod extrusion process with upper bound analysis method

The prediction of central bursting defects in the rod extrusion process through conical dies using the upper bound analysisis investigated. A kinematically admissible velocity field, including the radial and angular velocity components, is proposed. A newcriterion is presented to predict the occurrence of the central bursting defects. Parameter bobt, which represents the risk probability ofcracking, is proposed. It is calculated using the shape of the boundary at the entrance by minimizing the total power dissipationduring the extrusion process. When bobt is equal to or greater than bcr, central bursting occurs. Furthermore, the quantitativerelationships between central bursting defects and process parameters （semi die angle, reduction in area and frictional factor） arestudied. The results show that the central bursting defects are affected primarily by the reduction in area and the friction factor. Thepresented criterion is verified by comparing with the FEM simulation data and the results of the published paper.

LOWER BOUND LIMIT ANALYSIS OF THREE-DIMENSIONAL ELASTOPLASTIC STRUCTURES BY BOUNDARY ELEMENT METHOD

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.

Limit variation analysis of shallow rectangular tunnels collapsing with double-layer rock mass based on a three-dimensional failure mechanism

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.

Upper-bound limit analysis based on the natural element method

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.

Upper bound limit analysis of roof collapse in shallow tunnels with arbitrary cross sections under condition of seepage force

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.

Stability analysis for natural slope by kinematical approach

The stability of natural slope was analyzed on the basis of limit analysis. The sliding model of a kind of natural slope was presented. A new kinematically admissible velocity field for the new sliding model was constructed. The stability factor formulation by the upper bound theorem leads to a classical nonlinear programming problem, when the external work rate and internal energy dissipation were solved, and the constraint condition of the programming problem was given. The upper bound optimization problem can be solved efficiently by applying a nonlinear SQP algorithm, and stability factor was obtained, which agrees well with previous achievements.

Seismic response evaluation of base-isolated reinforced concrete buildings under bidirectional excitation

This paper reports on an investigation of the seismic response of base-isolated reinforced concrete buildings, which considers various isolation system parameters under bidirectional near-fault and far-fault motions. Three-dimensional models of 4-, 8-, and 12-story base-isolated buildings with nonlinear effects in the isolation system and the superstructure are investigated, and nonlinear response history analysis is carried out. The bounding values of isolation system properties that incorporate the aging effect of isolators are also taken into account, as is the current state of practice in the design and analysis of base-isolated buildings. The response indicators of the buildings are studied for near-fault and far-fault motions weight-scaled to represent the design earthquake （DE） level and the risk-targeted maximum considered earthquake （MCER） level. Results of the nonlinear response history analyses indicate no structural damage under DE-level motions for near-fault and far-fault motions and for MCER-level far-fault motions, whereas minor structural damage is observed under MCER- level near-fault motions. Results of the base-isolated buildings are compared with their fixed-base counterparts. Significant reduction of the superstructure response of the 12-story base-isolated building compared to the fixed-base condition indicates that base isolation can be effectively used in taller buildings to enhance performance. Additionally, the applicability of a rigid superstructure to predict the isolator displacement demand is also investigated. It is found that the isolator displacements can be estimated accurately using a rigid body model for the superstructure for the buildings considered.

Study on Failure Mechanism and Bearing Capacity of Three-Dimensional Rectangular Footing Subjected to Combined Loading

This paper presents two kinematic failure mechanisms of threc-dimensional rectangular footing resting on homogeneous undrained clay foundation under uniaxial vertical loading and uniaxial moment loading. The failure mechanism under vertical loading comprises a plane strain Prandti-type mechanism over the central part of the longer side, and the size of the mechanism gradually reduces at the ends of the longer side and over the shorter side as the corner of rectangular footing is being approached where the direction of soil motion remains normal to each corresponding side respectively. The failure mechanism under moment loading comprises a plane strain scoop sliding mechanism over the central part of the longer side, and the radius of scoop sliding mechanism increases linearly at the ends of the longer side. On the basis of the kinematic failure mechanisms mentioned above, the vertical ultimate bearing capacity and the ultimate bearing capacity against moment or moment ultimate bearing capacity are obtained by use of upper bound limit analysis theory. At the same time, numerical analysis results, Skempton＇ s results and Salgado et al. ＇s results are compared with this upper bound solution. It shows that the presented failure mechanisms and plastic limit analysis predictions are validated. In order to investigate the behaviors of undrained clay foundation beneath the rectangular footing subjected to the combined loadings, numerical analysis is adopted by virtue of the general-purpose FEM software ABAQUS, where the clay is assumed to obey the Mohr-Coulomb yielding criterion. The failure envelope and the ultimate bearing capacity are achieved by the numerical analysis results with the varying aspect ratios from length L to breadth B of the rectangular footing. The failure mechanisms of rectangular footing which are subjected to the combined vertical loading V and horizontal loading H （Vertical loading V and moment loading M, and horizontal loading H and moment loading M respectively are observed in the finite element analysis. ） is explained by use of the upper bound plasticity limit analysis theory. Finally, the reason of eccentricity of failure envelope in H-M loading space is given in this study, which can not be explained by use of the traditional ＇ swipe test＇.

Construction of improved rigid blocks failure mechanism for ultimate bearing capacity calculation based on slip-line field theory

Based on the slip-line field theory, a two-dimensional slip failure mechanism with mesh-like rigid block system was constructed to analyze the ultimate bearing capacity problems of rough foundation within the framework of the upper bound limit analysis theorem. In the velocity discontinuities in transition area, the velocity changes in radial and tangent directions are allowed. The objective functions of the stability problems of geotechnical structures are obtained by equating the work rate of external force to internal dissipation along the velocity discontinuities, and then the objective functions are transformed as an upper-bound mathematic optimization model. The upper bound solutions for the objective functions are obtained by use of the nonlinear sequential quadratic programming and interior point method. From the numerical results and comparative analysis, it can be seen that the method presented in this work gives better calculation results than existing upper bound methods and can be used to establish the more accurate plastic collapse load for the ultimate bearing capacity of rough foundation.

A LOWER BOUND APPROACH TO THE YIELD LOCI OF POROUS MATERIALS

A lower bound approach is proposed for the first time to solve the macroscopic yield loci of porous materials.The results are then compared with Gurson's upper bound yield loci and those of the experiments.It is shown that the present analysis is much more in accordance with the experimen- tal results than the Gurson's.

Flocking for swarm systems with fixed topology in a changing environment

This paper is mainly devoted to the flocking of a class of swarm with fixed topology in a changing environment. Firstly, the controller for each agent is proposed by employing the error terms between the state of the agent and the average state of its neighbors. Secondly, a sufficient condition for the swarm to achieve flocking is presented under assumptions that the gradient of the environment is bounded and the initial position graph is connected. Thirdly, as the environment is a plane, it is further proved that the velocity of each agent finally converges to the velocity of the swarm center although not one agent knows where the center of the group is. Finally, numerical examples are included to illustrate the obtained results.

Comparative study of material point method and upper bound limit analysis in slope stability analysis

The upper bound limit analysis(UBLA)is one of the key research directions in geotechnical engineering and is widely used in engineering practice.UBLA assumes that the slip surface with the minimum factor of safety(FSmin)is the critical slip surface,and then applies it to slope stability analysis.However,the hypothesis of UBLA has not been systematically verified,which may be due to the fact that the traditional numerical method is difficult to simulate the large deformation.In this study,in order to systematically verify the assumption of UBLA,material point method(MPM),which is suitable to simulate the large deformation of continuous media,is used to simulate the whole process of the slope failure,including the large-scale transportation and deposition of soil mass after slope failure.And a series of comparative studies are conducted on the stability of cohesive slopes using UBLA and MPM.The proposed study indicated that the slope angle,internal friction angle and cohesion have a remarkable effect on the slip surface of the cohesive slope.Also,for stable slopes,the calculation results of the two are relatively close.However,for unstable slopes,the slider volume determined by the UBLA is much smaller than the slider volume determined by the MPM.In other words,for unstable slopes,the critical slip surface of UBLA is very different from the slip surface when the slope failure occurs,and when the UBLA is applied to the stability analysis of unstable slope,it will lead to extremely unfavorable results.

Moving scanning emitter tracking by a single observer using time of interception:Observability analysis and algorithm

The target motion analysis（TMA） for a moving scanning emitter with known fixed scan rate by a single observer using the time of interception（TOI） measurements only is investigated in this paper.By transforming the TOI of multiple scan cycles into the direction difference of arrival（DDOA） model,the observability analysis for the TMA problem is performed.Some necessary conditions for uniquely identifying the scanning emitter trajectory are obtained.This paper also proposes a weighted instrumental variable（WIV） estimator for the scanning emitter TMA,which does not require any initial solution guess and is closed-form and computationally attractive.More importantly,simulations show that the proposed algorithm can provide estimation mean square error close to the Cramer-Rao lower bound（CRLB） at moderate noise levels with significantly lower estimation bias than the conventional pseudo-linear least square（PLS） estimator.