Numerical simulation of contact problems in vane machinery by a parametric quadratic programming method

Interference fits are widely used for connecting impeller and shaft assembly that are forced together slowly by pressing. The interference fit design ensures stable balance behavior and allows for positive contact between the impeller and shaft assembly throughout the range of operating speeds. In addition to maintaining radial contact, sufficient net radial interface pressure must remain in order to transmit torque when the rotational speed is very high. Therefore, the interference fit between the impeller and the shaft assembly is one of the most important factors influencing the performance of the turbo unit in the design of turbocharger compressor. A suitable fit tolerance needs to be considered in the structural design.A locomotive-type turbocharger compressor with 24 blades under combined centrifugal and interference fit loading is used for the analysis. The finite-element (FE) parametric quadratic programming (PQP) method developed based on the parametric variational principle (PVP) is used for the analysis of the stress distribution in the three-dimensional (3D) contact problem of impeller. The advantages of the parametric programming method compared with conventional approaches are that the penalty factors can be canceled and that solutions can be obtained directly without tedious iterative procedures such as the general incremental iterative method. To save time in the computation, a~multi-substructure technique is adopted for structural modeling. This not only simplifies the calculation, but also provides a convenient service for process computer-aided design (CAD) by means of FE simulation. The effects of the fit tolerance, coefficient of friction and rotational speed (centrifugal force), wall thickness of the shaft sleeve and the contact stress on the interference-fitting surfaces are studied in detail in the numerical computation. It is found that a nonuniform initial amount of interference in the structural design avoids the relative displacement generated and ensures uniformity of the contact stress. To assure quality of press-fitting, the amount of interference between the shaft sleeve and shaft should be strictly controlled to avoid the rapid increase of the contact stress. The numerical results demonstrate the high accuracy and good convergence of the algorithm presented here, which provides an effective approach that achieves more-reliable interference-fitted connections and more-precise assembly accuracy with lower manufacturing cost in the structural design.

[1]  A. Klarbring A mathematical programming approach to three-dimensional contact problems with friction , 1986 .

[2]  A. B. Turner,et al.  Development of a novel gas turbine driven centrifugal compressor , 2000 .

[3]  X. D. Fang,et al.  Design of interference fits via finite element method , 2000 .

[4]  Eugenio Dragoni,et al.  Optimal aspect ratio of interference fits for maximum load transfer capacity , 2005 .

[5]  Zhang Hong-wu,et al.  A combined parametric quadratic programming and iteration method for 3-D elastic-plastic frictional contact problem analysis , 1998 .

[6]  Peter Wriggers,et al.  A new algorithm for numerical solution of 3D elastoplastic contact problems with orthotropic friction law , 2004 .

[7]  Øyvind Bjørke,et al.  Computer Aided Tolerancing , 1989 .

[8]  Gu Yuanxian,et al.  A combined programming and iteration algorithm for finite element analysis of three-dimensional contact problems , 1995 .

[9]  Yi Qin,et al.  A new approach for the optimisation of the shrink-fitting of cold-forging dies , 2004 .

[10]  Stephen J. Wright Primal-Dual Interior-Point Methods , 1997, Other Titles in Applied Mathematics.

[11]  Bunyamin Aksakal,et al.  Stress analysis of interference fitted shaft–hub system under transient heat transfer conditions , 2004 .

[12]  William J.T. Daniel FLYWHEEL DESIGN BY THE FINITE ELEMENT METHOD , 1981 .

[13]  D. W. Meyer,et al.  A contact friction algorithm including nonlinear viscoelasticity and a singular yield surface provision , 1992 .

[14]  V. Ramamurti,et al.  Stress distribution in interference joints , 1994 .

[15]  Hong Wu Zhang,et al.  Quadratic Programming Method for Numerical Modeling of Elastoplastic Contact Problems: Theory, Software, Applications , 2004 .

[16]  Bernhard A. Schrefler,et al.  Global Constitutive Behavior of Periodic Assemblies of Inelastic Bodies in Contact , 2000 .

[17]  Y. Choi,et al.  A design method for cold backward extrusion using FE analysis , 2003 .

[18]  Masaru Nakazawa,et al.  Finite element incremental contact analysis with various frictional conditions , 1979 .

[19]  K. Sridhara,et al.  Numerical analysis and experimental verification of the radial growth of a turbocharger centrifugal compressor impeller , 1997 .

[20]  K. G. Murty,et al.  Complementarity problems , 2000 .

[21]  Changming Zhu,et al.  The solution of frictional contact problems using a finite element-mathematical programming method , 1994 .

[22]  E. A. Wilson,et al.  A Method for Determining the Surface Contact Stresses Resulting From Interference Fits , 1970 .

[23]  Joseph Edward Shigley,et al.  Standard Handbook of Machine Design , 2004 .

[24]  B. Dattaguru,et al.  Analysis of a cracked lug loaded by an interference fit pin , 1996 .