A Closed-Form Kinematic Model for Fiber-Reinforced Elastomeric Enclosures

Felt, Wyatt; David Remy, C.

doi:10.1115/1.4038220

In this work, we present a closed-form model, which describes the kinematics of fiber-reinforced elastomeric enclosures (FREEs). A FREE actuator consists of a thin elastomeric tube surrounded by reinforcing helical fibers. Previous models for the motion of FREEs have relied on the successive compositions of “instantaneous” kinematics or complex elastomer models. The model presented in this work classifies each FREE by the ratio of the length of its fibers. This ratio defines the behavior of the FREE regardless of the other parameters. With this ratio defined, the kinematic state of the FREE can then be completely described by one of the fiber angles. The simple, analytic nature of the model presented in this work facilitates the understanding and design of FREE actuators. We demonstrate the application of this model in an actuator design case study.

Introduction

Soft, fluid-driven actuators use structured compliance to create motion from the expansion of flexible volumes. Actuators based on these principles may, for example, contract along a central axis like biological muscles [1]. The McKibben muscle is a popular variety of such “pneumatic artificial muscles.” Each McKibben muscle consists of an elastomeric tube surrounded by a sleeve of braided helical fibers. The braid is made up of equal numbers of right-handed and left-handed helices of the same pitch. They were developed in the 1950 s by their namesake, Joseph Laws McKibben [2]. Since that time, they have been used in numerous applications such as legged robots and human assistive devices [3].

Bishop-Moser [4] extended the functional principle of McKibben muscles to a generalized class of cylindrical soft actuators known as fiber-reinforced elastomeric enclosures (FREEs). Like McKibben muscles, FREE actuators are formed from two sets of identical helical fibers. A “set” or “family” of fibers is a group of fibers characterized by the same angle with respect to the cylinder axis (e.g., α or β). In the McKibben muscle, the two sets maintain equal and opposite angles (α = −β). In a FREE, the fibers are wound with fiber angles selected by a designer to result in a desired behavior (Fig. 1). This choice of configuration permits complex actuated behaviors such as twisting while extending or contracting.

Fig. 1

View large Download slide

A FREE is a cylindrical soft pneumatic actuator. A FREE consists of an elastomeric tube surrounded by sets (or “families”) of identical helical fibers. The fibers in each family have the same angle with respect to the long axis (e.g., α) and the same unwound length (e.g., b_α). Shown here are examples of two-fiber-family FREEs. The families are, respectively, described by the angles α and β. The unwound length of the helical fibers in each family, b_α and b_β, remains constant, whereas the axial length l and diameter D change during actuation.

The kinematics of FREEs was initially developed only for small deformations. Krishnan et al. [5] described the motion of FREEs with instantaneous strain equations. These equations described the transformation of the FREE as a small change relative to its current configuration. To calculate the evolution of a FREE actuator over large deformations, the instantaneous strains were successively composed. This fiber-only kinematic model allowed the designer to consider the kinematics of the FREE fibers under tension without considering the elastomer.

Combined fiber-elastomer methods model the motion of specific FREE actuators with specified geometries and defined elastomeric properties. For instance, finite element solvers can be used to model FREEs [6]. Recent work has begun to explore the use of constitutive models to predict the motion of unloaded contracting, twisting, or bending fiber-reinforced actuators [7] and torsionally loaded fiber-reinforced actuators [8]. These constitutive methods rely on models of strain energy in the elastomers and fibers to relate internal pressures to predicted deformations for specific actuator geometries. The governing equations of Ref. [5] can also be adapted to consider the elastomer [9].

The model presented in this work allows the designer to consider the kinematics of FREE fibers independent of specific geometry and material choices. For instance, in an finite element analysis model, one would need to define the inner and outer diameter of the elastomer, the elastic material, a model for the elastic behavior, the length of the actuator, the two pitches of the fiber families, the fiber material, the number of reinforcing fibers in each family, the pressure ranges, etc. Once these choices have been made, a computationally intensive process must be conducted to evaluate the behavior of the specified FREE. After conducting a number of these numerical experiments, it is perhaps possible to extract some design heuristics from the observed behavior, but it is difficult to generalize the results.

In the model presented in this work, the behavior of the FREE actuators is described in terms of a single variable. This variable, η, is the ratio the lengths of the fibers in the two sets and is independent of the specific actuator geometry and materials. η is a design choice that is conserved over the course of actuation. The introduction of η allows this work to extend models developed for McKibben muscles [10] into the broader class of FREE actuators. The mathematical model presented in this work describes large-deformation kinematics of FREEs in closed-form and without the need to compose a succession of instantaneous strains.

The kinematic model presented here describes the length, diameter, rotation, and volume of FREE actuators as functions of the fiber angle β. The structure of this model provides a common language to describe every cylindrical, two-helical-fiber-set FREE by parameters that define the behavior (η) and state (β) of the actuator. The size and geometry of the actuator is given by the length b_β of the β fibers and the diameter D₀ (when β = −π/2). The simple, analytic nature of the model presented in this work facilitates the understanding and design of FREE actuators.

Kinematic Modeling

Assumptions and Definitions.

Our model assumes that the fiber-reinforced actuators are made from two sets of helical, inextensible fibers with respective angles α and β (Fig. 1). The fibers surround an elastomeric bladder containing the pressurized fluid. We assume that sufficient fibers are used to prevent the bladder from bulging between the gaps in the fibers. The elastomer is assumed to be infinitely extensible with negligible stiffness. The fibers are assumed to be always under tension from the internal pressure in the bladder. Because the individual fibers within the families are identical, the kinematics of only one fiber in each family need to be considered.

We also assume that the profile of the FREE actuator remains cylindrical. This approximates an unbent actuator away from its ends. At the actuator ends, the diameter tapers to match the fixed-diameter of the end [11]. The cylindrical assumption allows us to adapt simple helical formulas that have long been used to describe McKibben muscles [10]. The length of this cylinder is l and its diameter is D. All fibers wrap around it in a helical fashion. The individual fibers in their respective families behave identically. They all have the same axial length l and diameter D as the cylinder. Under these assumptions, we can relate the length of the cylinder l to the “unwound” length b of the fibers in a family via the cosine of their angle

l = b_{α} cos (α) = b_{β} cos (β)

(1)

Similarly, we can express the diameter of the cylinder D as a function of the fiber angles and fiber lengths

D = \frac{b_{α} sin (α)}{n_{α} π} = \frac{b_{β} sin (β)}{n_{β} π}

(2)

The diameter is additionally a function of the number n of times that the fiber circles the axis. For example, n_α = 0.75 signifies that, at the current diameter, the fibers of the α family circle the axis three-quarters of one time. The sign of n indicates the handedness of the helix (positive for right handed) and matches the sign of the corresponding fiber angle

\begin{matrix} sign (n_{α}) = sign (α) \\ sign (n_{β}) = sign (β) \end{matrix}

(3)

The cylinder diameter and length will change as the actuator volume expands. These changes will correspond to changes in the fiber angles α and β. The number n of fiber turns may also change as the actuator ends rotate relative to one another about the cylinder axis. The unwound length b of the fibers, however, remains constant.

Krishnan et al. [5] analyzed the instantaneous kinematics of FREEs with values of α and β between −90 deg and 90 deg. As they noted, however, the symmetry of FREEs makes this formulation redundant. We eliminate this redundancy, without loss of generality, by deliberately specifying which fiber family is labeled by α and which is labeled by β. We use α to describe the family with the greater or equal unwound fiber length b (i.e., b_α ≥ b_β, and thus

| α | \geq | β |

⁠). Furthermore, we can restrict our analysis to helices formed by the α fibers that are right handed and thus maintain a positive value of α. Because α is always positive, and α is strictly greater than β (β is negative when

| α | = | β |

⁠), both angles cannot be equal to zero simultaneously

\begin{matrix} 0 < α \leq \frac{π}{2} \\ - α \leq β < α \end{matrix}

(4)

The behavior of FREE actuators with left-handed (negative α) α fibers will simply be mirror-symmetric to the behavior modeled here (e.g., rotation that is modeled as counterclockwise will be clockwise in the actuator with left-handed α fibers).

Behavior Described by η, State by β.

The behavior of a FREE actuator is defined by the ratio η of the unwound lengths b of the fibers from the two sets

\begin{matrix} η = \frac{b_{β}}{b_{α}} \\ 0 < η \leq 1 \end{matrix}

(5)

Under the assumptions of this model, the angle α follows uniquely from the angle β. With η, a simple analytic relationship between the fiber angles can be developed from Eqs. and

α = {cos}^{- 1} (η cos (β))

(6)

Equation (6) provides a clear functional relationship between the angles of each FREE (Fig. 2) parameterized by η.

Fig. 2

The ratio 0<η=(bβ/bα)≤1 of the lengths of fibers in the two families determines the feasible angle combinations for a FREE actuator. The angle β can be used to specify the state of the actuator along the functional relationship of possible combinations. The locked manifold (dashed black line) is the set of angle combinations that maximize the volume of the actuator. When the volume of the actuator is increased by an internal pressure, the configuration of the actuator will advance toward the locked manifold. The gray regions of the figure are redundant and thus not specifically defined by this model. β is constrained to remain strictly less than α. Thus, when η = 1 (α = −β), β is constrained to be less than zero. This is indicated by the circle at the origin. The feasible region includes the lineα = −β but the line α = β is excluded. The values of η shownin this figure (and the others) are as follows: 1, 1516, 78, 1316, 34, 1116, 58, 916, 12, 716, 38, 516, 14, 316, 18, and 116.

View large Download slide

The ratio $0 < η = (b_{β} / b_{α}) \leq 1$ of the lengths of fibers in the two families determines the feasible angle combinations for a FREE actuator. The angle β can be used to specify the state of the actuator along the functional relationship of possible combinations. The locked manifold (dashed black line) is the set of angle combinations that maximize the volume of the actuator. When the volume of the actuator is increased by an internal pressure, the configuration of the actuator will advance toward the locked manifold. The gray regions of the figure are redundant and thus not specifically defined by this model. β is constrained to remain strictly less than α. Thus, when η = 1 (α = −β), β is constrained to be less than zero. This is indicated by the circle at the origin. The feasible region includes the lineα = −β but the line α = β is excluded. The values of η shownin this figure (and the others) are as follows: 1, $\frac{15}{16}, \frac{7}{8}, \frac{13}{16}, \frac{3}{4}, \frac{11}{16}, \frac{5}{8}, \frac{9}{16}, \frac{1}{2}, \frac{7}{16}, \frac{3}{8}, \frac{5}{16}, \frac{1}{4}, \frac{3}{16}, \frac{1}{8}$ ⁠, and $\frac{1}{16}$ ⁠.

The relationship between the fiber angles described by Eq. (6) has been observed previously but never defined in such an explicit form. Krishnan et al. [5] noted that each two-fiber-family FREE belongs to “a one-dimensional family of fiber angle configurations.” With the model presented here, it is now clear that the ratio η and Eq. (6) can be used to define this one-dimensional family of fiber angle configurations.

Size Described by b_β and D₀.

The length, diameter, and volume of a FREE actuator can be described as functions of η and β and variables that describe the dimensions of the particular actuator, b_β and D₀. The first of these, b_β, is the unwound length of the individual fibers in the β family (Figs. 1 and 3). The diameter D₀ is a standardizing measurement used for McKibben muscle actuators [10]. This quantity can also be defined for FREE actuators (Fig. 3). D₀ is calculated from Eq. . It is the diameter of the helices if β = −π/2 (and thus, α = π/2)

D_{0} = \frac{b_{α}}{n_{0, α} π} = - \frac{b_{β}}{n_{0, β} π}

(7)

where n₀ is the number of fiber turns at that diameter.

Fig. 3

When the helices are planar circular arcs, the diameter is D0. Shown here are the paths of an α fiber and a β fiber for η = 0.5, n0,α = 0.5 and n0,β = −0.25.

View large Download slide

When the helices are planar circular arcs, the diameter is D₀. Shown here are the paths of an α fiber and a β fiber for η = 0.5, n_0,_α = 0.5 and n_0,_β = −0.25.

Thus, each FREE actuator can be defined by the actuator-specific quantities η, b_β, and D₀. The state of the actuator is given by β. Examples of FREEs with the same values of b_β and D₀ (but different values of η and β) are shown in Fig. 4.

Fig. 4

Examples of FREEs with the same bβ and D0 but different η levels at various values of β. For the actuators shown here, D0 is half of bβ. Each row corresponds to a value of η and each column a value of β.

View large Download slide

Examples of FREEs with the same b_β and D₀ but different η levels at various values of β. For the actuators shown here, D₀ is half of b_β. Each row corresponds to a value of η and each column a value of β.

Calculating Rotation Δn, Diameter D, Length l, Volume V, and Surface Area A_surf.

The model presented in this work allows the length, rotation, diameter and volume of a FREE actuator to be described as functions of the actuator state β. The expression for the axial length l(β) (Fig. 5) is straightforward and is given by Eq.

l (β) = b_{β} cos (β)

(8)

Fig. 5

View large Download slide

The length and diameter of FREE actuators for various values of η across β. The length of the actuator is the cosine function scaled by b_β and is always greatest at β = 0 (l = b_β). The diameter scales linearly with D₀ and always increases with decreasing values of β. The maximum diameter is D₀.

Note that the length of the actuator scales linearly with the length b_β and that the maximum length achievable by the actuator is b_β.

The axial rotation of the actuator is designated by Δn and is conserved in the fiber turns n_α and n_β at D

\begin{matrix} n_{α} = n_{0, α} + Δ n \\ n_{β} = n_{0, β} + Δ n \end{matrix}

(9)

where Δn is zero when D = D₀. Equation leads to an expression for Δn

Δ n = - n_{0, β} \frac{sin (β) + sin (α)}{sin (α) - η sin (β)}

(10)

Note that there is no rotation when η = 1 (McKibben muscle)

\begin{matrix} η = 1 \to β = - α \to \\ sin (β) + sin (α) = sin (- α) + sin (α) = 0 \to \\ Δ n = 0 \end{matrix}

(11)

Substituting Eq. into Eq. gives the axial rotation Δn(β) in radians

Δ n (β) = - n_{0, β} \frac{\sqrt{1 - η^{2} {cos}^{2} (β)} + sin (β)}{\sqrt{1 - η^{2} {cos}^{2} (β)} - η sin (β)}

(12)

The rotation Δn scales linearly with $n_{0, β} = (- b_{β} / π D_{0})$ (Fig. 6).

Fig. 6

View large Download slide

The axial rotation of FREE actuators for various values of η across β. For η = 1 there is no rotation. The values shown here are normalized by n_0,_β. Note that $n_{0, β} = (- b_{β} / π D_{0})$ ⁠.

The diameter D(β) is found by substituting Eqs. , , , and into Eq. (Fig. 5)

D (β) = D_{0} \frac{\sqrt{1 - η^{2} {cos}^{2} (β)} - η sin (β)}{1 + η}

(13)

The diameter scales linearly with D₀.

Assuming that the thickness of the elastomer inside the fibers is negligible, the internal volume of the actuator is given by the volume V(β) of the cylinder contained within the fiber helices (Fig. 7)

V (β) = \frac{π}{4} D {(β)}^{2} l (β)

(14)

Fig. 7

View large Download slide

The volume of FREE actuators for various values of η across β. The volume values shown here are normalized by $D_{0}^{2} b_{β}$ ⁠. Each value of η corresponds to a unique angle β_LM that maximizes the cylinder volume.

Substituting Eqs. and into Eq. yields

V (β) = \frac{π}{4} D_{0}^{2} b_{β} \frac{cos (β) {(\sqrt{1 - η^{2} {cos}^{2} (β)} - η sin (β))}^{2}}{{(1 + η)}^{2}}

(15)

From Fig. 7, it is clear that the volume of each FREE type has a unique maximal point. For example, for the McKibben muscle case when η = 1 and α = −β, the volume is maximized when β ≈ −54.7 deg. This has long been known [10]. Here, the value of β that maximizes the volume is designated β_LM. Smaller values of η correspond to less negative values of β_LM. Because internal pressures drive the volume to expand, an actuator fabricated with a β value greater than β_LM will decrease in β under actuation. Similarly, actuators with β less than β_LM will increase in β under actuation. The derivative of the volume with respect to β is given by

\begin{matrix} \frac{d V}{d β} = - \frac{π}{4} D_{0}^{2} b_{β} \frac{{(γ - η sin (β))}^{2} (sin (β) γ + 2 η {cos}^{2} (β))}{{(1 + η)}^{2} γ} \\ γ = \sqrt{1 - η^{2} {cos}^{2} (β)} \end{matrix}

(16)

The maximum volume occurs when

d V / d β

is zero. Equation yields the following invertible expression relating η to the angle β_LM that maximizes the volume

\begin{matrix} \frac{d V}{d β} = 0 \to \\ η = \frac{1}{cos (β_{LM}) \sqrt{4 \cot^{2} (β_{LM}) + 1}} \\ - 2 {tan}^{- 1} (\sqrt{2 - \sqrt{3}}) \leq β_{LM} < 0 \end{matrix}

(17)

The α and β combinations that maximize FREE volume were described by Krishnan et al. [5] as a “locked manifold” (Fig. 2)

1 + 2 \cot (α_{LM}) \cot (β_{LM}) = 0

(18)

where the corresponding angles are designated with the subscript LM. “Locked” refers to the fact that internal pressure can no longer drive the actuator to deform because the volume is already maximized. Substituting Eq. (6) into Eq. (18) yields Eq. (17).

Example of Model Application: Actuator Design

Our model enables simple, closed-form design analyses for FREE actuators. As an example, this section applies the proposed model to design an actuator to meet specified kinematics.

Consider a FREE that is specified to contract from an unpressurized length l₁ = 5 cm to l₂ = 4 cm, while rotating a quarter-of-a-turn about its axis. The diameter of the actuator at the contracted state is to be D₂ = 2.5 cm. To ensure that the contracted configuration will be achievable through pressurization, the angle β₂ of the fibers in the contracted state is specified to be 10 deg greater than the angle β_LM that maximizes the volume.

By inspection, for an axially contracting actuator, the following must be true:

- 54.7 \deg < β_{LM} < β_{2} < β_{1} < 0

(19)

The initial unpressurized angle β₁ can be selected numerically such that the quarter-turn rotation constraint is satisfied

n_{1, β} - n_{2, β} = \frac{b_{β} sin (β_{1})}{D_{1} π} - \frac{b_{β} sin (β_{2})}{D_{2} π} = 0.25

(20)

where the rotation comes from the difference in the number of fiber turns given by inverting Eq. . The unknown values in Eq. are functions of the specified kinematics and β₁. b_β comes from inverting the length expression in Eq.

b_{β} = \frac{l_{1}}{cos (β_{1})}

(21)

The angle β₂ of the fibers in the contracted state is given by inverting Eq. and considering that β₂ must be negative due to Eq.

β_{2} = - {cos}^{- 1} (\frac{l_{2}}{b_{β}})

(22)

The angle β_LM is given by the constraint that β₂ be 10 deg greater than β_LM

β_{LM} = β_{2} - 10 \frac{π}{180}

(23)

which, by Eq. , can be used to calculate η. The value of D₀ is found by inverting Eq. with the values of D₂, β₂ and η

D_{0} = D_{2} \frac{1 + η}{\sqrt{1 - η^{2} {cos}^{2} (β_{2})} - η sin (β_{2})}

(24)

which allows D₁ to be calculated with Eq. , η and β₁.

The design specifications are achieved by fabricating an actuator with an unpressurized initial angle β₁ of −8.2 deg (Fig. 8, η = 0.714, b_β = 5.05 cm, D₀ = 3.4 cm).

Fig. 8

View large Download slide

The model can be used to design FREEs that meet specified kinematics. For instance, the FREE shown in this figure was designed to rotate a quarter turn while contracting 20%.

Discussion

The model presented in this work provides a closed-form framework for kinematic analysis and design of FREE actuators. The introduction of η and the analytical relationship between the fiber angles given in Eq. (6) is one of the major contributions of this work. Previously published fiber-only models have relied on “instantaneous” kinematics [5] to incrementally update fiber angles. To solve for large deformations with instantaneous kinematics, the nonlinear equations had to be iteratively solved and composed. The model presented in this work provides analytic functions describing the actuator rotation, length, diameter, and volume. These functions are parameterized by the kinematic state of the actuator given by the angle β. The simplicity of the presented model simplifies the design and understanding of FREEs. In this work, for example, we have shown how the model can be used to design a FREE that achieves desired kinematic behavior.

In addition to the closed-form kinematics, the model presented in this work has several improvements to previous FREE fiber-only models. The deliberate designation of the longer set of fibers with α allows the model presented in this work to describe FREEs with just a diagonal quadrant of the α-β coordinate frame. The ratio η leads to a simple parametric functional relationship between α and β (the first of its kind) allows the state of a FREE with a particular η value to be parameterized by a single angle β (rather than describing the state with potentially infeasible combinations of α and β). Previous work discovered that “every FREE belongs to a one-dimensional family of fiber angle configurations” [5]. This work provides the first analytic description of these configurations.

The model presented in this work shares the assumptions of previously published models fiber models [4,5]. The reformulation presented here is a simplification of the kinematic description in the previous models. Accordingly, the experimental verification of the previous models can be considered verification of the present work.

Like the models that have preceded it, the present model has limitations. External loading, for example, could buckle one or both of the fiber families. This would violate the assumption that the fibers are under tension. The model presented in this work is limited to rotation and/or length changes. Additional fiber families on FREEs can create planar [7] or helical [4,12] bends. The angles of the fibers in this model are constrained to be nonzero. So-called “straight-fiber” actuators are not governed by the equations presented here [13]. The governing equations of Ref. [5] have been adapted to include elastomer effects and noncylindrical deformation [9].

The model presented in this work does not take into account noncylindrical deformations or strain in the elastomers. The unmodeled elastomer strain will limit the motion of a FREE actuator to a small section of the possible fiber-angle-combinations defined by η. To account for the effects of the elastomer, a designer could take the insights from this fiber-only model and further explore them with fiber-elastomer models (e.g., see Refs. [6–9]).

The model presented in this work will facilitate the growing understanding of FREEs. The identification of the descriptor η allows the behavior of FREEs to be described independent of actuator-specific geometry. The linearly scaling functions make this behavior simple to predict and understand. As the understanding of FREEs grows, engineers will find new opportunities for these unique actuators to expand the functional ability of soft, fluid-driven systems.

Acknowledgment

The authors would like to acknowledge Joshua Bishop-Moser and Daniel Bruder for their thoughtful reviews of this work prior to submission and Audrey Sedal for helping the authors understand the constitutive models that have been used to describe FREEs in other works. This material is based upon work supported by the National Science Foundation Graduate Research Fellowship under Grant No. DGE 1256260. Any opinion, findings, and conclusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the National Science Foundation.

Funding Data

National Science Foundation Graduate Research Fellowship (Grant No. DGE 1256260).
Toyota Research Institute.

Nomenclature

b_α, b_β =: unwound length of the individual helical fibers in the family
D =: diameter of the cylindrical actuator
D₀ =: the diameter of the zero-height helical fibers when β = −π/2
l =: axial length of the cylindrical actuator
n_0,_α, n_0,_β =: the number of turns around the axis made by each fiber at the diameter D₀
n_α, n_β =: the number of turns around the axis made by each fiber at the current diameter
V =: volume of the cylindrical actuator
α =: angle of the family of helical fibers with equal or longer unwound length
α_LM, β_LM =: values of α and β on the locked manifold for a given value of η
β =: angle of the other family of helical fibers
γ =: the quotient of the sine of α and η
Δn =: the axial rotation between the actuator ends (β = −π/2 → Δn = 0)
η =: ratio of the lengths of the β fibers to that of the α fibers

References

1.

Daerden

,

F.

, and

Lefeber

,

D.

,

2002

, “

Pneumatic Artificial Muscles: Actuators for Robotics and Automation

,”

Eur. J. Mech. Environ. Eng.

,

47

(

1

), pp.

11

–

21

.http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.89.6717&rep=rep1&type=pdf

2.

Agerholm

,

M.

, and

Lord

,

A.

,

1961

, “

The Artificial Muscle of Mckibben

,”

Lancet

,

277

(

7178

), pp.

660

–

661

.

Google Scholar

Crossref

3.

Andrikopoulos

,

G.

,

Nikolakopoulos

,

G.

, and

Manesis

,

S.

,

2011

, “

A Survey on Applications of Pneumatic Artificial Muscles

,”

19th Mediterranean Conference on Control Automation

(

MED

), Corfu, Greece, June 20–23, pp.

1439

–

1446

.

4.

Bishop-Moser

,

J. L.

,

2014

, “

Design of Generalized Fiber-Reinforced Elasto-Fluidic Systems

,”

Ph.D. thesis

, University of Michigan, Ann Arbor, MI.https://deepblue.lib.umich.edu/handle/2027.42/107202

5.

Krishnan

,

G.

,

Bishop-Moser

,

J.

,

Kim

,

C.

, and

Kota

,

S.

,

2015

, “

Kinematics of a Generalized Class of Pneumatic Artificial Muscles

,”

ASME J. Mech. Rob.

,

7

(

4

), p.

041014

.

Google Scholar

Crossref

6.

Connolly

,

F.

,

Polygerinos

,

P.

,

Walsh

,

C. J.

, and

Bertoldi

,

K.

,

2015

, “

Mechanical Programming of Soft Actuators by Varying Fiber Angle

,”

Soft Rob.

,

2

(

1

), pp.

26

–

32

.

Google Scholar

Crossref

7.

Connolly

,

F.

,

Walsh

,

C. J.

, and

Bertoldi

,

K.

,

2017

, “

Automatic Design of Fiber-Reinforced Soft Actuators for Trajectory Matching

,”

Proc. Natl. Acad. Sci.

,

114

(

1

), pp.

51

–

56

.

Google Scholar

Crossref

8.

Sedal

,

A.

,

Bruder

,

D.

,

Bishop-Moser

,

J.

,

Vasudevan

,

R.

, and

Kota

,

S.

,

2017

, “

A Constitutive Model for Torsional Loads on Fluid-Driven Soft Robots

,”

ASME

Paper No. DETC2017-67970.

9.

Singh

,

G.

, and

Krishnan

,

G.

,

2017

, “

A Constrained Maximization Formulation to Analyze Deformation of Fiber Reinforced Elastomeric Actuators

,”

Smart Mater. Struct.

,

26

(

6

), p.

065024

.

Google Scholar

Crossref

10.

Chou

,

C.-P.

, and

Hannaford

,

B.

,

1996

, “

Measurement and Modeling of Mckibben Pneumatic Artificial Muscles

,”

IEEE Trans. Rob. Autom.

,

12

(

1

), pp.

90

–

102

.

Google Scholar

Crossref

11.

Singh

,

G.

, and

Krishnan

,

G.

,

2015

, “

An Isoperimetric Formulation to Predict Deformation Behavior of Pneumatic Fiber Reinforced Elastomeric Actuators

,”

IEEE/RSJ International Conference on Intelligent Robots and Systems

(

IROS

), Hamburg, Germany, Sept. 28–Oct. 2, pp.

1738

–

1743

.

12.

Gilbertson

,

M. D.

,

McDonald

,

G.

,

Korinek

,

G.

,

de Ven

,

J. D. V.

, and

Kowalewski

,

T. M.

,

2017

, “

Serially Actuated Locomotion for Soft Robots in Tube-Like Environments

,”

IEEE Rob. Autom. Lett.

,

2

(

2

), pp.

1140

–

1147

.

Google Scholar

Crossref

13.

Nakamura

,

T.

, and

Shinohara

,

H.

,

2007

, “

Position and Force Control Based on Mathematical Models of Pneumatic Artificial Muscles Reinforced by Straight Glass Fibers

,”

IEEE International Conference on Robotics and Automation

(

ICRA

), Rome, Italy, Apr. 10–14, pp.

4361

–

4366

.

A Closed-Form Kinematic Model for Fiber-Reinforced Elastomeric Enclosures

Introduction

Kinematic Modeling

Assumptions and Definitions.

Behavior Described by η, State by β.

Size Described by b_β and D₀.

Calculating Rotation Δn, Diameter D, Length l, Volume V, and Surface Area A_surf.

Example of Model Application: Actuator Design

Discussion

Acknowledgment

Funding Data

Nomenclature

References

Get Email Alerts

Cited By

ASME Journals

ASME Conference Proceedings

ASME eBooks

Resources

Opportunities

A Closed-Form Kinematic Model for Fiber-Reinforced Elastomeric Enclosures

Introduction

Kinematic Modeling

Assumptions and Definitions.

Behavior Described by η, State by β.

Size Described by bβ and D0.

Calculating Rotation Δn, Diameter D, Length l, Volume V, and Surface Area Asurf.

Example of Model Application: Actuator Design

Discussion

Acknowledgment

Funding Data

Nomenclature

References

Get Email Alerts

Cited By

Related Articles

Related Proceedings Papers

Related Chapters

ASME Journals

ASME Conference Proceedings

ASME eBooks

Resources

Opportunities

This Feature Is Available To Subscribers Only

Size Described by b_β and D₀.

Calculating Rotation Δn, Diameter D, Length l, Volume V, and Surface Area A_surf.