Published: 07 October 2016

The effect of dry friction forces on the process of dielectric wafer grinding

L. A. Igumnov1
V. S. Metrikin2
A. V. Grezina3
A. G. Panasenko4
1, 2Research Institute for mechanics, National Research Lobachevsky State University of Nizhni Novgorod, Nizhny Novgorod, Russia
3, 4National Research Lobachevsky State University of Nizhni Novgorod, Nizhny Novgorod, Russia
Corresponding Author:
V. S. Metrikin
Views 21
Reads 10
Downloads 1345


We have investigated possible motions of the holder on top of the polishing pad during the process of dielectric plate grinding taking into account the forces of dry friction about its axis. A mathematical model of the mechanical device has been elaborated to describe the process of dielectric wafer grinding. The model is in the form of a non-autonomous nonlinear system with a variable structure. The structure of the phase space of the dynamical system was investigated, the qualitative studies of the possible motion modes were carried out. The values of the geometrical and dynamic parameters that qualitatively and quantitatively influence the modes of the holder motion were obtained. It was found that the inclusion of dry friction forces on the axis of the free holder result in a periodic motion of the mechanism with long stops. We present the calculations of the parameters for the type 3PD-320 machine.

1. Introduction

The authors in [1, 2] focus on the study of the kinematics of the process and the quality of the abrasive machining of wafer surfaces. Here, we investigate the dynamics of the mechanical grinding of dielectric wafers pasted to a free holder, whose axes are subjected to frictional forces. Schematic diagram of the mechanism is shown in Fig. 1. The wafers to be ground are pasted to the bottom of the holder (a small circle with the center at O4), which is placed on the polishing pad (a large circle with its center at O3 ). The grinding is carried out by a complex relative motion of the holder and the polishing pad. The polishing pad rotates at angular velocity ω0. The motion of the holder’s axis (the center of O4) is determined by the operation of a crank-and-rod mechanism in which the crank O1B rotates at angular velocity ω. The crank O1B through the rod AB, which is hinged at points A and B, acts on the rigid structure AO2O4, so that the axis of the holder (O4) oscillates. The rotation of the holder is determined by the effect of the friction forces between the polishing pad and the wafers pasted to the holder, and by the friction torque on the axis of the holder.

Fig. 1Schematic diagram of the grinding mechanism

Schematic diagram of the grinding mechanism

1.1. Mathematical model

In the development of the mathematical model it is assumed that the relative material removal rate is proportional to the relative velocity of the piece being ground and the polishing pad (Preston’s law [3]).

Introducing the notation:

O1B=b, AB=c, AO2=d, O2O3=R0, O2O4=R1, O2O1=f, ψ=ωt,

and making a series of simple transformations, we obtain the formula for calculating the velocity Vr of the point M on the holder relative to the polishing pad in the form:



G=f2+d2+b2-c2,Φ1t=2fbsinψ, cosφ0=f2-c2+d22fd,
ρ=RR0, ρ1=1+ρρ-2cos(φ-φ0), ρ2=rR0,
tgϑ=f-bsinψbcosψ, tgμ=-sinφ0(ρ-cosφ-φ0)-1.

It follows from Preston’s law [3] that the relative rate of material removal at different points of the wafers on the holder S(r,γ)=Sρ2,γ is determined from the following Eq. (2):


Fig. 2 shows for f= 22.7; c= 22.5; d= 3.5; R0= 40; b= 0.4. in three-dimensional space [gammaγ,ro2ρ,V/R0(S)] the results of the calculations of the relative value of material removal on the wafers being ground depending on the coordinates of the wafer points: the angle γ, the radius ρ, ω at ω0=0.5c-1 for the mode of the holder fixed on its axis, using the MATLAB package. Let us now consider in more detail the process of grinding in the free mode of the holder. In this case, the holder’s rotation is determined by the action of the moment M1 of the frictional forces on the holder from the pad and the friction torque on the axis of the holder.

To determine Μ1, we choose a coordinate system x'y' parallel to the system xy with the origin that coincides with the center of the holder. In this system, the inertial forces acting on small parts of the holder are equal to ΔF=-ΔmW0 where Δm is a part of the mass, and W0 is the acceleration of the center of the holder. These forces do not generate torque, since the holder is disk-shaped and the wafers being ground are fixed, as a rule, symmetrically without changing the position of the center of gravity that coincides with the center of the disk. Thus we conclude that the moment Μ1 is only determined by the action of the friction forces between the holder (the wafers on the holder) and the polishing pad.

A detailed study of the friction forces associated with the process of polishing is a separate large task. Therefore, to simplify, we assume that the frictional force acting on the element of the holder with the area ds is equal to:


The admissibility of this representation of the frictional forces is based on the similarity of the abrasive with lubricant to a high viscosity liquid. Constant χ is proportional to its viscosity. Then we can write Eq. (4):


Fig. 2The dependence of the relative rate of material removal S on the coordinate points of the wafer

The dependence of the relative rate of material removal S on the coordinate points of the wafer

Going to the coordinate system ξ, ζ, η fixed relative to the holder (the axis ξ is directed along O4W, the axis ζ – along O4U, the axes ξ, ζ, η form a right-handed space) and calculating the vector value [r×νr] as:


it is possible to compute the moment relative to the axis η in the following way:


where ω3 is the angular speed of the holder in the coordinate system x and y. Here, r0 is the radius of the holder. Thus, in the chosen coordinate system we have:


where J is the moment of inertia of the holder with the wafers. Assume further that the force of friction on the axis is the Coulomb-Amontons friction and represent the moment in the following form:

ΜZ=Bsignω2, ω20, ΜZ-B,B, ω2=0.

Since ω3=ω1+ω2, where ω2 is the velocity of the holder relative to its axis, and ω1 is the angular velocity of the holder drive piece, which can be considered for the above scheme approximately equal to ω1=Asin(ωt+φ0), ω0 is the polishing pad velocity, the Eq. (6) can be rewritten as Eq. (7):

dω2dt+λω2=-Aωcosωt+φ0+λsinωt+φ0+λω0-ΜZω2J, λ=χπr042J.

Introducing the dimensionless time τ=ωt+φ0-Φ, coordinate x=(ω2-ω0)/ω0 and parameters k=λ/ω, F=Aλ2+ω2/(ωω0), a=B/(Jω0ω), Eq. (7) can be rewritten as follows Eq. (8):

x˙+kx=-Fcosτ-asignx+1, tgΦ=λω.

The phase space of Eq. (8) is two-dimensional in the coordinates (τ,x). It has a straight line Г (x=-1), which divides the plane into three subspaces X+(x>-1), X-(x<-1) and Г (x=-1). The motion of the image point in said subspaces is described by the Eq. (9):

x˙+kx=-Fcosτ-asignx+1, x-1, k+Fcosτa, x=-1.

Note that the motion of the image point in the subspace Г (x=-1) takes place in time intervals Δi[τi,τi+1] corresponding to the junction of the phase trajectories coming from the subspaces X-, X+.

Fig. 3 (a, b) shows the phase trajectories for different parameter values a= 9 (Fig. 3(a)) and a= 7 (Fig. 3(b)), respectively, at k= 0.71, F= 3.2.

Fig. 3Phase trajectories in the plane τ, x

Phase trajectories in the plane τ, x


Phase trajectories in the plane τ, x


What all the parameters have in common is the presence of parameter intervals Δi of the holder stops, which increase with a rise of the coefficient of static friction. The rotation with the stopping of the holder begins only from some values of a (see Fig. 3(b)). With the increase of the parameter k the intervals are shifted to the region of the shorter times, and the oscillation amplitude decreases. In the ideal case, when we can neglect friction on the axis of the holder (Μ2= 0), it follows from Eq. (6) that with the obvious relation (χr04π/2)/(Jω0)>0 we will obtain limtω3(t)=ω0. This means that the angular velocity of the holder ω2 relative to its axis is ω2=ω0-ω1(t). Note that for the chosen values of the parameters of the mechanism the relation ω1=φ˙0.25sin(ωt+φ0) holds true and hence ω2>0. Assuming Μ20 and that it can be approximated by the force of friction obeying the Coulomb-Amontons law, it follows from Eq. (6) that limtω3(t)=ω0-Δω, where Δω is a small constant, which is equal to Δω=2M2max/(χπr04). In this case, for the steady state mode, instead of the relationship Eq. (1) we will have the following Eq. (10):


2. Conclusions

Note that according to Eq. (10):

– In the ideal case (Μ2= 0) there is no dependence of the relative velocity Vr from the position of the point on the holder ρ2, that is, the removal of the material on the grinding disc will be uniform across the surface;

– The friction on the axis of the holder leads to a weak dependence Vr from ρ2 and thus to the non-uniform removal of the material;

– The dynamics of the mechanism is determined by the ratio of the friction torque on the axis and in the contact of the wafers with the polishing pad;

– There are three modes of the mechanism motion: free rotation, rotation with periodic stops of the holder; the holder does not rotate on its axis.


  • Khomich N. S., Lugovik А. Yu., Fedortsev Р. V., Korzun A. E., Kukhto P. V. Simulation of the kinematics of the process of magnetic abrasive polishing of wafers. BSTU Herald, Vol. 1, 2009, p. 33-38, (in Russian).
  • Gavrishuk Е. М., Komarovv N., Metrikin V. S., Panasenko А. G. Mathematical of wafer grinding on the machines of type 4PD-200 and 3PD-320. Bulletin of the Samara SC RАS, Vol. 13, Issue 1, 2011, p. 992-995, (in Russian).
  • Nanz G., Camille L. E. Modeling of chemical-mechanical polishing: a review. IEEE Transactions on Semiconductor Manufacturing, Vol. 8, Issue 4, 1995, p. 382-389.

About this article

25 May 2016
03 September 2016
07 October 2016
Mathematical models in engineering
mathematical model
phase space
regimes with a long stop

Development of mathematical model is financially supported by Russian Science Foundation, project No. 16-19-10237. Investigation of numerical results was financed within the framework of the base part of State Task of the Ministry of Education and Science Project, No. 2014/134 2226.