Simple shear test

[QuentinJouet]

VERONIQUE :

Mise en place d'une expérience de type cisaillement simple pour explorer la rhéologie de la banquise, en régime granulaire sec.
Échantillon: assemblage de floes à haute densité, de forme carrée (x,y)
Conditions limites: un bord inférieur fixe, un bord supérieur à vitesse imposée en x et pression (contrainte normale) imposée en y (pas de forçage surfacique!). Conditions périodiques aux bords latéraux.
*Il faudrait donc pouvoir imposer la vitesse d'un obstacle qui sert de bord ou de ''mur'' supérieur et aussi contrôler la contrainte normale (sigma_yy) appliquée sur ce mur.

[QuentinJouet]

@veronique.dansereau je ne suis pas sûr de comprendre "pression (contrainte normale) imposée en y (pas de forçage surfacique!)", y-a-t'il une vitesse en y ? ou la possibilité d'un déplacement en y du bord supérieur ?@veronique.dansereau je ne suis pas sûr de comprendre "pression (contrainte normale) imposée en y (pas de forçage surfacique!)", y-a-t'il une vitesse en y ? ou la possibilité d'un déplacement en y du bord supérieur ?

[QuentinJouet]

SILOUANE :

@veronique.dansereau j'ai fait un petit script pour générer un input et lancer un run avec comme entrée le nombre et la taille des grains. Les floes sont poussés par le courant vers en bas à droite et sont entourés par deux "bumpy obstacles". 1_shearTest_100f_r25

On peut envisager de rajouter dans floedyn un mode de forcage --fmodes 7 pour contraindre le déplacement d'un floe (celui du haut ici donc), et soit ajouter une liste de paramètres soit hard-coder des infos plus fines dans physical_data.hpp (on pourrait par exemple imposer directement une composante de vitesse, une force, bloquer la rotation, définir sur quels floes on applique ces contraintes et dans quel intervalle de temps, etc). Qu'en penses-tu ?

En revanche au niveau du post-traitement, quelles sorties te faudrait-il ? J'en ai discuté un peu avec Aurélien, si j'ai bien compris vous calculez des contraintes de cisaillement (dans floedyn les floes ne se déforment pas, donc on ne peut pas trop calculer des contraintes dessus...)@veronique.dansereau j'ai fait un petit script pour générer un input et lancer un run avec comme entrée le nombre et la taille des grains. Les floes sont poussés par le courant vers en bas à droite et sont entourés par deux "bumpy obstacles".

1_shearTest_100f_r25.mp4

On peut envisager de rajouter dans floedyn un mode de forcage --fmodes 7 pour contraindre le déplacement d'un floe (celui du haut ici donc), et soit ajouter une liste de paramètres soit hard-coder des infos plus fines dans physical_data.hpp (on pourrait par exemple imposer directement une composante de vitesse, une force, bloquer la rotation, définir sur quels floes on applique ces contraintes et dans quel intervalle de temps, etc). Qu'en penses-tu ?
En revanche au niveau du post-traitement, quelles sorties te faudrait-il ? J'en ai discuté un peu avec Aurélien, si j'ai bien compris vous calculez des contraintes de cisaillement (dans floedyn les floes ne se déforment pas, donc on ne peut pas trop calculer des contraintes dessus...)