PDF version: Notes on fiber biomechanics – Logan Thrasher Collins

**Elastic fiber models**

- For an elastic fiber in which a linear relationship between force and change in length is assumed, the force is given by F = k(L – L
_{0}). - To normalize for other elastic fibers with different starting lengths, this equation is divided by L
_{0}to give F = k(L/L_{0}– 1). It is common practice to represent L/L_{0}as a parameter λ (called the stretch ratio). - As such, F is found using the formula below. Note that the quantity λ – 1 is referred to as the strain.

- While linear models are often useful, many real fibers exhibit finite extensibility (a nonlinear phenomenon) after exceeding a certain critical strain value λ
_{c}. That is, the force necessary to extend the fiber farther after exceeding λ_{c}increases rapidly. Finite extensibility can be modeled using the following equation which divides k by a term dependent on λ and λ_{c}.

- To model a muscle, let L
_{0}represent the muscle’s length in its inactive state and L_{contracted}represent the muscle’s length in its contracted state. Unlike the spring, the contracted state is used as the reference length. The contraction stretch is described by the ratio λ_{contracted}= L_{contracted}/L_{0}while the stretch ratio remains as λ = L/L_{0}.

- If this muscle is contracted without carrying a load such that F = 0, then λ = λ
_{contracted}. If the muscle acquires a load and so must maintain a constant length equal to its original length L_{0}(to “hold the load steady”), then the force in the muscle is F = k(1/λ_{contracted}– 1). - To generalize this model for 3-dimensional space, the locations of the fiber’s endpoints A and B are used. The fiber’s length and orientation are given below.

- The following force vectors can act on point B and on point A. The stretch ratio is still λ = L/L
_{0}.

**Viscous fiber models**

- Purely viscous behavior (as with liquids) can be described 1-dimensionally using the equation below where c
_{η}is a damping coefficient.

- The normalized rate of deformation is equivalent to the above formula without the damping coefficient. In addition, the rate of deformation can be written in terms of the stretch ratio λ = L/L
_{0}.

- If one endpoint of a filament of fluid is moved with constant velocity, its position is given by x
_{B}= L_{0}+ vt. This means that the rate of deformation is v/x_{B}.

- Solving the above equation gives the following result. For a constant rate of elongation, the point x
_{B}must be displaced exponentially over time.

- Given endpoint displacements u
_{A}and u_{B}, the total displacement is ΔL = u_{B}– u_{A}. Using this quantity, the stretch λ and the strain ε can be written using the equations below.

- For the small strains (i.e. |ε| is much less than 1) that result from small stretches, some approximations can be made which reduce the force equation to the following form.

**Viscoelastic fiber models**

- Many biological materials exhibit viscoelastic behavior rather than elastic behavior. In viscoelastic systems, the force on a fiber with a constant length decreases over time and applying constant force causes the length to increase.
- The strain response of a fiber to force is given as ε(t). The creep function J(t) describes the fiber’s tendency to permanently deform as strain is applied. Many possible creep functions can be devised depending on the system. The creep function is related to the strain by a factor of force increase F
_{0}.

- Strain responses follow the principle of superposition. That is, if a force is applied at τ
_{1}and then another force is applied at τ_{2}, a total strain response can be expressed as a sum of the two individual strain responses.

- For an arbitrary history of applied forces, an integral formulation of strain response is used. In this equation, the change in force at each infinitesimal time interval is multiplied by the creep function.

- Similarly, the force resulting from an imposed strain history can be expressed as an integral where G is a function that describes the relaxation of the fiber with time (analogous to the creep function, but the opposite concept).

**Reference: **Oomens, C., Brekelmans, M., & Baaijens, F. (2009). *Biomechanics: Concepts and Computation*. Cambridge University Press.

** **