commit:2017:08_10
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
commit:2017:08_10 [2017/08/10 15:53] – [Divers] wautelet | commit:2017:08_10 [2017/08/11 16:43] (current) – wautelet | ||
---|---|---|---|
Line 1: | Line 1: | ||
- | ===== Commit 2017-06-07 ====== | + | ===== Commit 2017-08-11 ====== |
Ce commit est pour améliorer quelques routines liées au contact et pour mettre progressivement mes développements sur la version courante. | Ce commit est pour améliorer quelques routines liées au contact et pour mettre progressivement mes développements sur la version courante. | ||
- | ===== Bicubic Coons ===== | + | ===== Méthode d' |
- | J' | + | Suite à mes derniers développements sur le Lagrangien Augmenté, j' |
- | Si nous avons une surface plane constituée de patch de Coons, nous avons d'office une continuité des tangentes au niveau des frontières d'un patch à un autre. Lors que la surface n'est plus plane, il faut certaines conditions pour que les tangentes soient continues, ce qui n'est toujours aisé à remplir ! | + | Pour pouvoir utiliser les nouveaux schémas |
- | + | ||
- | Voici les commandes | + | |
< | < | ||
- | surfaceset | + | parameters[' |
- | | + | augLagAugmentation = AugLagExtrapolationAugmentation(alm) |
+ | augLagAugmentation.setExtrapolationMethod(parameters[' | ||
</ | </ | ||
- | Dans le cas de la boite carré, voici les différences que nous observons pour une pénalité | + | Parmi ces schémas, celui qui semble le plus efficace à l' |
- | < | + | ===== Robustesse des projections ===== |
- | Version avec patch de Coons bicubique : | + | |
- | User CPU : 1h6m52.28s | + | === Critère de divergence des itérations === |
- | Real CPU : 19m52.39s | + | |
- | + | ||
- | [TSC-STP] | + | |
- | [TSC-ITE] | + | |
- | [TSC-INW] | + | |
- | [TSC-EXW] | + | |
- | [TSC-CPU] | + | |
- | [TSC-REA] | + | |
- | [TSC-KER] | + | |
- | [TSC-MEM] | + | |
- | [TSC-EXT] | + | |
- | [TSC-EXT] | + | |
- | [TSC-EXT] | + | |
- | [TSC-EXT] | + | |
- | [TSC-EXT] | + | |
- | [TSC-EXT] | + | |
- | </ | + | |
- | < | + | Afin d' |
- | Version avec patch de Coons : | + | |
+ | === Critère de convergence === | ||
+ | |||
+ | J'ai constaté que le terme de droite dans les critères normés sur les angles étaient trop stricts lorsque le résidu était faible ou les normes de tangentes étaient trop faibles. J'ai ajouté tout simplement une opération maximale entre la valeur 1.0 et la valeur actuelle de la norme pour passer automatiquement d'une norme relative à une norme absolue. | ||
+ | |||
+ | === Inexact Line Search === | ||
+ | |||
+ | J'ai constaté que l’exécution du Line Search inexacte renvoie une valeur du pas proche de zéro (ou même zéro) suite à des erreurs d' | ||
+ | |||
+ | === Interface des tolérances In/Out pour les outils de contact === | ||
+ | |||
+ | Après avoir lancé une série de test sur le cluster sur la squareBox et sur le sRail, j'ai été assez surpris de mes résultats sur la sensibilité aux coefficients de pénalité. En affichant les zones potentiellement en contact, j'ai constaté que des projections à priori valides n' | ||
- | User CPU : 1h10m25.63s | + | Clairement là, il y a un problème suite à l' |
- | Real CPU : 21m6.52s | + | |
- | + | ||
- | [TSC-STP] | + | |
- | [TSC-ITE] | + | |
- | [TSC-INW] | + | |
- | [TSC-EXW] | + | |
- | [TSC-CPU] | + | |
- | [TSC-REA] | + | |
- | [TSC-KER] | + | |
- | [TSC-MEM] | + | |
- | [TSC-EXT] | + | |
- | [TSC-EXT] | + | |
- | [TSC-EXT] | + | |
- | [TSC-EXT] | + | |
- | [TSC-EXT] | + | |
- | [TSC-EXT] | + | |
- | </ | + | |
- | Pour montrer la discontinuité de la normale aux interfaces entre les patchs de Coons, il suffit d' | + | Sans tolérance Out, on observe que des points ne sont pas en contact |
- | Version avec patch de Coons | + | {{: |
- | {{: | + | Avec tolérance Out de 1E-3, on observe que ces points sont en contact sur le coin de la boîte. |
- | Version avec patch de Coons bicubique | + | {{: |
+ | Il reste malheureusement à choisir la valeur de cette tolérance par cas-test ou pour tous les cas tests (1E-6?). Il faut savoir que les tests du type SurroundednessTest2D ou SurroundednessTest3D utilise une tolérance de 1E-6 pour le IN et OUT. | ||
- | {{: | + | < |
+ | skinsetPunchTool = ContactTool(skinsetPunch) | ||
+ | skinsetPunchTool.setOuterTolerance(parameters[' | ||
+ | skinsetPunchTool.setInnerTolerance(parameters[' | ||
+ | </ | ||
+ | === Projection Informations === | ||
- | ===== Correction Patch de Coons ===== | + | Finalement, j'ai ajouté la possibilité |
- | Par la même occasion, j'ai corrigé une erreur dans les dérivées des tangentes pour les patchs de Coons. Elles sont utilisées | + | < |
+ | skinsetPunchTool = ContactTool(skinsetPunch) | ||
+ | skinsetPunchTool.setProjectionInfoVerbose(True) | ||
+ | </ | ||
+ | |||
+ | Ce qui donne par exemple | ||
+ | |||
+ | < | ||
+ | Projection Infos: 2897/3924 effective projections and 25/3924 failed projections | ||
+ | (average nbNRIt = 3.59372; average nbLSIt = 2.74042; average effective nbLSIt = 2.6236) | ||
+ | </ | ||
===== Divers ===== | ===== Divers ===== | ||
=== Visualisation === | === Visualisation === | ||
+ | |||
+ | J'ai corrigé quelques bugs : | ||
+ | |||
+ | * Affichage des territoires de contact pour l'auto contact ! | ||
+ | * Ajout d'un test sur la présence d' | ||
=== Cas test === | === Cas test === | ||
Line 92: | Line 88: | ||
< | < | ||
- | [a]:mtGeo\mtGeoBicubicCoons.cpp | + | [a]:mtGeo/ |
- | [a]:mtGeo\mtGeoBicubicCoons.h | + | [a]:mtGeo/ |
+ | [a]: | ||
+ | [a]: | ||
+ | [a]: | ||
+ | [a]: | ||
[r]: | [r]: | ||
</ | </ | ||
Line 100: | Line 100: | ||
< | < | ||
- | [a]:apps\qs\contactBicubicCoons.py | + | [a]: |
[r]: | [r]: | ||
</ | </ | ||
- | --- // | + | --- // |
commit/2017/08_10.txt · Last modified: 2017/08/11 16:43 by wautelet