Nowdays mechanical engineering application needs to know the fractional calculus show us the paper the experimental and theoretical results of the dynamic analysis of single degree of freedom systems, SDFS, are presented and compared among them. In this particular study the differential equation used to obtain the dynamic response of a simple degree of freedom system was determined from the Kelvin-Voigt fractional derivative rheological models, KVFRdm. Current practice in seismic engineering to analyze viscoelastic problems such as the one mentioned here, resort to ordinary differential equations, ODE. However, in the last decades, it has been found that the use of fractional differential equations, FDE, lead to more comprehensive evaluations of time-dependent problems such as the SDFS. In this paper the values of the mass, the stiffness and damping of the SDFS considered in this study, are deduced from the experimental results obtained from forced vibrations induced by the hydraulic shaking table testing of a single degree, HSTTSDF, developed and certified in the Engineering Institute in Geotechnical Coordination, as well as from their analytic interpretation. After realizing simulations of the SDFS using FDE the results were compared with those obtained from the HSTTSDF. It is worth to stress the fact that in this study only the damped element of the fractional differential equation, KVFRdm, was considered fractional. Accordingly, a coupled behaviour between the stiffness and damping is assumed; in other words, the fractionary KVFRdm has the ability of modeling non structured element such as polymers, soils, and others, with such versatility of varying the proportion stiffness-damping as a function of the fractionary exponent $\alpha$ varying within the interval $(0<\alpha <1)$. On the basis of the results included in this paper it can be concluded that using the FDE more trustable results are obtained, which might have a direct influence on the procedures presently carried out in the seismic engineering practice.

Recent research shows results of phenomena using ordinary differential equations and reproduces the phenomenon in a way somewhat far from reality. On the other hand, fractional calculus offers mathematical tools with great potential to simulate these complex phenomena. Readers interested in the theory, modeling, and applications of fractional calculation to various types of problems can consult the following references [1], [2],[9], [11] and others. From the results shown in these references that include a great diversity of topics, we can conclude that fractional differential equations offer the ability to simulate with greater precision and simplicity the phenomena of great complexity, contrary to the situation that occurs when trying to reproduce them using classical differential equations.

Fractional derivative model are widely applied to easier fitting of real process and phenomena studied in physics, engineering, population dynamics and others. The fractional calculus show us an alternative way to find the dynamic behaviour and the importance to create a real model like a SDFS to show the behaviour and compare with fractional solution.

The research tries to represent and model one of the famous S1GL models using materials and geometries in order to show the dynamic behaviour in the system and compare it with the differential equations integer, $\alpha =1$, or fractional in case the answer differs from the integer. References for this type of experimental research used were not found; however, there are articles related to the fractional analytical damped oscillator [3], [7], among others.

Fractional calculus treats fractional derivatives and integrals of arbitrary order. Ordinary calculation has integer derivatives:$\frac{d^n}{d{x}^n}$, $n=0,1,2,...$ If n is changed for an arbitrary quantity $\alpha$ then we have the notation $\frac{d^{\alpha }}{d{x}^{\alpha }}$ where $\alpha$ can be real: positive integer (ordinary calculation), positive fractions, negative integers, negative fractions or even complex numbers. It is important to note that “fractional” does not mean just fractions, but an arbitrary amount including integers, fractions and complex numbers in general. The notation of a fractional integral is denoted $D^{-n}$. When $\alpha$ has an arbitrary value, there is no physical interpretation for $D^{\alpha }$ or $D^{-\alpha }$.

The birth of the fractional calculus occurs when in a letter of L’Hôpital, dated September 30, 1695, addressed to Leibniz who was asked of the notation used in his publications for the umpteenth derivative of a function, $\frac{D^{n}f\left(x\right)}{D{x}^n}$, then, what will result if $n\ne 1$. Leibniz replied “an apparent paradox of something that will be extracted and one day will have“useful consequences”. In the case presented in this paper, the fractional exponent allows to reproduce with better approximation the dynamic response of systems of a degree of freedom than the integer solution.

The Hydraulic Shaking Table Testing of a Single Degree, HSTTSD, used in this research was developed for this purpose in the Geotechnical Coordination of the Engineering Institute UNAM, II-UNAM. The HSTTSD consists of a mechanical-hydraulic control system and a calibration-data acquisition program.

In general, the operation of the HSTTSD, consists of defining a model with a maximum weight of $1ton-force$, then the type of displacement is defined by harmonic (Figure 2) or random (Figure 3) function that is measured with a displacement sensor, LDT, and unidirectional accelerometer. The shaking table has a rectangular geometry $2.1m*2m$ longitudinal and transverse directions, respectively. It has a maximum peak-peak displacement of $18cm$, a frequency range of $0.25-10Hz$ and has a lightweight platform, with weight of $80kg$.

The fractional differential equation of motion is proposed from the fractional rheological model corresponding to the Kelvin-Voigt fractional derivative rheological models where the common damper is replaced by a fractional damper, triangular element in Figure 6, called “spring-pot or KVFRdm”. This represents a combination of springs and dampers in a ladder-type arrangement as shown in Figure 7.

The fractional differential equation of motion according to the dynamic balance of forces acting in the system: inertial forces, stiffness, damping and dynamic action. Thus, according to the rheological model of Figure 6 we have:

where m, c and k are the body mass, system damping coefficient and system stiffness, respectively, $F\left(t\right)$ is a dynamic excitation and $\alpha$ is the order of fractional derivative, $0<\alpha <1$.

The Caputo definition was used in this research. Further note that the initial conditions and the parameters is completely know. To know about the definition or computation of fractional derivatives can be found in a classical text [6] [8] [9] is recommended.

It is important to note that the term $\frac{d^{\alpha }x\left(t\right)}{d{t}^{\alpha }}$ represents a combination of velocity and displacement, so its units are a function of the fractional exponent [$\frac{1}{s^{\alpha }}$, s=seconds], which makes the fractional differential equation inhomogeneous from the point of view of its units. However, the parameter ${\lambda }^{(1-\alpha )}$ is entered in the $\frac{d^{\alpha }x\left(t\right)}{d{t}^{\alpha }}$ term of the fractional differential equation. Therefore the units of the parameter $\lambda$ are $s^{(-1)}$ ($\frac{1}{s};s=seconds$) and its value is equal to one, so it is ensured that the units of the term $\frac{d^{\alpha }x\left(t\right)}{d{t}^{\alpha }}$ do not change when the value of the exponent $\alpha$ changes; hence, the fractional differential equation remains homogeneous from the point of view of its units.

The SDFS is a mainly two-piece designed model by stainless steel, one is the bar and the other is the weight and both are fixed with a rope, then the bar is built-in into the base of the SDFS. The bar is composed of a bar of circular section, diameter $3.6mm$, height $13cm$ and the mass is a prism with a quadrangular section, the dimension by side $5cm$, height $3.63cm$ and the weight $0.7155kg$. The built-in SDFS is observed in Figure 8.

The ODEM and FDEM were solved using the procedure of the numerical methods with the platform of the software "Matlab" specifically the evaluation of the derivation was carried out using the application proposed in Simulink-Matlab by Nintiger [5].

The experiment was carried out in the HSTTSD and the experimental response was compared with the experimental solution obtained against the FDEM considering $\alpha \in (0.3,0.4,...,1)$ intervals 0.1. This is intended to visualize the approximation when using FDE applied to the solution of a SDFS dynamic response.

The experimental test consisted of using a sinusoidal displacement with frequencies of $3.0Hz$, amplitude of $1.15cm$ and duration of $33seconds$. Figure 11 shows the end of the displacement history $(32-33s)$ corresponding to $3.0Hz$. It should be noted that the HSTTSD displacements was defined as MTS_DISP.

Figure 12 show the acceleration history, it was obtained from the HSTTSD platform and the total acceleration of SDFS; the identification of these responses in the figure is HSTTSD_f3hz where the excitation has a frequency of $3.0Hz$. Likewise, the mass response was identified as Exp_f3M1H13_TA for the frequencies of $3.0Hz$.

Figure 13 shows the results of the SDFS response in terms of the history of accelerations for the case of excitation with a frequency of 3 $Hz$. The graph shows the experimental behaviour, Exp_f3M1H13_TA, of the total accelerations of the mass, the numerical result obtained with the fractional numerical results using the fractional differential equation with fractional exponents $\alpha \in (0.3,0.4,\dots ,1)$ in intervals of $0.1$, At_TF$=0.3$ to $1.0$. It is important to note, the graph included in Figure 13 only present the behaviour corresponding to the second 26 to 27 section of the total responses of $33s$, where the experimental values show greater discrepancies between the experimental and analytical values At_TF$=1.0$. The graph in this figure show the concordance between the experimental results and those calculated with the numerical solution of the fractional differential equation compared with the whole differential equation and while the fractional exponent value decreases, the severity of the relative (total) accelerations increases appreciably. Thus, the fractional exponent $\alpha$ is directly related to the total damping of the system. The graph show the value of the exponent $\alpha$ equal to $0.9$, A_Tf$=0.9$, reproduces with very good approximation to the experimental amplitudes resulting. Consequently, for this particular case, the fractional differential equation with the value of the exponent $\alpha =0.9$ allows to reproduce the experimental response of SDFS with a better approximation. Also, Figure 14 shows coincidences and discrepancies between the experimental response and the fractional model; however, it presents a very good agreement.

The way to evaluate the robustness of the value of the exponent $\alpha =0.9$ is optimal to reproduce the dynamic response of the SDFS. The responses monitored and obtained from the fractional solution with fractional exponents $\alpha$ between $0.3$ to $1.0$, using increments of $0.1$ are shown in Figure 13, where it indicates that the exponent $\alpha$ has an important influence on the response, it also shows that the $0.9$ value of the fractional exponent remains the one that allows the best approximation. Note, by increasing the frequency from $3.0Hz$, the behaviour of SDFS, the fractional exponent $\alpha =0.9$ continued to approach the experimental response better. Then, once the value of the exponent is defined, it retains its mechanical properties of SDFS and can continue to be used to calculate the response of the exponent but subject to similar excitations and with significantly different frequencies; the vibration characteristics of HSTTSD defects or incorrect mounting also contribute to generate of untoward vibration.

It is important to mention, the Figure 13 include the analytical integer solution obtained with the whole differential equation, A_Tf$=1.0$, and the analytical solutions calculated with the fractional differential equation solution, A_Tf$=0.9$, as well as the experimental measurements. This allows to directly compare the results. In general, the Figure 14 show the solution with the fractional differential equation is more attached to that recorded experimentally.

The results clearly show the direct effect of the $\alpha$ value on the system response. For $\alpha =0$, the system behaves as a body without damping, while for $\alpha =1$, the damping of the system depends on the speed, $\frac{dx\left(t\right)}{dt}$; so the body presents a behaviour as a viscous element. For intermediate values $0<\alpha <1$, complex responses are obtained since the term of the damping of the fractional differential equation jointly presents the characteristics of a stiffness and viscous body.

This research present an easy mathematical model capable to obtain the dynamic response of a SDFS of a fractional damping viscoelastic system, where the fractional derivative can be or any order. The experimental test presented above, its a justification to use a fractional calculus models of viscoelastic material behaviour is able to predict accurately the dynamic behaviour of SDFS during the test. The results of the simulation are compared with the experimental ones, concluding that the tendencies observed in the measurement have been reproduced. In this particular case, the fractional derivative is found to be or order $\alpha =0.9$.

To achieve competitive models it is necessary to create experimental models capable of representing the real world to compare, analyze and study its behaviour; in this particular case it was observed that the dynamic behaviour of S1GL should coincide with the whole ODEM; however, it did not happen that way, despite trying to make a perfect model and perform the experiment as best as possible. It is necessary to continue with experimental research and create models that are easy to apply and compare with reality using fractional calculus.

The need to generate results, from experimental results, and translate them to a theory is of vital importance. At present, this task is facilitated since it has had the support of numerical methods with which fractional differential equations can be solved, giving results similar to real systems without taking into account greater number of variables with more complex differential equations.

The research carried out allows to conclude that to use the fractional differential equations applied to truly representative models that are suitable to study and model the behaviour of the SDFS subjected to vibrations of any type. These conclusions can extend to problems in physical, engineering dependent on time such as propagation of waves in terrestrial media, fluid flow in porous media, consolidation of compressible soils, among others.