v3.04.305 SSLV305 - Calcul de la pression d’un réservoir sous séisme suivant l’Eurocode 8#
Résumé:
Ce cas test est destiné à valider la macro-commande DEFI_PRES_EC8.
Modélisation B#
But de la modélisation#
Cette modélisation a pour but de valider le calcul de la pression impulsive rigide sur le fond du réservoir, c’est à dire pour des rayons plus petits que le rayon du réservoir.
Résultats de référence#
Pour valider cela, nous disponons d’abaques issus de l’Eurocode 8 qui donnent trois courbes de pression en fonction du rapport \(r/R\). Les 3 courbes correspondent à trois valeurs différentes de \(\gamma = H/R\) avec \(H\) la hauteur d’eau et \(R\) le rayon du réservoir.
En abscisses, il s’agit de la pression impulsive normalisée à \(\rho_l R ZPA_h\), avec \(\rho_l\) la masse volumique de l’eau et \(ZPA_h\) l’accélération impulsive rigide horizontale (opérande ACCE_SP_H). Ces courbes ont été numérisées pour fournir des valeurs de références.
\(r/R\) |
\(P_i^{adim} ,\gamma=3\) |
\(P_i^{adim}, \gamma=1\) |
\(P_i^{adim}, \gamma=0.5\) |
0 |
0 |
0 |
0 |
0.1 |
0.10176 |
0.06342 |
0.01769 |
0.2 |
0.19911 |
0.12684 |
0.03687 |
0.3 |
0.29941 |
0.19321 |
0.05899 |
0.4 |
0.3997 |
0.25958 |
0.08554 |
0.5 |
0.5 |
0.33038 |
0.11504 |
0.6 |
0.60029 |
0.40265 |
0.15339 |
0.7 |
0.70058 |
0.48082 |
0.20058 |
0.8 |
0.80235 |
0.56342 |
0.25958 |
0.9 |
0.90412 |
0.65191 |
0.33038 |
1 |
1 |
0.75073 |
0.42477 |
Paramètres#
On conserve les paramètres de la modélisation A exepté le rayon que l’on fixe à \(R=H/\gamma\) pour les trois valeurs de \(\gamma\) à tester.
Grandeurs testées et résultats#
On teste la valeur de la pression normalisée pour les 11 rapports \(r/R\) du tableau ci-dessus et pour les 3 valeurs de \(\gamma\). Les valeurs n’ayant pas d’intérêt, on préfère donner la figure comparant les résultats avec les références.
Commentaires#
Les résultats sont très proches de la référence ce qui valide le calcul de la pression impulsive rigide sur le fond du réservoir.
Modélisation C#
But de la modélisation#
Cette modélisation a pour but de valider l’utilisation de la formule issue de DEFI_PRES_EC8 avec FORCE_COQUE_FO d’AFFE_CHAR_MECA.
Paramètres#
On conserve les paramètres de la modélisation A.
Détails#
Le maillage est un cylindre avec fond (et plafond) modélisé en éléments DKT. Le fond du réservoir est encastré. On applique la pression issue de DEFI_FONC_EC8. Le calcul est élastique.
Grandeurs testées et résultats#
Après s’être assuré que la déformée avait globalement la forme attendue, des tests de non-régression ont été ajoutés sur deux points au milieu de la hauteur d’eau, le premier d’azimut 0 et le second d’azimut 180°. Les valeurs n’ont pas d’importance ici car ce sont des tests de non-régression.
L’image ci-dessous montre la déformée du réservoir.