International Scholarly Research Notices

International Scholarly Research Notices / 2011 / Article

Research Article | Open Access

Volume 2011 |Article ID 309691 | https://doi.org/10.5402/2011/309691

A. H. El-Sinawi, "Modeling and Control of Cutter Vibration in the Presence of Random Distribution of Microhardness of Workpiece and Nonlinear Cutting Forces in Lathe Process", International Scholarly Research Notices, vol. 2011, Article ID 309691, 11 pages, 2011. https://doi.org/10.5402/2011/309691

Modeling and Control of Cutter Vibration in the Presence of Random Distribution of Microhardness of Workpiece and Nonlinear Cutting Forces in Lathe Process

Academic Editor: D. Zhou
Received02 Jun 2011
Accepted18 Jul 2011
Published18 Sep 2011

Abstract

This work presents a comprehensive approach to the control of tool's position, in the presence of machine tool structure vibration, nonlinear cutting force, and random tool vibration due to random distribution of microhardness of workpiece material. The controller is combination of Proportional and linear quadratic gaussian- (P-LQG-) type constructed from an augmented model of both tool-actuator dynamics and a nonlinear dynamic model relating tool displacement to cutting forces. The latter model is obtained using black-box system identification of experimental orthogonal cutting data in which tool displacement is the input and cutting force is the output. The controller is evaluated and its performance is demonstrated.

1. Introduction

Performance of machining processes is assessed by dimensional and geometrical accuracy as well as the surface texture of the part. Factors such as cutting force, workpiece vibration, machine-tool vibration, process instability, tool wear, and thermal deformation deteriorate this performance. Recent interest in high-speed machining and high manufacturing efficiency requires faster, higher bandwidth actuation and controllers that are more robust. Besides determining part functional behavior, the surface texture also plays a key role in the area of manufacturing process control [1]. The basis for using surface texture as means of process control is derived from the fact that a slight change in the manufacturing process manifests itself as a corresponding change in the resulting surface geometry. The roughness, waviness, lay, and flaws constitute the texture of the work piece surface [2]. Surface texture of the work piece is highly dependent on tool position during cutting which in turn is affected by cutting dynamics. Therefore, extensive research work has been done on modeling and analysis of tool work-piece interaction and the dynamics of the cutting process. Altintas et al. [3] have presented a cutting force model as a function of regenerative chip thickness, cutting speed, and the velocity and acceleration terms of the vibrations. In their work, amplitude and frequency of inner and outer vibration waves are generated on the chip by an instrumented fast tool servo-powered by a piezo actuator. They established a model in which the coefficients of dynamic cutting forces are identified and used in analyzing the effect of cutting speed, tool wear, vibration frequency, and wavelength on the chatter stability of a turning process. Ikua et al. [4] have studied the cutting force generated in ball end milling and presented a model for the radial and tangential forces as a function of depth of cut and other factors related to the geometry of the tool and workpieces.

Researchers have tried to establish quantitative relations between the surface finish and the cutting parameters of the depth of cut, feed, and cutting speed using methods such as wavelength decomposition surface roughness, wavelet analysis, and tool vibration [4ā€“9]. The cutter-workpiece interaction forces are assumed to cause relative displacement between the cutter and the work piece, which influences the surface generation mechanism. Analytical, experimental and mechanistic methods were used to predict the interaction forces. The analytical methods focused on establishing a relationship between the cutting force and the instantaneous uncut chip cross-section [10, 11] and the nonlinear mechanisms [12, 13]. The analytical models were not capable of predicting the dynamic forces accurately due to the secondary nonlinear effects that stems from the tool/workpiece interaction. The experimental methods include static, dynamic, and time-series methods. The static methods are assumed based on linear assumptions and linear cutting conditions [14]. The dynamic methods replaced the static methods, and the cutting process is assumed to be a combination of two independent actions, namely, the wave cutting action and the wave removing action [15]. The time-series method was formulated to identify the dynamic cutting force coefficients as well as the transfer functions of the three-dimensional dynamic cutting process directly from operating data [16]. The mechanistic modeling methods view the machining process as a combination of the chip load/cutting force relationship, cutting-tool geometry, machining conditions, and tool/work-piece displacement due to cutting forces. The workpiece dynamic behavior is constructed using either the distributed-parameter [17] or the lumped-parameter approaches [18]. Few researchers investigated the dynamic effect of the tool holder. Shawky and Elbestawi [19] concluded that the effect of tangential force is insignificant, and no model of the tool holder is needed in that direction. Their experimental work showed that modeling the tool by a second order system in each axis provides significant approximation of the dynamic behavior of a single-point tool during cutting [19].

Control of machining process includes three levels of control that one might encounter in a controller for machining process; namely, they are the servo control, the process control [20], and the supervisory control [21]. In the servo control process, the motion of the cutting tool is taken relative to that of the workpiece. The process control level is used to control the cutting forces and tool wear to maintain high production rates and good part quality. The highest level of control is the supervisory control and it directly measures product-related variables such as part dimensions and surface finish. Different approaches are used to correct for machining errors by means of dynamic error compensation. More recently, Tian et al. [22] have presented a methodology for modeling and control of a high-precision flexure-based mechanism for ultraprecision turning operation based on the position control of an auxiliary precision mechanism utilized on the turret of the conventional lathe to implement nanometer level in feed. El-Sinawi [23] has presented an optimal control method for controlling the tool position in orthogonal cutting in both feed are radial directions. Moradi et al. [24] have presented a robust control method for orthogonal turning process in the process was modeled as a single degree of freedom model that includes the effect of tool flank wear with a control input of the system being force applied to the tool provided by a piezoactuator. Huang et al. [25] have presented a tool wear detection based on cutting force monitoring.

In this work, a control scheme for the purpose of improving surface texture of turned surfaces through control of tool position in the presence of nonlinear cutting force is developed.

Improvement of surface quality is achieved via active positioning of the tool or cutter through the attenuation of cutting forces effect on the tool in both radial and tangential directions. The process is assumed to be stochastic due to both process and measurement noise. Experimental force-displacement data is used to construct a nonlinear dynamic model of dynamic cutting forces. Various nonlinear models are constructed using system identification techniques including ARX, nonlinear ARX, and Hammerstein-Wiener techniques [26]. The best nonlinear model obtained that closely fits experimental data is then linearized and constructed in state-space form to later utilize in the linear quadratic Gaussian controller. Simulation results will be used to verify the effectiveness of the proposed modeling and control approach and to enhance current understanding of the mechanisms responsible for generating both stochastic and deterministic components of the surface texture. Simulation work of the proposed modeling and control techniques will be based on actual parameters for actuators and the cutting process data, which is readily available in the literature [27].

2. Model Development

2.1. Nonhomogeneous Distribution of Microhardness

A statistical approach was proposed by Zhang and Kapoor [28] to formulate the nonhomogeneous distribution of microhardness which is based on random sampling processes in mixing of chemicals presented by Tucker [29]. A representative part is chosen from the work-piece material. The microstructural constituents as well as their microhardness values could be experimentally identified, say š»š‘–for š‘–=1,2,ā€¦,š‘›, where š‘› represents the number of measurement locations on the chosen part. Based on these data, the three statistical parameters of šœ‡š‘Ž, variance šœŽ2š‘Ž, and the correlation coefficient function šœŒ(š‘Ÿ), are calculated using the standard statistical formulas. Note that the parameter š‘Ÿ of the correlation coefficient function represents the distance in space between the ith and iā€²th measured locations. These three parameters provide quantitative information on the size, the shape, and the segregation of microstructures and depict the microhardness distribution present in the material being machined. They also form a basis for quantifying the effect of microconstituent scale on the bulk properties of hardness of the material being machined.

2.2. Sampling Process and Parameter Estimation

As shown in Figure 1, the geometric shape of each sample is a parallelepiped being šœ‹š·/š‘›š‘  long and a cross-sectional area equal to a product of depth of cut and feed. Due to the nonhomogeneity, each sample holds its own distinct ratio of microstructures. Therefore, šœ‡š‘  is different from sample to sample. The sample hardness means šœ‡š‘ š‘— for the jth sample can be obtained asšœ‡š‘ š‘—=1š‘š‘ š‘š‘ ī“š‘–=1š»š‘ ,(1) where š‘š‘  is the number of locations within a sample (a subset of measurement locations, š‘›) and š‘—=1,2,ā€¦,š‘›š‘ . The mean hardness value of all sample (šœ‡š‘ 1,šœ‡š‘ 2,ā€¦,šœ‡š‘ š‘›š‘ ) would be equal if the microstructures in the material were uniformly distributed (in a microscale). However, this is never the case in practice [28]. As a result, the sample hardness mean, šœ‡š‘ , differs from sample to sample. In fact, this variation reveals, from a microbase analysis, the main source of random excitation phenomenon observed during machining. When the cutting tool is cutting a sample (a portion of the circumference of the work-piece being cut) with a small value of mean hardness, the magnitude of the cutting force drops. On the other hand, the cutting force jumps up when the cutting tool is meeting a sample with a large value of mean hardness. Because the variation in the mean hardness value from sample to sample is of random nature, the dynamic variation of the cutting force is also random in nature. The larger the mean hardness value variation, the stronger the random excitation system present during machining.

Since the sample hardness mean šœ‡š‘ , the key parameter of the random excitation system, is a random variable, it is best dealt with statistically. It can be shown that the mean of the sample mean distribution is equal to the mean of its population distribution that is, šœ‡š‘ =šœ‡š‘Ž.. However, estimating the variance of the sample hardness mean distribution or the sample mean variance, using the standard statistical formulas, is not practical, because such an evaluation would require the subgrouping of the identified locations of microstructures into individual samples. Therefore, this procedure to evaluate šœŽ2š‘  would be extremely tedious if applied in practice. Furthermore, the main drawback of this procedure is that quantitative relationship between the sample mean variance and the geometric sample shape; namely, the three cutting parameters cannot be easily established. For this reason, the sample mean variance theory [29] is employed to establish this quantitative relationship. The sample mean variance theory stemmed from the need to characterize the composition uniformity of a chemical mixing process, a similar process of taking random samples from the population distribution and evaluating the sample mean values afterwards. The theory offers a mathematical approach to quantitatively evaluate the sample mean variance based on the knowledge of the mean, the variance, and the correlation coefficient function of the population distribution, from which the random samples are being taken.

To estimate the variance of the sample hardness mean distribution šœŽ2š‘ , we start by assuming the work-piece material to be a two-component mixture consisting of components A and B. This indicates that the cutting edge (tool) experience different microstructural portions during cutting. The analysis may be generalized into multiple-component mixtures by letting one component be A and all others (including free space) to be B. in this case the analysis must be performed (š‘›āˆ’1) times for a mixture of š‘› components. If samples with identical geometric shape (each has a volume š‘‰) are taken from the mixture, the average concentration of component A in an identical sample is given byšœ‡š‘ =1š‘‰ī€œš‘‰š‘Žš‘‘š‘£,(2) while the average of šœ‡š‘  will depend on the average of š‘Ž; that is,šœ‡š‘ =1š‘‰ī€œš‘‰š‘Žš‘‘š‘£.(3)

For convenience, we define the quantity š‘ as š‘=š‘Žāˆ’š‘Ž,(4) so that the deviation of šœ‡š‘  from the average for any particular sample is šœ‡āˆ’šœ‡š‘ =1š‘‰ī€œš‘‰š‘(š‘„)š‘‘š‘£.(5)

The square deviation of šœ‡š‘  is now ī€·šœ‡āˆ’šœ‡š‘ ī€ø2=1š‘‰2ī€œš‘‰ī€œš‘‰ā€²ī€·š‘„š‘(š‘„)š‘ī…žī€øš‘‘š‘£ī…žš‘‘š‘£,(6) the double integral is taken over all pairs of incremental volumes within the sample.

If the square deviation of šœ‡š‘  is averaged over all samples in the mixture, one gets the sample variance šœŽ2š‘ šœŽ2š‘ =ī€·šœ‡š‘ āˆ’šœ‡š‘ ī€ø2.(7)

Let š‘€ be the number of samples in the mixture. If we take š‘¦š‘– to be the spatial location of a reference point in the ith sample, and if š‘„ is referred to a local coordinate system with its origin at š‘¦š‘– and oriented with respect to the sample, then (6) and (7) can be combined to give šœŽ2š‘ =1š‘€š‘€ī“š‘–=1ī€œš‘‰ī€œš‘‰ā€²š‘ī€·š‘„+š‘¦š‘–ī€øš‘ī€·š‘„ī…ž+š‘¦š‘–ī€øš‘‘š‘£ī…žš‘‘š‘£.(8)

The order of the integration and summation can be reversed, bringing the summation inside the integral. Then, inspection of (5) shows that1š‘€š‘€ī“š‘–=1š‘ī€·š‘„+š‘¦š‘–ī€øš‘ī€·š‘„ī…ž+š‘¦š‘–ī€ø=šœŽ2š‘ŽšœŒ(š‘Ÿ),(9) where šœŽ2š‘Ž is the variance in concentration among all samples in the mixture, šœŒ(š‘Ÿ) is called the correlation coefficient function, and š‘Ÿ is the magnitude of the distance between the points š‘„ and š‘„ī…ž for||š‘Ÿ=š‘„āˆ’š‘„ī…ž||.(10) Equation (8) can now be expressed as the average of the correlation coefficient function over a single representative samplešœŽ2š‘ =1š‘‰2ī€œš‘‰ī€œš‘‰ā€²šœŽ2š‘ŽšœŒ(š‘Ÿ)š‘‘š‘£ī…žš‘‘š‘£.(11)

The outer integral in (11) specifies an incremental volume dv, for which the inner integral is evaluated. This inner integral can be expressed in spherical coordinates, ī€œš‘‰ā€²šœŽ2š‘ŽšœŒ(š‘Ÿ)š‘‘š‘£ī…ž=ī€œšœ™ī€œšœƒī€œš‘ŸšœŽ2š‘ŽšœŒ(š‘Ÿ)š‘Ÿ2sin(šœ™)š‘‘š‘Ÿš‘‘šœƒš‘‘šœ™,(12) with the limits on š‘Ÿ,šœƒ, and šœ™ chosen to limit the integral to points within the sample volume as shown in Figure 1. Note that š‘Ÿ is simply the distance between the volume increment š‘‘š‘£ī…ž and the increment dv specified by the outer integral. Let š‘„ be the location of this later incremental volume. To avoid the difficulty of specifying limits on šœƒ and šœ™ in (12), we define a function š‘Šāˆ—(š‘„,š‘Ÿ) as follows. Consider a spherical shell of radius š‘Ÿ and thickness dr, centered at š‘„. š‘Šāˆ—(š‘„,š‘Ÿ) is the fraction of the volume of that shell lying inside the sample as shown in Figure 3. Since the integrand in (12), we define a function of š‘Ÿ only, one can easily evaluate the integrals with respect to šœƒ and šœ™ when the integrand is multiplied by š‘Šāˆ—. That is, ī€œš‘‰ā€²šœŽ2š‘ŽšœŒ(š‘Ÿ)š‘‘š‘£ī…ž=ī€œāˆžš‘Ÿ=0šœŽ2š‘ŽšœŒ(š‘Ÿ)š‘Ÿ2sin(šœ™)4šœ‹š‘Ÿ2š‘‘š‘Ÿ.(13)

The infinite upper limit is used on š‘Ÿ, because š‘Šāˆ— will be zero when š‘Ÿ is so great that the entire shell lies outside the sample.

Substituting back into (11) and rearranging givesšœŽ2š‘ =4šœ‹šœŽ2š‘Žš‘‰ī€œāˆžš‘Ÿ=0šœŒ(š‘Ÿ)š‘Ÿ2ī‚ø1š‘‰ī€œš‘‰š‘Šāˆ—ī‚¹(š‘„,š‘Ÿ)š‘‘š‘£š‘‘š‘Ÿ.(14)

The terms within the square brackets in (14) are related to the sample geometry only. Define this part as the sample shape function š‘Š(š‘„,š‘Ÿ), because it can be integrated with respect to the sample volume separately.1š‘Š(š‘Ÿ)=š‘‰ī€œš‘‰š‘Šāˆ—(š‘„,š‘Ÿ)š‘‘š‘£,(15) by substituting (15) into (14), the formula to evaluate the sample mean variance is given byšœŽ2š‘ =4šœ‹šœŽ2š‘Žš‘‰ī€œāˆžš‘Ÿ=0šœŒ(š‘Ÿ)š‘Ÿ2š‘Š(š‘Ÿ)š‘‘š‘Ÿ.(16)

Substituting for š‘‰ from Figure 1, (16) becomesšœŽ2š‘ =4š‘›š‘ šœŽ2š‘Žī€œš·š‘“š‘‘āˆžš‘Ÿ=0šœŒ(š‘Ÿ)š‘Ÿ2š‘Š(š‘Ÿ)š‘‘š‘Ÿ.(17)

Note that (17) consists of two parts. One part is the integral, and the other is the constant term (4š‘›š‘ šœŽ2š‘Ž/š·š‘“š‘‘), where š‘›š‘ is the number of samples, š· is the workpiece diameter, š‘“ is the feed rate, and š‘‘ is the depth of cut. For a given work-piece material under different cutting parameter settings, the constant term decreases as any of the terms in the denominator increases. The integral term always increases as any of the three parameters f, d, and/or D increases, because the geometric shape factor š‘Š(š‘Ÿ) always increases (sometimes very slightly), and at the same time, the correlation coefficient function šœŒ(š‘Ÿ) remains unchanged as long as the material of the work-piece remains unchanged. Therefore, the magnitude of the sample mean variance šœŽ2š‘  is a compromise between these two terms. In fact, (17) serves as a mathematical model to predict values of the sample mean variance under different cutting parameter settings.

It is well known in statistics that the sampling distributions of sample means are normal even though the samples may be drawn from a nonnormal population such as uniform distributions. The mathematical justification of this relationship is the central limit theorem. Based on this theorem, the distribution of the sample hardness mean, šœ‡š‘ , is normal. The mean of the sample hardness means is expressed asšœ‡š‘ ī€ŗšœ‡=šøš‘ ī€»=šœ‡,(18) where šœ‡š‘Ž is the mean of the population distribution of hardness and the variance of the sample hardness means is given by (17).

Since the mean level of the dynamic cutting force is directly related to the mean value šœ‡š‘ , this indicates that the random excitation system (dynamic cutting force) can be mathematically describes by a normal distribution as shown in Figure 2.

The variance of the normal distribution šœŽ2š‘  represents the variation level o the dynamic cutting force about its mean level. The larger the variance šœŽ2š‘ , the more significant the šœ‡š‘  variation as well as the variation level of the dynamic cutting force. As the variance šœŽ2š‘  decreases, the variation level of the dynamic cutting force decreases accordingly. When the variance šœŽ2š‘  approaches zero, indicating that all sample mean hardness values are equal to each other, the random excitations caused by the cutting force diminishes.

Finally, the following very important conclusions can be drawn for the preceding sections.(1)The random excitation system can be mathematically described by a normal distribution.(2)The two parameters šœ‡š‘  and šœŽ2š‘  of the normal distribution fully describe the characteristics of the random excitation system. (3)The variation level of the dynamic cutting force represented by šœŽ2š‘ , is directly affected by the variation of the cutting parameters; see (17).

These conclusions are vital to the design of the proposed Kalman estimator, which will be the presented later on.

3. Nonlinear Force Model

Force-displacement data obtained from experimental cutting of a 6061 Aluminum workpiece using a carbide tool with 0Ā° rake angle and 7Ā°clearance angle, depth of cut of 0.4ā€‰mm, feed of 0.050ā€‰mm/rev, and spindle speed: 2200ā€‰rev/min, which are shown in Figures 5 and 6, have been used to construct a force displacement model using black-box system identification.

Four dynamic models with tool displacement as input and cutting force as output are constructed are indentified and compared as shown in Figure 3. In reference to Figure 3, the first model identified as LinMod1 is a linear ARX model of the form š‘“š‘(š‘”)+š‘Ž1š‘“š‘(š‘”āˆ’1)+š‘Ž2š‘“š‘(š‘”āˆ’2)+ā‹Æ+š‘Žšœš‘“š‘(š‘”āˆ’šœ)=š‘1š‘„š‘(š‘”)+š‘2š‘„š‘(š‘”āˆ’1)+ā‹Æ+š‘šœ‡š‘„š‘(š‘”āˆ’šœ‡āˆ’1)+š‘’(š‘”),(19) where š‘“š‘,š‘„š‘ are the cutting force and tool displacement in the corresponding direction, respectively. š‘“š‘(š‘”āˆ’1)+š‘“š‘(š‘”āˆ’2),ā€¦,š‘“š‘(š‘”āˆ’šœ),š‘„š‘(š‘”)+š‘„š‘(š‘”āˆ’1),ā€¦,š‘„š‘(š‘”āˆ’šœ‡āˆ’1) are delayed input and output variables called regressors. Linear ARX model predicts the output š‘“š‘ as a weighted sum of its regressors. šœ is the number of past output terms, while šœ‡ is the number of past input terms used to predict the current output š‘“š‘, and š‘’(š‘”) is a white noise sequence. The second and third models identified as LinMod3 and LinMod4 are constructed as follows.

LinMod3 is constructed using predictive error minimization or PEM with various constraints on the order or number of states of the state-space model. This method yields a discrete time-domain state-space dynamic model of the formš‘¦š‘„(š‘”+š‘‡š‘ )=š“š‘„(š‘”)+šµš‘¢(š‘”)+š¾š‘’(š‘”),(š‘”)=š¶š‘„(š‘”)+š·š‘¢(š‘”)+š‘’(š‘”),(20) in which š‘„ are the states, š‘¦ is output, š“, šµ, š¶, š¾, š· represent dynamics, input output, disturbance and feed through matrices, respectively [30]. Narx1 model is constructed using Nonlinear ARX model which is an extension of ARX model given in (1) except that the input-output mapping is nonlinear and of the form;š‘“š‘›ī€·š‘“(š‘”)=š¹š‘(š‘”āˆ’1),š‘“š‘(š‘”āˆ’2),ā€¦,š‘„š‘(š‘”),š‘„š‘ī€ø(š‘”āˆ’1),ā€¦,(21) where š‘“š‘› is the output and š¹ is a nonlinear mapping function, such as wavelet network, tree partition, and so forth, [31]. Nhw1 model shown in Figure 3 is constructed using Hammerstein-Wiener model; see [31] for further details. It is clear from Figure 3 that the Narx1 has the best fit to experimental data with 99.9% fit.

3.1. Linearization of the Nonlinear ARX Model

The nonlinear model that best fits input-output data is clearly the nonlinear ARX model mapped with a single-layer sigmoid function. The nonlinear mapping of input-output data done with š¹(š‘£)=šæ(š‘£āˆ’š‘Ÿ)+š‘‘+š‘”(š‘„(š‘£āˆ’š‘Ÿ)),(22) where š‘£ is a vector of the regressors. šæ(š‘£āˆ’š‘Ÿ)+š‘‘ is the output of the linear block, š‘‘ is a scalar offset, š‘”(š‘„(š‘£āˆ’š‘Ÿ)) is the output of the nonlinear function block, and š‘„ is a projection matrix and š‘Ÿ is the mean of regressors (š‘£). Function š‘”(š‘£)expressed as āˆ‘š‘”(š‘£)=š‘›š‘˜=1š›¼š‘˜šœ™(š›½š‘˜(š‘£āˆ’šœ‰š‘˜)) with šœ™(š‘ )=(š‘’š‘ +1)āˆ’1 is the nonlinearity estimator of the sigmoid function type. š›½š‘˜ is a row vector such that š›½š‘˜(š‘£āˆ’šœ‰š‘˜) is a scalar, š‘  is the Laplace variable, and š‘› is the number of sigmoid network units. The model obtained by this approach is of the form given in (3) with šœ=4, and šœ‡=4.

Steady state value of the step response can be used later on, in linearizing the model. See the appendix for linearization details.

Linearized state-space model of the radial force is Ģ‡š‘‹š‘Ÿ(š‘”)=š“š‘Ÿš‘‹š‘Ÿ(š‘”)+šµš‘Ÿš‘„š‘š‘Ÿ,š‘“š‘„=š¶š‘Ÿš‘‹š‘Ÿ(š‘”)+š·š‘Ÿš‘„š‘š‘Ÿ,(23) where Ģ‡š‘‹š‘Ÿ(š‘”) are the states of the radial force model, š“š‘Ÿ,šµš‘Ÿ,š¶š‘Ÿ, and š·š‘Ÿ are dynamic, input, output, and feed through matrices, respectively. š‘“š‘„,š‘„š‘š‘Ÿ are radial cutting force and radial tool displacement, respectively. The same procedure is carried out for the tangential cutting force model which yieldsĢ‡š‘‹š‘”(š‘”)=š“š‘”š‘‹š‘”(š‘”)+šµš‘”š‘„š‘š‘”,š‘“š‘¦=š¶š‘”š‘‹š‘”(š‘”)+š·š‘”š‘„š‘š‘”.(24)

3.2. Modeling of the Machining Process

Figure 4 shows the difference between tool post (passive) and proposed active tool fixture that will allow for implementation of the proposed controller. Figure 4(a) shows the conventional tool post rigidly attached to the machine tool structure, while Figure 4(b) shows the active tool holder platform equipped with two actuators placed between the tool and the machine tool structure. This will allow the control force to provide necessary manipulation of the tool to maintain a constant depth of cut.

During cutting, the tool is perturbed form nominal depth of cut by two inputs; namely, the dynamic cutting forces in both radial and tangential directions. š‘¤(š‘”) represents nonhomogenous micro-hardness distribution model presented in Section 2, producing the stochastic component of the cutting force.

State-space model of the tool-actuator assembly in the radial direction can be represented asĢ‡š‘„š‘Ÿ=š“š‘š‘Ÿš‘„š‘Ÿ+šµš‘š‘Ÿš¹š‘Žš‘Ÿ,š‘Œš‘š‘Ÿ=š¶š‘š‘Ÿš‘„š‘Ÿ+š·š‘š‘Ÿš¹š‘Žš‘Ÿ.(25) such that š“š‘š‘Ÿ=ī€ŗ01āˆ’š‘˜š‘Žš‘/š‘šāˆ’š‘š‘Žš‘/š‘šī€»,šµš‘š‘Ÿ=ī€ŗ01/š‘šī€»,š¶š‘š‘Ÿ=[10],š·š‘š‘Ÿ=0, where š‘„š‘Ÿ is vector of the tool states, namely, displacement and velocity. š‘š is the mass of the tool and its mount, š‘˜š‘Žš‘,š‘š‘Žš‘ are the elastic stiffness and damping coefficient of the actuator, respectively. š‘˜š‘Žš‘,š‘š‘Žš‘ can be obtained by an FR test of the tool-actuator assembly. Notice that the same procedure can be carried out for the tangential direction, since both actuators are orthogonal. Refer to [23, 27] for more information on the model development.

4. Control Strategy and the Controller Design

The control strategy will be centered on reducing the amplitude of the tool' dynamic displacement to zero for the purpose of maintaining a constant depth of cut and subsequently, a smooth surface texture of the workpiece. This requires minimization of the error between the desired and the actual position of the tool. The desired or nominal tool position is the one yielding a constant depth of cut. However, the existence of the dynamic cutting forces perturbs the tool form its nominal position, and thus, varies the depth of cut causing subsequent deterioration in surface texture. Attenuating the effect of the cutting force on the toolā€™s position, the error in the tool position will be minimized. This can be achieved by reacting on the tool with an equal but opposite forces through proper actuation. This task is difficult due to (a) the existence of process and measurement noise and (b) the nonlinearity of the dynamic cutting forces. To overcome these difficulties, the controller must be able to estimate the force needed to minimize the error in the tool position from noisy process using measurement data contaminated with noise. In addition to that, the controller has to track and manipulate the position of the tool effectively in the presence of nonlinear dynamic cutting forces. Added to all of the above, the controller has to maintain high stability and performance under various disturbance characteristics. To do so, the controller has to be constructed based on a model that takes into account the dynamics of the tool-actuator system as well as the cutting force dynamics. Therefore, The plant used for this purpose is a combination (i.e., augmentation) of the dynamic cutting forces and tool-actuator dynamics in both radial and tangential directions. The first step in designing the proposed active controller is the design of the LQG estimator and regulator gains šæš‘’ and š¾š‘…, respectively. The analysis presented here is in state-space such that the augmented system depicted by (23) and (25) is;š“š‘Ž=āŽ”āŽ¢āŽ¢āŽ£š“š‘š‘Ÿ0šµš‘Ÿš¶š‘š‘Ÿš“š‘ŸāŽ¤āŽ„āŽ„āŽ¦,šµš‘Ž=āŽ”āŽ¢āŽ¢āŽ£šµš‘š‘Ÿšµš‘Ÿš·š‘š‘ŸāŽ¤āŽ„āŽ„āŽ¦,š¶š‘Ž=ī€ŗš·š‘Ÿš¶š‘š‘Ÿš¶š‘Ÿī€»,š·š‘Ž=ī€ŗš·š‘Ÿš·š‘š‘Ÿī€»,(26) where the matrices š“š‘Ž, šµš‘Ž, š¶š‘Ž, and š·š‘Ž are, respectively, the dynamic, input, output, and direct transmission matrices of the augmented system. The four matrices in the foregoing are used to design the LQG estimator and regulator gain matrices šæš‘’ and š¾š‘… such that šæš‘’=āŽ”āŽ¢āŽ¢āŽ£ī€ŗšæš‘’š‘š‘Ÿī€»ī€ŗšæš‘’š‘Ÿī€»āŽ¤āŽ„āŽ„āŽ¦,š¾š‘…=āŽ”āŽ¢āŽ¢āŽ£ī€ŗš¾š‘…š‘š‘Ÿī€»ī€ŗš¾š‘…š‘Ÿī€»āŽ¤āŽ„āŽ„āŽ¦.(27)

This augmentation is necessary to incorporate the cutting force dynamics in the process of determining the optimal estimator and regulator gains. Equation (27) shows that šæš‘’ and š¾š‘… are both partitioned in two parts each, where []š‘š‘Ÿ corrects š‘„š‘š‘Ÿ (i.e., tool-actuator model states) and []š‘Ÿ corrects š‘‹š‘Ÿ (i.e., cutting force states). The LQG (virtual model of the system) should be subjected to all inputs that the actual plant is subjected to, including the control force. Figure 7 shows the control scheme implementation in the radial direction. The proportional controller is used to adjust for any variation in the cutting force due to variation of any or all of the cutting parameters.

5. Numerical Example

Force-displacement results obtained from experimental cutting of a 6061 Aluminum workpiece using a carbide tool with 0Ā° rake angle and 7Ā° clearance angle, depth of cut of 0.4ā€‰mm, feed of 0.050ā€‰mm/rev, and spindle speed: 2200ā€‰rev/min, shown in Figures 5 and 6, are used to construct a force-displacement model using black-box system identification. Linearized State-space model obtained by nonlinear ARX modeling presented here in canonical form isĢ‡š‘‹š‘ŸāŽ”āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ£āŽ¤āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ¦š‘‹(š‘”)=1.145āˆ’0.70.5651āˆ’0.443520000100000.50š‘ŸāŽ”āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ£000āŽ¤āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ¦š‘„(š‘”)+1024š‘š‘Ÿ,š‘“š‘„[]š‘‹=1500Ɨāˆ’266.2153.9āˆ’213.8460.7š‘Ÿ[0]š‘„(š‘”)+š‘š‘Ÿ.(28) ETREMA actuator with peak-to-peak excursion of 50Ɨ10āˆ’6š‘š is utilized as the active manipulator of the tool. The actuatorā€™s elastic stiffness and damping used in this study were taken as follows [27], š‘š=0.53kg for the moving tool-actuator assembly. š‘˜š‘Žš‘=14.6Ɨ106ā€‰N/m, and š‘š‘Žš‘=10ā€‰kg/s. Figure 8 shows the base force transmitted from the machine-tool structure to the cutter through the actuator. This force is due to the vibration of the machine tool at low resonance frequencies as indicated in [27]. Figure 9 shows the vibration of the cutter due to random distribution of microhardness according to the model developed in Section 2 of this paper and represented by š‘¤(š‘”) in Figure 7. The measurement noise š‘£(š‘”) is assumed to be random with zero mean and variation of Ā±10% of measurement values.

6. Results and Analysis

State-space model of the augmented system given in (26) is used to determine the LQG gains. LQG gains are highly dependent on weight matrices of the performance index. See [32] for further details on the selection of performance index weights.

The controller is implemented in both radial and tangential directions as outlined in Figure 7. However, results here are only for the radial direction (depth of cut direction), since the implementation in the feed direction is only a duplicate of the radial. Figure 10 shows the toolā€™s displacement in the radial direction with and without control, where it is clear from the latter that the controller has succeeded in reducing the dynamic tool-displacement by approximately 40%ā€“50%. The controller maintained excellent stability in the presence of process and measurement noise as shown by the bounded response of the cutter after control implementation. Control force in the radial direction is shown in Figure 11, which is well within the actuatorā€™s force capabilities. Cutting force as modeled by the nonlinear system identification shows that the cutting force has a significant magnitude especially when the random force due to nonhomogenous distribution of micro-hardness is taken into consideration. Figure 12 shows the cutting force to be reduced significantly when the tool displacement (i.e., vibration) is reduced as shown in Figure 10. The proportional gain used in simulation is š‘ƒ=1Ɨ106 and can be tuned for better performance of the control scheme with upper limit determined based on the saturation force of the actuator.

7. Conclusions

In this research work, active control of an orthogonal cutting process is presented. Nonlinear dynamic cutting force model is generated from actual machining data obtained in both radial and tangential directions. An LQG-based controller is developed based on an augmented model of both tool-actuator dynamics and a linearized model of the dynamic cutting force. The control objective is to eliminate the toolā€™s dynamic displacement and maintain a constant depth of cut. Simulation results have shown that the proposed control strategy has managed to significantly reduce the dynamic displacement of the tool with minimal force and calculations efforts even in the presence of significant randomness in the process and measurement as well as the nonlinearity of the cutting force. Experimental study is undergoing to verify the integrity and performance of the controller in real-time applications.

Appendix

Linearization of a Nonlinear System (Ljung [31])

Common use of time-varying linear systems is related to linearization of nonlinear systems around a certain ā€œnominalā€ trajectory. Therefore, for a nonlinear system described byš‘¦š‘„(š‘”+1)=š‘“((š‘„(š‘”),š‘¢(š‘”))+š‘Ÿ(š‘„(š‘”),š‘¢(š‘”))š‘¤(š‘”),(š‘”)=ā„Ž(š‘„(š‘”))+š‘š(š‘„(š‘”),š‘¢(š‘”))š‘£(š‘”).(A.1)

Assuming that the disturbance terms š‘¤(š‘”)andš‘£(š‘”) are white and small, with the nominal and disturbance-free (š‘¤(š‘”)ā‰”0,š‘£(š‘”)ā‰”0) behavior of the system corresponds to an input sequence š‘¢āˆ—(š‘”) and corresponding trajectoryš‘„āˆ—(š‘”). Neglecting nonlinear terms, the differences Ī”š‘„(š‘”)=š‘„(š‘”)āˆ’š‘„āˆ—ī€·š‘„(š‘”),Ī”š‘¦(š‘”)=š‘¦(š‘”)āˆ’ā„Žāˆ—ī€ø,(š‘”)Ī”š‘¢(š‘”)=š‘¢(š‘”)āˆ’š‘¢āˆ—(š‘”).(A.2)

are then subjected to Ī”š‘„(š‘”+1)=š¹(š‘”)Ī”š‘„(š‘”)+šŗ(š‘”)Ī”š‘¢(š‘”)+š‘¤(š‘”),Ī”š‘¦(š‘”)=š»(š‘”)Ī”š‘„(š‘”)+š‘£(š‘”).(A.3) such that, šœ•š¹(š‘”)=|||šœ•š‘„š‘“(š‘„,š‘¢)š‘„āˆ—(š‘”),š‘¢āˆ—(š‘”),šœ•šŗ(š‘”)=|||šœ•š‘¢š‘“(š‘„,š‘¢)š‘„āˆ—(š‘”),š‘¢āˆ—(š‘”),šœ•š»(š‘”)=|||šœ•š‘„ā„Ž(š‘„)š‘„āˆ—(š‘”).(A.4)

In view of the assumption of small disturbances, the cross terms with disturbance. The terms š‘¤(š‘”)and š‘£(š‘”)in (A.3) above are white disturbances with the following covariance properties:š‘…1(š‘”)=šøš‘¤(š‘”)š‘¤š‘‡(ī€·š‘„š‘”)=š‘Ÿāˆ—(š‘”),š‘¢āˆ—ī€ø(š‘”)šøš‘¤(š‘”)š‘¤š‘‡(š‘”)š‘Ÿš‘‡ī€·š‘„āˆ—(š‘”),š‘¢āˆ—ī€ø,š‘…(š‘”)2(š‘”)=šøš‘£(š‘”)š‘£š‘‡ī€·š‘„(š‘”)=š‘šāˆ—(š‘”),š‘¢āˆ—(ī€øš‘”)šøš‘£(š‘”)š‘£š‘‡(š‘”)š‘šš‘‡ī€·š‘„āˆ—(š‘”),š‘¢āˆ—(ī€ø,š‘…š‘”)12ī€·š‘„(š‘”)=š‘Ÿāˆ—(š‘”),š‘¢āˆ—ī€ø(š‘”)šøš‘¤(š‘”)š‘£š‘‡(š‘”)š‘šš‘‡ī€·š‘„āˆ—(š‘”),š‘¢āˆ—ī€ø.(š‘”)(A.5) The model is now linear in the vicinity of nominal trajectory.

Acknowledgment

The author acknowledges of the support of the American University of Sharjah.

References

  1. M. Shiraishi and H. Sumiya, ā€œImprovement of geometrical errors by surface roughness and tool position controls,ā€ in Proceedings of the Annual Meeting of the American Society of Mechanical Engineers, PED-vol. 45, pp. 9ā€“22, 1990. View at: Google Scholar
  2. M. Sander, A Practical Guide to the Assessment of Surface Texture, Mahr Perthen, Gottingen, Germany, 1991.
  3. Y. Altintas, M. Eynian, and H. Onozuka, ā€œIdentification of dynamic cutting force coefficients and chatter stability with process damping,ā€ CIRP Annalsā€”Manufacturing Technology, vol. 57, no. 1, pp. 371ā€“374, 2008. View at: Publisher Site | Google Scholar
  4. B. W. Ikua, H. Tanaka, F. Obata, and S. Sakamoto, ā€œPrediction of cutting forces and machining error in ball end milling of curved surfaces -I theoretical analysis,ā€ Precision Engineering, vol. 25, no. 4, pp. 266ā€“273, 2001. View at: Publisher Site | Google Scholar
  5. S. M. Pandit and M. S. Shunmugan, ā€œSignature of machine tool errors on surface texture,ā€ Annals of the CIRP, PED-vol. 45, pp. 63ā€“74, 1990. View at: Google Scholar
  6. W. S. Lin, B. Y. Lee, and C. L. Wu, ā€œModeling the surface roughness and cutting force for turning,ā€ Journal of Materials Processing Technology, vol. 108, no. 3, pp. 286ā€“293, 2001. View at: Publisher Site | Google Scholar
  7. O. B. Abouelatta and J. MĆ”dl, ā€œSurface roughness prediction based on cutting parameters and tool vibrations in turning operations,ā€ Journal of Materials Processing Technology, vol. 118, pp. 269ā€“277, 2001. View at: Publisher Site | Google Scholar
  8. S. Fu, B. Muralikrishnan, and J. Raja, ā€œEngineering surface analysis with different wavelet bases,ā€ Journal of Manufacturing Science and Engineering, Transactions of the ASME, vol. 125, no. 4, pp. 844ā€“852, 2003. View at: Publisher Site | Google Scholar
  9. K. A. Risbood, U. S. Dixit, and A. D. Sahasrabudhe, ā€œPrediction of surface roughness and dimensional deviation by measuring cutting forces and vibrations in turning process,ā€ Journal of Materials Processing Technology, vol. 132, pp. 203ā€“214, 2003. View at: Publisher Site | Google Scholar
  10. G. B. Boothroyd and W. A. Knight, Fundamentals of Machining and Machine Tools, Marcel Dekker, New York, NY, USA, 1989.
  11. K. F. Ehmann, S. G. Kapoor, R. E. DeVor, and I. Lazoglu, ā€œMachining process modeling: a review,ā€ Journal of Manufacturing Science and Engineering, Transactions of the ASME, vol. 119, no. 4, pp. 655ā€“663, 1997. View at: Google Scholar
  12. P. Albrecht, ā€œDynamics of the metal cutting process,ā€ Journal of Engineering for Industry, vol. 87, pp. 429ā€“441, 1965. View at: Google Scholar
  13. N. H. Hanna and S. A. Tobias, ā€œA theory of nonlinear regenerative chatter,ā€ Journal of Engineering for Industry, vol. 96, pp. 247ā€“253, 1974. View at: Google Scholar
  14. M. A. El Baradie, ā€œStatistical analysis of the dynamic cutting coefficients and machine tool stability,ā€ Journal of Engineering for Industry, vol. 115, no. 2, pp. 205ā€“214, 1993. View at: Google Scholar
  15. J. Peters and P. Vanherck, ā€œMachine tool stability tests and the incremental stiffness,ā€ Annals of the CIRP, vol. 17, pp. 225ā€“232, 1969. View at: Google Scholar
  16. S. M. Pandit, T. L. Subramanian, and S. M. Wu, ā€œModeling machine tool chatter by time series,ā€ Journal of Engineering for Industry, vol. 97, no. 1, pp. 211ā€“215, 1975. View at: Google Scholar
  17. M. U. Jen and E. B. Magrab, ā€œThe dynamic interaction of the cutting process, workpiece, and lathe's structure in facing,ā€ Journal of Manufacturing Science and Engineering, Transactions of the ASME, vol. 118, no. 3, pp. 348ā€“357, 1996. View at: Google Scholar
  18. L. Kops, M. Gould, and M. Mizrach, ā€œImproved analysis of the workpiece accuracy in turning based on the emerging diameter,ā€ Journal of Engineering for Industry, vol. 115, no. 3, pp. 253ā€“257, 1993. View at: Google Scholar
  19. A. M. Shawky and M. A. Elbestawi, ā€œAn enhanced dynamic model in turning including the effect of ploughing forces,ā€ Journal of Manufacturing Science and Engineering, Transactions of the ASME, vol. 119, no. 1, pp. 10ā€“20, 1997. View at: Google Scholar
  20. L. K. Daneshmend and H. A. Pak, ā€œModel reference adaptive control of feed force in turning,ā€ Journal of Dynamic Systems, Measurement and Control, vol. 108, no. 3, pp. 215ā€“222, 1986. View at: Google Scholar
  21. T. E. Bailey, D. M. Jenkins, A. D. Spence, and M. A. Elbestawi, ā€œIntegrated modeling for metal removal operations,ā€ in Proceedings of the ASME Dynamic Systems and Control Division (DSCD '96), vol. 58, pp. 191ā€“198, 1996. View at: Google Scholar
  22. Y. Tian, B. Shirinzadeh, and D. Zhang, ā€œA flexure-based mechanism and control methodology for ultra-precision turning operation,ā€ Precision Engineering, vol. 33, no. 2, pp. 160ā€“166, 2009. View at: Publisher Site | Google Scholar
  23. A. H. El-Sinawi, ā€œTwo-dimensional vibration suppression in turning using optimal control of the cutting tool,ā€ International Journal of Machining and Machinability of Materials, vol. 3, no. 1-2, pp. 91ā€“103, 2008. View at: Publisher Site | Google Scholar
  24. H. Moradi, M. R. Movahhedy, and G. R. Vossoughi, ā€œRobust control strategy for suppression of regenerative chatter in turning,ā€ Journal of Manufacturing Processes, vol. 11, no. 2, pp. 55ā€“65, 2009. View at: Publisher Site | Google Scholar
  25. S. N. Huang, K. K. Tan, Y. S. Wong, C. W. de Silva, H. L. Goh, and W. W. Tan, ā€œTool wear detection and fault diagnosis based on cutting force monitoring,ā€ International Journal of Machine Tools and Manufacture, vol. 47, no. 3-4, pp. 444ā€“451, 2007. View at: Publisher Site | Google Scholar
  26. Y. Zhu, ā€œEstimation of an N-L-N Hammerstein-Wiener model,ā€ Automatica, vol. 38, no. 9, pp. 1607ā€“1614, 2002. View at: Publisher Site | Google Scholar
  27. A. H. El-Sinawi and A. R. Kashani, ā€œImproving surface roughness in turning using a Kalman estimator-based feed forward control of tool's position,ā€ Journal of Materials Processing Technology, vol. 167, pp. 54ā€“61, 2005. View at: Google Scholar
  28. G. M. Zhang and S. G. Kapoor, ā€œDynamic generation of machined surface, part 1: description of random excitation system,ā€ Journal of Engineering for Industry, vol. 113, no. 2, pp. 137ā€“144, 1991. View at: Google Scholar
  29. C. L. Tucker, ā€œSample variance measurement of mixing,ā€ Chemical Engineering Science, vol. 36, no. 11, pp. 1829ā€“1839, 1981. View at: Google Scholar
  30. M. Farina and L. Piroddi, ā€œAn iterative algorithm for simulation error based identification of polynomial input-output models using multi-step prediction,ā€ International Journal of Control, vol. 83, no. 7, pp. 1442ā€“1456, 2010. View at: Publisher Site | Google Scholar
  31. L. Ljung, System Identification: Theory for the User, Prentice Hall, New York, NY, USA, 2nd edition, 1999.
  32. W. Gawronsky, Advanced Structural Dynamics and Active Control of Structures, Springer, New York, NY, USA, 2004.

Copyright Ā© 2011 A. H. El-Sinawi. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

795Ā Views | 705Ā Downloads | 1Ā Citation
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.