Recent advancements in additive manufacturing(AM)have revolutionized the design and production of complex engineering microstructures.Despite these advancements,their mathematical modeling and computational analysis r...Recent advancements in additive manufacturing(AM)have revolutionized the design and production of complex engineering microstructures.Despite these advancements,their mathematical modeling and computational analysis remain significant challenges.This research aims to develop an effective computational method for analyzing the free vibration of functionally graded(FG)microplates under high temperatures while resting on a Pasternak foundation(PF).This formulation leverages a new thirdorder shear deformation theory(new TSDT)for improved accuracy without requiring shear correction factors.Additionally,the modified couple stress theory(MCST)is incorporated to account for sizedependent effects in microplates.The PF is characterized by two parameters including spring stiffness(k_(w))and shear layer stiffness(k_(s)).To validate the proposed method,the results obtained are compared with those of the existing literature.Furthermore,numerical examples explore the influence of various factors on the high-temperature free vibration of FG microplates.These factors include the length scale parameter(l),geometric dimensions,material properties,and the presence of the elastic foundation.The findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the results of this research will have great potential in military and defense applications such as components of submarines,fighter aircraft,and missiles.展开更多
At the first time,the finite element method was used to model and analyze the free vibration and transient response of non-uniform thickness bi-directional functionally graded sandwich porous(BFGSP)skew plates.The who...At the first time,the finite element method was used to model and analyze the free vibration and transient response of non-uniform thickness bi-directional functionally graded sandwich porous(BFGSP)skew plates.The whole BFGSP skew-plates is placed on a variable visco-elastic foundation(VEF)in the hygro-thermal environment and subjected to the blast load.The BFGSP skew-plate thickness is permitted to vary non-linearly over both the length and width of the skew-plate,thereby faithfully representing the real behavior of the structure itself.The analysis is based on a four-node planar quadrilateral element with eight degrees of freedom per node,which is approximated using Lagrange Q_(4)shape function and C^(1)level non-conforming Hermite shape function based on refined higher-order shear deformation plate theory.The forced vibration parameters of the non-uniform thickness BFGSP skew-plate are fully determined using Hamilton's principle and the Newmark-βdirect integration technique.Accuracy of the calculation program is validated by comparing its numerical results with those from reputable sources.Furthermore,a thorough assessment is conducted to determine the impact of various parameters on the free and forced vibration responses of the non-uniform thickness BFGSP skew-plate.The findings of the paper may be used in the development of civil and military structures in situations that are prone to exceptional forces,such as explosions and impacts load.展开更多
The strategies that minimize the overall solution time of multiple linear systems in 3D finite element method (FEM) modeling of direct current (DC) resistivity were discussed. A global stiff matrix is assembled and st...The strategies that minimize the overall solution time of multiple linear systems in 3D finite element method (FEM) modeling of direct current (DC) resistivity were discussed. A global stiff matrix is assembled and stored in two parts separately. One part is associated with the volume integral and the other is associated with the subsurface boundary integral. The equivalent multiple linear systems with closer right-hand sides than the original systems were constructed. A recycling Krylov subspace technique was employed to solve the multiple linear systems. The solution of the seed system was used as an initial guess for the subsequent systems. The results of two numerical experiments show that the improved algorithm reduces the iterations and CPU time by almost 50%, compared with the classical preconditioned conjugate gradient method.展开更多
In an injection moulding process, the mould cooling s ystem is very important as an efficient and uniform cooling effect can improve b oth the productivity and part quality. Due to the complexity of the process, muc h...In an injection moulding process, the mould cooling s ystem is very important as an efficient and uniform cooling effect can improve b oth the productivity and part quality. Due to the complexity of the process, muc h research on the mould cooling analysis and cooling design optimization has bee n focused on the core and the cavity, excluding other systems of injection mould s. However, the runner system introduces a considerable amount of heat into the mould. In recent years, more and more hot runner systems are being applied in th e moulding industry to save material and decrease losses of injection pressure. This raises the need to include the hot runner system in the cooling analysis. I n this paper, a photo frame part was studied. The mould was built with a hot run ner system. Two thermal sensors were installed: one measures the temperature of lateral surface of hot runner nozzle; the other measures the plastic temperature from the core side. A pressure sensor was also installed to measure the pressur e of the core impression. Cooling analysis was performed using ABAQUS, ananalysi s software based on the Finite Element Method (FEM). The assembly including core , cavity and plastic part was modeled. Heat conduction from hot runner to cavity and from polymer melt to the mould and force convection on the cooling channel surfaces were studied. The natural convection between the ambient air and the ex terior mould surface was ignored. The simulations were adjusted with the experim ental results to find out the heat input from hot runner and its influence on mo uld cooling. Finally, the optimal cooling design and optimal injection condition were obtained.展开更多
Fluid-solid interaction problems have been studied q uite extensively in the past years. Rotor-bearing system is a typical example. Fluid field is changed under the exciting of rotor vibration. On the same ti me, a ne...Fluid-solid interaction problems have been studied q uite extensively in the past years. Rotor-bearing system is a typical example. Fluid field is changed under the exciting of rotor vibration. On the same ti me, a net force caused by fluid pressure exerts on rotor, which will change roto r vibration. So, the fluid-solid coupled analysis method must be used. Traditionally, numerical difference method was used to solve fluid problems. The coupled fluid-solid equation could not be set up based on the method. It is no t until finite element method was used in fluid dynamics area then can the coupl ed dynamics be researched. Recently many experimental, analytical and numerical studies have been used in the area . But in these investigations, it is a ssumed that the solid vibration could not be influenced by fluid. In the other w ords, the force exerted on solid from fluid was neglected in the papers. So, the models built were some kinds of semi-coupled model only. In this paper, the Galerkin finite-element method, two-dimension vibration equ ation of rigid body and Navier-Stokes equations are used to build a full-coupl ed fluid-solid model in rotor-bearing system. Some assumptions are taken: 1) In fluid equation, the nonlinear terms are relatively small and neglected. 2) The gravity takes no effect on this system. 3) The bearing and the rotor are long. Flow and leakage along the axis is neglec ted. 4) The fluid is a kind of Newtonian incondensable viscous fluid. 5) The rotor is considered to be a rigid body. Using the model established, we calculated all the examples given by paper , results show the error are less than 7%. So the full-coupled model is built c orrectly. Examples are given in the end of the paper. After analyzing the examples, we get some conclusions: 1) In rotor-bearing system, while being taken under two conditions that whether coupled method is taken or not, difference of pressure and vibration amplitude could reach 76% and 120%. Therefore coupled method must be taken to investigate fluid-solid system. 1) Amplitude of fluid pressure can be more or less influenced by rotor unbalance , gap, eccentricity and other factors. 2) By using coupling method, results show that the amplitudes of vibration and p ressure are greater than ignoring the method. It should be paid more attention t o.展开更多
A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element...A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.展开更多
The extended finite element method(XFEM) is a numerical method for modeling discontinuities within the classical finite element framework. The computation mesh in XFEM is independent of the discontinuities, such that ...The extended finite element method(XFEM) is a numerical method for modeling discontinuities within the classical finite element framework. The computation mesh in XFEM is independent of the discontinuities, such that remeshing for moving discontinuities can be overcome. The extended finite element method is presented for hydro-mechanical modeling of impermeable discontinuities in rock. The governing equation of XFEM for hydraulic fracture modeling is derived by the virtual work principle of the fracture problem considering the water pressure on crack surface. The coupling relationship between water pressure gradient on crack surface and fracture opening width is obtained by semi-analytical and semi-numerical method. This method simplifies coupling analysis iteration and improves computational precision. Finally, the efficiency of the proposed method for modeling hydraulic fracture problems is verified by two examples and the advantages of the XFEM for hydraulic fracturing analysis are displayed.展开更多
Texture evolution and inhomogeneous deformation of polycrystalline Cu during uniaxial compression are investigated at the grain scale by combining crystal plasticity finite element method(CPFEM) with particle swarm op...Texture evolution and inhomogeneous deformation of polycrystalline Cu during uniaxial compression are investigated at the grain scale by combining crystal plasticity finite element method(CPFEM) with particle swarm optimization(PSO) algorithm. The texture-based representative volume element(TBRVE) is used in the crystal plasticity finite element model, where a given number of crystallographic orientations are obtained by means of discretizing the orientation distribution function(ODF) based on electron backscattered diffraction(EBSD) experiment data. Three-dimensional grains with different morphologies are generated on the basis of Voronoi tessellation. The PSO algorithm plays a significant role in identifying the material parameters and saving computational time. The macroscopic stress–strain curve is predicted based on CPFEM, where the simulation results are in good agreement with the experimental ones. Therefore, CPFEM is a powerful candidate for capturing the texture evolution and clarifying the inhomogeneous plastic deformation of polycrystalline Cu. The simulation results indicate that the <110> fiber texture is generated finally with the progression of plastic deformation. The inhomogeneous distribution of rotation angles lays the foundation for the inhomogeneous deformation of polycrystalline Cu in terms of grain scale.展开更多
The aim of this work is to propose a 3D FE model of a theoretical assembling straight bevel gear pair to analyze the contact fatigue on the tooth surface and the bending fatigue in the tooth root. Based on the cumulat...The aim of this work is to propose a 3D FE model of a theoretical assembling straight bevel gear pair to analyze the contact fatigue on the tooth surface and the bending fatigue in the tooth root. Based on the cumulative fatigue criterion and the stress-life equation, the key meshing states of the gear pair were investigated for the contact fatigue and the bending fatigue. Then, the reliability of the proposed model was proved by comparing the calculation result with the simulation result. Further study was performed to analyze the variation of the contact fatigue stress and the bending fatigue stress under different loads. Furthermore, the roles of the driving pinion and the driven gear pair were evaluated in the fatigue life of the straight bevel gear pair and the main fatigue failure mode was determined for the significant gear. The results show that the fatigue failure of the driving pinion is the main fatigue failure for the straight bevel gear pair and the bending fatigue failure is the main fatigue failure for the driving pinion.展开更多
The motor and trailer cars of a high-speed train were modeled as a multi-rigid body system with two suspensions. According to structural characteristic of a slab track, a new spatial vibration model of track segment e...The motor and trailer cars of a high-speed train were modeled as a multi-rigid body system with two suspensions. According to structural characteristic of a slab track, a new spatial vibration model of track segment element of the slab track was put forward. The spatial vibration equation set of the high-speed train and slab track system was then established on the basis of the principle of total potential energy with stationary value in elastic system dynamics and the rule of ＂set-in-right-position＂ for formulating system matrices. The equation set was solved by the Wilson-θ direct integration method. The contents mentioned above constitute the analysis theory of spatial vibration of high-speed train and slab track system. The theory was then verified by the high-speed running experiment carried out on the slab track in the Qinghuangdao-Shenyang passenger transport line. The results show that the calculated results agree well with the measured rcsults, such as the calculated lateral and vertical rail displacements are 0.82 mm and 0.9 mm and the measured ones 0.75 mm and 0.93 mm, respectively; the calculated lateral and vertical wheel-rail forces are 8.9 kN and 102.3 kN and the measured ones 8.6 kN and 80.2 kN, respectively. The interpolation method, that is, the lateral finite strip and slab segment element, for slab deformation proposed is of simplification and applicability compared with the traditional plate element method. All of these demonstrate the reliability of the theory proposed.展开更多
Facing the problems lack of considering the non-uniform distribution of the static bias magnetic field and computing the panicle displacements in the simulation model of electromagnetic acoustic transducer （EMAT）, a...Facing the problems lack of considering the non-uniform distribution of the static bias magnetic field and computing the panicle displacements in the simulation model of electromagnetic acoustic transducer （EMAT）, a multi-field coupled model was established and the finite element method （FEM） was presented to calculate the entire transduction process. The multi-field coupled model included the static magnetic field, pulsed eddy current field and mechanical field. The FEM equations of the three fields were derived by Garlerkin FEM method. Thus, the entire transduction process of the EMAT was calculated through sequentially coupling the three fields. The transduction process of a Lamb wave EMAT was calculated according to the present model and method. The results show that, by the present method, it is valid to calculate the particle displacement under the given excitation signal and non-uniformly distributed static magnetic field. Calculation error will be brought about if the non-uniform distribution of the static bias magnetic field is neglected.展开更多
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.展开更多
In-service hydrocarbons must be transported at high temperature and high pressure to ease the flow and prevent the solidification of the wax fraction. The pipeline containing hot oil will expand longitudinally due to ...In-service hydrocarbons must be transported at high temperature and high pressure to ease the flow and prevent the solidification of the wax fraction. The pipeline containing hot oil will expand longitudinally due to the rise in temperature. If such expansion is resisted, for example by frictional effects over a kilometer or so of pipeline, compressive axial stress will be built up in the pipe-wall. The compressive forces are often so large that they induce vertical buckling of buffed pipelines, which can jeopardize the structural integrity of the pipeline. A typical initial imperfection named continuous support mode of submarine pipeline was studied. Based on this type of initial imperfection, the analytical solution of vertical thermal buckling was introduced and an elastic-plasticity finite element analysis （FEA） was developed. Both the analytical and the finite element methodology were applied to analyze a practice in Bohai Gulf, China. The analyzing results show that upheaval buckling is most likely to build up from the initial imperfection of the pipeline and the buckling temperature depends on the amplitude of initial imperfection. With the same amplitude of initial imperfection, the triggering temperature difference of upheaval buckling increases with covered depth of the pipeline, the soil strength and the friction between the pipeline and subsoil.展开更多
The finite element analysis （FEA） software Ansys was employed to study the stress state of the dies of both plane and non-plane parting face structures with uniform interference and the die of plane parting face str...The finite element analysis （FEA） software Ansys was employed to study the stress state of the dies of both plane and non-plane parting face structures with uniform interference and the die of plane parting face structure with non-uniform interference. Considering the symmetry of the die, a half gear tooth model of the two-ring assembled die with 2.5 GPa inner pressure was constructed. Four paths were defined to investigate the stress state at the bottom comer of the die where stress concentration was serious. FEA results show that the change of parting face from non-plane to plane can greatly reduce the stress at the teeth tips of the die so that the tip fracture is avoided. The interference structure of the die is the most important influencing factor for the stress concentration at the bottom comer. When non-uniform interference is adopted the first principal stress at the comer on the defined paths of the die is much lower than that with uniform interference. The bottom hole radius is another important influencing factor for the comer stress concentration. The first principal stress at the comer of the plane parting face die with non-uniform interference is reduced from 2.3 to 1.9 GPa when the hole radius increases from 12.5 to 16.0 mm. The optimization of the die structure increases the life of the die from 100 to 6 000 hits.展开更多
A new method for optimizing a butterfly-shaped linear ultrasonic motor was proposed to maximize its mechanical output. The finite element analysis technology and response surface methodology were combined together to ...A new method for optimizing a butterfly-shaped linear ultrasonic motor was proposed to maximize its mechanical output. The finite element analysis technology and response surface methodology were combined together to realize the optimal design of the butterfly-shaped linear ultrasonic motor. First, the operation principle of the motor was introduced. Second, the finite element parameterized model of the stator of the motor was built using ANSYS parametric design language and some structure parameters of the stator were selected as design variables. Third, the sample points were selected in design variable space using latin hypercube Design. Through modal analysis and harmonic response analysis of the stator based on these sample points, the target responses were obtained. These sample points and response values were combined together to build a response surface model. Finally, the simplex method was used to find the optimal solution. The experimental results showed that many aspects of the design requirements of the butterfly-shaped linear ultrasonic motor have been fulfilled. The prototype motor fabricated based on the optimal design result exhibited considerably high dynamic performance, such as no-load speed of 873 ram/s, maximal thrust of 27.5 N, maximal efficiency of 43%, and thrust-weight ratio of 45.8.展开更多
Efficient bolted joint design is an essential part of designing the minimum weight aerospace structures, since structural failures usually occur at connections and interface. A comprehensive numerical study of three-d...Efficient bolted joint design is an essential part of designing the minimum weight aerospace structures, since structural failures usually occur at connections and interface. A comprehensive numerical study of three-dimensional(3D) stress variations is prohibitively expensive for a large-scale structure where hundreds of bolts can be present. In this work, the hybrid composite-to-metal bolted connections used in the upper stage of European Ariane 5ME rocket are analyzed using the global-local finite element(FE) approach which involves an approximate analysis of the whole structure followed by a detailed analysis of a significantly smaller region of interest. We calculate the Tsai-Wu failure index and the margin of safety using the stresses obtained from ABAQUS. We find that the composite part of a hybrid bolted connection is prone to failure compared to the metal part. We determine the bolt preload based on the clamp-up load calculated using a maximum preload to make the composite part safe. We conclude that the unsuitable bolt preload may cause the failure of the composite part due to the high stress concentration in the vicinity of the bolt. The global-local analysis provides an efficient computational tool for enhancing 3D stress analysis in the highly loaded region.展开更多
A three-dimensional finite element simulation was carried out to investigate the effects of tunnel construction on nearby pile foundation.The displacement controlled model (DCM) was used to simulate the tunneling-indu...A three-dimensional finite element simulation was carried out to investigate the effects of tunnel construction on nearby pile foundation.The displacement controlled model (DCM) was used to simulate the tunneling-induced volume loss effects.The numerical model was verified based on the results of a centrifuge test and a set of parametric studies was implemented based on this model.There is good agreement between the trend of the results of the centrifuge test and the present model.The results of parametric studies show that the tunnelling-induced pile internal force and deformation depend mainly on the pile?tunnel distance,the pile length to tunnel depth ratio and the volume loss.Two different zones are separated by a 45° line projected from the tunnel springline.Within the zone of influence,the pile is subjected to tensile force and large settlement;whereas outside the zone of influence,dragload and small settlement are induced.It is also established that the impact of tunnelling on a pile group is substantially smaller as compared with a single pile in the same location with the rear pile in a group,demonstrating a positive pile group effect.展开更多
Crack is found to be a major distress that affects the performance of the epoxy asphalt pavement.An extended finite element method was proposed for investigating the fracture properties of the epoxy asphalt mixture.Fi...Crack is found to be a major distress that affects the performance of the epoxy asphalt pavement.An extended finite element method was proposed for investigating the fracture properties of the epoxy asphalt mixture.Firstly,the single-edge notched beam test was used to analyze the temperature effect and calculate the material parameters.Then,the mechanical responses were studied using numerical analysis.It is concluded that 5℃ can be selected as the critical temperature that affects the fracture properties,and numerical simulations indicate that crack propagation is found to significantly affect the stress state of the epoxy asphalt mixture.The maximum principal stress at the crack surface exhibits different trends at various temperatures.Numerical solution of stress intensity factor can well meet the theoretical solution,especially when the temperature is lower than 5℃.展开更多
文摘Recent advancements in additive manufacturing(AM)have revolutionized the design and production of complex engineering microstructures.Despite these advancements,their mathematical modeling and computational analysis remain significant challenges.This research aims to develop an effective computational method for analyzing the free vibration of functionally graded(FG)microplates under high temperatures while resting on a Pasternak foundation(PF).This formulation leverages a new thirdorder shear deformation theory(new TSDT)for improved accuracy without requiring shear correction factors.Additionally,the modified couple stress theory(MCST)is incorporated to account for sizedependent effects in microplates.The PF is characterized by two parameters including spring stiffness(k_(w))and shear layer stiffness(k_(s)).To validate the proposed method,the results obtained are compared with those of the existing literature.Furthermore,numerical examples explore the influence of various factors on the high-temperature free vibration of FG microplates.These factors include the length scale parameter(l),geometric dimensions,material properties,and the presence of the elastic foundation.The findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the results of this research will have great potential in military and defense applications such as components of submarines,fighter aircraft,and missiles.
文摘At the first time,the finite element method was used to model and analyze the free vibration and transient response of non-uniform thickness bi-directional functionally graded sandwich porous(BFGSP)skew plates.The whole BFGSP skew-plates is placed on a variable visco-elastic foundation(VEF)in the hygro-thermal environment and subjected to the blast load.The BFGSP skew-plate thickness is permitted to vary non-linearly over both the length and width of the skew-plate,thereby faithfully representing the real behavior of the structure itself.The analysis is based on a four-node planar quadrilateral element with eight degrees of freedom per node,which is approximated using Lagrange Q_(4)shape function and C^(1)level non-conforming Hermite shape function based on refined higher-order shear deformation plate theory.The forced vibration parameters of the non-uniform thickness BFGSP skew-plate are fully determined using Hamilton's principle and the Newmark-βdirect integration technique.Accuracy of the calculation program is validated by comparing its numerical results with those from reputable sources.Furthermore,a thorough assessment is conducted to determine the impact of various parameters on the free and forced vibration responses of the non-uniform thickness BFGSP skew-plate.The findings of the paper may be used in the development of civil and military structures in situations that are prone to exceptional forces,such as explosions and impacts load.
基金Projects(40974077,41164004)supported by the National Natural Science Foundation of ChinaProject(2007AA06Z134)supported by the National High Technology Research and Development Program of China+2 种基金Projects(2011GXNSFA018003,0832263)supported by the Natural Science Foundation of Guangxi Province,ChinaProject supported by Program for Excellent Talents in Guangxi Higher Education Institution,ChinaProject supported by the Foundation of Guilin University of Technology,China
文摘The strategies that minimize the overall solution time of multiple linear systems in 3D finite element method (FEM) modeling of direct current (DC) resistivity were discussed. A global stiff matrix is assembled and stored in two parts separately. One part is associated with the volume integral and the other is associated with the subsurface boundary integral. The equivalent multiple linear systems with closer right-hand sides than the original systems were constructed. A recycling Krylov subspace technique was employed to solve the multiple linear systems. The solution of the seed system was used as an initial guess for the subsequent systems. The results of two numerical experiments show that the improved algorithm reduces the iterations and CPU time by almost 50%, compared with the classical preconditioned conjugate gradient method.
文摘In an injection moulding process, the mould cooling s ystem is very important as an efficient and uniform cooling effect can improve b oth the productivity and part quality. Due to the complexity of the process, muc h research on the mould cooling analysis and cooling design optimization has bee n focused on the core and the cavity, excluding other systems of injection mould s. However, the runner system introduces a considerable amount of heat into the mould. In recent years, more and more hot runner systems are being applied in th e moulding industry to save material and decrease losses of injection pressure. This raises the need to include the hot runner system in the cooling analysis. I n this paper, a photo frame part was studied. The mould was built with a hot run ner system. Two thermal sensors were installed: one measures the temperature of lateral surface of hot runner nozzle; the other measures the plastic temperature from the core side. A pressure sensor was also installed to measure the pressur e of the core impression. Cooling analysis was performed using ABAQUS, ananalysi s software based on the Finite Element Method (FEM). The assembly including core , cavity and plastic part was modeled. Heat conduction from hot runner to cavity and from polymer melt to the mould and force convection on the cooling channel surfaces were studied. The natural convection between the ambient air and the ex terior mould surface was ignored. The simulations were adjusted with the experim ental results to find out the heat input from hot runner and its influence on mo uld cooling. Finally, the optimal cooling design and optimal injection condition were obtained.
文摘Fluid-solid interaction problems have been studied q uite extensively in the past years. Rotor-bearing system is a typical example. Fluid field is changed under the exciting of rotor vibration. On the same ti me, a net force caused by fluid pressure exerts on rotor, which will change roto r vibration. So, the fluid-solid coupled analysis method must be used. Traditionally, numerical difference method was used to solve fluid problems. The coupled fluid-solid equation could not be set up based on the method. It is no t until finite element method was used in fluid dynamics area then can the coupl ed dynamics be researched. Recently many experimental, analytical and numerical studies have been used in the area . But in these investigations, it is a ssumed that the solid vibration could not be influenced by fluid. In the other w ords, the force exerted on solid from fluid was neglected in the papers. So, the models built were some kinds of semi-coupled model only. In this paper, the Galerkin finite-element method, two-dimension vibration equ ation of rigid body and Navier-Stokes equations are used to build a full-coupl ed fluid-solid model in rotor-bearing system. Some assumptions are taken: 1) In fluid equation, the nonlinear terms are relatively small and neglected. 2) The gravity takes no effect on this system. 3) The bearing and the rotor are long. Flow and leakage along the axis is neglec ted. 4) The fluid is a kind of Newtonian incondensable viscous fluid. 5) The rotor is considered to be a rigid body. Using the model established, we calculated all the examples given by paper , results show the error are less than 7%. So the full-coupled model is built c orrectly. Examples are given in the end of the paper. After analyzing the examples, we get some conclusions: 1) In rotor-bearing system, while being taken under two conditions that whether coupled method is taken or not, difference of pressure and vibration amplitude could reach 76% and 120%. Therefore coupled method must be taken to investigate fluid-solid system. 1) Amplitude of fluid pressure can be more or less influenced by rotor unbalance , gap, eccentricity and other factors. 2) By using coupling method, results show that the amplitudes of vibration and p ressure are greater than ignoring the method. It should be paid more attention t o.
基金Project supported by the National Basic Research Program of China (973Project) (No.2002CB412709) and the National Natural Science Foundation of China (Nos.50278012,10272027,19832010)
文摘A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.
基金Project(2011CB013505)supported by the National Basic Research Program of ChinaProject(51279100)supported by the National Natural Science Foundation of China
文摘The extended finite element method(XFEM) is a numerical method for modeling discontinuities within the classical finite element framework. The computation mesh in XFEM is independent of the discontinuities, such that remeshing for moving discontinuities can be overcome. The extended finite element method is presented for hydro-mechanical modeling of impermeable discontinuities in rock. The governing equation of XFEM for hydraulic fracture modeling is derived by the virtual work principle of the fracture problem considering the water pressure on crack surface. The coupling relationship between water pressure gradient on crack surface and fracture opening width is obtained by semi-analytical and semi-numerical method. This method simplifies coupling analysis iteration and improves computational precision. Finally, the efficiency of the proposed method for modeling hydraulic fracture problems is verified by two examples and the advantages of the XFEM for hydraulic fracturing analysis are displayed.
基金Projects(51305091,51475101) supported by the National Natural Science Foundation of ChinaProject(20132304120025) supported by Specialized Research Fund for the Doctoral Program of Higher Education,China
文摘Texture evolution and inhomogeneous deformation of polycrystalline Cu during uniaxial compression are investigated at the grain scale by combining crystal plasticity finite element method(CPFEM) with particle swarm optimization(PSO) algorithm. The texture-based representative volume element(TBRVE) is used in the crystal plasticity finite element model, where a given number of crystallographic orientations are obtained by means of discretizing the orientation distribution function(ODF) based on electron backscattered diffraction(EBSD) experiment data. Three-dimensional grains with different morphologies are generated on the basis of Voronoi tessellation. The PSO algorithm plays a significant role in identifying the material parameters and saving computational time. The macroscopic stress–strain curve is predicted based on CPFEM, where the simulation results are in good agreement with the experimental ones. Therefore, CPFEM is a powerful candidate for capturing the texture evolution and clarifying the inhomogeneous plastic deformation of polycrystalline Cu. The simulation results indicate that the <110> fiber texture is generated finally with the progression of plastic deformation. The inhomogeneous distribution of rotation angles lays the foundation for the inhomogeneous deformation of polycrystalline Cu in terms of grain scale.
基金Project(51105287) supported by the National Natural Science Foundation of ChinaProject(2012BAA08003) supported by the Key Research and Development Project of New Products and New Technologies of Hubei Province, ChinaProject(2011-P05) supported by the State Key Laboratory of Materials Processing and Die & Mould Technology, Huazhong University of Science and Technology,China
文摘The aim of this work is to propose a 3D FE model of a theoretical assembling straight bevel gear pair to analyze the contact fatigue on the tooth surface and the bending fatigue in the tooth root. Based on the cumulative fatigue criterion and the stress-life equation, the key meshing states of the gear pair were investigated for the contact fatigue and the bending fatigue. Then, the reliability of the proposed model was proved by comparing the calculation result with the simulation result. Further study was performed to analyze the variation of the contact fatigue stress and the bending fatigue stress under different loads. Furthermore, the roles of the driving pinion and the driven gear pair were evaluated in the fatigue life of the straight bevel gear pair and the main fatigue failure mode was determined for the significant gear. The results show that the fatigue failure of the driving pinion is the main fatigue failure for the straight bevel gear pair and the bending fatigue failure is the main fatigue failure for the driving pinion.
基金Project(2007CB714706) supported by the National Basic Research Program of ChinaProject (50678176) supported by the National Natural Science Foundation of ChinaProject(NCET-07-0866) supported by the Program for New Century Excellent Talents in University
文摘The motor and trailer cars of a high-speed train were modeled as a multi-rigid body system with two suspensions. According to structural characteristic of a slab track, a new spatial vibration model of track segment element of the slab track was put forward. The spatial vibration equation set of the high-speed train and slab track system was then established on the basis of the principle of total potential energy with stationary value in elastic system dynamics and the rule of ＂set-in-right-position＂ for formulating system matrices. The equation set was solved by the Wilson-θ direct integration method. The contents mentioned above constitute the analysis theory of spatial vibration of high-speed train and slab track system. The theory was then verified by the high-speed running experiment carried out on the slab track in the Qinghuangdao-Shenyang passenger transport line. The results show that the calculated results agree well with the measured rcsults, such as the calculated lateral and vertical rail displacements are 0.82 mm and 0.9 mm and the measured ones 0.75 mm and 0.93 mm, respectively; the calculated lateral and vertical wheel-rail forces are 8.9 kN and 102.3 kN and the measured ones 8.6 kN and 80.2 kN, respectively. The interpolation method, that is, the lateral finite strip and slab segment element, for slab deformation proposed is of simplification and applicability compared with the traditional plate element method. All of these demonstrate the reliability of the theory proposed.
基金Project(10974115) supported by the National Natural Science Foundation of China
文摘Facing the problems lack of considering the non-uniform distribution of the static bias magnetic field and computing the panicle displacements in the simulation model of electromagnetic acoustic transducer （EMAT）, a multi-field coupled model was established and the finite element method （FEM） was presented to calculate the entire transduction process. The multi-field coupled model included the static magnetic field, pulsed eddy current field and mechanical field. The FEM equations of the three fields were derived by Garlerkin FEM method. Thus, the entire transduction process of the EMAT was calculated through sequentially coupling the three fields. The transduction process of a Lamb wave EMAT was calculated according to the present model and method. The results show that, by the present method, it is valid to calculate the particle displacement under the given excitation signal and non-uniformly distributed static magnetic field. Calculation error will be brought about if the non-uniform distribution of the static bias magnetic field is neglected.
基金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.
基金Project(51021004) supported by Innovative Research Groups of the National Natural Science Foundation of ChinaProject(40776055) supported by the National Natural Science Foundation of china+1 种基金Project(1002) supported by State Key Laboratory of Ocean Engineering Foundation, ChinaProject(NCET 11 0370) supported by the Program for New Century Excellent Talents in Universities of China
文摘In-service hydrocarbons must be transported at high temperature and high pressure to ease the flow and prevent the solidification of the wax fraction. The pipeline containing hot oil will expand longitudinally due to the rise in temperature. If such expansion is resisted, for example by frictional effects over a kilometer or so of pipeline, compressive axial stress will be built up in the pipe-wall. The compressive forces are often so large that they induce vertical buckling of buffed pipelines, which can jeopardize the structural integrity of the pipeline. A typical initial imperfection named continuous support mode of submarine pipeline was studied. Based on this type of initial imperfection, the analytical solution of vertical thermal buckling was introduced and an elastic-plasticity finite element analysis （FEA） was developed. Both the analytical and the finite element methodology were applied to analyze a practice in Bohai Gulf, China. The analyzing results show that upheaval buckling is most likely to build up from the initial imperfection of the pipeline and the buckling temperature depends on the amplitude of initial imperfection. With the same amplitude of initial imperfection, the triggering temperature difference of upheaval buckling increases with covered depth of the pipeline, the soil strength and the friction between the pipeline and subsoil.
基金Project(2006BAF04B06) supported by the National Key Technology R & D Program of ChinaProject(2005AA101B19) supported by the Key Technology R & D Program of Hubei Province, China
文摘The finite element analysis （FEA） software Ansys was employed to study the stress state of the dies of both plane and non-plane parting face structures with uniform interference and the die of plane parting face structure with non-uniform interference. Considering the symmetry of the die, a half gear tooth model of the two-ring assembled die with 2.5 GPa inner pressure was constructed. Four paths were defined to investigate the stress state at the bottom comer of the die where stress concentration was serious. FEA results show that the change of parting face from non-plane to plane can greatly reduce the stress at the teeth tips of the die so that the tip fracture is avoided. The interference structure of the die is the most important influencing factor for the stress concentration at the bottom comer. When non-uniform interference is adopted the first principal stress at the comer on the defined paths of the die is much lower than that with uniform interference. The bottom hole radius is another important influencing factor for the comer stress concentration. The first principal stress at the comer of the plane parting face die with non-uniform interference is reduced from 2.3 to 1.9 GPa when the hole radius increases from 12.5 to 16.0 mm. The optimization of the die structure increases the life of the die from 100 to 6 000 hits.
基金Projects(51275235, 50975135) supported by the National Natural Science Foundation of ChinaProject(U0934004) supported by the Natural Science Foundation of Guangdong Province, ChinaProject(2011CB707602) supported by the National Basic Research Program of China
文摘A new method for optimizing a butterfly-shaped linear ultrasonic motor was proposed to maximize its mechanical output. The finite element analysis technology and response surface methodology were combined together to realize the optimal design of the butterfly-shaped linear ultrasonic motor. First, the operation principle of the motor was introduced. Second, the finite element parameterized model of the stator of the motor was built using ANSYS parametric design language and some structure parameters of the stator were selected as design variables. Third, the sample points were selected in design variable space using latin hypercube Design. Through modal analysis and harmonic response analysis of the stator based on these sample points, the target responses were obtained. These sample points and response values were combined together to build a response surface model. Finally, the simplex method was used to find the optimal solution. The experimental results showed that many aspects of the design requirements of the butterfly-shaped linear ultrasonic motor have been fulfilled. The prototype motor fabricated based on the optimal design result exhibited considerably high dynamic performance, such as no-load speed of 873 ram/s, maximal thrust of 27.5 N, maximal efficiency of 43%, and thrust-weight ratio of 45.8.
基金Project(282522)supported by the European Union's Research and Innovation Funding Programme
文摘Efficient bolted joint design is an essential part of designing the minimum weight aerospace structures, since structural failures usually occur at connections and interface. A comprehensive numerical study of three-dimensional(3D) stress variations is prohibitively expensive for a large-scale structure where hundreds of bolts can be present. In this work, the hybrid composite-to-metal bolted connections used in the upper stage of European Ariane 5ME rocket are analyzed using the global-local finite element(FE) approach which involves an approximate analysis of the whole structure followed by a detailed analysis of a significantly smaller region of interest. We calculate the Tsai-Wu failure index and the margin of safety using the stresses obtained from ABAQUS. We find that the composite part of a hybrid bolted connection is prone to failure compared to the metal part. We determine the bolt preload based on the clamp-up load calculated using a maximum preload to make the composite part safe. We conclude that the unsuitable bolt preload may cause the failure of the composite part due to the high stress concentration in the vicinity of the bolt. The global-local analysis provides an efficient computational tool for enhancing 3D stress analysis in the highly loaded region.
文摘A three-dimensional finite element simulation was carried out to investigate the effects of tunnel construction on nearby pile foundation.The displacement controlled model (DCM) was used to simulate the tunneling-induced volume loss effects.The numerical model was verified based on the results of a centrifuge test and a set of parametric studies was implemented based on this model.There is good agreement between the trend of the results of the centrifuge test and the present model.The results of parametric studies show that the tunnelling-induced pile internal force and deformation depend mainly on the pile?tunnel distance,the pile length to tunnel depth ratio and the volume loss.Two different zones are separated by a 45° line projected from the tunnel springline.Within the zone of influence,the pile is subjected to tensile force and large settlement;whereas outside the zone of influence,dragload and small settlement are induced.It is also established that the impact of tunnelling on a pile group is substantially smaller as compared with a single pile in the same location with the rear pile in a group,demonstrating a positive pile group effect.
基金Project(50578038)supported by the National Natural Science Foundation of China
文摘Crack is found to be a major distress that affects the performance of the epoxy asphalt pavement.An extended finite element method was proposed for investigating the fracture properties of the epoxy asphalt mixture.Firstly,the single-edge notched beam test was used to analyze the temperature effect and calculate the material parameters.Then,the mechanical responses were studied using numerical analysis.It is concluded that 5℃ can be selected as the critical temperature that affects the fracture properties,and numerical simulations indicate that crack propagation is found to significantly affect the stress state of the epoxy asphalt mixture.The maximum principal stress at the crack surface exhibits different trends at various temperatures.Numerical solution of stress intensity factor can well meet the theoretical solution,especially when the temperature is lower than 5℃.
