Abstract

Compliant joints have significant advantages compared to rigid-body hinges due to a monolithic design and the absence of friction, which prevents effects like wear, backlash, and stick-slip behavior. However, the loading capability is often limited and the support stiffness generally decreases during rotation. A new design principle called closed form pressure balancing has been proposed as a solution to improve these limitations. By using an incompressible fluid as the main compliant element, the support stiffness becomes independent of rotation and buckling no longer limits the loading capability. This work analyzes the fundamental working principle behind closed form pressure balancing and introduces a 2D design model to determine stiffness properties. The design model is validated with a finite element model and used to construct an optimization strategy for optimum joint performance. Additionally, a conversion model and some practical considerations are presented for the transition to a 3D design model.

1 Introduction

Compliant hinges have introduced elastic deformation of flexural elements as an alternative to sliding and rolling contacts in rigid-body hinges. Due to their monolithic design, significant advantages are the absence of wear, backlash, stick-slip behavior, and the need for lubrication, allowing high repeatable motions with high precision [13]. General disadvantages include a limited range of motion due to stress concentrations [3] and a moving center of rotation (CoR), called axis drift [4]. Combinations of leaf flexures have been proposed as solutions to these limitations, like the cross-axis flexure [5], cartwheel hinge [3], X2-joint [6], butterfly hinge [7], trapezoidal hinge [8], and CR joints [4]. Other disadvantages, also applicable to these examples, are a limited support stiffness and risk of buckling when the hinges are loaded in compression [9]. Although these effects depend on the design, they often worsen during rotation as leaf flexures lose their support stiffness when deflected [10].

Specific design methods to improve support stiffness and load capability have been proposed. These include the use of inversion and isolation [9], compliant aided rolling contacts [11,12], pre-curved flexures [13], torsional stiffeners [10], and folded leaf springs [14]. However, most methods reintroduce contact surfaces, increase the kinematic complexity, or still suffer from a decreasing performance during rotation. Additionally, as these solutions are based on the use of solid leaf flexures, the design trade-off remains between thinner flexures for a lower rotation stiffness and larger range of motion, and thicker flexures for a higher support stiffness.

This relation between support and rotation stiffness limits the application range of compliant hinges. In positioning mechanisms for example, the decreasing support stiffness during rotation causes a relatively large footprint compared to the effective workspace [15]. Similarly, biomechanical applications for the replacement of joints loaded in compression, like ankles, hips, or spinal discs, remain difficult to design [1618]. Other examples that can benefit from axial stiff compliant joints are the connection between a piston and eccentric running track or universal joints in parallel kinematic positioning systems, like large ground based telescopes [1921]. Improving the support stiffness over a larger rotation could thus increase the range of applications.

Recent work [22] has proposed a new design principle that could improve the support stiffness, utilizing an enclosed incompressible fluid. A prototype has shown an increase in stiffness ratio of the axial support stiffness over the rotation stiffness [22], although the characteristic behavior is still unknown. This current paper presents simplified analytical models to estimate the behavior of this new conceptual principle named closed form pressure balancing. First, the working principle and embodiment are elaborated. An analytical 2D design model to determine characteristic stiffness behavior is presented and verified with a finite element method (FEM) model. The models are then used to analyze the movement of the CoR, study the influence of various design parameters, and develop an optimization method for designing. Finally, 2D to 3D conversion models are introduced that are used to translate the 2D results to a 3D design.

2 Design Model

The design principle analyzed in this work is the use of an incompressible fluid as main compliant element to be loaded under compressive loads [22]. Using a fluid fundamentally changes the working principles behind the axial and rotation stiffnesses. The liquid state allows fluid molecules to freely reorganize and thereby accommodate a change in shape of the volume. Neglecting inertial and viscous effects, this theoretically results in zero rotation stiffness. However, incompressibility makes the fluid resist any volume changes when loaded in compression, resulting in a significant axial stiffness. The first major advantage compared to flexures is that this working principle of axial stiffness is insensitive to the shape of the fluid volume, theoretically decoupling the axial and rotation stiffnesses. This could prevent a decrease in axial stiffness during rotation and eliminates the trade-off between a high axial stiffness and low rotation stiffness. This characteristic behavior, caused by the closed volume of incompressible fluid, is referred to as “closed form pressure balancing” [22]. However, performance will be influenced by the shape and stiffness of the elastic solid required to contain the fluid. The influence of this elastic body is elaborated on through the use of a simplified 2D analytical model.

2.1 Conceptual Design.

A cell is thus needed that optimally transfers the properties of the fluid to a physical joint. In this work, a cell with circular enclosures of thickness t0 is proposed as shown in Fig. 1. When vertically compressed, fluid is pushed sideways into the enclosure, causing an internal pressure. During rotation and shear, however, no pressure buildup is expected as the volume only needs to change shape, making the rotation and shear stiffnesses fully dependent on the enclosure’s bending stiffness. The circular shape ensures an even distribution of stresses in the material under internal pressure, whereas the curved geometry allows for a low bending stiffness and thus a low rotation stiffness. The rigid connectors are added to clamp the horizontal segments as they potentially provide additional compliance to the enclosures, which would reduce axial stiffness. The second major advantage of this embodiment compared to flexures is that compressive loads are converted to a tensile load in the enclosures, eliminating the risk of buckling. An incidental drawback, however, is a limited shear stiffness due to the low bending stiffness of the enclosure, which is further addressed in Sec. 5.3.

Fig. 1
Embodiment of the pressure balanced compliant joint (left) and the enclosure shape for different angles of θ (right)
Fig. 1
Embodiment of the pressure balanced compliant joint (left) and the enclosure shape for different angles of θ (right)
Close modal

The angle of the circular segment, indicated with θ in Fig. 1, is expected to influence both stiffness and range of motion. For smaller θ the enclosure length increases, which lowers the elongation stiffness and bending stiffness of the enclosures and thus decreases both axial and rotation stiffnesses of the joint. Furthermore, larger arc segments allow larger joint rotations while bending with distributed compliance. As shown in the top right of Fig. 1, smaller arc segments transition to lumped compliance or even require elongation for anti-clockwise rotations of the upper bar to bridge the gap between the endpoints of both connectors, which increases rotation stiffness. Smaller values of θ thus seem favorable, although it reduces axial stiffness. However, as the bending and elongation stiffnesses of the enclosures scale with t03 and t01 respectively, reducing t0 is considered most efficient to minimize rotation stiffness. Therefore, θ has been increased to π/2(rad) to improve axial stiffness while limiting the increase in rotation stiffness due to lumped compliance.

The following sections present two analytical models to predict axial, rotation, and shear stiffnesses. All models assume a linear elastic material and out-of-plane width b of 1 (m). The fluid is assumed incompressible, requiring its initial area A0 to remain constant. The horizontal parts of the elastic cell are assumed rigid by clamping of the rigid connectors in Fig. 1.

2.2 Axial Stiffness.

The axial compressive stiffness model uses two additional simplifications. First, the enclosures are assumed to maintain a circular shape due to the internal pressure buildup. Second, the enclosure’s bending stiffness is neglected compared to its elongation stiffness for inflation as only thin enclosures are considered to obtain a low rotation stiffness. This includes a free rotation of the enclosures at their connection point to the horizontal segments, as illustrated in Fig. 2.

Fig. 2
Enclosure showing the assumption of free rotation (left) and the inner wall geometric description in a deformed configuration for a compression displacement dh (right)
Fig. 2
Enclosure showing the assumption of free rotation (left) and the inner wall geometric description in a deformed configuration for a compression displacement dh (right)
Close modal
The modeling approach consists of two parts. The first part is a geometrical description of the required deformation of the inner wall of the enclosure, as schematically shown in the right of Fig. 2. A vertical displacement of the top causes fluid to move into the enclosures. The fluid surface area Ae in a single enclosure for a compression displacement dh is given by Eq. (1), where L0 is the length of the horizontal segments and h0 the initial height of the fluid as in Fig. 1.
Ae=(A0L0(h0dh))/2
(1)
The new geometry of the inner wall follows from two constraints. First, the opening of the inner wall should match with the center volume as given by Eq. (2), where Ra is the radius and θ the opening angle of the inner wall shown in Fig. 2. Second, the volume enclosed by the inner wall should comply with Ae, described in Eq. (3). Together this results in the parameters Ra and θ for every compression dh, from which the inner wall length c can be obtained using Eq. (4).
h0dh=2Rasin(θ/2)
(2)
Ae=πRa20.5Ra2(θsin(θ))
(3)
c=(2πθ)Ra
(4)
The second part of the model uses the derived geometry of the inner wall to determine the axial stiffness of the joint. The steps described in this part are schematically depicted in Fig. 3. Each enclosure is considered as a segment of a hollow cylinder subjected to an internal pressure, for which analytical models exist that describe its linear elastic expansion. As derived by Ref. [23], the radial displacement u of a material point at a center-distance r in a cylinder, subjected to an internal pressure p, is described by Eq. (5), where Sa and Sb are respectively the undeformed inner and outer radii of the cylinder, E the modulus of elasticity, and v the Poisson’s ratio.
u(Sa,Sb,p,r)=p(1+v)E(Sb2Sa2)((12v)Sa2r+Sa2Sb2r)
(5)
To implement Eq. (5), two additions to the model are required. First, it is noted in the left of Fig. 2 that the endpoints of the enclosure move closer together during compression. To account for this inward movement, the inflation of the enclosures is approached in steps, where in each step Eq. (5) is applied to the updated geometry. The second addition relates to the requirement of Eq. (5) that the initial internal pressure is zero. By using the stepwise approach for inflation, the initial pressure in subsequent modeling steps will be nonzero. Therefore, the geometry at the start of each step is projected back to a state of zero pressure before Eq. (5) is applied. Combined, this results in a sequence of computations for each inflation step. To determine the force Fi+1 for a compression dh,i+1 and the corresponding secant stiffness, this sequence is as follows, as numbered in Fig. 3:
  1. Using the geometrical analysis, the updated (i + 1) inner radius Ra,i+1, angle θi+1, and length of the inner wall ci+1 are determined based on the compression dh,i+1.

  2. With the pressure pi and both inner and outer radii Ra,i and Rb,i known from the previous step (i), Eq. (5) is applied twice to find the geometry of the enclosure from the previous step, with angle θi, for the case its internal pressure would be zero. This results in Eqs. (6) and (7) where both equations are solved simultaneously for Sa and Sb, respectively the inner and outer radii of the enclosure in step (i) under zero internal pressure.
    Ra,i=Sa+u(Sa,Sb,pi,Sa)
    (6)
    Rb,i=Sb+u(Sa,Sb,pi,Sb)
    (7)
  3. For the depressurized enclosure obtained in step 2, the radial displacement ua required to increase the inner wall length to ci+1 is determined using Eq. (8). The pressure pi+1 to achieve this displacement ua from the depressurized state can then be calculated with Eq. (5), using r = Sa. Substituting this pi+1 again in Eq. (5) using r = Sb yields the displacement ub of the outer wall.
    ua=ci+12πθiSa
    (8)
  4. From the displacements ua and ub, the enclosure thickness in the updated step (i + 1) is determined to find the new outer radius Rb,i+1 as described in Eq. (9).
    Rb,i+1=Ra,i+1+((Sb+ub)(Sa+ua))
    (9)
  5. The pressure pi+1 is combined with the elastic force in the enclosures to determine the compression force. As derived by Ref. [23], the tangential stress σθθ as function of position r is described by Eq. (10). Integrating this equation over the wall thickness in depressurized condition, defined by Sa and Sb, results in Eq. (11). Combined with the internal pressure on the top horizontal segment of the joint, the compression force Fi+1 is given by Eq. (12), where the elastic force Fθθ is added as a tensile force on the top horizontal segment as shown in Fig. 3.
    σθθ=p(Sa2Sb2Sa2+Sa2Sb2(Sb2Sa2)r2)
    (10)
    Fθθ,i+1=pi+1(bSb2Sa2)(SaSb2Sa3)
    (11)
    Fi+1=pi+1(bL0)+2cos(θi+1/2)Fθθ,i+1
    (12)
    The axial secant stiffness Ky then follows from Eq. (13).
    Ky=Fi+1dh,i+1
    (13)
Fig. 3
Flowchart for the calculation of a single inflation step in the enclosure and the corresponding compression force
Fig. 3
Flowchart for the calculation of a single inflation step in the enclosure and the corresponding compression force
Close modal

2.3 Rotation and Shear Stiffnesses.

For rotation or horizontal shearing of the upper rigid connector, no significant fluid flow from the center into the enclosures is expected. As this prevents the buildup of pressure, the rotation and shear stiffnesses depend only on the bending stiffness of the enclosures. To analyze this relatively large and non-linear motion, the pseudo-rigid-body (PRB) model in Fig. 4 is constructed, based on a PRB model specific for initially curved beams [24]. Although the basis-PRB model is derived for beams clamped on one side, it is fixed at both ends to the horizontal segments in Fig. 4 by using the symmetry of the enclosure and applying the prescribed angle ζ1 from Ref. [24] at both ends. With the initial enclosure angle of θ = π(rad), the parameters for the PRB model in Fig. 4 follow directly from Ref. [24] as summarized in Table 1. Note that Le in Table 1 indicates the initial enclosure length at the middle surface of the wall, opposed to c defined earlier as the enclosure length at the inner wall.

Fig. 4
Model definition of the PRB model used to analyze both the rotation and shear stiffnesses
Fig. 4
Model definition of the PRB model used to analyze both the rotation and shear stiffnesses
Close modal
Table 1

Parameters for the PRB model in Fig. 4 

ParameterSymbolValueUnit
Enclosure lengthLeπ(h0 + t0)/2(m)
Moment of inertiaIbt03/12(m4)
Length 1L10.1296 · Le(m)
Length 2L20.3513 · Le(m)
Stiffness 1K13.3546 · EI/Le(Nm)
Stiffness 2K22.4764 · EI/Le(Nm)
Fixed angle 1ζ10.2008(rad)
Initial angle 2ζ20.7854(rad)
Initial angle 3ζ31.1692(rad)
ParameterSymbolValueUnit
Enclosure lengthLeπ(h0 + t0)/2(m)
Moment of inertiaIbt03/12(m4)
Length 1L10.1296 · Le(m)
Length 2L20.3513 · Le(m)
Stiffness 1K13.3546 · EI/Le(Nm)
Stiffness 2K22.4764 · EI/Le(Nm)
Fixed angle 1ζ10.2008(rad)
Initial angle 2ζ20.7854(rad)
Initial angle 3ζ31.1692(rad)
To derive the rotation and shear stiffnesses, the principle of virtual work is applied to the motion of the upper rigid connector and torsion springs between the links. Due to a similarity in derivation, only equations for the right enclosure are shown here. First, the rotation of each spring θi in terms of the position [x3, y3] and orientation [α] of the upper rigid connector is required. To obtain these rotations, the coordinates [x1RD, y1RD] and [x1RU, y1RU] of the lower and upper springs with stiffness K1 in Fig. 4 are defined as in Eqs. (14)(17), where L0 is the length of the horizontal segments and L1 and ζ1 respectively indicate the length and initial angle of the links closest to the rigid connectors.
x1RD=12L0+L1cos(ζ1)
(14)
y1RD=L1sin(ζ1)
(15)
x1RU=x3+12L0cos(α)+L1cos(α+ζ1)
(16)
y1RU=y312L0sin(α)L1sin(α+ζ1)
(17)
With above coordinates, the orientations of the middle links of length L2 can be defined using the triangle in the bottom right of Fig. 4. The distance ZR between the top and bottom springs is given by Eq. (18). Using this distance ZR, the triangle orientation βR and outer angle γR follow from Eqs. (19) and (20).
ZR=(x1RDx1RU)2+(y1RDy1RU)2
(18)
βR=sin1(x1RDx1RUZR)
(19)
γR=sin1(ZR2L2)
(20)
The absolute orientation of the middle links (L2) can now be described both in terms of the triangle (βR, γR) and the angles between the links (ζi, θi, α). Combining both options results in the rotation of the springs described by Eqs. (21)and (23).
θ1RD=γR+βRζ1ζ2
(21)
θ1RU=γRβRζ1ζ2α
(22)
θ2R=θ1RDθ1RUα
(23)
Next, static equilibrium is determined by setting the total virtual work Wt for a virtual rotation δα to zero as in Eq. (24), where i is the summation over all springs and δθi the virtual spring rotation. To evaluate torque T, all δθi need to be expressed in terms of δα. In this model the CoR is assumed in the fluid center, as indicated in Fig. 4 and later verified in the FEM model. Thus, a virtual rotation δα also causes virtual displacements δx3 and δy3, leading to Eq. (25) for all δθi.
Wt,rotation=Tδαi=16Kiθiδθi=0
(24)
δθi=[θix3(h0+t02cosα)+θiy3(h0+t02sinα)+θiα]δα
(25)
The secant rotation stiffness Kα is now given by Eq. (26).
Kα=Tα=(i=16Kiθiδθi)/δαα
(26)
For shear stiffness the same derivation is used, where a horizontal force Fx is assumed to translate the upper rigid connector in horizontal direction. This results in Eq. (27) for the virtual work Wt,shear and Eq. (28) for the virtual rotations δθi.
Wt,shear=Fxδx3i=16Kiθiδθi=0
(27)
δθi=θix3δx3
(28)
The secant shear stiffness Ks is now given by Eq. (29).
Ks=Fxx3=(i=16Kiθiδθi)/δx3x3
(29)

3 Finite Element Model and Validation

To verify the analytical stiffness models and investigate further performance behavior, the pressure balanced joint is modeled in comsol multiphysics v5.5. The FEM model consists of a 2D geometry in which the horizontal segments are modeled as rigid domains, following the assumption in the analytical models. The bottom segment is given a fixed constraint while a prescribed displacement/rotation is imposed on the upper segment. The enclosures are modeled as linear elastic material with E = 1 (MPa) and ν = 0.49 (−) to simulate an incompressible silicon rubber, also used as reference material throughout the rest of this study. To implement the internal pressure buildup in the incompressible fluid, a pressure p is added as boundary load on the entire inner boundary of the cell. The magnitude of p is defined in a global equation as in Eq. (30), stating p is increased such that the inner area Ainner equals the initial area A0. A constant 105, the order of number of elements, is added in Eq. (30) to increase the emphasis of this constraint in the overall system of equations and to ensure the constraint to be met when the model converges. The inner area Ainner is evaluated using the divergence theorem as described in Eq. (31) [22], where x and y represent the position and nx and ny the components of the inwards pointing normal vector n at the inner boundary.
(AinnerA01)105=0
(30)
Ainner=Anxx+nyy2dA
(31)
After imposing a motion, the solid mechanics module determines the resulting deformations in the enclosure using non-linear geometry, a physics controlled “extremely fine” mesh, a relative tolerance of 1 · 10−6 and the fully coupled mumps solver. The resulting stiffness is obtained by evaluating the reaction force from the upper rigid domain.

For pure compression, the upper rigid domain is prescribed a downward translation u0y = −dh, with the horizontal translation (u0x) and rotation (ϕ0) fixated as no movement is expected in these degrees-of-freedom (DOF). For rotation however, the CoR is unknown. Therefore both translations are not prescribed, allowing the model to translate the upper rigid domain to a state of equilibrium. Similarly for shear, the vertical translation is not prescribed although rotation is again constrained.

3.1 Validation Method.

The stiffness from the analytical models is compared to the FEM model for a range of design parameters. Two parameter ratios have been identified that influence the design: (L0/h0) and (t0/h0). To independently vary these ratios without scaling the design, a normalized height h0 = 1 (m) is set. Although (L0/h0) influences the stiffness, as confirmed in Sec. 4.2, it does not alter the enclosure geometry which affects the validity of the analytical models. Therefore, only (t0/h0) ratios are considered here.

To analyze joint configurations able to reach state-of-the-art ranges of motion [14], a maximum ratio of (L0/h0) = 4 is allowed. This is based on the maximum rotation αmax, indicated in Fig. 4 and given by Eq. (32), before internal parts make contact assuming the CoR is in the fluid center.
αmax=2tan1(h0/L0)
(32)
Therefore, a range of [2; 4] for (L0/h0) is analyzed in this work, of which (L0/h0) = 3 is used for validation. The (t0/h0) ratio is chosen as [0.01; 0.10] to account for thicknesses from a general plate ratio up to one order thicker enclosures. For the axial stiffness, the compression is approached in n = 25 steps with a maximum compression ratio Cratio=5%, defined in Eq. (33). This Cratio is chosen such that the strain of the enclosure remains below 10% as a linear elastic material model is used [25]. Due to geometrical limitations of the PRB model, the rotation and shear stiffnesses are evaluated from the neutral position up to 5deg rotation and a translation dx with shear ratio Sratio=20% as in Eq. (34).
Cratio=dhh0=0.05
(33)
Sratio=dxh0=0.20
(34)
Figure 5 presents the comparison between the analytical and FEM model. Different combinations of (t0/h0) ratios are plotted to optimally present the encountered effects.
Fig. 5
Comparison between the analytical models (A) and the comsol model (C) for the (a) axial, (b) rotation, and (c) shear stiffnesses, using E = 1 (MPa), v = 0.49 (−), and h0 = 1 (m). The error markers indicate the error of the analytical model relative to comsol.
Fig. 5
Comparison between the analytical models (A) and the comsol model (C) for the (a) axial, (b) rotation, and (c) shear stiffnesses, using E = 1 (MPa), v = 0.49 (−), and h0 = 1 (m). The error markers indicate the error of the analytical model relative to comsol.
Close modal

4 Joint Behavior and Implementation

Both analytical and FEM models are used in this section to analyze the CoR and characteristic stiffness behavior and to introduce an optimization strategy. Additionally, a conversion model and some practical considerations are presented for 3D design implementations.

4.1 Center of Rotation.

Several methods to determine the CoR exist in literature [8,26]. Constructing the instantaneous center is considered most accurate and implemented as described in Ref. [26]. First, the displacements of the endpoints of the upper rigid connector are tracked in the FEM model during rotation. For each displacement, a perpendicular bisector is constructed. The intersection of the two bisectors then indicates the CoR.

From simulations with varying parameters, only the (L0/h0) ratio showed to significantly influence the axis drift. This is expected as the position of the upper rigid connector depends on the constant volume constraint and only the (L0/h0) ratio influences the volume shape and magnitude. For all parameters, the instantaneous CoR is initially located in the fluid center, as indicated in Fig. 4. When normalized to the length L0, the drift of the CoR during rotation for various (L0/h0) ratios is relatively constant as shown in Fig. 6(a), showing a predominant translation to the right. To interpret the effect of this axis drift, the normalized displacement of the upper rigid connector is shown in Fig. 6(b), evaluated at its midpoint (x3, y3). As the CoR appears to translate faster to the right, the rigid connector moves slightly upwards.

Fig. 6
Normalized displacement of (a) the CoR and (b) the upper midpoint (x3, y3) during a 5deg clockwise rotation, simulated for (t0/h0) = 0.05, h0 = 1 (m), E = 1 (MPa), and v = 0.49 (−)
Fig. 6
Normalized displacement of (a) the CoR and (b) the upper midpoint (x3, y3) during a 5deg clockwise rotation, simulated for (t0/h0) = 0.05, h0 = 1 (m), E = 1 (MPa), and v = 0.49 (−)
Close modal

4.2 Parametric Analysis and Optimization.

Using the analytical models, the influence of the design parameters on stiffness behavior is analyzed for h0 = 1 (m). Only the axial and rotation stiffnesses are considered here as the shear stiffness is inherently low. However, if shear deformations cannot be allowed in an intended application, measures are needed as will be suggested in Sec. 5. The analysis uses a compression of Cratio=1%, well within the linear range of the deformation curve of Fig. 5(a). For rotation stiffness, the analytical model was shown to be most accurate at 3.5deg rotation for (L0/h0) = 3. However, the vertical compression of the circular enclosure scales linearly with L0 for a constant rotation. Therefore, the optimum rotation αopt to determine rotation stiffness with a minimal error at other (L0/h0) ratios is corrected by Eq. (35) in this analysis.
αopt=3.5deg3(L0/h0)
(35)
The result of the stiffness analysis is shown in Fig. 7. The effect of the (t0/h0) ratio on both stiffnesses corresponds with the assumption in Sec. 2.1 that thickness t0 scales the axial and rotation stiffnesses with t01 and t03 respectively. Looking at the (L0/h0) ratio, both stiffnesses scale quadratic which is explained as follows. For a given compression dh, the fluid is pushed sideways into the enclosures, and with it the pressure that is generated, increases with L01. As this pressure acts on the upper surface also increased by L01, the axial stiffness scales with (L0/h0)2. For a given rotation, the vertical displacements of the circular enclosures, and with their reaction forces, increase with L01. With the moment-arms of these reaction forces also related to L01, the rotation stiffness scales with (L0/h0)2 as well.
Fig. 7
Parametric analysis on the (a) axial and (b) rotation stiffnesses for h0 = 1 (m), E = 1 (MPa) and v = 0.49 (−)
Fig. 7
Parametric analysis on the (a) axial and (b) rotation stiffnesses for h0 = 1 (m), E = 1 (MPa) and v = 0.49 (−)
Close modal
Additionally, the stiffness ratio Kratio given in Eq. (36) is analyzed as suggested by Ref. [4]. This non-dimensionless ratio indicates the joints performance as an axial stiff rotational joint because the unconstrained rotational DOF should have a relatively low stiffness compared to the constrained axial DOF [8]. The result of this analysis is shown in Fig. 8(a). The (L0/h0) ratio has no influence on Kratio as both axial and rotation stiffnesses scale with (L0/h0)2. Although decreasing thickness t0 significantly improves the stiffness ratio, it must be noted that it decreases the axial stiffness as well.
Kratio=KyKα
(36)
To decide on the trade-off in thickness t0, an optimization is proposed with an example shown in Fig. 8(b). The objective is to maximize Kratio for optimum performance as axial stiff rotational joint. A minimum axial stiffness (e.g., Ky > 5 · 105 (N/m)) should be defined as constraint for the application. Additionally, Eq. (32) can be used to constrain the maximum (L0/h0) ratio for a desired range of motion (e.g., αmax>30deg). As the objective increases for smaller (t0/h0) ratios, the optimum parameters are then located at the intersection of both constraints, as indicated in Fig. 8(b).
Fig. 8
(a) Analysis of the stiffness ratio for h0 = 1 (m), E = 1 (MPa), v = 0.49 (−) and (b) an example of an optimization strategy
Fig. 8
(a) Analysis of the stiffness ratio for h0 = 1 (m), E = 1 (MPa), v = 0.49 (−) and (b) an example of an optimization strategy
Close modal

4.3 2D to 3D Design Conversion.

So far a 2D embodiment has been analyzed for the characteristics of closed form pressure balancing, but extruding the 2D geometry does not result in a feasible design as fluid will leak from the front and back. Potential 3D embodiments are proposed in Fig. 9, including a rectangular design for a single-DOF hinge and a circular design for a double-DOF universal joint, using an axisymmetric topology. The double curvature of the enclosures in Fig. 9 indicates that a significant shearing motion is geometrically limited. Therefore, only designs with rotational DOF are considered useful.

Fig. 9
3D embodiments of the pressure balanced joint as (a) a single-DOF rotational joint and (b) double-DOF universal joint
Fig. 9
3D embodiments of the pressure balanced joint as (a) a single-DOF rotational joint and (b) double-DOF universal joint
Close modal
This step to 3D designs does influence the stiffness properties. The rotation stiffness increases due to the double curvature of the enclosure and stiffness from the front and back segments. As the double curvature depends on the 3D shape, a general correction factor for this increase is not considered possible. The axial stiffness on the other hand decreases as the fluid also moves into the added enclosure segments during compression, reducing the overall required inflation. This effect can be generally predicted as the enclosure inflation in 3D is similar to 2D, apart from a circumferential component that will be discussed in Sec. 5. The prediction uses the increase in enclosure length from the 2D to 3D embodiment, defined as the circumference of the rigid connector adjacent to an enclosure as indicated in Fig. 9. With all enclosure segments in parallel, the axial stiffness is expected to scale inversely proportional to this enclosure length. As will be shown with an example for both the rectangular and circular designs, two steps define the stiffness prediction:
  1. The out-of-plane width b in the 2D analytical model is scaled to equalize the top surface areas of the 2D and 3D design for an equal volume displacement during compression. For the examples, this results in Eqs. (37) and (38), where w is the width of the rectangular joint as in Fig. 9.
    brectangular=w
    (37)
    bcircular=πL0/4
    (38)
  2. A correction factor fcor is defined to predict the 3D axial stiffness from the 2D result as in Eq. (39). The factor is based on the ratio of the 2D (L2D) and 3D (L3D) enclosure lengths, where L2D only consists of two segments at the sides. For the examples, this results in Eqs. (40) and (41).
    Ky,3D=fcorKy,2D
    (39)
    fcor,rectangular=L2DL3D=2brectangular2(w+L0)
    (40)
    fcor,circular=L2DL3D=2bcircularπL0=12
    (41)
To validate the prediction, both designs were implemented in a 3D comsol model using the approach described in Sec. 3. Different however is the use of a 3D geometry, to model both radial and circumferential expansions, and the 3D volume integral from Eq. (42), to evaluate the inner volume Vinner for incompressibility. Simultaneously, the axial stiffnesses were predicted with the analytical steps presented above. The parameters used are h0 = 1 (m), (L0/h0) = 3, (t0/h0) = 0.05, and width w = 6 (m). The results, evaluated at a compression ratio Cratio=1%, are given in Table 2.
Vinner=Anxx+nyy+nzz3dA
(42)
Table 2

Comparison between the predicted and simulated 3D axial stiffness of a rectangular and circular 3D designs

Joint typePredicted KySimulated Ky
Rectangular2.73 · 106 (N/m)2.74 · 106 (N/m)
Circular8.02 · 105 (N/m)8.61 · 105 (N/m)
Joint typePredicted KySimulated Ky
Rectangular2.73 · 106 (N/m)2.74 · 106 (N/m)
Circular8.02 · 105 (N/m)8.61 · 105 (N/m)

4.4 Combined Loadcases.

The axial and rotation stiffnesses were separately analyzed in preceding sections for the characteristics of closed form pressure balancing. However, a more realistic load case will be one with an axial load during rotation, causing an internal pressure that pushes the enclosures into a circular shape. As this negatively influences their bending stiffness, an increase in rotation stiffness is expected for increasing axial loads. It is therefore important to identify its effect on performance.

Two simulations are compared for a rotation of the 3D circular design with its parameters from Sec. 4.3. First, the upper rigid domain is rotated over 10deg while it is free to translate vertically. Second, a constant vertical load Fload on the upper domain is simulated during rotation by prescribing a variable translation dz to it. This dz is defined in a global equation stating that its value is increased such that the vertical reaction force Fz equals Fload as in Eq. (43). For two axial loads, corresponding to a compression ratio of 1% and 2%, the resulting actuation torque is shown in Fig. 10.
(FzFload1)=0
(43)
Fig. 10
Effect of a combined loadcase with an axial load on the rotation stiffness in a 3D circular design, simulated for (L0/h0) = 3, (t0/h0) = 0.05, h0 = 1 (m), E = 1 (MPa), and v = 0.49 (−)
Fig. 10
Effect of a combined loadcase with an axial load on the rotation stiffness in a 3D circular design, simulated for (L0/h0) = 3, (t0/h0) = 0.05, h0 = 1 (m), E = 1 (MPa), and v = 0.49 (−)
Close modal

5 Discussion

5.1 Design Model.

The analytical design model presented in Sec. 2, able to estimate the onset axial and rotation stiffnesses within 10% error, provides a design tool to obtain initial design parameters for further optimization. However, three limitations exist. Application of the model is restricted to relative small compressions, because of the use of a linear elastic material, and to limited rotations, as the two middle links in the PRB model align at some point. Furthermore, the constraint of a constant inner volume cannot be implemented in the PRB model and only embodiments with specific angular segments θ can be analyzed for which parameters of the PRB model exist [24]. These last two are limitations as they prevent an analytical analysis of combined loadcases during rotation.

5.2 Finite Element Method Model and Validation.

Section 3 introduced a FEM modeling approach and validation of the analytical models, requiring two remarks. First, the use of Young’s modulus E = 1 (MPa) is not considered to influence the validation or performances presented, apart from stiffness magnitudes, as all stiffnesses scale linearly with E. Second, although no experimental validation has been presented in this paper, it is noted that a prototype has been tested in a case study in Ref. [27]. Also, a similar FEM modeling approach was experimentally validated in Ref. [22].

The remainder of this section discusses the validation results from Fig. 5. For the axial stiffness, a structural underestimation of the simulated stiffness is seen. This can be explained by neglecting the bending stiffness in both the enclosures and their connection points to the rigid connectors in the analytical model. At relative small compressions, this effect is more significant for thicker enclosures, as expected. However, at larger compressions the analytical model performs better for thicker enclosures. This is probably caused by a growing error in the adapted enclosure thickness (step 4 in Fig. 3), which is more significant for thin enclosures.

The analytical rotation stiffness shows an initial underestimation before rapidly increasing, which is explained by the PRB model. Links L1 next to the horizontal segments are fixed at an angle ζ1, so the required deformation is distributed over all torsion springs. In the FEM model however, the enclosures maintain more circular for a constant inner volume, resulting in larger enclosure deformations close to the rigid connectors and increased stiffness. At larger joint rotations, the PRB model approaches the singularity where the middle links align. As this requires large rotations in the torsion springs, the rotation stiffness rapidly increases.

The analytical shear stiffness shows a structural underestimation as well. Similar to the rotation stiffness, the constraint on incompressibility causes larger deformations close to the rigid connectors in the FEM model, resulting in a larger stiffness. Despite that, it is noted that the shear stiffness is significantly lower than the axial stiffness. This can be explained by the rolling motion of the enclosures during a horizontal movement of the upper rigid connector, best compared to the motion seen in tape loops [28].

5.3 Joint Performance.

The assumption on the central position of the CoR was verified in Sec. 4.1, although Fig. 6 showed that the CoR moves sideways during rotation. This axis drift, relative to the joints dimensions, is significant compared to most conventional compliant joints [29]. Additionally, the axis drift is expected to increase significantly for combined load cases as the compressive load creates a horizontal force component during rotation. Due to the relative low shear stiffness, this horizontal force can translate the upper part of the joint sideways and obstruct the joint’s functionality as a hinge. Adding stiffness elements to improve the shear stiffness is therefore recommended, as demonstrated in the case study in Ref. [27]. Here, an intersecting set of wires was added inside the cell between the top and bottom rigid connectors that created shear stiffness without affecting the rotation stiffness.

The performance in terms of stiffnesses was presented in Sec. 4.2. For this analysis, only secant stiffnesses have been used as Fig. 5 showed all stiffnesses to behave approximately linear. With respect to the stiffness ratio and optimization, two remarks are made. First, the stiffness ratio is not dimensionless and improves for smaller designs. Second however, stresses and failure have not been considered. When included, they will limit the minimum size and achievable stiffness ratio as stresses increase for smaller designs. Also with respect to the range of motion, it is recommended to analyze material behavior in future research.

5.4 Conversion to 3D Designs.

Section 4.3 described the negative effect of conversion to 3D designs on the axial and rotation stiffnesses and showed in Table 2 that the 3D axial stiffness can be estimated from the 2D model. The underestimation for the circular embodiment in Table 2 is the result of this particular 3D shape, that apart from the radial expansion as in the 2D design also expands in circumferential direction, which is not included in the prediction method. In the rectangular shape, this is less relevant as most segments only expand radially. Although no 3D prediction is presented for the rotation stiffness, it is noted that this stiffness is not required in the optimization for finding initial design parameters, as described in Sec. 4.2.

Another aspect for 3D implementations however, is the effect of combined load cases as presented in Sec. 4.4. Figure 10 confirms that compressive loads during a rotation significantly increase rotation stiffness. An additional explanation for this effect is a difference in deformation of the enclosure at opposite sides of the joint. On one side the rigid connector tips move apart resulting in an inward movement of the enclosure, while on the other side the tips move closer together resulting in an outward movement of the enclosure. With the pressure in both enclosures being equal, the shape difference causes a reaction torque on the upper rigid domain, increasing rotation stiffness. As this effect is inherent to the enclosures, it should be taken into account when designing for a maximum rotation stiffness.

6 Conclusion

This work presents the characteristics of closed form pressure balancing, a design approach to improve axial performance in compliant rotational joints. By using an incompressible fluid, the axial and rotational stiffnesses become decoupled which prevents an axial stiffness decrease during rotation. Additionally, compressive forces are converted to tensile loads which prevent buckling and thus increase loading capability. Analytical and FEM models developed are capable to determine initial stiffnesses of a 2D embodiment and provide a design tool to find optimum design parameters. With the proposed design, relative large ranges of motion can be achieved that only depend on joint length. Optimal joint performance in terms of stiffnesses is achieved at smaller enclosure thicknesses, although a minimum axial stiffness should be maintained. The center of rotation, initially at the center of the joint, shows a horizontal drift up to 8% relative to the joints length during rotation. However, a relative low shear stiffness makes the CoR susceptible to external loads, for which additional elements in parallel are recommended. Additionally, a conversion model is presented capable to predict the axial stiffness in 3D designs from the 2D model. Finally, axial loads are shown to significantly increase the rotation stiffness, which needs to be taken into account in the design process.

Conflict of Interest

There are no conflicts of interest.

Data Availability Statement

The datasets generated and supporting the findings of this article are obtainable from the corresponding author upon reasonable request.

References

1.
Howell
,
L. L.
,
2013
,
21st Century Kinematics
,
Springer
,
London
, pp.
189
216
.
2.
Lusty
,
C.
,
Bailey
,
N. Y.
, and
Keogh
,
P. S.
,
2017
, “
A Comparison of the Capacity for Precise Motion Control Between Hinged Joints and Flexure Joints
,”
Proceedings of the ASME 2017 International Mechanical Engineering Congress and Exposition. Volume 4A: Dynamics, Vibration, and Control
,
Tampa, FL
,
Nov. 3–9
, ASME, p. V04AT05A019.
3.
Kang
,
D.
, and
Gweon
,
D.
,
2013
, “
Analysis and Design of a Cartwheel-Type Flexure Hinge
,”
Precis. Eng.
,
37
(
1
), pp.
33
43
.
4.
Trease
,
B. P.
,
Moon
,
Y. -M.
, and
Kota
,
S.
,
2004
, “
Design of Large-Displacement Compliant Joints
,”
ASME J. Mech. Des.
,
127
(
4
), pp.
788
798
.
5.
Jensen
,
B. D.
, and
Howell
,
L. L.
,
2002
, “
The Modeling of Cross-Axis Flexural Pivots
,”
Mech. Mach. Theory
,
37
(
5
), pp.
461
476
.
6.
Martin
,
J.
, and
Robert
,
M.
,
2011
, “
Novel Flexible Pivot With Large Angular Range and Small Center Shift to be Integrated Into a Bio-Inspired Robotic Hand
,”
J. Intel. Mater. Syst. Struct.
,
22
(
13
), pp.
1431
1437
.
7.
Henein
,
S.
,
Spanoudakis
,
P.
,
Droz
,
S.
,
Myklebust
,
L. I.
, and
Onillon
,
E.
,
2003
, “
Flexure Pivot for Aerospace Mechanisms
,”
10th European Space Mechanisms and Tribology Symposium
,
San Sebastian, Spain
,
Sept. 24–26
, pp.
285
288
.
8.
Xu
,
P.
,
Jingjun
,
Y.
,
Guanghua
,
Z.
,
Shusheng
,
B.
, and
Zhiwei
,
Y.
,
2008
, “
Analysis of Rotational Precision for an Isosceles-Trapezoidal Flexural Pivot
,”
ASME J. Mech. Des.
,
130
(
5
), p.
052302
.
9.
Guérinot
,
A. E.
,
Magleby
,
S. P.
,
Howell
,
L. L.
, and
Todd
,
R. H.
,
2004
, “
Compliant Joint Design Principles for High Compressive Load Situations
,”
ASME J. Mech. Des.
,
127
(
4
), pp.
774
781
.
10.
Wiersma
,
D. H.
,
Boer
,
S. E.
,
Aarts
,
R. G. K. M.
, and
Brouwer
,
D. M.
,
2014
, “
Design and Performance Optimization of Large Stroke Spatial Flexures
,”
ASME J. Comput. Nonlinear. Dyn.
,
9
(
1
), p.
011016
.
11.
Cannon
,
J. R.
,
Lusk
,
C. P.
, and
Howell
,
L. L.
,
2005
, “
Compliant Rolling-Contact Element Mechanisms
,”
Proceedings of the ASME 2005 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference. Volume 7: 29th Mechanisms and Robotics Conference, Parts A and B
,
Long Beach, CA
,
Sept. 24–28
,
American Society of Mechanical Engineers Digital Collection
, pp.
3
13
.
12.
Jeanneau
,
A.
,
Herder
,
J.
,
Laliberté
,
T.
, and
Gosselin
,
C.
,
2004
, “
A Compliant Rolling Contact Joint and Its Application in a 3-DOF Planar Parallel Mechanism With Kinematic Analysis
,”
Proceedings of the ASME 2004 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference. Volume 2: 28th Biennial Mechanisms and Robotics Conference, Parts A and B
,
Salt Lake City, UT
,
Sept. 28–Oct. 2
,
American Society of Mechanical Engineers Digital Collection
, pp.
689
698
.
13.
Brouwer
,
D. M.
,
Meijaard
,
J. P.
, and
Jonker
,
J. B.
,
2009
, “
Elastic Element Showing Low Stiffness Loss at Large Deflection
,”
Proceedings of the 24th Annual Meeting of the American Society of Precision Engineering
,
Monterey, CA
, pp.
30
33
.
14.
Naves
,
M.
,
Aarts
,
R. G. K. M.
, and
Brouwer
,
D. M.
,
2019
, “
Large Stroke High Off-Axis Stiffness Three Degree of Freedom Spherical Flexure Joint
,”
Precis. Eng.
,
56
, pp.
422
431
.
15.
Folkersma
,
K. G. P.
,
Boer
,
S. E.
,
Brouwer
,
D. M.
,
Herder
,
J. L.
, and
Soemers
,
H. M. J. R.
,
2012
, “
A 2-DOF Large Stroke Flexure Based Positioning Mechanism
,”
Proceedings of the ASME 2012 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference. Volume 4: 36th Mechanisms and Robotics Conference, Parts A and B
,
Chicago, IL
,
Aug. 12–15
,
American Society of Mechanical Engineers Digital Collection
, pp.
221
228
.
16.
Nguyen
,
T. T.
,
Dao
,
T. -P.
, and
Huang
,
S. -C.
,
2017
, “
Biomechanical Design of a Novel Six DOF Compliant Prosthetic Ankle-Foot 2.0 for Rehabilitation of Amputee
,”
Proceedings of the ASME 2017 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference. Volume 5A: 41st Mechanisms and Robotics Conference
,
Cleveland, OH
,
Aug. 6–9
,
American Society of Mechanical Engineers Digital Collection
, p. V05AT08A013.
17.
Nguyen
,
T. T.
,
Le
,
H. G.
,
Dao
,
T. -P.
, and
Huang
,
S. -C.
,
2017
, “
Evaluation of Structural Behaviour of a Novel Compliant Prosthetic Ankle-Foot
,”
2017 International Conference on Mechanical, System and Control Engineering (ICMSC)
,
St. Petersburg, Russia
, pp.
58
62
.
18.
Halverson
,
P. A.
,
Howell
,
L. L.
, and
Bowden
,
A. E.
,
2009
, “
A Flexure-Based Bi-Axial Contact-Aided Compliant Mechanism for Spinal Arthroplasty
,”
Proceedings of the ASME 2008 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference. Volume 2: 32nd Mechanisms and Robotics Conference, Parts A and B
,
Brooklyn, NY
,
Aug. 3–6
,
American Society of Mechanical Engineers Digital Collection
, pp.
405
416
.
19.
Nijssen
,
J.
,
Kempenaar
,
A.
, and
Diepeveen
,
N.
,
2018
, “
Development of an Interface Between a Plunger and an Eccentric Running Track for a Low-Speed Seawater Pump
,”
Fluid Power Networks: Proceedings 11th International Fluid Power Conference
,
Aachen, Germany
,
Mar. 19–21
, pp.
370
379
.
20.
Sneed
,
R. C.
,
Cash
,
M. F.
,
Chambers
,
T. S.
, and
Janzen
,
P. C.
,
2010
, “
Six Degrees of Freedom, Sub-Micrometer Positioning System for Secondary Mirrors
,”
Proc. SPIE 7733, Ground-Based and Airborne Telescopes III
, Vol. 77332R, pp.
1002
1012
.
21.
Cash
,
M.
,
Bruch
,
D.
,
Jahn
,
B.
, and
Keas
,
P.
,
2008
, “
Practical Considerations of Joint Friction and Backlash in Large Ground-Based Telescope Secondary Optic Positioning Systems
,”
Proc. SPIE 7018, Advanced Optical and Mechanical Technologies in Telescopes and Instrumentation
, Vol. 70181T, pp.
652
663
.
22.
Gomes
,
R. F. P.
,
Nijssen
,
J. P. A.
, and
van Ostayen
,
R. A. J.
,
2020
, “
Design of a Compliant Hinge Based on Closed Form Pressure Balancing
,”
Proceedings of the ASME 2020 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference. Volume 10: 44th Mechanisms and Robotics Conference (MR)
,
Virtual, Online
,
Aug. 17–19
,
American Society of Mechanical Engineers Digital Collection
, p. V010T10A003.
23.
Bower
,
A. F.
,
2009
,
Applied Mechanics of Solids
,
CRC Press
,
Boca Raton, FL
.
24.
Venkiteswaran
,
V. K.
, and
Su
,
H. -J.
,
2018
, “
A Versatile 3R Pseudo-Rigid-Body Model for Initially Curved and Straight Compliant Beams of Uniform Cross Section
,”
ASME J. Mech. Des.
,
140
(
9
), p.
092305
.
25.
Martins
,
P. a. L. S.
,
Jorge
,
R. M. N.
, and
Ferreira
,
A. J. M.
,
2006
, “
A Comparative Study of Several Material Models for Prediction of Hyperelastic Properties: Application to Silicone-Rubber and Soft Tissues
,”
Strain
,
42
(
3
), pp.
135
147
.
26.
Linß
,
S.
,
Erbe
,
T.
,
Theska
,
R.
, and
Zentner
,
L.
,
2011
, “
The Influence of Asymmetric Flexure Hinges on the Axis of Rotation
,”
Proceedings of the 56th International Scientific Colloquium
,
Ilmenau University of Technology, Germany
,
Sept. 12–16
.
27.
Nijssen
,
J. P. A.
,
2022
, “Compliant Full Film Lubricated Bearings: Concept Design & Development,” Ph.D. thesis,
Delft University of Technology
,
Delft, Netherlands
.
28.
Vehar
,
C.
,
Kota
,
S.
, and
Dennis
,
R.
,
2008
, “
Closed-Loop Tape Springs as Fully Compliant Mechanisms: Preliminary Investigations
,”
Proceedings of the ASME 2004 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference. Volume 2: 28th Biennial Mechanisms and Robotics Conference, Parts A and B
,
Salt Lake City, UT
,
Sept. 28–Oct. 2
,
American Society of Mechanical Engineers Digital Collection
, pp.
1023
1032
.
29.
Farhadi Machekposhti
,
D.
,
Tolou
,
N.
, and
Herder
,
J. L.
,
2015
, “
A Review on Compliant Joints and Rigid-Body Constant Velocity Universal Joints Toward the Design of Compliant Homokinetic Couplings
,”
ASME J. Mech. Des.
,
137
(
3
), p.
032301
.