复合型裂纹张开位移和裂尖附近塑性区域的有限元分析

用ＣＡＤ技术和Ｃ语言知识，对受双向载荷含Ⅰ－Ⅱ复合型裂纹十字形板进行了计算机自动划分有限单元体网格并生成数据文件。采用四边形八节点等参元对含Ⅰ－Ⅱ复合型裂纹十字形板进行了弹塑性有限元分析，获得了裂尖附近塑性区形状。弹塑性有限元法计算的裂纹面张开位移与试验结果相吻合。

Prediction of the thermo-elastic properties of knitted structural composites using FEM

It is a very important and complex task to estimate the thermo-elasticproperties of a textile structural composite. In this paper, the finite element method (FEM) wasused for the prediction of the orthotropic thermo-elastic properties of a composite reinforced byglass fiber knitted fabric. In order to define the final 3-D configuration of the loop reinforcingstructure, the interactions between the adjacent loops, the large displacement and the contactelements without friction were considered. The values predicted were compared with the experimentalresults.

Time-dependent lateral response of pile embedded in elasto-plastic soil

A two-dimensional （2D） finite element analysis was carried out to assess the time-dependent behavior of single vertical pile embedded in elasto-plastic soil. The finite element analyses were carried out using the linear elastic model for the structure of the pile, while the Mohr-Coulomb model was used for representing the soil behavior surrounding the pile. The study includes cohesionless and cohesive soil to assess the lateral response of pile in the two types of soil. The whole geoteehnical model is suitable for problem of piles to determine the design quantities such as lateral deformation, lateral soil stress and its variation with time. The model is verified based on the results of published cases and there is good comparison between the results of published ease and the present simulation model. It is found that, the pile in cohesionless soil has more resistance in the rapid loading and less one in the long term loading. On the other hand, the pile in cohesive soil shows opposite behavior.

Stability analysis of a concrete gravity dam and its foundation

The stability of dams and their foundations is an important problem to which dam engineers have paid close attention over the years. This paper presented two methods to analyze the stability of a gravity dam and its foundation. The direct analysis method was based on a rigid limit equilibrium method which regarded both dam and the rock foundation as undeformable rigid bodies. In this method, the safety factor of potential sliding surfaces was computed directly. The second method, the indirect analysis method, was based on elasto-plastic theory and employs nonlinear finite element method (FEM) in the analysis of stresses and deformation in the dam and its foundation. The determination of the safety degree of the structure was based on the convergence and abrupt the change criterion. The results obtained showed that structures' constituent material behavior played an active role in the failure of engineered structures in addition to the imposed load.

A new numerical method for determining collapse load-carrying capacity of structure made of elasto-plastic material

Determination of collapse load-carrying capacity of elasto-plastic material is very important in designing structure. The problem is commonly solved by elasto-plastic finite element method （FEM）. In order to deal with material nonlinear problem involving strain softening problem effectively, a new numerical method-damped Newton method was proposed. The iterative schemes are discussed in detail for pure equilibrium models. In the equilibrium model, the plasticity criterion and the compatibility of the strains are verified, and the strain increment and plastic factor are treated as independent unknowns. To avoid the stiffness matrix being singularity or condition of matrix being ill, a damping factor a was introduced to adjust the value of plastic consistent parameter automatically during the iterations. According to the algorithm, the nonlinear finite element program was complied and its numerical example was calculated. The numerical results indicate that this method converges very fast for both small load steps and large load steps. Compared with those results obtained by analysis and experiment, the predicted ultimate bearing capacity from the proposed method is identical.

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.

A Comparison of Different Approaches in Numerical Modeling of Pavement Granular Material

Modeling pavement granular materials have played through an experimental or numerical approach to predict the a significant role in pavement design procedure. Modeling can be granular behavior during cyclic loading. Current design process in western Australia is based on linear elastic analysis of layers. The analysis is largely performed through a well-known program CIRCLY which is applied to model bound pavement material behavior. The KENLAYER is one of the common pavement software models used for pavement design in the United State which performs non-linear analysis for granular materials. Alternatively, a general finite element program such as ABAQUS can be used to model the complicated behavior ofmultilayer granular materials. This study is to compare results of numerical modeling with these three programs on a sample constructed pavement model. Moreover, a parametric study on the effects of Poisson ratio over the surface deflection of the flexible pavement has been conducted. It is found that increase in Poisson ratio of asphalt layer will increase the surface deflection while the increase in Poisson ratio of granular layers decreases the surface deflection.

Dynamic Analysis of Kineto-Elastic Beam System with Second-order Effect

Dynamic equations of motional flexible beam elements were derived considering second-order effect. Non-linear finite element method and three-node Euler-Bernoulli beam elements were used. Because accuracy is higher in non-linear structural analysis,three-node beam elements are used to deduce shape functions and stiffness matrices in dynamic equations of flexible elements. Static condensation method was used to obtain the finial dynamic equations of three-node beam elements. According to geometrical relations of nodal displacements in concomitant and global coordinate system,dynamic equations of elements can be transformed to global coordinate system by concomitant coordinate method in order to build the global dynamic equations. Analyzed amplitude condition of flexible arm support of a port crane,the results show that second-order effect should be considered in kinetic-elastic analysis for heavy load machinery of big flexibility.

Numerical simulation of the separation between concrete face slabs and cushion layer of Zipingpu dam during the Wenchuan earthquake

The Zipingpu concrete-faced rockfill dam(CFRD)experienced strong ground motion from the 2008 Wenchuan earthquake.Separation between concrete face slabs and the cushion layer was observed after the earthquake.The separation voids under the stage III slabs make up 55%of the total area of the stage III slabs.The observed maximum height of the separation voids was nearly 23 cm at the top of the stage III slabs.Separation voids were also observed locally below the top of stage II slabs near the left abutment,with a maximum height of 7 cm.In this study,a static and dynamic elasto-plastic finite element analysis on Zipingpu CFRD was conducted to capture the separation during the Wenchuan earthquake.The rockfill materials were described using a state-dependent elasto-plastic model that considered particle breakage.The model parameters of rockfill materials were obtained from feedback analysis.The numerical results were largely consistent with the field measurements during construction and after the Wenchuan earthquake.A three-dimensional state-dependent elasto-plastic model that can trace the separation and re-contact of a soil-structure interface was employed to investigate the interaction between concrete face slabs and a cushion layer.The analysis showed the distribution of separation voids observed in the Zipingpu CFRD has a close relationship to the water level and slab dislocations at the time of the earthquake.The phenomenon of the separation from the Wenchuan earthquake was successfully captured by the proposed numerical procedure.

Base force element method (BFEM) on complementary energy principle for linear elasticity problem

Using the concept of base forces as state variables,a new finite element method-the base force element method (BFEM) on complementary energy principle for linear elasticity problems is presented.Firstly,an explicit expression of compliance matrix for an element is derived through base forces by dyadic vectors.Then,the explicit control equations of finite element method of complementary energy principle are derived using Lagrange multiplier method.Thereafter,the base forces element procedure for linear elasticity is developed.Finally,several examples are analyzed to illustrate the reliability and accuracy of the formulation and the procedure.

Numerical and experimental analysis of the reusability of spring energized metal C rings

It is well known that static seals are usually designed to be used once. In this paper, we discuss the reusability of a type of static seal called the spring energized metal C ring, which is required to be reused 4-5 times without maintenance during the lifetime of the pump. A theoretical analysis based on the nonlinear finite element method and an experimental investigation of the reusability of the C ring are carried out. The effects of the elastic modulus of the ring material, the ring wall thickness, the amount of assembly interference between the spring and the C ring, the operating temperature and the compression ratio are discussed. The parameter S, which denotes the reaction force difference percentage, is determined to measure the reusability of the C ring. The results show that the ring wall thickness plays an important role in the reusability. Although the normal force of a thick wall decreases after the ring is reused eight times, the value is still higher than the first normal force of a thin wall. In contrast, the elastic modulus of the material has the lowest impact. Here, a pre-compression ratio of 18.2% not only produces a larger first normal force but also leads to a smaller reduction in the normal force when the ring is reused. A simulated experiment is carried out to investigate the influence of the coating of the C ring on its reusability. The results show that when plastic deformation occurs on the contact surface, the contact state of the C ring coating becomes a major factor that affects the seal performance. Aging treatment at high temperatures （250℃） reduces gized metal C ring, which leads to the deterioration in the reusability the binding force of the silver coating on the spring ener of the ring.

Submicron size-scale mapping of carbonate effective elastic properties from FIB-SEM images and finite element method

In this paper, an automatic unstructured focused ion beam （FIB） and scanning electron microscopy （SEM） images induced representative volume element （RVE） finite element （FE） method is developed to predict submicron scale carbonate rock effective Young＇s and bulk moduli and Poisson＇s ratio on parallel CPU-GPU platform. Based on high resolution-contrast surface morphology and internal fabric-texture structure images from carbonate rock specimen （covered 0.12-64 μm2 area and 8000 μm3 domain）, the cubic RVE FE models are constructed from different sites through Avizo with user-defined parameters Matlab coding. The effective Young＇s and bulk moduli and Poisson＇s ratio of the different RVEs and porosity and pore size are computed by using periodic boundary condition in the well-known FE software Abaqus. FE mesh sensitivity analysis has been conducted where all moduli converge to a certain constant value at larger FE mesh density. The effect of fabric-texture （pore size, shape, and distribution） on the elastic properties is discussed. The correlations between the computed effective elastic properties and pore size, porosity, RVE size have been established. The simulation results show that the effective Young＇s and bulk moduli and Poisson＇s ratio have strong anisotropic behavior and depend on RVE size, porosity and pore size. The RVE size, porosity and pore size are three independent factors in affecting of the effective elastic moduli, the effect mechanism of porosity and pore size is same while the effect mechanism of RVE size is difference.

HIGH ACCURACY ANALYSIS OF THE FINITE ELEMENT METHOD FOR NONLINEAR VISCOELASTIC WAVE EQUATIONS WITH NONLINEAR BOUNDARY CONDITIONS

The standard finite elements of degree p over the rectangular meshes are applied to solve a kind of nonlinear viscoelastic wave equations with nonlinear boundary conditions, and the superclose property of the continuous Galerkin approximation is derived without using the nonclassical elliptic projection of the exact solution of the model problem. The global superconvergence of one order higher than the traditional error estimate is also obtained through the postprocessing technique.

A finite element approach for seismic design of arch dam-abutment structures

We present a 3-D finite element (FE) approach to find the optimal distribution of seismic reinforcement force to secure high arch dam-abutment structures against certain earthquake actions. Nonlinear FE time history analysis is performed on the structure to find the seismic responses, using the associated elastic-perfectly plastic material description. The concept of plastic complementary energy is introduced to structural dynamics to quantify the structure's resistance against the seismic action throughout the time history and to indicate the critical moments when extreme extents of dynamic failure occur. Meanwhile the distributions of the unbalanced force at these critical moments reveal the dominant patterns of the dynamic failure. By the principle of minimum plastic complementary energy, the unbalanced force is just the counterforce of optimal reinforcement force to secure the self-unsupportable structure against the earthquake, which makes the seismic design more targeted and effective. Seismic design analysis is performed on Maji high arch dam-abutment structure. The results could to a large extent guide the seismic design, showing that several structural surfaces lying at the upper abutment are the most seismically vulnerable. This application indicates good applicability of this approach to large-scale projects.