===== Commit - LPx - 12 Novembre 2019 ===== AEJ : bug fix calcul du volume des triangles Axisym ===== AEJ : ===== ==== Calcul des volumes des éléments triangulaires en Axisym : ==== * Lorsque Yves a implémenté l'état plan généralisé (2.5D), il a malheureusement oublié de modifier 2 appels à la fonction 'computeVolume' (ce qui marchait vu les valeurs par défaut définie dans la définition de la fonction) SfIntegratorTemplate : double computeVolume(const std::vector Co, double thick = 1.0, bool isAxisym=false); * Le calcul du volume axisymétrique des triangleSide et QuadraticTriangleSide était donc faux (ce qui pose de gros problème pour les AEJ en axisymétrique), mais vu l'impossibilité de faire des monos-éléments AEJ (il faut des voisins), c'était passé sous le radar... * Le comportement des AEJ en hypo semblent déjà plus proches de celui en Etat plan de défo (c'est encore loin d'être parfait) === STEPS :=== ^ Test ^ Old ^ New ^ absDiff ^ relDiff (%) ^ |apps.aej.taylor2D - Number of steps [ .py][ .res] | 3753| 3617| -136| -3.6| |apps.aej.taylor2DChAna - Number of steps [ .py][ .res] | 1936| 841| -1095| -56.6| |apps.aej.taylor2DChNum - Number of steps [ .py][ .res] | 119| 123| 4| 3.4| |apps.aej.taylor2DHypo - Number of steps [ .py][ .res] | 5489| 11281| 5792| 105.5| |apps.aej.taylor2DHypoCHAna - Number of steps [ .py][ .res] | 1622| 206| -1416| -87.3| |apps.aej.taylor2DHypoCHNum - Number of steps [ .py][ .res] | 125| 111| -14| -11.2| === ITERATIONS : === ^ Test ^ Old ^ New ^ absDiff ^ relDiff (%) ^ |apps.aej.taylor2DChAna - Number of mech. iterations [ .py][ .res] | 6954| 3116| -3838| -55.2| |apps.aej.taylor2DChNum - Number of mech. iterations [ .py][ .res] | 359| 370| 11| 3.1| |apps.aej.taylor2DHypoCHAna - Number of mech. iterations [ .py][ .res] | 5627| 626| -5001| -88.9| |apps.aej.taylor2DHypoCHNum - Number of mech. iterations [ .py][ .res] | 385| 336| -49| -12.7| ===== Fichiers ajoutés/supprimés : ===== Added : Deleted : Moved : ===== Tests ajoutés/supprimés ===== added : deleted : Moved : --- //[[L.Papeleux@ulg.ac.be|Luc Papeleux]] 2019/11/12 //