Next Article in Journal
A New and Versatile Adjustable Rigidity Actuator with Add-on Locking Mechanism (ARES-XL)
Next Article in Special Issue
Investigating the Effect of a Mechanism Combined with a Speed-Increasing Gear and a Pneumatic Artificial Muscle
Previous Article in Journal
Efficient Structure-Based Models for the McKibben Contraction Pneumatic Muscle Actuator: The Full Description of the Behaviour of the Contraction PMA
Previous Article in Special Issue
Pneumatic Multi-Pocket Elastomer Actuators for Metacarpophalangeal Joint Flexion and Abduction-Adduction
 
 
Erratum published on 13 March 2018, see Actuators 2018, 7(1), 9.
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling the Static Force of a Festo Pneumatic Muscle Actuator: A New Approach and a Comparison to Existing Models

1
Department of Electrical Engineering and Computer Science, Control Systems Group, Technical University of Berlin, Einsteinufer 17, D-10587 Berlin, Germany
2
Department VII, Electrical Engineering, Mechatronics and Optometry, Beuth Hochschule Berlin University of Applied Science, Luxemburger Straße 10, D-13353 Berlin, Germany
*
Author to whom correspondence should be addressed.
Actuators 2017, 6(4), 33; https://doi.org/10.3390/act6040033
Submission received: 7 September 2017 / Revised: 15 October 2017 / Accepted: 20 October 2017 / Published: 2 November 2017
(This article belongs to the Special Issue Pneumatic Actuators)

Abstract

:
In this paper, a new approach for modeling the static force characteristic of Festo pneumatic muscle actuators (PMAs) will be presented. The model is physically motivated and therefore gives a deeper understanding of the Festo PMA. After introducing the new model, it will be validated through a comparison to a measured force map of a Festo DMSP-10-250 and a DMSP-20-300, respectively. It will be shown that the error between the new model and the measured data is below 4.4% for the DMSP-10-250 and below 2.35% for the DMSP-20-300. In addition, the quality of the presented model will be compared to the quality of existing models by comparing the maximum error. It can be seen that the newly introduced model is closer to the measured force characteristic of a Festo PMA than any existing model.

1. Introduction

In recent years, pneumatic muscle actuators (PMAs, Figure 1) were integrated in many robotic systems [1,2,3,4] and are still particularly favored, especially in applications that have to be lightweight and powerful.
As can be seen in [6], most robotic systems driven by PMAs, but also driven by other actuators, require an underlying torque controller. While integrating PMAs in robotic systems, the control of PMAs and especially the precise torque control of PMA-driven joints consequently becomes an essential feature. The modeling of PMAs, strongly coupled to model-based controller designs, is therefore still a much discussed topic.
On the one hand, PMAs have their own dynamics, like hysteresis [7], some thermodynamic effects [8,9] and can even be interpreted as a spring–damper combination [10]. On the other hand, it has been shown in [11,12,13] that controlling PMAs, with only a static force map approach, can lead to accurate and high-performance controller designs. The reason is that the PMA dynamics are dominated by the fluid dynamics of the air inside the PMA and the stream of mass characterizing the valve that controls the PMA [14,15]. Nevertheless, a model describing the static PMA force precisely is crucial for the accuracy of the torque controller. Finding the most accurate model to describe the static force characteristic is therefore one goal of this paper.
Different approaches for modeling the static force characteristic of PMAs can be found in literature. The most popular approach is based only on air compression [16]. Although this approach is valid for McKibben PMAs, it is shown by measurements in [14] that this approach does not hold for Festo PMAs. To improve the accuracy of the model presented in [16], the model has been extended in many different ways [8,9,17,18]. It has been shown that these extensions are leading to fewer errors between measured static force maps and the forces predicted by the models. Another approach for the specific use on Festo PMAs has been presented in [12]. The basic idea is the approximation of a Festo PMA as a piston with a virtual, pressure-dependent piston area and a spring that counteracts the expansion. A second approach for the specific use on Festo PMAs is presented in [19], where the static force characteristic is supposed to behave like a mechanical spring with variable stiffness. Because of uncertain parameters, all models, except the first presented in [16], have in common that they must be identified by minimization of the error between a measured static force map and the force defined by the model. By comparing the maximum error, it will be shown in this paper how the accuracy of existing models varies. Furthermore, a new approach for modeling the static force characteristic of Festo PMAs will be presented. The model is physically motivated and therefore gives a deeper understanding of the Festo PMA. After introducing the new model in Section 3 and an explicit discussion of existing models in Section 4.1, all models will be compared to a measured force map of a Festo DMSP-10-250 (DMSP-<initial inner diameter in millimeter>-<initial length in millimeter>) and a DMSP-20-300 in Section 4.3. It will be shown that the error between the new model and the measured data is below 4.4% for the DMSP-10-250 and below 2.35% for the DMSP-20-300. In addition, the quality of the presented model will be compared to the quality of existing models by comparing the maximum error.

2. The Static Force Characteristic

A PMA is characterized by its force map. The PMA force F PMA is a function of the muscle pressure p and its length L, as can be found in [5,17,20,21]. The force map (Figure 2) can be read in the following way: Starting with an unfastened PMA at initial conditions, which means initial length and atmospheric pressure, the PMA does not exert any force. Now putting pressure inside, the PMA gets shorter until it is fully contracted (about 25% [5]). The exerted force is still zero. While pulling the pressurized PMA back to initial length, the PMA will react with the maximum force. The exerted force can be varied between zero and a length-dependent upper limit by varying the PMA pressure [13].

3. A New Model of the Static Force Characteristic

PMAs are a combination of a flexible tube and two stiff aluminium connectors (see Figure 1). The membrane that the tube is made of is a combination of a stiff aramid fiber mesh and a flexible rubber that encloses the air inside the PMA. While putting pressure inside the PMA, the membrane expands and, due to the stiff aramid fibers, the PMA gets shorter. According to [17], the approach for the presented PMA model also is that the fiber length L Fiber stays constant and therefore only the membrane rubber deforms.

3.1. PMA Volume

Following the ideas of [14,15,17], the PMA volume can be approximated by a cylinder. The volume
V ( L , D ) = π 4 D 2 L
of a cylinder is a function of length L and diameter D. Cutting the PMA membrane open and flattening it to a plain (Figure 3), the dependency of the diameter on the length can be calculated by using the Pythagoras theorem [17]. Supposing that L Fiber is constant, the diameter equation
D ( L ) = L Fiber 2 L 2 n π
holds, whereby
L Fiber = L 0 cos Θ 0 and n = L 0 tan Θ 0 π D 0 .
The index 0 indicates the initial state of the PMA, with atmospheric pressure p 0 , initial length L 0 , inner diameter D 0 , membrane thickness H 0 and fiber angle Θ 0 . In [14], it is proven by measurements that the initial fiber angle is Θ 0 = 28.6 . This is the only value that can be found in literature. The initial membrane thickness, H 0 = 1.8 mm for both the DMSP-10-250 and the DMSP-20-300, can easily be calculated from the given initial inner diameter D 0 and the measurable outer initial diameter of the PMA rubber tube.
Inserting the functional dependency on the diameter back into the approximated volume (1), the PMA volume loses its dependency on the diameter and is only a function of the PMA length:
V ( L ) = L · L Fiber 2 4 π n 2 L 3 4 π n 2 .

3.2. Energy-Based Modeling

While pulling a contracted PMA, the muscle reacts with a force F PMA against the pulling direction. The virtual work of the PMA W PMA is given through
W PMA = F PMA · d L .
Furthermore, the virtual work of the PMA can be separated into two parts. On the one hand, the virtual work W VAE has to be done to change the included air volume. On the other hand, additional virtual work W Elast is necessary to change the potential energy of the elastic membrane rubber. Therefore,
W PMA = W VAE + W Elast
always holds.

3.2.1. Virtual Work of the Changing Air Volume

The virtual work needed to change the included air volume inside is given by
W VAE = p · d V .

3.2.2. Elastic Energy of the Membrane

While putting pressure inside the PMA, the actuator gets shorter and expands. The deformation of the membrane is a plane state strain and can be described by the strain in direction of the PMA length ε L and the PMA perimeter ε PE . Rotating the coordinate system by the membrane fiber angle Θ , the deformation is given by the strain in fiber direction ε AF and the direction of pure rubber ε RU that is perpendicular to the fiber direction (see Figure 4). Following the approach of a constant fiber length, the strain ε AF is always zero. This means that the PMA membrane only expands perpendicularly to the fibers.
For this paper, a one-dimensional state of stress is assumed and, therefore, according to Hooke’s law, the tension inside the rubber σ RU will be approximated by
σ RU = E RU ( L ) · ε RU ,
where the modulus of elasticity is supposed to be a function of the PMA length.
Following the approach of a one-dimensional membrane deformation, the strain
ε RU = L L 0 L 0 2 + D D 0 D 0 2
is given by the Pythagoras theorem.
Rotating the tension σ RU back to initial coordinates, it is possible to calculate the tension in length direction
σ L = σ RU · sin θ = E RU ( L ) · ε RU · sin θ = E RU ( L ) · L L 0 L 0
and the tension in perimeter direction of the PMA
σ PE = σ RU · cos θ = E RU ( L ) · ε RU · cos θ = E RU ( L ) · D D 0 D 0 .
Multiplying the tension with the edge surfaces of the unreeled membrane, the virtual work to deform the membrane in length direction is given through
W Elast L = σ L · H 0 · π · D ( L ) = F L · d L .
This approach neglects any effect of lateral contraction. The virtual work that is necessary to deform the membrane in perimeter direction can be calculated in an analog way:
W Elast PE = σ PE · H 0 · L · π = F PE · d D .
The negative sign in Equation (13) is necessary because an increase in length of the PMA has to result in a decreasing perimeter, and this is defined to be positive for this paper.

3.2.3. Summarizing the Energy

Inserting Equations (5), (7), (12) and (13) in Equation (6), the force of the PMA is given by
F Martens ( p , L ) = F PMA ( p , L ) = p · d V d L + F PE · d D d L F L .
It must be noted that Equation (14) is positive, in case the PMA exerts a pulling force.
Because the modulus of elasticity of the membrane rubber and the fiber angle cannot be measured directly, both terms are identified by minimization. The modulus of elasticity E RU used in Equations (10) and (11) is approximated by a polynomial of the third order:
E RU ( L ) = c 3 L 3 + c 2 L 2 + c 1 L + c 0 .
This order is chosen by experiment. It can be seen that higher order polynomials do not lead to better results. However, for lower order polynomials, the calculated error, defined later in this paper (Equation (23)), is much higher. Furthermore, the approximately known fiber angle Θ 0 is corrected by a constant d 0 :
Θ 0 corr = Θ 0 + d 0 .
Θ 0 corr is now used as the new, corrected fiber angle. In this paper, all parameters c j and d j ( j N 0 ) have to be identified through solving the minimization problem (22). The parameters c j and d j are calculated in such a way that the quadratic error between a measured force map and the forces predicted by the models is the smallest. A more detailed discussion of the optimization will be given later in this paper.

4. Validation and Comparison to Existing Models

4.1. Existing Models

Defining a valid model describing the static force characteristic of a PMA is still a much discussed topic. The first model that was introduced to describe a McKibben Muscle [16] was only based on the energy that is needed to change the inner PMA air volume. Effects of elasticity of the membrane material were fully neglected:
F McKibben ( p , L ) = p · d V d L = p · L Fiber 2 3 cos 2 Θ 1 4 π n 2 .
It is shown in [14] that the remaining error between the measured static force characteristic of a Festo PMA and (17) is not negligible. The reason for this is that the Festo PMA stores potential energy in its deformed membrane and the McKibben PMA does not.
A modified version of Equation (17) is presented in [8,9], where Equation (17) has been corrected by two additional factors c 0 and c 1 :
F Andrikopoulos ( p , L ) = p · c 0 · π D 0 2 4 3 tan 2 Θ 0 1 c 1 · L 0 L L 0 2 1 sin 2 Θ 0 .
Although the parameters c 0 and c 1 were calculated directly in [8,9], in this paper, both parameters will be identified by solving the minimization (22). A third improved, inspired by Equation (17) model variant is presented in [18]:
F Sarosi ( p , L ) = c 0 · p + c 1 e c 2 · ( L 0 L ) L 0 + p · d 0 · L 0 L L 0 + d 1 + d 2 .
A model for the specific use on a Festo PMA is presented in [12]. This model is a combination of the pressure-virtual-piston-area product p · A ( L ) and a length-dependent counter force F c ( L ) . The idea is that the PMA behaves like a combination of a pneumatic piston with a variable piston area and a mechanical spring that counteracts the expansion of the PMA:
F Hildebrandt ( p , L ) = p · A ( L ) F c ( L ) = p · j = 0 2 c j L j j = 0 3 d j L j + d 4 L 2 / 3 .
In [19], the static force characteristic of a Festo PMA is supposed to be equivalent to a mechanical spring, with a displacement-pressure-dependent spring stiffness k ( p , Δ L ) :
F Wickramatunge ( p , L ) = c 3 p 2 + c 2 p Δ L + c 1 Δ L 2 + c 0 = k ( p , Δ L ) · Δ L , with Δ L = L min ( L ) .

4.2. Test Rig

The static force maps of a Festo DMSP-10-250 and a DMSP-20-300 have been measured with the test rig depicted in Figure 5.
During the measurement, all screws are fastened and the PMA length is fixed. The pressure is varied between the length-dependent lower limit and the maximum pressure. While varying the pressure, the PMA force is measured by a 1D-forcesensor KD9363-0.5t and a GSV-1A measurement amplifier, both from ME-Meßsysteme GmbH (Hennigsdorf, Germany). Combining the pressure-dependent force characteristics for different lengths, it is easy to determine the static force characteristic as shown in Figure 2. To reduce the influence of measurement errors, the measurement process has been repeated ten times. The depicted force map therefore shows the mean values. Both measured force maps, for the Festo DMSP-10-250 and the DMSP-20-300, can be found in the appendix in Table A1 and Table A2, respectively.

4.3. Results

The measured static force maps are characterizing the force behavior of the Festo DMSP-10-250 and DMSP-20-300, respectively. First, the presented PMA models (14), (18)–(21) are identified by minimizing the quadratic error between the measured force map and the force map calculated by the model. The optimization has been solved by using the MATLAB (R2015b, The MathWorks, Inc., Natick, MA, USA) function fminsearch:
min p i L j F Measurement ( p i , L j ) F Model ( p i , L j ) 2 i = Number of pressure points j = Number of length points .
Because the optimization problem (22) is nonlinear, only a local, start point-dependent minimum can be found. The start point for all models was identified by iterative testing. It can be seen that good results—in the sense of a small error—can be achieved if the start points are chosen with respect to the physical meaning and within a proper range. According to the SI units— 1 m for lengths and 1 Pa = 1 × 10 5 bar for pressure—any parameter c j and d j ( j N 0 ) is set to 1 if it is only multiplied to a length-dependent factor and to 1 × 10 5 if it is multiplied by the pressure. The modulus of elasticity of rubber, necessary for Equation (14), is supposed to be in the area of 1 MPa , and therefore the chosen start point is set to 1 × 10 6 . The initial fiber angle Θ 0 in Equation (14) is supposed to be correct. The chosen parameter estimation’s start points for all models are given without units in Table 1.
To demonstrate the validity of the the estimated parameters c 0 3 and d 0 defining the new model (14), they are given in Table 2. As assumed, the values for the modulus of elasticity are in the range of MPa = 1 × 10 6 Pa , and the initial fiber angle is only slightly corrected by 5.73 for the DMSP-10-250 and 4.18 for the DMSP-20-300, respectively.
All identified models can now be used for calculating their own specific force maps. The error between the calculated and the measured force map
e r r = 100 · m a x ( | F Measurement F Model | ) m a x ( F Measurement )
is defined as the maximum difference between a measured and a calculated force point, normalized to the maximum measured force. Multiplied by 100, the error is given by percentage and is shown in Table 3.
It can be seen in Table 3 that the presented model is closer to the measured static force map of the Festo DMSP-10-250 and DMSP-20-300 than any existing model. The error of Equation (14) for the Festo DMSP-10-250 is smaller than 4.4% and, for the Festo DMSP-20-300, smaller than 2.35%. Furthermore, it can be seen that Equation (17) leads to a maximum error for both PMAs and does not seem to be accurate enough to describe the static force characteristic of a Festo PMA.

5. Conclusions

In this paper, a new approach for modeling the static force characteristic of Festo PMAs is presented. After a detailed derivation of the new model (14), the validity of the model is demonstrated by a comparison to a measured static force map of a Festo DMSP-10-250 and DMSP-20-300, respectively. It is shown that the maximum error between the„true“ force and the one predicted with Equation (14) is smaller than 4.4% for the Festo DMSP-10-250 and smaller than 2.35% for the Festo DMSP-20-300. Furthermore, it is shown that the presented model is closer to the measured data than any existing model that can be found in literature.

Acknowledgments

All presented research results have been supported by the German Federal Ministry for Economic Affairs and Energy (BMWi), especially The Central Innovation Programme for SMEs and the research project Development of a PMA-driven Exoskeleton for the Upper Body (project number ZF4007503). Furthermore, special thanks are extended to Michael Drummond, who developed the excellent test rig that has been used for measuring all static force maps.

Author Contributions

Mirco Martens was the main author of the manuscript and developed the presented model of a Festo pneumatic muscle actuator. Furthermore, Mirco Martens investigated the new model and all existing models under the aspect of accuracy by comparing them to the measured static force maps. Ivo Boblan was the main supervisor and adviser of the work and contributed to the manuscript by reviewing and revising its contents.

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

Appendix A. Measurement Data

Appendix A.1. Force Map of DMSP-20-300

Table A1. Measured force map of a Festo DMSP-20-300 in [ N ] .
Table A1. Measured force map of a Festo DMSP-20-300 in [ N ] .
Length [m]0.2960.2820.2670.250.2370.222
Pressure [Pa]
35,000000000
60,00054.900000
85,000113.400000
110,00017900000
135,000241.800000
146,250273.321.10000
157,500300.242.90000
168,75033170.30000
180,000355.693.70000
197,500402.513127.7000
215,000446.5165.261.6000
232,500495.9204.192.7000
250,000543.4244.7122000
281,250623.5309.7170.242.800
312,500708.8380.422180.600
343,750787.9443.7277.4118.700
375,000870.9518.9328157.600
450,0001065.2676.4462.3252.165.50
525,0001268.9837.5587.8352.3124.70
600,0001460.11004714.8443.5188.10
675,0001653.51164.3840.7535.8253.30

Appendix A.2. Force Map of DMSP-10-250

Table A2. Measured force map of a Festo DMSP-10-250 in [ N ] .
Table A2. Measured force map of a Festo DMSP-10-250 in [ N ] .
Length [m]0.250.2340.2190.2040.198
Pressure [Pa]
500000000
68,75032.16670000
132,50084.83330000
196,250139.50000
260,000195.66670000
285,000217.666717000
310,000239.666735.1667000
335,000261.666752.6667000
360,000283.833371.3333000
402,187321.5101.522.833300
444,375359132.333346.166700
486,562396.6667162.569.333300
528,750434.3333193.666792.166700
559,062461.6667216.3333109.333311.50
589,375488.3333238.1667126.333323.33330
619,688515.8333261142.6667350
650,000542.5282.8333158.6667460

References

  1. Boblan, I.; Maschuw, J.; Engelhardt, D.; Schulz, A.; Schwenk, H.; Bannasch, R.; Rechenberg, I. A human-like robot hand and arm with fluidic muscles: Modelling of a muscle driven joint with an antagonistic setup. In Proceedings of the 3rd International Symposium on Adaptive Motion in Animals and Machines, Ilmenau, Germany, 25–30 September 2005. [Google Scholar]
  2. Caldwell, D.G.; Tsagarakis, N.G.; Kousidou, S.; Costa, N.; Sarakoglou, I. “Soft” exoskeletons for upper and lower body rehabilitation—Design, control and testing. Int. J. Humanoid Robot. 2007, 4, 549–573. [Google Scholar] [CrossRef]
  3. Balasubramanian, S.; Wei, R.; Perez, M.; Shepard, B.; Koeneman, E.; Koeneman, J.; He, J. RUPERT: An exoskeleton robot for assisting rehabilitation of arm functions. In Proceedings of the 2008 Virtual Rehabilitation, Vancouver, BC, Canada, 25–27 August 2008. [Google Scholar]
  4. Noritsugu, T.; Tanaka, T. Application of rubber artificial muscle manipulator as a rehabilitation robot. IEEE/ASME Trans. Mechatron. 1997, 2, 259–267. [Google Scholar] [CrossRef] [Green Version]
  5. Festo AG. Datasheet—Fluidic Muscle DMSP/MAS; Festo AG: Esslingen am Neckar, Germany, 2016. [Google Scholar]
  6. Siciliano, B.; Sciavicco, L.; Villani, L.; Oriolo, G. Motion Control. In Robotics: Modelling, Planning and Control; Springer: London, UK, 2009; pp. 303–361. [Google Scholar]
  7. Lin, C.J.; Lin, C.R.; Yu, S.K.; Chen, C.T. Hysteresis modeling and tracking control for a dual pneumatic artificial muscle system using Prandtl-Ishlinskii model. Mechatronics 2015, 28, 35–45. [Google Scholar] [CrossRef]
  8. Andrikopoulos, G.; Nikolakopoulos, G.; Manesis, S. Incorporation of thermal expansion in static force modeling of pneumatic artificial muscles. In Proceedings of the 2015 23rd Mediterranean Conference on Control and Automation (MED), Torremolinos, Spain, 16–19 June 2015; pp. 414–420. [Google Scholar]
  9. Andrikopoulos, G.; Nikolakopoulos, G.; Manesis, S. Novel Considerations on static force modeling of pneumatic muscle actuators. IEEE/ASME Trans. Mechatron. 2016, 21, 2647–2659. [Google Scholar] [CrossRef]
  10. Kerscher, T.; Albiez, J.; Zollner, J.M.; Dillmann, R. Evaluation of the dynamic model of fluidic muscles using quick-release. In Proceedings of the First IEEE/RAS-EMBS International Conference on Biomedical Robotics and Biomechatronics, Pisa, Italy, 20–22 February 2006; pp. 637–642. [Google Scholar]
  11. Hildebrandt, A.; Sawodny, O.; Neumann, R.; Hartmann, A. A cascaded tracking control concept for pneumatic muscle actuators. In Proceedings of the 2003 European Control Conference (ECC), Cambridge, UK, 1–4 September 2003; pp. 2517–2522. [Google Scholar]
  12. Hildebrandt, A.; Sawodny, O.; Neumann, R.; Hartmann, A. Cascaded control concept of a robot with two degrees of freedom driven by four artificial pneumatic muscle actuators. In Proceedings of the 2005 American Control Conference, Portland, OR, USA, 8–10 June 2005. [Google Scholar]
  13. Martens, M.; Passon, A.; Boblan, I. A sensor-less approach of a torque controller for pneumatic muscle actuator driven joints. In Proceedings of the 2017 3rd International Conference on Control, Automation and Robotics (ICCAR), Nagoya, Japan, 22–24 April 2017; pp. 477–482. [Google Scholar]
  14. Boblan, I. Modellbildung und Regelung Eines Fluidischen Muskelpaares. Ph.D. Thesis, Technical University of Berlin, Berlin, Germany, 2010. [Google Scholar]
  15. Tóthová, M.; Pitel, J. Simulation of actuator dynamics based on geometric model of pneumatic artificial muscle. In Proceedings of the 2013 IEEE 11th International Symposium on Intelligent Systems and Informatics (SISY), Subotica, Serbia, 26–28 September 2013; pp. 233–237. [Google Scholar]
  16. Schulte, H. The Characteristics of the McKibben Artificial Muscle. In The Application of External Power in Prosthetics and Orthotics; National Academy of Sciences—National Research Council: Washington, DC, USA, 1961. [Google Scholar]
  17. Chou, C.P.; Hannaford, B. Static and dynamic characteristics of McKibben pneumatic artificial muscles. In Proceedings of the 1994 IEEE International Conference on Robotics and Automation, San Diego, CA, USA, 8–13 May 1994; pp. 281–286. [Google Scholar]
  18. Sárosi, J.; Biro, I.; Nemeth, J.; Cveticanin, L. Dynamic modeling of a pneumatic muscle actuator with two-direction motion. Mech. Mach. Theory 2015, 85, 25–34. [Google Scholar] [CrossRef]
  19. Wickramatunge, K.C.; Leephakpreeda, T. Study on mechanical behaviors of pneumatic artificial muscle. Int. J. Eng. Sci. 2010, 48, 188–198. [Google Scholar] [CrossRef]
  20. Kelasidi, E.; Andrikopoulos, G.; Nikolakopoulos, G.; Manesis, S. A survey on pneumatic muscle actuators modeling. In Proceedings of the 2011 IEEE International Symposium on Industrial Electronics, Subotica, Serbia, 8–10 September 2011; pp. 1263–1269. [Google Scholar]
  21. Tsagarakis, N.; Caldwell, D.G. Improved modelling and assessment of pneumatic muscle actuators. In Proceedings of the 2000 ICRA IEEE International Conference on Robotics and Automation, San Francisco, CA, USA, 24–28 April 2000. [Google Scholar]
Figure 1. Different Festo pneumatic muscle actuators (PMAs). The depicted connector form is called DMSP. See the data sheet [5] for further information.
Figure 1. Different Festo pneumatic muscle actuators (PMAs). The depicted connector form is called DMSP. See the data sheet [5] for further information.
Actuators 06 00033 g001
Figure 2. Measured static force map of a Festo DMSP-10-250 with an initial length of 250 mm and an initial inner diameter of 10 mm .
Figure 2. Measured static force map of a Festo DMSP-10-250 with an initial length of 250 mm and an initial inner diameter of 10 mm .
Actuators 06 00033 g002
Figure 3. Unreeled membrane of PMA.
Figure 3. Unreeled membrane of PMA.
Actuators 06 00033 g003
Figure 4. Schematic of the PMA membrane and the fiber angle.
Figure 4. Schematic of the PMA membrane and the fiber angle.
Actuators 06 00033 g004
Figure 5. Photo of the test rig that the characteristic force maps have been measured with.
Figure 5. Photo of the test rig that the characteristic force maps have been measured with.
Actuators 06 00033 g005
Table 1. Start points of parameter estimation for all models without units.
Table 1. Start points of parameter estimation for all models without units.
c 0 c 1 c 2 c 3 d 0 d 1 d 2 d 3 d 4
F Andrikopoulos 1 × 10 5 1 × 10 5 -------
F Wickramatunge 11 1 × 10 5 1 × 10 10 -----
F Hildebrandt 1 × 10 5 1 × 10 5 1 × 10 5 -11111
F Sarosi 1 × 10 5 11- 1 × 10 5 1 × 10 5 1--
F Martens 1 × 10 6 1 × 10 6 1 × 10 6 1 × 10 6 0----
Table 2. Estimated parameters for F Martens ; optimization start point is given in Table 1.
Table 2. Estimated parameters for F Martens ; optimization start point is given in Table 1.
c 0 c 1 c 2 c 3 d 0
DMSP-10-250 74.085 MPa 689.20 MPa / m 1.8370 GPa / m 2 848.79 MPa / m 3 5.73
DMSP-20-300 93.232 MPa 715.29 MPa / m 1.5483 GPa / m 2 502.95 MPa / m 3 4.18
Table 3. Force error of each model by percentage.
Table 3. Force error of each model by percentage.
err of F McKibben F Andrikopoulos F Wickramatunge F Hildebrandt F Sarosi F Martens
DMSP-10-25046.1%20.05%13.49%10.12%5.1%4.4%
DMSP-20-30030%13.04%8.2%5.75%3.59%2.35%

Share and Cite

MDPI and ACS Style

Martens, M.; Boblan, I. Modeling the Static Force of a Festo Pneumatic Muscle Actuator: A New Approach and a Comparison to Existing Models. Actuators 2017, 6, 33. https://doi.org/10.3390/act6040033

AMA Style

Martens M, Boblan I. Modeling the Static Force of a Festo Pneumatic Muscle Actuator: A New Approach and a Comparison to Existing Models. Actuators. 2017; 6(4):33. https://doi.org/10.3390/act6040033

Chicago/Turabian Style

Martens, Mirco, and Ivo Boblan. 2017. "Modeling the Static Force of a Festo Pneumatic Muscle Actuator: A New Approach and a Comparison to Existing Models" Actuators 6, no. 4: 33. https://doi.org/10.3390/act6040033

APA Style

Martens, M., & Boblan, I. (2017). Modeling the Static Force of a Festo Pneumatic Muscle Actuator: A New Approach and a Comparison to Existing Models. Actuators, 6(4), 33. https://doi.org/10.3390/act6040033

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop