Stochastic Systems: Modeling, Analysis, Synthesis, Control, and their Applications to EngineeringView this Special Issue
Research Article | Open Access
Claudio Floris, "Stochastic Stability of Damped Mathieu Oscillator Parametrically Excited by a Gaussian Noise", Mathematical Problems in Engineering, vol. 2012, Article ID 375913, 18 pages, 2012. https://doi.org/10.1155/2012/375913
Stochastic Stability of Damped Mathieu Oscillator Parametrically Excited by a Gaussian Noise
This paper analyzes the stochastic stability of a damped Mathieu oscillator subjected to a parametric excitation of the form of a stationary Gaussian process, which may be both white and coloured. By applying deterministic and stochastic averaging, two Itô’s differential equations are retrieved. Reference is made to stochastic stability in moments. The differential equations ruling the response statistical moment evolution are written by means of Itô’s differential rule. A necessary and sufficient condition of stability in the moments of order r is that the matrix of the coefficients of the ODE system ruling them has negative real eigenvalues and complex eigenvalues with negative real parts. Because of the linearity of the system the stability of the first two moments is the strongest condition of stability. In the case of the first moments (averages), critical values of the parameters are expressed analytically, while for the second moments the search for the critical values is made numerically. Some graphs are presented for representative cases.
The so-called Mathieu equation has attracted the attention of the scholars in the past, and many papers and books are available on this subject such as [1, 2]. The interest in this equation stems from two facts: first, its dynamics is very rich; second, it describes the vibrations of important mechanical systems such as a prismatic bar stretched or compressed by a sinusoidal axial force, a pendulum whose support is subjected to a sinusoidal motion, and an elliptic membrane [3, 4]. In general, attention is focused on the stability of the motion.
When the excitation acting on the Mathieu oscillator is a stochastic process, its response is a stochastic process too. Thus, the analyst is concerned with the problem of the statistical characterization of the response, and in the case of a system prone to instability the problem of the stochastic stability must be solved. Nevertheless, studies on the Mathieu systems with stochastic excitation are not numerous. Probably, the first authors who addressed the problem were Stratonovich and Romanovski . Analyses of the Mathieu-type equations with external random excitation were presented by Dimentberg  and by Cai and Lin .
A Mathieu system with random parametrical excitation was considered by Rong et al.  and by Xie  (see the references cited herein too). The first authors use the method of multiple scale to determine the equations of the amplitude and the phase of the response. The check for almost-sure stability is made by computing the largest Lyapunov exponent. On the contrary, in ([9, Chapter 7]) the method of stochastic averaging is used, and the stability in moments is analyzed. It recalled that there are different definitions of stochastic stability, for which the reader is referred to [10, 11]. According to the selected stability definition, different stability bounds are found. Thus, it is not surprising that the results of [8, 9] are different. However, the conclusions that are reached are opposite: in  it is claimed that the stochastic parametric excitation is always detrimental with respect to the stability, while in  it is affirmed that there are cases in which the stochastic excitation stabilizes the system.
In this paper, the Mathieu damped oscillator parametrically excited by a Gaussian noise is considered again. A suitable coordinate transformation is made, and deterministic and stochastic averaging are applied in such a way that the second-order motion equation is replaced by two first-order equations, are interpreted in Itô’s sense. As the Fokker-Planck-Kolmogorov equation for the response probability density function (PDF) is not analytically solvable, the stochastic stability in moments is considered as in [9, 12]: the differential equations governing the response moment evolution are written by means of the rules of Itô’s stochastic differential calculus [13–16]. A necessary and sufficient condition of stability in the moments of order is that the matrix of the coefficients of the ODE system governing them has negative real eigenvalues and complex eigenvalues with negative real parts. Because of the linearity of the system the stability of the first two moments is the strongest condition of stability so that only the first two moments are considered ([11, Chapter 7]). Analytical stability bounds are established for the first moments, while the stability of the second moments is searched numerically.
2. Statement of the Problem
Consider the Mathieu-type stochastic differential equation where is the ratio of critical damping, is the undamped pulsation of the oscillator, is a small parameter, β and Ω are the amplitude and the pulsation of the sinusoidal term, respectively, and F(t) is a stochastic stationary Gaussian process with zero mean. Because of the parametric excitation F(t) the response X(t) is a stochastic process too.
The problem of the statistical characterization of (2.3) will be tackled by using the stochastic averaging method [6, 11, 17–20]. However, in the case of a harmonic term like that in (2.3), the classic coordinate transformation , where and are stochastic processes, does not work well as the joint probability density function (PDF) of the response is not separable ([11, chapter 7]). A suitable coordinate transformation was proposed in [6, 17, 18], which reads as According to the principles of deterministic averaging for weakly nonlinear systems  it is assumed that
The expressions in (2.5) are not exact but only approximate. The exact expression of the first derivative is . According to the method of deterministic averaging  the second and fourth terms are neglected so that the derivative retains the same form as it would have if and were constant.
The method of stochastic averaging is applied to (2.6), (2.7). In the field of stochastic dynamics this method of analysis was proposed first by Stratonovich [17, 18] (see also [11, 19, 20]) as an extension to stochastic systems of the deterministic method by Bogoliubov and Mitropolsky  (see Mettler , too). Then, it was rigorously demonstrated by Khasminskii . It involves two phases that can be performed in either order: in the former, the deterministic terms that do not contain the forcing functions are averaged. In the second phase, the terms containing the forcing functions are reduced to Gaussian stationary white noises.
Now, the forcing terms must be worked out. In the classic stochastic averaging method it is required that the process F(t) is broadbanded: in order to remove this restriction, a different way will be followed. Stratonovich [17, Volume 1, Chapter 7] suggested that, given a stochastic process F(t), which may be even narrowbanded, the product is replaceable by a stationary Gaussian white noise (t) having the autocorrelation function . The constant is given by the following integral: , where and are the auto-correlation function and the power spectral density (PSD) of , respectively. However, the forcing terms in (2.8) contain the contributions too. These contributions are replaced by a stationary Gaussian white noise (t), whose intensity is computed by adapting to the present case the derivation of Stratonovich [17, Volume 1, Section 7.2]. We have By expanding the product of the trigonometric functions in the integral (2.9), we obtain: . In the last expression, the second and the third addenda, when averaged, give rise to zero, while the fifth term is a faster oscillatory one that according to Stratonovich can be neglected. Thus, the integral (2.9) is equal to The expansion of leads to the same result. In much the same way, it results that , that is, they are uncorrelated.
Performing these operations, (2.8) simplify into For the sake of simplicity, the white noises in the previous equations have unit intensities. More-over, the forcing terms in the second equation have the sign plus as the white noise processes have the same probabilistic characteristics ([24, chapter 6, section 6.3.2]).
In order to transform (2.11) into two Itô-type stochastic differential equations, the so-called Wong-Zakai-Stratonovich corrective terms must be added to the drift terms [25, 26]. They are computed according to the following formula: where in the present case n = m = 2 and , are elements of the diffusion matrix Hence, the two Itô stochastic differential equations that govern the problem are where and (t) and (t) are two standard Wiener processes (Brownian motion), for which the formal relationship holds.
From a theoretical point of view the probabilistic characterization of the random vector X = should be obtained by solving the so-called Fokker-Planck-Kolmogorov (FPK) equation in the joint PDF of the state variables [11, 17, 18, 24] where the summation rule with respect to a repeated index has been adopted, and ,, with the matrix being defined in (2.13), while the apex “t” denotes transpose. Unfortunately, as in many other cases the FPK equation (2.15) does not have an analytical solution because the excitation is a parametric one. Thus, another method has to be chosen in order to characterize the probabilistic characteristics of the response. The moment equation approach is used herein.
3. Moment Equation Approach
In (2.14) the excitation is multiplicative only, and there is not an external excitation. Thus, zero solution = = 0 for all t satisfies them, even if the multiplicative excitation is a stochastic process. We are concerned with the stability of zero solution. There exist different definitions of stochastic stability, for which the reader is referred to Chap. 6 of . Here, the stability of zero solution is analyzed in the moments. In fact, the stability of the first two moments is the strongest for linear autonomous systems under multiplicative Gaussian excitation [11, 27]. Thus, the stability of the first two moments only will be considered here. As the response of (2.14) is not Gaussian, in order to characterize it statistically, from a theoretical point of view the knowledge of the infinite hierarchy of the moments would be necessary. Nevertheless, the system of (2.14) is linearly parametric: it has been shown that the equations for the statistical moments of such a type of systems are a close set , and they can be solved in succession. Thus, the convergence of the moments of first and second order is a necessary condition for being stable the moments of higher orders.
In the field of dynamic stability analyses, a system is stable when it comes back to the initial configuration after being subjected to a small perturbation. Zero statistical moments correspond to a zero solution. When the zero solution is perturbed, the response moments are no longer zero: if the stochastic Mathieu oscillator is stable, once the perturbation is removed, they decay to zero, otherwise they grow without limits. The first step is to write the ordinary differential equations (ODEs) ruling the time evolution of the response moments. Since (2.14) is an Itô system, use is made of Itô’s differential rule [13–15], which reads as In (3.1) Ψ is a nonanticipating function of the state variables and the summation rule is adopted. It is recalled that (3.1) retains the terms of order dt only, as dB is of order. In order to write the moment equations, the appropriate nonanticipating function to be chosen is , where i, j are positive integers or zero with the constraint , with being the order of the moments to be computed.
By choosing with i + j = 1, the ODEs for the first order moments are where the dots mean derivatives with respect to time and . Analogously, the ODEs ruling the second moments are found to be
By inspecting (3.3), it is noted that the forcing terms are absent, which due to the fact that in (2.14) the excitation is purely parametric. Thus, (3.3) can be written in compact matrix form as where (t) is a vector collecting all the moments of order of the system states and is a matrix of the coefficients as shown in (3.3). The solution to (3.4) is where is a vector whose entries are the initial conditions for the moments. These constitute the perturbation to zero solution.
It is well known (e.g., see ) that, as grows, the response moments decay to zero whenever the matrix has negative real eigenvalues and complex eigenvalues with negative real parts. Otherwise, the moments increase without limits. Thus, the condition of stability in moments is that the eigenvalues of the matrix have negative real parts. In this way, the stochastic problem is led to the classic deterministic problem of studying the eigenvalues of a matrix. Since the matrix depends on the system parameters (2.1), there exist critical values of these quantities for which the real part of almost an eigenvalue is zero. Increasing them further, a real part becomes positive, and the moments grow to infinity.
In order to study the eigenvalues of , its characteristic equation is formed: where is a unit matrix having the same order as . Equation (3.6) must specialized for the order of the moments. As previously advanced, due to the linearity of the system and to the Gaussianity of the input, the moment stability analysis is limited to the first two moments.
Coming back to the first-order moments, from (3.2) the characteristic equation looks like where from (2.2) . The roots of (3.7) are In examining the eigenvalues given by (3.8), it is necessary to distinguish between the case of real eigenvalues and that of complex conjugate ones. The eigenvalues are real numbers when where for the sake of simplicity β is assumed to be positive. Otherwise, they are complex conjugate numbers. Clearly, in the latter case the stability condition for the first moments is Equation (3.10) requires that the oscillator is damped and the amount of damping depends on both the oscillator frequency and the intensity of the exciting noise. This requirement may be rather restrictive when is not small.
Now, let us consider the case in which (3.9) is satisfied. From (3.8) it is seen that there is passage to instability when the eigenvalue with the sign plus before the square root becomes zero, that is, Solving (3.11) with respect to , we obtain the critical value of this quantity: where β must be larger than the value in the right-hand side of (3.9). By comparing (3.12) with (3.10) it is found that the critical value given by the former is always smaller than that given by the latter. Keeping into account that , the critical intensity of the excitation depends on the form of power spectral intensity of this: in Section 4 the cases of white and coloured noises will be considered.
As regards the second moment stability, since the moment equations are three, (3.3), the characteristic equation associated with the matrix is of the third order. The roots of such a type of equation have analytical expressions. However, they are rather cumbersome as many parameters enter them. Thus, it has been preferred to proceed numerically. By using a computer algebra software a parameter is varied till one out of the three eigenvalues becomes zero: the corresponding value of the parameter is the critical one. Routh-Hurwitz criteria [4, 29] are not used. This is why in another study devoted to the stability of elastic columns with memory-dependent damping  it has been found that these criteria may overevaluate the critical values by a 30%.
4. Stability Analyses
The present section is devoted to the applications of the theory previously explained. Three cases are considered for the parametric excitation : (1) stationary Gaussian white noise; (2) coloured Gaussian process obtained by passing a stationary Gaussian white noise through a first-order linear filter; (3) coloured Gaussian process obtained by passing a stationary Gaussian white noise through a second-order linear filter. This choice makes the comparisons easier as the three stochastic processes have a common parameter that is the intensity σ2 of the white noise.
The power spectral density (PSD) of a stationary stochastic process is defined as where is the autocorrelation function of the process and . Since the autocorrelation function of a stationary Gaussian white noise can be expressed as ( is Dirac’s delta), its PSD is , constant on the whole real axis. Thus, for a white noise the first criterion of stability for the first-order moments (3.10) gives The critical value of σ2 increases linearly with the system ratio of critical damping ζ0 and is inversely proportional to system pulsation , but it does not depend on Ω.
Now, the white noise is passed through a first-order linear filter (the Langevin equation) to have the excitation , that is, where has unit strength. The PSD of is The critical value is found by solving (3.10) with respect to and taking the expression of this quantity into account, that is, The result of (4.5) is plotted in Figure 1 as a function of Ω for , , and . It is found that is an increasing function of Ω.
The third type of excitation is obtained by passing the white noise through a second-order linear filter, that is, In this case the PSD of isfor the first condition of stability is computed from The results of (4.8) are plotted in Figure 2 for , , and ; the top plot refers to and the bottom plot to . Both plots show resonance with a marked minimum when Ω equates .
As regards the second condition of stability, from (3.12) is plotted in Figure 3 for , . Since this condition is valid only when is larger than the right-hand side of (3.9), the plots are drawn for with . Equation (3.12) gives a result with physical meaning only when is positive. From the plots—the bottom is a detail—it is seen that the interval of validity is small, and it narrows as κ increases. All the curves have a peak for .
Before presenting the results of the stability analyses for the second-order moments, it is necessary to compare the three PSD functions used in the analyses: they are depicted in Figure 4.The parameter is common to the three PDFs, and it constitutes the strength of the white noise. The plots clearly show as the white noise is by far the most severe excitation that may excite a dynamic system, if the colored excitations are obtained by passing it through a linear filter. As a consequence of this only observation, it cannot produce stabilizing effects on a dynamic system. The Wiener process yielded by the first-order filter of (4.3) (the Langevin equation) is the second in order of severity. As regards the processes generated by the second-order filter of (4.6), they have less strength as the abscissae of the peaks move farther from the origin. These remarks allow to explain the results of the analyses keeping in mind that the excitation affects the moments equations through and . The following anal-yses are of two types: (1) search of the critical value of keeping β, the amplitude of the sinusoidal term in (2.1), constant; (2) search of the critical value of β keeping constant. For the sake of simplicity in all analyses the parameter is worth one. It is recalled that there is passage to instability when an eigenvalue becomes zero, which happens when equates or β equates .
For white noise excitation the plots of and as a function of Ω are in Figures 5 and 6, respectively, where = 2π and ζ0 = 0.05. Both and are constant with respect to Ω, being the stable regions below the straight lines. This result is not surprising: in fact, the excitation enters the moment equations through and , which are equal and do not vary with Ω in the case of white noise, whose PSD is constant on the whole frequency axis. increases as decreases (Figure 5), and has a similar trend (Figure 6).
Figure 7 shows as a function of time for , (top plot), and (bottom plot). The moment equations (3.3) are numerically integrated by means of a fourth-order Runge-Kutta method. The initial perturbation is given by means of the initial conditions . In the former case is lesser than , and starts from the prescribed value 0.1, then it decays to zero and does not show any oscillation. In the latter case is larger than , and grows without limit. The absence of oscillations is due to the fact that in applying the stochastic averaging the oscillatory terms are cancelled.
As regards the excitation given by the first-order filter of (4.3), the plots of and as a function of are shown in Figures 8 and 9, respectively, where , , and . Both quantities increase as increases. In fact, the PSD (4.4) is a monotonically decreasing func-tion of , which causes the excitation to diminish with . The curve in the plot of Figure 8 is obtained by keeping β equal to 0.1. The increase of is marked and reaches a 70% as Ω passes from 0.5π to 6.5π. The curve of refers to : in the same interval the increase is less marked and amounts to about 28%. The analyses are made in a discrete series of values of Ω obtaining a set of points in the plane (σ2, Ω) or (β, Ω). Then, the curves are traced by means of the routine Spline of MAPLE.
The plots for the excitation given by the second-order filter (4.6) are shown in Figures 10 and 11 and are obtained for this data set: , = 2π, , in Figure 10 for , in Figure 11 for , , black line , blue line . In both plots the curve deriving from is the higher since the excitation strength is smaller, as already explained. In both plots there is a marked valley for : since (Ω) assumes its largest value (see Figure 4), there is a kind of stochastic resonance , even if the system remains stable.
5. Summary and Conclusions
In this paper the issue of stability of the stochastic Mathieu oscillator is addressed. In order to solve the problem, a suitable coordinate transformation and stochastic averaging are applied to the original dynamic system. In this way, two first-order stochastic differential equations are obtained. Then, they are transformed into two Itô-type stochastic equations. For such a type of stochastic differential equations Itô’s differential rule is applicable and allows one to derive the ODEs ruling the time evolution of the response statistical moments.
In the dynamic system considered here the excitation is merely linearly parametric. Hence, the moments ODEs are linear and homogeneous, and the ODEs for the moments of order are uncoupled from the ODEs for moments of different order. Because of this feature the ODEs for the moments of order admit zero solution, which means that the system is at rest, and the parametric excitation does not disturb this state. The system is stable in moments when these return to zero after the application of an external perturbation; otherwise, it is unstable and the moments diverge. The necessary and sufficient condition for the stability of the moments of order is that all the eigenvalues of the matrix of the coefficients of the ODEs written for those moments have negative real part. The passage to the instability requires that an eigenvalue becomes zero. In this way, the stochastic problem is reduced to the deterministic problem of the study of the eigenvalues of a matrix.
The procedure developed in this paper is applied to three types of Gaussian stationary stochastic excitations: (1) white noise; (2) white noise filtered through a first-order filter; (3) white noise filtered through a second-order filter. In order to make it possible to compare the results, in cases (2.2) and (2.3) the source white noise is the same as that of case (1), and it is by far the more severe excitation (Figure 4).
As regards the stability of the first-order moments, it is substantially led by the system damping. Several parameters affect the stability of the second-order moments: since the characteristic equation of the matrix is of the third order, so that its roots have complicated expressions, it has been chosen to proceed numerically. The analyses were aimed at finding the critical values of the white noise strength σ2 and of the amplitude β of the sinusoidal term in the motion equation (2.1). Differently from the deterministic Mathieu oscillator the pulsation Ω of this term has little or no effects on the bounds of stability, which is due to the presence of damping. Vice versa, in the case of excitation deriving from a second-order filter there is a stochastic resonance when Ω equates the peak frequency of the PSD of the excitation. In any case, the critical values and are by far smaller for the white noise excitation since as observed it is the strongest excitation. Thus, the statement of [9, 31] that a white noise excitation may cause a stabilizing effect is not found true. The results of these papers agree with those of : a direct comparison of the results is not possible as the system studied in that reference is a little different and a quite different method of analysis is used there.
- R. Campbell, Théorie Générale de l’Equation de Mathieu et de Quelques Autres Equations Differentielles de la Mécanique, Masson, Paris, France, 1955.
- N. W. McLachlan, Theory and Application of Mathieu Functions, Dover Publications, Mineola, NY, USA, 1962.
- S. P. Timoshenko, Theory of Elastic Stability, Dover Publications, Mineola, NY, USA, 2009.
- L. Meirovitch, Methods of Analytical Dynamics, Dover Publications, Mineola, NY, USA, 2003.
- R. L. Stratonovic and Y. M. Romanovski, “Parametric effect of a random force on linear and nonlinear oscillatory systems,” in Non-Linear Transformation of Random Processes, P. I. Kuznetsov, R. L. Stratonovic, and V. I. Tikhonov, Eds., p. 341, Pergamon Press, New York, NY, USA, 1965.
- M. F. Dimentberg, Statistical Dynamics of Nonlinear and Time-Varying Systems, Research Studies Press, Taunton, UK, 1988.
- G. O. Cai and Y. K. Lin, “Nonlinearly damped systems under simultaneous broad-band and harmonic excitations,” Nonlinear Dynamics, vol. 6, no. 2, pp. 163–177, 1994.
- H. Rong, G. Meng, X. Wang, W. Xu, and T. Fang, “Invariant measures and Lyapunov exponents for stochastic Mathieu system,” Nonlinear Dynamics, vol. 30, no. 4, pp. 313–321, 2002.
- W.-C. Xie, Dynamic Stability of Structures, Cambridge University Press, New York, NY, USA, 2006.
- Y. K. Lin and G. Q. Cai, “Stochastic stability of nonlinear systems,” International Journal of Non-Linear Mechanics, vol. 29, no. 4, pp. 539–553, 1994.
- Y. K. Lin and G. Q. Cai, Probabilistic Structural Dynamics: Advanced Theory and Applications, McGraw-Hill, New York, NY, USA, 1995.
- C. Floris, “Stochastic stability of a viscoelastic column axially loaded by a white noise force,” Mechanics Research Communications, vol. 38, no. 1, pp. 57–61, 2011.
- K. Itô, “On stochastic differential equations,” Memoirs of the American Mathematical Society, vol. 4, pp. 289–302, 1951.
- K. Itô, “On a formula concerning stochastic differentials,” Nagoya Mathematical Journal, vol. 3, pp. 55–65, 1951.
- M. Di Paola, “Stochastic differential calculus,” in Dynamic Motion: Chaotic and Stochastic Behaviour, F. Casciati, Ed., vol. 340, pp. 29–92, Springer, Vienna, Austria, 1993.
- M. Grigoriu, Stochastic Calculus, Birkhäuser Boston, Boston, Mass, USA, 2002.
- R. L. Stratonovich, Topics in the Theory of Random Noise, vol. 1, Gordon and Breach, New York, NY, USA, 1963.
- R. L. Stratonovich, Topics in the Theory of Random Noise, vol. 2, Gordon and Breach, New York, NY, USA, 1967.
- J. B. Roberts and P.-T. D. Spanos, “Stochastic averaging: an approximate method of solving random vibration problems,” International Journal of Non-Linear Mechanics, vol. 21, no. 2, pp. 111–134, 1986.
- W. Q. Zhu, “Stochastic averaging methods in random vibration,” ASME Applied Mechanics Review, vol. 41, no. 5, pp. 189–197, 1988.
- N. N. Bogoliubov and Y. A. Mitropolsky, Asymptotic Methods in the Theory of Non-Linear Oscillations, Gordon and Breach, New York, NY, USA, 1961.
- E. Mettler, “Stability and vibration problems of mechanics under harmonic excitation,” in Dynamic Stability of Structures, G. Herrmann, Ed., Pergamon Press, Oxford, UK, 1967.
- R. Z. Khasminskii, “A limit theorem for the solution of differential equations with random right hand sides,” Theory of Probability and Applications, vol. 11, no. 3, pp. 390–405, 1966.
- T. T. Soong and M. Grigoriu, Random Vibration of Mechanical and Structural Systems, Prentice Hall, Englewood Cliffs, NJ, USA, 1993.
- E. Wong and M. Zakai, “On the relation between ordinary and stochastic differential equations,” International Journal of Engineering Science, vol. 3, no. 2, pp. 213–229, 1965.
- R. L. Stratonovich, “A new representation for stochastic integrals and equations,” SIAM Journal on Control and Optimization, vol. 4, no. 2, pp. 362–371, 1966.
- L. Arnold, “A formula connecting sample and moment stability of linear stochastic systems,” SIAM Journal on Applied Mathematics, vol. 44, no. 4, pp. 793–802, 1984.
- M. Di Paola, “Linear systems excited by polynomials of filtered Poisson pulses,” ASME Journal of Applied Mechanics, vol. 64, no. 3, pp. 712–717, 1997.
- N. G. Chetayev, The Stability of Motion, Pergamon Press, New York, NY, USA, 1971.
- H. Calisto and M. G. Clerc, “A new perspective on stochastic resonance in monostable systems,” New Journal of Physics, vol. 12, Article ID 113027, 2010.
- S. T. Ariaratnam and D. S. F. Tam, “Parametric random excitation of a damped Mathieu oscillator,” Journal of Applied Mathematics and Mechanics (ZAMM), vol. 56, pp. 449–452, 1976.
Copyright © 2012 Claudio Floris. 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.