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.展开更多
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>展开更多
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.展开更多
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 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.展开更多
Coriolis effects,encompassing the dilative,compressive,and deflective manifestations,constitute pivotal considerations in the centrifugal modelling of high-speed granular run-out processes.Notably,under the deflective...Coriolis effects,encompassing the dilative,compressive,and deflective manifestations,constitute pivotal considerations in the centrifugal modelling of high-speed granular run-out processes.Notably,under the deflective Coriolis condition,the velocity component parallel to the rotational axis exerts no influence on the magnitude of Coriolis acceleration.This circumstance implies a potential mitigation of the Coriolis force's deflective impact.Regrettably,extant investigations predominantly emphasize the dilative and compressive Coriolis effects,largely neglecting the pragmatic import of the deflective Coriolis condition.In pursuit of this gap,a series of discrete element method(DEM)simulations have been conducted to scrutinize the feasibility of centrifugal modelling for dry granular run-out processes under deflective Coriolis conditions.The findings concerning the deflective Coriolis effect reveal a consistent rise in the run-out distance by 2%–16%,a modest increase in bulk flow velocity of under 4%,and a slight elevation in average flow depth by no more than 25%.These alterations display smaller dependence on the specific testing conditions due to the granular flow undergoing dual deflections in opposing directions.This underscores the significance and utility of the deflective Coriolis condition.Notably,the anticipated reduction in error in predicting the final run-out distance is substantial,potentially reaching a 150%improvement compared to predictions made under the dilative and compressive Coriolis conditions.Therefore,the deflective Coriolis condition is advised when the final run-out distance of the granular flow is the main concern.To mitigate the impact of Coriolis acceleration,a greater initial height of the granular column is recommended,with a height/width ratio exceeding 1,as the basal friction of the granular material plays a crucial role in mitigating the deflective Coriolis effect.For more transverse-uniform flow properties,the width of the granular column should be as large as possible.展开更多
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.展开更多
The flexible structure of photonic crystal fibre not only offers novel optical properties but also brings some difficulties in keeping the fibre structure in the fabrication process which inevitably cause the optical ...The flexible structure of photonic crystal fibre not only offers novel optical properties but also brings some difficulties in keeping the fibre structure in the fabrication process which inevitably cause the optical properties of the resulting fibre to deviate from the designed properties. Therefore, a method of evaluating the optical properties of the actual fibre is necessary for the purpose of application. Up to now, the methods employed to measure the properties of the actual photonic crystal fibre often require long fibre samples or complex expensive equipments. To our knowledge, there are few studies of modeling an actual photonic crystal fibre and evaluating its properties rapidly. In this paper, a novel method, based on the combination model of digital image processing and the finite element method, is proposed to rapidly model the optical properties of the actual photonic crystal fibre. Two kinds of photonic crystal fibres made by Crystal Fiber A/S are modeled. It is confirmed from numerical results that the proposed method is simple, rapid and accurate for evaluating the optical properties of the actual photonic crystal fibre without requiring complex equipment.展开更多
An approach was presented to characterize the stress response of workpiece in hard machining, accounting for the effect of the initial workpiece hardness in addition to temperature, strain and strain rate on flow stre...An approach was presented to characterize the stress response of workpiece in hard machining, accounting for the effect of the initial workpiece hardness in addition to temperature, strain and strain rate on flow stress in this paper. AISI H13 die steel was chosen to verify this methodology. The proposed flow stress model demonstrates a good agreement with experimental data. Therefore, the proposed model can be used to predict the corresponding flow stress-strain response ofAISl H13 die steel with variation of the initial workpiece hardness in hard machining.展开更多
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.
基金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>
基金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.
基金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.
基金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.
基金supported by the National Natural Science Foundation of China(Grant Nos.42120104008 and 42307214)the Postdoctoral Fellowship Program of CPSF(Grant No.GZB20230620).
文摘Coriolis effects,encompassing the dilative,compressive,and deflective manifestations,constitute pivotal considerations in the centrifugal modelling of high-speed granular run-out processes.Notably,under the deflective Coriolis condition,the velocity component parallel to the rotational axis exerts no influence on the magnitude of Coriolis acceleration.This circumstance implies a potential mitigation of the Coriolis force's deflective impact.Regrettably,extant investigations predominantly emphasize the dilative and compressive Coriolis effects,largely neglecting the pragmatic import of the deflective Coriolis condition.In pursuit of this gap,a series of discrete element method(DEM)simulations have been conducted to scrutinize the feasibility of centrifugal modelling for dry granular run-out processes under deflective Coriolis conditions.The findings concerning the deflective Coriolis effect reveal a consistent rise in the run-out distance by 2%–16%,a modest increase in bulk flow velocity of under 4%,and a slight elevation in average flow depth by no more than 25%.These alterations display smaller dependence on the specific testing conditions due to the granular flow undergoing dual deflections in opposing directions.This underscores the significance and utility of the deflective Coriolis condition.Notably,the anticipated reduction in error in predicting the final run-out distance is substantial,potentially reaching a 150%improvement compared to predictions made under the dilative and compressive Coriolis conditions.Therefore,the deflective Coriolis condition is advised when the final run-out distance of the granular flow is the main concern.To mitigate the impact of Coriolis acceleration,a greater initial height of the granular column is recommended,with a height/width ratio exceeding 1,as the basal friction of the granular material plays a crucial role in mitigating the deflective Coriolis effect.For more transverse-uniform flow properties,the width of the granular column should be as large as possible.
基金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.
基金Project supported by the National Basic Research Program of China(Grant No.2010CB328206)the National Natural Science Foundation of China(Grant No.60977033)the Science and Technology Innovation Foundation for Excellent Doctors of Beijing Jiaotong University,China(Grant Nos.141055522 and 141060522)
文摘The flexible structure of photonic crystal fibre not only offers novel optical properties but also brings some difficulties in keeping the fibre structure in the fabrication process which inevitably cause the optical properties of the resulting fibre to deviate from the designed properties. Therefore, a method of evaluating the optical properties of the actual fibre is necessary for the purpose of application. Up to now, the methods employed to measure the properties of the actual photonic crystal fibre often require long fibre samples or complex expensive equipments. To our knowledge, there are few studies of modeling an actual photonic crystal fibre and evaluating its properties rapidly. In this paper, a novel method, based on the combination model of digital image processing and the finite element method, is proposed to rapidly model the optical properties of the actual photonic crystal fibre. Two kinds of photonic crystal fibres made by Crystal Fiber A/S are modeled. It is confirmed from numerical results that the proposed method is simple, rapid and accurate for evaluating the optical properties of the actual photonic crystal fibre without requiring complex equipment.
基金the Natural Science Foundation of Jiangxi Province(No.550067)the National Natural Science Foundation of China(No.50465003)Foundation of the State Key Laboratory of Plastic Forming Simulation and Die & Mould Technology(No.06-3)
文摘An approach was presented to characterize the stress response of workpiece in hard machining, accounting for the effect of the initial workpiece hardness in addition to temperature, strain and strain rate on flow stress in this paper. AISI H13 die steel was chosen to verify this methodology. The proposed flow stress model demonstrates a good agreement with experimental data. Therefore, the proposed model can be used to predict the corresponding flow stress-strain response ofAISl H13 die steel with variation of the initial workpiece hardness in hard machining.
文摘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.
