Uniaxial poroelastic tendon model with crimped fibre recruitment

Technische Samenvatting: Uniaxiaal Poroelastisch Tendon-model met Crimped Fibre Recruitment

Probleemstelling
Pezen zijn complexe zachte weefsels die voornamelijk bestaan uit collageenfibrillen, een niet-collageen matrix (NCM) en water (55–70%). Hoewel de mechanische rol van "crimped" (gekrulde) collageenfibrillen—waarbij fibrillen recht worden en pas belasting opnemen nadat ze zijn ontkruld—experimenteel goed gedocumenteerd is, schieten bestaande wiskundige modellen vaak tekort in het integreren van dit recruitmentgedrag binnen een poroelastisch kader. Eerdere poroelastische modellen hebben ofwel fibrillen volledig genegeerd door het vaste skelet te behandelen als een homogeen neo-Hookeaans materiaal, of hebben fibrillen opgenomen zonder rekening te houden met de progressieve recruitment door crimp. Deze studie adresseert de kloof in het begrip van hoe vezel-recruitment binnen een niet-lineair poroelastisch kader de respons van het weefsel op belasting en ontlasting beïnvloedt, specifiek met betrekking tot diffusiesnelheden, stress-rek-asymmetrie en hysteresis.

Methodologie
De auteurs ontwikkelden een eendimensionaal, niet-lineair poroelastisch model (het FIB-model) dat een pees representeert die begrensd wordt door spier en bot. Het vaste skelet wordt gemodelleerd als een composiet van:

1. Crimped Fibrils: Stijve, initieel gekrulde collageenfibrillen ingebed in een matrix. Zij dragen geen belasting totdat ze gestrekt (ontkruld) zijn. De recruitment wordt beheerst door een distributie van crimp-hoeken $\theta(r)$, waarbij fibrillen progressief recht worden naarmate het weefsel uitrekt.
2. Non-Collagenous Matrix (NCM): Een zachte, neo-Hookeaanse achtergrondmatrix die continu belasting draagt.

Het model gaat uit van incompressibele vaste korrels en vloeistof, waarbij de vloeistofstroom wordt beheerst door de wet van Darcy en een porositeitsafhankelijke permeabiliteit (Kozeny-Carman wet). De leidende vergelijking is een niet-lineaire diffusievergelijking voor porositeit, waarbij de diffusiecoëfficiënt $D$ afhankelijk is van de lokale stijfheid ($\partial s' / \partial \hat{\Phi}$).

Het FIB-model wordt vergeleken met een controlemodel (NH-model) waarbij het vaste skelet een puur neo-Hokaans materiaal is zonder vezel-recruitment, gebruikmakend van equivalente macroscopische Young-moduli. De studie onderzoekt twee randvoorwaarden bij het spierinterface:

• Belasting: Instantane toepassing van een constante trekspanning op een ontspannen pees.
• Ontlasting: Instantane vrijgave van de belasting van een voorverstrekte pees.

Simulaties werden uitgevoerd met een eindige volumemethode in de ruimte en een impliciete methode met variabele stapgrootte in de tijd.

Belangrijkste Bijdragen

• Integratie van Recruitment in Poroelasticiteit: De studie breidt een eerder ontwikkeld poroelastisch kader uit door expliciet een vezel-recruitmentregel op basis van fibril crimp op te nemen, wat een directe vergelijking mogelijk maakt tussen een vezelversterkt weefsel en een homogeen neo-Hokaans weefsel.
• Afleiding van een Specifieke Stresswet: De auteurs hebben een specifieke stresswet afgeleid voor het FIB-model die rekening houdt met de overgang van de gekrulde staat (nul belasting) naar de gespannen staat (lineaire belasting), waarbij een onderscheid wordt gemaakt tussen de stressbijdragen van de fibrillen en de NCM.
• Identificatie van Diffusie-Stijfheid Koppeling: Het werk benadrukt dat de diffusiecoëfficiënt niet constant is maar dynamisch varieert met fibril-recruitment. Naarmate fibrillen rekruteren, neemt de stijfheid toe, wat de diffusiesnelheid verhoogt; omgekeerd, wanneer fibrillen opnieuw krullen (re-crimp), daalt de stijfheid, waardoor de diffusie afneemt.

Resultaten

• Belastingsfase:

  • Zachtere Respons: Het FIB-model vertoont een zachtere stress-rek respons dan het NH-model. Hierdoor rekt de FIB-pees onder dezelfde toegepaste belasting verder uit (hogere steady-state rek/porositeit).
  • Langzamere Benadering van Steady State: Ondanks het grotere rekvermogen, nadert het FIB-model de steady state langzamer dan het NH-model. Dit komt omdat de initiële recruitmentfase gepaard gaat met een lage stijfheid (nabij nul diffusie), wat de voortplanting van de belasting door het weefsel vertraagt.
  • Ruimtelijke Heterogeniteit: Het FIB-model vertoont scherpe kantelpunten in stress-, porositeits- en fluxprofielen op de locatie waar de fibrillen net taut (gespannen) worden ($Z_c$), een kenmerk dat afwezig is in de vloeiende profielen van het NH-model.

• Ontlastingsfase:

  • Asymmetrie en Trage Relaxatie: Het FIB-model vertoont een significante asymmetrie tussen belasting en ontlasting. Tijdens de ontlasting, naarmate het weefsel ontspant, gaan de fibrillen opnieuw krullen (re-crimp), wat het weefsel doet verzachten. Deze verzachting vermindert de diffusiecoëfficiënt, wat een feedbackloop creëert die het relaxatieproces drastisch vertraagt. Het weefsel bereikt vaak niet volledig een ontspannen staat binnen de simulatietijd, waarbij een niet-verwaarloosbaar deel van de fibrillen nog steeds gespannen blijft.
  • Hysteresis: Deze asymmetrie manifesteert zich als een hysteresis-lus in de gemiddelde stress-rek curve. De hysteresis is het meest uitgesproken bij lagere toegepaste belastingen waar fibrillen vroeg in de ontlastingsfase opnieuw gaan krullen. Bij hogere belastingen blijven de fibrillen langer gespannen, wat zorgt voor een snellere initiële relaxatie en minder hysteresis, hoewel het systeem nog steeds niet volledig ontspant binnen de gestelde tijdspanne.

• Beschermingsmechanisme: Het model toont aan dat de aanwezigheid van gekrulde fibrillen de pees in staat stelt om grotere totale extensies te ondergaan zonder de fibrillen of de NCM excessief te belasten, wat fungeert als een natuurlijk beschermingsmechanisme voor structurele componenten.

Betekenis en Claims
Het artikel claimt dat het integreren van gekrulde vezel-recruitment in een poroelastisch kader zorgt voor een biologisch realistischere interpretatie van de peesmechanica dan homogene neo-Hokaanse modellen. De primaire betekenis ligt in het verklaren van de oorsprong van hysteresis en de belasting/ontlastings-asymmetrie, niet enkel door viscoelasticiteit, maar door de interactie tussen de diffuse aard van vloeistofstroom en de stijfheidsafhankelijke recruitment van vezels.

De auteurs merken op dat de FIB-modelparameters (crimp-hoek, fibril fractie, fibril modulus) directer meetbaar en biologisch interpreteerbaar zijn dan de fenomenologische constanten die vereist zijn voor neo-Hokaanse wetten. Het model reproduceert succesvol experimentele observaties van de 'toe regions' en hysteresis. De auteurs erkennen echter bescheiden de beperkingen, waaronder de aanname dat gekrulde fibrillen nul belasting dragen (wat strikt genomen niet waar kan zijn), de exclusie van vezel-glijden (fibre sliding), en de vereenvoudiging naar een eendimensionale geometrie. Het kader wordt gepresenteerd als een fundament voor toekomstig werk over peesschade (bijv. het breken van fibrillen) en hogere-dimensionale modellering.