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

Influence of Various Earth-Retaining Walls on the Dynamic Response Comparison Based on 3D Modeling

The present work aims to assess earthquake-induced earth-retaining(ER)wall displacement.This study is on the dynamics analysis of various earth-retaining wall designs in hollow precast concrete panels,reinforcement concrete facing panels,and gravity-type earth-retaining walls.The finite element(FE)simulations utilized a 3D plane strain condition to model full-scale ER walls and numerous nonlinear dynamics analyses.The seismic performance of differentmodels,which includes reinforcement concrete panels and gravity-type and hollowprecast concrete ER walls,was simulated and examined using the FE approach.It also displays comparative studies such as stress distribution,deflection of the wall,acceleration across the wall height,lateral wall displacement,lateral wall pressure,and backfill plastic strain.Three components of the created ER walls were found throughout this research procedure.One is a granular reinforcement backfill,while the other is a wall-facing panel and base foundation.The dynamic response effects of varied earth-retaining walls have also been studied.It was discovered that the facing panel of the model significantly impacts the earthquake-induced displacement of ER walls.The proposed analytical model’s validity has been evaluated and compared with the reinforcement concrete facing panels,gravity-type ER wall,scientifically available data,and American Association of State Highway and Transportation Officials(AASHTO)guidelines results based on FE simulation.The results of the observations indicate that the hollow prefabricated concrete ER wall is the most feasible option due to its lower displacement and high-stress distribution compared to the two types.The methodology and results of this study establish standards for future analogous investigations and professionals,particularly in light of the increasing computational capabilities of desktop computers.

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.

Vehicle–track–tunnel dynamic interaction:a finite/infinite element modelling method

The aim of this study is to develop coupled matrix formulations to characterize the dynamic interaction between the vehicle,track,and tunnel.The vehicle–track coupled system is established in light of vehicle–track coupled dynamics theory.The physical characteristics and mechanical behavior of tunnel segments and rings are modeled by the finite element method,while the soil layers of the vehicle–track–tunnel(VTT)system are modeled as an assemblage of 3-D mapping infinite elements by satisfying the boundary conditions at the infinite area.With novelty,the tunnel components,such as rings and segments,have been coupled to the vehicle–track systems using a matrix coupling method for finite elements.The responses of sub-systems included in the VTT interaction are obtained simultaneously to guarantee the solution accuracy.To relieve the computer storage and save the CPU time for the large-scale VTT dynamics system with high degrees of freedoms,a cyclic calculation method is introduced.Apart from model validations,the necessity of considering the tunnel substructures such as rings and segments is demonstrated.In addition,the maximum number of elements in the tunnel segment is confirmed by numerical simulations.

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.

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.

Deformation Stability of GH4033 Superalloy in the Hot Continuous Rolling Process Based on Dynamic Material Model and Finite Element Model

The flow stress behavior of GH4033 superalloy was determined by the hot compression tests at the temperatures of 1223-1473 K and the total strains of 0.6 with the strain rates of 0.001-30.0 s^(-1) by using cylindrical samples.The processing maps based on the dynamic material model(DMM)combined with the corresponding microstructure observations indicate the reasonable processing domain locating at the strain rates of 0.1-1.0 s^(-1) and the deformation temperature of 1273-1423 K.Meanwhile,the numerical simulation based on finite element model(FEM)described the variation of the effective strain,effective strain rate and the temperature for the core node,and unveiled the influence of the hot rolling parameters considering the initial temperature(T_(0))range of 1223-1473 K and the first-stand biting velocity(v_(0))range of 0.15-0.35 m·s^(-1).Furthermore,the deformation stability of GH4033 superalloy in the round rod hot continuous rolling(HCR)process is described and analyzed by coupling the three-dimensional(3-D)processing map,and the spatial trajectory lines were determined by the numerically simulated temperatures,the strains and the strain rates.Finally,the results show that the hot deformation stability of GH4033 can be achieved by the rolling process parameters located at T_(0)=1423 K and v_(0)=0.25 m·s^(-1).Additionally,the practical HCR processes as T_(0)=1423 K and v_(0)=0.15,0.25,0.35 m·s^(-1) were operated to verify the influence of the hot rolling parameters on the hot deformation stability by the microstructure observation of the final products.

PRINCIPAL COMPONENT DECOMPOSITION BASED FINITE ELEMENT MODEL UPDATING FOR STRAIN-RATE-DEPENDENCE NONLINEAR DYNAMIC PROBLEMS

Thin wail component is utilized to absorb impact energy of a structure. However, the dynamic behavior of such thin-walled structure is highly non-linear with material, geometry and boundary non-linearity. A model updating and validation procedure is proposed to build accurate finite element model of a frame structure with a non-linear thin-walled component for dynamic analysis. Design of experiments （DOE） and principal component decomposition （PCD） approach are applied to extract dynamic feature from nonlinear impact response for correlation of impact test result and FE model of the non-linear structure. A strain-rate-dependent non-linear model updating method is then developed to build accurate FE model of the structure. Computer simulation and a real frame structure with a highly non-linear thin-walled component are employed to demonstrate the feasibility and effectiveness of the proposed approach.

Multi-field dynamic modeling and numerical simulation of aluminum alloy resistance spot welding

In order to explore the influence of welding parameters and to investigate the Al alloy (AA) nugget formation process, a comprehensive model involving electrical-thermal-mechanical and metallurgical analysis was established to numerically display the resistance spot welding (RSW) process within multiple fields and understand the AA-RSW physics. A multi-disciplinary finite element method (FEM) framework and a empirical sub-model were built to analyze the affecting factors on weld nugget and the underlying nature of welding physics with dynamic simulation procedure. Specifically, a counter-intuitive phenomenon of the resistance time-variation caused by the transient inverse virtual variation (TIVV) effect was highlighted and analyzed on the basis of welding current and temperature distribution simulation. The empirical model describing the TIVV phenomenon was used for modifying the dynamic resistance simulation during the AA spot welding process. The numerical and experimental results show that the proposed multi-field FEM model agrees with the measured AA welding feature, and the modified dynamic resistance model captures the physics of nugget growth and the electrical-thermal behavior under varying welding current and fluctuating heat input.

Analysis of the Elastic-Plastic Dynamic Finite Element on a Jointed Rock Mass

Aim To study the elastic plastic dynamical constitutive relations about a jointed rock mass under explosion load and its computer simulation. Methods\ Stress history is taken into account and stresses will follow changes in time during a period of explosion load. According to the principle of static force balance, the corresponding nodal concentrated force is calculated and the nodal displacement is counted. The elastic plastic dynamic finite element equations are thus obtained. Results\ A finite element method is given for a jointed rock mass under explosion load. Conclusion\ The problem of large plastic deformation for jointed rock mass on blasting was efficiently resolved through dynamic finite element analysis and the range of damages by blasting simulated, and this pushes forward the problem to engineering practice.

Simulation Research of Fault Model of Detecting Rotor Dynamic Eccentricity in Brushless DC Motor Based on Motor Current Signature Analysis

基于Ansoft／Maxwell设置动态偏心故障，建立求解电机电感和磁链的有限元模型，通过仿真，证明了将感应电机动态偏心故障的特征频率经过简化后，同样适用于无刷直流电动机。基于Ansoft／Simplorer建立无刷直流电动机系统的仿真模型。在电机稳态运行下，对定子电流进行傅里叶分析，研究并建立基于定子电流监测动态偏心故障的仿真模型：动态偏心故障与特征频率的关系、动态偏心故障程度与特征频率幅值的关系。进而研究了无刷直流电动机稳态运行时转速波动对偏心故障监测的影响。仿真结果表明，转子偏心程度加大，特征频率的幅值增加。

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 model for arch bridge vibration dynamics considering effect of suspender length adjustment on geometry stiffness matrix

In this paper, we established a finite element (FEM) model to analyze the dynamic characteristics of arch bridges. In this model, the effects of adjustment to the length of a suspender on its geometry stiffness matrix are stressed. The FEM equations of mechanics characteristics, natural frequency and main mode are set up based on the first order matrix perturbation theory. Applicantion of the proposed model to analyze a real arch bridge proved the improvement in the simulation precision of dynamical characteristics of the arch bridge by considering the effects of suspender length variation.

Finite Element Simulation of In-Stent Restenosis with Tissue Growth Model

In this study, a finite element simulation of in-stent restenosis (ISR) is conducted to simulate the deployment and expansion of a stent in an occluded artery with a contact model and a mechanics-based growth model. A tissue growth model based on the multiplicative decomposition of deformation is applied to investigate the growth of the plaque and artery wall upon the stent’s implantation. Due to the high stresses at the contact points between the stent struts and the tissue, further tissue injury or restenosis is observed. The simulation results show that after the stent deployment, the von Mises stress is significantly larger in the plaque compared to the artery wall, especially in the region that is in contact with the stent. However, the growth of the plaque and artery tends to even out the stress concentration over time. The tissue growth is found to be more significant near the inner wall than the outer layer. A 0.77 mm restenosis is predicted, which agrees with published clinical observations. The features of the artery growth are carefully analyzed, and the underlying mechanism is discussed. This study is the first attempt to apply finite element analysis to artery restenosis, which establishes a framework for predicting ISR’s occurrence and severity. The results also provide insights into understanding the underlying mechanism of in-stent restenosis.

DYNAMIC MODELLING OF BAR-GEAR MIXED MULTIBODY SYSTEMS USING A SPECIFIC FINITE ELEMENT METHOD

A new dynamic model for mixed, flexible bar and gear multibody systems is developed based on a specific finite element method, and a new gear element is proposed. The gear element can take into account the time variant stiffness, the gear errors and mass unbalance. The model for geared multibody systems can couple the gear meshing and the flexibility of all contained components. The kinematic and dynamic analyses of the geared multibody systems are expounded and illustrated on an example composed of three gears, two bars and one slider.

An Integrated Dynamic Model of Ocean Mining System and Fast Simulation of Its Longitudinal Reciprocating Motion

An integrated dynamic model of China＇s deep ocean mining system is developed and the fast simulation analysis of its longitudinal reciprocating motion operation processes is achieved. The seafloor tracked miner is built as a three-dimensional single-body model with six-degree-of-freedom. The track-terrain interaction is modeled by partitioning the track-terrain interface into a certain number of mesh elements with three mutually perpendicular forces, including the normal force, the longitudinal shear force and the lateral shear force, acting on the center point of each mesh element. The hydrodynamic force of the miner is considered and applied. By considering the operational safety and collection efficiency, two new mining paths for the miner on the seafloor are proposed, which can be simulated with the established single-body dynamic model of the miner. The pipeline subsystem is built as a three-dimensional multi-body discrete element model, which is divided into rigid elements linked by flexible connectors. The flexible connector without mass is represented by six spring-damper elements. The external hydrodynamic forces of the ocean current from the longitudinal and lateral directions are both considered and modeled based on the Morison formula and applied to the mass center of each corresponding discrete rigid element. The mining ship is simplified and represented by a general kinematic point, whose heave motion induced by the ocean waves and the longitudinal and lateral towing motions are considered and applied. By integrating the single-body dynamic model of the miner and the multi-body discrete element dynamic model of the pipeline, and defining the kinematic equations of the mining ship, the integrated dynamic model of the total deep ocean mining system is formed. The longitudinal reciprocating motion operation modes of the total mining system, which combine the active straight-line and turning motions of the miner and the ship, and the passive towed motions of the pipeline, are proposed and simulated with the developed 3D dynamic model. Some critical simulation results are obtained and analyzed, such as the motion trajectories of key subsystems, the velocities of the buoyancy modules and the interaction forces between subsystems, which in a way can provide important theoretical basis and useful technical reference for the practical deep ocean mining system analysis, operation and control.

Numerical simulation analysis for deformation deviation and experimental verification for an antenna thin-wall parts considering riveting assembly with finite element method

In the process of thin-wall parts assembly for an antenna,the parts assembly deformation deviation is occurring due to the riveting assembly.In view of the riveting assembly deformation problems,it can be analyzed through transient and static simulation.In this work,the theoretical deformation model for riveting assembly is established with round head rivet.The simulation analysis for riveting deformation is carried out with the riveting assembly piece including four rivets,which comparing with the measuring points experiment results of riveting test piece through dealing with the experimental data using the point coordinate transform method and the space line fitting method.Simultaneously,the deformation deviation of the overall thin-wall parts assembly structure is analyzed through finite element simulation;and its results are verified by the measuring experiment for riveting assembly with the deformation deviation of some key points on the thin-wall parts.Through the comparison analysis,it is shown that the simulation results agree well with the experimental results,which proves the correctness and effectiveness of the theoretical analysis,simulation results and the given experiment data processing method.Through the study on the riveting assembly for thin-wall parts,it will provide a theoretical foundation for improving thin-wall parts assembly quality of large antenna in future.

Multiscale Hybrid-Mixed Finite Element Method for Flow Simulation in Fractured Porous Media

The multiscale hybrid-mixed(MHM)method is applied to the numerical approximation of two-dimensional matrix fluid flow in porous media with fractures.The two-dimensional fluid flow in the reservoir and the one-dimensional flow in the discrete fractures are approximated using mixed finite elements.The coupling of the two-dimensional matrix flow with the one-dimensional fracture flow is enforced using the pressure of the one-dimensional flow as a Lagrange multiplier to express the conservation of fluid transfer between the fracture flow and the divergence of the one-dimensional fracture flux.A zero-dimensional pressure(point element)is used to express conservation of mass where fractures intersect.The issuing simulation is then reduced using the MHM method leading to accurate results with a very reduced number of global equations.A general system was developed where fracture geometries and conductivities are specified in an input file and meshes are generated using the public domain mesh generator GMsh.Several test cases illustrate the effectiveness of the proposed approach by comparing the multiscale results with direct simulations.

Dynamic Characteristics of Inflated Torus Using Finite Element Analysis

In recent years, inflatable structures have been a subject of interest for space applications such as communication antenna, solar thermal propulsion and entry/landing systems. The inflatable structures characterized by high strength-to-mass ratios, minimal stowage volume, which makes them suitable for cost-effective large space structures. A typical example for the inflatable structure is the inflated torus which often used in order to provide structure support. In this study, our main focus is to understand the dynamic characteristics of an inflated torus in order to formulate an accurate mathematical model suitable for active vibration control applications. A commercial finite element package, ANSYS, is used to model the inflated torus. To verify the model the obtained frequencies and mode shapes are compared with the published results, which are derived using analytical approach, the verification shows a good agreement between the FEM and the analytical results. Based on the verified model, parametric study was investigated; the material thickness increase causes the natural frequencies decrease, while the increase of the inflation pressure simply results in stiffening the ring, which means that the natural frequency increased. The FEM analysis gives an easy and fast way for the vibration analysis of the structures compared with the complicated analytical solutions.

Study on Dynamic Mechanical Behavior of Al-Mg-Si Alloy

The dynamic mechanical behavior of Al-Mg-Si alloy was investigated under different strain rates by mechanical property and microstructure characterization,constitutive behavior analysis and numerical simulation in the present study.As the strain rate increases,the yield strength,ultimate tensile strength and elongation increase first,then remain almost constant,and finally increase.The alloy always exhibits a typical ductile fracture mode,not depending on the strain rate.However,as the strain rate increases,the number of dimples gradually increases.Tensile deformation can refine grains,however,the grain structure is slightly affected by the strain rate.An optimized Johnson-Cook constitutive equation was used to describe the mechanical behavior and obtained by fitting the true stress-strain curves.The parameter C was described by a function related to the strain rate.The fitting true stress-strain curves by the JC model agree very well with the experimental true stress-strain curves.The true stress-strain curves calculated by the finite element numerical simulation agree well with the experimental true stress-strain curves.