Next Article in Journal
A Weighted Solution Concept under Replicated Behavior
Next Article in Special Issue
Logarithm-Based Methods for Interpolating Quaternion Time Series
Previous Article in Journal
Design of Observer and Dynamic Output Feedback Control for Fuzzy Networked Systems
Previous Article in Special Issue
Using Regression Analysis for Automated Material Selection in Smart Manufacturing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Machine Learning-Based Parsimonious Models for Static Modeling of Fluidic Muscles in Compliant Mechanisms

by
Monika Trojanová
*,
Alexander Hošovský
and
Tomáš Čakurda
Department of Industrial Engineering and Informatics, Faculty of Manufacturing Technologies with a Seat in Presov, Technical University of Kosice, Bayerova 1, 080 01 Presov, Slovakia
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(1), 149; https://doi.org/10.3390/math11010149
Submission received: 2 November 2022 / Revised: 11 December 2022 / Accepted: 23 December 2022 / Published: 28 December 2022
(This article belongs to the Special Issue Recent Advances of Mathematics in Industrial Engineering)

Abstract

:
This paper uses computational intelligence and machine learning methods to describe experimental modeling performed to approximate the static characteristics of one type of fluidic muscle from the manufacturer FESTO for three different muscle sizes. For the experiments, measured data from the manufacturer and data from a real system (i.e., test device) were used. The measurements, which took place on the experimental equipment, were carried out in two stages (i.e., when the muscle was pressed and when the muscle was relaxed). The resulting measured characteristics were obtained by averaging two values at a given moment. MATLAB® software was used for simulations, in which four models were created: MLP, SVM, ANFIS, and a custom model (i.e., polynomial model). Given that most articles mainly interpret their results graphically when approximating characteristics, in this article, the outputs of the models are also compared with the measured data based on the SSE, NRMSE, SBC, and AIC performance indicators, enabling a more relevant and comprehensive overview of the performance of the individual models. The outputs of the best models described in this article reach an accuracy of 89.90% to 98.74% (all from the MLP group), depending on the muscle size, compared to real measured outputs.

1. Introduction

When designing machines and equipment, it is essential to know the operating parameters of all components and the conditions of the environment in which the equipment is to be applied. Drives are no exception. Conventional types of drive, such as hydraulic actuators or electric motors, have been available for a long time and are commonly used. On the other hand, unconventional types of drive (e.g., pneumatic artificial muscles (PAMs)), although they have been known for several years, continue to be of interest to researchers due to their properties. On the one hand, there are advantages that encourage scientists to try to implement them in practice [1,2], e.g., the possibility of producing greater forces, safe human–machine (robot) interaction, easy replacement in case of failure, etc. On the other hand, their unfavorable properties must be examined and described in detail before they are applied in practice, e.g., hysteresis, friction between fibers, and their static and dynamic characteristics. If a pneumatic artificial muscle is to be used to drive a machine or device, it is necessary to differentiate its parameters depending on the type of muscle. Pneumatic artificial muscles are a general name for a group of unconventional drives whose history dates back to the 1930s when the first pneumatic artificial muscle—the fluid-driven muscle—was invented. The originator of this muscle was Garasiev, who invented braided/netted or embedded membranes, stretching membranes, and rearranging membranes [3]. The aforementioned muscle belongs to one of the subgroups of PAMs: braided muscles. On the other hand, some of the most recent pneumatic artificial muscles are the muscles produced by FESTO or Electro Cam Corporation. A chronological overview of some PAMs that show signs of similarity to braided muscles in their construction or design can be found in Table 1. In addition to this subgroup of muscles, we subsequently distinguish pleated muscles [4,5], netted muscles [6,7,8], and embedded muscles [9,10,11,12]. The last subgroup contains muscles that are the result of design, especially in recent years. Compared to the original braided muscle, these are structures formed from materials such as silicone rubber [13], polyethylene tubes [14], or silicon substrates containing embedded Kevlar fibers [15].
Our object of interest is one type of unconventional drive, namely, pneumatic artificial muscles (specifically, fluidic muscles). A tensile action element formed by a rubber tube and an insert made of aramid fibers is hermetically sealed. After the supply of pressured air to the membrane, an oblong contraction occurs, resulting in a tensile force. The muscle develops the most significant tensile force at the beginning of the contraction; then, as the contraction increases, the tensile force decreases. The areas in which these drives are used are diverse, including clamping, vibrating, positioning, and the use of muscle as a pneumatic spring. Nevertheless, not all possibilities are explored. The reason for this is that, despite their many advantages, the fluidic muscles’ characteristics are nonlinear, and hysteresis is also present. A detailed investigation of the characteristics is essential before a control can be designed to meet the precision and safety requirements of technical equipment.
The static force model of a fluidic muscle is a very important element in the simulation of the behavior of PAM-actuated mechanisms. Most of the models are of a geometric type, relying on the geometric parameters of a muscle to predict the muscle force. Even though they have general characteristics, their reliance on geometric parameters can be problematic under certain conditions (e.g., real-time control) if some parameters are difficult to measure (or even to obtain). In such cases, empirical models may be more suitable, as most of them use a lower number of parameters, although their accuracy may be limited, especially in certain parts of their operational range. Machine learning models have the advantage of being very flexible and can be made quite accurate for the whole range of operating conditions. Nevertheless, their high number of parameters can result in very complex models, which are undesirable in practical scenarios. This paper analyzes the performance of selected machine learning models, with emphasis also being put on the evaluation of the models’ complexity to obtain parsimonious models with a good balance of accuracy and complexity (in terms of the number of model parameters). In addition, the performance of models for different types of fluidic muscles based on their diameter is compared. While the overall shape of the function is comparable for all three types and, therefore, scalable, the measured characteristics actually exhibit differences in output force, and the presumed scalability impairs the resulting model performance if higher accuracy is required.
Therefore, this article describes the results of an experimental study where the goal was to approximate the static characteristics of FESTO fluidic muscles from measured data using computational intelligence and machine learning methods. In total, three different sizes of one muscle type were assessed. The main goal of the research described in this article was to create models that interpret relationships between force, contraction, and pressure approximated by ANFIS, custom (polynomial), MLP, and SVM models for the muscle to the highest possible extent. This article is divided into seven parts, where the Introduction (Section 1) is followed by an overview of related works (Section 2). Subsequently, the object of interest (fluidic muscles) is described in Section 3, where we present the technical data of the selected DMSP muscles (i.e., DMSP-5, DMSP-20, DMSP-40), along with the experimental systems that were made at the authors’ workplace and in which these types of muscles were implemented for the purpose of driving the systems. Section 4 is devoted to describing the models chosen for the modeling phase. The most extensive part of the article is Section 5 (Results), where the resulting models and their outputs are described and graphically interpreted, together with the values of performance indicators for individual muscles. A summary of the results together with their discussion is contained in Section 6, followed by the Conclusion. Since the outputs of the created experimental models were statistically evaluated, Appendix A summarizes the information about four performance metrics (i.e., SSE, NRMSE, AIC, and SBC).

2. Related Works

If a group of braided muscles is considered—such as FESTO, McKibben, or Shadow muscles—the properties of such PAMs differ depending on their shape and the material(s) used. It follows from the above that the force that the muscles can generate is also influenced by several factors (e.g., the pressure applied to the muscle and the contraction ratio). A muscle that is under pressure is highly nonlinear. Therefore, controlling and regulating such a system is difficult. Hence, the effort to solve such problems, e.g., the nonlinearity of the muscle characteristics, which has been a subject of interest of researchers for several decades (e.g., [16,32,33,34]). Modeling and simulations give them space to be able to describe the system as accurately as possible, which is a fundamental prerequisite for further work with the system, e.g., control design, etc.
Analytical modeling is a frequently used approach, since the goal is to find a mathematical model describing the system without the need to obtain real measured data from the system. The source of information for analytical models is knowledge of the physics of ongoing processes and first principles. Experimental modeling is suitable if a real system is available, where it can obtain a representative dataset through measurement. At the same time, this type of modeling is less time-consuming compared to analytical modeling. However, many researchers also exercise the option of combining modeling methods when they choose one of the following within the analytical models [35,36]:
  • Models based on muscle geometry [37,38,39];
  • Models based on muscle properties [40,41];
  • Empirical models [42,43];
  • Phenomenological biomimetic/biomechanical models [44,45,46].
Knowledge of the static characteristics of PAMs is important when designing systems driven in the same way as this type of actuator. Several approaches can be applied to approximate muscle characteristics, e.g., an approximation of the characteristics using a parametric model of the muscle [36], an approximation based on maximum muscle strength [47], or an approximation using an exponential or polynomial function [48,49]. One of the latest articles [50] assessed the effect of the altitude of operation on the pulling force, presenting a numerical approximation of the static characteristics of pneumatic artificial muscles. In contrast, Ref. [51] investigates the static characteristics of pneumatic artificial muscles under isobaric, isotonic, and isometric conditions; at the same time, it presents the dynamic characteristics of the muscle, which interpret the response to a gradual change in pressure under a constant loading force. For testing, the authors developed a laboratory stand described in the article. However, not all models are suitable for approximating the characteristics of FESTO muscles. In [52], the authors state that the existing models of pneumatic artificial muscles are unsuitable for FESTO muscles, as their accuracy is low and, at the same time, they contain too many parameters in their settings. The model presented in the article made it possible to obtain static characteristics with an accuracy of 10%. At the same time, the work also deals with the dynamic modeling of PAM characteristics. The authors of [16] describe methods for determining the static characteristics of pneumatic artificial muscles of the FESTO and Shadow types. These characteristics were monitored under isotonic, isobaric, and isothermal conditions. In the research described in the article, an experimental approach was chosen, where data were obtained by measurements on a test stand and, thus, experimental characteristics were compared with simulated characteristics.
In Table 2, an overview of some related works that were devoted to modeling the characteristics of pneumatic artificial muscles is presented. In addition to the author, the year of publication of the paper, and the citation for the publication itself, the table lists the types of muscle that the authors investigated/described. These were mostly fluidic or McKibben muscles. The table also shows the models that were applied. Most of the publications used real measured data from the system to compare the calculated/simulated outputs of the models. At the same time, the table indicates the type of measurement (i.e., isotonic/isobaric/isometric).

3. The Object of Interest: Fluidic Muscle

To approximate their characteristics, fluidic muscles from the manufacturer FESTO were chosen. The manufacturer offers them in two variants: the MAS muscle type has a screwed connection, while the DMSP muscle type has press-fitted connections. Figure 1 shows the DMSP type of muscle, which falls into the category of single-acting traction drives. The standard sizes of DMSP muscles offered by the manufacturer are 5, 10, 20, and 40. In this study, three muscles of different sizes were considered: DMSP-40 (Figure 1a), DMSP-20 (Figure 1b), and DMSP-5 (Figure 1c). The designation of the muscles by size is derived from their internal diameter. The nominal muscle length is selectable within a specific range for each muscle size. The force also depends on the size of the muscle. An overview of the technical data of the researched muscles is provided in Table 3.
The muscles that were the object of this research were applied at the authors’ workplace in various experimental stands. These devices were assembled to investigate the characteristics and hysteresis of the fluidic muscles, etc.
Figure 2a shows a functional model of a 3-DOF planar robotic arm actuated with DMSP-20 (one pair of FESTO DMSP 20-250N-RM-CM) and DMSP-40 (one pair of FESTO DMSP 40-400N-RM-CM and one pair of FESTO DMSP 20-400N-RM-CM) fluidic muscles. The structure of the manipulator is composed of aluminum profiles, which allow modification of the structure (e.g., in case of the need to shorten/extend one of the arms). Rotational movement is achieved by the antagonistic engagement of three pairs of muscles, joints, and the gear belt subsystem. Electronic pressure regulators (MATRIX ERP50) regulate the pressure in the muscles, which also include pressure sensors. MATLAB/Simulink software is used to control the arm. Control signals are processed with a Humusoft MF624 DAQ PCI card.
A test rig consisting of a multiparallel 1-DOF arm actuated with three pairs of DMSP-5 fluidic muscles is shown in Figure 2b. The rig is made of aluminum profiles and uses six identical muscles with a length of 200 mm (DMSP 5-200N-RM-CM) to actuate a revolute joint (Bosch LF6S drive head with toothed belt and pulley) with a threaded rod arm and a 3 kg payload. The pressure in the muscles is controlled using two proportional pressure regulators with a 0–10 V control signal range (FESTO VPPM-6L-L-1-G18-0L10H-V1P-S1), and pressurized air is generated using a silent compressor (ABAC OS20P-1,5-6RM). The joint angle is sensed using a high-resolution rotary encoder (US Digital E3 with 3600 CPR). All signal processing and generation of control signals is provided by a PC (i7-8750H/8GB) with a DAQ (Humusoft MF644) connected to the computer using a Thunderbolt interface.
A schematic diagram of the test stand that was used for measuring the characteristics of the DMSP-5 muscle is shown in Figure 3. The test stand was designed to measure the muscle’s static characteristics under isobaric conditions (i.e., with constant pressure). The stand consisted of a structure on which the DMSP-5 muscle was attached (both ends of the muscle were firmly attached to the frame of the structure). Compressed air was let into the muscle by the compressor, which was regulated as needed (1–6 bar, with a step of 0.5 bar) using a proportional pressure regulator. With the help of a sliding mechanism, the required contraction of the muscle was ensured, which was measured using a digital height gauge. Therefore, the input value for the measurement was the required muscle contraction (calculated to the length of the muscle) under constant pressure in the muscle. The output monitored variable was the tensile force of the muscle recorded using the force sensor at a given pressure constant.

4. Approximations of the Muscle Characteristics: Theoretical Background

The following well-known and frequently used models were used to approximate the characteristics (i.e., the forces of FESTO fluidic muscles depending on the muscle contraction and pressure in the muscle):
  • MLP (multilayer perceptron network) model;
  • SVM (support-vector machine) model;
  • ANFIS (adaptive neuro-fuzzy inference system) model;
  • Custom (polynomial) model.
The following parts of this chapter describe the basic theoretical premises for the individual models in more detail.

4.1. MLP Model

Since the investigated systems described above represent nonlinear dynamic systems. The first model that was chosen for the approximation of muscle characteristics was the MLP model. The MLP network is one of the most well-known neural network models that allows the classification or approximation of nonlinear functions. The basic task of the network is to transform inputs uj (number of inputs = p) into outputs y with the help of activation functions ϕi—most often hyperbolic tangents (Equation (1)) or logistic functions (Equation (2)) [57,58]:
t a u h   ( x ) = 1 e x p ( 2 x ) 1 + e x p ( 2 x )
l o g i s t i c   ( x ) = 1 1 + e x p ( x )
The architecture consists of a minimum of one hidden layer with a minimum of one neuron called a perceptron (M = the number of perceptrons of the hidden layer), which is subsequently linked to a neuron on the output layer. The hidden layer parameters (hidden layer weight wij) are nonlinear in nature and determine the direction and position of the basis functions. Conversely, the output layer parameters (output layer weights wi) are linear and determine the amplitude of the functions. The output of the MLP network can then be expressed by the following equation [57,58]:
y = i = 0 M w i ϕ i ( j = 0 p w i j u j )
Derivations of the output y of the MLP network according to the weights of the hidden layer wij can be determined by Equation (4), which applies to the weight of the connection between the jth input and the ith neuron of the hidden layer. The function g(x) represents an activation function, e.g., tanh(x). After instituting this form of function into Equation (4), Equation (5) arises. It follows that the activation function ϕi is dependent on two parameters, namely, the inputs uj and the weights of the neurons on the hidden layer wij. The mentioned method represents the backpropagation algorithm for determining the gradients of the MLP network according to the weights [57,58].
y w i j = w i     d g ( x ) d x u j   ;   p r e   u 0 = 1 ,   i = 1 M ,   j = 0 p
y w i j = w i     ( 1 ϕ i 2 ) u j   ;   p r e   u 0 = 1

4.2. SVM Model

The SVM model is a very popular supervised machine learning algorithm used to solve classification and regression tasks, as it offers the best results in many cases. SVM is a method where the goal of the algorithm is to find a hypersurface that divides the data into classes, but with the condition that these classes are as far apart as possible. If the data can be divided by a linear straight line, the model is a linear SVM. However, if it is necessary to use a nonlinear function for data distribution, the model is a nonlinear SVM, which uses a combination of support vectors with a kernel function.
The method of support vectors is based on the essence of calculating the scalar products of two observations (Equation (6)), while Equation (7) applies to the linear classifier of support vectors [59,60]:
x i , x i = j = 1 p x i j x i j
f ( x ) = β 0 + i = 1 n α i x , x i
It follows from these equations that n parameters αi, i = 1, …, n are assumed. To estimate the parameters β0, α1, …, αn, it is necessary to determine the scalar products between all pairs of two observations. However, if these scalar products are replaced by a generalized function that quantifies the similarity of two observations, the linear kernel function (Equation (8)) arises. The task of the kernel function is to transform the data into a higher-dimensional space, where it is possible to find a hypersurface that can divide the data into classes. Within the framework of nonlinear kernel functions, the polynomial kernel (Equation (9)) is known, where d interprets the degree of the polynomial in the form of a whole positive number; along with the radial kernel (Equation (10)), where c is a positive constant [59,60]:
K ( x i , x i ) = j = 1 p x i j x i j
K ( x i , x i ) = ( 1 + j = 1 p x i j x i j ) d
K x i , x i = e x p ( c j = 1 p ( x i j x i j ) 2 )
If a support-vector classifier is combined with a nonlinear kernel, the result is a nonlinear SVM [59,60]:
f ( x ) = β 0 + i ϵ S n α i K ( x , x i )
where S is the set of support vectors.

4.3. ANFIS Model

The third model that was chosen for the approximation of the characteristics was ANFIS. Its name—adaptive neuro-fuzzy inference system—implies that it is a hybrid neuro-fuzzy system, where the fuzzy system is interpreted by a rule of the Takagi–Sugeno type and has the following form:
I f   x 1   i s   A 1 p   a n d   x 2   i s   A 2 p a n d   x n   i s   A n p   t h a n   o p = α 0 p + α 1 p x 1 +   α n p x n
where xi is the input linguistic variable, A i p is the fuzzy set of the ith input and the pth rule, op is the consequent of the pth rule, and α i p represents the parameters of the consequent.
ANFIS consists of five layers, with each layer playing a different role [61,62]:
LAYER 1—Fuzzification layer: O i 1 is the output of node i of the first layer, while each node is adaptive and has a membership function μ A i p ( x ) , for which Equation (13) applies, where a p ,   b p ,   c p is a set of parameters also called premise parameters:
μ A i p ( x i ) = 1 / ( 1 + ( ( x i c p a p ) 2 ) b p )
LAYER 2—Operational layer of T-norm: Each node is fixed and corresponds to one rule. The inputs are membership function values combined by T-norm, while the output wp is the result of all incoming signals:
w p = μ A i p ( x 1 ) μ A j p ( x 2 )                     i = 1 , 2 ;   j = 1 , 2 ;   p = 1 , , 4
LAYER 3—Normalization layer: The task is to determine the relative degree of fulfillment of each rule Rp, as follows:
w p ¯ = w p p = 1 4 w p
LAYER 4—Layer of consequences: Within this layer, the nodes are connected to exactly one node of the normalization layer. The task of the nodes is to calculate the output of the rule in the form of Equation (16), where the parameters α i p denote consequents:
O p 4 = w p ¯ o p = w p ¯ ( α 0 p + α 1 p x 1 + α 2 p x 2 )
LAYER 5—Aggregation layer: The task is to determine the final value of all inputs from the consequent layer:
O 5 = p = 1 4 w p ¯ ( α 0 p + α 1 p x 1 + α 2 p x 2 )  

5. Results

This chapter presents the results of the approximation of the static characteristics of FESTO muscles using the selected models of computational intelligence and machine learning presented in Section 4. Each of the models was created using one of the toolboxes of MATLAB version 2021b and assessed based on four performance metrics, as presented in Appendix A. In addition to the actual approximation results (also interpreted graphically), the used MATLAB toolboxes are briefly described at the beginning of the chapter, along with the essence of configuring the parameters of the individual models. In addition to using the toolboxes, several commands were used—especially for identifying the values of performance indicators of individual models and the graphical interpretation of model outputs, since individual toolboxes do not provide the same interpretation of outputs.
Depending on the muscle size, different datasets were available, which were subsequently used in creating, training, and testing the models. For the DMSP-20 and DMSP-40 muscles, the data were read from the static characteristics provided by the manufacturer. For the DMSP-40 muscle, the tensile strength of the muscle was read over the entire contraction range (not just the operating range). There was a contraction interval of −5% to 30% with a contraction step of 0.5% for individual pressures (0–6 bar). A total of 497 samples were thus obtained. For the DMSP-20 muscle, the tensile force data were read in the same way for individual pressure values of 0–6 bar, but the contraction interval was from −5% to 27% in steps of 0.5%. A total of 455 samples were obtained. However, the dataset for the DMSP-5 muscle was not read from the manufacturer’s specifications; rather, the data were measured using an experimental stand (Figure 3). Through the measurements, characteristic curves were obtained for pressures in the interval of 1–6 bar with a step of 0.5 bar for contractions in the interval of 0 to 20% with a step of 1%. Measurements were carried out in two stages (i.e., inflating and deflating the muscle). A pair of values was thus obtained for each contraction and pressure, and they were subsequently averaged. There were 11 characteristic curves with a total of 497 samples. The total number of samples for individual muscle dimensions was subsequently redistributed at a ratio of 70/15/15% for the training/validation/testing phases, respectively. The total number of samples and the sampling ratios for individual muscles are summarized in Table 4. Normalization of the data was performed implicitly before further processing within the models.

5.1. Neural Network Fitting and Configuration of the Created MLP Models

The MLP models of the neural network were created in the Neural Network Fitting toolbox, which is included in the machine learning and deep learning toolbox group within the MATLAB program. The advantage of this toolbox is that after importing the data and configuring the basic parameters of the network, it is possible to monitor the architecture of the network, but also to perform training/validation/testing without having to use another MATLAB environment. At the same time, the toolbox provides several options for displaying plots (for example, we monitored error histograms and regression plots).
For each muscle type (i.e., DMSP-5, DMSP-20, and DMSP-40), we created five MLP models according to the number of neurons in the hidden layer (1, 2, 5, 10, or 20). The layer size of these networks was in count 1. The number of parameters was also monitored for the individual created models. The input for the network was a two-row matrix with N samples, where each row represented one input parameter—namely, muscle pressure p and muscle contraction h. The output was a one-row matrix with N samples of the tensile force of muscle F. As a training algorithm, we used the Bayesian regularization algorithm.

5.2. Regression Learner App and Configuration of the Created SVM Models

The SVM models were created in the Regression Learner App which, like the previous toolbox, is classified in MATLAB’s machine learning and deep learning toolbox group. This toolbox serves as a comprehensive tool for the creation and subsequent training and testing of linear and nonlinear regression models (e.g., SVM, trees, neural network models, etc.). If the user is not sure which model is most suitable, they can choose several types of models and monitor their success/error rates based on the RMSE criterion. The results of the models can also be interpreted with different graphs, such as response plots, validation residual plots, or predicted vs. actual plots.
A total of six SVM models were created for each muscle. The individual SVM models differed from one another in the functions used: linear, quadratic, cubic, fine Gaussian, medium Gaussian, and coarse Gaussian. The input of the network was a matrix of size 3xN, where the first and second rows represented the inputs pressure and contraction, while the third row was the output, i.e., the tensile force of the muscle. N represents the number of samples in the given matrix. The toolbox allows the selection of the type of validation: cross-validation, holdout validation, or resubstitution validation. In this study, holdout validation was used. All three phases (i.e., training, validation, and testing) were implemented in the aforementioned toolbox. In addition to performance indicators, it is possible to monitor other parameters of the models in the overviews of the created models, such as kernel scale, bias, etc.

5.3. Neuro-Fuzzy Designer and Configuration of the Created ANFIS Models

The ANFIS models were created in Neuro-Fuzzy Designer, which belongs to the control system design and analysis toolbox group. This toolbox allows for designing, training, and validating ANFIS models by importing input/output training data. After importing the data, it enables the generation of an FIS structure (grid partition or subclustering options). For the membership function (MF), the toolbox provides a choice of settings. For its inputs, it is possible to select the number and type of the MFs. The MF type is chosen for its outputs. At the same time, it allows the structure of the created model to be displayed with inputs/outputs, inputs/outputs of the membership function, and rules. The structure also shows what logical operations are used. Before the actual training, it is necessary to choose the optimization method (hybrid or backpropagation), epochs, and tolerance. During training, the toolbox displays a training error. To test the models, they were exported to the workspace of the MATLAB environment and tested using the command evalfis.
Since there are several options to choose from when creating the FIS structure and configuring the membership function, a total of eight ANFIS models were created.
  • The first four models, labeled ANFIS 1–4, had the following features: an FIS structure created by grid partition; 4, 6, 8, and 10 input MFs, and Gaussian type; a constant output MF type; a hybrid optimization method; 0 error tolerance; 10 epochs.
  • The other four models, labeled ANFIS 5–8, had the following features: an FIS structure created using grid partition; 4, 6, 8, and 10 input MFs, and Gaussian type; a linear output MF type; a hybrid optimization method; 0 error tolerance; 10 epochs.
As the number of MFs used in the model structure increased, so did the total number of model parameters, as also shown in the overview tables themselves.

5.4. Curve Fitting Toolbox and the Configuration of the Custom Equation Models

Curve Fitting Toolbox is classified in MATLAB’s Math, Statistics, and Optimization section. The toolbox enables the analysis of imported data, finding functions for interpreting the input data (interpolant, lowess, polynomial, or custom equations), plotting curves and surfaces, and evaluating the resulting functions based on several performance metrics (e.g., SSE, R-squared, RMSE, etc.). The toolbox also provides modeling techniques such as smoothing, interpolation, or splines.
To compare the results of the MLP, SVM, and ANFIS models, a simple four-parameter model for the force function of the fluidic muscles was used. The model included the multiplication of two terms: one with the first independent variable (i.e., pressure), and the second with the second-degree polynomial function of the second independent variable (i.e., contraction). The model was named a custom model, taken for purposes of comparison from [63]. The name custom model is related to the Curve Fitting Toolbox used, where the custom equation option was selected in the toolbox when entering the adopted model. This model, commonly applied in the model-based control design of FM-actuated mechanisms mentioned in [63], has the form of Equation (18):
f ( x , y ) = x   ( a 1 y 2 + a 2 y + a 3 ) + a 4
The task of the toolbox was to find the values of the coefficients a 1   a 4 for the specified measured data. Therefore, one custom model was determined for each muscle. Testing of the created models, which each had four parameters, was carried out outside the toolbox environment.

5.5. Approximation of the Characteristics of PAMs: DMSP 40

When using muscles, but also when examining them, it is important to know their working ranges. The tensile force of the muscle F (N) is dependent on the contraction of the muscle, but the dependence is nonlinear and is different for different muscles depending on their size. The manufacturer indicates the permissible working range of the muscles and defines it with the help of limit curves.
Figure 4 shows the dependences of tensile force F (N) on contraction h (%) from the nominal muscle length for the DMSP-40 muscle type at constant pressures in the range from 0 to 6 bar, as measured and interpreted by the manufacturer. In the picture, four boundary curves can be seen, where boundary 1 expresses the minimum theoretical force at the maximum operating pressure (6000 N). Boundary 2 represents the maximum operating pressure, which is a value of 6 bar. Boundary 3 expresses the maximum deformation; for the DMSP-40 muscle, this is 25% of the nominal muscle length. Finally, boundary 4 is the maximum pretensioning, which for the DMSP-40 muscle is −5% of the nominal muscle length. The grey area in the graph expresses the permissible operating range for the given muscle.
By recording the data from the force characteristics of a given type of muscle, a set of input/output data used for creating models in MATLAB was created. Table 5, Table 6, Table 7 and Table 8 present overviews of the parameters of the created individual models, along with performance metrics comparing the individual models from the validation phase:
An overview of the results of the MLP models is presented in Table 5, which compares five MLP models that differ from one another in the number of neurons on the hidden layer (and, in turn, the number of parameters). Based on the NRMSE criterion, the model with 10 neurons (marked as MLP 10) was the best model, with 97.95% agreement with the reference model (provided by the manufacturer). From the theoretical starting points presented in the article, it is obvious that for the AIC and SBC metrics the value when the indicator’s value is ideal is unknown. These metrics become important if several models can be compared with one another. It is true that the lower the value of the AIC or SBC, the better the model. The MLP 10 model—which achieved the best evaluation within the NRMSE criterion and, subsequently, in the SSE—had the lowest values of the AIC and SBC indicators. Although this model does not have the lowest number of parameters, it is also not among the most complex models.
  • The second type of model that was chosen for approximation was the SVM model. The results of validation based on performance indicators are shown in Table 6. It is obvious that none of the created SVM models could achieve the same results as the best MLP models, since the best-rated model based on the NRMSE criterion (as well as other indicators) was a model with a fine Gaussian function (SVM 4), with 75.52% agreement with the measured data.
  • The third type of model was the ANFIS hybrid model. A total of eight ANFIS models were created, which differed in the type of membership function used and the number of function inputs. An overview of the parameters of the created ANFIS models, along with the results of the performance evaluation based on individual indicators after the validation phase for the DMSP-40 muscle, can be seen in Table 7. Based on the NRMSE, the best model was ANFIS 8, with a rate of agreement of 84.09%, with 10 MF inputs and a linear function type. However, this model did not fare the best in the evaluation by the AIC and SBC, as it was the model with the highest number of parameters, which was reflected in the values of these performance indicators. Therefore, the next considered model—with the same type of MF, but a different number of inputs and better AIC and SBC evaluations—could be the ANFIS 6 or ANFIS 5 models. Compared to the best MLP and SVM models, ANFIS ranked between these models in the NRMSE evaluation.
  • Since the best models among the MLP, SVM, and ANFIS models had a higher number of parameters, one more model was chosen for comparison: a model with four parameters, interpreted as a custom model. Table 8 shows the identified parameters of the model a 1   a 4 . At the same time, the model was assessed by performance metrics, where the NRMSE indicator achieved 53.18% compliance with the measured data.

5.6. Summary of the Model Results for DMSP-40 (Graphical Interpretation)

Figure 5 presents a summary of the model results interpreted graphically for the best models of the given groups, where the simulated force is presented as “Output” in correlation with the force measured by FESTO—marked as “Target”—for DMSP-40 within the test phase. The points labeled “Data” are data that were used for the testing phase. The regression line expresses the agreement between the targets and outputs. The more data type values lie on the regression line, the greater the correlation between the outputs and targets. From the results interpreted in Table 5, Table 6, Table 7 and Table 8, it is obvious that within the NRMSE criterion, from best to worst, the first model was the MLP model, followed by the ANFIS, SVM, and custom models. These results can be confirmed by the graphs in Figure 5. The dispersion of values for the custom model in the region of the regression line is considerable, while the dispersion of values for the MLP model is minimal. Within the dispersion of values, it can be assumed that these may be values on the edge or outside the operating range of the muscle. For the custom model, it can be seen that the predicted forces are lower than the measured output. Differences in performance probably result from differences in characteristics in a given area of the operating range.
Relations of force, contraction, and pressure approximated by the ANFIS, custom, MLP, and SVM models for the DMSP-40 model are shown in Figure 6. These are 3D surfaces that were created from the simulated outputs of the muscle tensile force parameter models based on the measured values of muscle contractions at constant pressures. As mentioned above, the best model for the DMSP-40 muscle was the MLP model, whose curves and corresponding surface can be seen in the lower left panel of Figure 6.

5.7. Approximation of the Characteristics of PAMs: DMSP 20

A similar interpretation of the dependence of the tensile force on the muscle contraction also applies to the DMSP-20 muscle, as shown graphically in Figure 7. Just as was the case for DMSP-40 (Figure 4), this is a graphical dependence measured and interpreted by the muscle manufacturer. In the image, it can be seen that the minimum theoretical force at the maximum operating pressure for DMSP-20 is 1500 N (Figure 7; boundary 1). The maximum operating pressure is a value of 6 bar (Figure 7: boundary 2). Moreover, the maximum deformation is the same as for DMSP-40, i.e., 25% of the nominal length of the muscle (Figure 7; boundary 3). However, the maximum pretensioning is somewhat smaller, at −4% of the nominal length of the muscle (Figure 7; boundary 4).
As for DMSP-40, the input/output data from the manufacturer’s operating range for DMSP-20 were read and used to create models, the results of which are interpreted in Table 9, Table 10, Table 11 and Table 12:
  • Table 9 contains the resulting performance indicators for five MLP models that were modeled under the same conditions as the MLP models for the DMSP-40 muscle. Not only within the NRMSE criterion, but also within the SSE, AIC, and SBC indicators, the MLP 10 model with 10 neurons and an NRMSE value of 89.90% agreement with the reference model emerged as the best model.
  • An overview of the parameters and results of the SVM models for the DMSP-20 muscle can be seen in Table 10. Among the total of six models created, the best-rated model based on NRMSE (63.45%) was the SVM 4 model, which interprets the SVM model with the fine Gaussian function. The values of the metrics SSE, AIC, and SBC also had the lowest values compared to the other models. Although this model showed the best results in terms of indicators, compared to the best MLP model it achieved a significantly worse result.
  • The best ANFIS model for the DMSP-20 muscle also achieved a similar evaluation within the NRMSE as achieved by the best SVM. As part of the overview of the results in Table 11, it can be seen that the best model was ANFIS 7—a model with eight MF inputs and a linear MF. This model has a higher number of parameters, which corresponds to worse ratings within the AIC and SBC indicators. However, within the NRMSE criterion, it achieved the best rating (65.74%). An alternative to the use of the ANFIS 7 model could be the ANFIS 5 or ANFIS 6 models, which achieve worse results within the NRMSE but are simpler models in terms of their number of parameters.
  • Similar to the DMSP-40 muscle, for the DMSP-20 muscle the models with a higher number of parameters emerged as the best in the evaluation. However, the custom model created based on Equation (18) achieved a very low value (only 20.52%) in the evaluation based on the NRMSE criterion. The parameters of the custom model, together with the values of other performance metrics, are listed in Table 12.

5.8. Summary of the Model Results for DMSP-20 (Graphical Interpretation)

In Figure 8, it is possible to see the simulated force (Output) in correlation with the force measured by FESTO (Target) for DMSP-20 within the test phase. The points labeled “Data” represent data that were used for the testing phase. The regression line expresses the agreement between the target and output. From the results interpreted in the figure, along with those presented in Table 9, Table 10, Table 11 and Table 12, it is clear that the MLP model performed best and the custom model performed the worst (as for the DMPS−40 muscle) within the NRMSE criterion. As with DMSP-40, it can be assumed that the dispersion values may be values on the edge or outside the operating range of the muscle. However, for the custom model, the predicted forces were significantly lower than the measured output. Again, performance differences probably resulted from characteristics in a given area of the operating range.
The obtained simulated outputs of the DMSP-20 muscle’s tensile force depending on the muscle contraction and the pressure in the muscle, which were approximated by the selected models, are graphically displayed for the best individual models in Figure 9. The individual points of the simulated characteristics (red points) create 3D surfaces, and if they are compared with the characteristics given by the manufacturer (Figure 7), it is clear that the MLP model achieved the best approximation. The custom model resulting from the evaluation of performance indicators, as well as from graphical interpretation, was not suitable for approximating the static characteristics of the DMSP-20 muscle.

5.9. Approximation of the Characteristics of PAMs: DMSP-5 Experimental Measured Data

Figure 10 shows the graphical dependence of tensile force on contraction for the DMSP-5 muscle, as reported by the muscle manufacturer. The theoretical force at the maximum operating pressure is much smaller than that of the previous muscles due to its technical parameters (140 N; designation 1 in Figure 10). The maximum operating pressure is the same as for the previous muscles. The maximum contraction is 5% lower (20% of the nominal length of the muscle), and the maximum pretensioning is significantly lower than that of the DMSP-20 and DMSP-40 muscles (−1% of the nominal muscle length).
However, data from the manufacturer were not used to approximate the characteristics. Experimental measurements were taken to determine the characteristics of the DMSP-5 muscle on the test device (Figure 3) described in Section 3. Since two sets of measurements were taken for each pressure value (1–6 bar, with a step of 0.5 bar; one measurement at 0–20% contraction; a second measurement at 20–0% contraction), these measurements were subsequently processed and adjusted by averaging the values. The results are the measured static characteristics of the DMSP-5 fluidic muscle, as presented in Figure 11. These measured characteristics subsequently served to create and verify the models.
Table 13, Table 14, Table 15 and Table 16 present the interpretation of the resulting models created for the DMSP-5 muscle:
  • Table 13 shows an overview of the parameters and indicators of the MLP models. Compared to the MLP models for DMSP-20 and DMSP-40, based on the NRMSE criterion, the model with 20 neurons (MLP 20) emerged as the best model, with an agreement of 98.74% with the measured data. However, selection based on the SSE, AIC, and SBC metrics would also confirm this selection, as the values of these indicators were the lowest for the model in question.
  • The variability in the results of the SVM models for DMSP-5 is visible in the values of the indicators shown in Table 14 (NRMSE ranging from 36.96% to 81.93%). The best model reached an NRMSE of 81.93% using a cubic function (SVM 3). If the best model was chosen based on the SSE, the choice would be the same. However, a change would occur based on the AIC and SBC, when the SVM 5 model with the medium Gaussian function and an NRMSE value of 80.76% would also be a suitable candidate.
  • The ANFIS models created and tested for the DMSP-5 muscle achieved much better results compared to the ANFIS models for DMSP-20 and DMSP-40. Moreover, in the comparison of the created models, the results of the NRMSE evaluation for the best ANFIS model (NRMSE = 96.31%) were similar to those of the best MLP model (NRMSE = 98.74%). In Table 15, it can be seen that the best ANFIS model was the ANFIS 6 model, which contains 6 MF inputs and a linear MF. The number of parameters in this model is not the lowest, but at the same time, it does not belong to the group of models with the highest number of parameters. This model achieved the best values in almost all indicators.
  • The last approximated function of the static characteristics of the DMSP-5 muscle based on the measured data was a custom model; the transcription of the general equation is identical to Equation (18), and the values of the model parameters are summarized in Table 16. Within the framework of the custom models created for individual muscles, they achieved the highest rating for NRMSE (72.86%). Nevertheless, their results were not comparable with those of the models listed above.

5.10. Summary of the Model Results for DMSP-5 (Graphical Interpretation)

Figure 12 shows an overview of the force simulated (Output) by the ANFIS, custom, MLP, and SVM models in correlation with the force measured (Target) for DMSP-5 in the test phase. From the graphs selected for the best models, it is possible to confirm what was presented in Table 13, Table 14, Table 15 and Table 16, i.e., that the best model for approximating static characteristics based on the NRMSE evaluation criterion is the MLP model. At the same time, the largest dispersion of values around the regression line (and, thus, the worst model) is the custom model.
The static characteristics of the DMSP-5 muscle approximated by the best models, along with their created 3D surfaces, can be seen in Figure 13. From the measured characteristics shown in Figure 11—which were used to create the models and, subsequently, to compare their outputs—it is obvious that the maximum tensile force at a pressure of 6 bar should reach 150 N. However, from the graphs shown in Figure 13, it can be seen that the custom model is not capable of approximating the given value; therefore, the custom model was the worst among the individual evaluation models (Table 13, Table 14, Table 15 and Table 16). On the other hand, the curves of the MLP and ANFIS models are very similar and almost identical to the measured output, as proven by the values of the NRMSE criterion.

6. Summary and Discussion

This final section contains a summary of the results together with their discussion and is divided into two parts: The first part summarizes the results of the approximation models for the individual muscles, always choosing the best from a given group of models. The results as described and graphically interpreted in Figure 14, Figure 15 and Figure 16 compare the models’ outputs with the force measured for the individual muscles. This section is followed by a discussion of the results based on the best models according to their correspondence with the measured data and based on their complexity (i.e., number of parameters). The discussion also summarizes the results of the models of all architectures for each of the modeled muscle types using selected performance metrics, which are interpreted graphically in the form of boxplots.
The most relevant graphical overview of the model outputs for the DMSP-40 muscle is presented in Figure 14, which summarizes all of the outputs of the best models in one graph together with the measured data. The measured tensile force is indicated by a black dashed line labeled F. The outputs of the individual models are shown by solid colored lines:
  • MLP model with 10 neurons and an NRMSE value of 97.95% (red line, marked as YS-mlp);
  • ANFIS model with 10 MF inputs and a linear MF, with an NRMSE value of 84.09% (green line, marked as YS-anfis);
  • SVM model with a fine Gaussian function and an NRMSE value of 75.52% (blue line, marked as YS-svm);
  • Custom model with an NRMSE value of 53.18% (yellow line, marked as YS-custom).
An overview of the model outputs, summarizing all of the outputs of the best models in one graph together with the measured data for the DMSP-20 muscle, is shown in Figure 15. The measured tensile force is indicated by the black dashed line marked F. The outputs of the individual models are shown by solid colored lines:
  • MLP model with 10 neurons and an NRMSE value of 89.90% (red line, marked as YS-mlp);
  • ANFIS model with 8 MF inputs and a linear MF, with an NRMSE value of 65.74% (green line, marked as YS-anfis);
  • SVM model with a fine Gaussian function and an NRMSE value of 63.45% (blue line, marked as YS-svm);
  • Custom model with an NRMSE value of 20.52% (yellow line, marked as YS-custom).
A comparison of the outputs of the best models in one graph, together with the experimentally obtained data for the DMSP-5 muscle, is shown in Figure 16. The measured tensile force is indicated by a black dashed line labeled F. The outputs of the individual models are shown by solid colored lines:
  • MLP model with 20 neurons and an NRMSE value of 98.74% (red line, marked as YS-mlp);
  • ANFIS model with 6 MF inputs and a linear MF, with an NRMSE value of 96.31% (green line, marked as YS-anfis);
  • SVM model with a cubic function and an NRMSE value of 81.93% (blue line, marked as YS-svm);
  • Custom model with an NRMSE value of 72.86% (yellow line, marked as YS-custom).
Taking into account the results of the measured force characteristics of the investigated muscle types (either supplied by the muscle manufacturer or determined during the experiments), it can be observed that the function has a similar form but with different scaling. However, simple scaling of the output of a single model trained with the data for one of the muscle diameters would not be sufficient to obtain a good performance when higher accuracy is required. The 3D plots of the function reveal two important areas that manifest themselves differently when using various types of muscle force function approximators. One of the areas is the area of zero-force values, which increases with the decrease in the value of muscle pressure. The other area includes the area of non-zero muscle force, i.e., the area where a muscle actually performs useful work and generates force. It is noteworthy that in the case of measurements carried out for DMSP-5, the distribution of the data was denser for set pressures because of the lower pressure step (50 kPa, compared to 100 kPa for the characteristics given by the manufacturer). In addition, the mean curves of two directions of length changes were used; therefore, the hysteresis of the force function was not considered.
As shown in the results for model performance, the area of zero force appears to be more demanding in terms of the required flexibility of the model for achieving higher accuracy. As a result of this, the performance of the simple custom model with four parameters is significantly reduced compared to the other models. While the inaccuracies of the muscle force model approximators in the zero-force area caused by the output reaching negative values could be handled using the saturation function with a lower bound of zero, this solution would not help for the output errors above zero. Moreover, the regression plots for test data shown in Figure 12 illustrate that both the SVM and custom models also have inferior performance in the area of higher force. Even though the approximation capabilities of ANFIS and MLP are naturally expected to be more powerful due to their higher numbers of parameters, it is important to obtain parsimonious models where the balance of model complexity and performance is taken into consideration. These can be compared using the AIC and SBC performance metrics which, based on the results, still favor the MLP models.
In Table 17, which summarizes the best results for each type of DMSP fluidic muscle, it is clear that the MLP model prevails as the best within the individual muscles. Nevertheless, let us analyze these results:
  • For the DMSP-40 muscle, the MLP 10 model achieved an AIC of 3678.1 and an SBC of 3850.6. The other models that approached these values were also from the MLP group. The SVM models started at AIC 6330 and SBC 6880 (SVM 4). The ANFIS models started at similar values—AIC 6190 and SBC 6450 (ANFIS 5). The custom model achieved values in these metrics that ranked it among the models with weaker results.
  • A similar scenario occurred for the DMSP-20 muscle, with the best AIC and SBC values of 3624.1 and 3793, respectively, for the 10-neuron MLP. Although the AIC and SBC values of the best model were similar to those for the DMSP-40 muscle, at the same time, the models that approached the best results were also from the group of MLP models, with only one change. The SVM, ANFIS, and custom models achieved lower values in the AIC and SBC parameters than for the DMSP-40 muscle. For example, the best SVM model started at AIC 5000 and SBC 5600 (SVM 4), while the best ANFIS model started at AIC 5160 and SBC 5290 (ANFIS 1). The values of these two metrics for the custom model were around 5400.
  • However, the most fundamental difference in the AIC and SBC metrics was observed for the DMSP-5 muscle. The metric values of the best MLP model, with 20 neurons, can be seen in Table 17. However, interesting results were also achieved by the ANFIS 6 model, with AIC 332 and SBC 787, as well as by the custom model, where the AIC reached 998.3944 and the SBC reached 1012.2. The best SVM model (SVM 5) started with AIC and SBC values of 971.5456 and 5830, respectively.
It is also interesting to generalize the performance of the models using the selected performance metrics for all muscle types in question (Figure 17 and Table 17). Generally, we were interested in the best results for the approximation of muscle force function. The boxplots in Figure 17 allow the capabilities of the selected architectures of the approximators to be evaluated. However, it must be noted that only the NRMSE metric includes normalization and, thus, makes the results comparable between the muscle types. The highest median values of the NRMSE metric, recorded for DMSP-5, were probably caused by the different method of obtaining data, as well as the slightly increased number of data samples. On the other hand, the best-performing function approximators had fit values in the order of 90% (Table 17). The AIC metrics show similar ranges of their 25th and 75th percentiles, but the unnormalized form causes different placement of their respective boxplots on the y-axis. In terms of model parsimony, the AIC and SBC criteria are more important because of their reflection of the complexity/fit balance. The boxplots of the AIC criterion also confirm that several models performed much better compared to others (i.e., low-value outliers). This applied to the MLP models, which dominated the other model architectures for this type of function.
It appears probable that the sigmoid activation function is especially suited for the shape of a DMSP muscle force function and can obtain a good approximation with the best complexity/fit balance, as reflected in the values of AIC. The summary results shown in Table 17 allow the best model approximator to be determined—in all three cases, it is the MLP architecture, with the 10-neuron structure for DMSP-20 and DMSP-40, and with the 20-neuron structure for DMSP-5. However, the difference in NRMSE between MLP-10 and MLP-20 is roughly 1.5%, which makes MLP-10 possibly the best choice for the approximation of DMSP muscles’ force function.
Table 18 shows the comparison of the relevant metrics for the best-performing model in this work (i.e., MLP10 for DMSP-20, which was the only muscle used in other references) with four other models used in other works.
It must be noted that the number of data points was not available and had to be estimated based on the used metrics (it was determined to be 417). It can be observed that the value of RMSE was the lowest for the case of MLP10; of course, as expected in the case of ML models, this was achieved with a higher number of parameters. Nevertheless, the value of AIC, which takes the balance between model accuracy and complexity into account, was still the lowest. This suggests that the high flexibility of ML models might be rewarded with a sufficiently high increase in accuracy to retain the beneficial balance of accuracy/complexity in the resulting models. In consequence, when high accuracy is of interest (e.g., in simulation models used for optimization or verification), the parsimonious ML models can be useful. In addition, the high flexibility of these models allows researchers to vary their structure, include additional inputs easily, and also achieve very high accuracies that other models might not be capable of [66].
The algorithms for training the ML models were sequential minimal optimization for the SVM models, a hybrid method (i.e., LS and backpropagation) for training the ANFIS models, and the Levenberg–Marquardt algorithm for training the MLP models. The SMO algorithm for SVM training is said to reduce the time complexity to O(n) [67]. Hybrid training of ANFIS has a dominant time complexity of O(n2) [68], while LM training of MLP has a time complexity of O(n3) [69]. It should be noted that these are the time complexities for the training of the models, which is expected to be carried out offline if their use in the simulation of PAM-driven mechanisms is assumed.

7. Conclusions

The aim of this research was to approximate the static characteristics of FESTO DMSP-type fluidic muscles using selected methods of computational intelligence and machine learning. The study, which analyzed the results of the experiments described in this article, can be divided into the following stages:
  • Reading of the static characteristics of the dependence of the tensile strength of the DMSP-20 and DMSP-40 muscles upon muscle contraction at constant pressures.
  • Experimental acquisition of a representative dataset on the dependence of the tensile force of the DMSP-5 muscle on the contraction of the muscle at constant pressures on the testing device developed by the authors.
  • Creation of MLP, SVM, ANFIS, and custom models in MATLAB and associated toolboxes for individual types of muscles.
  • Testing of the trained models on an independent dataset and performance evaluation of model outputs compared with real outputs, based on the performance indicators NRMSE, SSE, AIC, and SBC.
Based on the results, which are described in detail in the Results and Discussion section, it is obvious that the functions of the force characteristics of the muscles have a similar form, but with different scaling; therefore, each type of muscle mentioned in this article must be considered separately. Nevertheless, the results can be summarized in the following findings:
  • The flexibility of the created models was lower in the area of zero force (maximum contraction values), which reduced the performance of the models.
  • In general, models from the MLP group achieved the best results in all performance metrics.
  • The sigmoid activation function appears to be a suitable candidate for the muscle force function of the DMSP muscles; not only can it approximate with high values of fit, it also offers a compromise in terms of the complexity-to-fit ratio of the model.
  • Given that the aim of this article was to find models with not only high performance but also a good complexity/fit balance, the best model based on the results appears to be the MLP with 10 neurons for each of the mentioned types of DMSP muscle from the manufacturer FESTO.
Considering the above information, in the future, we would like to improve the best models—whose performance is reduced due to low flexibility in the region of zero force (negative values)—by using a saturation function with a lower limit defined at zero. By increasing the flexibility of models in the field of zero force and, at the same time, using the information and knowledge that we acquired during the processing of this article, we wish to ensure the optimal conditions for the creation of models of the mechanisms described in this article and to which DMSP muscles are applied.

Author Contributions

Conceptualization, A.H. and M.T.; methodology, A.H.; measuring data, A.H.; modeling in MATLAB, M.T.; formal analysis, M.T.; investigation, A.H. and M.T.; writing—original draft preparation, M.T., A.H. and T.Č.; writing—review and editing, M.T., A.H. and T.Č.; visualization, M.T. and T.Č.; project administration, M.T.; funding acquisition, A.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Slovak Research and Development Agency under contract No. APVV-19-0590—Modular Multifunctional Inspection Workplace using Computational Intelligence Techniques; and by the projects VEGA 1/0700/20—Object Recognition Techniques with Convolutional Neural Networks; and KEGA 055TUKE-4/2020—Transfer of Knowledge from Research on Digitization of Production Processes into Study Programs of the Faculty of Manufacturing Technologies, granted by the Ministry of Education, Science, Research, and Sport of the Slovak Republic.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Performance Metrics

To evaluate and compare the results of the models’ validation in the process of approximating the characteristics of fluid muscles, we used four performance indicators:
  • SSE—sum of square error;
  • NRMSE—normalized root-mean-square error;
  • SBC—Schwarz Bayesian information criterion;
  • AIC—Akaike information criterion.
These metrics can be divided into two groups: One group consists of performance indicators that evaluate the error rate of models (i.e., NRMSE and SSE). The second group includes those that are assessed based on the numbers of parameters of the individual models (i.e., AIC and SBC).

Appendix A.1. SSE and NRMSE

The sum of square error (SSE) is a performance indicator that expresses the error as the sum of the squares of the difference between the measured (observed) value yi for the ith sample and the model−predicted response value y ^ i for the ith sample, where i = 1 ,   2 ,   3 , ,   n . The following relationship applies to determine the value of the SSE performance indicator [70,71]:
S S E = i = 1 n ( y i y ^ i ) 2
In general, the lower the value of the SSE, the smaller the variability of the model compared to the measured values.
The normalized root-mean-square error (NRMSE) is a performance indicator that compares the degree of agreement between the output of the model for the ith sample and the measured value (output of the system) for the ith sample, where n is the number of samples. The value of the NRMSE indicator is determined based on the following condition [72,73]:
N R M S E = ( 1 j = 1 n [ y j y ^ j ] 2 j = 1 n [ y j 1 n j = 1 n y j ] 2 ) × 100 %
An NRMSE value of 100% indicates that the model output and the system output are identical.

Appendix A.2. SBC and AIC

The Schwarz Bayesian information criterion (SBC) is a performance indicator that evaluates which model from the collection of created candidate models M k 1   M k L is the most suitable for describing measurements n with a response y based on the available data. The essence is to find the candidate model that is closest to the real model g(y)). If we suppose that y is generated to an unknown density g(y), then for each model M k of the candidate models there is a k-dimensional parametric class of density function [74,75]:
( k ) = { f ( y | θ k ) : θ k Θ ( k ) }
The SBC is based on the empirical log-likelihood. Let L ( θ k | y ) be the designation for the likelihood of the density f ( y | θ k ) and θ ^ k be the designation of a vector of estimates, which was created by maximizing the likelihood L ( θ k | y ) over Θ ( k ) . Subsequently, it is possible to establish the SBC criterion for the candidate model M k as follows [74,75]:
S B C = 2 ln L ( θ ^ k | y ) + k   ln ( n )
In practice, if the values of the SBC for the individual-created models are known, then the model that achieves the lowest SBC value during evaluation is the closest to the real model. In MATLAB, it is possible to set the SBC criterion using the following equation:
S B C = n · log ( S S E / n ) + p · log ( n )
where n is the number of samples, p is the number of parameters of the given created model, and SSE is the sum of square error.
The Akaike information criterion (AIC), like the SBC, is a performance indicator that expresses the quality of the candidate model M k compared to the actual model g(y) with respect to the number of parameters of the given model. The condition for using criterion is that the model must be tested on another dataset. As for the SBC indicator, the following condition applies:
( k ) = { f ( y | θ k ) : θ k Θ ( k ) }
Since L ( θ k | y )   = f ( y | θ k ) and θ ^ k is a vector of estimates created by maximizing the likelihood L ( θ k | y ) over Θ ( k ) , it is possible to establish the AIC for the candidate model M k , as follows [76,77]:
A I C = 2 log f ( y | θ ^ k ) + 2 k  
Equation (A7) can be divided into two parts, which gives us two different statistics: The first is the concept of goodness-of-fit 2 log f ( y | θ ^ k ) , which reflects the agreement of the model. 2k is called the penalty term and is a bias correction. In practice, the AIC is evaluated in the same way as the SBC—the lower the value of the criterion, the closer the simulated model is to the real one. In MATLAB, the AIC can be determined as follows:
A I C = n · log ( S S E / n ) + 2 · p  
where the interpretation of the individual parameters is the same as for the SBC.

References

  1. Kalita, B.; Leonessa, A.; Dwivedy, S.K. A Review on the Development of Pneumatic Artificial Muscle Actuators: Force Model and Application. Actuators 2022, 11, 288. [Google Scholar] [CrossRef]
  2. Craddock, M.; Augustine, E.; Konerman, S.; Shin, M. Biorobotics: An Overview of Recent Innovations in Artificial Muscles. Actuators 2022, 11, 168. [Google Scholar] [CrossRef]
  3. Daerden, F.; Lefeber, D. Pneumatic Artificial Muscles: Actuators for Robotics and Automation. Eur. J. Mech. Environ. Eng. 2002, 47, 11–21. [Google Scholar]
  4. Daerden, F.; Lefeber, D. The Concept and Design of Pleated Pneumatic Artificial Muscles. Int. J. Fluid Power 2001, 2, 41–50. [Google Scholar] [CrossRef] [Green Version]
  5. Terryn, S.; Brancart, J.; Lefeber, D.; Van Assche, G.; Vanderborght, B. A Pneumatic Artificial Muscle Manufactured Out of Self-Healing Polymers That Can Repair Macroscopic Damages. IEEE Robot. Autom. Lett. 2018, 3, 16–21. [Google Scholar] [CrossRef]
  6. Yariott, J.M. Fluid Actuator. U.S. Patent US3645173A, 29 February 1972. [Google Scholar]
  7. Immega, G. Bellows Actuator. U.S. Patent US5181452A, 26 January 1993. [Google Scholar]
  8. Kukolj, M. Axially Contractable Actuator. U.S. Patent US4733603A, 29 March 1988. [Google Scholar]
  9. Morin, H.A. Elastic Diaphragm. U.S. Patent US2642091A, 16 June 1953. [Google Scholar]
  10. Han, K.; Kim, N.-H.; Shin, D. A Novel Soft Pneumatic Artificial Muscle with High-Contraction Ratio. Soft Robot. 2018, 5, 554–566. [Google Scholar] [CrossRef]
  11. Veale, A.J.; Xie, S.Q.; Anderson, I.A. Modeling the Peano Fluidic Muscle and the Effects of Its Material Properties on Its Static and Dynamic Behavior. Smart Mater. Struct. 2016, 25, 065014. [Google Scholar] [CrossRef]
  12. Wirekoh, J.; Valle, L.; Pol, N.; Park, Y.-L. Sensorized, Flat, Pneumatic Artificial Muscle Embedded with Biomimetic Microfluidic Sensors for Proprioceptive Feedback. Soft Robot 2019, 6, 768–777. [Google Scholar] [CrossRef] [Green Version]
  13. Oguntosin, V.; Akindele, A. Design and Characterization of Artificial Muscles from Wedge-like Pneumatic Soft Modules. Sens. Actuators A Phys. 2019, 297, 111523. [Google Scholar] [CrossRef]
  14. Greer, J.D.; Morimoto, T.K.; Okamura, A.M.; Hawkes, E.W. Series Pneumatic Artificial Muscles (SPAMs) and Application to a Soft Continuum Robot. IEEE Int. Conf. Robot Autom. 2017, 2017, 5503–5510. [Google Scholar] [CrossRef]
  15. Wirekoh, J.; Park, Y.-L. Design of Flat Pneumatic Artificial Muscles. Smart Mater. Struct. 2017, 26, 035009. [Google Scholar] [CrossRef]
  16. Takosoglu, J.E.; Laski, P.A.; Blasiak, S.; Bracha, G.; Pietrala, D. Determining the Static Characteristics of Pneumatic Muscles. Meas. Control 2016, 49, 62–71. [Google Scholar] [CrossRef]
  17. Gaylord, R.H. Fluid Actuated Motor System and Stroking Device. U.S. Patent US2844126A, 22 July 1958. [Google Scholar]
  18. Tagami, T.; Miyazaki, T.; Kawase, T.; Kanno, T.; Kawashima, K. Pressure Control of a Pneumatic Artificial Muscle Including Pneumatic Circuit Model. IEEE Access 2020, 8, 60526–60538. [Google Scholar] [CrossRef]
  19. Tondu, B.; Lopez, P. The McKibben Muscle and Its Use in Actuating Robot-arms Showing Similarities with Human Arm Behaviour. Ind. Robot. Int. J. 1997, 24, 432–439. [Google Scholar] [CrossRef]
  20. Jamwal, P.K.; Hussain, S.; Xie, S.Q. Dynamic Modeling of Pneumatic Muscles Using Modified Fuzzy Inference Mechanism. In Proceedings of the 2009 IEEE International Conference on Robotics and Biomimetics (ROBIO), Guilin, China, 19–23 December 2009; pp. 1451–1456. [Google Scholar]
  21. Yamamoto, Y.; Matsunaga, N.; Okajima, H. Robust Variable Stiffness Control of McKibben Type Pneumatic Artificial Muscle Arm by Using Multiple Model Error Compensators. In Proceedings of the 2017 17th International Conference on Control, Automation and Systems (ICCAS), Jeju, Korea, 18–21 October 2017; pp. 957–962. [Google Scholar]
  22. Ku, K.K.; Bradbeer, R. Static Model of the Shadow Muscle under Pneumatic Testing. In Proceedings of the 4th Regional Inter-University Postgraduate Electrical and Electronics Engineering Conference, Macao, China, 13–14 July 2006. [Google Scholar]
  23. Marcinčin, J.; Palko, A. Negative Pressure Artificial Muscle—An Unconventional Drive of Robotic and Handling Systems; Transactions of the University of Košice: Košice, Slovakia, 1993; pp. 350–354. [Google Scholar]
  24. Novák-Marcinčin, J.; Smrček, J. The UPAM-New Actuator for Advanced Robotic Devices. In Proceedings of the 25th International Symposium on Industrial Robots, Hannover, Germany, 20–27 April 1994; pp. 299–304. [Google Scholar]
  25. Novák-Marcinčin, J.; Smrček, J.; Niznik, J. Bioactuators-Most Efficient Actuators for Biomechanics. In Proceedings of the 7th International Conference on Measurement in Clinical Medicine (IMEKO TC-13)—Model Based Biomeasurements, Stara Lesna, Slovak Republic, 6–9 September 1995; pp. 369–371. [Google Scholar]
  26. Bergemann, D.; Lorenz, B.; Thallemer, A. Actuating Means. U.S. Patent US6349746B1, 26 February 2002. [Google Scholar]
  27. Chen, Y.; Sun, N.; Liang, D.; Qin, Y.; Fang, Y. A Neuroadaptive Control Method for Pneumatic Artificial Muscle Systems with Hardware Experiments. Mech. Syst. Signal Process. 2021, 146, 106976. [Google Scholar] [CrossRef]
  28. Sun, N.; Liang, D.; Wu, Y.; Chen, Y.; Qin, Y.; Fang, Y. Adaptive Control for Pneumatic Artificial Muscle Systems With Parametric Uncertainties and Unidirectional Input Constraints. IEEE Trans. Ind. Inform. 2020, 16, 969–979. [Google Scholar] [CrossRef]
  29. Cao, Y.; Huang, J.; Xiong, C. Single-Layer Learning Based Predictive Control with Echo State Network for Pneumatic Muscle Actuators-Driven Exoskeleton. IEEE Trans. Cogn. Dev. Syst. 2020, 13, 80–90. [Google Scholar] [CrossRef]
  30. Davis, D.; Carlson, J. Fluidic Actuator. U.S. Patent US20040107829A1, 10 June 2004. [Google Scholar]
  31. Davis, D.L.; Carlson, J.A. Fluidic Actuator. U.S. Patent US6868773B2, 22 March 2005. [Google Scholar]
  32. Kotkas, L.; Zhurkin, N.; Donskoy, A.; Zharkovskij, A. Design and Mathematical Modeling of a Pneumatic Artificial Muscle-Actuated System for Industrial Manipulators. Machines 2022, 10, 885. [Google Scholar] [CrossRef]
  33. Capace, A.; Cosentino, C.; Amato, F.; Merola, A. A Multistate Friction Model for the Compensation of the Asymmetric Hysteresis in the Mechanical Response of Pneumatic Artificial Muscles. Actuators 2019, 8, 49. [Google Scholar] [CrossRef] [Green Version]
  34. Kalita, B.; Dwivedy, S.K. Nonlinear Dynamic Response of Pneumatic Artificial Muscle: A Theoretical and Experimental Study. Int. J. Non-Linear Mech. 2020, 125, 103544. [Google Scholar] [CrossRef]
  35. 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, Gdansk, Poland, 27–30 June 2011; pp. 1263–1269. [Google Scholar]
  36. Piteľ, J.; Balara, M.; Hošovský, A.; Tóthová, M. Pneumatické Umelé Svaly: Modelovanie, Simulácia, Riadenie; Technická Univerzita v Košiciach: Košice, Slovakia, 2015; Volume 1. [Google Scholar]
  37. Chou, C.-P.; Hannaford, B. Measurement and Modeling of McKibben Pneumatic Artificial Muscles. IEEE Trans. Robot. Autom. 1996, 12, 90–102. [Google Scholar] [CrossRef]
  38. Tondu, B.; Lopez, P. Modeling and Control of McKibben Artificial Muscle Robot Actuators. IEEE Control Syst. Mag. 2000, 20, 15–38. [Google Scholar] [CrossRef]
  39. Doumit, M.; Fahim, A.; Munro, M. Analytical Modeling and Experimental Validation of the Braided Pneumatic Muscle. IEEE Trans. Robot. 2009, 25, 1282–1291. [Google Scholar] [CrossRef]
  40. Zhang, J.-F.; Yang, C.-J.; Chen, Y.; Zhang, Y.; Dong, Y.-M. Modeling and Control of a Curved Pneumatic Muscle Actuator for Wearable Elbow Exoskeleton. Mechatronics 2008, 18, 448–457. [Google Scholar] [CrossRef]
  41. Kothera, C.; Jangid, M.; Sirohi, J.; Wereley, N. Experimental Characterization and Static Modeling of McKibben Actuators. J. Mech. Des. 2009, 131, 091010-1–091010-10. [Google Scholar] [CrossRef]
  42. Wickramatunge, K.C.; Leephakpreeda, T. Study on Mechanical Behaviors of Pneumatic Artificial Muscle. Int. J. Eng. Sci. 2010, 48, 188–198. [Google Scholar] [CrossRef]
  43. Wickramatunge, K.C.; Leephakpreeda, T. Empirical Modeling of Dynamic Behaviors of Pneumatic Artificial Muscle Actuators. ISA Trans. 2013, 52, 825–834. [Google Scholar] [CrossRef]
  44. Rosen, J.; Fuchs, M.B.; Arcan, M. Performances of Hill-Type and Neural Network Muscle Models-toward a Myosignal-Based Exoskeleton. Comput. Biomed. Res. 1999, 32, 415–439. [Google Scholar] [CrossRef] [Green Version]
  45. Serres, J.L.; Reynolds, D.B.; Phillips, C.A.; Rogers, D.B.; Repperger, D.W. Characterisation of a Pneumatic Muscle Test Station with Two Dynamic Plants in Cascade. Comput. Methods Biomech. Biomed. Eng. 2010, 13, 11–18. [Google Scholar] [CrossRef]
  46. Pitel’, J.; Tóthová, M. Dynamic Modeling of PAM Based Actuator Using Modified Hill’s Muscle Model. In Proceedings of the 14th International Carpathian Control Conference (ICCC), Rytro, Poland, 26–29 May 2013; pp. 307–310. [Google Scholar]
  47. Tóthová, M.; Piteľ, J.; Hošovský, A.; Sárosi, J. Numerical Approximation of Static Characteristics of Mckibben Pneumatic Artificial Muscle. Int. J. Math. Comput. Simul. 2015, 9, 228–233. [Google Scholar]
  48. Piteľ, J.; Tóthová, M. Modelling of Pneumatic Muscle Actuator Using Hill’s Model with Different Approximations of Static Characteristics of Artificial Muscle. MATEC Web Conf. 2016, 76, 02015. [Google Scholar] [CrossRef] [Green Version]
  49. Sárosi, J. New Approximation Algorithm for the Force of Fluidic Muscles. In Proceedings of the 7th IEEE International Symposium on Applied Computational Intelligence and Informatics, Timisoara, Roamnia, 24–26 May 2012; pp. 229–233, ISBN 978-1-4673-1013-0. [Google Scholar]
  50. Das, B.; Dutta, P.P.; Dutta, P.P. Numerical Approximation of Performance of Pneumatic Artificial Muscles with Change in Altitude. In Control and Measurement Applications for Smart Grid; Suhag, S., Mahanta, C., Mishra, S., Eds.; Springer Nature: Singapore, 2022; pp. 219–228. [Google Scholar]
  51. Pietrala, D. The Characteristics of a Pneumatic Muscle. EPJ Web Conf. 2017, 143, 02093. [Google Scholar] [CrossRef]
  52. Donskoj, A.; Kotkas, L.; Salova, T.; Barilovich, V.; Girshov, V.; Mertens, K.; Cherkesov, G.; Donskoy, A.; Fokin, G.; Akhmetova, I. Modelling the Static and Dynamic Characteristics of Pneumatic Muscle. IOP Conf. Ser. Earth Environ. Sci. 2019, 337, 012042. [Google Scholar] [CrossRef]
  53. Wang, H.; Yan, H.; Wang, H.; Yang, Z.; Ni, Z.; Li, Z. Study on Static Characteristics of Pneumatic Muscles. MATEC Web Conf. 2018, 232, 04071. [Google Scholar] [CrossRef]
  54. Takosoglu, J. Static Characteristics of the New Artificial Pneumatic Muscle. EPJ Web Conf. 2022, 269, 01061. [Google Scholar] [CrossRef]
  55. FESTO Fluidic Muscle DMSP/MAS_Catalog. 2017. Available online: https://www.festo.com/cat/en-us_us/data/doc_enus/PDF/US/DMSP-MAS_ENUS.PDF (accessed on 7 November 2022).
  56. FESTO Fluidic Muscle DMSP/MAS_Catalogue. 2008. Available online: https://www.festo.com/rep/en_corp/assets/pdf/info_501_en.pdf (accessed on 7 November 2022).
  57. Nelles, O. Nonlinear System Identification: From Classical Approaches to Neural Networks and Fuzzy Models; Springer: Berlin/Heidelberg, Germany, 2001; ISBN 978-3-540-67369-9. [Google Scholar]
  58. Karray, F.O.; Silva, C.W.D. Soft Computing and Intelligent Systems Design: Theory, Tools and Applications, 1st ed.; Addison-Wesley: Harlow, UK; New York, NY, USA, 2004; ISBN 978-0-321-11617-8. [Google Scholar]
  59. Machová, K. Nové Trendy SU-Štatistický Prístup; Technická Univerzita v Košiciach: Košice, Slovakia, 2016; ISBN 978-80-553-2602-3. [Google Scholar]
  60. Rodríguez-Pérez, R.; Vogt, M.; Bajorath, J. Support Vector Machine Classification and Regression Prioritize Different Structural Features for Binary Compound Activity and Potency Value Prediction. ACS Omega 2017, 2, 6371–6379. [Google Scholar] [CrossRef] [Green Version]
  61. Jang, J.-S.R. ANFIS: Adaptive-Network-Based Fuzzy Inference System. IEEE Trans. Syst. Man Cybern. 1993, 23, 665–685. [Google Scholar] [CrossRef]
  62. Rusgiyono, A. Adaptive Neuro Fuzzy Inference System (ANFIS) Approach for Modeling Paddy Production Data in Central Java. J. Phys. 2018, 1217, 9. [Google Scholar] [CrossRef]
  63. Yang, H.; Yu, Y.; Zhang, J. Angle Tracking of a Pneumatic Muscle Actuator Mechanism under Varying Load Conditions. Control Eng. Pract. 2017, 61, 1–10. [Google Scholar] [CrossRef]
  64. 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, BioRob 2006, Pisa, Italy, 20–22 February 2006; pp. 637–642. [Google Scholar]
  65. Jouppila, V.; Ellman, A. Position Control of PWM-Actuated Pneumatic Muscle Actuator System. In Proceedings of the ASME 2011 International Mechanical Engineering Congress & Exposition, Denver, CO, USA, 11 November 2011; Volume 7, p. 12. [Google Scholar]
  66. Hošovský, A.; Balara, M. Pneumatic Artifical Muscle Force Function Approximation Using Anfis. J. Appl. Sci. Thermodyn. Fluid Mech. 2009, 3, 6. [Google Scholar]
  67. Chen, P.-H.; Fan, R.-E.; Lin, C.-J. A Study on SMO-Type Decomposition Methods for Support Vector Machines. IEEE Trans. Neural Netw. 2006, 17, 893–908. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Wardoyo, R.; Nur Afifa, L. Computing the Time Complexity of ANFIS Algorithm. Int. J. Adv. Res. Comput. Eng. Technol. (IJARCET) 2018, 7, 132–135. [Google Scholar]
  69. Farias, F.S.; Azevedo, R.A.; Rivera, E.C.; Herrera, W.E.; Filho, R.M.; Lima, L.P. Product Quality Monitoring Using Extreme Learning Machines and Bat Algorithms: A Case Study in Second-Generation Ethanol Production. In Computer Aided Chemical Engineering; Klemeš, J.J., Varbanov, P.S., Liew, P.Y., Eds.; 24 European Symposium on Computer Aided Process Engineering; Elsevier: Amsterdam, The Netherlands, 2014; Volume 33, pp. 955–960. [Google Scholar]
  70. Pham, H. A New Criterion for Model Selection. Mathematics 2019, 7, 1215. [Google Scholar] [CrossRef] [Green Version]
  71. Nainggolan, R.; Perangin-angin, R.; Simarmata, E.; Tarigan, A.F. Improved the Performance of the K-Means Cluster Using the Sum of Squared Error (SSE) Optimized by Using the Elbow Method. J. Phys. Conf. Ser. 2019, 1361, 012015. [Google Scholar] [CrossRef]
  72. The MathWorks, Inc. System Identification Toolbox User’s Guide. 2022. Available online: https://www.mathworks.com/help/pdf_doc/ident/ident_ug.pdf (accessed on 7 November 2022).
  73. Hošovsky, A.; Trojanova, M.; Pitel, J.; Svetlik, J. Identification of DMSP-5 Fluidic Muscle Dynamics Using Hammerstein Model. In Proceedings of the 2018 IEEE 12th International Symposium on Applied Computational Intelligence and Informatics (SACI), Timisoara, Romania, 17–19 May 2018; pp. 000319–000324. [Google Scholar]
  74. Neath, A.A.; Cavanaugh, J.E. The Bayesian Information Criterion: Background, Derivation, and Applications. WIREs Comput. Stat. 2012, 4, 199–203. [Google Scholar] [CrossRef]
  75. Chen, J.; Chen, Z. Extended Bayesian Information Critera for Model Selection with Large Model Spaces. Biometrika 2008, 95, 759–771. [Google Scholar] [CrossRef] [Green Version]
  76. Cavanaugh, J.E.; Neath, A.A. The Akaike Information Criterion: Background, Derivation, Properties, Application, Interpretation, and Refinements. WIREs Comput. Stat. 2019, 11, e1460. [Google Scholar] [CrossRef]
  77. Dridi, N.; Hadzagic, M. Akaike and Bayesian Information Criteria for Hidden Markov Models. IEEE Signal Process. Lett. 2019, 26, 302–306. [Google Scholar] [CrossRef]
Figure 1. Fluidic muscles with press-fitted connections: (a) DMSP-40, (b) DMSP-20, and (c) DMSP-5 [55].
Figure 1. Fluidic muscles with press-fitted connections: (a) DMSP-40, (b) DMSP-20, and (c) DMSP-5 [55].
Mathematics 11 00149 g001
Figure 2. Experimental system: (a) 3-DOF planar robotic arms actuated with DMSP-20 and DMSP-40 fluidic muscles; (b) 1-DOF arm with multiparallel connection of 6 DMSP-5 fluidic muscles.
Figure 2. Experimental system: (a) 3-DOF planar robotic arms actuated with DMSP-20 and DMSP-40 fluidic muscles; (b) 1-DOF arm with multiparallel connection of 6 DMSP-5 fluidic muscles.
Mathematics 11 00149 g002
Figure 3. Schematic of the test stand used for measuring the characteristics of the DMSP-5 muscle.
Figure 3. Schematic of the test stand used for measuring the characteristics of the DMSP-5 muscle.
Mathematics 11 00149 g003
Figure 4. Operating range for the DMSP-40 muscle [56].
Figure 4. Operating range for the DMSP-40 muscle [56].
Mathematics 11 00149 g004
Figure 5. Force simulated (Output) by the ANFIS, custom, MLP, and SVM models in correlation with the force measured by FESTO (Target) for DMSP-40: test phase.
Figure 5. Force simulated (Output) by the ANFIS, custom, MLP, and SVM models in correlation with the force measured by FESTO (Target) for DMSP-40: test phase.
Mathematics 11 00149 g005
Figure 6. Relations of force, contraction, and pressure approximated by the ANFIS, custom, MLP, and SVM models for the DMSP-40 muscle.
Figure 6. Relations of force, contraction, and pressure approximated by the ANFIS, custom, MLP, and SVM models for the DMSP-40 muscle.
Mathematics 11 00149 g006
Figure 7. Operating range for the DMSP-20 muscle [56].
Figure 7. Operating range for the DMSP-20 muscle [56].
Mathematics 11 00149 g007
Figure 8. Force simulated (Output) by the ANFIS, custom, MLP, and SVM models in correlation with the force measured by FESTO (Target) for DMSP-20: test phase.
Figure 8. Force simulated (Output) by the ANFIS, custom, MLP, and SVM models in correlation with the force measured by FESTO (Target) for DMSP-20: test phase.
Mathematics 11 00149 g008
Figure 9. Relations of force, contraction, and pressure approximated by the ANFIS, custom, MLP, and SVM models for the DMSP-20 muscle.
Figure 9. Relations of force, contraction, and pressure approximated by the ANFIS, custom, MLP, and SVM models for the DMSP-20 muscle.
Mathematics 11 00149 g009
Figure 10. Operating range for the DMSP-5 muscle [55].
Figure 10. Operating range for the DMSP-5 muscle [55].
Mathematics 11 00149 g010
Figure 11. Force characteristics of FESTO’s DMSP-5 muscle measured using the experimental setup shown in Figure 3.
Figure 11. Force characteristics of FESTO’s DMSP-5 muscle measured using the experimental setup shown in Figure 3.
Mathematics 11 00149 g011
Figure 12. Force simulated (Output) by the ANFIS, custom, MLP, and SVM models in correlation with the force measured (Target) for DMSP-5: test phase.
Figure 12. Force simulated (Output) by the ANFIS, custom, MLP, and SVM models in correlation with the force measured (Target) for DMSP-5: test phase.
Mathematics 11 00149 g012
Figure 13. Relations of force, contraction, and pressure approximated by the ANFIS, custom, MLP, and SVM models for the DMSP-5 muscle.
Figure 13. Relations of force, contraction, and pressure approximated by the ANFIS, custom, MLP, and SVM models for the DMSP-5 muscle.
Mathematics 11 00149 g013
Figure 14. Comparison of the model outputs with the force measured by FESTO for the DMSP-40 muscle.
Figure 14. Comparison of the model outputs with the force measured by FESTO for the DMSP-40 muscle.
Mathematics 11 00149 g014
Figure 15. Comparison of the model outputs with the force measured by FESTO for the DMSP-20 muscle.
Figure 15. Comparison of the model outputs with the force measured by FESTO for the DMSP-20 muscle.
Mathematics 11 00149 g015
Figure 16. Comparison of the model outputs with the measured force for the DMSP-5 muscle.
Figure 16. Comparison of the model outputs with the measured force for the DMSP-5 muscle.
Mathematics 11 00149 g016
Figure 17. Summary of the results of the models of all architectures for each of the modeled muscle types using the selected performance metrics.
Figure 17. Summary of the results of the models of all architectures for each of the modeled muscle types using the selected performance metrics.
Mathematics 11 00149 g017
Table 1. Chronological overview of some braided pneumatic artificial muscles [16].
Table 1. Chronological overview of some braided pneumatic artificial muscles [16].
Name of MuscleAuthorYear
of Invention
References
McKibben muscleGaylord1953[17,18,19]
Air muscleGreenhill, Shadow Robot Company1988[20,21,22]
Under-pressure artificial muscleMarcincin, Palko1993[23,24,25]
Actuating means, fluidic muscleBergemann, Lorenz, Thallemer, FESTO AG & Co.1999[26,27,28,29]
Fluidic actuatorDavis, Carlson, Electro Cam Corporation2003[30,31]
Table 2. Overview of the related works.
Table 2. Overview of the related works.
AuthorYearCiteType
of Muscle
ModelExperimental StandData Obtained
from Experimental Stand
IsotonicIsobaricIsometric
Doumit et al.2009[39]McKibbenGeometricYesNoNoYes
Kothera et al.2009[41]McKibbenGeometricYesYesNoYes
Sárosi2012[49]FluidicGeometricYesNoYesNo
Wickramatunge et al.2013[43]FluidicEmpiricalYesYesYesNo
Tóthová et al.2015[47]FluidicGeometric/empiricalNo---
Takosoglu et al.2016[16]Fluidic
McKibben
GeometricYesYesYesYes
Pietrala 2017[51]McKibben-YesYesYesYes
Wang et al.2018[53]FluidicGeometricYesNoNoNo
Donskoj et al.2019[52]FluidicGeometricNo---
Takosoglu 2022[54]McKibben-YesYesYesYes
Table 3. Overview of the technical data of the selected DMSP muscles [55,56].
Table 3. Overview of the technical data of the selected DMSP muscles [55,56].
Type of MuscleDMSP-5DMSP-20DMSP-40
Internal diameter (mm)52040
Nominal length (mm)30 … 100060 … 9000120 … 9000
Max. permissible pretensioning1% of nominal length4% of nominal length5% of nominal length
Max. permissible contraction20% of nominal length25% of nominal length25% of nominal length
Theoretical force (N)14015006000
Max. additional load (kg)580250
Operating pressure (bar)0 … 60 … 60 … 6
Ambient temperature (°C)−5 … +60−5 … +60−5 … +60
Table 4. Division ratios for training, validating, and testing the created models.
Table 4. Division ratios for training, validating, and testing the created models.
Type of MuscleTotal Number
of Samples
Data Division Ratio %Sample Division
DMSP-523170/15/15161/35/35
DMSP-2045570/15/15319/68/68
DMSP-4049770/15/15347/75/75
Table 5. Overview of the parameters of the created MLP models and performance metrics of the DMSP-40 muscle validation results.
Table 5. Overview of the parameters of the created MLP models and performance metrics of the DMSP-40 muscle validation results.
ModelMLP 1MLP 2MLP 5MLP 10MLP 20
Number of neurons1251020
Number of parameters59214181
NRMSE 60.3162.8794.7797.9590.36
SSE259,540,000227,080,0004,500,700689,79015,308,000
AIC6553.464954570.23678.15298.6
SBC6574.56532.94658.63850.65639.5
Table 6. Overview of the parameters of the created SVM models and performance metrics of the DMSP-40 muscle validation results.
Table 6. Overview of the parameters of the created SVM models and performance metrics of the DMSP-40 muscle validation results.
ModelSVM 1SVM 2SVM 3SVM 4SVM 5SVM 6
Kernel scale0.771620.771620.771620.351.45.7
Bias1.45 × 1031.16 × 1031.11 × 1031.67 × 1032.24 × 1033.84 × 103
Epsilon1.76 × 1021.76 × 1021.76 × 1021.76 × 1021.76 × 1021.76 × 102
FunctionLinearQuadraticCubicFine GaussianMedium GaussianCoarse Gaussian
Alpha3161379413273151
Beta200000
NRMSE46.2854.9452.6175.5261.6957.32
SSE4.75 × 1083.34 × 1083.70 × 1089.87 × 1072.42 × 1083.00 × 108
AIC7.48 × 1036.94 × 1036.91 × 1036.33 × 1036.65 × 1036.92 × 103
SBC8.82 × 1037.52 × 1037.30 × 1036.88 × 1036.96 × 1037.55 × 103
Table 7. Overview of the parameters of the created ANFIS models and performance metrics of the DMSP-40 muscle validation results.
Table 7. Overview of the parameters of the created ANFIS models and performance metrics of the DMSP-40 muscle validation results.
ModelANFIS 1ANFIS 2ANFIS 3ANFIS 4ANFIS 5ANFIS 6ANFIS 7ANFIS 8
Input MF4 4 6 6 8 8 10 10 4 4 6 6 8 8 10 10
MF typeConstantConstantConstantConstantLinearLinearLinearLinear
Parameters32609614064132224340
NRMSE63.2768.9973.9577.5275.6578.8381.0484.09
SSE2.22 × 1081.58 × 1081.12 × 1088.33 × 1079.76 × 1077.38 × 1075.92 × 1074.17 × 107
AIC6.53 × 1036.42 × 1036.32 × 1036.26 × 1036.19 × 1036.18 × 1036.26 × 1036.31 × 103
SBC6.67 × 1036.67 × 1036.72 × 1036.85 × 1036.45 × 1036.74 × 1037.20 × 1037.75 × 103
Table 8. Overview of the parameters of the created custom model and performance metrics of the DMSP-40 muscle validation results.
Table 8. Overview of the parameters of the created custom model and performance metrics of the DMSP-40 muscle validation results.
Model ParametersMetrics
a10.2629NRMSE53.18
a2−39.51 SSE3.6116 × 108
a3924.9 AIC6.7156 × 103
a494.56 SBC6.7325 × 103
Table 9. Overview of the parameters of the created MLP models and performance metrics of the DMSP-20 muscle validation results.
Table 9. Overview of the parameters of the created MLP models and performance metrics of the DMSP-20 muscle validation results.
ModelMLP 1MLP 2MLP 5MLP 10MLP 20
Number of neurons1251020
Number of parameters59214181
NRMSE39.8666.2081.6289.9088.74
SSE38,800,00012,255,0003,623,1001,093,7001,361,200
AIC5175.94659.54129.13624.13803.6
SBC5196.54696.64215.637934137.4
Table 10. Overview of the parameters of the created SVM models and performance metrics of the DMSP-20 muscle validation results.
Table 10. Overview of the parameters of the created SVM models and performance metrics of the DMSP-20 muscle validation results.
ModelSVM 1SVM 2SVM 3SVM 4SVM 5SVM 6
Kernel scale0.8351170.8351170.8351170.351.45.7
Bias331.7927240.2869221.3975367.2308460.73361055.019
Epsilon44.097244.097244.097244.097244.097244.0972
FunctionLinearQuadraticCubicFine GaussianMedium GaussianCoarse Gaussian
Alpha319150138145118194
Beta200000
NRMSE29.6222.8219.4663.4543.4730.46
SSE5.31 × 1076.39 × 1076.96 × 1071.43 × 1073.43 × 1075.19 × 107
AIC5.95 × 1035.69 × 1035.71 × 1035.00 × 1035.35 × 1035.69 × 103
SBC7.27 × 1036.31 × 1036.28 × 1035.60 × 1035.83 × 1036.49 × 103
Table 11. Overview of the parameters of the created ANFIS models and performance metrics of the DMSP-20 muscle validation results.
Table 11. Overview of the parameters of the created ANFIS models and performance metrics of the DMSP-20 muscle validation results.
ModelANFIS 1ANFIS 2ANFIS 3ANFIS 4ANFIS 5ANFIS 6ANFIS 7ANFIS 8
Input MF4 4 6 6 8 8 10 10 4 4 6 6 8 810 10
MF typeConstantConstantConstantConstantLinearLinearLinearLinear
Parameters32609614064132224340
NRMSE44.1955.0857.9059.7060.8064.7665.7465.41
SSE3.34 × 1072.16 × 1071.90 × 1071.74 × 1071.71 × 1071.33 × 1071.26 × 1071.28 × 107
AIC5.16 × 1035.22 × 1035.29 × 1035.38 × 1035.23 × 1035.36 × 1035.55 × 1035.78 × 103
SBC5.29 × 1035.47 × 1035.69 × 1035.95 × 1035.49 × 1035.91 × 1036.47 × 1037.18 × 103
Table 12. Overview of the parameters of the created custom model and performance metrics of the DMSP-20 muscle validation results.
Table 12. Overview of the parameters of the created custom model and performance metrics of the DMSP-20 muscle validation results.
Model ParametersMetrics
a1−0.1398 NRMSE20.52
a2−4.77 SSE6.7771 × 107
a3182.7 AIC5.4277 × 103
a477.43 SBC5.4441 × 103
Table 13. Overview of the parameters of the created MLP models and performance metrics of the DMSP-5 muscle validation results.
Table 13. Overview of the parameters of the created MLP models and performance metrics of the DMSP-5 muscle validation results.
ModelMLP 1MLP 2MLP 5MLP 10MLP 20
Number of neurons1251020
Number of parameters59214181
NRMSE80.9186.5695.4597.1198.74
SSE8.32 × 1034.12 × 103473.054190.327236.5044
AIC837.8817683.7842207.579137.2615−264.1914
SBC855.0938714.766279.8698178.400614.6444
Table 14. Overview of the parameters of the created SVM models and performance metrics of the DMSP-5 muscle validation results.
Table 14. Overview of the parameters of the created SVM models and performance metrics of the DMSP-5 muscle validation results.
ModelSVM 1SVM 2SVM 3SVM 4SVM 5SVM 6
Kernel scale1.2971891.2971891.2971890.351.45.7
Bias14.644995.5236595.47279219.6549131.0298658.20743
Epsilon2.1126762.1126762.1126762.1126762.1126762.112676
FunctionLinearQuadraticCubicFine GaussianMedium GaussianCoarse Gaussian
Alpha18013811912570156
Beta200000
NRMSE36.9676.7681.9371.4880.7655.03
SSE9.07 × 1041.23 × 1047.46 × 1031.86 × 1048.45 × 1034.62 × 104
AIC1.74 × 1031.66 × 1031.04 × 1031.26 × 103971.54561.54 × 103
SBC7.27 × 1036.31 × 1036.28 × 1035.60 × 1035.83 × 1036.49 × 103
Table 15. Overview of the parameters of the created ANFIS models and performance metrics of the DMSP-5 muscle validation results.
Table 15. Overview of the parameters of the created ANFIS models and performance metrics of the DMSP-5 muscle validation results.
ModelANFIS 1ANFIS 2ANFIS 3ANFIS 4ANFIS 5ANFIS 6ANFIS 7ANFIS 8
Input MF4 4 6 6 8 8 10 10 4 4 6 68 8 10 10
MF typeConstantConstantConstantConstantLinearLinearLinearLinear
Parameters32609614064132224340
NRMSE81.4087.2890.5293.2693.7196.3186.1489.97
SSE7.90 × 1033.70 × 1032.05 × 1031.04 × 1039.03 × 1023.10 × 1024.38 × 1032.30 × 103
AIC8.80 × 1027.60 × 1026.96 × 1026.27 × 1024.43 × 1023.32 × 1021.13 × 1031.21 × 103
SBC9.90 × 1029.67 × 1021.03 × 1031.11 × 1036.63 × 1027.87 × 1021.90 × 1032.38 × 103
Table 16. Overview of the parameters of the created custom model and performance metrics of the DMSP-5 muscle validation results.
Table 16. Overview of the parameters of the created custom model and performance metrics of the DMSP-5 muscle validation results.
Model ParametersMetrics
a10.07496NRMSE72.86
a2−2.533SSE1.6812 × 104
a325.39AIC998.3944
a4−16.72SBC1.0122 × 103
Table 17. Summary of the best results for each type of DMSP fluidic muscle.
Table 17. Summary of the best results for each type of DMSP fluidic muscle.
MetricsDMSP-40DMSP-20DMSP-5
MLP 10MLP 10MLP 20
AIC3678.13624.1−264.1914
SBC3850.6379314.6444
NRMSE97.95%89.90%98.74%
Table 18. Performance comparison of several state-of-art models with the best ML model.
Table 18. Performance comparison of several state-of-art models with the best ML model.
ModelAnalyticalMaximum ForceExponential FunctionPolynomial FunctionMLP 10
[47,64][47,65][47,49][47,48]
Number
of parameters
5862141
Number
of samples
417417417417455
RMSE81.22473209.2164.549.02791
SSE2,784,00093,740,00018,420,00010,990,0001,093,700
AIC3682.25154.74472.24286.83624.1
SBC3702.45186.94496.44371.53793
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Trojanová, M.; Hošovský, A.; Čakurda, T. Evaluation of Machine Learning-Based Parsimonious Models for Static Modeling of Fluidic Muscles in Compliant Mechanisms. Mathematics 2023, 11, 149. https://doi.org/10.3390/math11010149

AMA Style

Trojanová M, Hošovský A, Čakurda T. Evaluation of Machine Learning-Based Parsimonious Models for Static Modeling of Fluidic Muscles in Compliant Mechanisms. Mathematics. 2023; 11(1):149. https://doi.org/10.3390/math11010149

Chicago/Turabian Style

Trojanová, Monika, Alexander Hošovský, and Tomáš Čakurda. 2023. "Evaluation of Machine Learning-Based Parsimonious Models for Static Modeling of Fluidic Muscles in Compliant Mechanisms" Mathematics 11, no. 1: 149. https://doi.org/10.3390/math11010149

APA Style

Trojanová, M., Hošovský, A., & Čakurda, T. (2023). Evaluation of Machine Learning-Based Parsimonious Models for Static Modeling of Fluidic Muscles in Compliant Mechanisms. Mathematics, 11(1), 149. https://doi.org/10.3390/math11010149

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