commit:2015:07_13
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revisionLast revisionBoth sides next revision | ||
commit:2015:07_13 [2015/07/13 17:20] – wautelet | commit:2015:07_13 [2015/07/13 18:47] – wautelet | ||
---|---|---|---|
Line 5: | Line 5: | ||
==== Méthode de mise à jour du Lagrangien augmenté tangentiel ==== | ==== Méthode de mise à jour du Lagrangien augmenté tangentiel ==== | ||
- | $$\boldsymbol{\lambda^{(k, | + | La force tangentielle de contact est définie comme : |
- | J'ai ajouté une option pour calculer une longueur caractéristique du contact dans le critère géométrique normalisé. | + | $$\boldsymbol{f}_t^{(k, |
+ | où | ||
+ | * $k$ est l'indice de l' | ||
+ | * $i$ est l' | ||
+ | * $\boldsymbol{f}_t$ est la force tangentielle. | ||
+ | * $\boldsymbol{\lambda}_t$ est le vecteur lagrangien augmenté tangentiel. | ||
+ | * $\boldsymbol{g}_t$ est le vecteur gap tangentiel. | ||
+ | * $C_t$ est le coefficient de pénalité tangentiel. | ||
- | Par défaut | + | Puisque le plan de normal $\boldsymbol{n}^{(k)}$ où se trouve |
+ | $$\boldsymbol{\lambda}_t^{(k,i)} = \boldsymbol{T} ( \boldsymbol{n}^{(i)}, | ||
+ | où | ||
+ | * $k$ est l' | ||
+ | * $i$ est l' | ||
+ | * $\boldsymbol{n}$ est la normale à la surface maître évaluée au point correspondant | ||
+ | * $\boldsymbol{T}$ est l' | ||
+ | * $\boldsymbol{\lambda}_t$ est le vecteur lagrangien augmenté tangentiel. | ||
- | Malgré sa simplicité, cette approche peut poser un problème du point de vue du contact normal, où c'est l'interpénétration des corps en contact qui compte et dès lors, il faut définir une distance caractéristique normale à la surface de contact. Par exemple, | + | L' |
+ | |||
+ | - soit une projection : $$\boldsymbol{T} ( \boldsymbol{n}^{(i)}, \boldsymbol{n}^{(k)} )= \boldsymbol{I} - \boldsymbol{n}^{(i)} \otimes \boldsymbol{n}^{(i)} | ||
+ | - soit une rotation (Formule | ||
+ | | ||
+ | * $\cos{\phi} = \boldsymbol{n}^{(k)} \times \boldsymbol{n}^{(i)}$ | ||
+ | * $\sin{\phi} = ||\boldsymbol{n}^{(k)} | ||
+ | * $\boldsymbol{n} = \frac{\boldsymbol{n}^{(k)} | ||
+ | * $\tilde{ }$ est l'opérateur du produit vectoriel sous forme matricielle. | ||
+ | |||
+ | On constate que dans le cas d' | ||
+ | |||
+ | Dans Metafor, nous pouvons faire le choix entre ces deux méthodes (pour les matériaux de contact avec une composante tangentielle) : | ||
< | < | ||
- | augLagCriterion = AugLagNormalisedGeoCriterion(alm) | + | materset.define(1, AugLagStickingContactMaterial) |
- | augLagCriterion.setCharacteristicLengthMeth(ALM_GEO_LENGTH) | + | materset(1).put(TANGENTAUGLAGUPDATEMETHOD, |
</ | </ | ||
- | Elle est implémentée pour le contact 2D (état plan déformation et modélisation axisymétrique) et le contact en 3D. | + | OPTION = TALUM_ROTATION |
- | Elle a été testé dans le cas de la rupture des éléments. | + | |
+ | Par défaut, nous avons toujours la méthode TALUM_PROJECTION, | ||
+ | |||
+ | Finalement, la linéarisation cohérente | ||
- | <note important> | + | ==== Divers - Tutorials ==== |
+ | J'ai ajouté le jeux de données des projets de seconde session de cette année. | ||
===== Fichiers ajoutés/ | ===== Fichiers ajoutés/ | ||
< | < | ||
- | [a]:AugLagCharacteristicLengthMethod.h | + | [a]:mtContact/ |
- | [a]:AugLagCharacteristicLengthMethod.cpp | + | [a]:mtContact/ |
[r]: | [r]: | ||
</ | </ | ||
Line 36: | Line 66: | ||
< | < | ||
[r]: | [r]: | ||
- | [a]: | + | [a]:apps/ |
</ | </ | ||
- | --- // | + | --- // |
commit/2015/07_13.txt · Last modified: 2016/03/30 15:23 by 127.0.0.1