---
title: Fiber Reinforced Elastomeric Enclosures
url: https://www.emergentmind.com/topics/fiber-reinforced-elastomeric-enclosures-frees
type: topic
---

# Fiber Reinforced Elastomeric Enclosures

Searching arXiv for recent and foundational FREE papers to ground the article.
Fiber Reinforced Elastomeric Enclosures (FREEs) are soft, fluid-driven actuators consisting of an elastomeric tube reinforced by one or more families of helically wound fibers around a fluid-filled internal cavity. Under internal pressurization, the elastomer tends to expand, while the fibers impose geometric constraints that program coupled axial extension or contraction, twist, and radial expansion. Unlike rigid actuators, FREEs generate spatial forces and moments without imposing rigid kinematic constraints, which makes them particularly suitable for soft robots, continuum manipulators, exoskeletons, grippers, and modular parallel actuator arrangements [1805.00124][1902.00054].

## 1. Definition, structure, and distinguishing characteristics

In the formulation used across the cited literature, a FREE is a soft fluid-driven actuator built from an elastomeric cylindrical tube, a fiber reinforcement architecture, rigid end fittings, and an internal pressurized cavity. A common baseline design uses a single family of inextensible fibers wound helically at a prescribed fiber angle, denoted \(\Gamma\) or \(\alpha\), relative to the tube axis [1805.00124][1902.00054]. Other variants employ two helical families or an additional longitudinal strain-limiting fiber to convert extension into bending [2002.09577][2507.10121].

The defining mechanical feature of a FREE is that the fibers constrain the natural expansion of the elastomer in a geometrically programmed way. This coupling allows pressure-driven deformation modes including axial extension or contraction, twisting about the main axis, radial change, bending, coiling, or combinations thereof, depending on fiber layout and boundary conditions [1805.00124][2002.09577][2507.10121]. In the single-family cylindrical idealization, the principal design parameter is the fiber angle; changing this angle alters the balance between axial force and torsional moment, and more generally the relationship between pressure, deformation, and generated wrench [1805.00124].

FREEs differ from conventional rigid actuators in that they do not incorporate a rigid kinematic constraint such as a bearing-constrained shaft or lead screw. Instead, they can impart spatial forces and moments to attached bodies without restricting those bodies to the same directions of motion. This suggests a useful distinction between force generation and motion constraint: rigid actuators usually combine both in a single mechanism, whereas FREEs decouple them by remaining inherently compliant [1805.00124]. A plausible implication is that FREE arrays can be embedded into soft structures without introducing rigid joints while still producing multi-degree-of-freedom actuation.

The application space described in the literature is correspondingly broad. FREEs are reported as actuators in exoskeletons and wearable assistive devices, soft manipulators and tentacles, soft grippers and hands, locomoting soft robots and snake-like robots, vibration isolators, and parallel actuator arrangements [1902.00054]. In later assembly-level work, they are also treated as modular building blocks for soft continuum arms intended for assistive devices, agriculture, search applications, or surgery [2507.10121].

## 2. Geometry, kinematics, and deformation programming by fiber architecture

The simplest and most widely analyzed FREE geometry is a cylindrical tube with a single helical fiber family. Under the ideal assumptions of inextensible fibers, negligible wall thickness, cylindrical shape, and no bending or taper, the deformation state is commonly represented by generalized coordinates
\[
\vec{q}=
\begin{bmatrix}
\Delta l\\
\Delta \phi
\end{bmatrix},
\]
where \(\Delta l\) is change in length and \(\Delta \phi\) is twist about the main axis [1805.00124]. The associated generalized force vector is
\[
\vec{\tau}=
\begin{bmatrix}
F\\
M
\end{bmatrix},
\]
with \(F\) the axial force and \(M\) the torsional moment [1805.00124].

For a relaxed configuration with length \(L\), radius \(R\), and fiber angle \(\Gamma\), the fiber length and the number of fiber revolutions are given by
\[
B=\left|\frac{L}{\cos\Gamma}\right|,
\qquad
N=-\frac{L}{2\pi R}\tan\Gamma.
\]
These quantities remain fixed under the inextensible-fiber assumption [1805.00124]. The current length is
\[
l=L+\Delta l,
\]
and the fiber-length constraint yields the current radius
\[
r=\frac{B}{|2\pi N+\Delta\phi|\sqrt{1-\left(\frac{L+\Delta l}{B}\right)^2}}.
\]
The internal volume then follows as
\[
V(\vec{q})=\pi r^2 l
=\frac{\pi (L+\Delta l)\left(B^2-(L+\Delta l)^2\right)}{(2\pi N+\Delta\phi)^2}.
\]
Thus, under the ideal cylindrical model, length, twist, radius, and volume are not independent but are all determined by \(\vec{q}\) and the design parameters \((L,R,\Gamma)\) [1805.00124].

The fiber angle governs the qualitative actuation mode. Low \(|\Gamma|\), with fibers more aligned to the axis, favors axial contraction under pressurization and relatively smaller twist; high \(|\Gamma|\), with fibers more circumferential, allows more axial extension and often produces extension with twist [1805.00124]. In one formulation that includes pressure on the end caps, the ratio between axial force and twisting moment is
\[
\frac{F}{M'}=\frac{1-2\cot^2\Gamma}{2\cot\Gamma},
\qquad
M'=\frac{M}{r}.
\]
This establishes fiber angle as a direct programming variable for the force–moment balance [1805.00124].

A related geometric result appears in work on biologically inspired bending FREEs. For an extending FREE with relaxed fiber angle \(\alpha_0\), the behavior transitions at the “neutral angle”
\[
\arctan(\sqrt{2}) \approx 54.74^\circ.
\]
If \(|\alpha_0|<54.74^\circ\), the actuator is contracting under pressurization; if \(|\alpha_0|>54.74^\circ\), it is extending [2002.09577][2002.09577]. The same geometric threshold is discussed in constitutive studies of fiber-reinforced elastomers as the “magic angle” \(\alpha_m=\arctan\sqrt{2}\approx54.7^\circ\), the angle maximizing enclosed volume for fixed fiber length [1912.00989]. This suggests a deep connection between FREE actuation mode selection and angle-dependent anisotropic stress states in fiber-reinforced elastomers.

When a longitudinal strain-limiting fiber is added to an extending FREE, axial extension is converted into bending. Using the geometric relations
\[
L=\mathcal{R}\theta,\qquad
L_0=(\mathcal{R}-R)\theta,
\]
and the FREE radius–length relation
\[
R=\frac{\sqrt{B^2-L^2}}{2\pi N},
\]
the curvature is
\[
\mathcal{K}
=
2\pi N\left(1-\frac{L_0}{L}\right)\left(B^2-L^2\right)^{-1/2},
\]
and the maximum curvature becomes
\[
\mathcal{K}_{\max}
=
\frac{1}{R_0}
\left(\frac{\sin\alpha_0}{\sin54.74^\circ}\right)
\left(1-\frac{\cos\alpha_0}{\cos54.74^\circ}\right).
\]
In that framework, \(\mathcal{K}_{\max}\) is monotonically increasing in \(\alpha_0\) for \(\alpha_0\in[54.74^\circ,90^\circ)\), so larger fiber angles yield larger attainable curvature [2002.09577].

## 3. Force generation, fluid Jacobians, and pressure-to-wrench mappings

A central line of FREE research formulates actuation through a state-dependent fluid Jacobian. For a single cavity actuator with volume \(V=V(\vec{q})\), volumetric flow satisfies
\[
\dot{V}=\bar{J}_q(\vec{q})\,\dot{\vec{q}},
\qquad
\bar{J}_q(\vec{q})=\frac{\partial V}{\partial \vec{q}}.
\]
For the cylindrical single-family FREE with \(\vec{q}=[\Delta l,\Delta\phi]^T\),
\[
\bar{J}_q(\vec{q})=
\begin{bmatrix}
\displaystyle \frac{\pi(B^2-3(L+\Delta l)^2)}{(2\pi N+\Delta\phi)^2}
&
\displaystyle \frac{2\pi(L+\Delta l)\left((L+\Delta l)^2-B^2\right)}{(2\pi N+\Delta\phi)^3}
\end{bmatrix}.
\]
Using power balance between mechanical power and fluid power,
\[
\dot{\vec{q}}^{\,T}\vec{\tau}=p\dot{V},
\]
one obtains the pressure-to-force relation
\[
\vec{\tau}(\vec{q},p)=\bar{J}_q^T(\vec{q})\,p.
\]
For a given state, the generalized force direction is encoded in \(\bar{J}_q^T\), while the magnitude scales linearly with the scalar pressure \(p\) [1805.00124].

This model explicitly separates active fiber-generated forces from passive elastomeric forces:
\[
\vec{\tau}_{\text{total}}=\vec{\tau}+\vec{\tau}_{\text{elast}}.
\]
The pressure model concerns \(\vec{\tau}=\bar{J}_q^T p\), while passive elastomer response is treated separately or subtracted experimentally [1805.00124].

The same work generalizes the approach to arrays of \(n\) FREEs mounted in parallel to a common end effector. For actuator \(i\), the local generalized force \(\vec{\tau}_i=[F_i,M_i]^T\) is mapped into a 6D end-effector wrench by its attachment offset \(\vec{d}_i\) and axis direction \(\hat{a}_i\). The net wrench becomes
\[
\vec{f}(\vec{x},\vec{p})=\sum_{i=1}^n \bar{J}_{x,i}^T(\vec{x})\,p_i
=\bar{J}_x^T(\vec{x})\,\vec{p},
\]
where \(\vec{x}\) is end-effector state and \(\bar{J}_x(\vec{x})\) is the stacked fluid Jacobian in end-effector coordinates [1805.00124]. This formulation treats a parallel FREE assembly as a state-dependent linear map from pressure inputs to spatial wrench.

A closely related lumped model, developed for dynamic control of a single FREE, resolves pressure-induced axial force and torsional moment through thin-wall equilibrium. With linearized elastomer force and damping,
\[
F_e=-k_e s-c_e\dot{s},
\qquad
M_e=-k_t\phi-c_t\dot{\phi},
\]
the governing equations become
\[
F_1-k_e s-c_e\dot{s}+\pi r^2 P(1-2\cot^2\gamma)=m_1\ddot{s},
\]
\[
M_1-k_t\phi-c_t\dot{\phi}-2\pi r^3 P\cot\gamma=I_1\ddot{\phi}.
\]
Under low-pressure operating conditions for the studied actuator, simulations showed \(r\) and \(\gamma\) change less than \(2\%\), allowing a linearized control model with \(r\approx R\) and \(\gamma\approx\Gamma\) [1912.07380][2110.01591].

In more recent continuum-arm simulation work, FREE physics is embedded into a Cosserat rod model. There, the fiber-angle kinematics are written in terms of axial and radial stretches
\[
\lambda_1=\frac{l}{l_0},\qquad \lambda_2=\frac{r}{r_0},
\]
and current fiber angles satisfy
\[
\alpha = \tan^{-1}\left( \frac{\lambda_2}{\lambda_1 \tan(\alpha_0) + \dfrac{r}{l}\,\delta} \right),
\qquad
\beta = \tan^{-1}\left( \frac{\lambda_2}{\lambda_1 \tan(\beta_0) + \dfrac{r}{l}\,\delta} \right).
\]
The corresponding axial force and torsional couple are then taken from analytical FREE theory as functions of \(P\), \(r\), \(\alpha\), and \(\beta\), and applied as distributed loads to a Cosserat rod [2507.10121]. This suggests a hierarchy of models: Jacobian-based quasi-static wrench models, lumped dynamic end-cap models, and rod-based assembly models all use pressure-induced force transmission constrained by fiber geometry, but at different levels of abstraction.

## 4. Constitutive modeling and model classes

The constitutive modeling of FREEs spans several levels of fidelity. A comparative study evaluates three model classes for single fiber-family FREEs: a linear lumped-parameter model, a nonlinear continuum-mechanics model, and a neural network model [1902.00054]. The benchmark comprises 8 FREE designs, 22,880 pressure–kinematics configurations, and approximately 12,575 unbuckled data points [1902.00054].

The linear lumped-parameter model imposes an inextensible-fiber kinematic constraint and a thin-wall approximation, from which the volume
\[
V(\vec{q})=\pi\frac{l B^2-l^3}{(\Phi+\varphi)^2}
\]
and the fluid Jacobian
\[
J_V=
\begin{bmatrix}
\displaystyle \pi\frac{B^2-3l^2}{(\Phi+\varphi)^2}
&
\displaystyle -2\pi\frac{lB^2-l^3}{(\Phi+\varphi)^3}
\end{bmatrix}
\]
are derived [1902.00054]. With pressure-generated generalized force
\[
\vec{\tau}_{\text{fluid}}=J_V^T P
\]
and a linear wall stiffness
\[
\vec{\tau}_{\text{wall}}=-K\Delta\vec{q}
=
-\begin{bmatrix}
k_a & k_c\\
k_c & k_b
\end{bmatrix}
\begin{bmatrix}
l-L\\
\varphi
\end{bmatrix},
\]
the full model is
\[
\vec{\tau}=J_V^T P-K\Delta\vec{q}.
\]
This model uses only three experimentally identified scalars \(k_a,k_b,k_c\), but has a singularity at
\[
\varphi=-\Phi=-\frac{L}{R_o}\tan\Gamma,
\]
where the fibers become parallel to the axis [1902.00054].

The nonlinear continuum model treats the wall as an incompressible hyperelastic composite with finite strain. Using the kinematic mapping
\[
r=\sqrt{\frac{R^2-R_i^2}{\lambda_z}+r_i^2},\qquad
\theta=\Theta+Z\frac{\varphi}{L},\qquad
z=\lambda_z Z,\qquad
\lambda_z=\frac{l}{L},
\]
and a deformation gradient
\[
F=
\begin{bmatrix}
\displaystyle \frac{R}{r\lambda_z} & 0 & 0\\
0 & \displaystyle \frac{r}{R} & r\frac{\Phi}{L}\\
0 & 0 & \lambda_z
\end{bmatrix},
\]
it decomposes the Helmholtz free energy as
\[
\Psi_{\text{total}}=\Psi_{\text{isotropic}}+\Psi_{\text{anisotropic}},
\]
with
\[
\Psi_{\text{isotropic}}=\frac{C_1}{2}(I_1(F^T F)-3),
\qquad
\Psi_{\text{anisotropic}}=\frac{C_2}{2}(I_4(F^T F,\Gamma)-1)^2.
\]
Axial force and moment then follow from stress integration:
\[
F=-2\pi\int_{r_i}^{r_o}\sigma_{zz}\,r\,dr+\pi r_i^2 P,
\qquad
M=-2\pi\int_{r_i}^{r_o}\sigma_{\theta z}\,r^2\,dr.
\]
This model uses only two physical parameters, \(C_1\) and \(C_2\), corresponding to elastomer and fiber stiffness [1902.00054].

The neural network model maps \((l,\varphi,P,\Gamma,L,R_i,R_o)\) to \((F,M)\) using a shallow fully connected network with one hidden layer of 6 neurons, hyperbolic tangent activation, and 62 total parameters [1902.00054]. It achieves the best actuator-specific accuracy when trained and tested on the same sample, but does not generalize well across design space [1902.00054].

The comparison yields a characteristic trade-off. The neural net achieves the highest peak performance, while the first principles-based models generalize best across all actuator design parameters tested [1902.00054]. Quantitatively, the continuum model has the lowest average error across all training–test combinations, approximately \(18.1\%\), with moderate error range, whereas the lumped model has the lowest single-case error, \(2.19\%\), but an average error of approximately \(34.9\%\), and the neural network can reach diagonal errors around \(2.32\text{–}11.5\%\) but with errors up to approximately \(190\%\) under cross-design generalization [1902.00054]. A plausible implication is that constitutive structure is particularly valuable when FREE models are intended for new geometries rather than re-identification of one fabricated specimen.

Beyond actuator-level models, continuum constitutive work on fiber-reinforced elastomers supplies angle-dependent material laws relevant to FREE walls. A hyperelastic model of incompressible fiber-reinforced elastomer sheets proposes
\[
W = C_1 (I_1 - 3) + C_2 (I_1 - 3)^2
    + C_3 (I_1 - 3)^2 (I_4 - 1)
    + C_4 (I_4 - 1)^2 + C_5 (I_4 - 1)^3,
\]
with Cauchy stress
\[
\boldsymbol{\sigma}
= -p\mathbf{I} + 2 W_1 \mathbf{B} + 2 W_4\,\mathbf{m}\otimes\mathbf{m}.
\]
Experiments on PDMS composites with fiber orientations from \(0^\circ\) to \(90^\circ\) show that fibers are extensible and that the inextensibility constraint is invalid in those specimens [1912.00989]. This does not invalidate inextensible-fiber FREE models as ideal design tools, but it indicates that high-fidelity FREE wall simulations may require extensible-fiber constitutive laws rather than hard geometric constraints.

A further extension introduces finite-strain viscoelastic models for fiber-reinforced composites using multiplicative decomposition of the deformation gradient and fiber stretch. That framework includes Mooney–Rivlin or Neo–Hookean matrix models, Holzapfel-type fiber potentials, modified potentials for compressive buckling and slackness, and efficient time stepping for matrix and fiber viscoelasticity [1806.05948]. Although developed for pressurized multi-layer composite pipes, the formulation maps directly onto FREE-like structures under large deformation and offers a route to modeling creep, stress relaxation, and hysteresis [1806.05948].

## 5. Parallel combinations, modules, and force zonotopes

FREEs are especially important as components of parallel soft actuation systems because their wrenches superimpose without introducing rigid joints [1805.00124]. In a parallel array, each actuator contributes only along the half-ray generated by positive pressure. Because a single actuator cannot reverse its force direction without negative pressure, multi-dimensional controllability requires antagonistic arrangements [1805.00124].

The wrench set achievable at a given state \(\vec{x}\) under pressure bounds \(p_i\in[0,p_i^{\max}]\) is the force zonotope
\[
\mathcal{Z}(\vec{x})
=
\left\{
\bar{J}_x^T(\vec{x})\vec{p} : p_i\in[0,p_i^{\max}]
\right\}.
\]
Equivalently, with actuator direction vectors
\[
\mathbf{w}_i(\vec{x})=\bar{J}_{x,i}^T(\vec{x})\,p_i^{\max},
\]
\[
\mathcal{Z}(\vec{x})
=
\left\{
\sum_{i=1}^n \alpha_i \mathbf{w}_i(\vec{x}) : \alpha_i\in[0,1]
\right\},
\]
which is also the convex hull of the \(2^n\) extreme points generated by \(p_i\in\{0,p_i^{\max}\}\) [1805.00124]. The zonotope is used as a design and evaluation tool because it respects the directional nature of each actuator and the positivity of pressure bounds [1805.00124].

Several design implications follow directly. Coverage of a desired wrench set \(S\) requires \(S\subseteq\mathcal{Z}(\vec{x})\) at all relevant states. The zonotope also reveals whether an arrangement has sufficient antagonism to generate both signs of each wrench component, whether wrench capabilities are isotropic or anisotropic, and how actuation degenerates across the workspace as \(\bar{J}_x(\vec{x})\) changes [1805.00124]. The work argues that in an \(m\)-dimensional wrench space, at least \(m+1\) suitably arranged actuators are needed for full control under positive-pressure constraints [1805.00124].

An experimental example uses three FREEs in parallel, all aligned with the end-effector \(z\)-axis and attached symmetrically around it. The three actuators are one extending FREE with \(\Gamma=-85^\circ\) and two contracting FREEs with \(\Gamma=48^\circ\) and \(\Gamma=-48^\circ\), all with relaxed \(R=5\,\text{mm}\) and \(L=100\,\text{mm}\) [1805.00124]. This arrangement yields a well-balanced zonotope in the 2D wrench space \((F^{\hat{z}_e},M^{\hat{z}_e})\), enabling full 2D control with the minimal three actuators [1805.00124].

Module-level studies in finite elements examine assemblies of four FREEs arranged at the corners of a square and attached to a common end-effector [1912.07380][2110.01591]. In one LR configuration, diagonally opposite actuators have opposite handedness, allowing selective generation of extension, twist, or bending via different pressure patterns. In the reported cases, equal pressurization of all four FREEs produces pure axial motion, diagonal pairs produce net twist, and adjacent pairs produce approximate bending [1912.07380]. Finite element results further show that a 30-degree difference in winding angle dramatically alters the shape of the workspace generated by four FREEs assembled into a module [2110.01591].

A more recent assembly framework uses multiple FREEs as Cosserat rods connected in parallel and serial combinations. It studies a BR2 arm, composed of one bending FREE and two twisting FREEs in parallel, and a two-segment BR2-B3 or B3-BR2 arm with six FREEs total [2507.10121]. This work does not present a zonotope analysis, but it demonstrates the same general principle: distinct FREE types can be combined into modules that generate rich 3D bending, twisting, and stretching through coordinated pressure inputs [2507.10121].

## 6. Experimental validation, applications, and limitations

Experimental validation of FREE models spans force prediction, dynamic control, morphology emulation, and assembly-level reconstruction. In the foundational parallel-actuation study, a 2-DOF test platform evaluated 125 pressure triplets at each of four end-effector configurations for two independent sets of FREEs. The active wrench was isolated by subtracting the zero-pressure measurement,
\[
\vec{f}(\vec{x},\vec{p})=
\vec{f}_{\text{meas}}(\vec{x},\vec{p})
-
\vec{f}_{\text{meas}}(\vec{x},\vec{0}),
\]
and compared against \(\bar{J}_x^T(\vec{x})\vec{p}\) [1805.00124]. The reported prediction errors were less than \(1.5\) N force and less than \(8\times 10^{-3}\) Nm moment in RMSE [1805.00124].

The comparative actuator-modeling study generated benchmark data from eight FREE samples spanning a broad design space, with fiber angles
\[
15^\circ, 25^\circ, 36^\circ, 40^\circ, 50^\circ, 62^\circ, 73^\circ, 76^\circ
\]
and approximately \(L=90\text{–}129\) mm, \(R_i\approx4.77\) mm, \(R_o\approx6.13\text{–}6.74\) mm [1902.00054]. Each sample was tested over a grid of axial displacement, twist, and pressure, then all buckled configurations were excluded, leaving approximately \(55\%\) of the configurations as usable unbuckled data [1902.00054]. Buckling occurred in approximately \(45\%\) of tested configurations, with no simple pattern [1902.00054].

A dynamical control study demonstrates that a simple lumped-parameter model can support pressure-based rotational control of a single FREE. With a closed-loop controller embedded in the system, the lumped model predicts the actual rotational motion of a FREE with at most \(4\%\) error [2110.01591]. In thesis-level detail, step-input rotation experiments yield approximately \(4\%\) RMS deviation between simulation and experiment, while polynomial trajectory following yields approximately \(3\%\) RMS deviation [1912.07380]. This suggests that low-order models remain effective for control design even when they abstract away much of the wall mechanics.

Application studies illustrate how FREE geometry can be tuned to match biological motion. In a snake-mimic platform, extending-type FREEs with strain-limiting fibers were configured as head, midsection, and tail segments to emulate anti-predator thrashing in three snake genera. The curvature values of the fabricated soft-robotic head, midsection, and tail segments overlap with those exhibited by live snakes, and robot motion durations were less than or equal to those of snakes for all three genera [2002.09577]. That work uses fiber angles \(67.5^\circ\) and \(89^\circ\) in different modules, together with longitudinal strain limiters, to produce straight, kinked, curved, or coiled segments [2002.09577].

Assembly-level experiments validate rod-based simulation against optical motion capture and physics-informed reconstruction. In serial twist–twist and bend–bend assemblies, as well as a full BR2-B3 arm, errors in tip position, total twist, total bending, and total elongation are typically below \(10\%\) [2507.10121]. The same digital twin framework was used for interactive forward control of a six-FREE arm manipulating through a constrained opening, with successful retrieval in 20 out of 20 real-world trials using replayed pressure sequences [2507.10121].

Several limitations recur across the literature. Ideal cylindrical models fail when FREEs bend, buckle, kink, or taper [1805.00124]. Buckling was observed experimentally, especially for a \(-85^\circ\) FREE at higher pressures, reducing effective axial force transmission [1805.00124]. More generally, buckled states are excluded from many benchmark analyses because the standard models assume cylindrical, unbuckled geometries [1902.00054]. Passive elastomeric forces, hysteresis, fluid compressibility, rate dependence, and friction or fiber slip are also typically neglected or approximated [1805.00124][1902.00054][1912.07380].

A common misconception is that fiber geometry alone determines FREE behavior. The evidence is more specific. Fiber angle is the dominant kinematic programming parameter, but finite element results indicate that variations in the material properties of the elastic enclosure of a FREE are more significant than variations in fiber properties [2110.01591]. In module studies and single-actuator FEA, elastomer stiffness changes deformation much more strongly than fiber stiffness once the fibers are already much stiffer than the matrix [1912.07380][2110.01591]. Another misconception is that inextensible-fiber models are universally exact; planar constitutive experiments on fiber-reinforced elastomers show measurable fiber extensibility and angle-dependent fiber strain [1912.00989]. The idealized inextensible-fiber assumption is therefore best understood as a modeling abstraction that can be highly useful but is not universally literal.

The broader trajectory of FREE research points toward increasingly integrated models. Quasi-static fluid Jacobians support wrench design and controllability analysis [1805.00124]. Continuum and constitutive models support generalization across geometries [1902.00054]. Finite-strain viscoelastic and wave models indicate how damping, compressive buckling, and transient wall mechanics might be incorporated [1806.05948][1909.05627]. Rod-based assembly simulation enables digital twins for modular continuum arms [2507.10121]. This suggests a layered methodology for future work: use low-order geometric models for design intuition and control allocation, continuum constitutive models for design transfer and material selection, and assembly-level simulators for task-level behavior in complex soft robotic systems.

Source: https://www.emergentmind.com/topics/fiber-reinforced-elastomeric-enclosures-frees