AAA Abstract and Applied Analysis 1687-0409 1085-3375 Hindawi Publishing Corporation 871614 10.1155/2014/871614 871614 Research Article Fractional Order Models of Industrial Pneumatic Controllers Razminia Abolhassan 1 http://orcid.org/0000-0002-0286-7244 Baleanu Dumitru 2,3,4 Agarwal Praveen 1 Dynamical Systems & Control (DSC) Research Lab. Department of Electrical Engineering School of Engineering Persian Gulf University P.O. Box 75169, Bushehr Iran pgu.ac.ir 2 Department of Chemical and Materials Engineering Faculty of Engineering King Abdulaziz University P.O. Box 80204 Jeddah 21589 Saudi Arabia kau.edu.sa 3 Department of Mathematics and Computer Sciences Faculty of Arts and Sciences Çankaya University 06530 Ankara Turkey cankaya.edu.tr 4 Institute of Space Sciences Magurele-Bucharest Romania 2014 2422014 2014 17 12 2013 08 01 2014 24 2 2014 2014 Copyright © 2014 Abolhassan Razminia and Dumitru Baleanu. 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.

This paper addresses a new approach for modeling of versatile controllers in industrial automation and process control systems such as pneumatic controllers. Some fractional order dynamical models are developed for pressure and pneumatic systems with bellows-nozzle-flapper configuration. In the light of fractional calculus, a fractional order derivative-derivative (FrDD) controller and integral-derivative (FrID) are remodeled. Numerical simulations illustrate the application of the obtained theoretical results in simple examples.

1. Introduction

Fractional calculus is a powerful mathematical tool with a long history, but its application to engineering and modeling of physical systems has attracted much attention only in recent years . This theory generalizes the classical differentiation and integration into noninteger order ones. It has been found that in interdisciplinary fields, many systems can be described more accurately and more conveniently by fractional differential equations (FDEs). For instance, fractional derivatives have been widely used in the mathematical modeling of viscoelastic materials . The anomalous diffusion phenomena in nonhomogeneous media can be explained by noninteger derivative-based equations of diffusion . Another example for an element with fractional order model is fractance, which is an electrical circuit with noninteger order impedance and has a property that lies between resistance and capacitance . Moreover, it has been shown that the dynamical process of heat conduction can be modeled more adequately via fractional order calculus . In biology, the membranes of the biological cells are proven to have fractional order electrical conductance and are classified among noninteger order systems [6, 7]. In economics, it is known that some finance systems can display fractional order dynamics . For more examples of fractional order models see, for example, [9, 10] and the references therein. As an important industrial controller, fractional order versions of PID controllers have been considered in the literature . More theoretical problems and recent applications can be found in .

Pneumatic controllers are essential parts of the industrial automation systems containing several diaphragms and bellows . There have been a great development in low-pressure pneumatic controllers for industrial control systems in the past decades, and today they are used extensively in industrial processes. The advantages include safety and explosion-proof characteristic, simplicity, low cost, high compliance, and ease of maintenance . For instance, recently, pneumatic servo systems have been applied in several systems, such as industrial robots, rehabilitation tools, and medical and caregiver robots, aspherical glass molding machines, and precise positioning [14, 15].

Recently in  a model of a pneumatic vibration isolation system (PVIS) has been presented. However, a thorough study on the modeling of the pneumatic elements via fractional order differential equations has not been considered yet. Because of compressibility of the air, the control action in practical actuating valves may not be positive; that is, an error may exist in the valve-stem position. To overcome this imperfection, a more accurate model of pneumatic controllers is needed. Using fractional order dynamics for diaphragms, bellows, and other devices that may have sprig-property yields a more accurate model for such industrial systems than the classical integer order models. In particular, relaxation processes deviating from the classical exponential behavior are often encountered in the dynamics of complex materials. In many cases experimentally observed relaxation functions exhibit a stretched exponential decay [17, 18]. Such behavior may be seen in the stress relaxation of viscoelastic materials, such as polymers or critical gels, in the change carrier transport in amorphous semiconductors and in the attenuation of seismic waves . By this new approach, the available hysteresis in the diaphragms, which has been usually disregarded for simplicity in conventional approaches, can be easily modeled in the new approach. Fractional calculus allows for a rigorous and reliable modeling for such systems.

In this paper, in order to consider some real and nonideal features of a pneumatic structure, using a four-parameter fractional derivative Zener model for viscoelastic materials , we present a fractional order model for a nozzle-flapper-relay configuration. Then, under some simplified conditions, we provide a multiorder fractional derivative and integral controller, which considers the memory of the system in a compact form. The numerical simulations support the obtained results.

This paper is organized as follows. Section 2 provides the preliminary background from fractional calculus. In Section 3, the pneumatic and pressure systems are presented. Some illustrative examples are provided in Section 4. Finally, the conclusion remarks are given in Section 5.

2. Backgrounds

In this section, a brief background of fractional calculus is presented. The definition of the fractional integral is the extension of Cauchy formula for evaluating the integration. The q th order fractional integral of function x ( t ) with respect to t is defined by (1) J q x ( t ) = 1 Γ ( q ) 0 t ( t - s ) q - 1 x ( s ) d s , in which Γ ( q ) = 0 e - z z q - 1 d z , q > 0 is the Euler gamma function. In addition, there are some definitions for fractional derivatives such as Riemann-Liouville, Grunwald-Letnikov, and Caputo definitions . For example, the Riemann-Liouville (RL) fractional derivative of order q for function x ( t ) is defined by (2) D 0 RL t q x ( t ) = D m J ( m - q ) x ( t ) , where m - 1 < q < m , m Z + .

The Laplace transform of RL fractional derivative of order q of function x ( t ) is (3) L { D 0 RL t q x ( t ) } = s q X ( s ) - k = 0 m - 1 s k · D 0 RL t q - k - 1 x ( 0 ) , where m - 1 < q < m Z + . As can be seen from (3) for evaluating the Laplace transform of RL fractional derivative operator, the fractional order derivative of the function x ( t ) is needed as initial conditions. This is somewhat meaningless from the physical viewpoint. Therefore, the RL operator is useless in modeling of the physical systems.

Another definition for fractional derivative has been introduced by Caputo: (4) D 0 C t q x ( t ) : = D 0 RL t - ( m - q ) D m x ( t ) = 1 Γ ( m - q ) 0 t ( t - s ) m - q - 1 x ( m ) ( s ) d s , where m - 1 < q < m Z + .

In contrast to RL operator, the Laplace transform of the Caputo operator needs the integer order derivative of the function as the initial conditions, as shown by the following relation: (5) L { D 0 C t q x ( t ) } = s q X ( s ) - k = 0 m - 1 s q - k - 1 x ( k ) ( 0 ) , where m - 1 < q < m Z + .

3. Fractional Order Model of Pneumatic and Pressure Systems

Many industrial processes and pneumatic controllers involve the flow of a gas or air through connected pipelines and pressure vessels. Thus, it is logical to introduce a parameter for characterizing their specifications. Some quantities in pipelines and pressure vessels can be viewed as resistance ( R ) and capacitance ( C ) used in electrical circuits. Traditional definitions for these characteristics are as follows : (6) R = d ( Δ P ) d q , where d ( Δ P ) is a small change in the gas pressure difference between input and vessel gas and d q is a small change in the gas flow rate. Traditionally, the capacitance of the pressure vessel may be defined by (7) C = d m d p , where m and p are mass of gas stored in the vessel and gas pressure, respectively.

It can be shown that during the change of state of a polytropic process, from isothermal to adiabatic state, the capacitance is constant and can be obtained as  (8) C = V n R gas T , where n is the polytropic exponent, V , R gas , and T are the volume of vessel, gas constant, and absolute temperature, respectively. It is worth noting that in some practical cases, the polytropic exponent n is approximately 1.0–1.2 for gases in uninsulated metal vessels.

3.1. A Fractional Order Model for the Nozzle-Flapper-Relay Controller

A conventional apparatus in the industrial pneumatic control systems is the nozzle-flapper configuration. The nozzle-flapper amplifier converts displacement into a pressure signal. Since typical industrial process control systems require large output power to operate large pneumatic actuating valves, the power amplification of the nozzle-flapper amplifier is usually insufficient. To overcome this problem, a pneumatic relay can be connected to the nozzle-flapper. Therefore, a complete pneumatic amplifier is composed of two stages: nozzle-flapper as the first and pneumatic relay as the second amplifier stages. A schematic diagram of such a configuration is depicted in Figure 1.

Schematic diagram of a force-distance type of pneumatic proportional controller.

Assuming that the relationship between the variation in the nozzle back pressure p b and the variation in the nozzle-flapper distance x is linear, one can write (9) p b = k 1 x , where k 1 is a positive constant. We now develop a new fractional order model for the pneumatic relay containing a diaphragm as mentioned in Section 1 which considers hereditary property of the configuration. The elasticity and damping properties of the diaphragm cannot be measured directly because the shape of the diaphragm in the deflated condition is quite different from that in the inflated condition, and the deflated diaphragm is very difficult to handle due to its flexibility. Therefore, the stiffness of the diaphragm is obtained usually from the experimental data. A simple bleed type of pneumatic relay is depicted in Figure 2.

A simple model of bleed-type relay.

Following , a typical form of the four-parameter fractional derivative Zener model for viscoelastic materials can be written as (10) σ ( t ) + τ β d β d t β σ ( t ) = - E min ε ( t ) - E max τ β d β d t β ε ( t ) , where σ ( t ) is the stress, ε ( t ) is the strain, E max and E min are the maximum and minimum of elastic modulus, τ is the time constant, and β is the exponent of the fractional derivative. Based on the fact that for the pneumatic relay, the back pressure on the diaphragm is proportional to the stress and the displacement is proportional to the strain, (10) reduces to (11) p b ( t ) + τ new β d β d t β p b ( t ) = Ψ min z ( t ) + Ψ max τ new β d β d t β z ( t ) , where p b ( t ) is the back pressure acting on the top diaphragm of the pneumatic relay, z is the resulting displacement of the diaphragm from the equilibrium, Ψ max and Ψ min are proportional to the maximum and minimum stiffness of the diaphragm, respectively, and τ new is the new time constant which is a characteristic of the structure. Applying the Laplace transform on (11) yields (12) P b ( s ) Z ( s ) = Ψ min + Ψ max τ new β s β 1 + τ new β s β , where P b ( s ) and Z ( s ) are the Laplace transform of the back pressure signal p b ( t ) and the diaphragm displacement z ( t ) , respectively. The structure of the diaphragm is such that the relation between control pressure p c ( t ) and z ( t ) is linear; that is, (13) P c ( s ) = k 2 Z ( s ) , where k 2 is a positive constant and P c ( s ) is the Laplace transform of p c ( t ) . Using (9), (12), and (13) one can conclude that (14) P c ( s ) = k 1 k 2 1 + ( τ new s ) β Ψ min + Ψ max ( τ new s ) β    X ( s ) .

For the flapper, since there are two small movements ( e and y ) in opposite directions, one may add up the results of two movements into one displacement x . It can be easily shown that for the flapper movement the following relationship is held: (15) x ( t ) = b a + b e ( t ) - a a + b y ( t ) , where y is the displacement of the end of the flapper due to the bellows expansion or contraction. A usual approach is to consider the bellows action as a linear spring. However, as was developed for the diaphragm, for every practical bellows a simple fractional order model can be considered as follows which we call it a pseudo-spring: (16) A p c ( t ) = k s D 0 C t q y ( t ) , in which A is the effective area of the bellows and k s is the equivalent pseudo-spring constant, that is, the stiffness due to the action of the corrugated side of the bellows. Using (9), (14), (15), and (16), one can develop a block diagram for the nozzle-flapper-relay, as shown in Figure 3.

Block diagram for the fractional order nozzle-flapper-relay.

The transfer function between the control pressure P c ( s ) and the error signal E ( s ) can be obtained using Mason’s formula as follows: (17) P c ( s ) E ( s ) = ( b k s k 1 k 2 s q ( 1 + ( τ new s ) β ) ) × ( ( a + b ) k s s q ( Ψ min + Ψ max ( τ new s ) β ) + a A k 1 k 2 ( 1 + ( τ new s ) β ) ( Ψ min + Ψ max ( τ new s ) β ) ) - 1 .

As can be seen, this controller configuration yields a fractional order proper transfer function. In the following sections, by inserting a restrictor in the feedback path, a fractional order derivative and integral performance are achieved.

Remark 1.

It should be emphasized that in using transfer function (by the definition) all initial conditions must be set to zero. Since we have utilized the Caputo’s fractional derivative, all initial conditions are those which have been considered in the integer order cases. In other words, from (5), in computing the Laplace transform of the Caputo’s derivative, we need integer order derivatives of the signal; that is, L { D 0 C t q x ( t ) } = s q X ( s ) - k = 0 m - 1 s q - k - 1 x ( k ) ( 0 ) in which necessity of x ( k ) ( 0 ) is apparent. So, setting these initial conditions to zero in the development of transfer functions is naïve. Indeed, if the Riemann-Liouville derivative had been used, setting the initial conditions to zero does not have any meaning. Moreover, as there is not any superfluous element in all proposed pneumatic structures, one can conclude that the system is completely characterized by its transfer function which implies that the contributions of the consistent initial conditions can be seen in the response given by the transfer function. Any initial conditions can be regarded in the simulation by setting the desired values in the integrator parts which is an option in Simulink package.

3.2. Pneumatic Controller with Fractional Order Derivative Performance

Consider the pneumatic controller configuration shown in Figure 4, in which a restriction in the negative feedback path modifies the previous controller introduced in (17).

Pneumatic controller with a restriction in the feedback path.

Assuming a small positive step change in the error signal e , the control pressure p c changes almost instantaneously. The restriction R momentarily prevents the feedback bellows from sensing the pressure change p c . Thus, the feedback bellows does not respond momentarily, and the pneumatic actuating valve feels the full effect of the flapper movement. As the time goes on, the feedback bellows expands. The expansion of the bellows critically depends on the type of the material used in the construction of the bellows. The usual materials used in the practical bellows are not purely spring and exhibit something between elastic and semiplastic behavior. The differential equation that could describe this behavior under pressure changes is very complicated. Using a fractional order operator to describe such a pseudo-spring bock is logical, because the order of the operator gives an additional freedom to match the behavior. As can be observed in a practical implementation, the signal p c decays in time , and thus, the controller is of the derivative type or more accurately a proportional plus fractional order derivative. Indeed similar to the methodology discussed in Section 3.1, one can consider a fractional order model for this capacitance-resistance connection as 1 / ( R C s v + 1 ) .

Thus, a simplified block diagram for the controller which is called fractional order derivative controller can be developed as depicted in Figure 5.

A simplified block diagram for fractional order PD controller.

The transfer function of the control pressure with respect to the error signal is (18) P c ( s ) E ( s ) = ( b a + b k 1 k 2 ( 1 + ( τ new s ) β Ψ min + Ψ max ( τ new s ) β ) ) × ( 1 + k 1 k 2 ( 1 + ( τ new s ) β Ψ min + Ψ max ( τ new s ) β ) × ( a a + b ) ( A k s s q ) ( 1 R C s v + 1 ) ( 1 + ( τ new s ) β Ψ min + Ψ max ( τ new s ) β ) ) - 1 = ( k 1 k 2 b ( k s s q ) ( R C s v + 1 ) ( 1 + ( τ new s ) β ) ) × ( ( Ψ min + Ψ max ( τ new s ) β ) ( a + b ) × ( k s s q ) ( R C s v + 1 ) + k 1 k 2 a A ( 1 + ( τ new s ) β ) ) - 1 = ( K ( R C τ new β s q + v + β + R C s q + v + s q + τ new β s q + β ) ) × ( κ ( Ψ max R C τ new β s q + v + β + Ψ max τ new β s q + β + Ψ min R C s q + v + Ψ min s q τ new β ) + η + η τ new β s β ) - 1 in which K = k 1 k 2 k s b , κ = k s ( a + b ) , and η = k 1 k 2 a A .

Remark 2.

Note that if | k 1 k 2 ( ( 1 + ( τ new s ) β ) / ( Ψ min + Ψ max ( τ new s ) β ) ) ( a / ( a + b ) ) ( A / k s s q ) ( 1 / ( R C s v + 1 ) ) | 1 , then (19) P c ( s ) E ( s ) ( b a + b k 1 k 2 ( 1 + ( τ new s ) β Ψ min + Ψ max ( τ new s ) β ) ) × ( k 1 k 2 ( 1 + ( τ new s ) β Ψ min + Ψ max ( τ new s ) β ) × ( a a + b ) ( A k s s q ) ( 1 R C s v + 1 ) ( 1 + ( τ new s ) β Ψ min + Ψ max ( τ new s ) β ) ) - 1 = M ( R C s q + v + s q ) , where M = b k s ( a A ) - 1 . It can be easily seen that this transfer function is a pure fractional order derivative controller with two tunable orders. For this reason, we call it a fractional order D + D controller.

3.3. Pneumatic Controller with Fractional Order Integral Performance

Similar to the previous discussions in Section 3.2, we can develop a pneumatic controller with fractional order integral performance. Consider the controller configuration in Figure 6.

Pneumatic controller with a restriction in the feedback path which depicts an integral performance.

The bellows denoted by I is connected to the control pressure source without any restriction, though a restriction can be implemented for it. The bellows denoted by II is connected to the control pressure source through a restriction denoted by R . A small positive step change in the actuating error, e , will cause the back pressure in the nozzle to change instantaneously, and thus a change in the control pressure, p c , also occurs instantaneously. Due to the restriction of the valve in the path to bellows II, there will be a pressure drop across the valve. As the time goes, the air flows across the valve in such a way that the pressure change in bellows II reaches p c . Thus, bellows II will expand or contract as the time elapses causing to move the flapper an additional amount in the direction of the original displacement e . This will cause the back pressure, p b , in the nozzle to change continuously. In other words in addition to the direct action of the bellows I, a positive feedback effect emerges as the time goes on. Using the methodology employed in the previous section, it can be deduced that the obtained performance is inherently an integral action. However, because of the resulted fractional order behavior of the bellows and the diaphragm of the pneumatic relay, one can draw a block diagram as shown in Figure 7. The effect of positive feedback can be observed easily from the figure.

A simplified block diagram for fractional order PI controller.

In this block diagram we denote the fractional order model for the restrictor-bellows II with orders w . Note that the fractional order integral control action in the controller takes the form of slowly canceling the feedback that the proportional control originally provided.

Using Mason’s formula, one can develop the transfer function for the block diagram sketched in Figure 7: (20) P c ( s ) E ( s ) = ( b a + b k 1 k 2 ( 1 + ( τ new s ) β Ψ min + Ψ max ( τ new s ) β ) ) × ( 1 + k 1 k 2 ( 1 + ( τ new s ) β Ψ min + Ψ max ( τ new s ) β ) × ( a a + b ) ( A k s s q ) ( 1 - 1 R C s w + 1 ) ( 1 + ( τ new s ) β Ψ min + Ψ max ( τ new s ) β ) ) - 1 = ( b k 1 k 2 k s s q ( R C s w + 1 ) ( 1 + ( τ new s ) β ) ) × ( ( a + b ) k s s q ( R C s w + 1 ) ( Ψ min + Ψ max ( τ new s ) β ) + a A k 1 k 2 R C s w ( 1 + ( τ new s ) β ) ) - 1 .

Notice that if | k 1 k 2 ( ( 1 + ( τ new s ) β ) / ( Ψ min + Ψ max ( τ new s ) β ) ) ( a / ( a + b ) ) ( A / k s s q ) ( R C s w / ( R C s w + 1 ) ) | 1 , (20) can be simplified as follows: (21) P c ( s ) E ( s ) ( b a + b k 1 k 2 ( 1 + ( τ new s ) β Ψ min + Ψ max ( τ new s ) β ) ) × ( k 1 k 2 ( 1 + ( τ new s ) β Ψ min + Ψ max ( τ new s ) β ) × ( a a + b ) ( A k s s q ) ( 1 - 1 R C s w + 1 ) ( 1 + ( τ new s ) β Ψ min + Ψ max ( τ new s ) β ) ) - 1 = b k s s q ( R C s w + 1 ) a A R C s w = M s q ( 1 + 1 R C s w ) in which M = b k s ( a A ) - 1 . Note that if w > q , the transfer function indicates an integral plus derivative action, which can be called a fractional I + D controller.

Remark 3.

Using a similar method, one can develop a fractional PID controller for a nozzle-flapper-relay configuration.

4. Simulation Results

In this section some simulations are given for the developed theory in the previous sections.

Example 4.

Consider the following values for a nozzle-flapper amplifier discussed in Section 3.1: b = 9 ,    a = 1 ,    k 1 = 0.1 ,    k 2 = 0.1 , A = 2 ,    k s = 20 ,    τ new = 0.001 ,    Ψ min = 1 ,    Ψ max = 10 ,    β = 0.8 , q = 0.9 , 0.99 , 1 .

The step responses are depicted in Figure 8. As can be seen for larger β ’s the settling time will be better.

Step response of the nozzle-amplifier-relay (17) for various betas.

Remark 5.

The values of k 1 and k 2 have been chosen from the procedure and graphs given in . The tuning technique presented there is a rough approach to us for choosing the values of the pneumatic structures presented in this paper. We have used the algorithm and graphs proportionally to other parameters in the nozzle-flapper presented in our configurations. As the optimality of such values in the integer-order counterparts is a difficult problem (and in general is an open problem), its status in the fractional order cases is consequently more difficult and up to now there is no trend in choosing the parameters optimally.

Example 6.

Now consider the nozzle-flapper-relay studied in Section 3.2. The parameters are similar to those of Example 4, which additionally v = 0.94 and R C = 0.1 . Note that because of inherent hysteresis in the practical diaphragm and bellows, a delay is observed in the step responses which present a more accurate response than classical models using linear spring. Thus, considering fractional order model for the diaphragms and bellows is justified and is more consistent with practice.

Moreover the step changes in the error signal yield a differentiation performance for the apparatus designed in Figure 9.

Step response of the fractional order derivative nozzle-flapper-relay for (20).

Example 7.

As the final example consider the nozzle-flapper-relay studied in Section 3.3. The parameters are similar to those of Example 4, which additionally w = 0.94 and R C = 0.1 and k 1 = 2 , k 2 = 5 , and q = 0 . As can be seen from the figure if we set q = 0 , the configuration shown in Figure 6 is a proportional plus fractional integral controller. The step responses are proposed in Figure 10.

Step response of the fractional order derivative nozzle-flapper-relay for (21).

Remark 8.

We have employed the Simulink package to simulate these numerical examples. All standard blocks have been taken from its library. However, the fractional order derivative blocks have been generated via Oustaloup approximation frequency technique which is a routine and reliable method in simulating the fractional order control systems .

5. Conclusions

In this paper in the light of fractional calculus, some new models have been developed for three configurations of nozzle-flapper-relay amplifier. Indeed using a typical form of the four-parameter fractional derivative Zener model for viscoelastic materials, a fractional order model has been provided for the diaphragm of pneumatic relays and the spring property of the bellows. After developing block diagrams for each configuration, a fractional order derivative and integral controller with two freedoms in orders have been obtained. The fractional orders that are imposed in the models can be used as other adjustable parameters of the controllers which inherently depend upon the material of the bellows and diaphragms. Some numerical simulations are proposed to clarify the obtained results.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Cafagna D. Past and present—fractional calculus: a mathematical tool from the past for present engineers IEEE Industrial Electronics Magazine 2007 1 2 35 40 2-s2.0-36348976811 10.1109/MIE.2007.901479 Ogata K. Modern Control Engineering 2010 Prentice Hall Arkhincheev V. E. Anomalous diffusion in inhomogeneous media: some exact results Modelling, Measurement and Control A 1993 26 11 29 Krishna B. T. Studies on fractional order differentiators and integrators: a survey Signal Processing 2011 91 3 386 426 2-s2.0-78049324242 10.1016/j.sigpro.2010.06.022 ZBL1203.94035 Djordjevic V. D. Atanackovic T. M. Similarity solutions to nonlinear heat conduction and Burgers/Korteweg-deVries fractional equations Journal of Computational and Applied Mathematics 2008 222 2 701 714 2-s2.0-53549087275 10.1016/j.cam.2007.12.013 ZBL1157.35470 Cole K. S. Electric Conductance of Biological Systems 1993 New York, NY, USA Cold Spring Harbor Glockle W. G. Nonnenmacher T. F. A fractional calculus approach to self-similar protein dynamics Biophysical Journal 1995 68 1 46 53 2-s2.0-0028878140 10.1016/S0006-3495(95)80157-8 Laskin N. Fractional market dynamics Physica A 2000 287 3-4 482 492 MR1802416 10.1016/S0378-4371(00)00387-3 Ingman D. Suzdalnitsky J. Application of differential operator with servo-order function in model of viscoelastic deformation process Journal of Engineering Mechanics 2005 131 7 763 767 2-s2.0-23044492109 10.1061/(ASCE)0733-9399(2005)131:7(763) Debnath L. Recent applications of fractional calculus to science and engineering International Journal of Mathematics and Mathematical Sciences 2003 2003 54 3413 3442 MR2025566 10.1155/S0161171203301486 ZBL1036.26004 Hamamci S. E. Stabilization using fractional-order PI and PID controllers Nonlinear Dynamics 2008 51 1-2 329 343 2-s2.0-35448968971 10.1007/s11071-007-9214-5 ZBL1170.93023 Jun Yang X. Local Fractional Functional Analysis and Its Applications 2011 HongKong, China Asian Academic Publisher Limited Andersen B. W. The Analysis Design of Pneumatic Systems 1967 John Wiley & Sons Miyajima T. Fujita T. Sakaki K. Kawashima K. Kagawa T. Development of a digital control system for high-performance pneumatic servo valve Precision Engineering 2007 31 2 156 161 2-s2.0-33847153586 10.1016/j.precisioneng.2006.05.003 Nishioka Y. Suzumori K. Kanda T. Wakimoto S. Multiplex pneumatic control method for multi-drive system Sensors and Actuators A 2010 164 1-2 88 94 2-s2.0-78249263207 10.1016/j.sna.2010.09.004 Moon J.-H. Lee B.-G. Modeling and sensitivity analysis of a pneumatic vibration isolation system with two air chambers Mechanism and Machine Theory 2010 45 12 1828 1850 2-s2.0-77957259374 10.1016/j.mechmachtheory.2010.08.006 ZBL05829319 Atanackovic T. M. Pilipovic S. Zorica D. Forced oscillations of a body attached to a viscoelastic rod of fractional derivative type International Journal of Engineering Science 2013 64 54 65 MR3023956 10.1016/j.ijengsci.2012.12.006 Xu Z. Chen W. A fractional-order model on new experiments of linear viscoelastic creep of Hami Melon Computers & Mathematics with Applications 2013 66 677 681 MR3089376 10.1016/j.camwa.2013.01.033 Schiessel H. Metzler R. Blumen A. Nonnenmacher T. F. Generalized viscoelastic models: their fractional equations with solutions Journal of Physics A 1995 28 23 6567 6584 2-s2.0-0000902727 10.1088/0305-4470/28/23/012 ZBL0921.73154 Bentil S. A. Dupaix R. B. Exploring the mechanical behavior of degrading swine neural tissue at low strain rates via the fractional Zener constitutive model Journal of the Mechanical Behavior of Biomedical Materials 2014 30 83 90 10.1016/j.jmbbm.2013.10.020 Holm S. Näsholm S. P. Prieur F. Sinkus R. Deriving fractional acoustic wave equations from mechanical and thermal constitutive equations Computers & Mathematics with Applications 2013 66 621 629 MR3089370 10.1016/j.camwa.2013.02.024 Saedpanah F. Well-posedness of an integro-differential equation with positive type kernels modeling fractional order viscoelasticity European Journal of Mechanics A 2014 44 201 211 MR3145383 10.1016/j.euromechsol.2013.10.014 Podlubny I. Fractional Differential Equations 1999 New York, NY, USA Academic Oldham K. B. Spanier J. The Fractional Calculus, Mathematics in Science and Engineering 1974 Academic press Kilbas A. A. Srivastava H. M. Trujillo J. J. Theory and Applications of Fractional Differential Equations 2006 Elsevier Pritz T. Analysis of four-parameter fractional derivative model of real solid materials Journal of Sound and Vibration 1996 195 1 103 115 2-s2.0-0030575463 10.1006/jsvi.1996.0406 ZBL1235.34026 http://www.gearseds.com/documents.html Moon J.-H. Lee B.-G. Modeling and sensitivity analysis of a pneumatic vibration isolation system with two air chambers Mechanism and Machine Theory 2010 45 12 1828 1850 2-s2.0-77957259374 10.1016/j.mechmachtheory.2010.08.006 ZBL05829319 Oustaloup A. Levron F. Mathieu B. Nanot F. M. Frequency-band complex noninteger differentiator: characterization and synthesis IEEE Transactions on Circuits and Systems I 2000 47 1 25 39 2-s2.0-0033888365 10.1109/81.817385