Metafor

ULiege - Aerospace & Mechanical Engineering

User Tools

Site Tools


commit:2015:01_20

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision
Previous revision
commit:2015:01_20 [2015/01/21 09:52] carrettacommit:2015:01_20 [2016/03/30 15:23] (current) – external edit 127.0.0.1
Line 7: Line 7:
 Pour y parvenir, je vais procéder en 2 étapes. Tout d'abord, permettre à l'utilisateur d'utiliser une épaisseur hors plan - en état plan de déformation - différente de 1. Cette première étape est l'objet de ce commit. La deuxième consistera à prendre en compte la variation de cette grandeur au cours de la simulation. Pour y parvenir, je vais procéder en 2 étapes. Tout d'abord, permettre à l'utilisateur d'utiliser une épaisseur hors plan - en état plan de déformation - différente de 1. Cette première étape est l'objet de ce commit. La deuxième consistera à prendre en compte la variation de cette grandeur au cours de la simulation.
  
-La modification principale a été apportée dans le calcul du Jacobien (cfr.**SfIntegrationPointTemplate<Dimension2D>::computeJacobian(...)**). La composante 33 du Jacobien ne vaut donc plus nécessairement 1 en état plan de déformations mais cette valeur correspond désormais à l'épaisseur hors plan (hzz) définie par l'utilisateur. Pour prendre en compte ce paramètre dans le jeu de données, on l'introduit en argument de la fonction suivante:+La modification principale a été apportée dans le calcul du Jacobien (cfr.**SfIntegrationPointTemplate<Dimension2D>::computeJacobian(...)**). La composante 33 du Jacobien ne vaut donc plus nécessairement 1 en état plan de déformation, mais cette valeur correspond désormais à l'épaisseur hors plan (hzz) définie par l'utilisateur. Pour prendre en compte ce paramètre dans le jeu de données, on l'introduit en argument de la fonction suivante:
  
 <code> <code>
Line 17: Line 17:
 ==== Modifications apportées pour le calcul du contact ==== ==== Modifications apportées pour le calcul du contact ====
  
-Pour le cas du contact faisant appel à AREAINCONTACT, le calcul de l'aire associée aux noeuds a été modifié (avec l'aide de Gaëtan): en pratique, on multiplie la valeur utilisée jusqu'ici par la moyenne des épaisseurs hors plan des éléments incident aux noeuds de l'interaction de contact.+Pour le cas du contact faisant appel à AREAINCONTACT, le calcul de l'aire associée aux noeuds a été modifié (avec l'aide de Gaëtan): en pratique, on multiplie la valeur utilisée jusqu'ici par la moyenne des épaisseurs hors plan des éléments incidents aux noeuds de l'interaction de contact.
  
 ==== Modifications apportées dans l'ALE ==== ==== Modifications apportées dans l'ALE ====
  
-Sur les conseils de Luc, j'ai utilisé une valeur par défaut hzz = 0.0 dans les fonctions faisant appel au calcul du Jacobien (et pas 1.0) de manière à faire planter la simulation au cas où j'aurais oublié de passer l'épaisseur hors plan en argument qque part dans le code. C'était une bonne idée puisque ça m'évit de faire une erreur dans le calcul des intégrales sur les cellules utilisées en ALE. L'épaisseur était jusqu'ici supposée unitaire en Etat plan de déformation, on considère que c'est toujours le cas en passant explicitement une épaisseur de 1.0 dans les fonctions. Comme l'épaisseur hors plan ne varie pas au cours de la phase de transfert de données, le fait d'utiliser une épaisseur unitaire n'a pas d'incidence sur le résultat.+Sur les conseils de Luc, j'ai utilisé une valeur par défaut hzz = 0.0 dans les fonctions faisant appel au calcul du Jacobien (et pas 1.0) de manière à faire planter la simulation au cas où j'aurais oublié de passer l'épaisseur hors plan en argument qque part dans le code. C'était une bonne idée puisque ça m'évite de faire une erreur dans le calcul des intégrales sur les cellules utilisées en ALE. L'épaisseur était jusqu'ici supposée unitaire en Etat plan de déformation, on considère que c'est toujours le cas en passant explicitement une épaisseur de 1.0 dans les fonctions. Comme l'épaisseur hors plan ne varie pas au cours de la phase de transfert de données, le fait d'utiliser une épaisseur unitaire n'a pas d'incidence sur le résultat.
  
 Cette modification a été apportée dans la classe **mtGeoSide** ainsi que dans **SfIntegratorTemplate** (cfr. codes ci-dessous) Cette modification a été apportée dans la classe **mtGeoSide** ainsi que dans **SfIntegratorTemplate** (cfr. codes ci-dessous)
Line 75: Line 75:
 J'ai ajouté un répertoire appelé /gps/ (generalized plane strain) dans /apps/. Dans ce répertoire, j'ai ajouté les cas tests qui vérifient la validité des développements réalisés. J'ai ajouté un répertoire appelé /gps/ (generalized plane strain) dans /apps/. Dans ce répertoire, j'ai ajouté les cas tests qui vérifient la validité des développements réalisés.
  
-La premère série de cas tests consiste à exraser un bloc carré de 1 mm de coté. Trois types de chargement différents sur le bord supérieur sont testés : +La première série de cas tests consiste à éraser un bloc carré de 1 mm de coté. Trois types de chargement différents sur le bord supérieur sont testés : 
  
   - déplacement imposé   - déplacement imposé
commit/2015/01_20.1421830349.txt.gz · Last modified: 2016/03/30 15:22 (external edit)

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki