The microstructure and related property evolution induced by dynamic recrystallization(DRX)and static recrystallization(SRX)in thermo-mechanical process are two critical factors for the metal forming.The DRX and SRX a...The microstructure and related property evolution induced by dynamic recrystallization(DRX)and static recrystallization(SRX)in thermo-mechanical process are two critical factors for the metal forming.The DRX and SRX are determined by the grain level deformation and sequentially coupled.In order to fully capture the microstructure and mechanical property evolution,a crystal plasticity finite element based modelling method for DRX and SRX is proposed in the current work.The grain level deformation is calculated with crystal plasticity which is coupled with the recrystallization model straightforwardly,and both the grain deformation and microstructure evolution are updated simultaneously.The proposed method is validated with discontinuous DRX experiments and the effects of initial deformation conditions are well-captured.Two controversial mechanisms for recrystallization microstructure evolution,i.e.oriented nucleation and growth selection,are discussed in the current framework with the advantages of accurate grain level deformation and interaction predictions.Furthermore,the sequentially coupled DRX and SRX are modelled seamlessly in the current work which provides a critical method for fully integrated thermo-mechanical processes analysis.展开更多
Within the roadway advanced detection methods, DC resistivity method has an extensive application because of its simple principle and operation. Numerical simulation of the effect of focusing current on advanced detec...Within the roadway advanced detection methods, DC resistivity method has an extensive application because of its simple principle and operation. Numerical simulation of the effect of focusing current on advanced detection was carried out using a three-dimensional finite element method （FEM）, meanwhile the electric-field distribution of the point source and nine-point power source were calculated and analyzed with the same electric charges. The results show that the nine-point power source array has a very good ability to focus, and the DC focus method can be used to predict the aquifer abnormality body precisely. By comparing the FEM modelling results with physical simulation results from soil sink, it is shown that the accuracy of forward simulation meets the requirement and the artificial disturbance from roadway has no impact on the DC focus method.展开更多
The granule medium of discreteness is supposed to be continuous（Drucker-Prager model） in the existing finite element simulation analysis on the hot granule medium pressure forming（HGMF） process, so the granule med...The granule medium of discreteness is supposed to be continuous（Drucker-Prager model） in the existing finite element simulation analysis on the hot granule medium pressure forming（HGMF） process, so the granule medium may produce tensile stress in the process of pressure-transferring and flowing, which does not coincide with the reality. The analysis method, discrete element and finite element（DE-FE） coupling simulation, is proposed to solve the problem. The material parameters of simulation model are obtained by the pressure-transfer performance test of granule medium and the hot uniaxial tensile test of sheet metal. The DE-FE coupling simulation platform is established by adopting Visual Basic language. The features in the process that AA7075-T6 conical parts are formed by the HGMF process are analyzed and verified by the process test. The studies show that the results of DE-FE coupling simulation coincide well with the test results, which provides a new analysis method to solve the mechanics problem in the coupling of discrete and continuum.展开更多
A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underex...A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underexpanded axisymmetric jet. Several flow property distributions along the jet axis, including density, pres- sure and Mach number are obtained and the qualitative flowfield structures of interest are well captured using the proposed method, including shock waves, slipstreams, traveling vortex ring and multiple Mach disks. Two Mach disk locations agree well with computational and experimental measurement results. It indicates that the method is robust and efficient for solving the unsteady-state underexpanded axisymmetric jet.展开更多
As a new attempt,ball spinning was used to manufacture the nickel-titanium shape memory alloy（NiTi SMA） tube at elevated temperature.The NiTi bar with a nominal composition of Ni50.9Ti49.1（mole fraction,%） was sol...As a new attempt,ball spinning was used to manufacture the nickel-titanium shape memory alloy（NiTi SMA） tube at elevated temperature.The NiTi bar with a nominal composition of Ni50.9Ti49.1（mole fraction,%） was solution treated and was used as the original tube blank for ball spinning.Based on the variable temperature field and the constitutive equation,rigid-viscoplastic finite element method（FEM） was applied in order to simulate the ball spinning of NiTi SMA tube.The temperature field,the stress field,the strain field and the load prediction were obtained by means of FEM.FEM results reveal that there is a temperature increase of about 160 ℃ in the principal deformation zone of the spun part.It can be found from the stress fields and the strain fields that the outer wall of NiTi SMA tube is easier to meet the plastic yield criterion than the inner wall,and the plastic deformation zone is caused to be in a three-dimensional compressive stress state.The radial strain and the tangential strain are characterized by the compressive strain,while the axial strain belongs to the tensile strain.The variation of spinning loads with the progression of the ball is of great importance in predicting the stable flow of the spun part.展开更多
The micro quartz crystal tuning fork gyroscope is a new type of vibratory gyroscope. The gyroscope should be analyzed and simulated early in the design stage in order to offer reliable basis for design and to shorten ...The micro quartz crystal tuning fork gyroscope is a new type of vibratory gyroscope. The gyroscope should be analyzed and simulated early in the design stage in order to offer reliable basis for design and to shorten the period of development. Thus the vibratory characteristics of the gyroscope is simulated with the finite element method of coupled field. The optimum exciting frequency and the factors which influence the gyroscope sensitivity are determined. The method for adjusting the frequency deviation between driving and detecting modes is also proposed.展开更多
To analyze the stress state of steel orthotropic deck pavement and provide reference for the design of the overlay, the inner stress state and strain distribution of surfacing under the load of the deformation of the ...To analyze the stress state of steel orthotropic deck pavement and provide reference for the design of the overlay, the inner stress state and strain distribution of surfacing under the load of the deformation of the whole bridge structure and tyre load are analyzed by the finite element method of submodeling. Influence of surfacing modulus on the strain state of the overlay is analyzed for the purpose of the optimal design of the overlay structure. Analysis results show that the deformation of the whole bridge structure has no evident influence on the stress state of the overlay. The key factor of the overlay design is the transverse tensile strain in the overlay above the upper edge of web plate of rib. The stress state of the overlay is influenced evidently by the modulus of rigidity transform overlay. And the stress state of the overlay can be optimized and lowered by increasing the modulus and thickness of rigidity transform overlay, The fatigue test has been done to evaluate the fatigue performance and modulus of different deck pavement materials such as epoxy asphalt, SBS modified asphalt, rosphalt asphalt which can provide reference for deck pavement structure design.展开更多
Two full 3D steady mathematical models are developed by finite element method （FEM） to calcalate coupled physics fields. the electro-magnetic model is built and solved first and so is the fluid motion model with the...Two full 3D steady mathematical models are developed by finite element method （FEM） to calcalate coupled physics fields. the electro-magnetic model is built and solved first and so is the fluid motion model with the acquired electromagnetic force as source body forces in Navier-Stokes equations. Effects caused by the ferromagnetic shell, busbar system around, and open boundary problem as well as inside induced current were considered in terms of the magnetic field. Furthermore, a new modeling method is found to set up solid models and then mesh them entirely with so-called structuralized grids, namely hex-mesh. Examples of 75kA prebaked cell with two kinds of busbar arrangements are presented. Results agree with those disclosed in the literature and confirm that the coupled simulation is valid. It is also concluded that the usage of these models facilitates the consistent analysis of the electric field to magnetic field and then flow motion to the greater extent, local distributions of current density and magnetic flux density are very much dependent on the cell structure, the steel shell is a shield to reduce the magnetic field and flow pattern is two dimensional in the main body of the metal pad.展开更多
In this paper, a finite element model was developed for a turbine blade with thermal barrier coatings to investigate its failure behavior under cyclic thermal loading. Based on temperature and stress fields obtained f...In this paper, a finite element model was developed for a turbine blade with thermal barrier coatings to investigate its failure behavior under cyclic thermal loading. Based on temperature and stress fields obtained from finite element simulations, dangerous regions in ceramic coating were determined in terms of the maximum principal stress criterion. The results show that damage preferentially occurs in the chamfer and rabbet of a turbine blade with thermal barrier coatings and its thermal fatigue life decreases with the increase of thermal stress induced by high service temperature.展开更多
Recently,the finite element method(FEM) has been commonly applied in the engineering analysis of rotor dynamics.Gyroscopic moments,rotary inertia,transverse shear deformation and gravity can be included in computation...Recently,the finite element method(FEM) has been commonly applied in the engineering analysis of rotor dynamics.Gyroscopic moments,rotary inertia,transverse shear deformation and gravity can be included in computational models of rotor-bearing systems.In this paper,a finite element model and its solution method are presented for the calculation of the dynamics of dual rotor systems.A typical structure with two rotor shafts is discussed and the procedure for obtaining the coupling motion equations of the subsystems is illustrated.A computer program is developed to solve critical speeds and to simulate the transient motion.The influence of gyroscopic moments on co-rotation and counter-rotation is analyzed,and the effect of the speed ratio on critical speed is studied.The dynamic characteristics under different conditions of increasing speed during start-up are demonstrated by comparison with transient nodal displacements.The presented model provides a complete foundation for further investigation of the dynamics of dual rotor systems.展开更多
Research in different areas of orthopedic and trauma surgery requires a methodology that allows both a more economic approach and the ability to reproduce different situations in an easy way. Simulation models have be...Research in different areas of orthopedic and trauma surgery requires a methodology that allows both a more economic approach and the ability to reproduce different situations in an easy way. Simulation models have been introduced recently in bioengineering and could become an essential tool in the study of any physiological unity, regardless of its complexity. The main problem in modeling with finite elements simulation is to achieve an accurate reproduction of the anatomy and a perfect correlation of the different structures, in any region of the human body. Authors have developed a mixed technique, joining the use of a three-dimensional laser scanner Roland Picza captured together with computed tomography(CT) and 3D CT images, to achieve a perfect reproduction of the anatomy. Finite element(FE) simulation lets us know the biomechanical changes that take place after hipprostheses or osteosynthesis implantation and biological responses of bone to biomechanical changes. The simulation models are able to predict changes in bone stress distribution around the implant, so allowing preventing future pathologies. The development of a FE model of lumbar spine is another interesting application of the simulation. The model allows research on the lumbar spine, not only in physiological conditions but also simulating different load conditions, to assess the impact on biomechanics. Different degrees of disc degeneration can also be simulated to determine the impact on adjacent anatomical elements. Finally, FE models may be useful to test different fixation systems, i.e., pedicular screws, interbody devices or rigid fixations compared with the dynamic ones. We have also developed models of lumbar spine and hip joint to predict the occurrence of osteoporotic fractures, based on densitometric determinations and specific biomechanical models, including approaches from damage and fracture mechanics. FE simulations also allow us to predict the behavior of orthopedic splints applied to the correction of deformities, providing the recovering force-displacement and angle-moment curves that characterize the mechanical behavior of the splint in the overall range of movement.展开更多
Thermo-mechanical coupled finite element calculations were carried out to simulate the Gleeble compression of the samples of a titanium alloy(Ti60), and the results are analyzed and compared with the actual compressio...Thermo-mechanical coupled finite element calculations were carried out to simulate the Gleeble compression of the samples of a titanium alloy(Ti60), and the results are analyzed and compared with the actual compression tests conducted on a Gleeble 3800 thermo-mechanical simulator. The changes in temperature, stress and strain distribution in the samples and the source of error on the constitutive relations from Gleeble hot compression test were analyzed in detail. Both simulations and experiments showed that the temperature distribution in the specimen is not uniform during hot compression, resulting in significant deformation inhomogeneity and non-ignorable error in the flow stress strain relation,invalidating the uniform strain assumption commonly assumed when extracting the constitutive relation from Gleeble tests. Based on the finite element simulations with iterative corrections, we propose a scheme to refine the constitutive relations from Gleeble tests.展开更多
The evolution of microstructure in the drawing process of commercially pure aluminum wire （CPAW） does not only depend on the nature of materials, but also on the stress profile. In this study, the effect of stress p...The evolution of microstructure in the drawing process of commercially pure aluminum wire （CPAW） does not only depend on the nature of materials, but also on the stress profile. In this study, the effect of stress profile on the texture evolution of the CPAW was systematically investigated by combining the numerical simulation and the microstructure observation. The results show that the tensile stress at the wire center promotes the formation of 〈111〉 texture, whereas the shear stress nearby the rim makes little contribution to the texture formation. Therefore, the 〈111 〉 texture at the wire center is stronger than that in the surface layer, which also results in a higher microhardness at the center of the CPAW under axial loading.2017 Published by Elsevier Ltd on behalf of The editorial office of Journal of Materials Science ＆ Technology.展开更多
In counter-rotating electrochemical machining （CRECM）, a revolving cathode tool with hollow windows of various shapes is used to fabricate convex structures on a revolving part. During this process, the anode workpi...In counter-rotating electrochemical machining （CRECM）, a revolving cathode tool with hollow windows of various shapes is used to fabricate convex structures on a revolving part. During this process, the anode workpiece and the cathode tool rotate relative to each other at the same rotation speed. In contrast to the conventional schemes of ECM machining with linear motion of a block tool electrode, this scheme of ECM is unique, and has not been adequately studied yet. In this paper, the finite element method （FEM） is used to simulate the anode shaping process during CRECM, and the simulation process which involves a meshing model, a moving boundary, and a simulation algorithm is described. The simulated anode profiles of the convex structure at different processing times show that the CRECM process can be used to fabricate convex structures of various shapes with different heights. Besides, the variation of the inter-electrode gap indicates that this process can also reach a relative equilibrium state like that in conventional ECM. A rectangular convex and a circular convex are successfully fabricated on revolving parts. The experimental results indicate relatively good agreement with the simulation results. The proposed simulation process is valid for convex shaping prediction and feasibility studies as well.展开更多
To investigate the effects of billet geometry on the cold precision forging process of a helical gear, six different billet geometries were designed utilizing the relief-hole principle. And the influences of the bille...To investigate the effects of billet geometry on the cold precision forging process of a helical gear, six different billet geometries were designed utilizing the relief-hole principle. And the influences of the billet geometry on the forming load and the deformation uniformity were analyzed by three-dimensional （3D） finite element method （FEM） under the commercial software DEFORM 3D. The billet geometry was optimized to meet lower forming load and better deformation uniformity requirement. Deformation mechanism was studied through the distribution of flow velocity field and effective strain field. The forging experiments of the helical gear were successfully performed using lead material as a model material under the same process conditions used in the FE simulations. The results show that the forming load decreases as the diameter of relief-hole do increases, but the effect of do on the deformation uniformity is very complicated. The forming load is lower and the deformation is more uniform when do is 10 mm.展开更多
The finite element method （FEM） and particle image velocimetry （PIV） technique are utilized to get the flow field along the inlet passage, the chamber, the metering port and the outlet passage of spool valve at th...The finite element method （FEM） and particle image velocimetry （PIV） technique are utilized to get the flow field along the inlet passage, the chamber, the metering port and the outlet passage of spool valve at three different valve openings. For FEM numerical simulation, the stream function ψ-vorticity ω forms of continuity and Navier-Stokes equations are employed and FEM is applied to discrete the equations. Homemade simulation codes are executed to compute the values of stream function and vorticity at each node in the flow domain, then according to the correlation between stream function and velocity components, the velocity vectors of the whole field are calculated. For PIV experiment, pulse Nd： YAG laser is exploited to generate laser beam, cylindrical and spherical lenses are combined each other to produce 1.0 mm thickness laser sheet to illuminate the object plane, Polystyrene spherical particle with diameter of 30-50 μm is seeded in the fluid as a tracing particles, Kodak ES 1.0 CCD camera is employed to capture the images of interested, the images are processed with fast Fourier transform （FFT） cross-correlation algorithm and the processing results is displayed. Both results of numerical simulation and PIV experimental show that there are three main areas in the spool valve where vortex is formed. Numerical results also indicate that the valve opening have some effects on the flow structure of the valve. The investigation is helpful for qualitatively analyzing the energy loss, noise generating, steady state flow forces and even designing the geometry structure and flow passage.展开更多
Based on the theory of porous media, a general Gurtin variational principle for the initial boundary value problem of dynamical response of fluid-saturated elastic porous media is developed by assuming infinitesimal d...Based on the theory of porous media, a general Gurtin variational principle for the initial boundary value problem of dynamical response of fluid-saturated elastic porous media is developed by assuming infinitesimal deformation and incompressible constituents of the solid and fluid phase. The finite element formulation based on this variational principle is also derived. As the functional of the variational principle is a spatial integral of the convolution formulation, the general finite element discretization in space results in symmetrical differential-integral equations in the time domain. In some situations, the differential-integral equations can be reduced to symmetrical differential equations and, as a numerical example, it is employed to analyze the reflection of one-dimensional longitudinal wave in a fluid-saturated porous solid. The numerical results can provide further understanding of the wave propagation in porous media.展开更多
To study the hot deformation behavior of Mg-8.3 Gd-4.4 Y-1.5 Zn-0.8 Mn(wt%) alloy,hot compression tests were conducted using a Gleeble-3500 thermal simulator at temperatures ranging from 653 to773 K,true strain rates ...To study the hot deformation behavior of Mg-8.3 Gd-4.4 Y-1.5 Zn-0.8 Mn(wt%) alloy,hot compression tests were conducted using a Gleeble-3500 thermal simulator at temperatures ranging from 653 to773 K,true strain rates of 0.001-1 s^(-1),and a deformation degree of 60%.Results of hot compression experiments show that the flow stress of the alloy increases with the strain rate.The true stress-true strain curves are corrected by correcting the effect of temperature rise in the deformation process.Activation energy,Q,equal to 287380 J/mol and material constant,n,equal to 4.59 were calculated by fitting the true stress-true strain curves.Then,the constitutive equation was established and verified via finite element simulation.Results of the hot processing map show that the probability of material instability increases with the degree of deformation,which indicates that the material is not suitable for large deformation in a single pass.On the whole,the alloy is appropriate for multipass processing with small deformation and a suitable processing temperature and strain rate are 733 K and 0.01 s-1,respectively.展开更多
Materials with the same elastic modulus E and representative stress and strain （σr,εr） present similar indentation-loading curves, whatever the value of strain hardening exponent n. Based on this definition, a goo...Materials with the same elastic modulus E and representative stress and strain （σr,εr） present similar indentation-loading curves, whatever the value of strain hardening exponent n. Based on this definition, a good approach was proposed to extract the plastic properties or constitutive equations of metals from nanoindentation test combining finite element simulation. Firstly, without consideration of strain hardening, the representative stress was determined by varying assumed representative stress over a wide range until a good agreement was reached between the computed and experimental loading curves. Similarly, the corresponding representative strain was determined with different hypothetical values of strain hardening exponent in the range of 0-0.6. Through modulating assumed strain hardening exponent values to make the computed unloading curve coincide with that of the experiment, the real strain hardening exponent was acquired. Once the strain hardening exponent was determined, the initial yield stress ay of metals could be obtained by the power law constitution. The validity of the proposed methodology was verified by three real metals： AISI 304 steel, Fe andA1 alloy.展开更多
A thermo-mechanical coupling.finite element model was built to investigate the inertia friction welding of GH4169 bars. The remeshing and map solution techniques were adopted. Ttle whole welding process was investigat...A thermo-mechanical coupling.finite element model was built to investigate the inertia friction welding of GH4169 bars. The remeshing and map solution techniques were adopted. Ttle whole welding process was investigated by adopting an innovative heat generation model and the flywheel rotational speed measured via the experiment. The simulated evolution of axial shortening shows a good agreement with the experiment. In addition, extensive .strain concentration presents in the interface and flash, and the largest ,strain exists near the flash root. Moreover, an intere.sting thermal reflux phenomenon during the cooling stage was found.展开更多
基金supported by the National Natural Science Foundation of China(Nos.52105384 and U2141215).
文摘The microstructure and related property evolution induced by dynamic recrystallization(DRX)and static recrystallization(SRX)in thermo-mechanical process are two critical factors for the metal forming.The DRX and SRX are determined by the grain level deformation and sequentially coupled.In order to fully capture the microstructure and mechanical property evolution,a crystal plasticity finite element based modelling method for DRX and SRX is proposed in the current work.The grain level deformation is calculated with crystal plasticity which is coupled with the recrystallization model straightforwardly,and both the grain deformation and microstructure evolution are updated simultaneously.The proposed method is validated with discontinuous DRX experiments and the effects of initial deformation conditions are well-captured.Two controversial mechanisms for recrystallization microstructure evolution,i.e.oriented nucleation and growth selection,are discussed in the current framework with the advantages of accurate grain level deformation and interaction predictions.Furthermore,the sequentially coupled DRX and SRX are modelled seamlessly in the current work which provides a critical method for fully integrated thermo-mechanical processes analysis.
基金Project(41174103)supported by the National Natural Science Foundation of ChinaProject(20110162130008)supported by the PhD Program Foundation of Ministry of Education of ChinaProject(2011BAB04B08)supported by the National Key Technology R&D Program during the 12th Five-Year Plan of China
文摘Within the roadway advanced detection methods, DC resistivity method has an extensive application because of its simple principle and operation. Numerical simulation of the effect of focusing current on advanced detection was carried out using a three-dimensional finite element method （FEM）, meanwhile the electric-field distribution of the point source and nine-point power source were calculated and analyzed with the same electric charges. The results show that the nine-point power source array has a very good ability to focus, and the DC focus method can be used to predict the aquifer abnormality body precisely. By comparing the FEM modelling results with physical simulation results from soil sink, it is shown that the accuracy of forward simulation meets the requirement and the artificial disturbance from roadway has no impact on the DC focus method.
基金Projects(5130538651305385)supported by the National Natural Science Foundation of China+1 种基金Project(E2013203093)supported by the Natural Science Foundation of Hebei ProvinceChina
文摘The granule medium of discreteness is supposed to be continuous（Drucker-Prager model） in the existing finite element simulation analysis on the hot granule medium pressure forming（HGMF） process, so the granule medium may produce tensile stress in the process of pressure-transferring and flowing, which does not coincide with the reality. The analysis method, discrete element and finite element（DE-FE） coupling simulation, is proposed to solve the problem. The material parameters of simulation model are obtained by the pressure-transfer performance test of granule medium and the hot uniaxial tensile test of sheet metal. The DE-FE coupling simulation platform is established by adopting Visual Basic language. The features in the process that AA7075-T6 conical parts are formed by the HGMF process are analyzed and verified by the process test. The studies show that the results of DE-FE coupling simulation coincide well with the test results, which provides a new analysis method to solve the mechanics problem in the coupling of discrete and continuum.
文摘A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underexpanded axisymmetric jet. Several flow property distributions along the jet axis, including density, pres- sure and Mach number are obtained and the qualitative flowfield structures of interest are well captured using the proposed method, including shock waves, slipstreams, traveling vortex ring and multiple Mach disks. Two Mach disk locations agree well with computational and experimental measurement results. It indicates that the method is robust and efficient for solving the unsteady-state underexpanded axisymmetric jet.
基金Project(51071056) supported by the National Natural Science Foundation of ChinaProject(HEUCF121712) supported by the Fundamental Research Funds for the Central Universities of China
文摘As a new attempt,ball spinning was used to manufacture the nickel-titanium shape memory alloy（NiTi SMA） tube at elevated temperature.The NiTi bar with a nominal composition of Ni50.9Ti49.1（mole fraction,%） was solution treated and was used as the original tube blank for ball spinning.Based on the variable temperature field and the constitutive equation,rigid-viscoplastic finite element method（FEM） was applied in order to simulate the ball spinning of NiTi SMA tube.The temperature field,the stress field,the strain field and the load prediction were obtained by means of FEM.FEM results reveal that there is a temperature increase of about 160 ℃ in the principal deformation zone of the spun part.It can be found from the stress fields and the strain fields that the outer wall of NiTi SMA tube is easier to meet the plastic yield criterion than the inner wall,and the plastic deformation zone is caused to be in a three-dimensional compressive stress state.The radial strain and the tangential strain are characterized by the compressive strain,while the axial strain belongs to the tensile strain.The variation of spinning loads with the progression of the ball is of great importance in predicting the stable flow of the spun part.
文摘The micro quartz crystal tuning fork gyroscope is a new type of vibratory gyroscope. The gyroscope should be analyzed and simulated early in the design stage in order to offer reliable basis for design and to shorten the period of development. Thus the vibratory characteristics of the gyroscope is simulated with the finite element method of coupled field. The optimum exciting frequency and the factors which influence the gyroscope sensitivity are determined. The method for adjusting the frequency deviation between driving and detecting modes is also proposed.
文摘To analyze the stress state of steel orthotropic deck pavement and provide reference for the design of the overlay, the inner stress state and strain distribution of surfacing under the load of the deformation of the whole bridge structure and tyre load are analyzed by the finite element method of submodeling. Influence of surfacing modulus on the strain state of the overlay is analyzed for the purpose of the optimal design of the overlay structure. Analysis results show that the deformation of the whole bridge structure has no evident influence on the stress state of the overlay. The key factor of the overlay design is the transverse tensile strain in the overlay above the upper edge of web plate of rib. The stress state of the overlay is influenced evidently by the modulus of rigidity transform overlay. And the stress state of the overlay can be optimized and lowered by increasing the modulus and thickness of rigidity transform overlay, The fatigue test has been done to evaluate the fatigue performance and modulus of different deck pavement materials such as epoxy asphalt, SBS modified asphalt, rosphalt asphalt which can provide reference for deck pavement structure design.
基金the National High Technical Reasearch and Development Programme of China （No. 2003AA327140） the National Natural Science Foundation of China （No. 50374081）.
文摘Two full 3D steady mathematical models are developed by finite element method （FEM） to calcalate coupled physics fields. the electro-magnetic model is built and solved first and so is the fluid motion model with the acquired electromagnetic force as source body forces in Navier-Stokes equations. Effects caused by the ferromagnetic shell, busbar system around, and open boundary problem as well as inside induced current were considered in terms of the magnetic field. Furthermore, a new modeling method is found to set up solid models and then mesh them entirely with so-called structuralized grids, namely hex-mesh. Examples of 75kA prebaked cell with two kinds of busbar arrangements are presented. Results agree with those disclosed in the literature and confirm that the coupled simulation is valid. It is also concluded that the usage of these models facilitates the consistent analysis of the electric field to magnetic field and then flow motion to the greater extent, local distributions of current density and magnetic flux density are very much dependent on the cell structure, the steel shell is a shield to reduce the magnetic field and flow pattern is two dimensional in the main body of the metal pad.
基金supported by the National Natural Science Foundation of China(Nos.11002122,51172192,11272275 and 11002121)the Natural Science Foundation of Hunan Province, China(No.11JJ4003)the Key Project of Scientific Research Conditions in Hunan Province,China(No.2012TT2040)
文摘In this paper, a finite element model was developed for a turbine blade with thermal barrier coatings to investigate its failure behavior under cyclic thermal loading. Based on temperature and stress fields obtained from finite element simulations, dangerous regions in ceramic coating were determined in terms of the maximum principal stress criterion. The results show that damage preferentially occurs in the chamfer and rabbet of a turbine blade with thermal barrier coatings and its thermal fatigue life decreases with the increase of thermal stress induced by high service temperature.
文摘Recently,the finite element method(FEM) has been commonly applied in the engineering analysis of rotor dynamics.Gyroscopic moments,rotary inertia,transverse shear deformation and gravity can be included in computational models of rotor-bearing systems.In this paper,a finite element model and its solution method are presented for the calculation of the dynamics of dual rotor systems.A typical structure with two rotor shafts is discussed and the procedure for obtaining the coupling motion equations of the subsystems is illustrated.A computer program is developed to solve critical speeds and to simulate the transient motion.The influence of gyroscopic moments on co-rotation and counter-rotation is analyzed,and the effect of the speed ratio on critical speed is studied.The dynamic characteristics under different conditions of increasing speed during start-up are demonstrated by comparison with transient nodal displacements.The presented model provides a complete foundation for further investigation of the dynamics of dual rotor systems.
文摘Research in different areas of orthopedic and trauma surgery requires a methodology that allows both a more economic approach and the ability to reproduce different situations in an easy way. Simulation models have been introduced recently in bioengineering and could become an essential tool in the study of any physiological unity, regardless of its complexity. The main problem in modeling with finite elements simulation is to achieve an accurate reproduction of the anatomy and a perfect correlation of the different structures, in any region of the human body. Authors have developed a mixed technique, joining the use of a three-dimensional laser scanner Roland Picza captured together with computed tomography(CT) and 3D CT images, to achieve a perfect reproduction of the anatomy. Finite element(FE) simulation lets us know the biomechanical changes that take place after hipprostheses or osteosynthesis implantation and biological responses of bone to biomechanical changes. The simulation models are able to predict changes in bone stress distribution around the implant, so allowing preventing future pathologies. The development of a FE model of lumbar spine is another interesting application of the simulation. The model allows research on the lumbar spine, not only in physiological conditions but also simulating different load conditions, to assess the impact on biomechanics. Different degrees of disc degeneration can also be simulated to determine the impact on adjacent anatomical elements. Finally, FE models may be useful to test different fixation systems, i.e., pedicular screws, interbody devices or rigid fixations compared with the dynamic ones. We have also developed models of lumbar spine and hip joint to predict the occurrence of osteoporotic fractures, based on densitometric determinations and specific biomechanical models, including approaches from damage and fracture mechanics. FE simulations also allow us to predict the behavior of orthopedic splints applied to the correction of deformities, providing the recovering force-displacement and angle-moment curves that characterize the mechanical behavior of the splint in the overall range of movement.
基金the National Key Research and Development Program of China(2016YFB0701300,2011CB606404)the CAS informatization Program(XXH13506-304)+1 种基金the National Natural Science Foundation of China(51171195,50911130367)the CAS Shenyang Supercomputing Center
文摘Thermo-mechanical coupled finite element calculations were carried out to simulate the Gleeble compression of the samples of a titanium alloy(Ti60), and the results are analyzed and compared with the actual compression tests conducted on a Gleeble 3800 thermo-mechanical simulator. The changes in temperature, stress and strain distribution in the samples and the source of error on the constitutive relations from Gleeble hot compression test were analyzed in detail. Both simulations and experiments showed that the temperature distribution in the specimen is not uniform during hot compression, resulting in significant deformation inhomogeneity and non-ignorable error in the flow stress strain relation,invalidating the uniform strain assumption commonly assumed when extracting the constitutive relation from Gleeble tests. Based on the finite element simulations with iterative corrections, we propose a scheme to refine the constitutive relations from Gleeble tests.
基金financially supported by the State Grid Corporation of China (No. 52110416001z)the National Natural Science Foundation of China (No. 51331007)
文摘The evolution of microstructure in the drawing process of commercially pure aluminum wire （CPAW） does not only depend on the nature of materials, but also on the stress profile. In this study, the effect of stress profile on the texture evolution of the CPAW was systematically investigated by combining the numerical simulation and the microstructure observation. The results show that the tensile stress at the wire center promotes the formation of 〈111〉 texture, whereas the shear stress nearby the rim makes little contribution to the texture formation. Therefore, the 〈111 〉 texture at the wire center is stronger than that in the surface layer, which also results in a higher microhardness at the center of the CPAW under axial loading.2017 Published by Elsevier Ltd on behalf of The editorial office of Journal of Materials Science ＆ Technology.
基金supported by the Program for New Century Excellent Talents in University of China(NCET-10-0074)
文摘In counter-rotating electrochemical machining （CRECM）, a revolving cathode tool with hollow windows of various shapes is used to fabricate convex structures on a revolving part. During this process, the anode workpiece and the cathode tool rotate relative to each other at the same rotation speed. In contrast to the conventional schemes of ECM machining with linear motion of a block tool electrode, this scheme of ECM is unique, and has not been adequately studied yet. In this paper, the finite element method （FEM） is used to simulate the anode shaping process during CRECM, and the simulation process which involves a meshing model, a moving boundary, and a simulation algorithm is described. The simulated anode profiles of the convex structure at different processing times show that the CRECM process can be used to fabricate convex structures of various shapes with different heights. Besides, the variation of the inter-electrode gap indicates that this process can also reach a relative equilibrium state like that in conventional ECM. A rectangular convex and a circular convex are successfully fabricated on revolving parts. The experimental results indicate relatively good agreement with the simulation results. The proposed simulation process is valid for convex shaping prediction and feasibility studies as well.
基金Project(51105287)supported by the National Natural Science Foundation of China
文摘To investigate the effects of billet geometry on the cold precision forging process of a helical gear, six different billet geometries were designed utilizing the relief-hole principle. And the influences of the billet geometry on the forming load and the deformation uniformity were analyzed by three-dimensional （3D） finite element method （FEM） under the commercial software DEFORM 3D. The billet geometry was optimized to meet lower forming load and better deformation uniformity requirement. Deformation mechanism was studied through the distribution of flow velocity field and effective strain field. The forging experiments of the helical gear were successfully performed using lead material as a model material under the same process conditions used in the FE simulations. The results show that the forming load decreases as the diameter of relief-hole do increases, but the effect of do on the deformation uniformity is very complicated. The forming load is lower and the deformation is more uniform when do is 10 mm.
文摘The finite element method （FEM） and particle image velocimetry （PIV） technique are utilized to get the flow field along the inlet passage, the chamber, the metering port and the outlet passage of spool valve at three different valve openings. For FEM numerical simulation, the stream function ψ-vorticity ω forms of continuity and Navier-Stokes equations are employed and FEM is applied to discrete the equations. Homemade simulation codes are executed to compute the values of stream function and vorticity at each node in the flow domain, then according to the correlation between stream function and velocity components, the velocity vectors of the whole field are calculated. For PIV experiment, pulse Nd： YAG laser is exploited to generate laser beam, cylindrical and spherical lenses are combined each other to produce 1.0 mm thickness laser sheet to illuminate the object plane, Polystyrene spherical particle with diameter of 30-50 μm is seeded in the fluid as a tracing particles, Kodak ES 1.0 CCD camera is employed to capture the images of interested, the images are processed with fast Fourier transform （FFT） cross-correlation algorithm and the processing results is displayed. Both results of numerical simulation and PIV experimental show that there are three main areas in the spool valve where vortex is formed. Numerical results also indicate that the valve opening have some effects on the flow structure of the valve. The investigation is helpful for qualitatively analyzing the energy loss, noise generating, steady state flow forces and even designing the geometry structure and flow passage.
基金Project supported by the National Nattural Science Foundation of China(No.10272070)
文摘Based on the theory of porous media, a general Gurtin variational principle for the initial boundary value problem of dynamical response of fluid-saturated elastic porous media is developed by assuming infinitesimal deformation and incompressible constituents of the solid and fluid phase. The finite element formulation based on this variational principle is also derived. As the functional of the variational principle is a spatial integral of the convolution formulation, the general finite element discretization in space results in symmetrical differential-integral equations in the time domain. In some situations, the differential-integral equations can be reduced to symmetrical differential equations and, as a numerical example, it is employed to analyze the reflection of one-dimensional longitudinal wave in a fluid-saturated porous solid. The numerical results can provide further understanding of the wave propagation in porous media.
基金Project supported by the General Program of National Natural Science Foundation of China (51874062)。
文摘To study the hot deformation behavior of Mg-8.3 Gd-4.4 Y-1.5 Zn-0.8 Mn(wt%) alloy,hot compression tests were conducted using a Gleeble-3500 thermal simulator at temperatures ranging from 653 to773 K,true strain rates of 0.001-1 s^(-1),and a deformation degree of 60%.Results of hot compression experiments show that the flow stress of the alloy increases with the strain rate.The true stress-true strain curves are corrected by correcting the effect of temperature rise in the deformation process.Activation energy,Q,equal to 287380 J/mol and material constant,n,equal to 4.59 were calculated by fitting the true stress-true strain curves.Then,the constitutive equation was established and verified via finite element simulation.Results of the hot processing map show that the probability of material instability increases with the degree of deformation,which indicates that the material is not suitable for large deformation in a single pass.On the whole,the alloy is appropriate for multipass processing with small deformation and a suitable processing temperature and strain rate are 733 K and 0.01 s-1,respectively.
基金Project （51171125） supported by the National Natural Science Foundation of China Project （20110321051 ） supported by the Science and Technology Key Project of Shanxi Province, China
文摘Materials with the same elastic modulus E and representative stress and strain （σr,εr） present similar indentation-loading curves, whatever the value of strain hardening exponent n. Based on this definition, a good approach was proposed to extract the plastic properties or constitutive equations of metals from nanoindentation test combining finite element simulation. Firstly, without consideration of strain hardening, the representative stress was determined by varying assumed representative stress over a wide range until a good agreement was reached between the computed and experimental loading curves. Similarly, the corresponding representative strain was determined with different hypothetical values of strain hardening exponent in the range of 0-0.6. Through modulating assumed strain hardening exponent values to make the computed unloading curve coincide with that of the experiment, the real strain hardening exponent was acquired. Once the strain hardening exponent was determined, the initial yield stress ay of metals could be obtained by the power law constitution. The validity of the proposed methodology was verified by three real metals： AISI 304 steel, Fe andA1 alloy.
基金The work is supported by the National Natural Science Foundation of＂ China （51005180）.
文摘A thermo-mechanical coupling.finite element model was built to investigate the inertia friction welding of GH4169 bars. The remeshing and map solution techniques were adopted. Ttle whole welding process was investigated by adopting an innovative heat generation model and the flywheel rotational speed measured via the experiment. The simulated evolution of axial shortening shows a good agreement with the experiment. In addition, extensive .strain concentration presents in the interface and flash, and the largest ,strain exists near the flash root. Moreover, an intere.sting thermal reflux phenomenon during the cooling stage was found.
