| name | soft-tissue-mechanics |
| description | Soft tissue biomechanics — hyperelastic models (Fung exponential, Holzapfel-Gasser-Ogden anisotropic, Mooney-Rivlin), viscoelasticity (quasi-linear viscoelastic model, Prony series), failure criteria (tissue toughness, ultimate stress), tendon, ligament, cartilage, skin, arterial wall mechanics, preconditioning, toe region (collagen uncrimping), anisotropy from fiber orientation, and medical device design implications. |
| metadata | {"priority":7,"promptSignals":{"phrases":["soft tissue mechanics","tissue biomechanics","tendon mechanics","arterial wall","Fung exponential","Holzapfel"],"minScore":3}} |
Soft Tissue Biomechanics — Complete Skill
Soft Tissue Structure and Composition
Hierarchical Structure
Soft tissues (tendon, ligament, skin, arteries, cartilage):
Composed of: collagen fibers (primary load-bearing), elastin (large-deformation recovery), ground substance (proteoglycans + water), cells (fibroblasts, chondrocytes, SMCs)
Collagen hierarchy (tendon):
Tropocollagen molecules → microfibrils → fibrils (diameter 50–500 nm) → fascicles (100–500 μm) → tendon (mm scale)
Crimp angle: fibrils have sinusoidal crimp → responsible for toe region in stress-strain curve
Fiber orientation:
Tendons/ligaments: predominantly aligned fibers (transversely isotropic)
Arterial wall: two families of collagen fibers at ±angle to circumferential direction
Skin: random fiber network → nearly isotropic in-plane
Cartilage: zone-dependent orientation (superficial: parallel; middle: random; deep: perpendicular)
Mechanical Behavior
Non-Linear Stress-Strain Curve
Three distinct regions:
- Toe region (0–3% strain): collagen fibers uncrimp; low stiffness; dominated by elastin and ground substance
σ ≈ 0 to 1 MPa; tangent modulus E_toe = 0.01–0.5 MPa
- Heel region (transition): fibers progressively aligning and loading
- Linear region: crimped fibers fully straightened; fibers carry load; steeper slope
Tangent modulus E_lin = 1–2 GPa (tendon); 100–500 MPa (ligament); 0.5–10 MPa (skin)
- Failure: progressive fiber breakage; ultimate stress and failure (abrupt or gradual)
Tissue ultimate properties:
| Tissue | σ_ult [MPa] | ε_ult [%] | E_linear [MPa] | Notes |
|---|
| Tendon (patellar) | 50–100 | 8–15 | 1,000–2,000 | High strength; low compliance |
| Ligament (ACL) | 15–40 | 15–20 | 100–400 | More compliant than tendon |
| Aorta (circumferential) | 1–2 | 50–80 | 0.5–1.5 | Large deformation elastic |
| Skin | 5–30 | 20–80 | 15–150 | Anisotropic; highly variable |
| Cartilage (compression) | 0.5–5 | 10–20 | 0.1–1.0 (aggregate) | Biphasic; fluid-solid |
Hyperelastic Constitutive Models
Fung Exponential Model (Soft Tissues)
1D exponential model:
P(λ) = A × (exp(B × (λ² - 1)) - 1) [A, B = material constants; λ = stretch ratio]
σ_eng = P; σ_true = P × λ [Cauchy/true stress]
Full Fung exponential (2D):
W = c/2 × (exp(Q) - 1) [c = material constant; Q = quadratic form in strains]
Q = a₁×E₁₁² + a₂×E₂₂² + 2a₃×E₁₁×E₂₂ + a₄×E₁₂² [Green-Lagrange strains; a_i = material constants]
Stress: S_ij = ∂W/∂E_ij = c × exp(Q) × ∂Q/∂E_ij
Fitting: uniaxial tension + biaxial inflation data; least squares minimize ||P_model - P_data||²
Typical values (aorta): c = 1.0–10 kPa; a₁ = 0.5–2.0; a₂ = 0.5–2.0; a₃ = 0.1–1.0
Holzapfel-Gasser-Ogden (HGO) Anisotropic Model
Best model for fibered soft tissues (arteries, ligaments):
W = W_isotropic + W_anisotropic
W_iso = μ/2 × (I₁ - 3) [Neo-Hookean for matrix/ground substance; μ = shear modulus of matrix]
W_aniso = Σ_{α=1,2} k₁/(2k₂) × {exp(k₂ × (I₄α - 1)²) - 1} [two fiber families; I₄α = invariant for fiber direction]
[k₁ [kPa] = fiber stiffness parameter; k₂ = dimensionless exponential shape; I₄α = (a₀α · C · a₀α) = stretch² in fiber direction]
Fiber direction invariants:
I₄₁ = a₀₁ · C · a₀₁ = λ₁² [fiber 1; C = right Cauchy-Green; a₀α = reference fiber direction unit vector]
I₄₂ = a₀₂ · C · a₀₂ = λ₂² [fiber 2; symmetric about tissue axis]
Stress: S_aniso = 2 × Σ k₁ × (I₄α - 1) × exp(k₂×(I₄α-1)²) × a₀α⊗a₀α [2nd Piola-Kirchhoff]
Typical HGO parameters (human aortic adventitia):
μ = 10 kPa; k₁ = 5 kPa; k₂ = 10; fiber angles ±41.9° from circumferential direction
Dispersion of fibers (von Mises distribution — HGO-C model):
κ = 0: perfectly aligned fibers; κ = 1/3: isotropic distribution
I₄α_modified = κ×I₁ + (1-3κ)×I₄α [accounts for fiber dispersion; Gasser 2006]
κ = 0.226 for human aortic adventitia
Mooney-Rivlin for Soft Tissue
W = C₁₀(I₁-3) + C₀₁(I₂-3) — sometimes used for simple soft tissue (skin, fat)
Less accurate than Fung/HGO for fibered tissues; no toe region representation
C₁₀ = 0.1–5 kPa; C₀₁ = 0.1–2 kPa (soft hydrated tissues)
Quasi-Linear Viscoelasticity (QLV — Fung 1972)
Theory
Stress relaxation: apply step strain → stress decays with time
Creep: apply step stress → strain increases with time
QLV model:
σ(ε, t) = G(t) ∗ ∂σ_e(ε)/∂t [convolution of elastic response σ_e(ε) with reduced relaxation function G(t)]
Reduced relaxation function (Prony series):
G(t) = G_∞ + Σ G_i × exp(-t/τ_i) [G_∞ = equilibrium modulus ratio; G_i, τ_i = Prony coefficients]
Constraint: G_∞ + Σ G_i = 1
Continuous spectrum:
G(t) = 1 + C × [E₁(t/τ₂) - E₁(t/τ₁)] [Fung 1972 original; C, τ₁, τ₂ from relaxation test]
E₁ = exponential integral; typical: C = 0.025–0.2; τ₁ = 0.01 s; τ₂ = 100 s (broad relaxation spectrum)
Hysteresis:
Preconditioning: first several cycles shift mean response → stabilized behavior
Hysteresis energy = area within stress-strain loop = mechanical energy dissipated per cycle
Preconditioning
Preconditioning: repeated cyclic loading before testing to obtain repeatable mechanical response
Mechanism: progressive fiber re-arrangement and internal stress redistribution
Number of cycles: typically 10–20 cycles until response stabilizes
Effect: stiffness increases, toe region shortens (fibers become less crimped)
Design implication: use preconditioned properties for device design if tissue is chronically loaded (arteries, ligaments); use unconditioned for single-event failure (trauma prediction)
Tissue-Specific Notes
Arterial Wall
Structural layers (adventitia, media, intima):
Adventitia: predominantly collagen; prevents overstretch; stiff at high P
Media: SMC + collagen + elastin; elastic behavior; controls vasoactive response
Intima: thin (endothelium); fluid barrier; negligible structural role (increases in disease)
In-vivo state: arteries are under residual stress (opening angle test); residual stretch λ_z = 1.05–1.10 axially; circumferential stretch λ_θ = 1.3–1.7 at physiological pressure
Laplace: P × r / t = σ_hoop [circumferential wall stress; P ≈ 120 mmHg = 16 kPa; r/t ≈ 10 → σ = 160 kPa]
Cartilage (Biphasic Theory — Mow)
Two phases: solid matrix (collagen + proteoglycan); fluid (water, electrolytes = 70–80% total volume)
Biphasic equations: σ_total = σ_solid - φ_f × p × I [p = fluid pressure; φ_f = fluid fraction]
Permeability k: controls fluid exudation under compression: k = k₀ × (1 + e/e₀)^M [e = void ratio]
Aggregate modulus H_A = E(1-ν)/[(1+ν)(1-2ν)]; typical H_A = 0.5–1.0 MPa; k = 10⁻¹⁵ m²/Pa·s
Tendon Fatigue
Stress-life (S-N) for tendon:
σ_fatigue_limit ≈ 20–30% σ_ult (tissue; much lower than metal)
Damage accumulates with cycling: microtearing at ε > 4% in vitro
Medical Device Design Implications
Suture anchors, stents, orthopedic implants: must accommodate large tissue deformations without stress shielding or fretting
Catheter design: must navigate through arteries (R_min ≈ 5–10 mm for coronary; compliance matching)
Tissue engineering scaffolds: match mechanical properties of target tissue (stiffness, anisotropy, viscoelasticity)
Finite element simulation: use HGO or Fung model; incompressible (J=1 for soft tissue up to 10% strain)
Standards and References
| Standard | Scope |
|---|
| ASTM F2150 | Characterization of hydrogels for tissue-engineered constructs |
| ISO 10993-18 | Chemical characterization of materials in medical devices |
| Fung "Biomechanics: Mechanical Properties of Living Tissues" | Primary reference textbook |
| Holzapfel "Nonlinear Solid Mechanics" | HGO model development |
| Mow & Ratcliffe "Structure and Function of Articular Cartilage" | Biphasic cartilage theory |
Output
Provide: tissue type (tendon/ligament/artery/skin/cartilage; species; anatomical location; in-vivo vs. ex-vivo), loading mode (uniaxial tension/compression/biaxial/inflation; strain rate [/s]; temperature [°C]), material model selected (Fung/HGO/Mooney-Rivlin; basis for selection — fiber orientation data? IVUS/histology?), model parameters (c, a_i for Fung; μ, k₁, k₂, κ, fiber angle for HGO; from literature or fit to provided data; R² of fit), toe-to-linear transition (ε_toe [%]; E_toe [MPa]; E_linear [MPa]; transition stretch λ_transition), failure properties (σ_ult [MPa]; ε_ult [%]; failure mode: fiber rupture/matrix failure), viscoelastic properties (if needed: G_∞; Prony terms G_i/τ_i; hysteresis [%] from loop area), preconditioning protocol (N_cycles; stabilized properties used), medical device implication (stress shielding check; interface stress [kPa]; compliance mismatch factor), FEM recommendation (ABAQUS with UMAT for HGO; incompressible hyperelastic; mesh convergence for stress concentration), and reference (Fung 1993 textbook; Holzapfel 2000; tissue-specific measurement paper).
