A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic gr...A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic growth was simulated from undercooled nickel melt under the forced flow. The simulation results show that the asymmetry behavior of the dendritic growth is caused by the forced flow. When the flow velocity is less than the critical value, the asymmetry of dendrite is little influenced by the forced flow. Once the flow velocity reaches or exceeds the critical value, the controlling factor of dendrite growth gradually changes from thermal diffusion to convection. With the increase of the flow velocity, the deflection angle towards upstream direction of the primary dendrite stem becomes larger. The effect of the dendrite growth on the flow field of the melt is apparent. With the increase of the dendrite size, the vortex is present in the downstream regions, and the vortex region is gradually enlarged. Dendrite tips appear to remelt. In addition, the adaptive finite element method can reduce CPU running time by one order of magnitude compared with uniform grid method, and the speed-up ratio is proportional to the size of computational domain.展开更多
Advanced fiber reinforced polymer composites have been increasingly applied to various structural components. One of the important processes to fabricate high performance laminated composites is an autoclave assisted ...Advanced fiber reinforced polymer composites have been increasingly applied to various structural components. One of the important processes to fabricate high performance laminated composites is an autoclave assisted prepreg lay-up. Since the quality of laminated composites is largely affected by the cure cycle, selection of an appropriate cure cycle for each application is important and must be optimized. Thus, some fundamental model of the consolidation and cure processes is necessary for selecting suitable parameters for a specific application. This article is concerned with the ＂flow-compaction＂ model during the autoclave processing of composite materials. By using a weighted residual method, two-dimensional finite element formulation for the consolidation process of thick thermosetting composites is presented and the corresponding finite element code is developed. Numerical examples, including comparison of the present numerical results with one-dimensional and twodimensional analytical solutions, are given to illustrate the accuracy and effectiveness of the proposed finite element formulation. In addition, a consolidation simulation of AS4/3501-6 graphite/epoxy laminate is carded out and compared with the experimental results available in the literature.展开更多
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.展开更多
A mixed finite element formulation for viscoelastic flows is derived in this paper, in which the FIC （finite incremental calculus） pressure stabilization process and the DEVSS （discrete elastic viscous stress split...A mixed finite element formulation for viscoelastic flows is derived in this paper, in which the FIC （finite incremental calculus） pressure stabilization process and the DEVSS （discrete elastic viscous stress splitting） method using the Crank-Nicolson-based split are introduced within a general framework of the iterative version of the fractional step algorithm. The SU （streamline-upwind） method is particularly chosen to tackle the convective terms in constitutive equations of viscoelastic flows. Thanks to the proposed scheme the finite elements with equal low-order interpolation approximations for stress-velocity-pressure variables can be successfully used even for viscoelastic flows with high Weissenberg numbers. The XPP （extended Pom-Pom） constitutive model for describing viscoelastic behaviors is particularly integrated into the proposed scheme. The numerical results for the 4：1 sudden contraction flow problem demonstrate prominent stability, accuracy and convergence rate of the proposed scheme in both pressure and stress distributions over the flow domain within a wide range of the Weissenberg number, particularly the capability in reproducing the results, which can be used to explain the ＂die swell＂ phenomenon observed in the polymer injection molding process.展开更多
<span style="font-family:Verdana;">The present study deals with the unsteady dynamics of cavitation around the NACA 0015 hydrofoil in a channel. A finite element model is proposed to solve the governin...<span style="font-family:Verdana;">The present study deals with the unsteady dynamics of cavitation around the NACA 0015 hydrofoil in a channel. A finite element model is proposed to solve the governing equations of momentum and mass conservation. Turbulent flows around the hydrofoil are described by the Prandtl-Kolmogorov model. The cavitation phenomenon is modeled through a mixture model involving liquid and vapor flows and the Zwart-Gerber-Belamri (ZGB) model is considered to evaluate the transport of the water vapor fraction. The variational finite element model formulation includes the mixing of the characteristic method and the finite element. Also, at the open sides of the channel flow, an open boundary condition is imposed. Numerical experiments are performed for cavitation numbers 0.8 and 0.4. The presented model predicts the essential features of unsteady cavitating flows, the generation of vapor cavities, the time-dependent oscillations of the variables and the presence of vortical flow structures associated to vapor volume concentrations during the shedding process.</span>展开更多
Numerical modelling of geological processes, such as mantle convection, flow in porous media, and geothermal heat transfer, has become quite common with the increase in computing and the availability of usable softwar...Numerical modelling of geological processes, such as mantle convection, flow in porous media, and geothermal heat transfer, has become quite common with the increase in computing and the availability of usable software. Today modelling these dynamical processes entails the solving of the governing equations involving the mass, momentum, energy and chemical transport. These equations represent partial differential equations and must be solved on powerful enough computers because they require sufficient spatial and temporal resolution to be useful. We describe here the salient and outstanding features of the SEPRAN software package, developed in the Netherlands, as a case study for a robust and user-friendly soft- ware, which the geological community can utilize in handling many thermal-mechanical-chemical problems found in geology, which will include geothermal situations, where many types of partial differential equations must be solved at the same time with thermodynamical input parameters.展开更多
The comprehensive tire building and shaping processes are investigated through the finite element method(FEM)in this article.The mechanical properties of the uncured rubber from different tire components are investiga...The comprehensive tire building and shaping processes are investigated through the finite element method(FEM)in this article.The mechanical properties of the uncured rubber from different tire components are investigated through cyclic loading-unloading experiments under different strain rates.Based on the experiments,an elastoviscoplastic constitutive model is adopted to describe themechanical behaviors of the uncured rubber.The distinct mechanical properties,including the stress level,hysteresis and residual strain,of the uncured rubber can all be well characterized.The whole tire building process(including component winding,rubber bladder inflation,component stitching and carcass band folding-back)and the shaping process are simulated using this constitutive model.The simulated green tire profile is in good agreement with the actual profile obtained through 3D scanning.The deformation and stress of the rubber components and the cord reinforcements during production can be obtained fromthe FE simulation,which is helpful for judging the rationality of the tire construction design.Finally,the influence of the parameter“drum width”is investigated,and the simulated result is found to be consistent with the experimental observations,which verifies the effectiveness of the simulation.The established simulation strategy provides some guiding significance for the improvement of tire design parameters and the elimination of tire production defects.展开更多
The powder compaction simulations were performed to demonstrate deformation behavior of particles and estimate the effect of different punch speeds and particle diameters on the relative density of powder by a multi-p...The powder compaction simulations were performed to demonstrate deformation behavior of particles and estimate the effect of different punch speeds and particle diameters on the relative density of powder by a multi-particle finite element model(MPFEM). Individual particle discretized with a finite element mesh allows for a full description of the contact mechanics. In order to verify the reliability of compaction simulation by MPFEM, the compaction tests of porous aluminum with average particle size of 20 μm and 3 μm were performed at different ram speeds of 5, 15, 30 and 60 mm/min by MTS servo-hydraulic tester. The results show that the slow ram speed is of great advantage for powder densification in low compaction force due to sufficient particle rearrangement and compaction force increases with decrease in average particle size of aluminum.展开更多
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.展开更多
This paper models hot compression testing using a dilatometer in loading mode.These small-scale tests provide a high throughput at low cost,but are susceptible to inhomogeneity due to friction and temperature gradient...This paper models hot compression testing using a dilatometer in loading mode.These small-scale tests provide a high throughput at low cost,but are susceptible to inhomogeneity due to friction and temperature gradients.A novel method is presented for correcting the true stress-strain constitutive response over the full range of temperatures,strain-rates and strain.The nominal response from the tests is used to predict the offset in the stress-strain curves due to inhomogeneity,and this stress offsetΔσis applied piecewise to the data,correcting the constitutive response in one iteration.A key new feature is the smoothing and fitting of the flow stress data as a function of temperature and strain-rate,at multiple discrete strains.The corrected model then provides quantitative prediction of the spatial and temporal variation in strain-rate and strain throughout the sample,needed to correlate the local deformation conditions with the microstructure and texture evolution.The study uses a detailed series of 144 hot compression tests of a Zr-Nb alloy.While this is an important wrought nuclear alloy in its own right,it also serves here as a test case for modelling the dilatometer for hot testing of high temperature alloys,particularly those with dualα-βphase microstructures(such as titanium alloys).展开更多
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.展开更多
The occlusal design plays a decisive role in the fabrication of dental restorations.Dentists and dental technicians depend on mechanical simulations of mandibular movement that are as accurate as possible,in particula...The occlusal design plays a decisive role in the fabrication of dental restorations.Dentists and dental technicians depend on mechanical simulations of mandibular movement that are as accurate as possible,in particular,to produce interference-free yet chewing-efficient dental restorations.For this,kinetic data must be available,i.e.,movements and deformations under the influence of forces and stresses.In the present study,so-called functional data were collected from healthy volunteers to provide consistent information for proper kinetics.For the latter purpose,biting and chewing forces,electrical muscle activity and jaw movements were registered synchronously,and individual magnetic resonance tomograms(MRI)were prepared.The acquired data were then added to a large complex finite element model of the complete masticatory system using the functional information obtained and individual anatomical geometries so that the kinetics of the chewing process and teeth grinding could be realistically simulated.This allows developing algorithms that optimize computer-aided manufacturing of dental prostheses close to occlusion.In this way,a failure-free function of the dental prosthesis can be guaranteed and its damage during usage can be reduced or prevented even including endosseous implants.展开更多
Compaction processes are one the most important par ts of powder forming technology. The main applications are focused on pieces for a utomotive, aeronautic, electric and electronic industries. The main goals of the c...Compaction processes are one the most important par ts of powder forming technology. The main applications are focused on pieces for a utomotive, aeronautic, electric and electronic industries. The main goals of the compaction processes are to obtain a compact with the geometrical requirements, without cracks, and with a uniform distribution of density. Design of such proc esses consist, essentially, in determine the sequence and relative displacements of die and punches in order to achieve such goals. A.B. Khoei presented a gener al framework for the finite element simulation of powder forming processes based on the following aspects; a large displacement formulation, centred on a total and updated Lagrangian formulation; an adaptive finite element strategy based on error estimates and automatic remeshing techniques; a cap model based on a hard ening rule in modelling of the highly non-linear behaviour of material; and the use of an efficient contact algorithm in the context of an interface element fo rmulation. In these references, the non-linear behaviour of powder was adequately desc ribed by the cap plasticity model. However, it suffers from a serious deficiency when the stress-point reaches a yield surface. In the flow theory of plasticit y, the transition from an elastic state to an elasto-plastic state appears more or less abruptly. For powder material it is very difficult to define the locati on of yield surface, because there is no distinct transition from elastic to ela stic-plastic behaviour. Results of experimental test on some hard met al powder show that the plastic effects were begun immediately upon loading. In such mater ials the domain of the yield surface would collapse to a point, so making the di rection of plastic increment indeterminate, because all directions are normal to a point. Thus, the classical plasticity theory cannot deal with such materials and an advanced constitutive theory is necessary. In the present paper, the constitutive equations of powder materials will be discussed via an endochronic theory of plasticity. This theory provides a unifi ed point of view to describe the elastic-plastic behaviour of material since it places no requirement for a yield surface and a ’loading function’ to disting uish between loading an unloading. Endochronic theory of plasticity has been app lied to a number of metallic materials, concrete and sand, but to the knowledge of authors, no numerical scheme of the model has been applied to powder material . In the present paper, a new approach is developed based on an endochronic rate independent, density-dependent plasticity model for describing the isothermal deformation behavior of metal powder at low homologous temperature. Although the concept of yield surface has not been explicitly assumed in endochronic theory, it is shown that the cone-cap plasticity yield surface (Fig.1), which is the m ost commonly used plasticity models for describing the behavior of powder materi al can be easily derived as a special case of the proposed endochronic theory. Fig.1 Trace of cone-cap yield function on the meridian pl ane for different relative density As large deformation is observed in powder compaction process, a hypoelastic-pl astic formulation is developed in the context of finite deformation plasticity. Constitutive equations are stated in unrotated frame of reference that greatly s implifies endochronic constitutive relation in finite plasticity. Constitutive e quations of the endochronic theory and their numerical integration are establish ed and procedures for determining material parameters of the model are demonstra ted. Finally, the numerical schemes are examined for efficiency in the model ling of a tip shaped component, as shown in Fig.2. Fig.2 A shaped tip component. a) Geometry, boundary conditio n and finite element mesh; b) density distribution at final stage of展开更多
The images of granular ore media were captured by X-ray CT scanner. Combined with digital image processing and finite element techniques, the three-dimensional geometrical model, which represents the realistic pore st...The images of granular ore media were captured by X-ray CT scanner. Combined with digital image processing and finite element techniques, the three-dimensional geometrical model, which represents the realistic pore structure of the media, was constructed. With this model, three dimensional pore scale fluid flow among particles was simulated. Then the distributions of fluid flow velocity and pressure were analyzed and the hydraulic conductivity was calculated. The simulation results indicate the fluid flow behaviors are mainly dominated by the volume and topological structure of pore space. There exist obvious preferential flow and leaching blind zones simultaneously in the medium. The highest velocities generally occur in those narrow pores with high pressure drops. The hydraulic conductivity obtained by simulation is the same order of magnitude as the laboratory test result, which denotes the validity of the model. The pore-scale and macro-scale are combined and the established geometrical model can be used for the simulations of other phenomena during heap leaching process.展开更多
The hot deformation behavior of extruded AZ80 magnesium alloy was investigated using compression tests in the temperature range of 250–400 °C and strain rate range of 0.001–1.000 s–1. The 3D power dissipation ...The hot deformation behavior of extruded AZ80 magnesium alloy was investigated using compression tests in the temperature range of 250–400 °C and strain rate range of 0.001–1.000 s–1. The 3D power dissipation map was developed to evaluate the hot deformation mechanisms and determine the optimal processing parameters. Two domains of dynamic recrystallization were identified from the 3D power dissipation map, with one occurring in the temperature and strain rate range of 250–320 °C and 0.001–0.010 s–1and the other one occurring in the temperature and strain rate range of 380–400 °C and 0.001–0.003 s–1. In order to delineate the regions of flow instability, Prasad’s instability criterion, Murty’s instability criterion and Gegel’s stability criteria were employed to develop the 3D instability maps. Through microstructural examination, it is found that Prasad’s and Murty’s instability criteria are more effective than Gegel’s stability criteria in predicting the flow instability regions for extruded AZ80 alloy. Further, the 3D processing maps were integrated into finite element simulation and the predictions of the simulation are in good agreement with the experimental results.展开更多
A new coupled finite element formulation is proposed to calculate a conformation tensor model in two complex flows： a planar contraction flow and a planar flow around a symmetrically placed cylinder. The components o...A new coupled finite element formulation is proposed to calculate a conformation tensor model in two complex flows： a planar contraction flow and a planar flow around a symmetrically placed cylinder. The components of conformation tensor are first computed together with the velocity and pressure to describe the change of morphology of polymer chain coils in flow fields. Macroscopic quantities of viscoelastic flow are then calculated based on the conformation tensor. Comparisons between the numerical simulations and experiments for stress patterns and velocity profiles are carried out to prove the validity of the method.展开更多
In this paper, a unified internal state variable(ISV) model for predicting microstructure evolution during hot working process of AZ80 magnesium alloy was developed. A novel aspect of the proposed model is that the in...In this paper, a unified internal state variable(ISV) model for predicting microstructure evolution during hot working process of AZ80 magnesium alloy was developed. A novel aspect of the proposed model is that the interactive effects of material hardening, recovery and dynamic recrystallization(DRX) on the characteristic deformation behavior were considered by incorporating the evolution laws of viscoplastic flow, dislocation activities, DRX nucleation and boundary migration in a coupled manner. The model parameters were calibrated based on the experimental data analysis and genetic algorithm(GA) based objective optimization. The predicted flow stress, DRX fraction and average grain size match well with experimental results. The proposed model was embedded in the finite element(FE) software DEFORM-3 D via user defined subroutine to simulate the hot compression and equal channel angular extrusion(ECAE) processes. The heterogeneous microstructure distributions at different deformation zones and the dislocation density evolution with competitive deformation mechanisms were captured.This study can provide a theoretical solution for the hot working problems of magnesium alloy.展开更多
An emerging real-time ground compaction and quality control, known as intelligent compaction(IC), has been applied for efficiently optimising the full-area compaction. Although IC technology can provide real-time asse...An emerging real-time ground compaction and quality control, known as intelligent compaction(IC), has been applied for efficiently optimising the full-area compaction. Although IC technology can provide real-time assessment of uniformity of the compacted area, accurate determination of the soil stiffness required for quality control and design remains challenging. In this paper, a novel and advanced numerical model simulating the interaction of vibratory drum and soil beneath is developed. The model is capable of evaluating the nonlinear behaviour of underlying soil subjected to dynamic loading by capturing the variations of damping with the cyclic shear strains and degradation of soil modulus. The interaction of the drum and the soil is simulated via the finite element method to develop a comprehensive dataset capturing the dynamic responses of the drum and the soil. Indeed, more than a thousand three-dimensional(3D) numerical models covering various soil characteristics, roller weights, vibration amplitudes and frequencies were adopted. The developed dataset is then used to train the inverse solver using an innovative machine learning approach, i.e. the extended support vector regression, to simulate the stiffness of the compacted soil by adopting drum acceleration records. Furthermore, the impacts of the amplitude and frequency of the vibration on the level of underlying soil compaction are discussed.The proposed machine learning approach is promising for real-time extraction of actual soil stiffness during compaction. Results of the study can be employed by practising engineers to interpret roller drum acceleration data to estimate the level of compaction and ground stiffness during compaction.展开更多
基金Projects(51161011,11364024)supported by the National Natural Science Foundation of ChinaProject(1204GKCA065)supported by the Key Technology R&D Program of Gansu Province,China+1 种基金Project(201210)supported by the Fundamental Research Funds for the Universities of Gansu Province,ChinaProject(J201304)supported by the Funds for Distinguished Young Scientists of Lanzhou University of Technology,China
文摘A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic growth was simulated from undercooled nickel melt under the forced flow. The simulation results show that the asymmetry behavior of the dendritic growth is caused by the forced flow. When the flow velocity is less than the critical value, the asymmetry of dendrite is little influenced by the forced flow. Once the flow velocity reaches or exceeds the critical value, the controlling factor of dendrite growth gradually changes from thermal diffusion to convection. With the increase of the flow velocity, the deflection angle towards upstream direction of the primary dendrite stem becomes larger. The effect of the dendrite growth on the flow field of the melt is apparent. With the increase of the dendrite size, the vortex is present in the downstream regions, and the vortex region is gradually enlarged. Dendrite tips appear to remelt. In addition, the adaptive finite element method can reduce CPU running time by one order of magnitude compared with uniform grid method, and the speed-up ratio is proportional to the size of computational domain.
基金The project supported by the National Natural Science Foundation of China (10272037)The English text was polished by Yunming Chen.
文摘Advanced fiber reinforced polymer composites have been increasingly applied to various structural components. One of the important processes to fabricate high performance laminated composites is an autoclave assisted prepreg lay-up. Since the quality of laminated composites is largely affected by the cure cycle, selection of an appropriate cure cycle for each application is important and must be optimized. Thus, some fundamental model of the consolidation and cure processes is necessary for selecting suitable parameters for a specific application. This article is concerned with the ＂flow-compaction＂ model during the autoclave processing of composite materials. By using a weighted residual method, two-dimensional finite element formulation for the consolidation process of thick thermosetting composites is presented and the corresponding finite element code is developed. Numerical examples, including comparison of the present numerical results with one-dimensional and twodimensional analytical solutions, are given to illustrate the accuracy and effectiveness of the proposed finite element formulation. In addition, a consolidation simulation of AS4/3501-6 graphite/epoxy laminate is carded out and compared with the experimental results available in the literature.
基金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 National Natural Science Foundation of China (10672033,10590354,90715011 and 10272027)the National Key Basic Research and Development Program (2002CB412709)
文摘A mixed finite element formulation for viscoelastic flows is derived in this paper, in which the FIC （finite incremental calculus） pressure stabilization process and the DEVSS （discrete elastic viscous stress splitting） method using the Crank-Nicolson-based split are introduced within a general framework of the iterative version of the fractional step algorithm. The SU （streamline-upwind） method is particularly chosen to tackle the convective terms in constitutive equations of viscoelastic flows. Thanks to the proposed scheme the finite elements with equal low-order interpolation approximations for stress-velocity-pressure variables can be successfully used even for viscoelastic flows with high Weissenberg numbers. The XPP （extended Pom-Pom） constitutive model for describing viscoelastic behaviors is particularly integrated into the proposed scheme. The numerical results for the 4：1 sudden contraction flow problem demonstrate prominent stability, accuracy and convergence rate of the proposed scheme in both pressure and stress distributions over the flow domain within a wide range of the Weissenberg number, particularly the capability in reproducing the results, which can be used to explain the ＂die swell＂ phenomenon observed in the polymer injection molding process.
文摘<span style="font-family:Verdana;">The present study deals with the unsteady dynamics of cavitation around the NACA 0015 hydrofoil in a channel. A finite element model is proposed to solve the governing equations of momentum and mass conservation. Turbulent flows around the hydrofoil are described by the Prandtl-Kolmogorov model. The cavitation phenomenon is modeled through a mixture model involving liquid and vapor flows and the Zwart-Gerber-Belamri (ZGB) model is considered to evaluate the transport of the water vapor fraction. The variational finite element model formulation includes the mixing of the characteristic method and the finite element. Also, at the open sides of the channel flow, an open boundary condition is imposed. Numerical experiments are performed for cavitation numbers 0.8 and 0.4. The presented model predicts the essential features of unsteady cavitating flows, the generation of vapor cavities, the time-dependent oscillations of the variables and the presence of vortical flow structures associated to vapor volume concentrations during the shedding process.</span>
基金CMG and Geochemistry programs of the National Science foundationThe Netherlands Research Center for Integrated Solid Earth Science (ISES) project ME-2.7
文摘Numerical modelling of geological processes, such as mantle convection, flow in porous media, and geothermal heat transfer, has become quite common with the increase in computing and the availability of usable software. Today modelling these dynamical processes entails the solving of the governing equations involving the mass, momentum, energy and chemical transport. These equations represent partial differential equations and must be solved on powerful enough computers because they require sufficient spatial and temporal resolution to be useful. We describe here the salient and outstanding features of the SEPRAN software package, developed in the Netherlands, as a case study for a robust and user-friendly soft- ware, which the geological community can utilize in handling many thermal-mechanical-chemical problems found in geology, which will include geothermal situations, where many types of partial differential equations must be solved at the same time with thermodynamical input parameters.
基金funded by the NationalNatural Science Foundation of China (Nos.11902229,11502181)the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant Nos.XDB22040502,XDC06030200).
文摘The comprehensive tire building and shaping processes are investigated through the finite element method(FEM)in this article.The mechanical properties of the uncured rubber from different tire components are investigated through cyclic loading-unloading experiments under different strain rates.Based on the experiments,an elastoviscoplastic constitutive model is adopted to describe themechanical behaviors of the uncured rubber.The distinct mechanical properties,including the stress level,hysteresis and residual strain,of the uncured rubber can all be well characterized.The whole tire building process(including component winding,rubber bladder inflation,component stitching and carcass band folding-back)and the shaping process are simulated using this constitutive model.The simulated green tire profile is in good agreement with the actual profile obtained through 3D scanning.The deformation and stress of the rubber components and the cord reinforcements during production can be obtained fromthe FE simulation,which is helpful for judging the rationality of the tire construction design.Finally,the influence of the parameter“drum width”is investigated,and the simulated result is found to be consistent with the experimental observations,which verifies the effectiveness of the simulation.The established simulation strategy provides some guiding significance for the improvement of tire design parameters and the elimination of tire production defects.
基金supported by a grant-in-aid for the National Core Research Center Program from the Ministry of Education Science & Technology,Koreathe Korea Science & Engineering Foundation (No.R15-2006-022-03003-0)
文摘The powder compaction simulations were performed to demonstrate deformation behavior of particles and estimate the effect of different punch speeds and particle diameters on the relative density of powder by a multi-particle finite element model(MPFEM). Individual particle discretized with a finite element mesh allows for a full description of the contact mechanics. In order to verify the reliability of compaction simulation by MPFEM, the compaction tests of porous aluminum with average particle size of 20 μm and 3 μm were performed at different ram speeds of 5, 15, 30 and 60 mm/min by MTS servo-hydraulic tester. The results show that the slow ram speed is of great advantage for powder densification in low compaction force due to sufficient particle rearrangement and compaction force increases with decrease in average particle size of aluminum.
基金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.
基金Light Form,a UK Engineering and Physical Sciences Research Council(EPSRC)programme grant(EP/R001715/1)。
文摘This paper models hot compression testing using a dilatometer in loading mode.These small-scale tests provide a high throughput at low cost,but are susceptible to inhomogeneity due to friction and temperature gradients.A novel method is presented for correcting the true stress-strain constitutive response over the full range of temperatures,strain-rates and strain.The nominal response from the tests is used to predict the offset in the stress-strain curves due to inhomogeneity,and this stress offsetΔσis applied piecewise to the data,correcting the constitutive response in one iteration.A key new feature is the smoothing and fitting of the flow stress data as a function of temperature and strain-rate,at multiple discrete strains.The corrected model then provides quantitative prediction of the spatial and temporal variation in strain-rate and strain throughout the sample,needed to correlate the local deformation conditions with the microstructure and texture evolution.The study uses a detailed series of 144 hot compression tests of a Zr-Nb alloy.While this is an important wrought nuclear alloy in its own right,it also serves here as a test case for modelling the dilatometer for hot testing of high temperature alloys,particularly those with dualα-βphase microstructures(such as titanium alloys).
基金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.
基金We acknowledge the support of the German Research Foundation Grant Nos.SCHM 2456/5-1 and SCHW 307/30-1together with funding for the project initial phase from the German Federal Ministry for Economy and Technology Grant No.KF 2875101WM.(Bundesministerium für Wirtschaft und Technologie)according to a decision of the German Bundestag.
文摘The occlusal design plays a decisive role in the fabrication of dental restorations.Dentists and dental technicians depend on mechanical simulations of mandibular movement that are as accurate as possible,in particular,to produce interference-free yet chewing-efficient dental restorations.For this,kinetic data must be available,i.e.,movements and deformations under the influence of forces and stresses.In the present study,so-called functional data were collected from healthy volunteers to provide consistent information for proper kinetics.For the latter purpose,biting and chewing forces,electrical muscle activity and jaw movements were registered synchronously,and individual magnetic resonance tomograms(MRI)were prepared.The acquired data were then added to a large complex finite element model of the complete masticatory system using the functional information obtained and individual anatomical geometries so that the kinetics of the chewing process and teeth grinding could be realistically simulated.This allows developing algorithms that optimize computer-aided manufacturing of dental prostheses close to occlusion.In this way,a failure-free function of the dental prosthesis can be guaranteed and its damage during usage can be reduced or prevented even including endosseous implants.
文摘Compaction processes are one the most important par ts of powder forming technology. The main applications are focused on pieces for a utomotive, aeronautic, electric and electronic industries. The main goals of the compaction processes are to obtain a compact with the geometrical requirements, without cracks, and with a uniform distribution of density. Design of such proc esses consist, essentially, in determine the sequence and relative displacements of die and punches in order to achieve such goals. A.B. Khoei presented a gener al framework for the finite element simulation of powder forming processes based on the following aspects; a large displacement formulation, centred on a total and updated Lagrangian formulation; an adaptive finite element strategy based on error estimates and automatic remeshing techniques; a cap model based on a hard ening rule in modelling of the highly non-linear behaviour of material; and the use of an efficient contact algorithm in the context of an interface element fo rmulation. In these references, the non-linear behaviour of powder was adequately desc ribed by the cap plasticity model. However, it suffers from a serious deficiency when the stress-point reaches a yield surface. In the flow theory of plasticit y, the transition from an elastic state to an elasto-plastic state appears more or less abruptly. For powder material it is very difficult to define the locati on of yield surface, because there is no distinct transition from elastic to ela stic-plastic behaviour. Results of experimental test on some hard met al powder show that the plastic effects were begun immediately upon loading. In such mater ials the domain of the yield surface would collapse to a point, so making the di rection of plastic increment indeterminate, because all directions are normal to a point. Thus, the classical plasticity theory cannot deal with such materials and an advanced constitutive theory is necessary. In the present paper, the constitutive equations of powder materials will be discussed via an endochronic theory of plasticity. This theory provides a unifi ed point of view to describe the elastic-plastic behaviour of material since it places no requirement for a yield surface and a ’loading function’ to disting uish between loading an unloading. Endochronic theory of plasticity has been app lied to a number of metallic materials, concrete and sand, but to the knowledge of authors, no numerical scheme of the model has been applied to powder material . In the present paper, a new approach is developed based on an endochronic rate independent, density-dependent plasticity model for describing the isothermal deformation behavior of metal powder at low homologous temperature. Although the concept of yield surface has not been explicitly assumed in endochronic theory, it is shown that the cone-cap plasticity yield surface (Fig.1), which is the m ost commonly used plasticity models for describing the behavior of powder materi al can be easily derived as a special case of the proposed endochronic theory. Fig.1 Trace of cone-cap yield function on the meridian pl ane for different relative density As large deformation is observed in powder compaction process, a hypoelastic-pl astic formulation is developed in the context of finite deformation plasticity. Constitutive equations are stated in unrotated frame of reference that greatly s implifies endochronic constitutive relation in finite plasticity. Constitutive e quations of the endochronic theory and their numerical integration are establish ed and procedures for determining material parameters of the model are demonstra ted. Finally, the numerical schemes are examined for efficiency in the model ling of a tip shaped component, as shown in Fig.2. Fig.2 A shaped tip component. a) Geometry, boundary conditio n and finite element mesh; b) density distribution at final stage of
基金Projects (50934002, 51074013, 51104100) supported by the National Natural Science Foundation of China
文摘The images of granular ore media were captured by X-ray CT scanner. Combined with digital image processing and finite element techniques, the three-dimensional geometrical model, which represents the realistic pore structure of the media, was constructed. With this model, three dimensional pore scale fluid flow among particles was simulated. Then the distributions of fluid flow velocity and pressure were analyzed and the hydraulic conductivity was calculated. The simulation results indicate the fluid flow behaviors are mainly dominated by the volume and topological structure of pore space. There exist obvious preferential flow and leaching blind zones simultaneously in the medium. The highest velocities generally occur in those narrow pores with high pressure drops. The hydraulic conductivity obtained by simulation is the same order of magnitude as the laboratory test result, which denotes the validity of the model. The pore-scale and macro-scale are combined and the established geometrical model can be used for the simulations of other phenomena during heap leaching process.
基金financially supported by the National Science and Technology Major Project of China (No. 2012ZX04012011)the National Natural Science Foundation of China (No. 51375306)+1 种基金the National Basic Research Program of China (No. 2011CB012903)the Shanghai Science and Technology Innovation Action Plan (No. 14521100600)
文摘The hot deformation behavior of extruded AZ80 magnesium alloy was investigated using compression tests in the temperature range of 250–400 °C and strain rate range of 0.001–1.000 s–1. The 3D power dissipation map was developed to evaluate the hot deformation mechanisms and determine the optimal processing parameters. Two domains of dynamic recrystallization were identified from the 3D power dissipation map, with one occurring in the temperature and strain rate range of 250–320 °C and 0.001–0.010 s–1and the other one occurring in the temperature and strain rate range of 380–400 °C and 0.001–0.003 s–1. In order to delineate the regions of flow instability, Prasad’s instability criterion, Murty’s instability criterion and Gegel’s stability criteria were employed to develop the 3D instability maps. Through microstructural examination, it is found that Prasad’s and Murty’s instability criteria are more effective than Gegel’s stability criteria in predicting the flow instability regions for extruded AZ80 alloy. Further, the 3D processing maps were integrated into finite element simulation and the predictions of the simulation are in good agreement with the experimental results.
基金This work was financially supported by the National Natural Science Foundation of China(Nos.20204007 and 50390090)the Doctoral Foundation of National Education Committee of China(No.20030248008)the 863 Project of China(No.2002AA336120).
文摘A new coupled finite element formulation is proposed to calculate a conformation tensor model in two complex flows： a planar contraction flow and a planar flow around a symmetrically placed cylinder. The components of conformation tensor are first computed together with the velocity and pressure to describe the change of morphology of polymer chain coils in flow fields. Macroscopic quantities of viscoelastic flow are then calculated based on the conformation tensor. Comparisons between the numerical simulations and experiments for stress patterns and velocity profiles are carried out to prove the validity of the method.
基金funding supported by National Natural Science Foundation of China(No.52175285)Beijing Municipal Natural Science Foundation(No.3182025)+1 种基金National Defense Science and Technology Rapid support Project(No.61409230113)Scientific and Technological Innovation Foundation of Shunde Graduate School,USTB and Fundamental Research Funds for the Central Universities(No.FRFBD-20-08A,FRF-TP-20-009A2)。
文摘In this paper, a unified internal state variable(ISV) model for predicting microstructure evolution during hot working process of AZ80 magnesium alloy was developed. A novel aspect of the proposed model is that the interactive effects of material hardening, recovery and dynamic recrystallization(DRX) on the characteristic deformation behavior were considered by incorporating the evolution laws of viscoplastic flow, dislocation activities, DRX nucleation and boundary migration in a coupled manner. The model parameters were calibrated based on the experimental data analysis and genetic algorithm(GA) based objective optimization. The predicted flow stress, DRX fraction and average grain size match well with experimental results. The proposed model was embedded in the finite element(FE) software DEFORM-3 D via user defined subroutine to simulate the hot compression and equal channel angular extrusion(ECAE) processes. The heterogeneous microstructure distributions at different deformation zones and the dislocation density evolution with competitive deformation mechanisms were captured.This study can provide a theoretical solution for the hot working problems of magnesium alloy.
文摘An emerging real-time ground compaction and quality control, known as intelligent compaction(IC), has been applied for efficiently optimising the full-area compaction. Although IC technology can provide real-time assessment of uniformity of the compacted area, accurate determination of the soil stiffness required for quality control and design remains challenging. In this paper, a novel and advanced numerical model simulating the interaction of vibratory drum and soil beneath is developed. The model is capable of evaluating the nonlinear behaviour of underlying soil subjected to dynamic loading by capturing the variations of damping with the cyclic shear strains and degradation of soil modulus. The interaction of the drum and the soil is simulated via the finite element method to develop a comprehensive dataset capturing the dynamic responses of the drum and the soil. Indeed, more than a thousand three-dimensional(3D) numerical models covering various soil characteristics, roller weights, vibration amplitudes and frequencies were adopted. The developed dataset is then used to train the inverse solver using an innovative machine learning approach, i.e. the extended support vector regression, to simulate the stiffness of the compacted soil by adopting drum acceleration records. Furthermore, the impacts of the amplitude and frequency of the vibration on the level of underlying soil compaction are discussed.The proposed machine learning approach is promising for real-time extraction of actual soil stiffness during compaction. Results of the study can be employed by practising engineers to interpret roller drum acceleration data to estimate the level of compaction and ground stiffness during compaction.
