Three-Dimensional Finite Element Numerical Simulation and Physical Experiment for Magnetism-Stress Detecting in Oil Casing

The casing damage has been a big problem in oilfield production. The current detection methods mostly are used after casing damage, which is not very effective. With the rapid development of China's offshore oil industry, the number of offshore oil wells is becoming larger and larger. Because the cost of offshore oil well is very high, the casing damage will cause huge economic losses. What's more, it can also bring serious pollution to marine environment. So the effective methods of detecting casing damage are required badly. The accumulation of stress is the main reason for the casing damage. Magnetic anisotropy technique based on counter magnetostriction effect can detect the stress of casing in real time and help us to find out the hidden dangers in time. It is essential for us to prevent the casing damage from occurring. However, such technique is still in the development stage. Previous studies mostly got the relationship between stress and magnetic signals by physical experiment, and the study of physical mechanism in relative magnetic permeability connecting the stress and magnetic signals is rarely reported. The present paper uses the ANSYS to do the three-dimensional finite element numerical simulation to study how the relative magnetic permeability works for the oil casing model. We find that the quantitative relationship between the stress' s variation and magnetic induction intensity's variation is: Δδ =K* ΔB, K = 8.04×109, which is proved correct by physical experiment.

Numerical analysis of high‑speed railway slab tracks using calibrated and validated 3D time‑domain modelling

Concrete slabs are widely used in modern railways to increase the inherent resilient quality of the tracks,provide safe and smooth rides,and reduce the maintenance frequency.In this paper,the elastic performance of a novel slab trackform for high-speed railways is investigated using three-dimensional finite element modelling in Abaqus.It is then compared to the performance of a ballasted track.First,slab and ballasted track models are developed to replicate the full-scale testing of track sections.Once the models are calibrated with the experimental results,the novel slab model is developed and compared against the calibrated slab track results.The slab and ballasted track models are then extended to create linear dynamic models,considering the track geodynamics,and simulating train passages at various speeds,for which the Ledsgard documented case was used to validate the models.Trains travelling at low and high speeds are analysed to investigate the track deflections and the wave propagation in the soil,considering the issues associated with critical speeds.Various train loading methods are discussed,and the most practical approach is retained and described.Moreover,correlations are made between the geotechnical parameters of modern high-speed rail and conventional standards.It is found that considering the same ground condition,the slab track deflections are considerably smaller than those of the ballasted track at high speeds,while they show similar behaviour at low speeds.

Three-dimensional forward modeling for magnetotelluric sounding by finite element method

A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.

Three-Dimensional Simulations of RESET Operation in Phase-Change Random Access Memory with Blade-Type Like Phase Change Layer by Finite Element Modeling

An optimized device structure for reducing the RESET current of phase-change random access memory （PCRAM） with blade-type like （BTL） phase change layer is proposed. The electrical thermal analysis of the BTL cell and the blade heater contactor structure by three-dimensional finite element modeling are compared with each other during RESET operation. The simulation results show that the programming region of the phase change layer in the BTL cell is much smaller, and thermal electrical distributions of the BTL cell are more concentrated on the TiN/GST interface. The results indicate that the BTL cell has the superiorities of increasing the heating efficiency, decreasing the power consumption and reducing the RESET current from 0.67mA to 0.32mA. Therefore, the BTL cell will be appropriate for high performance PCRAM device with lower power consumption and lower RESET current.

Autologous nerve graft repair of different degrees of sciatic nerve defect:stress and displacement at the anastomosis in a three-dimensional finite element simulation model

In the repair of peripheral nerve injury using autologous or synthetic nerve grafting, the mag- nitude of tensile forces at the anastomosis affects its response to physiological stress and the ultimate success of the treatment. One-dimensional stretching is commonly used to measure changes in tensile stress and strain; however, the accuracy of this simple method is limited. There- fore, in the present study, we established three-dimensional finite element models of sciatic nerve defects repaired by autologous nerve grafts. Using PRO E 5.0 finite element simulation software, we calculated the maximum stress and displacement of an anastomosis under a 5 N load in 10-, 20-, 30-, 40-mm long autologous nerve grafts. We found that maximum displacement increased with graft length, consistent with specimen force. These findings indicate that three-dimensional finite element simulation is a feasible method for analyzing stress and displacement at the anas- tomosis after autologous nerve grafting.

Finite Element Analysis in Combination with Perfectly Matched Layer to the Numerical Modeling of Acoustic Devices in Piezoelectric Materials

The characterization of finite length Surface Acoustic Wave (SAW) and Bulk acoustic Wave (BAW) resonators is addressed here. The Finite Element Analysis (FEA) induces artificial wave reflections at the edges of the mesh. In fact, these ones do not contribute in practice to the corresponding experimental response. The Perfectly Matched Layer (PML) method, allows to suppress the boundary reflections. In this work, we first demonstrate the basis of PML adapted to FEA formalism. Next, the results of such a method are depicted allowing a discussion on the behavior of finite acoustic resonators.

FINITE ELEMENTAN ALYSIS FOR THE UNSTEADY NEARSHORECIR CULATION DUE TO WAVE－CURRENT INTER－ACTION（I）──NUMERICAL MODEL

In lhis paper, a numerical model for predicting the unsteady nearshore circulationdue to wave-current interaction was proposed. In addition io the traditional continuity,momentunm and energy. equations, the dispersion and refraction relations were includedin the governing equalions. Moreover, the effects of lateral shears, wind, radiation andbottom stresses were analysed in the governing equalions. Therefore, we expect thatthis model may more completely and exactly reflect the law of ware-currentinteraction. In part (II) we will adopt the selective lumping two-step explicit finite elementmethod to solve the model, and some examples will be presented.

Evaluation value of three-dimensional finite element model analysis for bone mineral density and bone metabolism activity in patients with osteoporosis

Objective: To study the evaluation value of three-dimensional finite element model analysis for bone mineral density (BMD) and bone metabolism activity in patients with osteoporosis. Methods: A total of 218 patients who were diagnosed with osteoporosis in the hospital between February 2014 and January 2017 were collected as observation group, and 100 healthy volunteers who received physical examination in the hospital during the same period were selected as normal control group. The femoral head of the two groups was analyzed by three-dimensional finite element model, and the femoral head BMD levels and serum bone metabolism index contents were measured. Pearson test was used to evaluate the evaluation value of femoral head three-dimensional finite element model for osteoporosis. Results: The cancellous bone and cortical bone Von Mises stress value of observation group were lower than those of normal control group, and femoral neck BMD value of observation group was lower than that of normal control group;serum bone metabolism index BGP content was lower than that of normal control group while NBAP, TRACP-5b and CTX-1 contents were higher than those of normal control group. Pearson test showed that the cancellous bone and cortical bone Von Mises stress value of patients with osteoporosis were directly correlated with BMD value and bone metabolism index contents. Conclusion: The three-dimensional finite element model analysis resultsof patients with osteoporosis can objectively reflect the femoral headBMD value and bone metabolism activity, and is a reliable way to evaluate the risk of long-term fractures.

Numerical analysis of rock fracturing by gas pressure using the extended finite element method

High energy gas fracturing is a simple approach of applying high pressure gas to stimulate wells by gen- erating several radial cracks without creating any other damages to the wells. In this paper, a numerical algorithm is proposed to quantitatively simulate propagation of these fractures around a pressurized hole as a quasi-static phenomenon. The gas flow through the cracks is assumed as a one-dimensional transient flow, governed by equations of conservation of mass and momentum. The fractured medium is modeled with the extended finite element method, and the stress intensity factor is calculated by the simple, though sufficiently accurate, displacement ex- trapolation method. To evaluate the proposed algorithm, two field tests are simulated and the unknown parameters are determined through calibration. Sensitivity analyses are performed on the main effective parameters. Considering that the level of uncertainty is very high in these types of engineering problems, the results show a good agreement with the experimental data. They are also consistent with the theory that the final crack length is mainly determined by the gas pressure rather than the initial crack length produced by the stress waves.

Three-dimensional finite element analysis on effects of tunnel construction on nearby pile foundation

A three-dimensional finite element simulation was carried out to investigate the effects of tunnel construction on nearby pile foundation.The displacement controlled model (DCM) was used to simulate the tunneling-induced volume loss effects.The numerical model was verified based on the results of a centrifuge test and a set of parametric studies was implemented based on this model.There is good agreement between the trend of the results of the centrifuge test and the present model.The results of parametric studies show that the tunnelling-induced pile internal force and deformation depend mainly on the pile?tunnel distance,the pile length to tunnel depth ratio and the volume loss.Two different zones are separated by a 45° line projected from the tunnel springline.Within the zone of influence,the pile is subjected to tensile force and large settlement;whereas outside the zone of influence,dragload and small settlement are induced.It is also established that the impact of tunnelling on a pile group is substantially smaller as compared with a single pile in the same location with the rear pile in a group,demonstrating a positive pile group effect.

Benchmarking of two three-dimensional numerical models in time/space domain to predict railway-induced ground vibrations

In the last 30 years,the scientific community has developed and proposed different models and numerical approaches for the study of vibrations induced by railway traffic.Most of them are formulated in the frequency/wave number domain and with a 2.5D approach.Three-dimensional numerical models formulated in the time/space domain are less frequently used,mainly due to their high computational cost.Notwithstanding,these models present very attractive characteristics,such as the possibility of considering nonlinear behaviors or the modelling of excess pore pressure and non-homogeneous and non-periodic geometries in the longitudinal direction of the track.In this study,two 3D numerical approaches formulated in the time/space domain are compared and experimentally validated.The first one consists of a finite element approach and the second one of a finite difference approach.The experimental validation in an actual case situated in Carregado(Portugal)shows an acceptable fitting between the numerical results and the actual measurements for both models.However,there are some differences among them.This study therefore includes some recommendations for their use in practical soil dynamics and geotechnical engineering.

Numerical Investigation of Thermal Behavior of CNC Machine Tool and Its Effects on Dimensional Accuracy of Machined Parts

The dimensional accuracy of machined parts is strongly influenced by the thermal behavior of machine tools (MT). Minimizing this influence represents a key objective for any modern manufacturing industry. Thermally induced positioning error compensation remains the most effective and practical method in this context. However, the efficiency of the compensation process depends on the quality of the model used to predict the thermal errors. The model should consistently reflect the relationships between temperature distribution in the MT structure and thermally induced positioning errors. A judicious choice of the number and location of temperature sensitive points to represent heat distribution is a key factor for robust thermal error modeling. Therefore, in this paper, the temperature sensitive points are selected following a structured thermomechanical analysis carried out to evaluate the effects of various temperature gradients on MT structure deformation intensity. The MT thermal behavior is first modeled using finite element method and validated by various experimentally measured temperature fields using temperature sensors and thermal imaging. MT Thermal behavior validation shows a maximum error of less than 10% when comparing the numerical estimations with the experimental results even under changing operation conditions. The numerical model is used through several series of simulations carried out using varied working condition to explore possible relationships between temperature distribution and thermal deformation characteristics to select the most appropriate temperature sensitive points that will be considered for building an empirical prediction model for thermal errors as function of MT thermal state. Validation tests achieved using an artificial neural network based simplified model confirmed the efficiency of the proposed temperature sensitive points allowing the prediction of the thermally induced errors with an accuracy greater than 90%.

Fracture of films caused by uniaxial tensions:a numerical model

Surface cracks are commonly observed in coatings and films.When structures with coatings are subject to stretching,opening mode cracks are likely to form on the surface,which may further lead to other forms of damage,such as interfacial delamination and substrate damage.Possible crack forms include cracks extending towards the interface and channeling across the film.In this paper,a two-dimensional numerical model is proposed to obtain the structural strain energy at arbitrary crack lengths for bilayer structures under uniaxial tension.The energy release rate and structural stress intensity factors can be obtained accordingly,and the effects of geometry and material features on fracture characteristics are investigated,with most crack patterns being confirmed as unstable.The proposed model can also facilitate the analysis of the stress distribution in periodic crack patterns of films.The results from the numerical model are compared with those obtained by the finite element method(FEM),and the accuracy of the theoretical results is demonstrated.

Kinematic source model for simulation of near-fault ground motion field using explicit finite element method

This paper briefly reviews the characteristics and major processes of the explicit finite element method in modeling the near-fault ground motion field. The emphasis is on the finite element-related problems in the finite fault source modeling. A modified kinematic source model is presented, in which vibration with some high frequency components is introduced into the traditional slip time function to ensure that the source and ground motion include sufficient high frequency components. The model presented is verified through a simple modeling example. It is shown that the predicted near-fault ground motion field exhibits similar characteristics to those observed in strong motion records, such as the hanging wall effect, vertical effect, fling step effect and velocity pulse effect, etc.

Finite Element Convergence for State-Based Peridynamic Fracture Models

We establish the a priori convergence rate for finite element approximations of a class of nonlocal nonlinear fracture models.We consider state-based peridynamic models where the force at a material point is due to both the strain between two points and the change in volume inside the domain of the nonlocal interaction.The pairwise interactions between points are mediated by a bond potential of multi-well type while multi-point interactions are associated with the volume change mediated by a hydrostatic strain potential.The hydrostatic potential can either be a quadratic function,delivering a linear force–strain relation,or a multi-well type that can be associated with the material degradation and cavitation.We first show the well-posedness of the peridynamic formulation and that peridynamic evolutions exist in the Sobolev space H2.We show that the finite element approximations converge to the H2 solutions uniformly as measured in the mean square norm.For linear continuous fi nite elements,the convergence rate is shown to be Ct Δt+Csh2/ε2,where𝜖is the size of the horizon,his the mesh size,and Δt is the size of the time step.The constants Ct and Cs are independent of Δt and h and may depend on ε through the norm of the exact solution.We demonstrate the stability of the semi-discrete approximation.The stability of the fully discrete approximation is shown for the linearized peridynamic force.We present numerical simulations with the dynamic crack propagation that support the theoretical convergence rate.

Numerical study of fatigue damage of asphalt concrete using cohesive zone model

In order to investigate the fatigue behavior of asphalt concrete, a new numerical approach based on a bi-linear cohesive zone model （CZM） is developed. Integrated with the CZM, a fatigue damage evolution model is established to indicate the gradual degradation of cohesive properties of asphalt concrete under cyclic loading. Then the model is implemented in the finite element software ABAQUS through a user-defined subroutine. Based on the proposed model, an indirect tensile fatigue test is finally simulated. The fatigue lives obtained through numerical analysis show good agreement with laboratory results. Fatigue damage accumulates in a nonlinear manner during the cyclic loading process and damage initiation phase is the major part of fatigue failure. As the stress ratio increases, the time of the steady damage growth stage decreases significantly. It is found that the proposed fatigue damage evolution model can serve as an accurate and efficient tool for the prediction of fatigue damage of asphalt concrete.

Numerical Modeling on Suspended Sediment Transportation in the Hangzhou Bay

In this paper the characteristics of tidal flow and seasonal variation of seidment content in the Hangzhou Bay and their affecting factors are studied. Field investigations and data analysis indicate that the sediment movement is mainly influenced by the Yangtze estuary and the sediment of the Yangtze estuary is induced by wind wave and tidal flow. Owing to the variation of dynamic conditions, the instantaneous sediment content is controlled by tidal flow, wind wave, depth of water and tidal range synthetically. A sediment content relationship formula is established with related factors. A non-equilibrium 2-dimensional numerical model of suspended sediment transportation is set up, and the finite element method is applied. The computation results of the model is in accordance with field data.

Stability assessment of landslide-prone road cut rock slopes in Himalayan terrain:A finite element method based approach

Large-scale slope destabilization could be aggravated due to swift urbanization and ever-rising demands of geoengineering projects such as dams,tunnels,bridges and widening roads.National Highway-58 connects Delhi to Badrinath in India,which passes through complex geomorphological and geological terrain and often encounters cut slopes susceptible to slope failures.In the present investigation,a detailed geotechnical appraisal is conducted along the road cut slopes from Rishikesh to Devprayag in the Himalayas.Twenty vulnerable road cut slopes were demarcated for detailed slope stability analysis using Phase2D finite element modeling simulator.Nonlinear generalized Hoek-Brown(GHB)criterion was adopted for stability analyses.Out of 20 slopes,five slopes(S6,S7,S18,S19 and S20)are unstable with factor of safety(FoS)less than or equal to 1,and thus needs immediate attention.The FoS values of four slopes(S2,S9,S13 and S17)lie between 1 and 1.3,i.e.marginally stable,and slopes S1,S3,S4,S5,S8,S10,Sll,S12,S14,S15 and S16 are stable.Mohr-Coulomb(MC)criterion was also adopted to compare the slope stability analysis with GHB criterion.The FoS calculated from GHB criterion is close to that using MC criterion for lower values of FoS whereas for higher values,the difference is marked.For the jointed rock in the Himalayan region,the nonlinear GHB criterion gives better results as compared to MC criterion and matches with the prevailing field conditions.Accordingly,some suggestions are proposed to strengthen the stability of cut slopes.

Numerical failure analysis of a continuous reinforced concrete bridge under strong earthquakes using multi-scale models

Previous failure analyses of bridges typically focus on substructure failure or superstructure failure separately. However, in an actual bridge, the seismic induced substructure failure and superstructure failure may influence each other. Moreover, previous studies typically use simplified models to analyze the bridge failure; however, there are inherent defects in the calculation accuracy compared with using a detailed three-dimensional （3D） finite element （FE） model. Conversely, a detailed 3D FE model requires more computational costs, and a proper erosion criterion of the 3D elements is necessary. In this paper, a multi-scale FE model, including a corresponding erosion criterion, is proposed and validated that can significantly reduce computational costs with high precision by modelling a pseudo-dynamic test of an reinforced concrete （RC） pier. Numerical simulations of the seismic failures of a continuous RC bridge based on the multi-scale FE modeling method using LS-DYNA are performed. The nonlinear properties of the bridge, various connection strengths and bidirectional excitations are considered. The numerical results demonstrate that the failure of the connections will induce large pounding responses of the girders. The nonlinear deformation of the piers will aggravate the pounding damages. Furthermore, bidirectional earthquakes will induce eccentric poundingsto the girders and different failure modes to the adjacent piers.