v5.01.109 SDND109 - Loi de comportement CHOC_ENDO, en dynamique non-linéaire#

Résumé:

Ce test valide en dynamique non-linéaire la loi de comportement CHOC_ENDO dédiée aux discrets.

Solutions de référence#

Modélisation A#

Chargements#

Le chargement est un déplacement imposé au nœud \(\mathrm{PT}3\) .

\(\mathit{Depl}(t)=0.25\ast \sin(2\pi t)+0.10\ast \sin(2\pi 1.5t)+0.10\ast \sin(2\pi 3t)\)

Le nœud \(\mathrm{PT}2\) est bloqué.

Résultats de référence#

Le déplacement est imposé par le chargement. La vérification se fait en comparant la réponse du discret par rapport aux courbes données dans les paramètres matériaux.

Incertitude sur la solution#

Modélisation B#

Chargements#

Le chargement est un déplacement imposé au nœud \(\mathrm{PT}1\) , c’est le même que celui de la modélisation \(A\) . Le nœud \(\mathrm{PT}2\) est bloqué.

Résultats de référence#

Le déplacement est imposé par le chargement. La vérification se fait en comparant la réponse du discret par rapport aux courbes données dans les paramètres matériaux.

Incertitude sur la solution#

Comparaison graphique entre la réponse et les données matériaux.

Modélisation C#

Chargements#

Le chargement est un déplacement imposé au nœud \(\mathrm{PT}1\) , c’est le même que celui desmodélisations \(A\) et \(B\) . Le nœud \(\mathrm{PT}2\) est bloqué.

Résultats de référence#

Le déplacement est imposé par le chargement. La vérification se fait en comparant la réponse du discret par rapport aux courbes données dans les paramètres matériaux.

Incertitude sur la solution#

Comparaison graphique entre la réponse et les données matériaux.

Modélisation D#

Chargements#

Le chargement est un déplacement imposé au nœud \(\mathrm{PT}1\) , c’est le même que celui de la modélisation \(A\) , \(B\) et \(C\) . Le nœud \(\mathrm{PT}2\) est bloqué. La seule différence est dans la définition du matériau où CRIT_AMOR = ’EXCLUS’.

Résultats de référence#

Le déplacement est imposé par le chargement. La vérification se fait en comparant la réponse du discret par rapport aux courbes données dans les paramètres matériaux.

Incertitude sur la solution#

Comparaison graphique entre la réponse et les données matériaux.

Modélisation A#

Caractéristiques de la modélisation#

Pour le 1er calcul le pas de temps est de \(0.005\mathrm{s}\) (■). Pour le 2nd calcul le pas de temps est de \(0.010\mathrm{s}\) (■).

Grandeurs testées et résultats#

Les grandeurs testées sont le déplacement venant du chargement, l’effort existant dans le discret ainsi que la variable de plasticité cumulée.

La figure ci-dessous présente les courbes issues des 2 simulations. À ces résultats sont superposés les courbes issues des paramètres matériaux. Les décharges en pointillés ( - - - - ) correspondent aux raideurs définies dans les paramètres du matériau.

Les points (●) correspondent aux tests réalisés dans le fichier de commandes.

../../../../_images/1000020100000DBC000008F81EC175354FF64115.png

Modélisation B#

Caractéristiques de la modélisation#

Pour le 1er calcul le pas de temps est de \(0.005\mathrm{s}\) (■). Pour le 2nd calcul le pas de temps est de \(0.010\mathrm{s}\) (■).

Grandeurs testées et résultats#

Les grandeurs testées sont le déplacement venant du chargement, l’effort existant dans le discret ainsi que la variable de plasticité cumulée.

La figure ci-dessous présente les courbes issues des 2 simulations. À ces résultats sont superposés les courbes issues des paramètres matériaux. Les décharges en pointillés ( - - - - ) correspondent aux raideurs définies dans les paramètres du matériau.

Les points (●) correspondent aux tests réalisés dans le fichier de commandes.

../../../../_images/1000020100000DBC000008F8A1FACA8EF78A4AE8.png

Modélisation C#

Caractéristiques de la modélisation#

Grandeurs testées et résultats#

Les grandeurs testées sont le déplacement venant du chargement, l’effort existant dans le discret ainsi que la variable de plasticité cumulée.

La figure ci-dessous présente les courbes issues des 2 simulations. À ces résultats sont superposés les courbes issues des paramètres matériaux. Les décharges en pointillés ( - - - - ) correspondent aux raideurs définies dans les paramètres du matériau.

Les points (●) correspondent aux tests réalisés dans le fichier de commandes.

../../../../_images/1000020100000F42000009F7D304B62B93606BB1.png

Le décalage entre l’effort et la valeur théorique vient de la prise en compte de l’amortissement qui “réduit” le seuil. L’effet de l’amortissement est pris en compte dans le calcul du seuil. Dans la définition du matériau DIS_CHOC_ENDO: CRIT_AMOR = ’INCLUS’.

Ce décalage n’a aucune raison de rester constant au cours du calcul, il vaut \(\text{Amor\_Nor}(p).V\) .

Modélisation D#

Caractéristiques de la modélisation#

Grandeurs testées et résultats#

Les grandeurs testées sont le déplacement venant du chargement, l’effort existant dans le discret ainsi que la variable de plasticité cumulée.

La figure ci-dessous présente les courbes issues des 2 simulations. À ces résultats sont superposés les courbes issues des paramètres matériaux. Les décharges en pointillés ( - - - - ) correspondent aux raideurs définies dans les paramètres du matériau.

Les points (●) correspondent aux tests réalisés dans le fichier de commandes.

../../../../_images/1000020100000F48000009FC39BB9F211D87CD0F.png

Par rapport à la modélisation C, il n’existe plus de décalage entre l’effort et le seuil.L’effet de l’amortissement n’est plus pris en compte dans le calcul du seuil. Dans la définition du matériau DIS_CHOC_ENDO: CRIT_AMOR = ’EXCLUS’.

Synthèse des résultats#

On trouve les résultats attendus.