Truth and Myths about 2D Tensegrity Trusses

Featured Application: We present an evaluation of 2D structural models of specific truss structures named tensegrities. Proper understanding of this unusual structure (according to the proposed definition and methods of analysis) is necessary for specific applications in building and bridge construction, deployable structures, and smart structures with monitoring, self-control, self-repair, and active control. New and future potential applications of the concept are tensegrity inspired metamaterials with unique mechanical properties. Abstract: The concept of tensegrity is understood in many ways. This term is often improperly used for structures that have some, but not necessarily the key, tensegrity properties. The concept of tensegrity systems is misused in reference to both mathematical models and completed engineering structures. The aim of the study is to indicate which of the plane (2D) trusses presented in the literature are erroneously classiﬁed as tensegrities. Singular value decomposition of the compatibility matrix and spectral analysis of the stiffness matrix with the effect of self-equilibrated forces is used for the analysis. A new precise deﬁnition of tensegrity trusses is proposed and implemented.


Introduction
Tensegrity structures consist of compression (struts) and tensile (cables) components which stabilize each other through pre-stress. The idea of tensegrity was first described about 60 years ago, but no precise, nor general, definition has currently been agreed upon. This thereby paves the way for confusion and misuse of the term in technical applications. The term is often used for trusses that have little to do with the tensegrity principle as defined by the creators of this idea.
The first patented definitions of tensegrity are very general. Fuller [1] described tensegrity systems as "islands of compression in a sea of tension". Snelson patented the system [2] which he called "Continuous Tension, Discontinuous Compression Structures." Emmerich [3] added the condition of self-stress state: Tensegrity structures consist of rods mounted in such a way that the struts remain physically isolated in a continuous set of cables. All these elements must be very close to one another and interconnected by tensional pre-stress. No external supports or anchors are required. It makes a structure as strong as the self-supporting structure, hence the name self-stress state. The most widely accepted definition of tensegrity is formulated by Pugh [4], in which he combines the two previous statements: "A tensegrity system is established when a set of discontinuous compression components interacts with a set of continuous tensile components to define a stable volume in space." This definition was narrowed by Motro [5,6]: "A tensegrity system is a system in a stable self-equilibrated state comprising a discontinuous set of compressed components inside a continuum of tensioned components." Hanaor [7] describes tensegrity structures as "internally The concept of tensegrity systems is misused in reference to both mathematical models and completed engineering structures where the terms "tensegrity bridges" or "tensegrity domes" are applied to the cable-strut models. The proposed evaluation method will allow proper understanding of these unique structures and ensure the correct use of the term in all applications.

Definition of Tensegrity
With regard to structural mechanics, the presence of self-stress states (understood as the system of self-equilibrated normal forces which satisfies homogeneous equations of equilibrium) and infinitesimal mechanisms which are stabilized by these states are the most important features of tensegrity structures. In this paper, considerations are limited to "pin-jointed plane structures" and the definition of tensegrity formulated by Kasprzak [20] is used. Based on the definitions known from the literature, tensegrity structures demonstrate the following characteristic features: • T-the structure is a truss, • S-there is a self-stress state, • M-there is an infinitesimal mechanism stiffened by the self-stress state, • D-the extremities of compressed components (struts) do not touch each other and struts constitute a discontinuous set, • I-the set of compressed components is included inside the set of tensile components, • C-tensile elements have no rigidity in compression-these are cables.
The characteristics above are used to classify structures as "pure tensegrities" or "structures with tensegrity features".
"Pure tensegrity" is represented by a structure that closely meets all the requirements arising from the definition of tensegrity and defining its specific properties. A tensegrity structure is a truss in which there is an infinitesimal mechanism stiffened by a self-stress state. Components in compression create a discontinuous system inside a continuous system of components in tension, which have no compressive rigidity. The structures that have all of the six features (T + S + M + D + I + C) are classified in the "pure tensegrity" group.
"Structures with tensegrity features" include a wider range of systems in which some properties of tensegrity can be ignored, e.g. a truss in a self-stress state, with tensile elements having no compressive rigidity, demonstrating at least one of the following features: The presence of an infinitesimal mechanism stabilized by self-stress state, the discontinuity of the system of compressed components, and/or the use of compressed elements inside a set of tension components. Structures assigned to this group fulfil three obligatory criteria (T + S + C) and have at least one of the features: M, D, or I.
From an engineering perspective, it is very important that the structures have all six features. Only pure tensegrities possess all the benefits of tensegrity structures, including the ability to control their stiffness and natural frequency by modifying the self-stress level.
Singular value decomposition of the compatibility matrix is used to determine the states of self-stress (S) and mechanisms (M). The nonlinear analysis with the use of geometric stiffness matrix, which takes into account the impact of the self-stress state, allows establishing whether the mechanism is infinitesimal or finite. The other characteristic features of tensegrity structures like D, I, and C are determined indirectly by the identified self-stress state.

Mathematical Description
The plane truss finite element of Young modulus E e , cross-sectional area A e and lengths L e (Figure 1a) is used to model the components of tensegrity structures. In global coordinate system (x, y) a truss element is described respectively by the compatibility matrix B e ∈ R 1×4 and by the geometric stiffness matrix K Ge ∈ R 4×4 , taking into account the impact of the self-stress state S e : where: c = The qualitative analysis focuses on n-membered (e = 1, 2, . . . , n) two-dimensional trusses, described by elasticity matrix E(∈ R n×n ): with m degrees of freedom q ∈ R m×1 : The compatibility matrix B(∈ R n×m ) and the geometric stiffness matrix K G (∈ R m×m ) for trusses are determined using the finite element formalism [49][50][51]: where C e ∈ R 4×m is a Boolean matrix. It is assumed that the number of global degrees of freedom of each element q i (i = 1, 2, 3, 4) corresponds to the number of global nodes of elements n 1 , n 2 as shown in Figure 1b. Consequently, the non-zero elements of C e for each element can be expressed as C iq i = 1.
The qualitative analysis focuses on -membered ( = 1,2, … , ) two-dimensional trusses, described by elasticity matrix (∈ ℝ × ): with degrees of freedom (∈ ℝ × ): The compatibility matrix (∈ ℝ × ) and the geometric stiffness matrix (∈ ℝ × ) for trusses are determined using the finite element formalism [49][50][51]: where (∈ ℝ × ) is a Boolean matrix. It is assumed that the number of global degrees of freedom of each element (i = 1, 2, 3, 4) corresponds to the number of global nodes of elements , as shown in Figure 1b. Consequently, the non-zero elements of for each element can be expressed as = 1. Equilibrium equations for trusses can be written in a stress form: where (∈ ℝ × ) is the normal force vector and (∈ ℝ × ) is the load vector, or via nodal displacements: where (∈ ℝ × ) is the linear stiffness matrix: To take into account the non-linearity of the equilibrium equations and the influence of the prestress (self-stress) in the structure, the geometric stiffness matrix (4)2 is introduced. Then the system of displacement equations (6) takes the form: A mathematical description of tensegrity-like pin-jointed structures can be found in [51][52][53][54][55][56].
The qualitative analysis of the set of equations for the structure can be done through singular value decomposition [57][58][59][60][61][62][63] of the compatibility matrix [64][65][66]. The singular value decomposition of matrix is a factorization in the form: Equilibrium equations for trusses can be written in a stress form: where S ∈ R n×1 is the normal force vector and P ∈ R m×1 is the load vector, or via nodal displacements: where K L (∈ R m×m ) is the linear stiffness matrix: To take into account the non-linearity of the equilibrium equations and the influence of the prestress (self-stress) in the structure, the geometric stiffness matrix (4) 2 is introduced. Then the system of displacement equations (6) takes the form: A mathematical description of tensegrity-like pin-jointed structures can be found in [51][52][53][54][55][56]. The qualitative analysis of the set of equations for the structure can be done through singular value decomposition [57][58][59][60][61][62][63] of the compatibility matrix B [64][65][66]. The singular value decomposition of matrix B is a factorization in the form: where Y(∈ R n×n ) and X(∈ R m×m ) are the orthogonal matrices and N(∈ R n×m ) is rectangular diagonal matrix. The orthogonal matrices Y and X as well as matrix N are related to eigenvectors and eigenvalues of the following problems: The solutions to the (10) and (11) eigen-problems, can be expressed as: The eigenvalues matrices M(∈ R n×n ) and L(∈ R m×m ) are dependent on the matrix N of the decomposition (9): One can notice that the product BB T can be considered as a matrix of symmetrized equations of equilibrium (5) with non-negative eigenvalues. Zero eigenvalues (if any) are related to the non-zero solution of homogeneous equations called a self-stress, or more precisely, self-equilibrated normal forces that satisfies homogeneous equations of equilibrium. The self-stress can be considered an eigenvector related to zero eigenvalue. In a similar way, the product B T B can be considered as a particular form of a linear stiffness matrix (7) with unit elasticity matrix E ≡ I. The eigenvalues are non-negative. Zero eigenvalues (if any) are related to the finite or infinitesimal mechanisms, which can be considered an eigenvector related to zero eigenvalue. Finite mechanisms (rigid-body motions) refer to the motions that do not change the distance between any pair of nodes. Infinitesimal mechanisms describe the local geometrical variability in the range of small displacements. Generally, the information from the null-space analysis alone is insufficient to determine the difference.
To establish whether the mechanism is infinitesimal or finite, it is necessary to apply nonlinear analysis with the use of geometric stiffness matrix (8), which is possible if the self-stress exists. If the eigenvalues of (K L + K G ) are positive numbers, the mechanism is infinitesimal. It means the structure is stable. Zero eigenvalues are related to finite mechanisms. A negative eigenvalue represents instability of the structure.
Full solutions to the eigen-problem of the stiffness matrix, taking into account the geometric stiffness matrix: can be expressed in the condensed form: where: are the matrices of eigenvalues and of normalized eigenvectors, respectively.

Results
The aim of the study is to indicate which of the plane (2D) trusses presented in the literature are wrongly classified as tensegrities. To do this, the qualitative analysis of selected trusses was carried out. Singular value decomposition of compatibility matrix B and spectral analysis of stiffness matrix with the effect of self-equilibrated forces (K L + K G ) were used for the analysis. Since the characteristic features of tensegrity structures do not depend on geometrical and material characteristics, it was assumed that the constant longitudinal stiffness E e A e = 1 and constant lengths a = 1. Eighteen plane tensegrity structures were analyzed. The results were presented in two ways (mathematical and graphical)-for selected structures the following matrices were determined: For the purpose of this paper, the plane (2D) trusses defined in the literature as "tensegrities" were classified into three groups: "Structures with tensegrity features", "pure tensegrities" and "non tensegrities". Figure 2a shows the first structure under analysis. It is the most popular structure in the literature [8,12,[15][16][17][18]37,41,42,64,[67][68][69][70][71][72][73][74][75][76][77][78], usually referred to as "the simplest 2D tensegrity structure" [77]. It is a truss (T) with six elements (n = 6) and five degrees of freedom (d.o.f.) (m = 5)-q = q 3 q 4 q 5 q 6 q 7 T . For this truss the compatibility matrix B ∈ R 6×5 is: There is no zero eigenvalue in L, which means that no mechanism is identified. The zero eigenvalue in M is responsible for the existence of the self-stress state (S): In this case, all eigenvalues of the matrix (K L + K G ) are positive: which means the structure is stable. For this truss the system of compression members (thick line) is discontinuous (D) and included inside the set of components in tension (I), which are cables (C). However, because of the absence of the mechanism, the truss shares only some features with a tensegrity.
If unsupported "X" truss is analyzed (Figure 2b), there are eight d.o.f. (m = 8) and the compatibility matrix B ∈ R 6×8 is as follows: Matrices L ∈ R 8×8 and M ∈ R 6×6 are respectively: The zero eigenvalue in is responsible for the existence of the self-stress state (S) defined by the same eigenvector as for the supported truss: In this case, three eigenvalues of the matrix ( + ) are zero: The zero eigenvalue in M is responsible for the existence of the self-stress state (S) defined by the same eigenvector as for the supported truss: In this case, three eigenvalues of the matrix (K L + K G ) are zero: The truss shown in Figure 4d has three zero eigenvalues in matrix M responsible for the existence of three self-stress states (S). These states need to be considered together. Superposition and normalization of self-stress states were performed, as shown in Figure 4d. For this truss the system of compression members is inside the set of tensioned elements (I), which are cables (C) but it does not satisfy the condition of discontinuity. According to the Skelton's definition [13], it is a class 2 tensegrity system-a system with two compressive elements connected at the nodes.  For these trusses there are no zero eigenvalues in matrix , which means that the system lacks the mechanisms. From the perspective of structural mechanics, the mechanism is one of the most important tensegrity features-the uniqueness of tensegrity structures is based on the fact that states of self-stress stabilize existing mechanisms. Because of the absence of the mechanism, these trusses For these trusses there are no zero eigenvalues in matrix L, which means that the system lacks the mechanisms. From the perspective of structural mechanics, the mechanism is one of the most important tensegrity features-the uniqueness of tensegrity structures is based on the fact that states of self-stress stabilize existing mechanisms. Because of the absence of the mechanism, these trusses are only structures with tensegrity features. The term "tensegrity" used for the above examples in the references [6,8,12,[15][16][17][18]20,37,41,42,64,[67][68][69][70][71][72][73][74][75][76][77][78][79][80][81][82][83][84][85][86] is not sufficiently accurate.
The structure in Figure 5a represents Similar but more practical example (to be used as a computational model of stadium or hall roofs) is presented in Figure 6. ] . This example is developed by extending the truss showed in Figure 5a. Here, the compatibility matrix is rectangular (∈ ℝ × ). The matrices (∈ ℝ × ) and (∈ ℝ × ) are respectively: which means that mechanisms are balanced by the self-stress state-the mechanisms are infinitesimal. For this truss the system of compression members (thick line) is discontinuous (D) and included inside the set of elements in tension (I), which are cables (C). This is a pure tensegrity with two mechanisms stabilized by one self-stress. A selection of pure tensegrities found in the literature is shown in Figure 8: thus the mechanism is balanced by the self-stress state-the mechanism is infinitesimal (describes the local geometrical variability). For this truss, the system of compression members (thick red line) is discontinuous (D) and included inside the set of elements in tension (I), which are cables (C). This is a pure tensegrity with one infinitesimal mechanism stabilized by one self-stress state. Similar but more practical example (to be used as a computational model of stadium or hall roofs) is presented in Figure 6. The structure is a truss (T) with 15 elements (n = 15) and 16 d.o.f.
(m = 16) [20,77,85,87]-q = q 3 q 4 q 5 q 6 q 7 q 8 q 9 q 10 q 13 q 14 q 15 q 16 q 17 q 18 T . This example is developed by extending the truss showed in Figure 5a. Here, the compatibility matrix is rectangular B ∈ R 15×16 . The matrices L ∈ R 16×16 and M ∈ R 15×15 are respectively: (thick red line) is discontinuous (D) and included inside the set of components in tension (I), which are cables (C).
For the first two cases (Figure 8a,b) there is one zero eigenvalue in matrix responsible for the existence of one mechanism (M). For the last two cases (Figure 8c,d) two mechanisms are identified. All mechanisms (Figure 9) are infinitesimal because eigenvalues of ( + ) are positive. The analyzed trusses are pure tensegrities with infinitesimal mechanisms stabilized by one self-stress. According to the features presented above the tensegrities presented in Figure 8 are promising as unit cells for the applications of 2D metamaterials with unusual mechanical properties [30]. (thick red line) is discontinuous (D) and included inside the set of components in tension (I), which are cables (C). For the first two cases (Figure 8a,b) there is one zero eigenvalue in matrix responsible for the existence of one mechanism (M). For the last two cases (Figure 8c,d) two mechanisms are identified. All mechanisms (Figure 9) are infinitesimal because eigenvalues of ( + ) are positive. The analyzed trusses are pure tensegrities with infinitesimal mechanisms stabilized by one self-stress. According to the features presented above the tensegrities presented in Figure 8 are promising as unit cells for the applications of 2D metamaterials with unusual mechanical properties [30]. For all the trusses there is one zero eigenvalue in matrix M which is responsible for the existence of the one self-stress state (S) shown in Figure 8. For these trusses the system of compression members (thick red line) is discontinuous (D) and included inside the set of components in tension (I), which are cables (C).
For the first two cases (Figure 8a,b) there is one zero eigenvalue in matrix L responsible for the existence of one mechanism (M). For the last two cases (Figure 8c,d) two mechanisms are identified. All mechanisms (Figure 9) are infinitesimal because eigenvalues of (K L + K G ) are positive. The analyzed trusses are pure tensegrities with infinitesimal mechanisms stabilized by one self-stress.
According to the features presented above the tensegrities presented in Figure 8 are promising as unit cells for the applications of 2D metamaterials with unusual mechanical properties [30].
(c) (d) Figure 8. Values of self-stress state for four pure tensegrity trusses (a-d). An interesting, as well as instructive, example is presented in Figure 10a. It is the simplest truss (T) with two elements ( = 2) and two d.o.f. ( = 2) [6,8,20,64,77] = [ ] . Here, the compatibility matrix is square (∈ ℝ × ), thus the matrices and (∈ ℝ × ) are the same: The zero eigenvalue in and is responsible for the existence of the mechanism (M) ( Figure  10b An interesting, as well as instructive, example is presented in Figure 10a. It is the simplest truss which indicates that the mechanism is balanced by the self-stress state. This structure does not satisfy the conditions of discontinuous compression members and the inclusion of compression elements within the set of elements in tension because there are no compression elements-the self-stress is balanced by reactions of the support. The truss is not a tensegrity because of lack of C, D, and I features. On the other hand, this is a good and instructive example to show the intuitive infinitesimal mechanism stabilized by also intuitive self-stress.  (Figure 11d).
For the first three examples (Figure 11a to c), there is no zero eigenvalue in the matrix which means that the system lacks the mechanism. There are three zero eigenvalues in the matrix , which are responsible for the existence of three different self-stress states (S). Two of them lead to the positive eigenvalues of the matrix ( + ), and the third self-stress gives one negative eigenvalue. Several combinations of the three identified self-stress states are possible. It is impossible to clearly define which elements are compressed and which are tensioned. This truss is not a tensegrity because of the absence of both the mechanism and D and I features.
For the last case ( Figure 11d) the absence of zero eigenvalues in matrices and indicates the lack of mechanisms and self-stress states, respectively. This truss is not a tensegrity. The zero eigenvalue in L and M is responsible for the existence of the mechanism (M) (Figure 10b):   (Figure 11a to c), there is no zero eigenvalue in the matrix which means that the system lacks the mechanism. There are three zero eigenvalues in the matrix , which are responsible for the existence of three different self-stress states (S). Two of them lead to the positive eigenvalues of the matrix ( + ), and the third self-stress gives one negative eigenvalue. Several combinations of the three identified self-stress states are possible. It is impossible to clearly define which elements are compressed and which are tensioned. This truss is not a tensegrity because of the absence of both the mechanism and D and I features.
For the last case ( Figure 11d) the absence of zero eigenvalues in matrices and indicates the lack of mechanisms and self-stress states, respectively. This truss is not a tensegrity.  (Figure 12b).
For the first truss, there is no zero eigenvalue in the matrix which means that the system lacks the mechanism. There is one zero eigenvalue in the matrix which is responsible for the existence of the one self-stress state (S), as shown in Figure 11a. This structure does not satisfy the condition of discontinuous compression members (thick red line) and the inclusion of compression components inside the set of tension components, which are cables (C). In [36] this truss is called "a twodimensional tensegrity structure", but because of the absence of both the mechanism and D and I features, the truss is not a tensegrity.
A change of support (Figure 12b) caused one zero eigenvalue in the matrix responsible for the For the first three examples (Figure 11a-c), there is no zero eigenvalue in the matrix L which means that the system lacks the mechanism. There are three zero eigenvalues in the matrix M, which are responsible for the existence of three different self-stress states (S). Two of them lead to the positive eigenvalues of the matrix (K L + K G ), and the third self-stress gives one negative eigenvalue. Several combinations of the three identified self-stress states are possible. It is impossible to clearly define which elements are compressed and which are tensioned. This truss is not a tensegrity because of the absence of both the mechanism and D and I features.
For the last case ( Figure 11d) the absence of zero eigenvalues in matrices L and M indicates the lack of mechanisms and self-stress states, respectively. This truss is not a tensegrity.
For the first truss, there is no zero eigenvalue in the matrix which means that the system lacks the mechanism. There is one zero eigenvalue in the matrix which is responsible for the existence of the one self-stress state (S), as shown in Figure 11a. This structure does not satisfy the condition of discontinuous compression members (thick red line) and the inclusion of compression components inside the set of tension components, which are cables (C). In [36] this truss is called "a twodimensional tensegrity structure", but because of the absence of both the mechanism and D and I features, the truss is not a tensegrity.
A change of support (Figure 12b) caused one zero eigenvalue in the matrix responsible for the existence of one mechanism (Figure 12c). The self-stress state is the same as for previous truss ( Figure  11b), but one eigenvalue of the matrix ( + ) is negative. This means structural instability-the mechanism is not stabilized by the self-stress state. The truss is not a tensegrity.

Summary and Conclusions
The qualitative evaluation of truss structures featured in this paper was used for tensegrity form finding. New, precise definition was proposed. Characteristic features of tensegrity structures were used to classify the structures described in the literature as "pure tensegrities" or "structures with For the first truss, there is no zero eigenvalue in the matrix L which means that the system lacks the mechanism. There is one zero eigenvalue in the matrix M which is responsible for the existence of the one self-stress state (S), as shown in Figure 11a. This structure does not satisfy the condition of discontinuous compression members (thick red line) and the inclusion of compression components inside the set of tension components, which are cables (C). In [36] this truss is called "a two-dimensional tensegrity structure", but because of the absence of both the mechanism and D and I features, the truss is not a tensegrity.
A change of support ( Figure 12b) caused one zero eigenvalue in the matrix L responsible for the existence of one mechanism (Figure 12c). The self-stress state is the same as for previous truss (Figure 11b), but one eigenvalue of the matrix (K L + K G ) is negative. This means structural instability-the mechanism is not stabilized by the self-stress state. The truss is not a tensegrity.

Summary and Conclusions
The qualitative evaluation of truss structures featured in this paper was used for tensegrity form finding. New, precise definition was proposed. Characteristic features of tensegrity structures were used to classify the structures described in the literature as "pure tensegrities" or "structures with tensegrity features". Only in the case of pure tensegrities (structures with all six characteristic features) is it possible to control truss properties, such as stiffness and natural frequency, by changing the level of self-stress.
The existing form-finding methods determine the self-equilibrated configuration in the absence of externally applied loads, but they are not always define other characteristic features of tensegrity structures. The identification of the self-stress state in a truss is not sufficient to call it a tensegrity structure. From the point of view of structural mechanics, one of the most important tensegrity features is the mechanism-the uniqueness of tensegrity structures is based on the fact that self-stress states stabilize existing infinitesimal mechanisms. The absence of the mechanism brings the trusses into the group of structures with tensegrity features.
In this paper, singular value decomposition of the compatibility matrix was used to determine the states of self-stress and mechanisms. The nonlinear analysis with the use of geometric stiffness matrix, which takes into account the impact of the self-stress state was used to establish whether the mechanism was infinitesimal or finite. The other characteristic features of tensegrity structures, that satisfy the condition of discontinuous compression members and the inclusion of compression components inside the set of tension components (cables) were determined indirectly by the identified self-stress state.
The compatibility matrix and the stiffness matrix with the effect of self-equilibrated forces were determined using the finite element method. The geometrical and material characteristics were not taken into account because the characteristic features of tensegrity structures do not depend on them.
The proposed method of evaluation allowed a precise understanding of this unusual structure and a way to specify the correct use of this term in these applications. The results are summarized in Table 1. Some 2D pure tensegrity structures (according to the proposed definition) exist, but most of the trusses are incorrectly named in the literature. The term "tensegrity" is misused and overused in the literature. Three pure tensegrities can be distinguished, with one infinitesimal mechanism stabilized by one self-stress. Three other more complex structures are pure tensegrities with two infinitesimal mechanisms stabilized by one self-stress.      In conclusion, tensegrity is overused in the literature on the subject, in particular, in the field of 2D models, where intuitive understanding of the specific features should be the easiest. Examples of pure tensegrities shown in Figures 5a and 6, and their extensions, can be recommended as 2D computational models for use in the coverings for stadiums and halls. The schemes from Figure 8 can be the basis for the construction of 2D metamaterials with specific mechanical features. The unit cells presented in Figures 5a and 8 are recommended as metamaterial unit cells. It should be noted that tensegrity trusses with computational schemes presented in Figures 5a, 6 and 8 can be successfully applied as deployable or smart engineering structures with monitoring, self-control, self-repair, and active control.
Author Contributions: Introduction was prepared by P.O. and W.G. Definition of tensegrity and mathematical description were written by P.O., J.K. and W.G. Results were obtained by P.O. and J.K. The analysis of the results and conclusions were written by P.O. and W.G.
Funding: This research received no external funding.