Next Article in Journal
Null Controllability of Hilfer Fractional Stochastic Differential Inclusions
Next Article in Special Issue
A Review on Fractional-Order Modelling and Control of Robotic Manipulators
Previous Article in Journal
Solvability of Nonlinear Impulsive Generalized Fractional Differential Equations with (p,q)-Laplacian Operator via Critical Point Theory
Previous Article in Special Issue
Fractional-Order Impedance Control for Robot Manipulator
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Using a Fully Fractional Generalised Maxwell Model for Describing the Time Dependent Sinusoidal Creep of a Dielectric Elastomer Actuator

Laboratory for Robotisation, Faculty of Mechanical Engineering, University of Maribor, 2000 Maribor, Slovenia
*
Author to whom correspondence should be addressed.
Fractal Fract. 2022, 6(12), 720; https://doi.org/10.3390/fractalfract6120720
Submission received: 26 September 2022 / Revised: 30 November 2022 / Accepted: 2 December 2022 / Published: 4 December 2022
(This article belongs to the Special Issue Applications of Fractional-Order Calculus in Robotics)

Abstract

:
Actuators made of dielectric elastomers are used in soft robotics for a variety of applications. However, due to their mechanical properties, they exhibit viscoelastic behaviour, especially in the initial phase of their performance, which can be observed in the first cycles of dynamic excitation. A fully fractional generalised Maxwell model was derived and used for the first time to capture the viscoelastic effect of dielectric elastomer actuators. The Laplace transform was used to derive the fully fractional generalised Maxwell model. The Laplace transform has proven to be very useful and practical in deriving fractional viscoelastic constitutive models. Using the global optimisation procedure called Pattern Search, the optimal parameters, as well as the number of branches of the fully fractional generalised Maxwell model, were derived from the experimental results. For the fully fractional generalised Maxwell model, the optimal number of branches was determined considering the derivation order of each element of the branch. The derived model can readily be implemented in the simulation of a dielectric elastomer actuator control, and can also easily be used for different viscoelastic materials.

1. Introduction

Dielectric elastomer actuators, also called DEAs, can have different structures. The first DEAs had parallel or cylindrical structures [1,2]. From then on, many different structures were developed: spring roll actuators, helical actuators, and stacked actuators where the actuators are stacked on top of each other [3,4,5]. All these actuators have the same basic structure and activation principle. Thus, regardless of whether the DEAs are parallel or cylindrical, their activation principle is the same. Their structure is like that of a capacitor. It has an upper and a lower electrode and an elastic dielectric. Some special rules apply to both the electrodes and the dielectric. The dielectric must be elastic. The electrodes must be conductive, and their mechanical structure must be the same or similar to that of the dielectric [6]. If their mechanical properties are different, the DEA may lose its mobility and efficiency. Usually, conductive pastes of carbon are chosen for the electrodes.
The material properties of the selected elastomer VHB 4910 are viscoelastic, which means that the material has both viscous and elastic properties. When the DEA is subjected to a sinusoidal load, it exhibits sinusoidal creep behaviour. When a high voltage is applied, coupling occurs between the electrical and mechanical forces. Suo et al. were the first researchers to describe the electromechanical coupling in detail in a DEA [7]. In the work of Gu et al. [8], the viscoelastic behaviour was captured using continuum mechanics and nonequilibrium thermodynamics, where the Helmholtz free energy along with the Gent model was used to derive the constitutive equations. The work shows a good match between the experimental and simulation results. However, extensive knowledge in the above areas is required, and the model cannot be integrated easily into the control simulations. The Prandtl–Ishlinskii model and the modified Prandtl–Ishlinskii model are used in the work of Zuo et al. [9,10]. The work is based on the fourth order polynomial to describe the asymmetric behaviour, and fixed weights with thresholds of the play operators are used to describe the rate dependence, while the second order derivative of the input voltage is introduced into the fourth order polynomial to describe the peak-to-peak shifts that depend on the frequencies. This is a complex process that results in a good match between the experimental data and the calculated responses. The Prandtl–Ishlinskii model can be used as a feedforward compensator. Wissler et al. [11] used the Prony series theory to model the time dependent viscoelastic behaviour of the elastomer used in soft actuators. The drawback of this method is the inability to use it in the control simulations. A standard linear solid rheological model which satisfies the thermodynamic consistency was also used to capture rate-dependant behaviour of the soft materials in the work of Zhao et al. [12]. This model could capture the rate dependent mechanical behaviour of soft materials. However, the introduction of the internal variables in the modelling process significantly increases the complexity of the proposed constitutive model and the number of governing equations.
Fractional calculus has proved to be a powerful mathematical tool in the approach used to modelling the time dependent mechanical behaviour of viscoelastic materials. Its advantages lie in the remarkable reduction in model parameters where fractional orders of derivatives are used. This can be especially seen in the field of fractional viscoelasticity, including soils [13], polymers [14,15] and construction materials [16]. Xu et al. [17] used fractional constitutive models of a fractional Kelvin–Voigt model, a fractional Maxwell model and a fractional Poynting–Thomson model, which is a springpot connected in series with two parallel springpots. In the work of Barretta et al. [18], hereditariness and nonlocality of bending problems were presented with the help of fractional operators. Time-dependent hereditary behaviour, which is typical of viscoelastic materials, has been modelled with the help of a springpot, a fractional Kelvin–Voigt and a fractional Maxwell model.
To capture the sinusoidal creep behaviour, the fully fractional generalised Maxwell model is derived in Section 2. The generalised fractional Maxwell model has already been used in the work of Luo et al. [19] to determine the storage and loss modulus, but not with the Laplace transform, nor to derive the governing equation of the DEA. Our model is the first to describe the material behaviour of the DEA using the Laplace transform. The fully fractional generalised Maxwell model was used because it does not require complex knowledge of continuum mechanics and thermodynamics. It can readily be derived using the Laplace transform to obtain the transfer function of the system. Using the lsim function in the Matlab software, the excitation voltage with three different frequencies can be implemented in the transfer function of the system, and the response of the DEA can be calculated easily. An optimisation procedure using the Pattern Search global optimisation solver was used to derive the material parameters based on the experimental results. Section 3 shows the results of the optimisation procedures. Section 4 summarises the conclusions.

2. Materials and Methods

2.1. Principles of the DEA

The structure of the DEA is similar to that of a capacitor. It has an upper and lower conductive electrode and an elastic dielectric. The structure of the DEA can be seen in Figure 1. Usually, elastic materials are chosen for dielectrics, such as VHB tapes. VHB 4910 tape was used in this study [20]. A high voltage DC is required to activate the DEAs. When the voltage is applied, the electrical force also known as the Maxwell force is generated between the upper and lower conductive electrodes, causing the DEA to contract vertically and expand longitudinally, as elastomer is known to be incompressible.
The Maxwell force can be calculated as
F M a x = F e l = ε 0 ε r E 2 l 1 l 2 = ε 0 ε r V l 3 2 l 1 l 2 N
ε 0   absolute   permittivity ε r = 4.7   relative   permittivity E   electric   field V   voltage   l 1 , l 2 , l 3 dimensions   of   the   DEA
The relative permittivity was chosen as 4.7 according to [8]. A conductive paste was chosen from BareConductive® [21].

2.2. Derivation of the Fully Fractional Generalised Maxwell Model

To derive the fully fractional generalised Maxwell model, one needs to use fractional derivatives. Fractional derivatives are derivatives that are not restricted to positive integers, but can be any real number. There are three definitions of fractional derivatives, namely, Riemann–Liouville, Caputo, and Gruenwald–Letnikov [22]. The Gruenwald–Letnikov definition is used, namely
D t p a f ( t ) = lim h 0 h p r = 0 m ( 1 ) r p r f ( t r h )
a , t   time   limits m   integer   order   of   derivation / integration p = p > 0 ,   derivation p < 0 ,   integration
Because it can handle the fractional derivation and integration easily. It is also suitable for the numerical calculation of fractional derivatives.
The generalised Maxwell model is used to describe constitutive models with viscoelastic effects. It consists of a spring connected in parallel to the branches of the Maxwell elements. The Maxwell element consists of a spring connected in series with a dashpot [15,16]. To convert this into a fully fractional generalised Maxwell model, all elements are replaced by the so-called springpot element [23,24]. The springpot element has two parameters, one of which represents the material properties and the other the derivative order. The derivative order is bounded between 0 and 1, since there is no known physical interpretation above 1 [24]. Figure 2 shows both models.
c i , 1 , c i , 2 material   properties   of   springpot α i , 1 , α i , 2   order   of   derivation k i   spring   constants m   mass   of   weight n   number   of   branches   with   the   Maxwell   element F e l   electrical   force   applied F g   weight   of   the   mass
The Laplace transform with Laplace operator is used to derive the governing equation of motion for the DEA using the fully fractional generalised Maxwell model. The Laplace transform turns the derivative into a multiplication and the integration into a division. The fully fractional generalised Maxwell model describes the material behaviour. The electric force is calculated using Equation (1). Since the weight is used to hold the DEA in the stretched position, its inertia must be included in the governing equation of motion. All branches of the fully fractional generalised Maxwell model are subject to the same displacement, as shown in Figure 3.
Each branch has two displacements because it has two springpots. The electric force and the force generated by the weight are distributed jointly between the branches of the fully fractional generalised Maxwell model. The governing equation of motion is calculated as
F = F 1 + i = 2 n + 1 F i + m x ¨ .
Displacement of each branch with the fractional Maxwell element is calculated as
x = x 2 , 1 + x 2 , 2 x = x 3 , 1 + x 3 , 2 . . . x = x i , 1 + x i , 2 n x = i = 2 n + 1 x i , 1 + i = 2 n + 1 x i , 2
The force in the first branch is calculated as
F 1 = c 1 , 1 d α 1 , 1 x d t α 1 , 1 .
The force in each branch containing two springpots is the same in each springpot, which is calculated as
F 2 = c 2 , 1 d α 2 , 1 x 2 , 1 d t α 2 , 1 = c 2 , 2 d α 2 , 2 x 2 , 2 d t α 2 , 2 F 3 = c 3 , 1 d α 3 , 1 x 3 , 1 d t α 3 , 1 = c 3 , 2 d α 3 , 2 x 3 , 2 d t α 3 , 2 . . . i = 2 n + 1 F i = i = 2 n + 1 c i , 1 d α i , 1 x i , 1 d t α i , 1 = i = 2 n + 1 c i , 2 d α i , 2 x i , 2 d t α i , 2
Expressing x i , 1 and x i , 2 from Equation (6) with the fractional integration, one gets
i = 2 n + 1 F i = i = 2 n + 1 c i , 1 d α i , 1 x i , 1 d t α i , 1 / d α i , 1 d t α i , 1 i = 2 n + 1 d α i , 1 F i d t α i , 1 = i = 2 n + 1 c i , 1 x i , 1 x i , 1 = i = 2 n + 1 1 c i , 1 d α i , 1 F i d t α i , 1 ; x i , 2 = i = 2 n + 1 1 c i , 2 d α i , 2 F i d t α i , 2 ;
Inserting the results from Equation (7) into Equation (4), one gets
n x = i = 2 n + 1 1 c i , 1 d α i , 1 F i d t α i , 1 + i = 2 n + 1 1 c i , 2 d α i , 2 F i d t α i , 2 .
Force   F i should be expressed from Equation (8). Laplace transformation is used, since force F i is part of the fractional integration with different orders of integration. Performing Laplace transformation on Equation (8) and expressing F i , one gets
i = 2 n + 1 F i ( s ) = n x ( s ) i = 2 n + 1 1 c i , 1 s α i , 1 + 1 c i , 2 s α i , 2 .
Performing the Laplace transform on Equation (3) and inserting Equation (9), one gets
F ( s ) = c 1 s α 1 x ( s ) + n x ( s ) i = 2 n + 1 1 c i , 1 s α i , 1 + 1 c i , 2 s α i , 2 + m s 2 x ( s ) .
Rearranging Equation (10) to get the transfer function, one gets
x ( s ) F ( s ) = 1 m s 2 + c 1 , 1 s α 1 , 1 + n i = 2 n + 1 1 c i , 1 s α i , 1 + 1 c i , 2 s α i , 2 .
Equation (11) represents the transfer function of the governing equations of motion for the fully fractional generalised Maxwell model with n branches, where the influence of the electric force and the force of the weight of the mass 316 g are combined in the equation. The add-on FOMCON is required to implement Equation (11) in Matlab [25]. Equation (11) is written in the relation force–displacement. To convert Equation (11) into the relation stress–strain, one needs to use
E = σ ε = F L A Δ l = k L A k = E A L
and
σ = η d α ε d t α η = σ d α ε d t α = F A d α Δ l L d t α = F A 1 L d α Δ l d t α = c L A c = η A L .

2.3. Experiments

The elastomer VHB 4910 was used to set up the experiment. The elastomer was cut to the initial dimensions of 49 mm × 50 mm × 1 mm, with only a 10 mm wide area used to expand the elastomer to the initial dimension, while the rest was used for clamping. The initial dimensions for the sinusoidal force excitation were 49 mm × 60 mm × 0.16 mm. The active area to which the conductive paste was applied was 40 mm × 60 mm.
Three different frequencies with an amplitude of 6 kV DC were used for the sinusoidal voltage excitation. The three frequencies were F 1 = 1 13 Hz , F 2 = 1 7 Hz , and F 3 = 1 5 Hz . These frequencies were chosen so that they were not multiples of each other. In this way, the frequencies are not associated with a common factor. The experiment lasted 155 s to capture the sinusoidal creep behaviour of the DEA. The displacement of the DEA was measured using a Wenglor laser sensor [26]. The experimental setup is shown in Figure 4.

2.4. Optimisation

To obtain the material parameters of the fully fractional generalised Maxwell model, the experimental results were optimised with the model. For the optimisation procedure, the global optimisation solver Pattern Search was chosen, which is also known as a direct search method in Matlab software. The flow of the optimisation procedure is shown in Figure 5. First, the user must specify the number of branches for the fully fractional generalised Maxwell model and the initial parameters for the model.
The electric and mechanical forces are calculated as the sum of the Maxwell force given by Equation (1) and the weight of the dead mass. In Equation (1), V is replaced with the calculation of the appropriate voltage regarding the frequency in use. Using the Matlab function lsim, the response of the transfer function of the fully fractional generalised Maxwell model to the electric and mechanical forces can be calculated easily. The responses for all frequencies were calculated and compared with the experimental results using the least squares method. The efficiency of the fully fractional generalised Maxwell model with a different number of branches was calculated using the method of the coefficient of determination, also known as the R 2 method, calculated as
R 2 = 1 S S r e s S S t o t .
S S t o t = i ( y i y ¯ ) 2   Total   sum   of   squares .   S S r e s = i ( y i f i ) 2   Residual   sum   of   squares .   y i   Measured   data . y ¯   Average   data . f i   Calculated   data .
After the optimisation procedure was completed, the best results were recorded and the final value of R 2 was calculated for each frequency. At the end, the average value of R 2 was calculated as R m e a n 2 for the whole frequency range. The number of branches was chosen to be between 1 and 5. Finally, each frequency was optimised individually for the model, to compare the results when only one frequency was optimised to the results where the whole frequency range was optimised. It was examined which the essential parameters of the model were for each frequency. Table 1 shows all symbols used in the research.

3. Results

After the optimisation procedure was set up using Matlab software and the Pattern Search global optimisation algorithm, the R m e a n 2 values shown in Table 2 were obtained. As can be seen, the number of branches increased the accuracy of the model, but not significantly. It can also be seen that adding more than three branches did not affect the efficiency of the model. The initial parameters used in the optimisation, as well as the optimal parameters obtained from the optimisation, are shown in Table 3.
From Table 3, in some branches the order of derivation does not change from the initial values. The derivation orders of 1 represent dashpot elements, and if both elements of the branch have the order of 1, they can be combined into one dashpot element. The same applies for spring elements if the order of derivation equals 0.
Figure 6 shows the experimental and calculated results from the optimised parameters of the fully fractional generalised Maxwell model for one to five branches. Adding more than two branches did not improve the efficiency of the model. The best matching between the data is seen for the middle frequency of 1/7 Hz, where the matching was 88%. For the frequency of 1/5 Hz, the worst matching was achieved of only 12.9%.
The optimisation of the model for an individual frequency was performed after the optimisation for the whole frequency range was carried out. From Table 2, it is seen that two additional branches of the fully fractional Maxwell elements are the most optimised. In Table 4, the initial and optimised parameters of only two additional branches are shown for each individual frequency for the fully fractional generalised Maxwell model. Each frequency demands its own material parameters, as well as different topology of the fully fractional generalised Maxwell model. Figure 7 shows matching between the experimental and calculated results if each frequency was optimised individually for the model.

4. Discussion

The fully fractional generalised Maxwell model was used for the first time to describe the material behaviour of a DEA. Table 2 shows that the number of branches contributes to the effectiveness of the model. However, the contribution of each additional branch was small. Increasing the number of branches up to two increased the effectiveness. Surprisingly, adding more branches did not improve the effectiveness of the model. It can be seen from Table 3 that when adding more than two branches, the optimised parameters were chosen such that the model could be reduced to the fully fractional generalised Maxwell model with only two branches. This was possible because the two α parameters within a branch were 1, which represents dashpots. Dashpots in series can be reduced to one dashpot. The reduced and rearranged model can be seen in Figure 8.
Concluding points can be itemized as follow:
  • The number of fully fractional Maxwell elements slightly affected the effectiveness of the model.
  • Adding more than two branches did not increase the effectiveness of the model.
  • The fully fractional Maxwell model was reduced to the model seen in Figure 8.
  • The middle frequency of 1/7 Hz had the best agreement of 0.88 between data.
  • Optimising each frequency individually drastically improved the overall agreement between data to 0.745.
  • Optimising each frequency individually has a drawback since each frequency requires its own material parameters.
  • Topology optimisation cannot be included into the Pattern Search algorithm.
The reduced and rearranged model has two springpots. The first springpot has the order of derivative equal to 0.2 which means it has 80% characteristics of a spring and only 20% characteristics of a damper. The other springpot has the order of derivative equal to 0.52 which means that nearly half of its characteristics are those of a spring, and half those of a damper.
The fully fractional generalised Maxwell model best describes the frequency of 1/7 Hz, where the match between the experimental and modelled response was 0.88, which was a good match. However, at the highest frequency of 1/5 Hz, the match between the data was the lowest and was only 0.12. The average match between the data over all three frequency ranges was 0.533, which is less than the methods used in the work of Gu et al. [8] and Zuo et al. [9] where agreement between data was more than 0.9. If only data from frequency of 1/7 Hz were compared to the data from the work of Gu and Zuo, then our method proves very efficient since it is much easier to implement it and use it in the control.
The initial and optimised parameters for the optimisation of the individual frequency with the fully fractional generalised Maxwell model are listed in Table 3. The match between the experimental and calculated responses was improved. The average match increased to 0.745, which is a good match and can be easily compared with the work of Gu and Zuo. However, the material parameters of the fully fractional generalised Maxwell model were different for each frequency, which was a drawback, since one would like to have universal material parameters for the entire frequency range.
Finding the optimal number of branches for the model cannot be included in the optimisation, because changing the number of branches changed the number of lower and upper bounds on the model. This is something that cannot be included in the optimisation solver for the Pattern Search. This can only be done by observing the material parameters of the individual branches within the model manually.

5. Conclusions

The match between experimental and calculated results was lower for the whole frequency range than in the work of others. Observing only the middle frequency, the method would be easily compared with the work of others. However, the fully fractional generalised Maxwell model can be derived and implemented easily, and the responses of the model can be determined quickly. Implementation of the model in simulation control is straightforward. The proposed method can easily be used on other materials with viscoelastic behaviour. In future work, topology optimisation could be included into the optimisation procedure.

Author Contributions

Conceptualisation, methodology, software, validation, visualisation, writing—review and editing, T.K., R.B. and J.G. All authors have read and agreed to the published version of the manuscript.

Funding

Research was co-financed by the Slovenian Research Agency under the Research Programme no. P2-0157, which operates within the framework of the University of Maribor, Faculty of Mechanical Engineering.

Data Availability Statement

The data that support the findings of this study are available upon reasonable request from the authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Youn, J.-H.; Jeong, S.M.; Hwang, G.; Kim, H.; Hyeon, K.; Park, J.; Kyung, K.-U. Dielectric Elastomer Actuator for Soft Robotics Applications and Challenges. Appl. Sci. 2020, 10, 640. [Google Scholar] [CrossRef] [Green Version]
  2. Bar-Cohen, Y. Electroactive Polymer (EAP) Actuators as Artificial Muscles: Reality, Potential, and Challenges; SPIE Press: Bellingham, WA, USA, 2001. [Google Scholar]
  3. Pei, Q.; Rosenthal, M.; Stanford, S.; Prahlad, H.; Pelrine, R. Multiple-degrees-of-freedom electroelastomer roll actuators. Smart Mater. Struct. 2004, 13, N86–N92. [Google Scholar] [CrossRef]
  4. Carpi, F.; Migliore, A.; Serra, G.; De Rossi, D. Helical dielectric elastomer actuators. Smart Mater. Struct. 2005, 14, 1210–1216. [Google Scholar] [CrossRef] [Green Version]
  5. Rui, Z.; Patrick, L.; Andreas, K.; Gabor, M.K. Spring roll dielectric elastomer actuators for a portable force feedback glove. In Proceedings of Smart Structures and Materials 2006: Electroactive Polymer Actuators and Devices (EAPAD), San Diego, CA, USA, 26 February–2 March 2006. [Google Scholar]
  6. Rosset, S.; Shea, H.R. Flexible and stretchable electrodes for dielectric elastomer actuators. Appl. Phys. A 2013, 110, 281–307. [Google Scholar] [CrossRef] [Green Version]
  7. Suo, Z.; Zhao, X.; Greene, W.H. A nonlinear field theory of deformable dielectrics. J. Mech. Phys. Solids 2008, 56, 467–486. [Google Scholar] [CrossRef]
  8. Gu, G.G.; Zhu, J.; Zhu, L.; Zhu, X. Modeling of Viscoelastic Electromechanical Behavior in a Soft Dielectric Elastomer Actuator. IEEE Trans. Robot. 2017, 133, 1263–1271. [Google Scholar] [CrossRef]
  9. Zou, J.; Gu, G. Modeling the Viscoelastic Hysteresis of Dielectric Elastomer Actuators with a Modified Rate-Dependent Prandtl–Ishlinskii Model. Polymers 2018, 10, 525. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Zou, J.; Gu, G. Dynamic modeling of dielectric elastomer actuators with a minimum energy structure. Smart Mater. Struct. 2019, 28, 085039. [Google Scholar] [CrossRef]
  11. Wissler, M.; Mazza, E. Mechanical behavior of an acrylic elastomer used in dielectric elastomer actuators. Sens. Actuators A Phys. 2007, 134, 494–504. [Google Scholar]
  12. Zhao, X.; Koh, S.; Suo, Z. Nonequilibrium Thermodynamics of Dielectric Elastomers. Int. J. Appl. Mech. 2011, 3, 203–217. [Google Scholar] [CrossRef]
  13. Xiang, G.; Yin, D.; Cao, C.; Gao, Y. Fractional description of creep behavior for fiber reinforced concrete: Simulation and parameter study. Constr. Build. Mater. 2022, 318, 126101. [Google Scholar] [CrossRef]
  14. Meng, R.; Yin, D.; Drapaca, C.S. Variable-order fractional description of compression deformation of amorphous glassy polymers. Comput. Mech. 2019, 64, 163–171. [Google Scholar] [CrossRef]
  15. Gao, Y.; Zhao, B.; Yin, D.; Yuan, L. A general fractional model of creep response for polymer materials: Simulation and model comparison. J. Appl. Polym. Sci. 2022, 139, 51577. [Google Scholar] [CrossRef]
  16. Su, X.; Chen, W.; Xu, W.; Liang, Y. Non-local structural derivative Maxwell model for characterizing ultra-slow rheology in concrete. Constr. Build. Mater. 2018, 190, 342–348. [Google Scholar]
  17. Xu, H.; Jiang, X. Creep constitutive models for viscoelastic materials based on fractional derivatives. Comput. Math. Appl. 2017, 73, 1377–1384. [Google Scholar] [CrossRef]
  18. Barretta, R.; Marotti de Sciarra, F.; Pinnola, F.P.; Vaccaro, M.S. On the nonlocal bending problem with fractional hereditariness. Meccanica 2022, 57, 807–820. [Google Scholar] [CrossRef]
  19. Luo, D.; Chen, H.-S. A new generalized fractional Maxwell model of dielectric relaxation. Chin. J. Phys. 2017, 55, 1998–2004. [Google Scholar] [CrossRef]
  20. 3M™. 3M™ VHB™ Tape Speciality Tapes. Available online: https://multimedia.3m.com/mws/media/986695O/3m-vhb-tape-specialty-tapes.pdf (accessed on 12 January 2021).
  21. Conductive, B. Electric Paint. Available online: https://www.bareconductive.com/shop/electric-paint-50ml/ (accessed on 11 July 2020).
  22. Podlubny, I. Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications; Elsevier Science: Amsterdam, The Netherlands, 1998. [Google Scholar]
  23. Mainardi, F. Fractional Calculus and Waves in Linear Viscoelasticity: An Introduction to Mathematical Models; Imperial College Press: London, UK, 2010; pp. 51–63. [Google Scholar]
  24. Mainardi, F.; Spada, G. Creep, relaxation and viscosity properties for basic fractional models in rheology. Eur. Phys. J. Spec. Top. 2011, 193, 133–160. [Google Scholar] [CrossRef]
  25. FOMCON. FOMCON Toolbox for Matlab. Available online: http://fomcon.net/ (accessed on 10 December 2020).
  26. Wenglor. YP05MGV80. Available online: https://www.google.com/url?sa=t&rct=j&q=&esrc=s&source=web&cd=1&ved=2ahUKEwihnaWa0bbjAhXiB50JHWepBx0QFjAAegQIBhAC&url=http%3A%2F%2Fwww.shintron.com.tw%2Fproimages%2Ftakex%2Fwenglor%2FYP05MGV80.pdf&usg=AOvVaw0hZdH3GNiqDLpm9BclkSVw (accessed on 15 July 2020).
Figure 1. Structure of the parallel DEA: (a) initial state and (b) activated state.
Figure 1. Structure of the parallel DEA: (a) initial state and (b) activated state.
Fractalfract 06 00720 g001
Figure 2. Two types of models used: (a) Fully fractional model and (b) Generalised model.
Figure 2. Two types of models used: (a) Fully fractional model and (b) Generalised model.
Fractalfract 06 00720 g002
Figure 3. Displacement of the elements.
Figure 3. Displacement of the elements.
Fractalfract 06 00720 g003
Figure 4. Experimental set up of the DEA.
Figure 4. Experimental set up of the DEA.
Fractalfract 06 00720 g004
Figure 5. Optimisation procedure.
Figure 5. Optimisation procedure.
Fractalfract 06 00720 g005
Figure 6. Experimental vs. calculated responses of the DEA for different numbers of branches (n) in the fully fractional generalised Maxwell model.
Figure 6. Experimental vs. calculated responses of the DEA for different numbers of branches (n) in the fully fractional generalised Maxwell model.
Fractalfract 06 00720 g006
Figure 7. Responses of individually optimised frequencies with the generalised fractional Maxwell model.
Figure 7. Responses of individually optimised frequencies with the generalised fractional Maxwell model.
Fractalfract 06 00720 g007
Figure 8. Reduced and rearranged optimised fully fractional generalised Maxwell model with two additional branches.
Figure 8. Reduced and rearranged optimised fully fractional generalised Maxwell model with two additional branches.
Fractalfract 06 00720 g008
Table 1. Nomenclature table with SI units.
Table 1. Nomenclature table with SI units.
SymbolUnitMeaning
A m 2 Area.
α i , 1 ,   α i , 2 /Order or fractional derivation of the springpots.
a , t s Time limits of the fractional derivation.
c i , 1 ,   c i , 2 N s α i , 1 / m ,   N s α i , 2 / m Material properties of the springpots.
E kg / ( ms 2 ) Modul of elasticity.
ε o F / m Absolute permittivity.
ε r /Relative permittivity.
ε 1Strain.
F M a x N Maxwell force.
F e l N Electrical force.
F i N Force in individual branch.
f i mm Calculated data.
i /Current number of fractional Maxwell element.
k N / m Spring constant.
l 1 , l 2 , l 3 m Dimensions of the DEA.
L m Initial length.
l m Displacement.
m /Integer order of derivation by the definition.
m kg Mass of weight.
n /Number of fractional Maxwell elements.
η Ns / m 2 Viscosity.
p /Fractional order of derivation by the definition.
R 2 /Coefficient of determination.
R m e a n 2 /Mean value of coefficient of determination.
s /Laplace operator.
S S t o t /Total sum of squares.
S S r e s /Residual sum of squares.
σ N / m 2 Stress.
V V Voltage.
x i m Displacement of individual branch.
y i mm Measured data.
y ¯ mm Averaged measured data.
Table 2. R m e a n 2 value of the fully fractional generalised Maxwell model for one to five branches.
Table 2. R m e a n 2 value of the fully fractional generalised Maxwell model for one to five branches.
Fully Fractional Generalised Maxwell Model
Number of Branches
R 2 m e a n
n= 10.5456
n= 20.5489
n = 30.5456
n = 40.5456
n = 50.5456
Table 3. Initial and optimised parameters of the fully fractional generalised Maxwell model for one to five branches.
Table 3. Initial and optimised parameters of the fully fractional generalised Maxwell model for one to five branches.
Param.n= 1 α 1 , 1 α 2 , 1 α 2 , 2
Initial 0.211
Optimised 0.20.0021
Parameters c 1 , 1 c 2 , 1 c 2 , 2
Initial 500500500
Optimised 62.9520.0360.142
Param.n= 2 α 1 , 1 α 2 , 1 α 2 , 2 α 3 , 1 α 3 , 2
Initial 0.21111
Optimised 0.2110.5230.046
Param. c 1 , 1 c 2 , 1 c 2 , 2 c 3 , 1 c 3 , 2
Initial 500500500500500
Optimised 62.74047.958173.21430.3810.215
Param.n= 3 α 1 , 1 α 2 , 1 α 2 , 2 α 3 , 1 α 3 , 2 α 4 , 1 α 4 , 2
Initial 0.2111111
Optimised 0.21110.00211
Param. c 1 , 1 c 2 , 1 c 2 , 2 c 3 , 1 c 3 , 2 c 4 , 1 c 4 , 2
Initial 500500500500500500500
Optimised 62.9500.2670.4640.1400.0880.0360.237
Param.n= 4 α 1 , 1 α 2 , 1 α 2 , 2 α 3 , 1 α 3 , 2 α 4 , 1 α 4 , 2 α 5 , 1 α 5 , 2
Initial 0.211111111
Optimised 0.20.00210.002110.00210.002
Param. c 1 , 1 c 2 , 1 c 2 , 2 c 3 , 1 c 3 , 2 c 4 , 1 c 4 , 2 c 5 , 1 c 5 , 2
Initial 500500500500500500500500500
Optimised 62.9500.1030.0360.3210.0360.0940.1570.0360.225
Param.n= 5 α 1 , 1 α 2 , 1 α 2 , 2 α 3 , 1 α 3 , 2 α 4 , 1 α 4 , 2 α 5 , 1 α 5 , 2 α 6 , 1 α 6 , 2
Initial 0.21111111111
Optimised 0.20.00211110.00210.00211
Param. c 1 , 1 c 2 , 1 c 2 , 2 c 3 , 1 c 3 , 2 c 4 , 1 c 4 , 2 c 5 , 1 c 5 , 2 c 6 , 1 c 6 , 2
Initial 500500500500500500500500500500500
Optimised 62.9500.0970.3150.0360.3560.0360.0970.0360.0950.0360.285
Table 4. Initial and optimised parameters of the fully fractional generalised Maxwell model for each individual frequency.
Table 4. Initial and optimised parameters of the fully fractional generalised Maxwell model for each individual frequency.
ParametersF = 1/13 Hz
n = 3
α 1 , 1 α 2 , 1 α 2 , 2 α 3 , 1 α 3 , 2
Initial 0.21111
Optimised 0.17951111
Parameters c 1 , 1 c 2 , 1 c 2 , 2 c 3 , 1 c 3 , 2 R 2
Initial 5005005005005000.658
Optimised 51.7730.0020.0020.0020.002
ParametersF = 1/7 Hz
n = 3
α 1 , 1 α 2 , 1 α 2 , 2 α 3 , 1 α 3 , 2
Initial 0.21111
Optimised 0.1880.2550.0460.2250.880
Parameters c 1 , 1 c 2 , 1 c 2 , 2 c 3 , 1 c 3 , 2 R 2
Initial 5005005005005000.907
Optimised 57.696569.3460.003406.7440.479
ParametersF = 1/5 Hz
n = 3
α 1 , 1 α 2 , 1 α 2 , 2 α 3 , 1 α 3 , 2
Initial 0.21111
Optimised 0.21110.002
Parameters c 1 , 1 c 2 , 1 c 2 , 2 c 3 , 1 c 3 , 2 R 2
Initial 5005005005005000.665
Optimised 73.13319.199188.074305.285191.367
R 2 m e a n 0.743
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Karner, T.; Belšak, R.; Gotlih, J. Using a Fully Fractional Generalised Maxwell Model for Describing the Time Dependent Sinusoidal Creep of a Dielectric Elastomer Actuator. Fractal Fract. 2022, 6, 720. https://doi.org/10.3390/fractalfract6120720

AMA Style

Karner T, Belšak R, Gotlih J. Using a Fully Fractional Generalised Maxwell Model for Describing the Time Dependent Sinusoidal Creep of a Dielectric Elastomer Actuator. Fractal and Fractional. 2022; 6(12):720. https://doi.org/10.3390/fractalfract6120720

Chicago/Turabian Style

Karner, Timi, Rok Belšak, and Janez Gotlih. 2022. "Using a Fully Fractional Generalised Maxwell Model for Describing the Time Dependent Sinusoidal Creep of a Dielectric Elastomer Actuator" Fractal and Fractional 6, no. 12: 720. https://doi.org/10.3390/fractalfract6120720

APA Style

Karner, T., Belšak, R., & Gotlih, J. (2022). Using a Fully Fractional Generalised Maxwell Model for Describing the Time Dependent Sinusoidal Creep of a Dielectric Elastomer Actuator. Fractal and Fractional, 6(12), 720. https://doi.org/10.3390/fractalfract6120720

Article Metrics

Back to TopTop