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 :math:`r/R`.
Les 3 courbes correspondent à trois valeurs différentes de :math:`\gamma = H/R` avec :math:`H` la hauteur d'eau et :math:`R` le rayon du réservoir.

.. image:: images/abaques_pi.png
    :width: 420
    :height: 536

En abscisses, il s'agit de la pression impulsive normalisée à :math:`\rho_l R ZPA_h`, avec :math:`\rho_l` la masse volumique de l'eau et :math:`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.

.. csv-table::

    ":math:`r/R`", ":math:`P_i^{adim} ,\gamma=3`", ":math:`P_i^{adim}, \gamma=1`", ":math:`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 à :math:`R=H/\gamma` pour les trois valeurs de :math:`\gamma` à tester.


Grandeurs testées et résultats
------------------------------

On teste la valeur de la pression normalisée pour les 11 rapports :math:`r/R` du tableau ci-dessus et pour les 3 valeurs de :math:`\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.

.. image:: images/figure_pir_norma.png
    :width: 948
    :height: 577


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.