Uniaxial poroelastic tendon model with crimped fibre recruitment

Sintesi Tecnica: Modello Poroelastico Uniaxiale di Tendine con Reclutamento di Fibre Crimpatura

Definizione del Problema
I tendini sono tessuti molli complessi composti principalmente da fibrille di collagene, una matrice non collagenica (NCM) e acqua (55–70%). Sebbene il ruolo meccanico delle fibrille di collagene "crimpatura" (crimped)—dove le fibrille si raddrizzano e reclutano carico solo dopo essere state raddrizzate—sia ampiamente documentato sperimentalmente, gli attuali modelli matematici spesso non riescono a integrare questo comportamento di reclutamento all'interno di un quadro poroelastico. I precedenti modelli poroelastici hanno o ignorato completamente le fibrille, trattando lo scheletro solido come un materiale neo-Hookeano omogeneo, o hanno incluso le fibre senza tenere conto del reclutamento progressivo dovuto alla crimpatura. Questo studio affronta la lacuna nella comprensione di come il reclutamento delle fibre all'interno di un quadro poroelastico non lineare influenzi la risposta del tessuto alla sollecitazione e al rilascio del carico, specificamente riguardo ai tassi di diffusione, l'asimmetria stress-strain e l'isteresi.

Metodologia
Gli autori hanno sviluppato un modello poroelastico non lineare monodimensionale (il modello FIB) che rappresenta un tendine delimitato da muscolo e osso. Lo scheletro solido è modellato come un composito di:

1. Fibrille Crimpatura (Crimped Fibrils): Fibrille di collagene rigide, inizialmente crimpatura, immerse in una matrice. Esse non sopportano alcun carico finché non vengono raddrizzate (uncrimped). Il reclutamento è governato da una distribuzione dell'angolo di crimpatura $\theta(r)$, dove le fibrille si raddrizzano progressivamente man mano che il tessuto si estende.
2. Matrice Non Collagenica (NCM): Una matrice di fondo morbida, neo-Hookeana, che sopporta il carico continuamente.

Il modello assume grani solidi e fluidi incompressibili, con il flusso del fluido governato dalla legge di Darcy e dalla permeabilità dipendente dalla porosità (legge di Kozeny-Carman). L'equazione governante è un'equazione di diffusione non lineare per la porosità, dove il coefficiente di diffusione $D$ dipende dalla rigidità locale ($\partial s' / \partial \hat{\Phi}$).

Il modello FIB è confrontato con un modello di controllo (modello NH) in cui lo scheletro solido è un materiale puramente neo-Hookeano senza reclutamento di fibre, utilizzando moduli di Young macroscopici equivalenti. Lo studio investiga due condizioni al contorno all'interfaccia muscolare:

• Carico (Loading): Applicazione istantanea di uno stress tensile costante a un tendine rilassato.
• Scarico (Unloading): Rilascio istantaneo del carico da un tendine pre-stirato.

Le simulazioni sono state eseguite utilizzando un metodo a volumi finiti nello spazio e un metodo implicito a passo temporale variabile nel tempo.

Contributi Chiave

• Integrazione del Reclutamento nella Poroelasticità: Lo studio estende un precedente quadro poroelastico per includere esplicitamente una regola di reclutamento delle fibre basata sulla crimpatura del collagene, consentendo un confronto diretto tra un tessuto rinforzato con fibre e un tessuto neo-Hookeano omogeneo.
• Derivazione di una Legge di Stress Su Misura: Gli autori hanno derivato una specifica legge di stress per il modello FIB che tiene conto della transizione dallo stato di crimpatura (carico zero) allo stato di tensione (carico lineare), distinguendo i contributi di stress delle fibrille e della NCM.
• Identificazione dell'Accoppiamento Diffusione-Rigidità: Il lavoro evidenzia che il coefficiente di diffusione non è costante ma varia dinamicamente con il reclutamento delle fibre. Man mano che le fibre vengono reclutate, la rigidità aumenta, aumentando il tasso di diffusione; viceversa, mentre le fibre tornano in stato di crimpatura, la rigidità diminuisce, riducendo la diffusione.

Risultati

• Fase di Carico:

  • Risposta più Morbida: Il modello FIB mostra una risposta stress-strain più morbida rispetto al modello NH. Di conseguenza, sotto lo stesso carico applicato, il tendine FIB si estende maggiormente (maggiore deformazione/porosità a regime).
  • Approccio più Lento allo Stato Stazionario: Nonostante si estenda di più, il modello FIB si avvicina allo stato stazionario più lentamente rispetto al modello NH. Ciò accade perché la fase iniziale di reclutamento comporta una bassa rigidità (diffusione vicina allo zero), rallentando la propagazione della risposta del carico attraverso il tessuto.
  • Eterogeneità Spaziale: Il modello FIB presenta punti di svolta netti nei profili di stress, porosità e flusso nella posizione in cui le fibrille stanno diventando tese ($Z_c$), una caratteristica assente nei profili fluidi del modello NH.

• Fase di Scarico:

  • Asimmetria e Rilassamento Lento: Il modello FIB mostra una significativa asimmetria tra carico e scarico. Durante lo scarico, mentre il tessuto si rilassa, le fibrille riacquistano la crimpatura, causando un ammorbidimento del tessuto. Questo ammorbidimento riduce il coefficiente di diffusione, creando un ciclo di feedback che rallenta drasticamente il processo di rilassamento. Il tessuto spesso non riesce a raggiungere uno stato completamente rilassato entro il tempo di simulazione, lasciando una frazione non trascurabile di fibrille ancora tese.
  • Isteresi: Questa asimmetria si manifesta come un ciclo di isteresi nella curva media stress-strain. L'isteresi è più pronunciata ai carichi inferiori, dove le fibrille riacquistano la crimpatura precocemente nella fase di scarico. Ai carichi più elevati, le fibrille rimangono tese più a lungo, permettendo un rilassamento iniziale più veloce e una riduzione dell'isteresi, sebbene il sistema non raggiunga comunque il rilassamento completo entro il tempo imposto.

• Meccanismo Protettivo: Il modello dimosta che la presenza di fibrille crimpatura permette al tendine di subire estensioni complessive maggiori senza sollecitare eccessivamente né le fibrille né la NCM, agendo come un meccanismo di protezione naturale per le componenti strutturali.

Significatività e Rivendicazioni
L'articolo sostiene che l'integrazione del reclutamento delle fibre crimpatura in un quadro poroelastico fornis la interpretazione meccanica dei tendini più biologicamente realistica rispetto ai modelli neo-Hookeani omogenei. La principale significatività risiede nello spiegare l'origine dell'isteresi e dell'asimmetria carico/scarico non solo attraverso la viscoelasticità, ma attraverso l'interazione tra la natura diffusiva del flusso dei fluidi e la rigidità dipendente dal reclutamento delle fibre.

Gli autori osservano che i parametri del modello FIB (angolo di crimpatura, frazione di fibrille, modulo delle fibrille) sono più direttamente misurabili e biologicamente interpretabili rispetto alle costanti fenomenologiche richieste dalle leggi neo-Hookeane. Il modello riproduce con successo le osservazioni sperimentali delle regioni "toe" e dell'isteresi. Tuttavia, gli autori riconoscono con modestia alcuni limiti, tra cui l'assunzione che le fibrille crimpatura non sopportino alcun carico (il che potrebbe non essere strettamente vero), l'esclusione dello scorrimento delle fibre e la semplificazione in una geometria monodimensionale. Il framework è presentato come una base per lavori futuri sul danno tendineo (ad esempio, rottura delle fibrille) e sulla modellazione tridimensionale.