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 too...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.展开更多
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 ...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.展开更多
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 thr...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.展开更多
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...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).展开更多
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 whi...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.展开更多
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 fi...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.展开更多
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...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.展开更多
The metal cutting process is accompanied by complex stress field,strain field,temperature field.The comprehensive effects of process parameters on chip morphology,cutting force,tool wear and residual stress are comple...The metal cutting process is accompanied by complex stress field,strain field,temperature field.The comprehensive effects of process parameters on chip morphology,cutting force,tool wear and residual stress are complex and inter-connected.Finite element method(FEM)is considered as an effective method to predict process variables and reveal microscopic physical phenomena in the cutting process.Therefore,the finite element(FE)simulation is used to research the conventional and micro scale cutting process,and the differences in the establishment of process variable FE simulation models are distinguished,thereby improving the accuracy of FE simulation.The reliability and effectiveness of FE simulation model largely depend on the accuracy of the simulation method,constitutive model,friction model,damage model in describing mesh element,the dynamic mechanical behavior of materials,the tool-chip-workpiece contact process and the chip formation mechanism.In this paper,the FE models of conventional and micro process variables are comprehensively and up-to-date reviewed for different materials and machining methods.The purpose is to establish a FE model that is more in line with the real cutting conditions,and to provide the possibility for optimizing the cutting process variables.The development direction of FE simulation of metal cutting process is discussed,which provides guidance for future cutting process modeling.展开更多
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-dimensio...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.展开更多
In this paper,a stabilized finite element technique,actualized by streamline upwind Petrov-Galerkin(SUPG)stabilized method and three-step finite element method(FEM),for large eddy simulation(LES)is developed to predic...In this paper,a stabilized finite element technique,actualized by streamline upwind Petrov-Galerkin(SUPG)stabilized method and three-step finite element method(FEM),for large eddy simulation(LES)is developed to predict the wind flow with high Reynolds numbers.Weak form of LES motion equation is combined with the SUPG stabilized term for the spatial finite element discretization.An explicit three-step scheme is implemented for the temporal discretization.For the numerical example of 2D wind flow over a square rib at Re=4.2×105,the Smagorinsky's subgrid-scale(SSGS)model,the DSGS model,and the DSGS model with Cabot near-wall model are applied,and their results are analyzed and compared with experimental results.Furthermore,numerical examples of 3D wind flow around a surface-mounted cube with different Reynolds numbers are performed using DSGS model with Cabot near-wall model based on the present stabilized method to study the wind field and compared with experimental and numerical results.Finally,vortex structures for wind flow around a surface-mounted cube are studied by present numerical method.Stable and satisfactory results are obtained,which are consistent with most of the measurements even under coarse mesh.展开更多
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 estab...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.展开更多
A model suitable for describing the mechanical response of thin elastic objects is proposed to simulate the deformation of guide wires in minimally invasive interventions. The main objective of this simulation is to p...A model suitable for describing the mechanical response of thin elastic objects is proposed to simulate the deformation of guide wires in minimally invasive interventions. The main objective of this simulation is to provide doctors an opportunity to rehearse the surgery and select an optimal operation plan before the real surgery. In this model the guide wire is discretized with the multi-body representation and its elastic energy derivate from elastic theory is a polynomial function of the nodal displacements. The vascular structure is represented by a tetrahedron mesh extended from the triangular mesh of the artery, which can be extracted from the patient＇s CT image data. The model applies the energy decline process of the conjugate gradient method to the deformation simulation of the guide wire. Experimental results show that the polynomial relationship between elastic energy and nodal displacements tremendously simplifies the evaluation of the conjugate gradient method and significantly improves the model＇s efficiency. Compared with models depending on an explicit scheme for evaluation, the new model is not only non-conditionally stable but also more efficient. The model can be applied to the real-time simulation of guide wire in a vascular structure.展开更多
Vaporizing foil actuator spot welding method is used in this paper to join magnesium alloy AZ31 and uncoated high-strength steel DP590,which are typically considered as un-weldable due to their high physical property ...Vaporizing foil actuator spot welding method is used in this paper to join magnesium alloy AZ31 and uncoated high-strength steel DP590,which are typically considered as un-weldable due to their high physical property disparities,low mutual solubility,and the lack of any intermetallic phases.Characterization results from scanning electron microscopy(SEM)and high-resolution transmission electron microscopy(HRTEM)of the weld interface indicate that the impact creates an Mg nanocrystalline interlayer with abundant Fe particles.The interlayer exhibits intact bonding with both DP590 and AZ31 substrates.To investigate the fundamental bond formation mechanisms at the interface,a finite element(FE)-based process simulation is first performed to calculate the local temperature and deformation at the interface under the given macroscopic experimental condition.Taking the FE results at the interface as inputs,molecular dynamics(MD)simulations are conducted to study the interlayer formation at the Mg/Fe interface during the impact and cooling.The results found a high velocity shearing-induced mechanical mixing mechanism that mixes Mg/Fe atoms at the interface and creates the interlayer,leading to the metallurgical bond between Mg/steel alloys.展开更多
Cables composed of rare-earth barium copper oxide(REBCO)tapes have been extensively used in various superconducting devices.In recent years,conductor on round core(CORC)cable has drawn the attention of researchers wit...Cables composed of rare-earth barium copper oxide(REBCO)tapes have been extensively used in various superconducting devices.In recent years,conductor on round core(CORC)cable has drawn the attention of researchers with its outstanding current-carrying capacity and mechanical properties.The REBCO tapes are wound spirally on the surface of CORC cable.Under extreme loadings,the REBCO tapes with layered composite structures are vulnerable,which can lead to degradation of critical current and even quenching of superconducting devices.In this paper,we simulate the deformation of CORC cable under external loads,and analyze the damage inside the tape with the cohesive zone model(CZM).Firstly,the fabrication and cabling of CORC are simulated,and the stresses and strains generated in the tape are extracted as the initial condition of the next step.Then,the tension and bending loads are applied to CORC cable,and the damage distribution inside the tape is presented.In addition,the effects of some parameters on the damage are discussed during the bending simulations.展开更多
A method of constructing three-dimensional process model for the punching cartridge cases is presented based on DEFORM simulation analysis. Using DEFORM software,the finite element simulation models for the punching a...A method of constructing three-dimensional process model for the punching cartridge cases is presented based on DEFORM simulation analysis. Using DEFORM software,the finite element simulation models for the punching and forming process of cartridge cases are established,and the corresponding simulation result model of each intermediate procedure is obtained by continuously performing the forming process simulation. The simulation model cannot annotate size and process information due to poor interface between DEFORM software and CAD software. Thus,a 3D annotation module is developed with secondary development technology of UG NX software. Consequently,the final process model with dimension and process information is obtained. Then,with the current 3D process management system,the 3D punching and forming process design of cartridge cases can be completed further. An example is also provided to illustrate that the relative error between the simulation process model and the physical model is less than 2%,which proves the validity and reliability of the proposed method in this study.展开更多
Cased explosives generate highly energetic fragments as their casing breaks up. Due to the complexity of casing fragment related behavior such as embedment, perforation and ricochet, it may be insufficient to use equi...Cased explosives generate highly energetic fragments as their casing breaks up. Due to the complexity of casing fragment related behavior such as embedment, perforation and ricochet, it may be insufficient to use equivalent triangular pressure loading in fragment impact simulations. This simplified method may over- or under-predict the target response. Recently, a procedure using contact techniques has been proposed to overcome such difficulties. It has been shown that the new method has the inherent capability in modeling the multi-piece and multi-hit fragment impact problems in a more realistic way. To investigate the applicability of the proposed method to simulations involving multi-layer penetration, the selected problems of fragment impact on multi-layer targets are described in this paper. It is demonstrated that this method is capable of predicting the complicated multi-layer structural response caused by fragment impact and penetration. Modeling procedures and some technical issues are also discussed.展开更多
Live bone inherently responds to applied mechanical stimulus by altering its internal tissue composition and ultimately biomechanical properties, structure and function. The final formation may structurally appear inf...Live bone inherently responds to applied mechanical stimulus by altering its internal tissue composition and ultimately biomechanical properties, structure and function. The final formation may structurally appear inferior by design but complete by function. To understand the loading response, this paper numerically investigated structural remodeling of mature sheep femur using evolutionary structural optimization method (ESO). Femur images from Computed Tomography scanner were used to determine the elastic modulus variation and subsequently construct finite element model of the femur with stiffest elasticity measured. Major muscle forces on dominant phases of healthy sheep gait were imposed on the femur under static mode. ESO was applied to progressively alter the remodeling of numerically simulated femur from its initial to final design by iteratively removing elements with low strain energy density (SED). The computations were repeated with two different mesh sizes to test the convergence. The elements within the medullary canal had low SEDs and therefore were removed during the optimization. The SEDs in the remaining elements varied with angle around the circumference of the shaft. Those elements with low SED were inefficient in supporting the load and thus fundamentally explained how bone remodels itself with less stiff inferior tissue to meet load demand. This was in line with the Wolff’s law of transformation of bone. Tissue growth and remodeling process was found to shape the sheep femur to a mechanically optimized structure and this was initiated by SED in macro-scale according to traditional principle of Wolff’s law.展开更多
It is well known that a SMPS （switched-mode power supply） is easy to produce strong EMI （electromagnetic interference） and fails in EMC （electromagnetic compatibility） test for its far field radiation exceeds th...It is well known that a SMPS （switched-mode power supply） is easy to produce strong EMI （electromagnetic interference） and fails in EMC （electromagnetic compatibility） test for its far field radiation exceeds the limits between 30-200 MHz. Based on asymmetry line antenna theory, a novel far field CM （common mode） radiation model, including an equivalent driving source, radiation structure and some key influence factors, is identified and built up for a small flyback power supply. Radiation characteristics of this model are predicted by using Ansoft HFSS software and the model effectiveness is verified by experiment. In the end, the radiation role of some key factors, such as the length of output cable, common mode impedance of AC grid, layout of cable and reflected ground, are studied using simulation in detail.展开更多
Fast mechanical switches(FMSs)are critical components of DC circuit breakers(DCCBs),which require the switch action time to break to a sufficient distance within 3 ms in the DC line breaking scenario,while ensuring a ...Fast mechanical switches(FMSs)are critical components of DC circuit breakers(DCCBs),which require the switch action time to break to a sufficient distance within 3 ms in the DC line breaking scenario,while ensuring a long service life.The breaking mechanism significantly affects the current interruption capability of DCCBs.The operation of the repulsion mechanism,along with the morphology of the arc and its transformation within the interrupter chamber,collectively influence the breaking performance of the FMSs.This paper presents a comprehensive analysis of the FMSs,which serves as the pivotal component of controlled resonance combination circuit breakers(CRCBs).This study establishes a multi physics coupling simulation analysis method based on the equivalent circuit of repulsion mechanism discharge,combined with electromagnetic and solid mechanics fields.By constructing a full cycle magnetohydrodynamic particle arc(MHP)model and using a combined simulation of Finite Element joint model(FEJM),the evolution law of arc characteristics during the superimposed current interruption process was systematically explored.The focus was on analysing the physical process of the zero zone of the superimposed arc,the multi physics field coupling relationship of the arc,and the interaction mechanism with external characteristic parameters.Further combining with optimisation design methods,the effectiveness of the model was verified through experiments,FEJM provides comprehensive technical support for effectively reflecting the stress issues of core components during the breaking process of FMS and can provide accurate theoretical references for the optimisation design of mechanical motion components in FMS.It also accurately represents the arc extinguishing process during the breaking of FMS and provides a convenient method for the selection and design of circuit parameters for the entire circuit breaker.展开更多
The constitutive relation of powder material was derived based on the assumption that metal powder is a kind of elasto-plastic material, complying with an elliptical yield criterion. The constitutive integration algor...The constitutive relation of powder material was derived based on the assumption that metal powder is a kind of elasto-plastic material, complying with an elliptical yield criterion. The constitutive integration algorithm was discussed. A way to solve the elastic strain increment in each iteration step during elasto-plastic transition stage was formulated. Different integration method was used for elastic and plastic strain. The relationship between model parameters and relative density was determined through experiments. The model was implemented into user-subroutines of Marc. With the code, computer simulations for compaction process of a balancer were performed. The part is not axisymmetric and requires two lower punches and one upper punch to form. The relative density distributions of two design cases, in which different initial positions of the punches were set, were obtained and compared. The simulation results indicate the influence of punch position and movement on the density distribution of the green compacts.展开更多
基金Supported by National Natural Science Foundation of China(Grant No.51005158)National Science and Technology Major Project of the Ministry of Science and Technology of China(Grant No.2013ZX04008-011-02)
文摘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.
基金Supported by the Strategic Priority Research Program of the Chinese Academy of Sciences under Grant No XDA09020402the National Integrate Circuit Research Program of China under Grant No 2009ZX02023-003+1 种基金the National Natural Science Foundation of China under Grant Nos 61261160500,61376006,61401444 and 61504157the Science and Technology Council of Shanghai under Grant Nos 14DZ2294900,15DZ2270900 and 14ZR1447500
文摘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.
基金Project(51675100)supported by the National Natural Science Foundation of ChinaProject(2016ZX04004008)supported by the National Numerical Control Equipment Major Project of ChinaProject(6902002116)supported by the Foundation of Certain Ministry of China
文摘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.
基金supported by the National Natural Science Foundation of China(Grant Nos.51805064,51701034)the Scientific and Technological Research Program of Chongqing Municipal Education Commission(Grant Nos.KJQN201801137,KJ1600922)+1 种基金the Basic and Advanced Research Project of Chongqing Science and Technology Commission(Grant Nos.cstc2017jcyj AX0062,cstc2018jcyj AX0035)the Chongqing University Key Laboratory of Micro/Nano Materials Engineering and Technology(Grant Nos.KFJJ2003)
文摘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).
文摘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.
文摘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.
文摘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.
基金supported by the National Natural Science Foundation of China(No.52175393)。
文摘The metal cutting process is accompanied by complex stress field,strain field,temperature field.The comprehensive effects of process parameters on chip morphology,cutting force,tool wear and residual stress are complex and inter-connected.Finite element method(FEM)is considered as an effective method to predict process variables and reveal microscopic physical phenomena in the cutting process.Therefore,the finite element(FE)simulation is used to research the conventional and micro scale cutting process,and the differences in the establishment of process variable FE simulation models are distinguished,thereby improving the accuracy of FE simulation.The reliability and effectiveness of FE simulation model largely depend on the accuracy of the simulation method,constitutive model,friction model,damage model in describing mesh element,the dynamic mechanical behavior of materials,the tool-chip-workpiece contact process and the chip formation mechanism.In this paper,the FE models of conventional and micro process variables are comprehensively and up-to-date reviewed for different materials and machining methods.The purpose is to establish a FE model that is more in line with the real cutting conditions,and to provide the possibility for optimizing the cutting process variables.The development direction of FE simulation of metal cutting process is discussed,which provides guidance for future cutting process modeling.
文摘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.
基金Project supported by the National Natural Science Foundation of China(No.51078230)the Research Fund for the Doctoral Program of Higher Education of China(No.200802480056)the Key Project of Fund of Science and Technology Development of Shanghai(No.10JC1407900),China
文摘In this paper,a stabilized finite element technique,actualized by streamline upwind Petrov-Galerkin(SUPG)stabilized method and three-step finite element method(FEM),for large eddy simulation(LES)is developed to predict the wind flow with high Reynolds numbers.Weak form of LES motion equation is combined with the SUPG stabilized term for the spatial finite element discretization.An explicit three-step scheme is implemented for the temporal discretization.For the numerical example of 2D wind flow over a square rib at Re=4.2×105,the Smagorinsky's subgrid-scale(SSGS)model,the DSGS model,and the DSGS model with Cabot near-wall model are applied,and their results are analyzed and compared with experimental results.Furthermore,numerical examples of 3D wind flow around a surface-mounted cube with different Reynolds numbers are performed using DSGS model with Cabot near-wall model based on the present stabilized method to study the wind field and compared with experimental and numerical results.Finally,vortex structures for wind flow around a surface-mounted cube are studied by present numerical method.Stable and satisfactory results are obtained,which are consistent with most of the measurements even under coarse mesh.
基金Projects (11202125, 61175038) supported by the National Natural Science Foundation of China
文摘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.
文摘A model suitable for describing the mechanical response of thin elastic objects is proposed to simulate the deformation of guide wires in minimally invasive interventions. The main objective of this simulation is to provide doctors an opportunity to rehearse the surgery and select an optimal operation plan before the real surgery. In this model the guide wire is discretized with the multi-body representation and its elastic energy derivate from elastic theory is a polynomial function of the nodal displacements. The vascular structure is represented by a tetrahedron mesh extended from the triangular mesh of the artery, which can be extracted from the patient＇s CT image data. The model applies the energy decline process of the conjugate gradient method to the deformation simulation of the guide wire. Experimental results show that the polynomial relationship between elastic energy and nodal displacements tremendously simplifies the evaluation of the conjugate gradient method and significantly improves the model＇s efficiency. Compared with models depending on an explicit scheme for evaluation, the new model is not only non-conditionally stable but also more efficient. The model can be applied to the real-time simulation of guide wire in a vascular structure.
基金sponsored by the US Department of Energy,Office of Vehicle Technology,under a prime contract with Oak Ridge National Laboratory(ORNL)。
文摘Vaporizing foil actuator spot welding method is used in this paper to join magnesium alloy AZ31 and uncoated high-strength steel DP590,which are typically considered as un-weldable due to their high physical property disparities,low mutual solubility,and the lack of any intermetallic phases.Characterization results from scanning electron microscopy(SEM)and high-resolution transmission electron microscopy(HRTEM)of the weld interface indicate that the impact creates an Mg nanocrystalline interlayer with abundant Fe particles.The interlayer exhibits intact bonding with both DP590 and AZ31 substrates.To investigate the fundamental bond formation mechanisms at the interface,a finite element(FE)-based process simulation is first performed to calculate the local temperature and deformation at the interface under the given macroscopic experimental condition.Taking the FE results at the interface as inputs,molecular dynamics(MD)simulations are conducted to study the interlayer formation at the Mg/Fe interface during the impact and cooling.The results found a high velocity shearing-induced mechanical mixing mechanism that mixes Mg/Fe atoms at the interface and creates the interlayer,leading to the metallurgical bond between Mg/steel alloys.
基金Project supported by the National Natural Science Foundation of China(Nos.U2241267,1217215511872195)。
文摘Cables composed of rare-earth barium copper oxide(REBCO)tapes have been extensively used in various superconducting devices.In recent years,conductor on round core(CORC)cable has drawn the attention of researchers with its outstanding current-carrying capacity and mechanical properties.The REBCO tapes are wound spirally on the surface of CORC cable.Under extreme loadings,the REBCO tapes with layered composite structures are vulnerable,which can lead to degradation of critical current and even quenching of superconducting devices.In this paper,we simulate the deformation of CORC cable under external loads,and analyze the damage inside the tape with the cohesive zone model(CZM).Firstly,the fabrication and cabling of CORC are simulated,and the stresses and strains generated in the tape are extracted as the initial condition of the next step.Then,the tension and bending loads are applied to CORC cable,and the damage distribution inside the tape is presented.In addition,the effects of some parameters on the damage are discussed during the bending simulations.
基金Supported by the National Defense Basic Scientific Research Project(A1020131011)
文摘A method of constructing three-dimensional process model for the punching cartridge cases is presented based on DEFORM simulation analysis. Using DEFORM software,the finite element simulation models for the punching and forming process of cartridge cases are established,and the corresponding simulation result model of each intermediate procedure is obtained by continuously performing the forming process simulation. The simulation model cannot annotate size and process information due to poor interface between DEFORM software and CAD software. Thus,a 3D annotation module is developed with secondary development technology of UG NX software. Consequently,the final process model with dimension and process information is obtained. Then,with the current 3D process management system,the 3D punching and forming process design of cartridge cases can be completed further. An example is also provided to illustrate that the relative error between the simulation process model and the physical model is less than 2%,which proves the validity and reliability of the proposed method in this study.
文摘Cased explosives generate highly energetic fragments as their casing breaks up. Due to the complexity of casing fragment related behavior such as embedment, perforation and ricochet, it may be insufficient to use equivalent triangular pressure loading in fragment impact simulations. This simplified method may over- or under-predict the target response. Recently, a procedure using contact techniques has been proposed to overcome such difficulties. It has been shown that the new method has the inherent capability in modeling the multi-piece and multi-hit fragment impact problems in a more realistic way. To investigate the applicability of the proposed method to simulations involving multi-layer penetration, the selected problems of fragment impact on multi-layer targets are described in this paper. It is demonstrated that this method is capable of predicting the complicated multi-layer structural response caused by fragment impact and penetration. Modeling procedures and some technical issues are also discussed.
文摘Live bone inherently responds to applied mechanical stimulus by altering its internal tissue composition and ultimately biomechanical properties, structure and function. The final formation may structurally appear inferior by design but complete by function. To understand the loading response, this paper numerically investigated structural remodeling of mature sheep femur using evolutionary structural optimization method (ESO). Femur images from Computed Tomography scanner were used to determine the elastic modulus variation and subsequently construct finite element model of the femur with stiffest elasticity measured. Major muscle forces on dominant phases of healthy sheep gait were imposed on the femur under static mode. ESO was applied to progressively alter the remodeling of numerically simulated femur from its initial to final design by iteratively removing elements with low strain energy density (SED). The computations were repeated with two different mesh sizes to test the convergence. The elements within the medullary canal had low SEDs and therefore were removed during the optimization. The SEDs in the remaining elements varied with angle around the circumference of the shaft. Those elements with low SED were inefficient in supporting the load and thus fundamentally explained how bone remodels itself with less stiff inferior tissue to meet load demand. This was in line with the Wolff’s law of transformation of bone. Tissue growth and remodeling process was found to shape the sheep femur to a mechanically optimized structure and this was initiated by SED in macro-scale according to traditional principle of Wolff’s law.
文摘It is well known that a SMPS （switched-mode power supply） is easy to produce strong EMI （electromagnetic interference） and fails in EMC （electromagnetic compatibility） test for its far field radiation exceeds the limits between 30-200 MHz. Based on asymmetry line antenna theory, a novel far field CM （common mode） radiation model, including an equivalent driving source, radiation structure and some key influence factors, is identified and built up for a small flyback power supply. Radiation characteristics of this model are predicted by using Ansoft HFSS software and the model effectiveness is verified by experiment. In the end, the radiation role of some key factors, such as the length of output cable, common mode impedance of AC grid, layout of cable and reflected ground, are studied using simulation in detail.
基金State Grid Corporation Headquarters Science and Technology Project,Grant/Award Number:5500-20220110A-1-1-ZN。
文摘Fast mechanical switches(FMSs)are critical components of DC circuit breakers(DCCBs),which require the switch action time to break to a sufficient distance within 3 ms in the DC line breaking scenario,while ensuring a long service life.The breaking mechanism significantly affects the current interruption capability of DCCBs.The operation of the repulsion mechanism,along with the morphology of the arc and its transformation within the interrupter chamber,collectively influence the breaking performance of the FMSs.This paper presents a comprehensive analysis of the FMSs,which serves as the pivotal component of controlled resonance combination circuit breakers(CRCBs).This study establishes a multi physics coupling simulation analysis method based on the equivalent circuit of repulsion mechanism discharge,combined with electromagnetic and solid mechanics fields.By constructing a full cycle magnetohydrodynamic particle arc(MHP)model and using a combined simulation of Finite Element joint model(FEJM),the evolution law of arc characteristics during the superimposed current interruption process was systematically explored.The focus was on analysing the physical process of the zero zone of the superimposed arc,the multi physics field coupling relationship of the arc,and the interaction mechanism with external characteristic parameters.Further combining with optimisation design methods,the effectiveness of the model was verified through experiments,FEJM provides comprehensive technical support for effectively reflecting the stress issues of core components during the breaking process of FMS and can provide accurate theoretical references for the optimisation design of mechanical motion components in FMS.It also accurately represents the arc extinguishing process during the breaking of FMS and provides a convenient method for the selection and design of circuit parameters for the entire circuit breaker.
基金Porject(50325516) supported by the National Natural Science Foundation of China Project(CG2003-GA005) supported by China Education and Research Grid (China Grid) Project(003019) supported by the Natural Science Foundation of Guangdong Province, China
文摘The constitutive relation of powder material was derived based on the assumption that metal powder is a kind of elasto-plastic material, complying with an elliptical yield criterion. The constitutive integration algorithm was discussed. A way to solve the elastic strain increment in each iteration step during elasto-plastic transition stage was formulated. Different integration method was used for elastic and plastic strain. The relationship between model parameters and relative density was determined through experiments. The model was implemented into user-subroutines of Marc. With the code, computer simulations for compaction process of a balancer were performed. The part is not axisymmetric and requires two lower punches and one upper punch to form. The relative density distributions of two design cases, in which different initial positions of the punches were set, were obtained and compared. The simulation results indicate the influence of punch position and movement on the density distribution of the green compacts.
