The main objective of the paper is to design a model reference adaptive controller (MRAC) with improved transient performance. A modification to the standard direct MRAC called fuzzy modified MRAC (FMRAC) is used in the paper. The FMRAC uses a proportional control based Mamdani-type fuzzy logic controller (MFLC) to improve the transient performance of a direct MRAC. The paper proposes the application of real-coded genetic algorithm (RGA) to tune the membership function parameters of the proposed FMRAC offline so that the transient performance of the FMRAC is improved further. In this study, a GA based modified MRAC (GAMMRAC), an FMRAC, and a GA based FMRAC (GAFMRAC) are designed for a coupled tank setup in a hybrid tank process and their transient performances are compared. The results show that the proposed GAFMRAC gives a better transient performance than the GAMMRAC or the FMRAC. It is concluded that the proposed controller can be used to obtain very good transient performance for the control of nonlinear processes.
1. Introduction
Industrial applications of liquid level control are found in food processing, beverage, dairy, filtration, effluent treatment, nuclear power generation plants, pharmaceutical industries, water purification systems, industrial chemical processing, boilers, and automatic liquid dispensing and replenishment devices [1]. Traditionally, proportional-integral-derivative (PID) controllers are used for this purpose. Ziegler-Nichols and Cohen-Coon methods are widely used to tune a PID controller [2]. In these methods, the PID controller parameters are tuned for a particular operating point around which the process can be considered linear [3]. The controller would give suboptimal performance whenever the operating point is shifted out of the linearized region or when the plant parameters change due to environment and ageing. So, the PID controller parameters need to be retuned. In order to have a good performance despite the environment-induced changes in the process parameters, the controller must adapt to changes in the plant dynamics. As the level process is nonlinear, [4, 5] adaptive controllers are more suitable than fixed-parameter controllers.
Gain scheduling, dual control, self-tuning regulation (STR), and MRAC are the strategies used in adaptive control [5]. In gain scheduling, there is no feedback to make necessary changes in the controller parameter values in the event of the schedule being incorrect [5]. Dual control is too complicated to be used for real-life problems [5]. Nonlinear nature of a system warrants the use of adaptive controllers such as MRAC and STR [3]. In STR, even a small error in process parameter estimation can lead to large changes in parameters resulting in possible oscillation of process variables. Model reference adaptive controller (MRAC) has been used in the control of nonlinear processes [4, 6–8]. Goodwin and Mayne [9] established global convergence and asymptotic properties of a direct adaptive controller for continuous time stochastic linear systems. In an MRAC, the process output y takes some finite time to converge with reference model output ym.
Improving the transient performance of the MRAC has been a point of research for a long time. A modified traditional MRAC proposed by [10] has shown improved transient and steady-state performance. Miller and Davison [11] proposed a controller, which comprises an LTI compensator together with a switching mechanism, which can give an arbitrarily good transient and steady-state performance for a single-input single-output linear time-invariant plant. Asan Mohideen et al. [4] proposed a GAMMRAC that uses a PID controller along with an MRAC to improve the transient performance of the MRAC in controlling the level of a hybrid tank process. The authors have proposed an FMRAC, which uses a Mamdani-type FLC to improve the performance of MRAC [12]. The fuzzy if-then rules and the membership functions are designed based on a proportional controller. The shape and the parameters of the fuzzy membership functions are designed by using the knowledge of the operator of the plant. In order to improve the performance of the FMRAC further, the membership function parameters can be fine-tuned by using intelligent optimization algorithms like particle swarm optimization (PSO), ant colony, bacterial foraging, and GA. Chang [13] applied RGA for system identification and offline tuning of a PID controller for a system whose structure is assumed to be known previously. Hu et al. [14] used GA to optimize a fuzzy PID controller. Valarmathi et al. [15] used RGA for system identification and controller tuning for a pH process. Asan Mohideen et al. [4] used RGA for system identification and controller tuning for level control. RGA overcomes some pitfalls of the binary-coded GA such as Hamming cliff problems [16], redundancy of codes [17], time-consuming encoding and decoding, large memory requirement, and low precision due to quantization error [4, 15].
This paper proposes RGA to fine-tune the membership functions of the FMRAC offline so that the transient and steady-state performance is further improved. The resultant proposed controller is referred to as GAFMRAC in the paper. A GAMMRAC [4], an FMRAC [12], and the proposed GAFMRAC are designed and implemented to control liquid level in a coupled tank process in this study.
2. Fuzzy Logic Controller
Fuzzy logic provides a simple method to design nonlinear controllers based on heuristic approach [18]. It has been used in control of nonlinear processes such as surge tank [18], continuous stirred tank reactor [19], level process in steam generators [20], and force control in end-milling [21].
The membership functions that are most commonly used in practice are the triangular and the trapezoidal membership functions.
The mathematical representation of a triangular membership function denoted by [abc] is
(1)Triangle(x:a,b,c)={0x<a(x-a)(b-a)a≤x<b(c-x)(c-b)b≤x≤c0x>c.
The mathematical representation of a trapezoidal membership function denoted by [abcd] is
(2)Trapezoid(x:a,b,c,d)={0x<a(x-a)(b-a)a≤x<b1b≤x<c(d-x)(d-c)c≤x≤d0x>d,
where the variables a, b, c, and d are the parameters of the membership functions. The Mamdani-type FLC (MFLC) is the most common in practice [22] and is used in this paper. There are a few defuzzification methods such as max membership principle, centroid method, weighted average method, and min max membership method. The centroid method is the most physically appealing of all the defuzzification methods [22]. This paper uses the centroid method of defuzzification. The output of the defuzzifier is given as the controller output to the process.
3. Fuzzy Modified Model Reference Adaptive Controller
The FMRAC is shown in Figure 1 [12]. It employs a proportional control based MFLC to improve the transient performance of a direct MRAC. The MFLC has the tracking error e of the MRAC as its input, and the defuzzified output of the MFLC is added to the MRAC output to generate the process input u. The output of the MFLC represents volumetric flow rate in liters per hour (lph). The tracking error e is given as
(3)e=y-ym,
where y is the output of the process and ym is that of the reference model. The controller output of the FMRAC is
(4)u=θ1uc-θ2y+output(MFIS).
If the cost function J(θ) that is to be minimized is taken as e2/2, change in the value of controller parameters with respect to time as per the MIT rule [5] is given in
(5)dθdt=-γ∂J∂θ=-γ∂e∂θ,
where ∂e/∂θ is the sensitivity derivative of the system, γ is the adaptation gain, and e is the tracking error. The adaptation gains that are used in updating the values of θ1 and θ2 are γ1 and γ2, respectively. Very small values are assigned for γ1 and γ2 initially and then they are gradually increased so that y asymptotically converges with ym in the shortest possible time [4]. The cost function, J(θ)=e2/2, is used in this paper.
Fuzzy modified model reference adaptive controller.
For systems with positive steady-state gain, when e is positive, the MFLC output will reduce the controller output. This will reduce y and hence e. Similarly, when e is negative, the MFLC output will increase the controller output. This contribution from the MFLC increases the response speed of the MRAC. In order to improve the performance of the FMRAC further, the membership function parameters can be fine-tuned by using RGA. The genetic algorithm [23] is a global search and optimization method based on Darwinian theory of evolution. GA is quite successful in locating the regions of the search space containing the global optimum, but not the global optimum itself [24]. To avoid getting trapped in local minima, the procedure explained in [4] can be followed. By following this procedure, optimal or near-optimal solutions can be ensured.
4. Tuning of FMRAC Using RGA
The setup in which RGA is used to find the optimal values of fuzzy membership function parameters of the MFLC that is part of the FMRAC is shown in Figure 2. The command signal is applied to the reference model and the controller. The error and the plant output are given to GA, and it returns the optimal values of membership function parameters. The parameters are updated in each generation until the stopping criterion is met.
Controller tuning using RGA.
The objective function to be minimized is the mean square error (MSE) given as
(6)MSE=(1N)∑k=0N(y(k)-ym(k))2,
where y(k) is the output of the model at the instant k and ym(k) is the output of the reference model at the instant k.
The GA maximizes the fitness function given by
(7)F=1(1+MSE).
Maximizing “F” in effect minimizes the objective function MSE. In the denominator, a “1” is added to make sure that the denominator value never becomes zero. Tournament selection is used in this application.
5. Results and Discussion
In this section, a GAMMRAC [4], an FMRAC, and a GAFMRAC are designed for the coupled tank setup used by [4] and their performances are compared.
5.1. Process Model
Figure 3 shows the coupled tank setup. It is identified as an overdamped second-order system with zero delay. The transfer function of the identified model is
(8)G(s)=0.181630000s2+700s+1.
Coupled tank setup.
The time constant in the transfer function is scaled by a factor of 50 in [4]. The flow rate F1 and the liquid level L2 are the manipulated variable and the controlled variable, respectively. The operating region of the coupled tank process is from 18.5 cm to 41 cm with the corresponding influent flow rate ranging from 615 liters per hour (lph) to 740 lph. The maximum value of F1 is 880 lph. This control element constraint is incorporated in the simulation experiments. In this paper, zero level represents 30 cm and zero influent flow represents 677 lph and, from this point onwards, all the values of F1 and L2 are given in relative values.
5.2. Controller Tuning
In this section, a GAMMRAC, an FMRAC, and a GAFMRAC are designed for the coupled tank system. The interaction between the two tanks makes the coupled tank process more sluggish than the coupled tank system without interaction. The time constant of a single tank level process is found to be approximately 200 seconds during the initial study of the hybrid tank process. Hence, a model representing a coupled tank system without interaction is selected as the reference model initially. Some other reference models are also tried to show that the proposed controllers work well even if the reference model and the process model are very different. The transfer function of the first reference model is
(9)GM1(s)=0.140000s2+400s+1.
The adaptation rules for the MRAC parameters θ1 and θ2 are
(10)dθ1dt=-γ1e(400s+140000s2+400s+1)uc,(11)dθ2dt=-γ2e(400s+140000s2+400s+1)y,
where γ1 and γ2 are the adaptation gains for θ1 and θ2, respectively.
5.2.1. GA Based Modified MRAC
In the GAMMRAC, the controller output u [4] is
(12)u=θ1uc-θ2y-(Kpe+Ki∫edt+Kddedt),
where Kp is the proportional gain, Ki is the integral gain, and Kd is the differential gain. A GAMMRAC is designed as per (8), (9), (10), (11), and (12). RGA is used to optimize the above said parameters. A six-digit floating-point number represents each of them. The steps referred to in Section 3 are followed to ensure the optimality of the solution. Different combinations of RGA parameters are tried and random restart is carried out twenty times. The average of the individual runs represents the solution. Maximum fitness is obtained with the assignment of the values 100, 20, 0.8, and 0.02 to number of generations, population size, crossover rate, and mutation rate, respectively. The stopping criterion used here is the number of generations.
Figure 4 shows the pattern of convergence of the best and average fitness values. At the nineteenth generation, both fitness values converge for the first time. After 31 generations, the best fitness changes very little for many more generations. Hence, the estimate can be said to be optimal or near optimal. The optimal values of Kp, Ki, and Kd are 999, 15, and 124, respectively.
Convergence of fitness in GAMMRAC.
5.2.2. Fuzzy Modified MRAC
An FMRAC is designed as shown in Figure 1. The inference used in this work is “max min” method, where “min” operation is used for “and” conjunction and “max” is used for “or” conjunction. Implication and aggregation are “min” and “max” methods, respectively. Centroid defuzzification method is used. Table 1 shows the design parameters of the MFLC used in the study. Input and output variables are assigned ranges of [-1010] and [-200200], respectively. Operating range of the process and final control element constraint on maximum flow F1 determine the ranges of the input and output variables, respectively. Each variable range is divided into five fuzzy subsets as shown in Table 1. The fuzzy subsets of the input variable e are large negative [a1b1c1d1]; small negative [a2b2c2]; zero [a3b3c3]; small positive [a4b4c4]; and large positive [a5b5c5d5]. The fuzzy subsets of the output variable flow rate are negative high [a6b6c6d6]; negative medium [a7b7c7d7]; zero [a8b8c8]; positive medium [a9b9c9d9]; and positive high [a10b10c10d10]. Triangular and trapezoidal membership functions are used as shown in the table. The total number of parameters is 36. The values of the adaptation gains γ1 and γ2 used are 0.0000004 and 0.000005, respectively. Increasing the values of the adaptation gains has resulted in an overshoot and damped oscillation in the response during simulation study. As zero overshoot is desired in this paper, larger values are not used for the adaptation gains even though larger adaptation gains reduce the rise time.
Design parameters of the MFLC.
Variable
Range
Fuzzy subset
Membership function (MF)
Parameters of MF
Input
[-1010]
Large negative
Trapezoidal
[-14.5-10-2.5-0.005]
Small negative
Triangular
[-2.5-0.0050]
Zero
Triangular
[-0.00500.005]
Small positive
Triangular
[00.0052.5]
Large positive
Trapezoidal
[0.0052.51014.5]
Output
[-200200]
Negative high
Trapezoidal
[-290-200-150-120]
Negative medium
Trapezoidal
[-150-120-200]
Zero
Triangular
[-20020]
Positive medium
Trapezoidal
[020120150]
Positive high
Trapezoidal
[120150200290]
Fuzzy rules are framed based on the following discussion. When e is large positive the controller output has to be reduced fast in order to have faster response. This is achieved by moving the control valve towards fully closed condition. This will reduce error fast and, as it becomes smaller, the control valve position moves towards the new steady-state position gradually so that y converges with ym asymptotically. Similarly, when e is large negative, the controller output is increased fast by moving the control valve towards fully open condition. As error moves towards zero, the control valve position also moves gradually from fully open to the new steady-state position. If e is zero, then control valve position needs no change.
The fuzzy if-then rules for the FMRAC are as follows.
If error is large positive, then flow rate is negative high.
If error is small positive, then flow rate is negative medium.
If error is zero, then flow rate is zero.
If error is small negative, then flow rate is positive medium.
If error is large negative, then flow rate is positive high.
5.2.3. GA Based Fuzzy Modified MRAC
RGA is used to optimize the membership function parameters of the FMRAC that is designed in the previous subsection. The steps referred to in Section 3 are followed to ensure the optimality of the solution. Since there are 36 parameters, each solution has 36 dimensions. A six-digit floating-point number represents each membership function parameter. Maximum fitness is obtained with the assignment of the values 100, 50, 0.8, and 0.02 to number of generations, population size, crossover rate, and mutation rate, respectively. The stopping criterion used here is the number of generations.
Figure 5 shows the pattern of convergence of the best fitness value and the average fitness value. At the twenty-ninth generation, the average fitness approaches the best fitness and since then they stay close to each other. After 50 generations, the best fitness does not change considerably for many more generations. Hence, the estimate can be said to be optimal or near-optimal. The membership function parameters of the FMRAC are substituted by the optimal values found by the application of RGA, and the resultant controller is called GAFMRAC. The optimal values of the membership function parameters are shown in Table 2.
Optimal parameters of the MFLC by using RGA.
Variable
Range
Fuzzy subset
Membership function (MF)
Optimal parameters of MF
Input
[-1010]
Large negative
Trapezoidal
[-10-4.91-2.58-0.00501]
Small negative
Triangular
[-2.57-0.00430]
Zero
Triangular
[-0.0049500.00573]
Small positive
Triangular
[00.004022.52]
Large positive
Trapezoidal
[0.00592.455.0910]
Output
[-200200]
Negative high
Trapezoidal
[-290-202.28-149.95-120.62]
Negative medium
Trapezoidal
[-150.27-119-18.860]
Zero
Triangular
[-15.05015.31]
Positive medium
Trapezoidal
[015.67120.15150.37]
Positive high
Trapezoidal
[119.8150.5203.5290]
Convergence of fitness in GAFMRAC.
5.3. Performance Analysis of the Proposed Controller
In order to carry out performance analysis of different controllers, a step input of size 50 is applied to the command signal uc and the responses are shown in Figure 6. The figure consists of the responses of the reference model, GAMMRAC, FMRAC, and GAFMRAC. Figure 7 shows the zoomed in view of Figure 6 in order to be able to appreciate the difference among the transient responses of the three controllers. The proposed controller of GAFMRAC performs the best as far as the performance criteria are concerned. The performance criteria are calculated for all the three controllers and shown in Table 3. Rise time Tr, settling time Ts, and mean square error (MSE) are the three performance indices considered in the paper. In this context, Ts is the time it takes for the process output y to reach within ±5% of its final value and stay there. Rise time Tr is the time it takes for the plant output y to reach 90% of the final value from 10% of the same. The main objective is to make the process follow the reference model as closely as possible with minimal MSE.
Comparison of step responses of different controllers.
Performance index
Reference model
GAMMRAC
FMRAC
GAFMRAC
Tr (s)
671.6
671.7
671.6
671.6
Ts (s)
1166.8
1170
1169.8
1166.8
MSE
—
1.6e-6
1.3e-7
5.1e-8
Response of the three controllers for a step change in set point.
A zoomed in view of the response of the three controllers shown in Figure 6.
The shortest rise time of 671.6 seconds enables the proposed controller to immediately adapt to any change in the reference level of the process. Usually, the shorter rise time would result in larger overshoot and accompanying oscillations. However, the output of the proposed controller has no overshoot at all and it has the lowest value of MSE. It has reduced the MSE by 60.4 percent when compared to FMRAC.
During the step response analysis, the process inputs in the cases of all the three controllers remain well within the bounds and are shown in Figure 8. The outputs of all the controllers oscillate initially before settling at a final value. Table 4 shows the statistical analysis of the controller outputs. The maximum and the standard deviation are the least in the case of GAFMRAC. Thus, it can be considered better than the other two controllers in terms of stringent requirements of the final control element. It shows that the proposed controller can be implemented to control the hybrid tank process.
Statistical analysis of different controller outputs.
Numerical descriptor
GAMMRAC
FMRAC
GAFMRAC
Maximum (lph)
52.6
43.6
43.5
Mean (lph)
29.2
29.2
29.2
Standard deviation (lph)
4.6
4.0
3.9
Process input for a step input of 50 to the command signal.
In order to check the robustness of the proposed controller, different types of command signals and different reference models are applied and the results are discussed below. Figure 9 shows the responses of the three controllers to a sinusoidal command signal when GM1(s) is used as the reference model. The values of MSE for the GAMMRAC, FMRAC, and GAFMRAC are 5.04e-7, 3.79e-8, and 2.26e-8, respectively. The GAFMRAC gives 40.4 percent and 95.5 percent reduction in MSE when compared with GAFMRAC and GAMMRAC. Figure 10 shows the responses of the three controllers to a step input followed by a sinusoidal signal. The values of MSE for the GAMMRAC, FMRAC, and GAFMRAC are 2.54e-6, 4.65e-7, and 1.84e-7, respectively. In this case, the GAFMRAC gives 60.4 percent and 92.8 percent reduction in MSE when compared with GAFMRAC and GAMMRAC.
Response of the three controllers to a sinusoidal command signal.
Response of the three controllers to a step input followed by a sinusoidal input.
Two different reference models whose parameters are very different from that of the process are used in order to test the robustness of the proposed controller. Figure 11 shows the responses of the three controllers to a sinusoidal command signal when a reference model whose transfer function is given in (13) is used. The values of MSE for GAMMRAC, FMRAC, and GAFMRAC are 6.11e-8, 1.61e-8, and 8.12e-9, respectively. The reference model has larger steady-state gain and smaller time constant when compared to that of the process model:
(13)GM2(s)=0.338000s2+600s+1.
Response of the three controllers to a sinusoidal command signal when GM2(s) is used as the reference model.
Figure 12 shows the responses of the three controllers to a square wave command signal when a first order transfer function given in (14) is used as the reference model. The values of MSE for MRAC, FMRAC, and GAFMRAC are 9.69e-5, 1.89e-5, and 1.82e-5, respectively. The reduction in MSE in the case of the GAFMRAC is 3.8% when compared to FMRAC:
(14)GM3(s)=0.2500s+1.
Response of the three controllers for a square wave command signal when GM3(s) is used as the reference model.
From the analysis of the results, it can be said that the proposed GAFMRAC gives very good transient and steady-state performance for any type of command signal. It is also shown that even when the reference model parameters are very different from the process model parameters, the proposed controller gives stable response and the process shows improved transient and steady-state performance when compared to the performance of GAMMRAC and FMRAC. Even if the plant parameters change due to aging and any other external disturbances, and consequently the process parameters and the reference model parameters become very different from each other, the process will faithfully follow the reference model.
6. Conclusions
The coupled tank process used in this study is a nonlinear process. Hence, an adaptive controller like MRAC is preferred compared to a PID controller in the study. For the coupled tank process, a GAMMRAC, an FMRAC, and a GAFMRAC are designed and their performances are compared. The application of RGA to fine-tune the membership function parameters of the FMRAC has given even better transient and steady-state performance than that of the FMRAC. The simulation results establish that the proposed GAFMRAC reduces MSE better when compared to FMRAC and GAMMRAC. It has also been shown that the proposed controller performs very well for different types of command signals like step, sinusoid, square, and step followed by a sinusoid. It has also performed very well, even when the reference model parameters are very different from the process model parameters indicating the robustness of the proposed design. The designed GAFMRAC can be used for the control of other nonlinear processes as well. Future work aims to adapt the values of the adaptation gains γ1 and γ2 by using fuzzy logic in order to obtain improved transient performance.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
PanH.WongH.KapilaV.de QueirozM. S.Experimental validation of a nonlinear backstepping liquid level controller for a state coupled two tank system2005131274010.1016/j.conengprac.2003.12.0192-s2.0-4344671584HagglundT.1994Triangle Park, NC, USAInstrument Society of America ResearchStephanopoulosG.1984Englewood Cliffs, NJ, USAPrentice HallAsan MohideenK.SaravanakumarG.ValarmathiK.DevarajD.RadhakrishnanT. K.Real-coded Genetic Algorithm for system identification and tuning of a modified Model Reference ADAptive Controller for a hybrid tank system20133763829384710.1016/j.apm.2012.08.019MR30205352-s2.0-84872600240AstromK. J.WittenmarkB.1989Addison-WesleyCartesD.WuL.Experimental evaluation of adaptive three-tank level control200544228329310.1016/S0019-0578(07)60181-52-s2.0-19544373312LiuT.-H.HsuH.-H.Adaptive controller design for a synchronous reluctance motor drive system with direct torque control20071581582410.1049/iet-epa:200700562-s2.0-34548007865MillerD. E.MansouriN.Model reference adaptive control using simultaneous probing, estimation, and control20105592014202910.1109/TAC.2010.2042983MR27224742-s2.0-77956533674GoodwinG. C.MayneD. Q.Continuous-time stochastic model reference adaptive control199136111254126310.1109/9.100934MR11304952-s2.0-0026258670DattaA.IoannouP. A.Performance analysis and improvement in model reference adaptive control199439122370238710.1109/9.362856MR13375632-s2.0-0028694983MillerD. E.DavisonE. J.An adaptive controller which provides an arbitrarily good transient and steady-state response1991361688110.1109/9.62269MR10842472-s2.0-0025755076MohideenK. A.ValarmathiK.Fuzzy modified model reference adaptive controller for improved transient responseProceedings of the International Conference on Power, Energy and Control (ICPEC '13)February 2013Dindigul, India45445710.1109/ICPEC.2013.65276992-s2.0-84880914046ChangW.Nonlinear system identification and control using a real-coded genetic algorithm200731354155010.1016/j.apm.2005.11.0242-s2.0-33750810348HuB.MannG. K. I.GosineR. G.New methodology for analytical and optimal design of fuzzy PID controllers19997552153910.1109/91.7979772-s2.0-0033204112ValarmathiK.DevarajD.RadhakrishnanT. K.Real-coded genetic algorithm for system identification and controller tuning20093383392340110.1016/j.apm.2008.11.0062-s2.0-64249083040EschelmanL. J.SchafferJ. D.Real-coded genetic algorithms and interval-schemata19932San Mateo, Calif, USAMorgan Kaufman187202HerreraF.LozanoM.VerdegayJ. L.Tackling real-coded genetic algorithms: operators and tools for behavioural analysis199812426531910.1023/A:10065049011642-s2.0-0032136585BanerjeeS.ChakrabartyA.MaityS.ChatterjeeA.Feedback linearizing indirect adaptive fuzzy control with foraging based on-line plant model estimation20111143441345010.1016/j.asoc.2011.01.0162-s2.0-79954588609SalehiS.ShahrokhiM.Adaptive fuzzy approach for H∞ temperature tracking control of continuous stirred tank reactors20081691101110810.1016/j.conengprac.2007.12.005EliasiH.DaviluH.MenhajM. B.Adaptive fuzzy model based predictive control of nuclear steam generators2007237666867610.1016/j.nucengdes.2006.08.0072-s2.0-33846917163ZuperlU.CusF.MilfelnerM.Fuzzy control strategy for an adaptive force control in end-milling2005164-1651472147810.1016/j.jmatprotec.2005.02.1432-s2.0-17844394962RossT. J.20052ndJohn Wiley & SonsGoldbergD. E.1989New York, NY, USAAddison-WesleyRajasekaranS.Vijayalakshmi PaiG. A.2003New Delhi, IndiaPrentice Hall