Page 3 - Cutting tool temperature prediction method using analytical model for end milling
P. 3

1790                                                                                  W. Baohai et al.




















                                                                        Fig. 2  Stress distribution on rake face.


            Fig. 1  End milling diagram and approximation heating model.
                                                                           Z  t  Z  L x  Z  L y
                                                                          a
                                                               Tðx; y; z; tÞ¼         Gðx; y; z; t; x p ; y ; 0; DÞ
                                                                                                 p
                                                                          k  0  0  0
            cutting action, the conductor consists of insert and cutter has
                                                                                 p
                                                                                      p
                                                                           Qðx p ; y ; sÞdy dx p ds         ð4Þ
            steady temperature. With the cutting action perform, it begins
            to be heated by heat source over insert surface. According to  The magnitude of heat flux is strongly affected by the fric-
            heat transfer, the conductor like this can be assumed as a  tion along the tool-chip interface. However, there are two dif-
            semi-infinite body in analysis temperature change. Therefore,  ferent friction states in this interface in the process of end
            as an approximation, the insert is modeled as a semi-infinite  milling, respectively called stick zone and slip zone, as shown
            rectangular corner heated by a time-varying heat flux on the  in Fig. 2.
            tool-chip interface, as shown in Fig. 1.              At the beginning of chip formation, the workpiece is
              In this paper, with the heat radiation and temperature effect  squeezed from the tool and in turn forms a chip. As a result,
            on the thermal properties of the tool material neglected, the  the shear stress e in the area that the chip closely contacts with
            governing equation and boundary condition for the tempera-  the tool tip is constant k chip . Moreover, there is very high nor-
            ture T of the insert are set as                    mal stress r in this area. The chip is able to separate from the
                                                               workpiece on the condition that the shear stress is increasingly
            (
               2    1  @T
              r T ¼                                           promoted. At this moment, the shear stress matches the normal
                    a  @t                                ð1Þ
                @T                                             stress by the friction coefficient l. In this study, the heat flux Q
                 @t
              k  ¼ Qðx; y; tÞ z ¼ 0; 0 6 x 6 L x ; 0 6 y 6 L y
                                                               is considered to come from the uniform heat source and the
            where a is the thermal diffusivity, t is the time, k the thermal  liner heat source based on the following empirical equations:
            conductivity, Q the heat flux, and L x and L y are the dimen-
                                                               Qðx; y; sÞ¼ q ðxÞqðsÞ                         ð5Þ
                                                                         1
            sions of the heat source. Except the interface between the tool
            and chip, the rest surface of tool is regarded as insulated and   1         0 6 x 6 bL x
            the initial temperature is set to be 0 °C.         q ðxÞ¼                                        ð6Þ
                                                                1
              The green function can be used to solve the above heat con-  ðL x  xÞ=½L x ð1  bފ  bL x 6 x 6 L x
                         5
            duction equation. The green function G for the temperature  where b is length ratio of stick zone over slip zone, q ðxÞ
                                                                                                            1
            of a semi-infinite corner excited by an instantaneous heat  reflects the type of heat source distribution, and qðsÞ is the
            source at time s at the surface point ðx ¼ x p ; y ¼ y ; z ¼ 0Þ is  source heat flux varying with time.
                                                   p
                                                                  Eqs. (5) and (6) then will be utilized to solve the tempera-
            Gðx; y; z; x p ; y ; 0; DÞ                         ture equation. First, the inner double integral of the green
                      p
                                "             !
                                           2
                  2        z 2       ðx þ x p Þ              function GR needs to be calculated:
              ¼   p ffiffiffi 3  exp  2   exp  2
                ðD pÞ      D             D                     GRðx; y; z; L x ; L y ; DÞ
                              !# "             !
                             2                2                    Z  L x  Z  L y
                     ðx  x p Þ       ðy þ y Þ
                                             p
               þ exp              exp                           ¼         Gðx; y; z; t; x p ; y ; 0; DÞdx p dy p
                                                                                       p
                         D 2              D 2                       0   0
                                                                                       "             !
                             2  !#                                    2       z 2   Z  L x  ðx þ x p Þ 2
                     ðy  y Þ                                   ¼        exp          exp
                           p
               þ exp                                     ð2Þ         p ffiffiffi 3   D 2              D 2
                        D 2                                        ð D pÞ            0
                                                                                 !#        "             !
                                                                                2       Z  L y          2
                                                                         ðx  x p Þ             ðy þ y Þ
                                                                                                       p
                                                                   þ exp                    exp
                p ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi                                           2      dx p            2
            D ¼ 2  aðt  sÞ                              ð3Þ                D            0          D
                                                                                2  !#
            where D is a parameter with units of length, and its physical  ðy  y Þ
                                                                               p
                                                                   þ exp           dy                        ð7Þ
            meaning is characteristic dimension for the penetration of tem-  D 2     p
            perature field at time t resulted from instantaneous heat source
            at time s. According to Osakis’ theory, 18  the solution of the  An error function erfðxÞ is introduced in the inner double
            temperature in a semi-infinite corner is            integral calculation:
   1   2   3   4   5   6   7