Inverse Eigenvalue Problem For Mass-Spring-Inerter Systems 111The work was supported by the National Natural Science Foundation of China under Grants 11925109 and 11688101.
Abstract
This paper has solved the inverse eigenvalue problem for “fixed-free” mass-chain systems with inerters. It is well known that for a spring-mass system wherein the adjacent masses are linked through a spring, the natural frequency assignment can be achieved by choosing appropriate masses and spring stiffnesses if and only if the given positive eigenvalues are distinct. However, when we involve inerters, multiple eigenvalues in the assignment are allowed. In fact, arbitrarily given a set of positive real numbers, we derive a necessary and sufficient condition on the multiplicities of these numbers, which are assigned as the natural frequencies of the concerned mass-spring-inerter system.
1 Introduction
Natural frequency, an inherent attribute of mechanical vibration systems, has attracted wide attention for its importance. In particular, purposefully allocating the natural frequencies to some pre-specified values provides an effective way to induce or evade resonance (see [2], [23]). This naturally raises the inverse eigenvalue problem (IEP), that is, to construct a vibration system whose natural frequencies, or mathematically known as eigenvalues, are given beforehand.
A well-known result on this problem is due to [8] and [17], which is addressed for mass-spring systems. Observe that in such a basic system, the adjacent masses are linked merely by a spring. Therefore, the IEP turns out to be the construction of a Jacobi matrix with its eigenvalues being assigned to a set of specified positive numbers. Borrowing the tools for Jacobi matrices, [8] and [17] assert that the IEP is solvable if and only if the given positive eigenvalues are distinct. Later on, a various of inverse problems on Jacobi matrices and Jacobi operators are investigated as well. But when dampers are taken into consideration, the matrices associated with masses, spring stiffnesses and damping coefficients in mass-spring-damper systems are no longer Jacobi matrices and the quadratic inverse eigenvalue problem (QIEP) is put forth. Nevertheless, most of the literatures on mass-spring-damper systems focus on distinct eigenvalues assignment [3, 13, 1, 16].
Interestingly, the IEP admits multiple eigenvalues, if we introduce a mechanical element called the inerter. This new mechanical device can simulate masses by changing inertance. It was theoretically first studied by [20], completing the analogy between electrical and mechanical networks (see Figure 1.2). Through physical realization, inerters have been applied to many engineering fields such as vibration isolators, landing gears, train suspensions, building vibration control, and so on [6, 22, 14, 21, 11]. In a mass-spring-inerter system, the neighbouring masses are linked by a parallel combination of a spring and an inerter. As a starting point, we restrict our interest in this paper to “fixed-free” systems. The term “fixed-free” means one end of the mass-chain system is attached to the ground while the other end is hanging free, as shown in Figure 1.2.
The free vibration equation of such a mass-spring-inerter system is described by
where and
(1) |
(2) |
(3) |
Here, real numbers , , for stand for the masses, spring stiffnesses and inertances. Unlike mass-spring systems, the well-studied Jacobi matrix theory cannot illuminate the IEP for mass-spring-inerter systems since the inertial matrix in (3) is a tridiagonal matrix. Recently, [10] found that inerters render the multiple eigenvalues possible for a mass-chain system. It showed that the multiplicity of a natural frequency must fulfill . Beyond that, little is known for the multiple eigenvalue case.
The purpose of this paper is to solve the IEP for mass-spring-inerter systems, where the eigenvalues are arbitrarily specified to positive real numbers. We deduce a necessary and sufficient condition for this assignment on the multiplicities of the given numbers. With the proposed critical criterion, the set structure of the given real numbers will be intuitively clear for the natural frequency assignment. Our construction further implies that masses of the system can be arbitrarily fixed beforehand for the assignment, where is the amount of the distinct assigned eigenvalues. More precisely, our construction is carried out by only adjusting massess, spring stiffnesses and inertances. It degenerates to the claim that, if the pre-specified eigenvalues are all distinct (), the IEP can be worked out by recovering and , whereas is fixed arbitrarily. This claim is exactly the main result of [20], which demonstrates an advantage of using inerters in the mass-fixed situation. Unfortunately, not all the natural frequency assignments are realizable by merely adjusting spring stiffnesses and inertances. An example of five-degree-of-freedom system in this paper shows that there exist some restrictive relationships between masses and given eigenvalues.
The organization of this paper is as follows. In Section 2, we state the main result by deducing a necessary and sufficient condition of the IEP for mass-spring-inerter systems, while the proofs are included in Sections 3 and 4. Conclusions are drawn in Section 5.


2 Main Result
The natural frequencies of a mass-spring-inerter system are completely determined by the eigenvalues of matrix pencil , where are defined by (1), (2) and (3), respectively. So, with a slight abuse of language, we will not distinguish the term “eigenvalues” from the “natural frequencies” in this article. We now raise our problem.
Problem 1. Arbitrarily given a set of real numbers , is it possible to recover matrices in (1), (2) and (3) by choosing , and for , so that the eigenvalues of matrix pencil are exactly , ?
Both [10] and [17] offered a positive answer to Problem 1 for the special case where the eigenvalues are all distinct. But the general situation should involve multiple eigenvalues, which is covered by the following theorem.
Theorem 2.1.
Remark 2.1.
Theorem 2.1 completely solves Problem 1 by providing the critical criterion (5). As indicated later (see Proposition 4.1 for details), when (5) holds, the recover of the relevant matrices allows a total of masses being taken arbitrarily, where is the number of distinct eigenvalues given beforehand. Particularly, for , each mass can be taken any fixed quantity in advance, as proved in [10]. However, when , it is generally impossible to achieve the natural frequency assignment with all the masses arbitrarily given. Example 2.1 suggests a restrictive relation between the masses and eigenvalues.
3 Proof of the necessity of Theorem 2.1.
This section is devoted to proving the necessity of Theorem 2.1, which is relatively easier than the argument for sufficiency. We begin by expressing in terms of a recursive sequence of polynomials. First, write
For , let , and be some matrices defined analogously as and in (1)–(3), respectively, but with order instead of . Next, denote as the determinant of , . Let and the leading principal minor of of order , where . So, to calculate , we only need to treat .
Remark 3.1.
For each , since , the Gershgorin’s circle theorem indicates that both and are positive definite matrices and so do their leading principal submatrices. Then, it follows that the roots of and are all real and positive (see [8, Theorem 1.4.3] ).
To facilitate the subsequent analysis, we introduce the following definition.
Definition 3.1.
Let and be two polynomials with degree , where . Suppose and both have distinct real roots, which are denoted by and , respectively. We say , if their leading coefficients are of the same sign and
The proof depends on a simple observation below.
Lemma 3.1.
The polynomials and satisfy
(7) |
with and .
Proof.
By the definition of , it is trivial that . For , expanding the leading principal minor of of order by cofactors of the th row shows . Furthermore, the expansion of by cofactors of the th row yields
as desired.
Lemma 3.2.
Let and be two polynomials that , then for any , .
Proof.
Since , considering Definition 3.1, we let for some and let and be the roots of and , respectively. Clearly,
Without loss of generality, assume the leading coefficients of and are both positive. Then,
This implies that for each , there is a root of falling in interval . Observing that the degree of is , the result follows immediately.
We present an important property enjoyed by sequence .
Lemma 3.3.
Suppose for some ,
(8) |
then
Moreover,
(i) if and , then
and ;
(ii) if and , then
and ;
(iii) if , then and .
Proof.
First, according to the definitions of , it is apparent that and for each . Let for some integer . Recalling (8), denote the roots of and by and , respectively. These roots fulfill
(9) |
In addition, the second equation of (7) implies . Now, we prove this lemma by discussing three cases.
(i) and . For this case, there exists some such that . We shall evaluate the sign of at . In fact, according to the second equation of (7),
Note that the leading coefficient of is positive, the definitions of and read
(10) |
which, together with (9), yields that for ,
and
Therefore,
This means that for each , there exists exactly one root of between and , and hence .
Moreover, , so
and
Applying Lemma 3.2 to the first equation of (7), we thus deduce
(12) |
Now, becasue of , it follows from (12) that
(ii) and . We first assume for some and evaluate the sign of at points and . Since Lemma 3.1 shows
Similarly, by Lemma 3.1, (9) and (10),
So, for each with , there is a root of falling in , and the rest two roots of lie in and , respectively. This infers by noting that . Hence, and applying Lemma 3.2 to the first equation of (7), one has
Then, and consequently
As for the situations where , and , an analogous treatment can be employed.
(iii) . We also first calculate the sign of at the roots of and . As before,
(16) |
Now, and the leading coefficient of is positive, if number is sufficiently large, it is evident that and
So, each interval in contains exactly one root of , which concludes
It is ready to prove the necessity of Theorem 2.1. To this end, we introduce some notations. Let be a polynomial whose roots are all real and for every . Denote as the multiplicity of root and for a real number , define .
The proof of the necessity of Theorem 2.1: First, in view of (4), we know that are the distinct roots of with multiplicities . To proceed the argument, note that Lemma 3.1 gives , then by using Lemma 3.3,
(21) |
Particularly, it turns out that all the roots of are distinct. So for all and then
(22) |
Now, by (21) and Lemma 3.3, for each ,
which yields that for all ,
(26) | |||||
On the other hand, given and , if and , Lemma 3.3 (i) indicates . Moreover, is a root of , it follows that . Otherwise, by virtue of Lemma 3.3, at least one of the following cases will happen:
(i) if and , then
(ii) if and , then
;
(iii) if , then .
All the above three cases lead to . So,
(29) |
Clearly, . By (26) and (29), it can be derived inductively that
Together with (22), the above inequality shows that for each ,
which completes the proof.
4 Proof of sufficiency of Theorem 2.1.
The sufficiency of Theorem 2.1 is a direct consequence of the following proposition.
Proposition 4.1.
Given real numbers and a polynomial with and , if for each , then there exist some , , , and distinct indices such that and
We now begin the construction of the required mass-chain system for Proposition 4.1. That is, to find a sequence of and indices such that and , , where . Taking account to [10, Theorem 4], we take . Moreover, denote , where
(30) |
Evidently, . We reorder the elements of by .
An important observation of Section 3 is that (22) implies every element in is equal to some . This could help us to design a rule to determine which . It is the key idea of our proof, so we offer an example to elaborate it.
Example 4.1.
Take and in Proposition 4.1. Then, we get sets as shown in Figure 4.4. Note that by (2)–(3), for all , and are located in the secondary diagonals of matrices K and B, respectively. We now introduce a matrix and assign the elements in the secondary diagonal the values taken from sets (see Figure 4.4). Specifically, for , we treat the first elements of by skipping and letting for . So, , as illustrated in Fig. 3. Repeat this procedure for elements with . For each with , we assign a value equivalent to . As a result, , , , .


In general, the rule to determine is summarized as follows:
(31) |
The proof of Proposition 4.1 thus will be completed in three steps.
Step 1: For each with and , we assign a value taken from in the light of (31).
So, the cardinal of set is .
Step 2: Rewrite the elements of by . Let for .
Step 3: Based on the above steps, we compute for each . Then, take some suitable parameters , so that if ,
(34) |
otherwise, for ,
(37) |
The construction of will be achieved by an induction method from to .
To proceed the proof, we first derive some technical lemmas.
Lemma 4.1.
Let and be two polynomials of degree that , where
,
and . Then, has real roots satisfying for and . Moreover, the following two statements hold:
(i) when , for any , if
(38) |
then
(39) |
(ii) when , for any , (39) holds provided that
(40) |
Proof.
Note that indicates and for each ,
which means has a root in . Further, since and , holds for all sufficiently large . On the other hand,
so has a root in . Clearly, .
The subsequent parts focus on Step 3 of the construction, whose key idea is to select some appropriate candidates for the roots of . We set these roots as , , where ,
(42) |
Next, define
(43) |
as well as
(44) |
Remark 4.1.
The above series of constants are repeatedly used in the next two lemmas (Lemmas 4.2–4.3), whose proofs are contained in Appendix C. Both the two lemmas concern the following polynomials
(49) |
whose roots and satisfy
(50) |
and for each ,
(51) |
Lemma 4.2.
Let and be two polynomials defined by (49)–(51) with
in (50).
For any given constants satisfying and
, the following two statements hold.
(i) If , then there exist two monic polynomials and with distinct roots and , respectively, satisfying and for ,
(52) |
In addition, for some numbers and with , and fulfill
(55) |
(ii) If , then there are some numbers , and with such that polynomials and satisfy equation (55).
Lemma 4.3.
Given with and
, the following two statements hold.
(i) For polynomials and defined by (49)–(51) with
in (50), there exist two monic polynomials and with distinct roots and , respectively, such that and (52) holds for all . In addition, for some numbers and with , and fulfill
(58) |
(ii) For with and , there are some numbers and with such that polynomials and satisfy equation (58).
Lemma 4.4.
Let and be two polynomials and be two numbers satisfying
(59) |
Then, there exist some monic polynomials , and some sequences of numbers , such that for each , the following two properties hold:
(i) and ;
(ii) if , then
(62) |
otherwise, for ,
(66) |
where and is defined by (30).
Proof.
For , we construct a series of numbers , and polynomials , on the basis of , and , , according to the following strategies:
We shall use the induction method to show that either strategy a) or strategy b) can be implemented for each . First, let . Observe that , it is easy to compute
Hence, . In addition, since (59) and (69) yield and , by applying Lemma 4.2(i) with , , , , , we can find some numbers with and
and two monic polynomials , such that (62) holds. So, strategy a) applies and both (i) and (ii) are true for these .
Now, assume that we have constructed the required , , and for some by following either strategy a) or strategy b), so that properties (i) and (ii) hold for . Considering Lemmas 4.2 and 4.3, we write with and with , . Here, for each ,
(72) |
Furthermore, if , then for each and ,
(75) |
Recall that and , (75) implies that for ,
(76) |
Moreover, by (75) again,
then a straightforward calculation leads to
(77) |
Note that by (48) in Remark 4.1,
so by virtue of (76) and (77),
(78) |
We remark that no matter strategy a) or b) applies, it always infers
and hence
(80) |
We now verify that at least one of the strategies a) and b) is valid for . It is discussed by two cases.
Case 1: .
Because of (72), we estimate directly by
Now, , it thus gives . So combining (78) and (80), it shows that the assumptions of Lemma 4.2 are fulfilled. Therefore, Lemma 4.2 is applicable and strategy a) works. We thus conclude properties (i) and (ii) for by Lemma 4.2 and
Case 2: for some . If , then , and hence is constructed by strategy b). In view of Lemma 4.3, the maximal root of is equivalent to , which shows
If , then is constructed by strategy a). Observe that
so by (78), . On the other hand, if we suppose for some , then because of (30). As a consequence,
(81) |
and hence
(82) |
We thus conclude (82) is always true in strategy a). By (80),
This combining with (78) and (80) verifies the assumptions of Lemma 4.3. Hence, Lemma 4.3 applies and strategy b) can be implemented. So, properties (i) and (ii) hold for due to Lemma 4.3 and
The induction is completed.
Proof of Proposition 4.1.
Let be defined in Lemma 4.4, then we can construct a series of numbers , and some monic polynomials , . Note that , set
(85) |
We shall see that meet the requirement of Proposition 4.1.
In fact, define a sequence of polynomials as follows:
(89) |
and let
(90) |
Clearly, The rest of the proof is to check whether satisfy the recursive formula in Lemma 3.1.
5 Concluding remarks.
The emergence of inerters in engineering brings some new phenomena in the study of inverse problems. Particularly, it enables a mass-chain system to possess multiple eigenvalues. This paper has solved the IEP for the “fixed-free” case, where the real numbers for eigenvalue assignment can be taken arbitrarily positive. Another common situation is that the both ends of the system are fixed at the wall. For such “fixed-fixed” systems, the construction cannot follow readily from the method developed in Section 4. To address this issue, a more detailed analysis is required and it would be our next work.
Appendix A Proof of Example 2.1
For each , let and be two monic polynomials such that
where and are the leading coefficients of and , respectively. Apparently, and
(94) |
Moreover, denote and as the roots of and . Define
We first assert that cannot contain any adjacent natural numbers. Otherwise, suppose there is a number such that . So, and by Lemma 3.3,
(95) |
Since the definition of indicates , then by (95),
Hence, , which contradicts to (94) that .
We in fact have derived , which together with Lemma 3.2 implies and . On the other hand, (22) and (26) in Section 3 infer , and hence . Observe that or , it then follows that . So, and by Lemma 3.3(i), and . For the roots of , they can be discussed similarly by using Lemma 3.3 repeatedly. Indeed, and lead to and , respectively. So,
In addition, indicates and
(97) |
At last, since , . By and (26), which immediately gives . So, . If , then (97) infers that . Therefore,
We thus summarize
(102) |
and
(105) |
This means and as a consequence, by Lemma 3.3,
(108) |
and
(111) |
With the above properties, we can also present the relationship between the roots of and for . In fact, when , and by (21), (102) and (105), . Then, (21), Lemma 3.2, (102) and (111) imply
(112) |
Now, we prove (6) by using reduction to absurdity. Suppose
(113) |
First, it is evident that for each , (102), (108) and (111) yield
(114) | |||||
In particular, when , (21) further implies
(115) |
Note that by comparing the leading coefficients of the polynomials in (108)–(111), we assert that for each with ,
(116) |
Finally, let us complete the proof by considering the following four cases.
Case 1: . First, we estimate for . Let , then by (102)–(105). Therefore, noting that , (21) and (116) lead to
(117) | |||||
So, if , the above inequality reduces to
(118) |
We next treat . Since (21) and (102) imply and , then
which is equivalent to Further, by (111) again, , so (112) shows
(119) | |||||
(120) |
As for , similar to (119), we compute by (112) that
(121) | |||||
As a result, by (118), (120) and (121), the following inequality always holds:
and thus (114) and (115) yield
(122) |
We proceed to the calculation of . Note that , then analogous to (117), we can prove which together with (114) and (115) leads to
(123) |
Now, combining (122) and (123), (113) derives
This is no other than
Consequently,
which contradicts to
the definition of in (113).
Case 2: . Then, (111) infers
In this case, and then (111) becomes
(126) |
Since and , comparing the leading coefficients of the polynomials in (126), we calculate
Consequently, by (21) and (112),
(127) |
and hence
Therefore, by (113),
which infers . This is impossible due to the definition of in (113).
Case 3: and . Then, and a similar argument as (127) indicates
It contradicts to (113) that .
Case 4: and . Then, (127) holds and according to (114),
By (117)–(120) and (123) in Case 1, we can demonstrate
The rest of the proof thus keeps the same as that in Case 1.
Appendix B Proof of Remark 4.1
Appendix C Proofs of Lemmas 4.2–4.3
Proof of Lemma 4.2.
Let . We first take a number satisfying
(128) |
Such indeed exists because of (49)–(51), which yield and
Then, let
(134) |
we shall show that all the above defined numbers and polynomials fulfill our requirements.
Observe that and by (50)–(51),
then . As a result, by (128),
(135) |
Since (128) indicates , (135) means is a well-defined monic polynomial of degree . Recall that and , applying Lemma 4.1 to polynomials and shows
(136) |
Therefore,
(137) |
which give
(138) |
Consequently, by (134) and (135),
(139) |
So far, we have verified that and satisfy the condition of Lemma 4.2. For these numbers, the first equality of (55) follows directly from (134). Considering (136), if the second inequality of (52) holds when , then in (134) will be the exact polynomial desired for both statements (i) and (ii).
Next, we check in (134).
The definition of infers
,
hence (138) implies that is a well-defined monic polynomial. We discuss this part by considering two cases.
(i) . In this case, fulfills the second equality of (55) by (134).
(ii) . Taking account to (137) and (138), we apply
Lemma 4.1 to polynomials and , then
(140) |
so the roots of are distinct. Further, we can verify the second equality of (55) from (134) again.
Now, it remains to show the second inequality of (52) for when . Combining (136) and (140), it gives . Fix an index , we compute
(141) |
Note that by (137) and by (136) and (140), (141) immediately leads to
(142) |
We are going to employ Lemma 4.1 to estimate term in (142). For this, denote . Recall that , then by (51), for each ,
which means . As a result, by (44),
So, by applying Lemma 4.1 to polynomials and , we conclude
(143) |
Then, , which together with (51) indicates
(144) | |||||
Next, we deal with in (142) for . As a matter of fact, by (51), for any ,
(145) | |||||
Proof of Lemma 4.3.
(i) Let and set
(152) |
By (49) and (51), it is clear that . Observe that
then is a well-defined monic polynomial of Now, (51) indicates , by applying Lemma 4.1 to polynomials and , it follows that the first roots of satisfy
(153) |
Moreover, the definition of in (152) shows
which yields Hence, .
Next, let
(158) |
where
(161) |
with
Therefore,
Note that , then (158) gives . As a result, by ,
So, is a well-defined monic polynomial of degree . Since (153) means that , by applying Lemma 4.1 to polynomials and , we deduce and
(162) |
Observe that plugging (152) and (158) into (58) immediately shows the validity of (58).
At last, we show the second inequality of (52) for . Fix an index . By plugging into (58), we obtain which equals to
(163) |
Note that for all because of (153). As for , by (51) and (153),
which together with (48) yields
Furthermore, since , (42), (48) and (51) lead to
Consequently, it follows from (162) and (163) that
Now, we prove for each . If ,
By applying Lemma 4.1 to polynomials and , it infers
When ,
and applying Lemma 4.1 to polynomials and shows
So both cases result in
which implies that for each . Then,
This finishes the proof of statement (i).
(ii) Let , ,
and ,
we can directly compute that and satisfy equation (58).
References
- [1] Bai Z 2007 Symmetric tridiagonal inverse quadratic eigenvalue problems with partial eigendata Inverse Problems 24 015005
- [2] Beeby S P, Tudor M J and White N M 2006 Energy harvesting vibration sources for microsystems applications Meas. Sci. Technol. 17 175–195
- [3] Cai Y F and Xu S F 2009 On a quadratic inverse eigenvalue problem Inverse Problems 25 085004
- [4] Chen M Z Q, Hu Y, Huang L and Chen G 2014 Influence of inerter on natural frequencies of vibration systems J. Sound Vib. 333 1874-1887
- [5] Chu M T 1998 Inverse eigenvalue problems SIAM Rev. 40 139
- [6] Dong X, Liu Y and Chen M Z Q 2015 Application of inerter to aircraft landing gear suspension Proceedings of the 34th Chinese Control Conference, Hangzhou, China 2066-2071
- [7] Gladwell G M L 1997 Inverse vibration problems for finite-element models Inverse Problems 13 311
- [8] Gladwell G M L 2005 Inverse problems in vibration second ed. Springer Science Business Media
- [9] Gladwell G M L 2006 Minimal mass solutions to inverse eigenvalue problems Inverse Problems 22 539
- [10] Hu Y L, Chen M Z Q and Smith M C 2018 Natural frequency assignment for mass-chain systems with inerters Mech. Syst. Signal Process. 108 126–139
- [11] Lazar I F, Neild S A and Wagg D J 2014 Using an inerter-based device for structural vibration suppression Earthq. Eng. Struct. Dynam. 43 1129–1147
- [12] Liang Y, Yamaura H and Ouyang H 2017 Active assignment of eigenvalues and eigen-sensitivities for robust stabilization of friction-induced vibration Mech. Syst. Signal Process. 90 254–267
- [13] Lin M M, Dong B and Chu M T 2010 Inverse mode problems for real and symmetric quadratic models Inverse Problems 26 065003
- [14] Matamoros-Sanchez A Z and Goodall R M 2015 Novel mechatronic solutions incorporating inerters for railway vehicle vertical secondary suspensions Veh. Syst. Dynam. 53 113–136
- [15] Mottershead J E and Ram Y M 2006 Inverse eigenvalue problems in vibration absorption: Passive modification and active control Mech. Syst. Signal Process. 20 5–44
- [16] Nylen P 1999 Inverse eigenvalue problem: existence of special mass-damper-spring systems Linear Algebr. Appl. 297 107–132
- [17] Nylen P and Uhlig F 1997 Inverse eigenvalue problem associated with spring-mass systems Linear Algebr. Appl. 254 409–425
- [18] Olsson P and Lidström P 2007 Inverse structural modification using constraints J. Sound Vib. 303 767–779
- [19] Ouyang H and Zhang J 2015 Passive modifications for partial assignment of natural frequencies of mass-spring systems Mech. Syst. Signal Process. 50 214–226
- [20] Smith M C 2002 Synthesis of mechanical networks: the inerter IEEE Trans. Autom. Control 47 1648–1662
- [21] Wang F C, Hong M F and Chen C W 2010 Building suspensions with inerters Proc. IMechE Part C: J. Mech. Eng. Sci. 224 1605–1616
- [22] Wang F C and Liao M K 2010 The lateral stability of train suspension systems employing inerters Veh. Syst. Dynam. 48 619–643
- [23] Yu Y H and Naganathan N G 2001 A literature review of automotive vehicle engine mounting systems Mech. Mach. Theory 36 123–142