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.

DEFORMATION ANALYSIS OF SHEET METAL SINGLE-POINT INCREMENTAL FORMING BY FINITE ELEMENT METHOD SIMULATION

Single-point incremental forming （SPIF） is an innovational sheet metal forming method without dedicated dies, which belongs to rapid prototyping technology. In generalizing the SPIF of sheet metal, the deformation analysis on forming process becomes an important and useful method for the planning of shell products, the choice of material, the design of the forming process and the planning of the forming tool. Using solid brick elements, the finite element method（FEM） model of truncated pyramid was established. Based on the theory of anisotropy and assumed strain formulation, the SPIF processes with different parameters were simulated. The resulted comparison between the simulations and the experiments shows that the FEM model is feasible and effective. Then, according to the simulated forming process, the deformation pattern of SPIF can be summarized as the combination of plane-stretching deformation and bending deformation. And the study about the process parameters＇ impact on deformation shows that the process parameter of interlayer spacing is a dominant factor on the deformation. Decreasing interlayer spacing, the strain of one step decreases and the formability of blank will be improved. With bigger interlayer spacing, the plastic deformation zone increases and the forming force will be bigger.

Flow simulation considering adsorption boundary layer based on digital rock and finite element method

Due to the low permeability of tight reservoirs,throats play a significant role in controlling fluid flow.Although many studies have been conducted to investigate fluid flow in throats in the microscale domain,comparatively fewer works have been devoted to study the effect of adsorption boundary layer(ABL)in throats based on the digital rock method.By considering an ABL,we investigate its effects on fluid flow.We build digital rock model based on computed tomography technology.Then,microscopic pore structures are extracted with watershed segmentation and pore geometries are meshed through Delaunay triangulation approach.Finally,using the meshed digital simulation model and finite element method,we investigate the effects of viscosity and thickness of ABL on microscale flow.Our results demonstrate that viscosity and thickness of ABL are major factors that significantly hinder fluid flow in throats.

Modeling and Simulation of Valve Cycle in Vein Using an Immersed Finite Element Method

A vein model was established to simulate the periodic characteristics of blood flow and valve deformation in blood-induced valve cycles.Using an immersed finite element method which was modified by a ghost fluid technique,the interaction between the vein and blood was simulated.With an independent solid solver,the contact force between vein tissues was calculated using an adhesive contact method.A benchmark simulation of the normal valve cycle validated the proposed model for a healthy vein.Both the opening orifice and blood flow rate agreed with those in the physiology.Low blood shear stress and maximum leaflet stress were also seen in the base region of the valve.On the basis of the healthy model,a diseased vein model was subsequently built to explore the sinus lesions,namely,fibrosis and atrophy which are assumed stiffening and softening of the sinus.Our results showed the opening orifice of the diseased vein was inversely proportional to the corresponding modulus of the sinus.A drop in the transvalvular pressure gradient resulted from the sinus lesion.Compared to the fibrosis,the atrophy of the sinus apparently improved the vein deformability but simultaneously accelerated the deterioration of venous disease and increased the risk of potential fracture.These results provide understandings of the normal/abnormal valve cycle in vein,and can be also helpful for the prosthesis design.

Seismic Wavelet Analysis Based on Finite Element Numerical Simulation

The practice of exploration and production has proved that explosives are excited in different surrounding rocks and the seismic wavelets collected have different characteristics. In this paper, by establishing a numerical model of the explosion in the well, using finite element analysis technology for numerical simulation, the simulation calculated the stress structure in the near-source area of the earthquake excitation, and extracted the seismic wavelet. The results show that the simulation seismic wavelet characteristics of different thin interbedded sand and mudstone structures have changed significantly. Through excitation simulation, the amplitude and spectrum information of seismic wavelets can be compared and analyzed, and the excitation parameters can be optimized. .

Calculation of effective temperature for pavement rutting using numerical simulation methods

In order to predict the long-term rutting of asphalt pavement, the effective temperature for pavement rutting is calculated using the numerical simulation method. The transient temperature field of asphalt pavement was simulated based on actual meteorological data of Nanjing. 24-hour rutting development under a transient temperature field was calculated in each month. The rutting depth accumulated under the static temperature field was also estimated and the relationship between constant temperature parameters was analyzed. Then the effective temperature for pavement rutting was determined based on the rutting equivalence principle. The results show that the monthly effective temperature is above 40 t in July and August, while in June and September it ranges from 30 to 40 Rutting development can be ignored when the monthly effective temperature is less than 30 t. The yearly effective temperature for rutting in Nanjing is around 38. 5 t. The long-term rutting prediction model based on the effective temperature can reflect the influences of meteorological factors and traffic time distribution.

Treatment of discontinuous interface in liquid-solid forming with extended finite element method

Extended finite element method(XFEM) is proposed to simulate the discontinuous interface in the liquid-solid forming process.The discontinuous interface is an important phenomenon happening in the liquid-solid forming processes and it is difficult to be simulated accurately with conventional finite element method(CFEM) because it involves solid phase and liquid phase simultaneously.XFEM is becoming more and more popular with the need of solving the discontinuous problem happening in engineering field.The implementation method of XFEM is proposed on Abaqus code by using UEL(user element) with the flowchart.The key is to modify the element stiffness in the proposed method by using UEL on the platform of Abaqus code.In contrast to XFEM used in the simulation of solidification,the geometrical and physical properties of elements were modified at the same time in our method that is beneficial to getting smooth interface transition and precise analysis results.The analysis is simplified significantly with XFEM.

Radiation heat transfer model for complex superalloy turbine blade in directional solidification process based on finite element method

For the sake of a more accurate shell boundary and calculation of radiation heat transfer in the Directional Solidification(DS) process, a radiation heat transfer model based on the Finite Element Method(FEM)is developed in this study. Key technologies, such as distinguishing boundaries automatically, local matrix and lumped heat capacity matrix, are also stated. In order to analyze the effect of withdrawing rate on DS process,the solidification processes of a complex superalloy turbine blade in the High Rate Solidification(HRS) process with different withdrawing rates are simulated; and by comparing the simulation results, it is found that the most suitable withdrawing rate is determined to be 5.0 mm·min^(-1). Finally, the accuracy and reliability of the radiation heat transfer model are verified, because of the accordance of simulation results with practical process.

An explicit method for numerical simulation of wave equations

In this paper, a method to develop a hierarchy of explicit recursion formulas for numerical simulation in an irregular grid for scalar wave equations is presented and its accuracy is illustrated via 2-D and 1-D models. Approaches to develop the stable formulas which are of 2M-order accuracy in both time and space with Mbeing a positive integer for regular grids are discussed and illustrated by constructing the second order （M= 1） and the fourth order （M = 2） recursion formulas.

An explicit method for numerical simulation of wave equations: 3D wave motion

In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer, and is illustrated by establishing second order （M=1） recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.

Sintering zone prediction in direct metal laser sintering by finite element method

A three-dimensional finite element thermal model in direct metal laser sintering(DMLS) including the effect of powder-to-solid transition were established to predict sintering zone, which benefited the determination of suitable process parameters in DMLS. The nonlinear transient model of the metals thermal conductivity for powder-to-solid transition was developed. The model uses solid thermal properties of material in liquid-phase zone, transitional ones in sintering or sintered zone and powder ones in unsintered zones of powder bed to predict, respectively. Sintering zone boundary was estimated by maximum temperature history profile. Experiments were carried out using multi-component Cu-based metal powder. Compared experimental and predicted results, the mean error of sintering depth and width are 7.8% and 14.4%, respectively, which confirms the accuracy of the FEM prediction.

MULTI-MATERIAL NUMERICAL SIMULATION OF MOVING SHOCK INTERACTING WITH CONSECUTIVE BUBBLES

Numerical simulations are performed on the interface with large deformation induced by the interaction between a moving shock and two consecutive bubbles. The high performance of the level set method for multi-material interfaces is demonstrated. Discontinuous Galerkin finite element method is used to solve Euleri- an equations. And the fifth-order weighted essentially non-oscillatory （WENO） scheme is used to solve the level set equation for capturing multi-material interfaces. The ghost fluid method is used to deal with the interfacial boundary condition. Results are obtained for two bubble interacting with a moving shock. The contours of the constant density and the pressure at different time are given. In the computational domain, three different cases are considered, i.e. two helium bubbles, a helium bubble followed by an R22 bubble in the direction of the moving shock, and an R22 bubble followed by a helium bubble. Computational results indicate that multi-mate- rial interfaces can be properly captured by the level set method. Therefore, for problems involving the flow of three different materials with two different interfaces, each interface separating two different materials can be similarly handled.

COUPLED NUMERICAL SIMULATION ON COLD ROLLER'S TEMPERATURE FIFLD-PHASE TRANSFORMATION - STRESS FIELD DURING ITS QUENCHING PSOCESS

Complicated changns occur inside the steel parts during quenching process. The abruptly changed boundary conditions make the temperature field,micro - structure and stress field change dramatically in very short time, and these variables take a contact interactions in the whole process. In this paper, a three dimensional non - linear mathematical model for queeching process has been founded and the numerical simulation on temperature field,microstructre and stress field has been realized.In the FEM analysis, the incremental iteration method is used to deal with such complicated nonlinear as boundary nonlinear, physical property nonlinear,transformation nonlinear etc.The effect of stress on transformation kinetics has been considered in the calculation of microstructure. In the stress field anal- ysis,a thermo- elasto - plastic model has been founded, which considers such factors as transforma- tion strain,transformation plastic strain, themal strain and the effect of temperature and transforma- tion on mechanical propertier etc. The transient temperature field, microstructure distribution and stress field of the roller on any time can be displayed vividly,and the cooling curve and the changes of stress on any position can also be given.

Review of Numerical Simulation of TGO Growth in Thermal Barrier Coatings

Thermally grown oxide(TGO)is a critical factor for the service life of thermal barrier coatings(TBC).Numerical simulations of the growth process of TGO have become an effective means of comprehensively understanding the progressive damage of the TBC system.At present,technologies of numerical simulation to TGO growth include two categories:coupled chemical-mechanical methods and mechanical equivalent methods.The former is based on the diffusion analysis of oxidizing elements,which can describe the influence of bond coat(BC)consumption and phase transformation in the growth process of TGO on the mechanical behavior of each layer of TBC,and has high accuracy for the thickness evolution of TGO,but they cannot describe the lateral growth of TGO and the rumpling phenomenon induced.The latter focuses on describing the final stress and strain state after the growth of a specific TGO rather than the complete growth processes of TGO.Based on the measured TGO thickness growth curve,simulations of thickening and lateral growth can be achieved by directly applying anisotropic volumetric strain to oxidized elements and switching elements properties from the BC to the TGO.

Study on Numerical Simulation for Control of Winding Process of Thin Wall Spiral Tube

Being aimed at the inside wall wrinkling and sinking phenomenon of palladium-yttrium alloy thin wall spiral tube used for preparation of high purity hydrogen, extraction of hydrogen isotope, and purification and separation of hydrogen in the winding process, this article analyzed the reasons for above phenomena, established a numerical simulation model of winding process of above tube, using elastic-plastic Finite Element method analyzed the max. tensile stress and max. compression stress and their locations, thereby provides a theory base for the control of working forming course of thin wall spiral tube.

The Numerical Simulation of Flow Field of Jet System

In this paper, the internal fluid motion of a jet system is described by the Navier Stokes mechanics equations. For the simulation of the motion, the penalty function finite element method is used, and the velocity vectors and stream function curves are obtained. Using the Prandtl theory, this paper derives the free jet velocity and the jet bunch width in a half-space, the latter of which is amended by experiment. The results obtained in this paper are applied to micro-type high pressure water jet cleaner and the ejector of rocket engine.

Numerical simulation of single roller melt spinning for NdFeB alloy based on finite element method

The numerical simulation model of single roller melt spinning for rapid quenching process of NdFeB alloy was built,and the vacuum chamber,cooling roller and sample were taken into account as a system.The existing mature technology was in order to verify the correctness of simulation.The rapid quenching ribbons with different roll speeds were used as the simulation objects.The results of the numerical simulation and experiments show that the validity of the model has been testified and the reasons of the formation of complete quenching ribbons and by-product have been explained.The experimental thickness of the ribbons is proportional to the theoretical thickness.In the same spray condition,with the roll speed increasing,the thickness decreases linearly.At the speed range of25-30 m·s^(-1),the simulated calculation date is close to the experimental date,which can be considered as an ideal technological parameter.

Numerical Simulation of a Multi-Frequency Resistivity Logging-While-Drilling Tool Using a Highly Accurate and Adaptive Higher-Order Finite Element Method

A novel,highly efficient and accurate adaptive higher-order finite element method(hp-FEM)is used to simulate a multi-frequency resistivity loggingwhile-drilling(LWD)tool response in a borehole environment.Presented in this study are the vector expression of Maxwell’s equations,three kinds of boundary conditions,stability weak formulation of Maxwell’s equations,and automatic hpadaptivity strategy.The new hp-FEM can select optimal refinement and calculation strategies based on the practical formation model and error estimation.Numerical experiments show that the new hp-FEM has an exponential convergence rate in terms of relative error in a user-prescribed quantity of interest against the degrees of freedom,which provides more accurate results than those obtained using the adaptive h-FEM.The numerical results illustrate the high efficiency and accuracy of the method at a given LWD tool structure and parameters in different physical models,which further confirm the accuracy of the results using the Hermes library(http://hpfem.org/hermes)with a multi-frequency resistivity LWD tool response in a borehole environment.

Numerical simulation of hydraulic fracturing and associated microseismicity using finite-discrete element method

Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid injection, which creates an interconnected fracture network and increases the hydrocarbonproduction. Meanwhile, microseismic （MS） monitoring is one of the most effective approaches to evaluatesuch stimulation process. In this paper, the combined finite-discrete element method （FDEM） isadopted to numerically simulate HF and associated MS. Several post-processing tools, includingfrequency-magnitude distribution （b-value）, fractal dimension （D-value）, and seismic events clustering,are utilized to interpret numerical results. A non-parametric clustering algorithm designed specificallyfor FDEM is used to reduce the mesh dependency and extract more realistic seismic information.Simulation results indicated that at the local scale, the HF process tends to propagate following the rockmass discontinuities; while at the reservoir scale, it tends to develop in the direction parallel to themaximum in-situ stress. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.