流体动压滑动轴承-转子系统动力特性分析

研究了非线性径向轴承支承的转子系统的动力行为。引入变分约束原理修正了流体润滑的Reynolds方程的变分形式 ,在几乎不增加计算量的情况下 ,求解了具有Reynolds边界的流体润滑问题 ,使得非线性油膜力及其Jacobian矩阵同时计算完成并且具有协调一致的精度。运用Newton Raphson方法在求得转子平衡点的同时求得了轴承的动力学系数。将预估 校正机理和Newton Raphson方法相结合给出了流体动压滑动轴承 转子系统Hopf分岔点所对应线性失稳转速的计算方法。运用打靶法并结合Floquet理论计算分析了流体动压滑动轴承 转子系统的非线性不平衡周期响应及其稳定性。数值结果表明上述方法不但节约了计算量 ,而且具有很高的精度。

Study on the FEM design of reinforced earth retaining wall with geogrid

At present,limit equilibrium method is often adopted in the design of reinforced earth retaining wall. Geotechnical engineers home and abroad have done a lot of work to improve the traditional calculation methods in recent years,while there are lots of defects. This paper first identifies the location of failure surface and safety factor through the finite element program of PLAXIS and then analyses the influencing factors of the stability of reinforced earth retaining wall with geogrid. The authors adopt strength reduction FEM (finite element method)in the design and stability analysis of reinforced earth retaining wall and have achieved some satisfying results. Without any assumptions,the new design method can automatically judge the failure mode of reinforced earth retaining wall,consider the influence of axial tensile stiffness of the reinforcement stripe on the stability of retaining wall,identify reasonable distance and length of the reinforcement stripe,and choose suitable parameters of reinforcement stripe,including strength,stiffness and pseudo-friction coefficient which makes the design optimal. It is proved through the calculation examples that this method is more reasonable,reliable and economical in the design of reinforced earth retaining wall.

Stability and reinforcement analysis of rock slope based on elasto-plastic finite element method

The rigid body limit equilibrium method(RBLEM) and finite element method(FEM) are two widely used approaches for rock slope's stability analysis currently. RBLEM introduced plethoric assumptions; while traditional FEM relied on artificial factors when determining factor of safety(FOS) and sliding surfaces. Based on the definition of structure instability that an elasto-plastic structure is not stable if it is unable to satisfy simultaneously equilibrium condition, kinematical admissibility and constitutive equations under given external loads, deformation reinforcement theory(DRT) is developed. With this theory, plastic complementary energy(PCE) can be used to evaluate the overall stability of rock slope, and the unbalanced force beyond the yield surface could be the identification of local failure. Compared with traditional slope stability analysis approaches, the PCE norm curve to strength reduced factor is introduced and the unbalanced force is applied to the determination of key sliding surfaces and required reinforcement. Typical and important issues in rock slope stability are tested in TFINE(a three-dimensional nonlinear finite element program), which is further applied to several representatives of high rock slope's stability evaluation and reinforcement engineering practice in southwest of China.

Analysis and research of long-term stability about the gob site in Northern Shaanxi Province

Based on the analysis of the basic characteristics for the gob site in Northern Shaanxi Province and the room and pillar mining way about coal mine, the variety rule of the coal beds below the site was studied by the using of FEM during the process of coal mining. The statuses of the stresses and strains and the varieties of the plastic area were simulated in the whole rock and coal pillars. The characters of stresses and deformation of the gob area of the coalmine were analyzed and evaluated after the site built in weathering. Moreover, the long-term stability of the gob area was predicted. As a result, the deformation of the gob area under the site is not been finished, and there is the danger that the gob site will collapsing as a whole; therefore, relative measures must be taken.

Coupled Numerical Analysis of the Stability Behaviour of Unsaturated Soil Slopes Under Rainfall Conditions

The stability behaviour of unsaturated soil slopes under rainfall conditions is investigated via a parametric finite element analysis, which is a fully coupled flow and deformation approach linked to a dynamic programming technique for determining the minimum factor of safety as well as its corresponding critical slip surface based on the stress fields from the numerical computation. The effects of rainfall features, soil strength parameters and permeability properties on slope stability are studied. The analyses revealed that the soil matric suction decreased during rainfall, especially in slopes with high permeability and/or with high suction angles of unsaturated soils. The influence of rainfall conditions on such slopes is quite obvious, and soil suction drops rapidly, which leads to a consequent quick reduction in the factor of safety.

Application of material strength reserve method to the stability analysis of earth dam of feilaixia multipurpose project

The material strength reserve method is practical in the study of the stability and failure mechanism of earth dam by analysing the development of failure zone of different shear strength parameters of the earth mass of the dam. The stability in the concrete dam and ensemble architecture has got general application while analysing. In combination with Feilaixia Multipurpose Project, application of this method to earth dam stability analysis was studied by plane Finite Element Method(FEM) for the first time. Through plane FEM, we can get the failure mechanism of earth dam and appraise to the security, for operating and managing put forward some reference suggestions.

Stress and Displacement Distribution of Soft Clay Slope with 2D and 3D Elastoplastic Finite Element Method

Based on elastoplastic model, 2D and 3D finite element method （FEM） are used to calculate the stress and displacement distribution in the soft clay slope under gravity and uniform load at the slope top. Stability analyses indicate that 3D boundary effect varies with the stress level of the slope. When the slope is stable, end effect of 3D space is not remarkable. When the stability decreases, end effect occurs; when the slope is at limit state, end effect reaches maximum. The energy causing slope failure spreads preferentially along y-z section, and when the failure resistance capability reaches the limit state, the energy can extend along x-axis direction. The 3D effect of the slope under uniform load on the top is related to the ratio of load influence width to slope height, and the effect is remarkable with the decrease of the ratio.

ERROR IMPROVEMENT OF TEMPORAL DISCRETIZATION IN TDFEM FOR ELECTROMAGNETIC ANALYSIS

Integral method is employed in this paper to alleviate the error accumulation of differential equation discretization about time variant t in Time Domain Finite Element Method (TDFEM) for electromagnetic simulation. The error growth and the stability condition of the presented method and classical central difference scheme are analyzed. The electromagnetic responses of 2D lossless cavities are investigated with TDFEM; high accuracy is validated with numerical results presented.

Dynamic response of jetting and cementing bucket platform to wave and current loading

The jetting and cementing bucket platform （JCBP） is a new type offshore oil-drilling platform. This paper aims to establish an analysis method for calculating the dynamic response of this platform. Based on the theory of elastic half space, the dynamic stiffness and damping of the platform＇ s foundation were obtained and attached to the end of the platform＇ s main jackets as a boundary condition. Then using finite element method （FEM）, the dynamic response of the platform due to wave and current loading was calculated. Furthermore, the whole platform＇ s finite element model was established and the dynamic response of the platform was calculated. The numerical results demonstrate that the present method by the usage of elastic half space theory and FEM is simple and it is reliable and efficient to calculate dynamic behavior of the platform in response to wave and current loading.

Application of the strength reduction FEM to diversion tunnel plug stability

Firstly,the common design principles for diversion tunnel plug are generalized,and two kinds of numerical analysis methods are discussed.Then the strength reduction FEM is introduced in numerical model analysis and the design steps of the plug's length are illustrated.During the progress to determine the plug's length,the equivalent plastic strain on the potential slip surface is assumed as the flag to tell the failure against sliding,and the plug stability is overall estimated from the plastic zone range and connectivity.

A Finite Element Model for the Instability Analysis of Flexible Pipes Tensile Armor Wires

The instability of the tensile armor wire of flexible pipes is a failure mode associated with deep and ultra-deep water applications. Real compressive forces acting on the pipe are necessary to trigger this process. The loss of stability may be divided into two distinct processes, according to the main direction of the wire＇s displacement： radial or lateral instability. This study aims at proposing a numerical tool for predicting lateral and radial critical buckling loads for the tensile armor wires of flexible pipes. A simple finite element model, based on springs and beams elements, was developed in ABAQUS~ to deal with this problem in an efficient and reliable manner. A parametric study was conducted concerning the behavior of the critical load when the laying angle, the initial curvature and the total pipe length are varied. The results were consistent with previously published literature data and analytical expressions, proving its applicability to pipe engineering projects. It also has the advantage of approaching the problem three-dimensionally, which allows further modelling modifications, such as including friction effects.

Mechanical Stability Analysis of BPM Support at SSRF

This paper concentrates on investigating the mechanical stability of a Beam Position Monitor (BMP)support prototype of the Shanghai Synchrotron Radiation Facility(SSRF)project.Both finite element analysis and vibration measurements have been performed.Inconsistent results between the simulations and experiments motivate us to study three connections between the support and the ground:ground bolt(used in the initial design),part grout and full grout(proposed in the later research).After changing the connection, the first eigenfrequency is increased from 20.2Hz to 50.2Hz,and the ratio of the integrated RMS displacement (4—50Hz)is decreased from 4.36 to 1.23 in the lateral direction.The mechanical stability is improved greatly.

Finite element reliability analysis of slope stability

The method of nonlinear finite element reliability analysis (FERA) of slope stability using the technique of slip surface stress analysis (SSA) is studied. The limit state function that can consider the direction of slip surface is given, and the formula-tions of FERA based on incremental tangent stiffness method and modified Aitken accelerating algorithm are developed. The limited step length iteration method (LSLIM) is adopted to calculate the reliability index. The nonlinear FERA code using the SSA technique is developed and the main flow chart is illustrated. Numerical examples are used to demonstrate the efficiency and robustness of this method. It is found that the accelerating convergence algorithm proposed in this study proves to be very efficient for it can reduce the iteration number greatly, and LSLIM is also efficient for it can assure the convergence of the iteration of the reliability index.

Finite element analysis for wellbore stability of transversely isotropic rock with hydraulic-mechanical-damage coupling

The finite element analysis （FEA） technology by hydraulic-mechanical-damage （HMD） coupling is proposed in this paper for wellbore stability analysis of transversely isotropic rock, developed basing on the recently established FEA technology for iso- tropic rock. The finite element （FE） solutions of numerical wellbore model, damage tensor calculation and Pariseau strength criterion for transversely isotropic rock are developed for researching the wellbore failure characteristics and computing the collapse and fracture pressure of laminated rock as shale reservoirs. The classic Blot constitutive for rock as porous medium is introduced to establish a set of FE equations coupling with elastic solid deformation and seepage flow. To be in accord with the inclined wellbore situation, the coordinate transformation for global, wellbore, in-situ stress and transversely isotropic for- mation coordinate systems is established for describing the in-situ stress field and the results in laminated rock. To be in accord with the practical situation, a three-dimensional FIE model is developed, in which several other auxiliary technologies are com- prehensively utilized, e.g., the typical Weibull distribution function for heterogeneous material description and adaptive tech- nology for mesh refinement. The damage tensor calculation technology for transversely isotropic rock are realized from the well-developed continuum damage variable of isotropic rock. The rock is subsequently developed into a novel conceptual and practical model considering the stress and permeability with the damage. The proposed method utilizing Parisean strength cri- terion fully reflects the strength parameters parallel or perpendicular to bedding of the transversely isotropic rock. To this end, an effective and reliable numerically three-step FEA strategy is well established. Numerical examples are given to show that the proposed method can establish efficient and applicable FE model and be suitable for analyzing the state of pore pressure and stress surrounding wellbore, furthermore to demonstrate the effectiveness and reliability of the instability analysis of wellbore failure region and the safe mud weight computation for collapse and fracture pressure of transversely isotropic rock.

Finite element analysis for inclined wellbore stability of transversely iso-tropic rock with HMCD coupling based on weak plane strength criterion

The finite element analysis(FEA) technology by hydraulic-mechanical-chemical-damage(HMCD) coupling is proposed in this paper for inclined wellbore stability analysis of water-sensitive and laminated rock, developed basing on the recently established FEA technology for transversely isotropic rock with hydraulic-mechanical-damage(HMD) coupling. The chemical activity of the drilling fluid is considered as phenomenological hydration behavior, the moisture content and parameters of rock considering hydration could be determined with time. The finite element(FE) solutions of numerical wellbore model considering the chemical activity of drilling fluid, damage tensor calculation and weak plane strength criterion for transversely isotropic rock are developed for researching the wellbore failure characteristics and computing the time-dependent collapse and fracture pressure of laminated rock as shale reservoirs. A three-dimensional FE model and elastic solid deformation and seepage flow coupled equations are developed, and the damage tensor calculation technology for transversely isotropic rock are realized by introducing effect of the hydration and the stress state under the current load. The proposed method utilizing weak plane strength criterion fully reflects the strength parameters in rock matrix and weak plane. To the end, an effective and reliable numerically three-step FEA strategy is well established for wellbore stability analysis. Numerical examples are given to show that the proposed method can establish efficient and applicable FE model and be suitable for analyzing the timedependsolutions of pore pressure and stresses, and the evolution region considering the hydration surrounding wellbore,furthermore to compute the collapse cycling time and the safe mud weight for collapse and fracture pressure of transversely isotropic rock.

Stability analysis of the temperature cracks in Xiaowan arch dam

In the Xiaowan arch dam there are massive temperature cracks nearly parallel to the dam axis. Obviously, whether the cracks may spread or not during the water storage process is one of the crucial factors for the safety of a dam. In this paper, a new type of crack element, in which the contact component is implicitly included into the concrete component, is proposed to simulate the effects of the existing cracks. The crack element is proved by numerical example to share the merits of both conventional contact elements and joint elements. With a finite element model of the cracked arch dam together with its rock foundation established, the transient displacement and stress fields of the dam are obtained. The complicated rock foundation, the construction process of the arch dam, the massive cracks, the transient temperature field, as well as the water storage process have been taken into consideration in the simulation. In addition to the global model, several sub-models for typical crack tips are also generated with finer elements placed around the tips. Thus, more accurate displacement and stress distribution are obtained by simultaneous sub-model simulation. Based on the calculation of stress intensity factor for crack tips by extension method, the temperature cracks in the Xiaowan arch dam are finally proved to be stable.

Stability analysis and a priori error estimate of explicit Runge-Kutta discontinuous Galerkin methods for correlated random walk with density-dependent turning rates

In this paper,we analyze the explicit Runge-Kutta discontinuous Galerkin(RKDG)methods for the semilinear hyperbolic system of a correlated random walk model describing movement of animals and cells in biology.The RKDG methods use a third order explicit total-variation-diminishing Runge-Kutta(TVDRK3)time discretization and upwinding numerical fluxes.By using the energy method,under a standard CourantFriedrichs-Lewy(CFL)condition,we obtain L2stability for general solutions and a priori error estimates when the solutions are smooth enough.The theoretical results are proved for piecewise polynomials with any degree k 1.Finally,since the solutions to this system are non-negative,we discuss a positivity-preserving limiter to preserve positivity without compromising accuracy.Numerical results are provided to demonstrate these RKDG methods.