Metafor

ULiege - Aerospace & Mechanical Engineering

User Tools

Site Tools

You are here: metafor.ltas.ulg.ac.be » Commits » 2014 » Commit 2014-03-19
Trace: • Commit 2015-11-11 • Commit 2017-09-05 • Commit 2015-10-10 • Commit 2015-12-11 • sidebar
commit:2014:03_19

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:2014:03_19 [2014/03/18 19:24] – [Cas Test] wauteletcommit:2014:03_19 [2016/03/30 15:23] (current) – external edit 127.0.0.1
Line 5: Line 5:
 ==== Ajout des matériaux de contact thermomécanique collant et avec frottement selon Tresca ==== ==== Ajout des matériaux de contact thermomécanique collant et avec frottement selon Tresca ====
  
-J'ai ajouté le matériau thermomécanique collant et avec frottement selon Tresca. +J'ai ajouté le matériau thermomécanique collant et avec frottement selon Tresca. Par la même occasion, j'ai documenté un peu plus tous les matériaux thermomécaniques de contact :   
 + 
 + 
 +=== TmFrictionlessContactMaterial === 
 + 
 +== Description == 
 + 
 +Loi de contact thermomécanique sans frottement.  
 + 
 +Le flux de chaleur $q_{N}$ normal à l'interface de contact (sortant de la surface esclave) est déterminée par  
 + 
 +$$q_{N} = h_c \left(p_{N} \right) \left(T^{S} - T^{M}\left(\bf{\xi}^{S}\right)\right), $$ 
 + 
 +où 
 +  * $p_{N}$ est la pression de contact,  
 +  * $T^{S}$ est la température du nœud esclave,  
 +  * $T^{M}\left(\bf{\xi}^{S}\right)$ est la température du point sur la surface maître qui correspond à la projection la plus proche du nœud esclave sur la surface maître, 
 +  * $h_c$ est le coefficient de résistance thermique de conduction. 
 + 
 +Le coefficient de résistance thermique de conduction $h_c$ est modélisée par  
 + 
 +$$h_c \left(p_{N} \right) = h_{c0} \left(\frac{p_{N}}{H_v}\right)^{w}, $$       
 + 
 +où  
 +     * $H_v$ la dureté du matériau selon Vickers, 
 +     * $w$ est un exposant,  
 +     * $h_{c0}$ le coefficient de résistance thermique de conduction nominale. 
 + 
 +== Paramètres == 
 + 
 +^                Nom                   ^     Code Metafor     ^  Type de dépendance  ^ 
 +| Pénalité normale                     |   ''PEN_NORMALE''    |        ''TM''        | 
 +| Profondeur de détection du contact   |    ''PROF_CONT''     |          -           | 
 +| Coefficient d'échange nominal        |  ''CTM_H_NOMINAL''   |          -           | 
 +| Exposant de la loi d'échange         |  ''CTM_EXPONENT_E''  |          -           | 
 +| Dureté du matériau                   |   ''CTM_HARDNESS''   |          -           | 
 +| Type de contact                      |    ''TYPE_CONT''     |          -           | 
 + 
 +:!: Cette loi n'est pas testée en 3D :!: 
 + 
 +=== TmStickingContactMaterial === 
 + 
 +== Description == 
 + 
 +Loi de contact thermomécanique collant.  
 + 
 +== Paramètres == 
 + 
 +^                Nom                   ^     Code Metafor     ^  Type de dépendance  ^ 
 +| Pénalité normale                     |   ''PEN_NORMALE''    |        ''TM''        | 
 +| Pénalité tangente                    |   ''PEN_TANGENT''    |        ''TM''        | 
 +| Profondeur de détection du contact   |    ''PROF_CONT''     |          -           | 
 +| Coefficient d'échange nominal        |  ''CTM_H_NOMINAL''   |          -           | 
 +| Exposant de la loi d'échange         |  ''CTM_EXPONENT_E''  |          -           | 
 +| Dureté du matériau                   |   ''CTM_HARDNESS''   |          -           | 
 +| Type de contact                      |    ''TYPE_CONT''     |          -           | 
 + 
 +:!: Cette loi n'est pas testée en 3D :!: 
 + 
 +=== TmCoulombContactMaterial === 
 + 
 +== Description == 
 + 
 +Loi de contact thermomécanique avec frottement de Coulomb. 
 + 
 +== Paramètres == 
 + 
 +^                Nom                   ^     Code Metafor     ^  Type de dépendance  ^ 
 +| Pénalité normale                     |   ''PEN_NORMALE''    |        ''TM''        | 
 +| Pénalité tangente                    |   ''PEN_TANGENT''    |        ''TM''        | 
 +| Profondeur de détection du contact   |    ''PROF_CONT''     |          -           | 
 +| Coefficient de frottement statique   |  ''COEF_FROT_STA''   |        ''TM''        | 
 +| Coefficient de frottement dynamique  |  ''COEF_FROT_DYN''   |        ''TM''        | 
 +| Coefficient d'échange nominal        |  ''CTM_H_NOMINAL''   |          -           | 
 +| Exposant de la loi d'échange         |  ''CTM_EXPONENT_E''  |          -           | 
 +| Dureté du matériau                   |   ''CTM_HARDNESS''   |          -           | 
 +| Type de contact                      |    ''TYPE_CONT''     |          -           | 
 + 
 +:!: Cette loi n'est pas testée en 3D :!: 
 + 
 +=== TmTrescaContactMaterial === 
 + 
 +== Description == 
 + 
 +Loi de contact thermomécanique avec frottement de Tresca. 
 + 
 +== Paramètres == 
 + 
 +^                Nom                   ^     Code Metafor     ^  Type de dépendance  ^ 
 +| Pénalité normale                     |   ''PEN_NORMALE''    |        ''TM''        | 
 +| Pénalité tangente                    |   ''PEN_TANGENT''    |        ''TM''        | 
 +| Profondeur de détection du contact   |    ''PROF_CONT''     |          -           | 
 +| Coefficient de frottement statique   |  ''COEF_FROT_STA''   |        ''TM''        | 
 +| Coefficient de frottement dynamique  |  ''COEF_FROT_DYN''   |        ''TM''        | 
 +| Coefficient d'échange nominal        |  ''CTM_H_NOMINAL''   |          -           | 
 +| Exposant de la loi d'échange         |  ''CTM_EXPONENT_E''  |          -           | 
 +| Dureté du matériau                   |   ''CTM_HARDNESS''   |          -           | 
 +| Type de contact                      |    ''TYPE_CONT''     |          -           | 
 + 
 +Le seuil est généralement calculé par $m\,\sigma_0\,/\sqrt{3}$ où //m// est le coefficient de frottement de Tresca et $\sigma_0$ est la limite élastique en traction du matériau. 
 + 
 +:!: Cette loi n'est pas testée en 3D :!:
  
 ==== Bug - Éléments de contact thermomécanique ==== ==== Bug - Éléments de contact thermomécanique ====
  
-J'ai corrigé un bug +L'aire de contact d'échange thermique n'est pas la même que l'aire de contact associée à chaque nœud esclave dans le cas de la méthode area in contact (cf. Thèse de Laurent Adam). Donc j'ai dû adapter le calcul de la pression de contact intervenant dans le coefficient de résistance thermique de conduction (voir début du commit). De fait, on peut utiliser les éléments de contact thermomécanique sans area in contact, dès lors on divise la force de contact par l'aire de contact d'échange thermique pour obtenir une pression de contact "équivalente". 
  
 ==== Divers - Gestion des itérations mécaniques/thermiques ==== ==== Divers - Gestion des itérations mécaniques/thermiques ====
Line 30: Line 131:
  
 Il est intéressant de désactiver ce critère d'arrêt, lorsque l'on désire comparer des méthodes de calcul sur des cas tests académiques.  Il est intéressant de désactiver ce critère d'arrêt, lorsque l'on désire comparer des méthodes de calcul sur des cas tests académiques. 
- + 
 +{{ :commit:2014:convergence.png?400 |}} 
  
 ==== Cas Test ====  ==== Cas Test ==== 
Line 41: Line 143:
 <code> <code>
 [r]: [r]:
-[a]:+[a]:mtMaterials/boundaries/TmStickingContactMaterial.h 
 +[a]:mtMaterials/boundaries/TmStickingContactMaterial.cpp 
 +[a]:mtMaterials/boundaries/TmTrescaContactMaterial.h 
 +[a]:mtMaterials/boundaries/TmTrescaContactMaterial.cpp
 </code> </code>
  
Line 48: Line 153:
 <code> <code>
 [r]: [r]:
-[a]:apps/tutorials/cylinder.py +[a]:apps/contactTm/RdEplBlockFrictionLess.py  
 +[a]:apps/contactTm/RdEplBlockSticking.py  
 +[a]:apps/contactTm/RdEplBlockCoulomb.py  
 +[a]:apps/contactTm/RdEplBlockTresca.py 
 </code> </code>
  
  
  --- //[[gwautelet@ulg.ac.be|Gaëtan WAUTELET]] 2014/03/18//  --- //[[gwautelet@ulg.ac.be|Gaëtan WAUTELET]] 2014/03/18//
commit/2014/03_19.1395167057.txt.gz · Last modified: (external edit)
Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki