Journal of Control Science and Engineering

Journal of Control Science and Engineering / 2012 / Article
Special Issue

Adaptive Control Theory and Applications

View this Special Issue

Research Article | Open Access

Volume 2012 |Article ID 846458 |

DongBin Lee, Peiman Naseradinmousavi, C. Nataraj, "Nonlinear Dynamic Model-Based Adaptive Control of a Solenoid-Valve System", Journal of Control Science and Engineering, vol. 2012, Article ID 846458, 13 pages, 2012.

Nonlinear Dynamic Model-Based Adaptive Control of a Solenoid-Valve System

Academic Editor: Lili Ma
Received06 Dec 2011
Revised08 Mar 2012
Accepted30 Mar 2012
Published27 Jun 2012


In this paper, a nonlinear model-based adaptive control approach is proposed for a solenoid-valve system. The challenge is that solenoids and butterfly valves have uncertainties in multiple parameters in the nonlinear model; various kinds of physical appearance such as size and stroke, dynamic parameters including inertia, damping, and torque coefficients, and operational parameters especially, pipe diameters and flow velocities. These uncertainties are making the system not only difficult to adjust to the environment, but also further complicated to develop the appropriate control approach for meeting the system objectives. The main contribution of this research is the application of adaptive control theory and Lyapunov-type stability approach to design a controller for a dynamic model of the solenoid-valve system in the presence of those uncertainties. The control objectives such as set-point regulation, parameter compensation, and stability are supposed to be simultaneously accomplished. The error signals are first formulated based on the nonlinear dynamic models and then the control input is developed using the Lyapunov stability-type analysis to obtain the error bounded while overcoming the uncertainties. The parameter groups are updated by adaptation laws using a projection algorithm. Numerical simulation results are shown to demonstrate good performance of the proposed nonlinear model-based adaptive approach and to compare the performance of the same solenoid-valve system with a non-adaptive method as well.

1. Introduction

In order to achieve advanced automation [1] in systems such as marine vessels or ship-based machinery system [2], solenoid actuators and valves are often used [3] to increase survivability and capability. One typical type of actuator driven by solenoids is shown in Figure 1, which is operated by the electromagnetic force. The electric-driven solenoid valve system [4] and its sophisticated control can provide high levels of automation in large systems. The useful function of the solenoid-valve, once an electrical signal (current or voltage) is applied, is to activate a mechanical motion such as displacement or rotation via the solenoid magnetic forces and torques. The proportional solenoids normally require integrated electronics for controlling the plunger to give such a signal. Hydrodynamic torque of a butterfly valve comprises the core knowledge of fluid valve system design [5], and it is known that most of the valves in real systems have strongly nonlinear characteristics between the force and displacement [6, 7]. The use of an intelligent approach [8, 9] such as adaptive, robust, optimal, or nonlinear control of the actuator-valve machinery systems will benefit a wide spectrum of nonlinear systems, compensating for nonlinearities [10] and dynamic characteristics. This approach will not only decrease the amount of cost and casualties but also improve the performance of the mechatronic system. To investigate the particular application, it is important to emphasize the nonlinear dynamic modeling analysis of such actuator-valve systems because the accuracy and reliability of these systems depend highly on the mathematical system modeling [11] and its validation. In [12], the authors developed and analyzed the nonlinear dynamic model of a solenoid-valve system; the reader is also referred to [13, 14] for recent modeling and analysis of solenoid actuators.

This paper will focus on model-based nonlinear adaptive control of an actuator-butterfly valve. The solenoid-valve system is described based on the exact model knowledge of the system. Figure 1 shows the integrated system, which consists of an electric-driven solenoid and a butterfly valve. The valve operates by solenoids that use a magnetic coil to move a movable plunger connected with the valve stem by means of a gear train and linkage. The control input is designed by substituting the current signal from the model of the electromagnetic force, pulling the plunger, and then controlling the angular position of the butterfly valve. The system has uncertainties in multiple parameters in the dynamic model, which requires the system to continuously adjust to the environment and consequently requires adaptation for sustainability and capability. The integrated system is highly nonlinear in addition to its parameter uncertainties. Hence, an adaptation law is proposed [8, 9, 17] and an adaptive control method is developed for the solenoid-valve system in multiple parametric uncertainties. A closed-loop stable controller is designed for the set-point trajectory tracking by introducing a Lyapunov-based stability analysis [9] based on the error signals of the nonlinear solenoid-valve system. The numerical results in the simulation are used for initial verification and performance evaluation.

2. Model-Based Nonlinear System

2.1. System Model

The dynamic equations of motion of the plunger and butterfly valve are given by [15] š‘šĢˆš‘„+šµ1Ģ‡š‘„+š‘˜š‘„=š¹magāˆ’š¹š‘,š½Ģˆš‘Ž+šµ2Ģ‡š›¼=š‘Ÿš‘”š¹š‘āˆ’š‘‡tot,(1) where š‘„(š‘”) is the displacement of the solenoid plunger, š›¼(š‘”) is the angle of butterfly disk, and š‘‡tot(š‘”) is the sum of the hydrodynamic and bearing torques expressed as š‘‡tot=š‘‡š‘+š‘‡ā„Ž. The bearing torque is given as š‘‡š‘=(šœ‹/8)šœ‡š·š‘ š·2š‘Ī”š‘ƒš‘£š¶š‘…(š›¼), where the bearing toque coefficient, š¶š‘…(š›¼)=ī”š¶2šæ+š¶2š·=ī”(1.14āˆšsin((š›¼/90)3180))2+(3āˆšcosš›¼)2, is obtained from the valve modeling (see [5, 15]) and the two subterms, lifting force š¶šæ=1.14āˆšsin((š›¼/90)3180) and drag force š¶š·=3āˆšcosš›¼, are nonlinear functions of the valve angle rotation š›¼(š‘”). The hydrodynamic torque š‘‡ā„Ž is obtained by reviewing three-dimensional hydrodynamic torque coefficient based on [7, 16] as š‘‡ā„Ž=(8/3šœ‹)šœŒš·3š‘š‘‰2š‘‚š‘‡š‘(š›¼)[(š‘‰š½/š‘‰š‘‚)(š›¼)]2, where both š‘‡š‘(š›¼) and (š‘‰š½/š‘‰š‘‚)(š›¼) depend on the closing angle (š›¼) of the butterfly valve and š·3š‘ term is a nonlinear term according to the pipe size. Solving the two equations in (1) with the contact force, š¹š‘(š‘”), and substituting the magnetic force, š¹mag(š‘”), into the equation yields āŽ›āŽœāŽš‘š+š½š‘Ÿ2š‘”āŽžāŽŸāŽ š‘Ÿš‘”Ģˆš‘„+āŽ›āŽœāŽšµ1+šµ2š‘Ÿ2š‘”āŽžāŽŸāŽ š‘Ÿš‘”Ģ‡š‘„+š‘˜š‘Ÿš‘”š‘„=š‘Ÿš‘”š¶2š‘22ī€·š¶1+š¶2š‘„ī€ø2š‘–2āˆ’š‘‡tot,(2) where the magnetic force š¹mag=(š¶2š‘2/2(š¶1+š¶2š‘„)2)š‘–2 used to lift the plunger of the solenoid actuator. The actuator is a current-controlled solenoid [15], proportional to the square of the current š‘–(š‘”), and š¶1, š¶2 are reluctances of the magnetic paths, obtained from the geometry of electric actuator [4, 15]. It is assumed that the pinion and the valve are moving at the same speed, that is, the gear ratio is 1ā€‰:ā€‰1 and š‘„(š‘”)=š‘Ÿš‘”š›¼ is the simple geometric relationship between the displacement of the pinion and the valve angle. Hence, the current source š‘–2(š‘”) is substituted for designing the closed-loop control input, š‘¢(š‘”), and then the following equation is obtained: āŽ›āŽœāŽš‘š+š½š‘Ÿ2š‘”āŽžāŽŸāŽ š‘Ÿš‘”Ģˆš‘„+āŽ›āŽœāŽšµ1+šµ2š‘Ÿ2š‘”āŽžāŽŸāŽ š‘Ÿš‘”Ģ‡š‘„+š‘˜š‘Ÿš‘”š‘„=š‘Ÿš‘”š¶2š‘22ī€·š¶1+š¶2š‘„ī€ø2š‘¢āˆ’š‘‡tot.(3)

For the subsequent controller design, multiplying (3) with the inverse term of the control input, 2(š¶1+š¶2š‘„)2/(š‘Ÿš‘”š¶2š‘2), yields a compact form of the dynamic equation as š‘€(š‘„,šœƒ)Ģˆš‘„+š¶(š‘„,šœƒ)Ģ‡š‘„+š·(š‘„,šœƒ)š‘„=š‘¢āˆ’šµī€·š‘‡š‘+š‘‡ā„Žī€ø,(4) where šœƒ describes a lumped expression of parameters obtained from (3) and (4), where each parameter is shown in Table 2 and the substituted terms are defined as follows: š‘€(š‘„,šœƒ)=āŽ›āŽœāŽš‘š+š½š‘Ÿ2š‘”āŽžāŽŸāŽ š‘Ÿš‘”šµ,šµ=2ī€·š¶1+š¶2š‘„ī€ø2š‘Ÿš‘”š¶2š‘2,š¶(š‘„,šœƒ)=āŽ›āŽœāŽšµ1+šµ2š‘Ÿ2š‘”āŽžāŽŸāŽ š‘Ÿš‘”šµ,š·(š‘„,šœƒ)=š‘˜š‘Ÿš‘”šµ.(5)

2.2. Error Signals Formulation

The following set-point control approach is used. Let š‘„š‘‘(š‘”) define the set-point trajectory and then the error can be defined as š‘’ā‰”š‘„š‘‘āˆ’š‘„,Ģ‡š‘’=Ģ‡š‘„š‘‘āˆ’Ģ‡š‘„,Ģˆš‘’=Ģˆš‘„š‘‘āˆ’Ģˆš‘„,(6) where Ģ‡š‘„š‘‘(š‘”) and Ģˆš‘„š‘‘(š‘”) are the first and second time derivatives of š‘„š‘‘(š‘”), which are assumed to be bounded. Premultiplying Ģˆš‘’(š‘”) in the last error signal of (6) with š‘€(š‘„,šœƒ) yields š‘€(š‘„,šœƒ)Ģˆš‘’=š‘€(š‘„,šœƒ)Ģˆš‘„š‘‘āˆ’š‘€(š‘„,šœƒ)Ģˆš‘„. Substituting š‘€(š‘„,šœƒ)Ģˆš‘„ in (4) into the above equation produces š‘€(š‘„,šœƒ)Ģˆš‘’=š‘€(š‘„,šœƒ)Ģˆš‘„š‘‘+š¶(š‘„,šœƒ)Ģ‡š‘„+š·(š‘„,šœƒ)š‘„āˆ’š‘¢+šµī€·š‘‡š‘+š‘‡ā„Žī€ø.(7)

A filtered error signal and its derivative are defined as š‘Ÿā‰”Ģ‡š‘’+šœ†1š‘’,Ģ‡š‘Ÿ=Ģˆš‘’+šœ†1Ģ‡š‘’,(8) where šœ†1āˆˆā„œ+ is a positive adjustable control gain. Multiplying (8) with š‘€(š‘„,šœƒ) and then substituting for š‘€(š‘„,šœƒ)Ģˆš‘’(š‘”) in (7) yields š‘€(š‘„,šœƒ)Ģ‡š‘Ÿ=š‘€(š‘„,šœƒ)Ģˆš‘„š‘‘+š¶(š‘„,šœƒ)Ģ‡š‘„+š·(š‘„,šœƒ)š‘„āˆ’š‘¢+šµī€·š‘‡š‘+š‘‡ā„Žī€ø+š‘€(š‘„,šœƒ)šœ†1Ģ‡š‘’+1š‘Ÿ2š‘”š‘’āˆ’1š‘Ÿ2š‘”š‘’,(9) where the last term š‘’(š‘”)/š‘Ÿ2š‘” is added and subtracted for further development of the control design based on Lyapunovā€™s method.

3. Lyapunov-Based Adaptive Feedback Control

Let š‘‰(š‘”) be a Lyapunov candidate function š‘‰=12ī‚€š‘Ÿš‘‡š‘€š‘Ÿ+š‘’š‘‡š‘’+š‘’š‘‡š›¼š‘’š›¼+ī‚Ī˜š‘‡Ī“āˆ’1ī‚Ī˜ī‚,(10) where the last term of the Lyapunov candidate function, Ī“=š›¾š¼š‘Ć—š‘, is a constant diagonal matrix with the gain value š›¾, š¼š‘Ć—š‘ is a š‘Ć—š‘ identity matrix, and the parameter estimation error, ī‚Ī˜, is defined as ī‚Ī˜=Ī˜āˆ’īĪ˜, where Ī˜āˆˆā„œš‘ is a known constant parameter vector and īĪ˜āˆˆā„œš‘ is the estimated constant parameter vector (see Table 1). Differentiating (10) yields Ģ‡ā€Œš‘‰=š‘Ÿš‘‡š‘€Ģ‡š‘Ÿ+12š‘Ÿš‘‡Ģ‡ā€Œš‘€š‘Ÿ+š‘’š‘‡Ģ‡š‘’+1š‘Ÿ2š‘”š‘’š‘‡Ģ‡š‘’āˆ’ī‚Ī˜š‘‡Ī“āˆ’1Ģ‡ā€ŒīĪ˜,(11) where the time derivative of the inertia matrix is obtained as Ģ‡ā€Œš‘€=(š‘š+š½/š‘Ÿ2š‘”)(4(š¶1+š¶2š‘„)/š‘2)Ģ‡š‘„, where Ģ‡ā€Œšµ=(4(š¶1+š¶2š‘„)/š‘Ÿš‘”š‘2)Ģ‡š‘„, š‘’š‘‡š›¼Ģ‡š‘’š›¼=(1/š‘Ÿ2š‘”)š‘’š‘‡Ģ‡š‘’ as š‘’š›¼ā‰”š›¼š‘‘āˆ’š›¼=š‘„š‘‘/š‘Ÿš‘”āˆ’š‘„/š‘Ÿš‘”=š‘’/š‘Ÿš‘”, Ģ‡š‘’š›¼ā‰”Ģ‡š›¼š‘‘āˆ’Ģ‡š›¼=Ģ‡š‘„š‘‘/š‘Ÿš‘”āˆ’Ģ‡š‘„/š‘Ÿš‘”=Ģ‡š‘’/š‘Ÿš‘”, in which the error signals of the valve angle, š‘’š›¼(š‘”), can be defined using the geometric relationship and Ģ‡ā€Œī‚Ī˜(š‘”)=āˆ’Ģ‡ā€ŒīĪ˜ comes from the definition of ī‚Ī˜.


š‘Š 1 0 1 Ģˆ š‘„ š‘‘ š‘Š 1 1 3 š¶ š‘… ( š›¼ ) š‘„
š‘Š 1 0 2 Ģˆ š‘„ š‘‘ š‘„ š‘Š 1 1 4 š‘‡ š‘ ( š›¼ ) [ š‘‰ š½ ( š›¼ ) / š‘‰ š‘‚ ] 2 š‘„ 2
š‘Š 1 0 3 Ģˆ š‘„ š‘‘ š‘„ 2 š‘Š 1 1 5 š¶ š‘… ( š›¼ ) š‘„ 2
š‘Š 1 0 4 Ģ‡ š‘„ š‘Š 1 1 6 Ģ‡ š‘„ Ģ‡ š‘’
š‘Š 1 0 5 Ģ‡ š‘„ š‘„ š‘Š 1 1 7 Ģ‡ š‘„ š‘’
š‘Š 1 0 6 Ģ‡ š‘„ š‘„ 2 š‘Š 1 1 8 š‘„ Ģ‡ š‘„ Ģ‡ š‘’
š‘Š 1 0 7 š‘„ š‘Š 1 1 9 š‘„ Ģ‡ š‘„ š‘’
š‘Š 1 0 8 š‘„ 2 š‘Š 1 2 0 š‘’
š‘Š 1 0 9 š‘„ 3 š‘Š 1 2 1 šœ† 1 Ģ‡ š‘’
š‘Š 1 1 0 š‘‡ š‘ ( š›¼ ) [ š‘‰ š½ ( š›¼ ) / š‘‰ š‘‚ ] 2 š‘Š 1 2 2 šœ† 1 š‘„ Ģ‡ š‘’
š‘Š 1 1 1 š¶ š‘… ( š›¼ ) š‘Š 1 2 3 šœ† 1 š‘„ 2 Ģ‡ š‘’
š‘Š 1 1 2 š‘‡ š‘ ( š›¼ ) [ š‘‰ š½ ( š›¼ ) / š‘‰ š‘‚ ] 2 š‘„ ā€”ā€”

ī Ī˜ 1 0 1 ī‚Š š‘€ š‘  2 ī š¶ 2 1 / ( ī š¶ 2 ī š‘ 2 ) ī Ī˜ 1 1 3 Ģ‚ ā€Œ š‘‡ š‘ 1 4 ī š¶ 1 / ( Ģ‚ š‘Ÿ š‘” ī š‘ 2 )
ī Ī˜ 1 0 2 ī‚Š š‘€ š‘  4 ī š¶ 1 / ī š‘ 2 ī Ī˜ 1 1 4 Ģ‚ ā€Œ š‘‡ ā„Ž 1 2 ī š¶ 2 / ( Ģ‚ š‘Ÿ š‘” ī š‘ 2 )
ī Ī˜ 1 0 3 ī‚Š š‘€ š‘  2 ī š¶ 2 / ī š‘ 2 ī Ī˜ 1 1 5 Ģ‚ ā€Œ š‘‡ š‘ 1 2 ī š¶ 2 / ( Ģ‚ š‘Ÿ š‘” ī š‘ 2 )
ī Ī˜ 1 0 4 ī š¶ š‘  2 ī š¶ 2 1 / ( ī š¶ 2 ī š‘ 2 ) ī Ī˜ 1 1 6 ī‚Š š‘€ š‘  2 ī š¶ 1 / ī š‘ 2
ī Ī˜ 1 0 5 ī š¶ š‘  4 ī š¶ 1 / ī š‘ 2 ī Ī˜ 1 1 7 ī‚Š š‘€ š‘  2 ī š¶ 1 / ī š‘ 2
ī Ī˜ 1 0 6 ī š¶ š‘  2 ī š¶ 2 / ī š‘ 2 ī Ī˜ 1 1 8 ī‚Š š‘€ š‘  2 ī š¶ 2 / ī š‘ 2
ī Ī˜ 1 0 7 Ģ‚ ā€Œ š‘˜ 2 ī š¶ 2 1 / ( ī š¶ 2 ī š‘ 2 ) ī Ī˜ 1 1 9 ī‚Š š‘€ š‘  2 ī š¶ 2 / ī š‘ 2
ī Ī˜ 1 0 8 Ģ‚ ā€Œ š‘˜ 4 ī š¶ 1 / ī š‘ 2 ī Ī˜ 1 2 0 1 / š‘Ÿ 2 š‘”
ī Ī˜ 1 0 9 Ģ‚ ā€Œ š‘˜ 2 ī š¶ 2 / ī š‘ 2 ī Ī˜ 1 2 1 ī‚Š š‘€ š‘  2 ī š¶ 2 1 / ( ī š¶ 2 ī š‘ 2 )
ī Ī˜ 1 1 0 Ģ‚ ā€Œ š‘‡ ā„Ž 1 2 ī š¶ 2 1 / ( Ģ‚ š‘Ÿ š‘” ī š¶ 2 ī š‘ 2 ) ī Ī˜ 1 2 2 ī‚Š š‘€ š‘  4 ī š¶ 1 / ī š‘ 2
ī Ī˜ 1 1 1 Ģ‚ ā€Œ š‘‡ š‘ 1 2 ī š¶ 2 1 / ( Ģ‚ š‘Ÿ š‘” ī š¶ 2 ī š‘ 2 ) ī Ī˜ 1 2 3 ī‚Š š‘€ š‘  2 ī š¶ 2 / ī š‘ 2
ī Ī˜ 1 1 2 Ģ‚ ā€Œ š‘‡ ā„Ž 1 4 ī š¶ 1 / ( Ģ‚ š‘Ÿ š‘” ī š‘ 2 ) ā€”ā€”

š‘š š½ š‘Ÿ š‘” šµ 1 šµ 2 š‘˜ š‘ Ī” š‘ƒ š‘£
0.1 1 . 0 4 š‘’ āˆ’ 6 1 š‘’ āˆ’ 2 1020 4 š‘’ 2 8 . 8 š‘’ 3 0.5

š¶ 1 š¶ 2 šœŒ š· š‘ š‘‰ š‘‚ šœ‡ š· š‘  š‘
1 . 5 7 š‘’ 6 6 . 3 2 š‘’ 8 1 š‘’ 3 5 ī…ž ī…ž āˆ¼ 8 ī…ž ī…ž

3.1. Design of Control Input

Substituting š‘€(š‘„,šœƒ)Ģ‡š‘Ÿ(š‘”) into (11) yields Ģ‡ā€Œš‘‰=š‘Ÿš‘‡āŽ›āŽœāŽš‘€Ģˆš‘„š‘‘+š¶Ģ‡š‘„+š·š‘„āˆ’š‘¢+šµī€·š‘‡š‘+š‘‡ā„Žī€ø+š‘€šœ†1Ģ‡š‘’+1š‘Ÿ2š‘”š‘’āŽžāŽŸāŽ āˆ’(Ģ‡š‘’+š‘’)š‘‡š‘’š‘Ÿ2š‘”+12š‘Ÿš‘‡Ģ‡ā€Œš‘€š‘Ÿ+š‘’š‘‡Ģ‡š‘’+1š‘Ÿ2š‘”š‘’š‘‡Ģ‡š‘’āˆ’ī‚Ī˜š‘‡Ī“āˆ’1Ģ‡ā€ŒīĪ˜,(12) where the last term in (9) premultiplied by š‘Ÿ(š‘”) came out of the parenthesis in (12) and is used for the definition of š‘Ÿ(š‘”) in (8). Then, combining the parameterized terms in (12) and substituting them into š‘ŠĪ˜ yields š‘ŠĪ˜=š‘€Ģˆš‘„š‘‘+š¶Ģ‡š‘„+š·š‘„+šµī€·š‘‡š‘+š‘‡ā„Žī€ø+š‘€šœ†1Ģ‡š‘’+1š‘Ÿ2š‘”š‘’+12Ģ‡ā€Œš‘€š‘Ÿ,(13) where š‘Š(Ģˆš‘„š‘‘,Ģ‡š‘„,š‘„,š‘Ÿ,š›¼,Ģ‡š‘’,š‘’)āˆˆā„œ1Ɨš‘ is a known regression vector, which is shown in the left side of Table 1 via the process given later (see (17)) and Ī˜ as the nominal value of the lumped parameter vector. Rearranging (12) produces Ģ‡ā€Œš‘‰=š‘Ÿš‘‡(š‘ŠĪ˜āˆ’š‘¢)āˆ’š‘’š‘‡š‘’š‘Ÿ2š‘”+š‘’š‘‡Ģ‡š‘’āˆ’ī‚Ī˜š‘‡Ī“āˆ’1Ģ‡ā€ŒīĪ˜,(14) where Ģ‡š‘’š‘‡š‘’/š‘Ÿ2š‘” is canceled in the last second term in (12), having the opposite sign because they are scalar, Ģ‡š‘’š‘‡š‘’=š‘’š‘‡Ģ‡š‘’.

The control input can be designed based on Lyapunov stability analysis, making Ģ‡ā€Œš‘‰ negative definite to be shown in the end, as š‘¢=š‘ŠīĪ˜+š‘˜1š‘Ÿ+š‘’,(15) where š‘Ÿ(š‘”) is a feedback error term, š‘˜1 is a positive constant as the control gain, š‘’(š‘”) is another feedback error term added to cancel the term having the opposite sign, š‘’š‘‡Ģ‡š‘’, outside the parenthesis by utilizing the definition of š‘Ÿ(š‘”) given in (8), and š‘ŠīĪ˜ captures the uncertainties associated with the elements of š‘€, š¶, š·, šµ(š‘‡š‘1+š‘‡ā„Ž1), š‘€šœ†1, 1/š‘Ÿ2š‘”, and Ģ‡ā€Œš‘€, which is defined as š‘ŠīĪ˜=ī‚Šš‘€Ģˆš‘„š‘‘+īš¶Ģ‡š‘„+īš·š‘„+īšµī‚µĢ‚ā€Œš‘‡š‘1š¶š‘…(š›¼)+Ģ‚ā€Œš‘‡ā„Ž1š‘‡š‘(š›¼)ī‚øš‘‰š½š‘‰š‘‚(š›¼)ī‚¹2ī‚¶+ī‚Šš‘€šœ†1Ģ‡š‘’+1Ģ‚š‘Ÿ2š‘”š‘’+12Ģ‡ā€Œī‚Šš‘€š‘Ÿ,(16) where the estimated parameter sets are given as ī‚Šš‘€(š‘„,šœƒ)=(īš‘š+Ģ‚ā€Œš½/Ģ‚š‘Ÿ2š‘”)Ģ‚š‘Ÿš‘”īšµ, where īšµ=2(īš¶1+īš¶2š‘„)2/Ģ‚š‘Ÿš‘”īš¶2īš‘2, īš¶(š‘„,šœƒ)=(īšµ1+īšµ2/Ģ‚š‘Ÿ2š‘”)Ģ‚š‘Ÿš‘”īšµ, īš·(š‘„,šœƒ)=Ģ‚ā€Œš‘˜Ģ‚š‘Ÿš‘”īšµ, Ģ‚ā€Œš‘‡š‘1=(šœ‹/8)īšœ‡īš·š‘ īš·2š‘Ī”īš‘ƒš‘£, and Ģ‚ā€Œš‘‡ā„Ž1=(8/3šœ‹)Ģ‚šœŒīš·3š‘īš‘‰2š‘‚. In order to develop the estimate parameter vector īĪ˜ in (16), we need to first define the regression and the estimate terms. Thus, the first term, ī‚Šš‘€Ģˆš‘„š‘‘, in (16) can be defined as ī‚Šš‘€Ģˆš‘„š‘‘=āŽ›āŽœāŽīš‘š+Ģ‚ā€Œš½Ģ‚š‘Ÿ2š‘”āŽžāŽŸāŽ Ģ‚š‘Ÿš‘”2ī‚€īš¶1+īš¶2š‘„ī‚2Ģ‚š‘Ÿš‘”īš¶2īš‘2Ģˆš‘„š‘‘=ī‚Šš‘€š‘ 2īš¶21īš¶2īš‘2Ģˆš‘„š‘‘+ī‚Šš‘€š‘ 4īš¶1īš‘2š‘„Ģˆš‘„š‘‘+ī‚Šš‘€š‘ 2īš¶2īš‘2š‘„2Ģˆš‘„š‘‘=ī‚ƒĢˆš‘„š‘‘š‘„Ģˆš‘„š‘‘š‘„2Ģˆš‘„š‘‘ī‚„ī„æī…€ī…€ī…€ī…€ī…€ī…€ī…€ī…€ī…€ī…€ī…€ī…ƒī…€ī…€ī…€ī…€ī…€ī…€ī…€ī…€ī…€ī…€ī…€ī…Œš‘Š101āˆ¼š‘Š103āŽ”āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ¢āŽ£ī‚Šš‘€š‘ 2īš¶21īš¶2īš‘2ī‚Šš‘€š‘ 4īš¶1īš‘2ī‚Šš‘€š‘ 2īš¶2īš‘2āŽ¤āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ„āŽ¦ī„æī…€ī…€ī…€ī…€ī…€ī…€ī…€ī…ƒī…€ī…€ī…€ī…€ī…€ī…€ī…€ī…ŒīĪ˜101āˆ¼īĪ˜103,(17) where ī‚Šš‘€š‘ =(īš‘š+Ģ‚ā€Œš½/Ģ‚š‘Ÿ2š‘”),īš¶š‘ =(īšµ1+īšµ2/Ģ‚š‘Ÿ2š‘”) (this is shown in Table 1), and Ģ‚š‘Ÿš‘”, īš¶1, or īš¶2 are canceled. š‘Š101āˆ¼š‘Š103 are the measurable regression terms and īĪ˜101āˆ¼īĪ˜103 is the estimated parameters, defined in Table 1, respectively. Similarly to (17), the rest of the terms in (16) are also given in Table 1.

3.2. Online Adaptation Laws for Parameter Updates

The following is constructed to define the known upper and lower bounds but with an unknown parameter of īĪ˜(š‘”) in the sense that īĪ˜_š‘—ā‰¤īĪ˜š‘—(š‘”)ā‰¤īĪ˜š‘—,(18) where īĪ˜š‘—(š‘”) are the estimated parameters as shown in Table 1 and īĪ˜_š‘— and īĪ˜š‘— are the lower and upper bounds of the estimated parameters, respectively, which will be set to the amount of percentage of their true values. The vector Ģ‡ā€ŒīĪ˜š‘—(š‘”) is designed to update using a projection-based algorithm as Ģ‡ā€ŒīĪ˜š‘—=Projī‚†Ī“š‘Šš‘‡š‘Ÿ,īĪ˜š‘—ī‚‡,(19) where Proj{ā‹…} is the projection operator [8] and each lumped parameter is adaptively updated using the adaptation laws [17] for online estimation of unknown parameter as follows: Projī‚†Ģ‡ā€ŒīĪ˜š‘—ī‚‡=Projī‚†Ī“š‘Šš‘‡š‘Ÿ,īĪ˜š‘—ī‚‡=āŽ§āŽŖāŽŖāŽØāŽŖāŽŖāŽ©Ī“š‘Šš‘‡š‘ŸifīĪ˜š‘—>īĪ˜_š‘—,īĪ˜š‘—<īĪ˜š‘—,Ī“š‘Šš‘‡š‘ŸifīĪ˜š‘—=īĪ˜_š‘—,ifĪ“š‘Šš‘‡š‘Ÿ>0,Ī“š‘Šš‘‡š‘ŸifīĪ˜š‘—=īĪ˜š‘—,ifĪ“š‘Šš‘‡š‘Ÿā‰¤0,0elsewhere.(20) Thus, substituting š‘¢(š‘”) in (15) into Ģ‡ā€Œš‘‰ in (13) yields Ģ‡ā€Œš‘‰=š‘Ÿš‘‡š‘Šī‚Ī˜āˆ’š‘Ÿš‘‡š‘˜1š‘Ÿāˆ’(Ģ‡š‘’+š‘’)š‘‡š‘’āˆ’š‘’š‘‡š‘’š‘Ÿ2š‘”+š‘’š‘‡Ģ‡š‘’āˆ’ī‚Ī˜š‘‡Ī“āˆ’1Ģ‡ā€ŒīĪ˜=āˆ’š‘Ÿš‘‡š‘˜1š‘Ÿāˆ’āŽ›āŽœāŽ1+1š‘Ÿ2š‘”āŽžāŽŸāŽ š‘’š‘‡š‘’+ī‚Ī˜š‘‡ī‚†š‘Šš‘‡š‘Ÿāˆ’Ī“āˆ’1Ģ‡ā€ŒīĪ˜ī‚‡.(21) Here, š‘Šī‚Ī˜ is defined as š‘Šī‚Ī˜=ī‚‹š‘€ī€·Ģˆš‘„š‘‘+šœ†1Ģ‡š‘’ī€ø+ī‚š¶Ģ‡š‘„+ī‚š·š‘„+ī‚šµĢƒā€Œš‘‡ā„Ž1š‘‡š‘(š›¼)š‘‰š½š‘‰š‘‚(š›¼)+ī‚šµĢƒā€Œš‘‡š‘1š¶š‘…(š›¼)+Ģƒš‘Ÿš‘”š‘’+12ī‚‹Ģ‡ā€Œš‘€š‘Ÿ,(22) where ī‚‹š‘€=š‘€āˆ’ī‚Šš‘€, ī‚š¶=š¶āˆ’īš¶, ī‚š·=š·āˆ’īš·, ī‚šµĢƒā€Œš‘‡ā„Ž1=šµš‘‡ā„Ž1āˆ’īšµĢ‚ā€Œš‘‡ā„Ž1, ī‚šµĢƒā€Œš‘‡š‘1=šµš‘‡š‘1āˆ’īšµĢ‚ā€Œš‘‡š‘1, Ģƒš‘Ÿš‘”=1/š‘Ÿ2š‘”āˆ’1/Ģ‚š‘Ÿ2š‘”, and ī‚‹Ģ‡ā€Œš‘€=Ģ‡ā€Œš‘€āˆ’ī‚ŠĢ‡ā€Œš‘€, in which Ģ‚ā€Œš‘‡ā„Ž1=(8/3šœ‹)Ģ‚šœŒīš·3š‘īš‘‰2š‘‚ and Ģ‚ā€Œš‘‡š‘1=(šœ‹/8)īšœ‡īš·š‘ īš·2š‘Ī”īš‘ƒš‘£. Actually, ī‚Ī˜(š‘”) can be expressed, for example, using the definitions of īĪ˜ in (16), as follows: owing to the subsequent adaptation law, the time derivative of š‘‰(š‘”) yields a negative definite function except the origin and upper bound by Ģ‡ā€Œš‘‰ā‰¤āˆ’š‘˜1ā€–š‘Ÿā€–2āˆ’āŽ›āŽœāŽ1+1š‘Ÿ2š‘”āŽžāŽŸāŽ ā€–š‘’ā€–2,(23) which can be written as Ģ‡ā€Œš‘‰ā‰¤āˆ’š‘˜2ā€–š‘§ā€–2,(24) where š‘˜2=min{š‘˜1,(1+1/š‘Ÿ2š‘”)} is a positive constant and š‘§=ī€ŗš‘Ÿš‘‡,š‘’š‘‡ī€»š‘‡.(25) Using Barbalatā€™s lemma [18], the set-point tracking error ā€–š‘§(š‘”)ā€–ā†’0, thus ā€–š‘Ÿ(š‘”)ā€–ā†’0 and ā€–š‘’(š‘”)ā€–ā†’0 as š‘”ā†’āˆž.

Remark 1. According to the analysis from (11) to (25) of Ģ‡ā€Œš‘‰, the property of š‘‰, and the control law of (15) with the parameter updates of (19) and the projection-based method of update laws of (20), it is straightforward to derive a conclusion that the tracking error vector š‘§(š‘”) in (25) is driven to zero. Thus, the set-point errors š‘Ÿ(š‘”), š‘’(š‘”), and š‘’š›¼(š‘”) also vanish and the parameter estimation error vector ī‚Ī˜ in (10) is bounded where īĪ˜, defined after (10), is bounded due to the projection-based update method and the constant known parameter, Ī˜. Owing to the bounds of š‘Ÿ(š‘”), š‘’(š‘”) in (8), Ģ‡š‘’(š‘”) is bounded and then š‘„(š‘”), resulting in š›¼(š‘”), and Ģ‡š‘„(š‘”) are bounded, respectively, where all desired trajectories such as š‘„š‘‘(š‘”) and Ģ‡š‘„š‘‘(š‘”) are assumed to be bounded. š‘€(ā‹…), šµ(ā‹…), š¶(ā‹…), and š·(ā‹…) matrices in (4) are bounded because šœƒ (due to Ī˜) and š‘„(š‘”) are bounded, and Ģ‡ā€Œšµ and Ģ‡ā€Œš‘€ after (11) are also bounded owing to the bounds of Ģ‡š‘„(š‘”). š‘‡tot is thus bounded because š›¼(š‘”) is bounded. Hence, š‘Š(ā‹…) and the control input š‘¢(š‘”) are bounded and, thus, the current is bounded. This leads to the boundedness of Ģˆš‘„(š‘”) in the dynamic model given in (4), which enables the boundedness of Ģˆš‘’(š‘”) in (6), and then the set-point tracking error dynamics Ģ‡š‘Ÿ(š‘”) in (8) is bounded. Therefore, we can conclude that all signals are bounded.

Remark 2. š‘ŠīĪ˜ in (16) where the known regression š‘Š(Ģˆš‘„š‘‘,Ģ‡š‘„,š‘„,š‘Ÿ,š›¼,Ģ‡š‘’,š‘’)āˆˆā„œ1Ɨ23 terms and the parameter estimates īĪ˜āˆˆā„œ23 are given in this system as in Table 1.

4. Simulation Results

Based on the dynamic model in (4), the numerical simulation is performed to verify the proposed controllers with consistently changing parameter values. The parameter values can be divided into two categories: operational values such as š·š‘ and š‘‰š‘‚ and uncertain values such as šµ1,šµ2, and šœ‡. After the flow velocity š‘‰š‘‚ is kept on 3.7 [m/s], the pipe diameter š·š‘ shown in Table 2 is used to vary from 5 inches to 8 inches as š·š‘={5.0,6.0,7.0,8.0}[in], where š·š‘={0.1270,0.1524,0.1778,0.2032}[m], by assuming that the pipe or transmission lines are different according to their applications. With each š·š‘ size, the variations of all the parameters are set to 30% for the simulation results provided here. The determined parameter vectors š‘‡š‘(š›¼) and [š‘‰š½/š‘‰š‘‚](š›¼) of the butterfly valve model are borrowed from [16] for the simulation given as look-up tables from the experimental data. The amount of upper and lower variation of the unknown parameter sets is 30% of their real values. The control gains were chosen selectively as š›¾=10, š‘˜1=250, and šœ†1=1.0 for all cases. MATLAB and Simulink are used for the simulation. The desired set-point distance of the solenoid, š‘„š‘‘(š‘”), is given as š“ā‹…(1āˆ’š‘’āˆ’šµš‘”), where š“=0.0148 and šµ=5.

A typical parameter set for this simulation is given by Table 2. Figures 2 and 3 show the rotation of the valve angle and the actual displacement of the plunger, respectively. This adaptive control approach shows better results compared to the results obtained from the previous research [15] using nonadaptive method, which are shown in Figures 4 and 5 for the displacement of the plunger and the angle of the butterfly valve. The developed mathematical model is the same as that of adaptive method based on the nonlinear models and the new specific approach in this paper is that the error signals are formulated by introducing the desired trajectory to reducing the displacement error in the adaptive method while overcoming the complicated parametric uncertainties. And also the response of the adaptive method shows that the motion process is more consistent in both variables, showing overdamped phenomena and short rising times, and then it quickly reaches the steady state. The actual movement is smoother and faster while all the parameters are varying 30% in each flow velocity.

Figures 6 and 7 are the filtered error š‘Ÿ(š‘”), defined in (8), and set-point error š‘’(š‘”), given in (6), respectively. Figure 8 shows the rate of the displacement to Figure 3. From Figure 8, all figures are presented without the notation of the adaptive approach. Figure 9 shows the control input š‘¢(š‘”) of the solenoid actuator for each pipe diameter, š·š‘, designed in (15) with the nonlinear adaptive controller by substituting the square of the current, š‘–2(š‘”). The electromagnetic force š¹mag(š‘”) given in (2) and the total torque š‘‡š‘”(š‘”) given in (1) by summing up the hydrodynamic and bearing torques are plotted in Figures 10 and 11, respectively. The torque coefficient š‘‡š‘(š‘”) shown in Figure 12 and the ratio of input and output jet velocities shown in Figure 13 are changed for every 5āˆ˜ of the butterfly angle. As the strokes are increased by the control input, the angles get larger and then accordingly the value of the inlet jet velocity increases, which affects the slower motion of the strokes and angles but the variables (strokes and angles) reach the steady state and the ratio as well.

The challenge is that most parameter terms in Table 1 are combined and lumped together due to the dynamic model in the presence of parameter uncertainty and the model is complicated owing to the control objectives, set-point regulation, and parametric adaptation. As given in the right side of Table 1, the unknown bounded parameter estimate vector īĪ˜āˆˆā„œ23 is shown in Figures 14, 15, 16, 17, 18, and 19. It can be seen that the parameters such as 101, 110, 111, 113, 116, 118, and 121 are quickly updated in the form of premultiplication by its regression term, š‘Š and go to steady state. Thus, they are more parametric-centric terms and related to the filtered error š‘Ÿ(š‘”). The parameters such as 102, 103, 104, 105, 117, 119, and 122 are updated according to the motion of the angle and stroke of the solenoid-butterfly valve system. Some parameters such as 106, 112, 114, and 115 are not updated but other parameters such as 107, 108, and 120 keep getting updated. Any parameters changing their values must be related to the control objective because these terms are incorporated into the control input in the form of the estimates.

5. Conclusion and Future Work

For developing advanced automation systems such as ship-based hydraulic systems, a typical solenoid-butterfly valve, which is driven by electromagnetic, fluid mechanics, and hydrodynamic forces and torques, is chosen as a continuation of previous research and an adaptive stable control approach with adaptation laws is developed accounting for uncertainties in multiple parameters on the nonlinear dynamic model. A stable adaptive controller of the solenoid-valve system is designed positioning the angle of the butterfly valve via a Lyapunov-based approach. The approach yields bounded error while adapting to the environment in the presence of complex uncertainties such as different physical appearances, uncertain parameters, operational characteristics, and parametric nonlinear dynamic models.

The parameter estimation for the unknown bounded parameters is performed using a projection algorithm whose output yields the upper and lower bounds. Numerical simulation is used to verify the performance of the proposed approach to show its effectiveness by comparing to the same dynamic model without adaptation from the previous research; when compared to the nonadaptive method, the responses of the plunger displacement and the rotating angle are steadier, smoother, and faster. Future work will be focused on demonstrating the results of hardware-in-the-loop or experiments for the nonlinear solenoid valve system as well as applying the suggested adaptive method based on Lyapunov-based control approach to the real-world system. Further research on developing control techniques using robust or optimal method would be continued to overcome nonlinearities such as hysteresis or nonlinear dynamics.


šµ1, šµ2:Damping coefficients of the solenoid and butterfly valve, respectively (Ns/m)
š¶1, C2:Reluctances of the magnetic paths obtained from plunger geometry of solenoid actuator
š¶š‘…(š›¼):Bearing torque coefficient as a function of the valve angle (š›¼)
š·š‘ :Stem diameter (m)
š·š‘:Pipe diameter (inch or m)
š¹mag: Magnetic force (N)
š¹š‘: Contact force or resultant force š¹š‘Ÿ(N)
š‘–:Current of solenoid actuator (A)
š½:Inertia moment (kgm2)
š‘˜:Spring stiffness (N/m)
š‘š:Mass of solenoid plunger (kg)
šœ‡:Friction coefficient of bearing area
š‘:Number of turns of the coil
Ī”š‘ƒš‘£:Valve differential pressure (psi)
š‘:Number of estimates
š‘Ÿš‘”:Radius of pinion gear (m)
šœŒ:Fluid mass density (kg/m3)
š‘¢:Control input (A2)
š‘‡š‘(š›¼):Hydrodynamic torque coefficient as a function of the valve angle (š›¼)
š‘‡š‘, š‘‡ā„Ž:Bearing and hydrodynamic torques, respectively (Nm)
š‘‰š½, š‘‰š‘‚, (š‘‰š½/š‘‰š‘‚)(š›¼)āˆ¶Jet velocity, mean flow velocity, and their ratio as a function of the valve angle (š›¼), respectively.


This research is supported by the Office of Naval Research (N00014-08-1-0435), which the authors gratefully acknowledge. Thanks are in particular due to Mr. Anthony Seman III. The authors would also like to thank Dr. Stephen Mastro and Mr. Frank Ferrese of Naval Surface Warfare Center (NSWC, Philadelphia) for help with many aspects of the paper. They are grateful to the anonymous reviewer for critical comments, which led to substantial improvement of the paper.


  1. A. Seman, ā€œAdaptive automation for machinery control,ā€ in Proceedings of the Office of Naval Research (ONR) Presentation, 2007. View at: Google Scholar
  2. R. Hughes, S. Balestrini, K. Kelly, N. Weston, and D. Mavris, ā€œModeling of an integrated reconfigurable intelligent system (IRIS) for ship design,ā€ in Proceedings of the ASNE Ships & Ship Systems Symposium, 2006. View at: Google Scholar
  3. P. Tich, V. Mark, P. Vrba, and M. Pechoucek, ā€œChilled water system automation,ā€ Rockwell Automation Case Study Report, 2005. View at: Google Scholar
  4. J. R. Brauer, Magnetic Actuators and Sensors, Wiley IEEE Press, Hoboken, NY, USA, 2006.
  5. Z. Leutwyier and H. Dalton, ā€œA CFD study of the flow field, resultant force, and aerodynamic torque on a symmetric disk butterfly valve in a compressible fluid,ā€ Journal of Pressure Vessel Technology, vol. 130, Article ID 021302, pp. 1ā€“10, 2008. View at: Google Scholar
  6. T. Sarpkaya, ā€œOblique impact of a bounded stream on a plane lamina,ā€ Journal of the Franklin Institute, vol. 267, no. 3, pp. 229ā€“242, 1959. View at: Google Scholar
  7. T. Sarpkaya, ā€œTorque and cavitation characteristics of butterfly valves,ā€ Journal of Applied Mechanics, vol. 29, pp. 511ā€“518, 1961. View at: Google Scholar
  8. M. Krstic, I. Kanellakopoulos, and P. Kokotovic, Nonlinear and Adaptive Control Design, Series on Adaptive and Learning Systems for Signal Processing, John Wiley & Sons, New York, NY, USA, 1st edition, 1995.
  9. D. B. Lee, C. Nataraj, and P. Naseradinmousavi, ā€œNonlinear model-based adaptive control of a solenoid-valve system,ā€ in Proceedings of the ASME Dynamic Systems and Control Conference (DSCC '10), pp. 1ā€“8, Boston, Mass, USA, September 2010. View at: Google Scholar
  10. J. Pohl, M. Sethson, P. Krus, and J.-O. Palmberg, ā€œModelling and validation of a fast switching valve intended for combustion engine valve trains,ā€ Proceedings of the Institution of Mechanical Engineers. Part I: Journal of Systems and Control Engineering, vol. 216, no. 2, pp. 105ā€“116, 2002. View at: Publisher Site | Google Scholar
  11. C. Nataraj, Vibration of Mechanical Systems, Cengage, 1st edition, 2011.
  12. P. Naseradinmousavi and C. Nataraj, ā€œNonlinear mathematical modeling of butterfly valves driven by solenoid actuators,ā€ Applied Mathematical Modelling, vol. 35, no. 5, pp. 2324ā€“2335, 2011. View at: Publisher Site | Google Scholar
  13. M. K. Zavarehi, P. D. Lawrence, and F. Sassani, ā€œNonlinear modeling and validation of solenoid-controlled pilot-operated servovalves,ā€ IEEE/ASME Transactions on Mechatronics, vol. 4, no. 3, pp. 324ā€“334, 1999. View at: Publisher Site | Google Scholar
  14. R. R. Chladny, C. R. Koch, and A. F. Lynch, ā€œModeling automotive gas-exchange solenoid valve actuators,ā€ IEEE Transactions on Magnetics, vol. 41, no. 3, pp. 1155ā€“1162, 2005. View at: Publisher Site | Google Scholar
  15. C. Nataraj and P. Mousavi, ā€œNonlinear analysis of solenoid actuators and butterfly valve systems,ā€ in Proceedings of the 14th International Ship Control Systems Symposium, pp. 1ā€“8, Ottawa, Canada, September 2009. View at: Google Scholar
  16. J. Y. Park and M. K. Chung, ā€œStudy on hydrodynamic torque of a butterfly valve,ā€ Journal of Fluids Engineering, Transactions of the ASME, vol. 128, no. 1, pp. 190ā€“195, 2006. View at: Publisher Site | Google Scholar
  17. J. B. Pomet and L. Praly, ā€œAdaptive nonlinear regulation: estimation from the Lyapunov equation,ā€ IEEE Transactions on Automatic Control, vol. 37, no. 6, pp. 729ā€“740, 1992. View at: Publisher Site | Google Scholar
  18. H. K. Khalil, Nonlinear Systems, Prentice Hall, Upper Saddle River, NJ, USA, 3rd edition, 2002.

Copyright © 2012 DongBin Lee et al. 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

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.