Study on Postseismic Impact of Wenchuan Earthquake on the Sichuan-Yunnan Region Based on Three-Dimensional Viscoelastic Finite Element Method

Based on the lateral segmentation and vertical stratification characteristics of the crustal medium in Sichuan-Yunnan region,and the asymmetry of the static dislocations on the coseismic fault plane of the Wenchuan M_S8.0 earthquake,we built a three-dimensional viscoelastic finite element model of the crust in the Sichuan-Yunnan region.The postseismic impact of the Wenchuan M_S8.0 earthquake on the Sichuan-Yunnan region was studied.The results show that:①The far-field horizontal deformation caused by the viscoelastic relaxation of the medium in the 10 years after the earthquake is about 0-20 mm within the Sichuan-Yunnan diamond-shaped block,which has a greater influence on north side and smaller on south side.②In the 10 years after the earthquake,the far-field vertical deformation caused by the viscoelastic relaxation effect of the medium is small,and it shows an increase of about 0-4 mm in most areas of the Sichuan-Yunnan diamond-shaped block.③The Xianshuihe fault and the eastern segment of the East Kunlun fault,which are close to the seismogenic fault,show a high gradient on deformation fields after the earthquake.④In order to compare with the strong earthquake activity in the SichuanYunnan region after the Wenchuan earthquake,the horizontal stress state and the Coulomb failure stress change of the active block boundary are also calculated.From the spatial distribution of the coseismic and postseismic displacement field,the fault activity characteristics reflected by the stress state and the stress loading of the fault layer reflected by the Coulomb failure stress change,there is a certain correlation with the spatial distribution of strong earthquake activity in this region.

Hermite Finite Element Method for Vibration Problem of Euler-Bernoulli Beam on Viscoelastic Pasternak Foundation

Viscoelastic foundation plays a very important role in civil engineering. It can effectively disperse the structural load into the foundation soil and avoid the damage caused by the concentrated load. The model of Euler-Bernoulli beam on viscoelastic Pasternak foundation can be used to analyze the deformation and response of buildings under complex geological conditions. In this paper, we use Hermite finite element method to get the numerical approximation scheme for the vibration equation of viscoelastic Pasternak foundation beam. Convergence and error estimation are rigourously established. We prove that the fully discrete scheme has convergence order O(τ2+h4), where τis time step size and his space step size. Finally, we give four numerical examples to verify the validity of theoretical analysis.

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.

Three-dimensional Finite Element Modeling of a Maxillary Premolar Tooth Based on the Micro-CT Scanning: A Detailed Description

This study describes the details of how to construct a three-dimensional （3D） finite element model of a maxillary first premolar tooth based on micro-CT data acquisition technique, MIMICS soft- ware and ANSYS software. The tooth was scanned by micro-CT, in which 1295 slices were obtained and then 648 slices were selected for modeling. The 3D surface mesh models of enamel and dentin were created by MIMICS （STL file）. The solid mesh model was constructed by ANSYS. After the material properties and boundary conditions were set, a loading analysis was performed to demonstrate the ap- plicableness of the resulting model. The first and third principal stresses were then evaluated. The re- suits showed that the number of nodes and elements of the finite element model were 56 618 and 311801, respectively. The geometric form of the model was highly consistent with that of the true tooth, and the deviation between them was ~）.28%. The loading analysis revealed the typical stress patterns in the contour map. The maximum compressive stress existed in the contact points and the maximum tensile stress existed in the deep fissure between the two cusps. It is concluded that by using the micro-CT and highly integrated software, construction of the 3D finite element model with high quality will not be difficult for clinical researchers.

Validation and application of three-dimensional discontinuous deformation analysis with tetrahedron finite element meshed block

In the last decade, three dimensional discontin- uous deformation analyses （3D DDA） has attracted more and more attention of researchers and geotechnical engineers worldwide. The original DDA formulation utilizes a linear displacement function to describe the block movement and deformation, which would cause block expansion under rigid body rotation and thus limit its capability to model block de- formation. In this paper, 3D DDA is coupled with tetrahe- dron finite elements to tackle these two problems. Tetrahe- dron is the simplest in the 3D domain and makes it easy to implement automatic discretization, even for complex topol- ogy shape. Furthermore, element faces will remain planar and element edges will remain straight after deformation for tetrahedron finite elements and polyhedral contact detection schemes can be used directly. The matrices of equilibrium equations for this coupled method are given in detail and an effective contact searching algorithm is suggested. Valida- tion is conducted by comparing the results of the proposed coupled method with that of physical model tests using one of the most common failure modes, i.e., wedge failure. Most of the failure modes predicted by the coupled method agree with the physical model results except for 4 cases out of the total 65 cases. Finally, a complex rockslide example demon- strates the robustness and versatility of the coupled method.

Finite element modeling of the viscoelastic responses of the eye during microvolumetric changes

A linear viscoelastic finite element model was built to investigate factors that influenced the intraocular pressure (IOP) elevations due to micro-volumetric changes in the eye at three different rates. The viscoelastic properties of the cornea and the sclera, including the instantaneous modulus, equilibrium modulus, and relaxation time constants, parametrically varied to examine their effects on IOP elevations at different rates of volumetric changes. The simulated responses were in good agreement with the previously reported experimental results obtained from porcine globes, showing the general trend of higher IOP elevations at faster rates. The simulations showed that all viscoelastic properties influenced the profile of the dynamic IOP due to volumetric changes, and the relative significance of a specific parameter was highly dependent on the rate of change.

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.

A mixed finite element scheme for viscoelastic flows with XPP model

A mixed finite element formulation for viscoelastic flows is derived in this paper, in which the FIC （finite incremental calculus） pressure stabilization process and the DEVSS （discrete elastic viscous stress splitting） method using the Crank-Nicolson-based split are introduced within a general framework of the iterative version of the fractional step algorithm. The SU （streamline-upwind） method is particularly chosen to tackle the convective terms in constitutive equations of viscoelastic flows. Thanks to the proposed scheme the finite elements with equal low-order interpolation approximations for stress-velocity-pressure variables can be successfully used even for viscoelastic flows with high Weissenberg numbers. The XPP （extended Pom-Pom） constitutive model for describing viscoelastic behaviors is particularly integrated into the proposed scheme. The numerical results for the 4：1 sudden contraction flow problem demonstrate prominent stability, accuracy and convergence rate of the proposed scheme in both pressure and stress distributions over the flow domain within a wide range of the Weissenberg number, particularly the capability in reproducing the results, which can be used to explain the ＂die swell＂ phenomenon observed in the polymer injection molding process.

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.

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.

Simulation of three-dimensional tension-induced cracks based on cracking potential function-incorporated extended finite element method

In the finite element method,the numerical simulation of three-dimensional crack propagation is relatively rare,and it is often realized by commercial programs.In addition to the geometric complexity,the determination of the cracking direction constitutes a great challenge.In most cases,the local stress state provides the fundamental criterion to judge the presence of cracks and the direction of crack propagation.However,in the case of three-dimensional analysis,the coordination relationship between grid elements due to occurrence of cracks becomes a difficult problem for this method.In this paper,based on the extended finite element method,the stress-related function field is introduced into the calculation domain,and then the boundary value problem of the function is solved.Subsequently,the envelope surface of all propagation directions can be obtained at one time.At last,the possible surface can be selected as the direction of crack development.Based on the aforementioned procedure,such method greatly reduces the programming complexity of tracking the crack propagation.As a suitable method for simulating tension-induced failure,it can simulate multiple cracks simultaneously.

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.

Alternating Direction Finite Volume Element Methods for Three-Dimensional Parabolic Equations

This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes, one is analogous to Douglas finite difference scheme with second-order splitting error, the other two schemes have third-order splitting error, and the last one is an extended LOD scheme. The L2 norm and H1 semi-norm error estimates are obtained for the first scheme and second one, respectively. Finally, two numerical examples are provided to illustrate the efficiency and accuracy of the methods.

Finite Element Modelling of Car Seat with Hyperelastic and Viscoelastic Foam Material Properties to Assess Vertical Vibration in Terms of Acceleration

Primary objective of automobile seats is to offer adequate level of safety and comfort to the seated human occupant, primarily against vibration. Ideally, any sort of automotive seat is constructed by mechanical framework, cushion, backrest and headrest. The frame structures are made of metallic alloys, while the cushion, backrest and headrest are made of polyurethane foam material. During the design phase of automotive seat, the greatest challenge is to assign realistic material properties to foam material;as it is non-linear in nature and exhibit hysteresis at low level stress. In this research paper, a car seat has been modelled in finite element environment by implementing both hyperelastic and viscoelastic material properties to polyurethane foam. The car seat has been excited with the loads due to car acceleration and human object and the effects of vibration in terms of vertical acceleration at different locations have been measured. The aims of this simulation study are to establish a car seat with the foam material properties as accurately as possible and provide a finite element set up of car seat to monitor the vertical acceleration responses in a reasonable way. The RMS acceleration values for headrest, backrest and cushion have been found to be 0.91 mm/sec2, 0.54 mm/sec2 and 0.47 mm/sec2, respectively, which showed that the car seat foam can effectively be modelled through combined hyperelastic and viscoelastic material formulations. The simulation outputs have been validated through real life testing data, which clearly indicates that this computerized simulation technique is capable of anticipating the acceleration responses at different car seat segments in a justified way.

Three-Dimensional Thermo-Elastic-Plastic Finite Element Method Modeling for Predicting Weld-Induced Residual Stresses and Distortions in Steel Stiffened-Plate Structures

The objective of the present paper is to develop nonlinear finite element method models for predicting the weld-induced initial deflection and residual stress of plating in steel stiffened-plate structures. For this purpose, three-dimensional thermo-elastic-plastic finite element method computations are performed with varying plate thickness and weld bead length (leg length) in welded plate panels, the latter being associated with weld heat input. The finite element models are verified by a comparison with experimental database which was obtained by the authors in separate studies with full scale measurements. It is concluded that the nonlinear finite element method models developed in the present paper are very accurate in terms of predicting the weld-induced initial imperfections of steel stiffened plate structures. Details of the numerical computations together with test database are documented.

CHEBYSHEV PSEUDOSPECTRAL-HYBRID FINITE ELEMENT METHOD FOR THREE-DIMENSIONAL VORTICITY EQUATION

In this paper,Chebyshev pseudospectral-finite element schemes are proposed for solving three dimensional vorticity equation.Some approximation results in nonisotropic Sobolev spaces are given.The generalized stability and the convergence are proved strictly.The numerical results show the advantages of this method.The technique in this paper is also applicable to other three-dimensional nonlinear problems in fluid dynamics.

THE CHARACTERISTIC FINITE ELEMENT ALTERNATING-DIRECTION METHOD AND ANALYSIS FOR THREE-DIMENSIONAL NUMERICAL RESERVOIR SIMULATION

Petroleum science has made remarkable progress in organic geochcmistry and in the research into the theories of petroleum origin, its transport and accumulation. In estimating the oil-gas resources of a basin, the knowledge of its evolutionary history and especially the numerical computation of fluid flow and the history of its changes under heat is vital. The mathematical model can be described as a coupled system of nonlinear partial differentical equations with initial-boundary value problems. This thesis, from actual conditions such as the effect of fluid compressibility and the three-dimensional characteristic of large-scale science-engineering computation, we put forward a kind of characteristic finite element alternating-direction schemes and obtain optimal order estimates in L^2 norm for the error in the approximate assumption.

Three-dimensional anisotropic induced polarization modeling using finite element method

With the development of geophysical exploration technology,the anisotropy of underground media has got more and more attention.At present,there are few studies on the anisotropy of the induced polarization method.This article explores the effect of anisotropy on the underground media of the induced polarization method under three-dimensional complex terrain.The research work transforms the underground electric field control equation into a variational problem,and use the unstructured finite element method to construct a large linear equation system for solving electric potentials.By the sparse matrix compression technique and symmetric successive over-relaxation preconditioned conjugate gradient algorithm(SSOR-PCG)to solve the equation system.Finally,the article uses the classic central gradient array method to obtain the forward apparent polarizability value.The calculation results of the model find that different anisotropic conditions will significantly affect the forward results which show a strong directional correlation,revealing the importance of considering anisotropy in practical work.

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.

Three-dimensional finite element analysls for external midface distraction after different osteotomy in patients with cleft lip and palate

Objective To study three - dimensional finite element analysis for external midface distraction after different osteotomy in patients with cleft lip and palate ( Clp) . Methods Three - dimensional Fem models of Le Fort Ⅰ,Ⅱand Ⅲ,osteotomy in Clp patients were estabolished. External midface distraction were simulated. An anteriorly and inferiorly directed 900 g force was