3D finite element method （FEM） simulation of groundwater flow during backward erosion piping

Backward erosion piping is an important failure mechanism for cohesive water retaining structures which are founded on a sandy aquifer. At present, the prediction models for safety assessment are often based on 2D assumptions. In this work, a 3D numerical approach of the groundwater flow leading to the erosion mechanism of backward erosion piping is presented and discussed. Comparison of the 2D and 3D numerical results explicitly demonstrates the inherent 3D nature of the piping phenomenon. In addition, the influence of the seepage length is investigated and discussed for both piping initiation and piping progression. The results clearly indicate the superiority of the presented 3D numerical model compared to the established 2D approach. Moreover, the 3D numerical results enable a better understanding of the complex physical mechanism involved in backward erosion piping and thus can lead to a significant improvement in the safety assessment of water retaining structures.

Propagations of Rayleigh and Love waves in ZnO films/glass substrates analyzed by three-dimensional finite element method

Propagation characteristics of surface acoustic waves（SAWs） in ZnO films/glass substrates are theoretically investigated by the three-dimensional（3D） finite element method. At first, for（11ˉ20） ZnO films/glass substrates, the simulation results confirm that the Rayleigh waves along the [0001] direction and Love waves along the [1ˉ100] direction are successfully excited in the multilayered structures. Next, the crystal orientations of the ZnO films are rotated, and the influences of ZnO films with different crystal orientations on SAW characterizations, including the phase velocity, electromechanical coupling coefficient, and temperature coefficient of frequency, are investigated. The results show that at appropriate h/λ, Rayleigh wave has a maximum k^2 of 2.4% in（90°, 56.5°, 0°） ZnO film/glass substrate structure; Love wave has a maximum k^2 of 3.81% in（56°, 90°, 0°） ZnO film/glass substrate structure. Meantime, for Rayleigh wave and Love wave devices, zero temperature coefficient of frequency（TCF） can be achieved at appropriate ratio of film thickness to SAW wavelength. These results show that SAW devices with higher k^2 or lower TCF can be fabricated by flexibly selecting the crystal orientations of ZnO films on glass substrates.

STUDY ON THREE-DIMENSIONAL FINITE BODIES CONTAINING CRACKS USING THE FINITE ELEMENT METHOD OF LINES

The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good accuracy can be obtained with relatively coarse girds. In particular, application to the tension specimen shows very good agreement with the evaluation of stress intensity factors, which is better than the results of other methods. This implies a considerable potential for using this method in the 3D analysis of finite geometry solids and suggests a possible extension of this technique to nonlinear material behavior.

Finite Element Simulation of Flexible Roll Forming with Supplemented Material Data and the Experimental Verification

Flexible roll forming is a promising manufacturing method for the production of variable cross section products. Considering the large plastic strain in this forming process which is much larger than that of uniform deformation phase of uniaxial tensile test, the widely adopted method of simulating the forming processes with non-supplemented material data from uniaxial tensile test will certainly lead to large error. To reduce this error, the material data is supplemented based on three constitutive models. Then a finite element model of a six passes flexible roll forming process is established based on the supplemented material data and the original material data from the uniaxial tensile test. The flexible roll forming experiment of a B pillar reinforcing plate is carried out to verify the proposed method. Final cross section shapes of the experimental and the simulated results are compared. It is shown that the simulation calculated with supplemented material data based on Swift model agrees well with the experimental results, while the simulation based on original material data could not predict the actual deformation accurately. The results indicate that this material supplement method is reliable and indispensible, and the simulation model can well reflect the real metal forming process. Detailed analysis of the distribution and history of plastic strain at different positions are performed. A new material data supplement method is proposed to tackle the problem which is ignored in other roll forming simulations, and thus the forming process simulation accuracy can be greatly improved.

3D Eddy Current and Temperature Field Analysis of Large Hydro-generators in Leading Phase Operations

As a common practice,a large hydro-generator will operate in leading phase conditions to absorb the reactive power of the power grid.However,the accurate and precise prediction of the leading phase operation capacity of a large hydro-generator has always been a formidable challenge to engineers and academicians because it is extremely hard to compute the eddy currents and losses as well as the local overheating in the pressure plate and finger.To address this problem,a full three dimensional(3D)finite element model and method of the coupled eddy current and temperature fields in the end region of a large hydro-generator are developed.The equivalent medium parameters used in the computations are comprehensively discussed.Moreover,some numerically based solution methodologies for accurate computation of the field and armature currents under different leading phase conditions are proposed.Numerical results on the coupled eddy current and temperature fields in the end regions of a 250 MW hydro-generator confirm positively the feasibility of the present work.

Lateral bearing characteristics of subsea wellhead assembly in the hydrate trial production engineering

Conductor and suction anchor are the key equipment providing bearing capacity in the field of deep-water drilling or offshore engineering,which have the advantages of high operation efficiency and short construction period.In order to drill a horizontal well in the shallow hydrate reservoir in the deep water,the suction anchor wellhead assembly is employed to undertake the main vertical bearing capacity in the second round of hydrate trial production project,so as to reduce the conductor running depth and heighten the kick-off point position.However,the deformation law of the deep-water suction anchor wellhead assembly under the moving load of the riser is not clear,and it is necessary to understand the lateral bearing characteristics to guide the design of its structural scheme.Based on 3D solid finite element method,the solid finite element model of the suction anchor wellhead assembly is established.In the model,the seabed soil is divided into seven layers,the contact between the wellhead assembly and the soil is simulated,and the vertical load and bending moment are applied to the wellhead node to simulate the riser movement when working in the deep water.The lateral bearing stability of conventional wellhead assembly and suction anchor wellhead assembly under the influence of wellhead load is discussed.The analysis results show that the bending moment is the main factor affecting the lateral deformation of the wellhead string;the anti-bending performance from increasing the outer conductor diameter is better than that from increasing the conductor wall thickness;for the subsea wellhead,the suction anchor obviously improves the lateral bearing capacity and reduces the lateral deformation.The conduct of the suction anchor wellhead assembly still needs to be lowered to a certain depth that below the maximum disturbed depth to ensure the lateral bearing stability,Thus,a method for the minimum conductor running depth for the suction anchor wellhead assembly is developed.The field implementations show that compared with the first round of hydrate trial production project,the conductor running depth is increased by 9.42 m,and there is no risk of wellhead overturning during the trial production.The method for determining the minimum conductor running depth in this paper is feasible and will still play an important role in the subsequent hydrate exploration and development.