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.

../../../../_images/abaques_pi.png

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.

../../../../_images/figure_pir_norma.png

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.

../../../../_images/sslv305c.png