Thèse APPROCHE DÉTERMINISTE DU SÉCHAGE DES AVIVÉS ...

Page-111-
L’outil numérique
En revanche, dans le cadre d’un séchage dissymétrique, la première hypothèse est maintenue, mais
le rayon de courbure est laissé variable pour qu’il puisse intervenir dans l’équilibre mécanique de
l’échantillon.
Début
- Maillage de l’épaisseur de la planche
- Conditions initiales et aux limites (transferts / mécanique)
- Actualisation du temps : i+ 1 i = + Δ
t t t
- Discrétisation des équations de transferts couplées
( g )
fi+ 1 X , T, P = 0
- Résolution du système par la méthode de Newton-Raphson
X ( x j ), T ( x j ), Pg ( x j ) , ( x j ) j:1→ Ncv
- Sauvegarde
Programme Principal : TransPore
- Changement des conditions aux limites
- Appel du Programme de Mécanique
- Test : Fin du pas de temps ?
Fin
nœuds du maillage
Figure III.11. Structure du code TransPore et du sous programme de mécanique.
(*) Le choix du pas Δ t est délicat car s’ il est grand, la convergence vers la solution du système sera difficile. En
revanche, s’il est petit, la résolution sera très lente. Par conséquent, TransPore a un pas de temps variable dont la
valeur s’adapte à la vitesse de convergence de la méthode de résolution du système d’équations.
¢¡¤£
- Discrétisation des équations traduisant l’équilibre mécanique
A⋅ X = B
t t
i+ 1 Final >
¥§¦©¨
Sous-Programme : Mécanique du séchage
- Résolution du système linéaire par la méthode des gradients Bi-Conjugués
Champ de déplacements et courbure ( u ', α, β ,1 ρ )
- Calcul des contraintes et des déformations liées à ce champ de déplacement.
- Sauvegarde et mémorisation des contraintes vues par le matériau à l’instant t i+
1
Δ
t