Research and Application of a Multi-Field Co-Simulation Data Extraction Method Based on Adaptive Infinitesimal Element

Regarding the spatial profile extraction method of a multi-field co-simulation dataset,different extraction directions,locations,and numbers of profileswill greatly affect the representativeness and integrity of data.In this study,a multi-field co-simulation data extractionmethod based on adaptive infinitesimal elements is proposed.Themultifield co-simulation dataset based on related infinitesimal elements is constructed,and the candidate directions of data profile extraction undergo dimension reduction by principal component analysis to determine the direction of data extraction.Based on the fireworks algorithm,the data profile with optimal representativeness is searched adaptively in different data extraction intervals to realize the adaptive calculation of data extraction micro-step length.The multi-field co-simulation data extraction process based on adaptive microelement is established and applied to the data extraction process of the multi-field co-simulation dataset of the sintering furnace.Compared with traditional data extraction methods for multi-field co-simulation,the approximate model constructed by the data extracted from the proposed method has higher construction efficiency.Meanwhile,the relative maximum absolute error,root mean square error,and coefficient of determination of the approximationmodel are better than those of the approximation model constructed by the data extracted from traditional methods,indicating higher accuracy,it is verified that the proposed method demonstrates sound adaptability and extraction efficiency.

Study on hot deformation behavior of homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy using a combination of strain-compensated Arrhenius constitutive model and finite element simulation method

Isothermal hot compression experiments were conducted on homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy to investigate hot deformation behavior at the temperature range of 673-773 K and the strain rate range of 0.001-1 s^(-1)by using a Gleeble-1500D thermo mechanical simulator.Metallographic characterization on samples deformed to true strain of 0.70 illustrates the occurrence of flow localization and/or microcrack at deformation conditions of 673 K/0.01 s^(-1),673 K/1 s^(-1)and 698 K/1 s^(-1),indicating that these three deformation conditions should be excluded during hot working of homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy.Based on the measured true stress-strain data,the strain-compensated Arrhenius constitutive model was constructed and then incorporated into UHARD subroutine of ABAQUS software to study hot deformation process of homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy.By comparison with measured force-displacement curves,the predicted results can describe well the rheological behavior of homogenized Mg-8.5Gd-4.5Y-0.8Zn-0.4Zr alloy,verifying the validity of finite element simulation of hot compression process with this complicated constitutive model.Numerical results demonstrate that the distribution of values of material parameters(α,n,Q and ln A)within deformed sample is inhomogeneous.This issue is directly correlated to the uneven distribution of equivalent plastic strain due to the friction effect.Moreover,at a given temperature the increase of strain rate would result in the decrease of equivalent plastic strain within the central region of deformed sample,which hinders the occurrence of dynamic recrystallization(DRX).

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. .

Multi-physical fields distribution in billet during helical electromagnetic stirring:A numerical simulation research

Electromagnetic stirring is one of the widely applied techniques to modify the quality of casting billets.Different from conventional rotate stirring,the helical stirring is more professional in assisting multi-dimensional flow of molten metal and eliminating solidification defects.In this study,the single-winding helical stirring(SWHS)was introduced,offering advantages such as smaller volume and lower electromagnetic shielding compared to traditional helical stirring methods.Following a comprehensive numerical simulation,the stirring parameters of SWHS were adjusted to yoke inclination angle of 43°and frequency of 12 Hz.The higher electromagnetic force and flow velocity in drawing direction,as well as the lower temperature gradient induced by the SWHS,are positive factors for homogeneous solidification of billet.The experimental results on Al-8%Si alloy and 0.4%C-1.1%Mn steel demonstrate that compared to rotate stirring,the SWHS process can induce better billet quality and is more effective in accelerating the equiaxed expansion and reducing element segregation.The SWHS process can enhance the equiaxed ratio of the billet by 58.3%and reduce segregation degree of carbon element by 10.97%.Consequently,SWHS holds great promise as a potential approach for improving the quality of continuous casting billets.

Random vibration analysis of FGM plates subjected to moving load using a refined stochastic finite element method

The article introduces a finite element procedure using the bilinear quadrilateral element or four-node rectangular element(namely Q4 element) based on a refined first-order shear deformation theory(rFSDT) and Monte Carlo simulation(MCS), so-called refined stochastic finite element method to investigate the random vibration of functionally graded material(FGM) plates subjected to the moving load.The advantage of the proposed method is to use r-FSDT to improve the accuracy of classical FSDT, satisfy the stress-free condition at the plate boundaries, and combine with MCS to analyze the vibration of the FGM plate when the parameter inputs are random quantities following a normal distribution. The obtained results show that the distribution characteristics of the vibration response of the FGM plate depend on the standard deviation of the input parameters and the velocity of the moving load.Furthermore, the numerical results in this study are expected to contribute to improving the understanding of FGM plates subjected to moving loads with uncertain input parameters.

COUPLED SIMULATION OF 3D ELECTRO-MAGNETO-FLOW FIELD IN HALL-HEROULT CELLS USING FINITE ELEMENT METHOD

Two full 3D steady mathematical models are developed by finite element method （FEM） to calcalate coupled physics fields. the electro-magnetic model is built and solved first and so is the fluid motion model with the acquired electromagnetic force as source body forces in Navier-Stokes equations. Effects caused by the ferromagnetic shell, busbar system around, and open boundary problem as well as inside induced current were considered in terms of the magnetic field. Furthermore, a new modeling method is found to set up solid models and then mesh them entirely with so-called structuralized grids, namely hex-mesh. Examples of 75kA prebaked cell with two kinds of busbar arrangements are presented. Results agree with those disclosed in the literature and confirm that the coupled simulation is valid. It is also concluded that the usage of these models facilitates the consistent analysis of the electric field to magnetic field and then flow motion to the greater extent, local distributions of current density and magnetic flux density are very much dependent on the cell structure, the steel shell is a shield to reduce the magnetic field and flow pattern is two dimensional in the main body of the metal pad.

Applications of finite element simulation in orthopedic and trauma surgery

Research in different areas of orthopedic and trauma surgery requires a methodology that allows both a more economic approach and the ability to reproduce different situations in an easy way. Simulation models have been introduced recently in bioengineering and could become an essential tool in the study of any physiological unity, regardless of its complexity. The main problem in modeling with finite elements simulation is to achieve an accurate reproduction of the anatomy and a perfect correlation of the different structures, in any region of the human body. Authors have developed a mixed technique, joining the use of a three-dimensional laser scanner Roland Picza captured together with computed tomography(CT) and 3D CT images, to achieve a perfect reproduction of the anatomy. Finite element(FE) simulation lets us know the biomechanical changes that take place after hipprostheses or osteosynthesis implantation and biological responses of bone to biomechanical changes. The simulation models are able to predict changes in bone stress distribution around the implant, so allowing preventing future pathologies. The development of a FE model of lumbar spine is another interesting application of the simulation. The model allows research on the lumbar spine, not only in physiological conditions but also simulating different load conditions, to assess the impact on biomechanics. Different degrees of disc degeneration can also be simulated to determine the impact on adjacent anatomical elements. Finally, FE models may be useful to test different fixation systems, i.e., pedicular screws, interbody devices or rigid fixations compared with the dynamic ones. We have also developed models of lumbar spine and hip joint to predict the occurrence of osteoporotic fractures, based on densitometric determinations and specific biomechanical models, including approaches from damage and fracture mechanics. FE simulations also allow us to predict the behavior of orthopedic splints applied to the correction of deformities, providing the recovering force-displacement and angle-moment curves that characterize the mechanical behavior of the splint in the overall range of movement.

RESEARCH ON THE SELECTION OF FRICTION MODELS IN THE FINITE ELEMENT SIMULATION OF WARM EXTRUSION

During the process of finite element simulation of precision warm forging, the selection of friction models has a direct effect on the precision accuracy of finite element simulation results. Among all the factors which influence the selection of friction models, the distribution rule of normal stress at the tool-workpiece interface is a key one. To find out the distribution rule of normal stress at the tool-workpiece interface, this paper has made a systematic research on three typical plastic deformation processes: forward extrusion, backward extrusion, and lateral extrusion by a method of finite element simulation. Then on the base of synthesizing and correcting traditional friction models, a new general friction model which is fit for warm extrusion is developed at last.

Finite element simulation of inertia friction welding of superalloy bars

A thermo-mechanical coupling.finite element model was built to investigate the inertia friction welding of GH4169 bars. The remeshing and map solution techniques were adopted. Ttle whole welding process was investigated by adopting an innovative heat generation model and the flywheel rotational speed measured via the experiment. The simulated evolution of axial shortening shows a good agreement with the experiment. In addition, extensive .strain concentration presents in the interface and flash, and the largest ,strain exists near the flash root. Moreover, an intere.sting thermal reflux phenomenon during the cooling stage was found.

Finite element simulation on the deep drawing of titanium thin-walled surface part

The deep drawing of titanium thin-walled surface part was simulated based on a self-developed three-dimensional finite element model. After an investigation on forming rules, a virtual orthogonal experimental design was adopted to determine the significance of processing parameters, such as die radius, blank holder force, and friction coefficient, on the forming process. The distributions of thickness and equivalent plastic strain of the drawn part were evaluated. The results show that die radius has a relative major influence on the deep drawing process, followed by friction coefficient and blank holder force.

THREE-DIMENSIONAL FINITE ELEMENT SIMULATION OF TOTAL KNEE JOINT IN GAIT CYCLE

Based on CT scanning pictures from a volunteer＇s knee joint, a three-dimensional finite element model of the healthy human knee joint is constructed including complete femur, tibia, fibular, patellar and the main cartilage and ligaments. This model was validated using experimental and numerical results obtained from other authors. The pressure distribution of contact surfaces of knee joint are calculated and analyzed under the load action of ‘heel strike＇, ‘single limb stance＇ and ‘toe-off＇. The results of the gait cycle are that the contact areas of medial cartilage are larger than that of lateral cartilage; the contact force and contact areas would grow larger with the load increasing; the pressure of lateral meniscus is steady, relative to the significant variation of peak pressure in medial meniscus; and the peak value of contact pressure on all components are usually found at about 4570 of the gait cycle.

Simulation-based Estimation of Thermal Behavior of Direct Feed Drive Mechanism with Updated Finite Element Model

Linear motors generate high heat and cause significant deformation in high speed direct feed drive mechanisms.It is relevant to estimate their deformation behavior to improve their application in precision machine tools.This paper describes a method to estimate its thermal deformation based on updated finite element（FE）model methods.Firstly,a FE model is established for a linear motor drive test rig that includes the correlation between temperature rise and its resulting deformation.The relationship between the input and output variables of the FE model is identified with a modified multivariate input/output least square support vector regression machine.Additionally,the temperature rise and displacements at some critical points on the mechanism are obtained experimentally by a system of thermocouples and an interferometer.The FE model is updated through intelligent comparison between the experimentally measured values and the results from the regression machine.The experiments for testing thermal behavior along with the updated FE model simulations is conducted on the test rig in reciprocating cycle drive conditions.The results show that the intelligently updated FE model can be implemented to analyze the temperature variation distribution of the mechanism and to estimate its thermal behavior.The accuracy of the thermal behavior estimation with the optimally updated method can be more than double that of the initial theoretical FE model.This paper provides a simulation method that is effective to estimate the thermal behavior of the direct feed drive mechanism with high accuracy.

Simulation of near-fault bedrock strong ground-motion field by explicit finite element method

Based on presumed active fault and corresponding model, this paper predicted the near-fault ground motion filed of a scenario earthquake （Mw=6 3/4 ） in an active fault by the explicit finite element method in combination with the source time function with improved transmitting artificial boundary and with high-frequency vibration contained. The results indicate that the improved artificial boundary is stable in numerical computation and the predicted strong ground motion has a consistent characteristic with the observed motion.

FINITE ELEMENT MODEL AND SIMULATION OF ROTARY FORGING OF A DISC

The rotary forging process of a disc is simulated by 3-D finite element method.The motion of the rotary the is described as the combination of a revolution round the machine axis and a spin round the rotary die axis. Therefore, the workpiece can be loaded and unloaded partly and cyclically by the cone surface of the rotary the continuously, according with the practical rotary forging process. From the siumulation rasults, the causes of center-thinning during rotary forging of discs are that the locally loading of rotary die made the workpiece center get high radial and tangential tensile stresses, and then the shortening in axial direction and the elongating in tangential and radial direction occur continuous- ly.

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.

Finite element simulation for mechanical response of surface mounted solder joints under different temperature cycling

Nonlinear finite element simulation for mechanical response of surface mounted solder joint under different temperature cycling was carried out. Seven sets of parameters were used in order to evaluate the influence of temperature cycling profile parameters. The results show that temperature cycling history has significant effect on the stress response of the solder joint. Based on the concept of relative damage stress proposed by the authors, it is found that enough high temperature holding time is necessary for designing the temperature cycling profile in accelerated thermal fatigue test.

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.

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.

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.

Evaluation of Guardrail Posts Installed in Asphalt Mow Strips by Static Finite Element Simulation

Asphalt mow strips are typically used as vegetation barriers around guardrail posts in the design of roadside safety structures. Asphalt mow strips have historically been modeled as a rigid layer in simulations;this assumption results in significant ground level restraint on the guardrail post. However, experiments have shown that asphalt rupture should be considered in the analysis of the response of guardrail posts embedded in mow strips. The present study investigates the effect of asphalt material properties and mow strip geometry on guardrail post performance using finite element simulations. Numerical simulations are performed and correlated with results from static experiments and material testing. The test simulations and experimental results are used to evaluate the response of guardrail posts with various mow strip designs to predict the level of restraint from the asphalt layer. The model is then used to investigate the effects of asphalt material properties and mow strip geometry on the overall performance of the system. The results demonstrate that including asphalt rupture in numerical simulations is essential in accurately predicting the behavior of guardrail posts installed in asphalt mow strips. In addition, mow strip geometry along with asphalt material properties significantly affect the guardrail post response.