Les trois variables d'états sont T pour la température,  M pour l'humidité et U pour les déplacements. D'abord les 3 équations régissant le phénomène:

\[  C_T \frac{\partial T}{\partial t} - \nabla \cdot ( K_T \nabla T) = 0  \]      (1)
\[  C_M \frac{\partial M}{\partial t} - \nabla \cdot ( K_M \nabla M + K_{TM} \nabla T) = 0  \]      (2)
\[  \rho \frac{\partial^2 U}{\partial t^2} - \nabla \cdot (E(\epsilon - \beta_M \Delta M - \beta_T \Delta T)) = 0  \]      (3)




\[  C_T = C_T(\rho,T,M)  \]      \[  C_M = C_M(\rho,T,M)  \]

\[  K_T = K_T(\rho,T,M)  \]      \[  K_M = K_M(\rho,T,M)  \]      \[  K_{TM} = K_{TM}(\rho,T,M)  \]

\[  E = E(\rho,T,M)  \]     \[  \beta_T = \beta_T(\rho,T,M)  \]      \[  \beta_M = \beta_M(\rho,T,M)  \]


Pour la première partie du processus de fabrication on applique une pression et de la chaleur à la plaque de MDF pour un court laps de temps (20 secondes). Si on note  \[  \Gamma_h  \] et \[  \Gamma_b  \] les faces supérieures et inférieure (respectivement) du panneau,\[  \Gamma_{sx}  \] et \[  \Gamma_{sy}  \] les rives sur lesquels on applique la symétrie des variables d'état et \[  \Gamma_0  \] le point (0,0,h/2) pour lequel on fixe les déplacements verticaux (exclus les mouvements de corps rigide) on a:

conditions initiales:  \[  T = T_0  \]      \[  M = M_0  \]      \[  U = 0  \]      \[  \sigma = 0  \]      \[  \sigma = E(\epsilon - \beta_M \Delta M - \beta_T \Delta T)  \]

conditions aux bords: \[  \sigma_n = P  \]      \[  T = T_P  \]    sur \[  \Gamma_h  \] et \[  \Gamma_b  \]      UZ = 0 sur \[  \Gamma_0  \]       \[  U_x = 0  \]      sur \[  \Gamma_{sx}  \]      \[  U_y = 0  \]     sur \[  \Gamma_{sy}  \]        


Cas 1



Temps initial t0
0
Temps final tf
20 sec
Température sur les faces sup. et inf. TP
150 C
Température initiale et ambiante T0 = Tinfini
20 C
Humidité initiale M0
6.5
Humidité ambiante Minfini
12.5
Déplacements initiaux U0 (0,0,0)
Pression sur les faces supérieure et inférieures P
4.5e+6 Pa
densité
540
B1 = B2
absorption 
1.7e-4
dessorption
3.4e-4
B3
absorption 
7.4e-3
dessorption 5.4e-3
E1 = E2
7.4e8
M <= 6.6
7.4e8 - 4.58333e7(M-6.6)
6.6 < M < 9.0
6.3e8
M >= 9.0
E3
1.7e7
M <= 6.6
1.7e7 - 2.08333e6(M-6.6)
6.6 < M < 9.0
1.2e7
M >= 9.0
G1
2.8e8
M <= 6.6
2.8e8 - 1.666e7(M-6.6)  
6.6 < M < 9.0
2.4e8
M >= 9.0
G2 = G3
8.4e7
M <= 6.6
8.4e7 - 5.41666e6(M-6.6)
6.6 < M < 9.0
7.12e7
M >= 9.0
Nu1
0.3
Nu3
0.2
D1 Coefficient de diffusion
1e-9
M <= 6.6
1e-9 - 3.7541e-10(M-6.6)
6.6 < M < 9.0
9.9e-11
M >= 9.0
KM
0
0
0
0
0
0
0
0
D1
KTM
1000D1
0
0
0
1000D1
0
0
0
1000D1
KT
0.185
0
0
0
0.185
0
0
0
0.185
BM
B1
0
0
0
B1
0
0
0
B3
BT
0 (tenseur)
CM
0.01*densité
CT
densité