About this Journal Submit a Manuscript Table of Contents
International Journal of Differential Equations
VolumeΒ 2010Β (2010), Article IDΒ 846107, 16 pages
http://dx.doi.org/10.1155/2010/846107
Research Article

On the Selection and Meaning of Variable Order Operators for Dynamic Modeling

School of Engineering, University of California, P.O. Box 2039, Merced, CA 95344, USA

Received 4 August 2009; Accepted 8 October 2009

Academic Editor: NikolaiΒ Leonenko

Copyright Β© 2010 Lynnette E. S. Ramirez and Carlos F. M. Coimbra. 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.

Abstract

We review the application of differential operators of noninteger order to the modeling of dynamic systems. We compare all the definitions of Variable Order (VO) operators recently proposed in literature and select the VO operator that has the desirable property of continuous transition between integer and non-integer order derivatives. We use the selected VO operator to connect the meaning of functional order to the dynamic properties of a viscoelastic oscillator. We conclude that the order of differentiation of a single VO operator that represents the dynamics of a viscoelastic oscillator in stationary motion is a normalized phase shift. The normalization constant is found by taking the difference between the order of the inertial term (2) and the order of the spring term (0) and dividing this difference by the angular phase shift between acceleration and position in radians (πœ‹), so that the normalization constant is simply 2/πœ‹.

1. Introduction

The integer order differential operators of classical calculus (such as the first or second order derivatives) are familiar to anyone who has an active interest in understanding dynamic systems. These differential operators are used to formulate models that accurately describe the majority of physical phenomena and are ubiquitous in the mathematical description of dynamic behavior. However effective these integer order differential operators are in general, there are more complex systems that are better characterized by dynamic behavior that lies in between the normal integer order description. A case in point is the so-called β€œviscoelastic” behavior, which has characteristics of both elastic (order zero) and viscous (order one) elements. It is thus natural to assume that differential operators of noninteger order, such as a 0.25, 0.50, or 0.75 would provide a convenient mathematical description to analyze these intermediate behaviors. The study of these noninteger differential operators falls under the general subject of what became known as Fractional Calculus, though the orders studied are not strictly limited to rational numbers.

A further generalization of the concept of noninteger order derivatives that is applicable to more complex systems is that of a derivative of varying order. One can find systems where the order of dynamics associated with each element is a function of time or frequency or position or any derivative of the position vector [1–3]. The objective of this work is to identify the most appropriate definition of a variable-order (VO) operator for modeling dynamic systems and to assign to the order of the derivative a physical meaning that will facilitate the understanding of its use in problems of vibration and control. First, we compare all VO operator definitions recently proposed in literature in order to select a definition that is better suited for modeling purposes. We then use the familiar example of viscoelastic harmonic oscillators to connect the order of the derivative to the dynamic properties of the oscillators.

Unlike ordinary derivatives, noninteger order derivatives are integrodifferential operators with either a power-law (in the case of fractional derivatives) or a variable-exponent (for VO derivatives) kernel. Thus, noninteger derivatives are nonlocal by definition and are ideally suited for modeling systems characterized by nonlocal (or memory-laden) behavior. Multiple definitions for a fractional derivative have been proposed, and there is no straightforward geometric or physical interpretation for the meaning of a fractional derivative, although a few interpretations have been proposed [4]. The subject of Fractional Calculus has developed rapidly, especially since in the last four decades, with quite a few recent books dedicated exclusively to the subject (see, e.g., [4–8]). Some applications that involve fractional derivatives include particle motion at small but finite Reynolds numbers [9–12], viscoelastic constitutive equations [13–16], transport dynamics with anomalous diffusion [17], and fractional order control systems [4, 18].

The concept of a VO operator is a much more recent development and is less widely-known. The order of the derivative/integral is allowed to vary over the domain of interest, thus they are suitable for modeling systems with evolving dynamics. Such systems include deformation of viscoelastic materials [19–22] and the mechanics of variable viscoelastic oscillators [1, 3]. Similar to the state of affairs in Fractional Calculus, multiple definitions of a VO derivative have been suggested [1, 2, 19, 23, 24], each preferred for different reasons. Since the kernel of the VO operators have a variable-exponent, analytical solutions to VO differential equations (VODEs) are more difficult to obtain, and have not been the focus of much attention. However, numerical solutions for the VODEs that have been formulated with the various VO operators have been developed [1–3]. In addition, rather than seeking a constant or multiple constants for the order of the derivative in the model that would accurately represent the data, in the VO case a function needs to be determined through mathematical and/or physical arguments. Previous applications of VO derivatives have either explicitly chosen the function [1–3] or found an approximation numerically [19, 20, 22]. This complication can be lessened if the functional form of the order is determined through physical arguments as in [21], where a model based on statistical mechanics was developed to describe the compression of a viscoelastic material.

The aim of this work is to compare the Variable Order operator definitions that have been proposed and to select the operator with the characteristics that are critical for the success of a dynamic model. We compare the various VO operators based on a very simple criteria: the VO operator must return the correct fractional derivative that corresponds to the argument of the functional order. In other words, if the argument π‘ž(𝑑) is equal to, say, 0, 0.50 or 1, then the function itself (zero-order derivative), the half-derivative, or the first derivative must be returned as the output of the operator. We will see that only two definitions previously proposed satisfy this elementary requirement. Of the two operator definitions that satisfy this property, one is more efficient from the numerical standpoint, and is therefore adopted in the remainder of this work. The appropriate operator then is used to study the somewhat familiar problem of a harmonically forced oscillator with viscoelastic damping. The goal of this second part of this work is to illustrate how a familiar problem in dynamics can be used to understand the meaning of a VO operator, and to understand how the dynamics in this familiar problem is affected by the physical parameters of the system using a VO analysis.

The next section presents an overview of the various VO operator definitions and a brief comparison of the VO operators applied to a harmonic and other bounded function. Subsequent to selecting the operator, we propose a VO model for the harmonically forced oscillator with viscoelastic damping of order 𝑝 (0<𝑝<1) and conduct a stationary analysis that yields a very concrete meaning to the order of the operator.

2. Variable Order Operators

The VO operator definitions that have been proposed are either direct extensions of the fractional calculus definitions or generalizations that arise from Laplace or Fourier transformations. In the direct extension approach, the constant exponent in the fractional operator is replaced with a function. For example, a VO integral is defined in [23] as

π’Ÿβˆ’π‘ž(𝑑)𝑐+1𝑓(𝑑)=Ξ“[]ξ€œπ‘ž(𝑑)𝑑𝑐(π‘‘βˆ’πœŽ)π‘ž(𝑑)βˆ’1𝑓(𝜎)d𝜎.(2.1) When π‘ž(𝑑)=𝛼= constant, then the 𝛼th-order fractional integral is recovered. Other definitions can be formulated by changing the form of the argument of the exponent to be π‘ž=π‘ž(𝜎) or π‘ž=π‘ž(π‘‘βˆ’πœŽ)and considering the Gamma function under the integral sign [24, 25]:

1D0βˆ’π‘ž(𝑑)ξ€œπ‘“(𝑑)=𝑑0(π‘‘βˆ’πœŽ)π‘ž(𝑑)βˆ’1Ξ“[]π‘ž(𝑑)𝑓(𝜎)d𝜎,(2.2)2D0βˆ’π‘ž(𝑑)ξ€œπ‘“(𝑑)=𝑑0(π‘‘βˆ’πœŽ)π‘ž(𝜎)βˆ’1Ξ“[]π‘ž(𝜎)𝑓(𝜎)d𝜎,(2.3)3D0βˆ’π‘ž(𝑑)ξ€œπ‘“(𝑑)=𝑑0(π‘‘βˆ’πœŽ)π‘ž(π‘‘βˆ’πœŽ)βˆ’1Ξ“[]π‘ž(π‘‘βˆ’πœŽ)𝑓(𝜎)d𝜎.(2.4) In the cases above, the lower terminal is set equal to 0, and it is assumed that 𝑓(0)=0 for 𝑑<0. Since (2.3) and (2.4) involve the variable of integration within the exponent, then this implies memory in the order, with the past states having a stronger effect on the order for definition (2.4) [24]. Also, the full convolution form of (2.4) enables use of the convolution properties to study the operator.

Similarly, a VO derivative definition can be obtained by directly substituting π‘ž=π‘ž(𝑑) in the Riemann-Liouville fractional derivative definition [23] valid for 0<π‘ž<1:

Dπ‘ž(𝑑)𝑐+1𝑓(𝑑)=Ξ“[]d1βˆ’π‘ž(𝑑)ξ€œd𝑑𝑑𝑐𝑓(𝜎)(π‘‘βˆ’πœŽ)π‘ž(𝑑)d𝜎.(2.5) Further definitions are obtained by taking π‘š first-order derivatives of the integrals defined in (2.2)–(2.4)

D0π‘ž(𝑑)d𝑓(𝑑)=π‘šdπ‘‘π‘šξ€œπ‘‘0(π‘‘βˆ’πœŽ)π‘šβˆ’π‘ž(𝑑,𝜎)βˆ’1Ξ“[π‘ž](𝑑,𝜎)𝑓(𝜎)d𝜎,(2.6) where π‘šβˆ’1<π‘ž(𝑑)<π‘š. Caputo-type VO operators are defined by taking the derivative of the function under the integrals

D0π‘ž(𝑑)ξ€œπ‘“(𝑑)=𝑑0(π‘‘βˆ’πœŽ)π‘šβˆ’π‘ž(𝑑,𝜎)βˆ’1Ξ“[]π‘“π‘ž(𝑑,𝜎)(π‘š)(𝜎)d𝜎,(2.7) where 𝑓(π‘š) denotes the π‘šth integer order derivative of 𝑓(𝑑). VO operators based on other fractional derivative definition forms have also been proposed. Samko and Ross [23] introduce a VO operator based on the Marchaud fractional derivative:

Dπ‘π‘ž(𝑑)𝑓(𝑑)=𝑓(𝑑)Ξ“[]1βˆ’π‘ž(𝑑)(π‘‘βˆ’π‘)π‘ž(𝑑)+π‘ž(𝑑)Ξ“[]ξ€œ1βˆ’π‘ž(𝑑)𝑑𝑐𝑓(𝑑)βˆ’π‘“(𝜎)(π‘‘βˆ’πœŽ)1+π‘ž(𝑑)d𝜎,(2.8) where 0<Reπ‘ž(𝑑)<1.

Using a different approach, Coimbra [1] begins with the Laplace transform of the Caputo operator to obtain a VO operator definition. Treating π‘ž(𝑑) as a running parameter and inverting back into the time domain yields the following definition valid for π‘ž(𝑑)<1:

πƒπ‘ž(𝑑)1𝑓(𝑑)=Ξ“[]ξ€œ1βˆ’π‘ž(𝑑)𝑑0+(π‘‘βˆ’πœŽ)βˆ’π‘ž(𝑑)𝑓(1)(𝜎)d𝜎+(𝑓(0+)βˆ’π‘“(0βˆ’))π‘‘βˆ’π‘ž(𝑑)Ξ“[].1βˆ’π‘ž(𝑑)(2.9) Extension of definition (2.9) to values of π‘ž(𝑑) larger than unity is possible as long as the higher-order derivatives are defined and integrable. Although this definition began with the Laplace transform of the Caputo operator, it is not strictly a Caputo generalization because of the addition of the initial condition term. As a result, Dπ‘ž(0)𝑓(0)β‰ 0 for any function 𝑓(𝑑) as is the the case for a Caputo-based operator. Soon et al. [3] show that a properly weighted sum of fractional order derivatives terms approximates this single term VO operator when a large number of terms are used, which implies convergence of both the operator and the numerical method.

Although definitions (2.8) and (2.9) were defined independently through different methods, they are similar. After integrating (2.8) by parts and simplifying, we arrive at

1Ξ“[]ξ€œ1βˆ’π‘ž(𝑑)𝑑𝑐𝑓(1)(𝜎)(π‘‘βˆ’πœŽ)π‘ž(𝑑)𝑓d𝜎+(𝑐)Ξ“[]1βˆ’π‘ž(𝑑)(π‘‘βˆ’π‘)π‘ž(𝑑).(2.10) The difference between (2.9) and (2.10) lies in the term that is evaluated at the lower terminal. For situations in which the system is assumed to be in dynamical equilibrium for 𝑑<0 such that 𝐷2𝑓(𝑑<0)=0 for any 𝑓(𝑑), definitions (2.9) and (2.10) return the same result. However, if 𝑓(0) is a true constant, such that 𝑓(0βˆ’)=𝑓(0+)=π‘Ž, then (2.9) would return 0 for the derivative, whereas (2.10) will not.

In summary, there are a total of nine VO operator definitions to be compared:

1D0π‘ž(𝑑)1𝑓(𝑑)=Ξ“[]d1βˆ’π‘ž(𝑑)ξ€œd𝑑𝑑0𝑓(𝜎)(π‘‘βˆ’πœŽ)βˆ’π‘ž(𝑑)d𝜎,(2.11)2D0π‘ž(𝑑)𝑓(𝑑)=𝑓(𝑑)Ξ“[]1βˆ’π‘ž(𝑑)(𝑑)π‘ž(𝑑)+π‘ž(𝑑)Ξ“[]ξ€œ1βˆ’π‘ž(𝑑)𝑑0𝑓(𝑑)βˆ’π‘“(𝜎)(π‘‘βˆ’πœŽ)1+π‘ž(𝑑)d𝜎,(2.12)3βˆ’5D0π‘ž(𝑑)𝑓d(𝑑)=ξ‚΅ξ€œd𝑑𝑑0(π‘‘βˆ’πœŽ)π‘ž(𝑑,𝜎)βˆ’1Ξ“[]π‘“ξ‚Άπ‘ž(𝑑,𝜎)(𝜎)d𝜎,(2.13)6βˆ’8D0π‘ž(𝑑)ξ€œπ‘“(𝑑)=𝑑0(π‘‘βˆ’πœŽ)π‘ž(𝑑,𝜎)βˆ’1Ξ“[π‘ž]𝑓(𝑑,𝜎)(1)(𝜎)d𝜎,(2.14)9D0βˆ’π‘ž(𝑑)1𝑓(𝑑)=Ξ“[]ξ€œ1βˆ’π‘ž(𝑑)𝑑0+(π‘‘βˆ’πœŽ)βˆ’π‘ž(𝑑)𝑓(1)((𝜎)d𝜎+𝑓(0+)βˆ’π‘“(0βˆ’))π‘‘βˆ’π‘ž(𝑑)Ξ“[],1βˆ’π‘ž(𝑑)(2.15) where π‘ž(𝑑,𝜎) in definitions (2.13) and (2.14) signify the three arguments: π‘ž(𝑑,𝜎)=π‘ž(𝑑), π‘ž(𝑑,𝜎)=π‘ž(𝜎), and π‘ž(𝑑,𝜎)=π‘ž(π‘‘βˆ’πœŽ). Each of the above definitions is defined for real derivative orders between 0 and 1. The value at all the lower terminals is set to 0, since we are interested in applying the operators to physical processes not necessarily at steady state (we examine a stationary problem in the next section). As is the case with fractional derivatives, there is no single VO derivative (or integral) definition that is widely considered to be the β€œcorrect" definition. Samko and Ross prefer definition (2.12) because the operator retains the symmetry on power functions that is found in the case of constant orders, that is:

Dπ‘ž(𝑑)𝑐++(π‘‘βˆ’π‘)𝛼=Ξ“[]𝛼+1Ξ“[]𝛼+1βˆ’π‘ž(𝑑)(π‘‘βˆ’π‘)π›Όβˆ’π‘ž(𝑑)(2.16) for 0<Reπ‘ž(𝑑)<1 and 𝛼>βˆ’1 [6]. Lorenzo and Hartley prefer the full convolution VO integral definition (2.4) because it satisfies the index rule for certain functions [25] and is time-invariant [24].

The approach chosen here is to determine which operator when acting upon a function returns the fractional derivative of the function at the corresponding time. This is an important characteristic from the aspect of physical modeling because it signifies that the operator yields a continuous transition of all orders of differentiation between integer orders. Thus, a smooth transition from zero-order dynamics to first-order dynamics is possible. For a qualitative comparison of the the VO operators, we look at the π‘ž(𝑑)=𝑑 derivative of two bounded functions: sin(2πœ‹π‘‘), and erfc(𝑑). The π‘ž(𝑑)=𝑑 derivative of both functions is computed numerically using a product trapezoidal rule to evaluate the convolution integrals [3, 26]. Plots of the 𝑑 derivative of sin(2πœ‹π‘‘) and erfc(𝑑) are shown in Figures 1-2. The 0, 0.25, 0.50, 0.75, and 1st-order derivatives of the sinusoidal function are shown for comparison since both the Riemann-Liouville and Caputo fractional derivatives return the same result. We show only the 0- and 1st-order derivatives of the erfc function since the different fractional derivative definitions yield different results. Also, note that Samko's Marchaud-based definition and Coimbra's definition coincide for these conditions, and therefore only definition (2.15) is plotted.

fig1
Figure 1: Plots of the 𝑑 derivatives of sin(2πœ‹π‘‘). The points are the 𝛼=0,0.25,0.50,0.75, and 1st-order fractional derivatives at 𝑑=𝛼. (a) Definition (2.13) with π‘ž(𝑑,𝜎)=π‘ž(𝑑) (thin line), π‘ž(𝑑,𝜎)=π‘ž(𝜎) (medium line), and π‘ž(𝑑,𝜎)=π‘ž(π‘‘βˆ’πœŽ) (thick line). Note that none of the definitions match the 1st-order derivative at 𝑑=1. (b) Caputo-type VO operators (2.14) with π‘ž(𝑑,𝜎)=π‘ž(𝑑) (thin line), π‘ž(𝑑,𝜎)=π‘ž(𝜎) (medium line), and π‘ž(𝑑,𝜎)=π‘ž(π‘‘βˆ’πœŽ) (thick line). Note that the variant of (2.14) with argument 𝑑 matches the corresponding 𝛼 fractional derivatives at 𝑑=𝛼. (c) Definition (2.11) (thick line) and Coimbra's operator (2.15) (thin line). The 𝑑-derivative defined by (2.15) is equivalent to the corresponding fractional derivatives at all the points.
fig2
Figure 2: Plots of the 𝑑 derivatives of erfc(𝑑). The points are the 0- and 1st-order derivatives at 𝑑=0 and 𝑑=1. Also shown for comparison is the 𝑑 derivative from Coimbra's operator (2.15) (dashed line). (a) Riemann-Liouville type operator (2.13) with π‘ž(𝑑,𝜎)=π‘ž(𝑑) (thin line), π‘ž(𝜎) (medium line), π‘ž(π‘‘βˆ’πœŽ) (thick line). (b) Caputo-type definition (2.14) with π‘ž(𝑑,𝜎)=π‘ž(𝑑) (thin line), π‘ž(𝜎) (medium line), π‘ž(π‘‘βˆ’πœŽ) (thick line). In this case the operator (2.14) with exponent π‘ž(𝑑) does not match Coimbra's operator, but is equivalent to the first derivative at 𝑑=1. (c) Definition (2.11) (thin line).

As expected, the Caputo-based VO operators have values of 0 at 𝑑=0, similar to the fractional derivative case. They do not return the zeroth-order derivative at 𝑑=0 for any functions that deviate from this (see Figure 2). Definition (2.14) is equivalent to (2.15) for the sine function, since 𝑓(𝑑) is continuous at 𝑑=0 and π‘ž(𝑑) is only a function of 𝑑. The Riemann-Liouville and Caputo type operators with the argument (π‘‘βˆ’πœŽ), also are shown to be equivalent because 𝑓(0)=0. This is analogous to the similarity of the Riemann-Liouville and Caputo fractional derivatives of functions when 𝑓(0)(0)=0 for 0β‰€π‘ž<1 [4]. Definitions (2.12) and (2.15) are the only operators that have the desirable property of returning the corresponding π‘žπ‘‘β„Ž order fractional derivative of π‘₯(𝑑) when π‘ž(𝑑)=𝑝 for both the sine and erfc funcitons. However, the convergence of (2.12) to that of (2.15) is slower due to the stronger singularity that must be evaluated in the convolution integral. Also, in the case of a true constant function 𝑓(𝑑) where 𝑓(0βˆ’)=𝑓(0+), (2.12) would not return 0 as the derivative, so we conclude that definition (2.15) is preferable for modeling dynamic systems.

Note that the operator defined in (2.15) is dynamically consistent with the causal behavior of the initial conditions. In other words, when π‘₯(𝑑) is a true constant from βˆ’βˆž to the initial time (𝑑=0+), the operator in (2.15) returns zero for all values of π‘ž(𝑑). However, if 𝑓(𝑑) is not continuous between 𝑑=0βˆ’ and 𝑑=0+, the operator returns the appropriate Heaviside contribution to the integral value of π’Ÿπ‘ž(𝑑)𝑓(𝑑). In accordance with this causal definition, we take the value of the physical variable 𝑓(𝑑) to be identically null from βˆ’βˆž to 0βˆ’ as a representation of dynamic equilibrium. A nonzero initial condition is treated as a Heaviside function at 𝑑=0, and therefore included in the second term of the definition of the operator (2.15).

Through the direct comparison of the various proposed VO operator definitions, we selected the operator that has fundamental characteristics that are desirable for physical modeling. Definition (2.15) proposed in [1] represents a continuous transition between the integer order derivatives, and returns a zero value for the derivative of a function that is constant from βˆ’βˆž<𝑑<∞; so we select it as the most appropriate definition. Now with the chosen operator, we proceed to connect the behavior of a VO operator with a physical quantity that is characteristic of a selected memory-laden system.

3. Stationary Analysis for Viscoelastic Oscillators

One of the drawbacks of VO modeling is that to date there is no clear physical understanding of what a VO derivative represents. The objective of this section is to illustrate how a variable order differential operator may be used to understand a familiar problem: the stationary analysis of a constant order viscoelastic oscillator. We will use a single VO operator of order π‘ž to replace multiple terms of constant order differential operators, including the viscoelastic term of constant order 𝑝. We seek an analytical expression for π‘ž to examine the effects of the parameters of the system on the dynamics of the oscillator. The exact expression for π‘ž(πœ”) is obtained from the stationary analysis of the problem. In order words, we look for two functions π‘ž that would allows us to replace the multiterm differential equation describing the steady motion of the oscillator with a single-term variable order equation. Note that the objective of this section is not to rehash the analysis of the constant order viscoelastic oscillator, rather we seek to find meaning for a VO operator by comparing its order with the physical parameters in a well-known problem. The reader who is interested in the details of the dynamics of constant order viscoelastic oscillators should consult [4, 27, 28].

The equation of motion for the constant order viscoelastic oscillator is π‘šπ·2π‘₯(𝑑)+𝑐𝐷𝑝π‘₯(𝑑)+π‘˜π·0π‘₯(𝑑)=𝐹(𝑑),(3.1) where 𝐹(𝑑)=𝐹0cos(Ω𝑑). When 𝑝=1, the system is an oscillator with viscous damping, and for 0<𝑝<1, the system is said to be characterized by viscoelastic damping. The equation is recast in dimensionless form using the following scaled parameters:

π‘₯Μƒπ‘₯=𝐿𝑐,̃𝑑=π‘‘πœ”π‘›,ΩΩ=πœ”π‘›,(3.2) where 𝐿𝑐 is a characteristic length (or amplitude of the motion) and πœ”π‘›=βˆšπ‘˜/π‘š is the undamped natural frequency of the system. The dimensionless equation of motion is

𝐷2𝐷̃π‘₯+πœ‰π‘ξ‚πΉΜƒπ‘₯+Μƒπ‘₯=π‘œξ‚€ξ‚Ξ©Μƒπ‘‘ξ‚cos,(3.3) where 𝐹0=𝐹0/π‘šπœ”2𝑛𝐿𝑐 is the dimensionless amplitude of the forcing, and πœ‰ is the damping ratio and is

πœ‰=π‘π‘˜π‘/2βˆ’1π‘šπ‘/2.(3.4)

Since we are concerned only with stationary behavior, πœ‚ and π‘ž are not functions of time, so we look for a relationship such as πœ‚ξ‚€ξ‚ξ‚π·π‘,Ξ©,πœ‰π‘žξ‚΅π‘,Ω,πœ‰ξ€·π·Μƒπ‘₯=2+πœ‰π·π‘+𝐷0𝐹̃π‘₯=Re0ξ‚€iΩ̃𝑑.exp(3.5) The stationary VO derivative of 𝑒iΩ̃𝑑 is

π·π‘ž(𝑝,Ω,πœ‰)𝐴𝑒iΩ̃𝑑=ξ‚€iξ‚Ξ©ξ‚π‘ž(𝑝,Ω,πœ‰)𝐴𝑒iΩ̃𝑑,(3.6) where definition (2.15) with a lower terminal of βˆ’βˆž is used since we are dealing with a stationary problem where the initial conditions are irrelevant. Rewriting (3.6) with (3.3) as the stationary solution Ω̃exp(i𝑑) yields

πœ‚ξ‚€iξ‚Ξ©ξ‚π‘ž(𝑝,Ω,πœ‰)Ω=1βˆ’2+ξ‚€iξ‚Ξ©ξ‚π‘πœ‰.(3.7) We now equate the real and imaginary parts of the above equation to arrive at

π‘žξ‚€ξ‚ξ‚π‘,Ξ©,πœ‰=tanβˆ’1βŽ›βŽœβŽœβŽπœ‰ξ‚Ξ©π‘sin((πœ‹/2)𝑝)Ω1βˆ’2+ξ‚Ξ©π‘βŽžβŽŸβŽŸβŽ 2πœ‰cos((πœ‹/2)𝑝)πœ‹,πœ‚ξ‚€ξ‚ξ‚=1𝑝,Ξ©,πœ‰ξ‚Ξ©π‘žξ‚™ξ‚€ξ‚Ξ©1βˆ’22+ξ‚€ξ‚Ξ©π‘πœ‰ξ‚2Ω+21βˆ’2ξ‚ξ‚Ξ©π‘πœ‰cos(πœ‹/2)𝑝.(3.8) The stationary solution for (3.5) can be written as ̃Ω̃̃π‘₯(𝑑)=𝐴cos(π‘‘βˆ’πœ™) where the amplitude 𝐴 and phase shift πœ™ are

𝐹𝐴=0Ω1βˆ’22+ξ‚€ξ‚Ξ©π‘πœ‰ξ‚2Ω+21βˆ’2Ω𝑝,πœ‰cos(πœ‹/2)π‘πœ™=tanβˆ’1βŽ›βŽœβŽœβŽπœ‰ξ‚Ξ©π‘sin((πœ‹/2)𝑝)Ω1βˆ’2+ξ‚Ξ©π‘βŽžβŽŸβŽŸβŽ .πœ‰cos((πœ‹/2)𝑝)(3.9) The procedure used to obtain the functional form for π‘ž and πœ‚ is easily extendable to systems that consist of multiple viscoelastic terms. For 𝑛 terms, the expression for π‘ž is

π‘ž=tanβˆ’1βŽ›βŽœβŽœβŽβˆ‘π‘›π‘˜=1πœ‰π‘˜ξ‚Ξ©π‘π‘˜sin(πœ‹/2)π‘π‘˜ξ‚€ξ‚Ξ©1βˆ’2+βˆ‘π‘›π‘˜=1πœ‰π‘˜ξ‚Ξ©π‘π‘˜cos(πœ‹/2)π‘π‘˜βŽžβŽŸβŽŸβŽ 2πœ‹.(3.10) Similarly, the expression for πœ‚ is

1πœ‚=ξ‚Ξ©π‘žξ„Άξ„΅ξ„΅βŽ·ξƒ©ξ‚Ξ©1βˆ’2+π‘›ξ“π‘˜=1πœ‰π‘˜ξ‚Ξ©π‘π‘˜πœ‹cos2π‘π‘˜ξƒͺ2+ξƒ©π‘›ξ“π‘˜=1πœ‰π‘˜ξ‚Ξ©π‘π‘˜πœ‹sin2π‘π‘˜ξƒͺ2.(3.11)

Equations (3.8)-(3.9) clearly show that π‘ž represents a scaled phase shift and πœ‚ is the scaled ratio of the amplitude of the forcing to the response. Thus, for the stationary solution of the oscillator the order of the derivative is connected to a physical quantity that is characteristic of the system. The multiterm equation in frequency that represents the stationary motion of a viscoelastic operator can be replaced by a single term parametric operator in frequency, where both the order of the derivative and the scaling function πœ‚ have physical meaning. A phase shift of πœ‹/2 implies that the response is proportional to the velocity and hence to the 1st-order derivative, while a phase shift of πœ‹ implies that it is proportional to the acceleration or to the 2nd-order derivative. The order of the VO operator that captures the whole dynamics of the systems is thus naturally connected with the phase shift between the response and the forcing.

Plots of π‘ž and πœ‚ versus Ω for damping orders of 𝑝=0.25,0.50,0.75, and 1 and various values of the damping ratio are shown in Figures 3, 4, 5, and 6. Also shown are the maps of π‘ž and πœ‚.

fig3
Figure 3: Plots of π‘ž and πœ‚ for 𝑝=0.25. (a) ξ‚π‘ž(Ξ©,πœ‰) versus Ω for πœ‰=0.1,0.5,1,1.5, and 2 (going from light to dark). Note that π‘ž aysmptotes to 2 as Ω increases. (b) Map of ξ‚π‘ž(Ξ©,πœ‰) where lower values are represented by darker colors. When Ω=1, π‘ž=0.25 for any value of πœ‰. (c) ξ‚πœ‚(Ξ©,πœ‰) versus Ω for πœ‰=0.1,0.5,1,1.5, and 2 (going from light to dark). (d) Map of ξ‚πœ‚(Ξ©,πœ‰). The lowest values of πœ‚ occur near Ω=1 and for small damping ratio.
fig4
Figure 4: Plots of π‘ž and πœ‚ for 𝑝=0.50. (a) ξ‚π‘ž(Ξ©,πœ‰) versus Ω for πœ‰=0.1,0.5,1,1.5, and 2 (going from light to dark). The asymptotic behavior is similar to the case when 𝑝=0.25. (b) Map of ξ‚π‘ž(Ξ©,πœ‰) where lower values are represented by darker colors. When Ω=1, π‘ž=0.50 for any value of πœ‰. (c) ξ‚πœ‚(Ξ©,πœ‰) versus Ω for πœ‰=0.1,0.5,1,1.5, and 2 (going from light to dark). (d) Map of ξ‚πœ‚(Ξ©,πœ‰).
fig5
Figure 5: Plots of π‘ž and πœ‚ for 𝑝=0.75. (a) ξ‚π‘ž(Ξ©,πœ‰) versus Ω for πœ‰=0.1,0.5,1,1.5, and 2 (going from light to dark). (b) Map of ξ‚π‘ž(Ξ©,πœ‰) where lower values are represented by darker colors. When Ω=1, π‘ž=0.75 for any value of πœ‰. (c) ξ‚πœ‚(Ξ©,πœ‰) versus Ω for πœ‰=0.1,0.5,1,1.5, and 2 (going from light to dark). (d) Map of ξ‚πœ‚(Ξ©,πœ‰).
fig6
Figure 6: Plots of π‘ž and πœ‚ for 𝑝=1, signifying the oscillator with viscous damping. (a) ξ‚π‘ž(Ξ©,πœ‰) versus Ω for πœ‰=0.1,0.5,1,1.5, and 2 (going from light to dark). (b) Map of ξ‚π‘ž(Ξ©,πœ‰) where lower values are represented by darker colors. When Ω=1, π‘ž=0.75 for any value of πœ‰. (c) ξ‚πœ‚(Ξ©,πœ‰) versus Ω for πœ‰=0.1,0.5,1,1.5, and 2 (going from light to dark). (d) Map of ξ‚πœ‚(Ξ©,πœ‰).

The expression for π‘ž reveals regions in which the three terms of the original equation of motion are dominant. For example, regions in which the order of the derivative is near 2 (such as systems with lower damping ratio and higher forcing frequency) are primarily dominated by the inertial term. Lower values of damping ratio and orders of viscoelastic damping reach the asymptotic value more quickly, suggesting that the change in the dynamics (and also the phase shift) is more sensitive to changes in frequency in those cases. The dependence of the order of the derivative on the damping ratio also changes when Ω=1, corresponding to the case where the driving frequency is the same as the natural frequency of the system. When Ω<1, then the systems with higher damping ratios have a higher value for π‘ž. Once Ω>1 then the behavior is switched with the systems with higher damping ratio having lower values for the order. For all cases, when Ω=1, then π‘ž=𝑝 for any damping ratio.

A scaled behavior identical to π‘ž is shown in the plots of πœ‚. The normalized amplitude ratio reaches an asymptotic value of πœ‚=1 as Ω increases. Thus, at higher frequencies the amplitude ratio is proportional to ξ‚Ξ©βˆ’π‘, and the order of the damping and the damping ratio do not have any effect on the amplitude of the stationary motion. Similar to the case for π‘ž, πœ‚ approaches the asymptotic value more quickly for small viscoelastic damping orders 𝑝 and lower damping ratios. Also for the cases with damping order 𝑝<1, the peak amplitude response shifts to higher frequencies for increasing damping ratio. For increasing 𝑝 and damping ratio, the peak response begins to flatten out. The maps of πœ‚ show that as the order of the viscoelastic damping increases, the regions where the amplitude of the response is damped increase.

4. Conclusions

This work advances our understanding of the use of variable order (VO) differential operators in dynamics in two substantial ways. First, we compare several definitions proposed in literature, and select the most suitable definition based on a few criteria: (1) the VO operator must be able to return all intermediate values between 0 and 1 that correspond to the argument of the order of differentiation, (2) the VO operator must be effectively evaluated numerically, and (3) all derivatives of a true constant (a function that is constant from βˆ’βˆž to +∞) must be zero. The operator defined in (2.15) satisfies these criteria for modeling dynamic systems. We then proceed to illustrate the meaning of a variable order of differentiation by analyzing a familiar dynamical problem (the stationary analysis of a viscoelastic oscillator). We determine that the order of differentiation for a single operator describing all dynamic elements in the stationary equation of motion (mass, damping and spring) is equal to the normalized phase shift. The normalization is easily understood as the quantity that transforms the maximum phase shift between acceleration and position (πœ‹) into the maximum difference between the order of acceleration (2) and the order of position (0). The normalization constant is thus 2/πœ‹, and the variable order differentiation is seen as being just the normalized phase shift in this problem, which gives us a straightforward interpretation of the meaning of variable orders of differentiation in dynamic systems.

References

  1. C. F. M. Coimbra, β€œMechanics with variable-order differential operators,” Annalen der Physik, vol. 12, no. 11-12, pp. 692–703, 2003. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  2. D. Ingman and J. Suzdalnitsky, β€œControl of damping oscillations by fractional differential operator with time-dependent order,” Computer Methods in Applied Mechanics and Engineering, vol. 193, no. 52, pp. 5585–5595, 2004. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  3. C. M. Soon, C. F. M. Coimbra, and M. H. Kobayashi, β€œThe variable viscoelasticity oscillator,” Annalen der Physik, vol. 14, no. 6, pp. 378–389, 2005. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  4. I. Podlubny, Fractional Differential Equations, vol. 198 of Mathematics in Science and Engineering, Academic Press, San Diego, Calif, USA, 1999. View at MathSciNet
  5. K. B. Oldham and J. Spanier, The Fractional Calculus: Theory and Applications of Differentiation and Integration to Arbitrary Order, vol. 11, Academic Press, London, UK, 1974. View at MathSciNet
  6. S. G. Samko, A. A. Kilbas, and O. I. Marichev, Integraly i proizvodnye drobnogo poryadka i nekotorye ikh prilozheniya, “Nauka i Tekhnika”, Minsk, Belarus, 1987. View at MathSciNet
  7. K. S. Miller and B. Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, John Wiley & Sons, New York, NY, USA, 1993. View at MathSciNet
  8. A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Theory and Applications of Fractional Differential Equations, vol. 204 of North-Holland Mathematics Studies, Elsevier, Amsterdam, The Netherlands, 2006. View at MathSciNet
  9. C. F. M. Coimbra and R. H. Rangel, β€œGeneral solution of the particle momentum equation in unsteady Stokes flows,” Journal of Fluid Mechanics, vol. 370, pp. 53–72, 1998. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  10. C. F. M. Coimbra and R. H. Rangel, β€œSpherical particle motion in harmonic Stokes flows,” AIAA Journal, vol. 39, no. 9, pp. 1673–1682, 2001. View at Publisher Β· View at Google Scholar
  11. C. F. M. Coimbra and M. H. Kobayashi, β€œOn the viscous motion of a small particle in a rotating cylinder,” Journal of Fluid Mechanics, vol. 469, pp. 257–286, 2002. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  12. L. E. S. Ramirez, E. A. Lim, C. F. M. Coimbra, and M. H. Kobayashi, β€œOn the dynamics of a spherical scaffold in rotating bioreactors,” Biotechnology and Bioengineering, vol. 84, no. 3, pp. 382–389, 2003. View at Publisher Β· View at Google Scholar Β· View at PubMed
  13. R. L. Bagley and P. J. Torvik, β€œA theoretical basis for the application of fractional calculus to viscoelasticity,” Journal of Rheology, vol. 27, no. 3, pp. 201–210, 1983. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  14. R. L. Bagley and P. J. Torvik, β€œOn the fractional calculus model of viscoelastic behavior,” Journal of Rheology, vol. 30, no. 1, pp. 133–155, 1986. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  15. N. Heymans and J.-C. Bauwens, β€œFractal rheological models and fractional differential equations for viscoelastic behavior,” Rheologica Acta, vol. 33, no. 3, pp. 210–219, 1994. View at Publisher Β· View at Google Scholar
  16. T. Pritz, β€œFive-parameter fractional derivative model for polymeric damping materials,” Journal of Sound and Vibration, vol. 265, no. 5, pp. 935–952, 2003. View at Publisher Β· View at Google Scholar
  17. R. Metzler and J. Klafter, β€œThe random walk's guide to anomalous diffusion: a fractional dynamics approach,” Physics Reports, vol. 339, no. 1, pp. 1–77, 2000. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  18. I. Podlubny, β€œFractional-order systems and PIλDμ-controllers,” IEEE Transactions on Automatic Control, vol. 44, no. 1, pp. 208–214, 1999. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  19. D. Ingman, J. Suzdalnitsky, and M. Zeifman, β€œConstitutive dynamic-order model for nonlinear contact phenomena,” Journal of Applied Mechanics, vol. 67, no. 2, pp. 383–390, 2000. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  20. D. Ingman and J. Suzdalnitsky, β€œApplication of differential operator with servo-order function in model of viscoelastic deformation process,” Journal of Engineering Mechanics, vol. 131, no. 7, pp. 763–767, 2005. View at Publisher Β· View at Google Scholar
  21. L. E. S. Ramirez and C. F. M. Coimbra, β€œA variable order constitutive relation for viscoelasticity,” Annalen der Physik, vol. 16, no. 7-8, pp. 543–552, 2007. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  22. D. Ingman and J. Suzdalnitsky, β€œResponse of viscoelastic plate to impact,” Journal of Vibration and Acoustics, vol. 130, no. 1, Article ID 011010, 8 pages, 2008. View at Publisher Β· View at Google Scholar
  23. S. G. Samko and B. Ross, β€œIntegration and differentiation to a variable fractional order,” Integral Transforms and Special Functions, vol. 1, no. 4, pp. 277–300, 1993. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  24. C. F. Lorenzo and T. T. Hartley, β€œVariable order and distributed order fractional operators,” Nonlinear Dynamics, vol. 29, no. 1–4, pp. 57–98, 2002. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  25. C. F. Lorenzo and T. T. Hartley, β€œInitialization, conceptualization, and application in the generalized fractional calculus,” NASA Technical Publication 98-208415, Lewis Research Center, NASA, Cleveland, Ohio, USA, 1998.
  26. K. Diethelm, N. J. Ford, A. D. Freed, and Y. Luchko, β€œAlgorithms for the fractional calculus: a selection of numerical methods,” Computer Methods in Applied Mechanics and Engineering, vol. 194, no. 6–8, pp. 743–773, 2005. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet
  27. Y. A. Rossikhin and M. V. Shitikova, β€œApplications of fractional calculus to dynamic problems of linear and nonlinear hereditary mechanics of solids,” Applied Mechanics Reviews, vol. 50, no. 1, pp. 15–67, 1997.
  28. Y. A. Rossikhin and M. V. Shitikova, β€œApplication of fractional derivatives to the analysis of damped vibrations of viscoelastic single mass systems,” Acta Mechanica, vol. 120, no. 1–4, pp. 109–125, 1997. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH Β· View at MathSciNet