Differences

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

--- commit:2017:04_22 [2017/04/25 22:52] – [Lagrangien Augmenté] wautelet
+++ commit:2017:04_22 [2017/04/26 14:50] (current) – wautelet
@@ Line 1: / Line 1: @@
-===== Commit 2017-04-22 ======
+===== Commit 2017-04-26 ======
 Ce commit est pour améliorer quelques routines liées au contact et pour mettre progressivement mes développements sur la version courante.
@@ Line 17: / Line 17: @@
 === Méthode interfacée ===
-Après un petit refactoring, pour ne pas changer, j'ai scindé en deux les méthodes de Newton Raphson Modifiée : CpuModifiedNRStiffnessUpdateCriterion et PeriodicModifiedNRStiffnessUpdateCriterion. Évidemment, la méthode par défaut est NRStiffnessUpdateCriterion. Toutes ces méthodes dérivent de StiffnessUpdateCriterionBase et ont pour but de décider si nous devons mettre à jour la matrice de raideur tangente lors du calcul de la correction des positions.
+Après un petit refactoring, pour ne pas changer, j'ai scindé en deux les méthodes de Newton Raphson Modifiée : CpuModifiedNRStiffnessUpdateCriterion et PeriodicModifiedNRStiffnessUpdateCriterion. Par ailleur, la méthode par défaut pour la stratégie de mise à jour de la matrice de raideur tangente est la méthode de Newton Raphson NRStiffnessUpdateCriterion. Toutes ces méthodes dérivent de StiffnessUpdateCriterionBase et ont pour but de décider si nous devons mettre à jour la matrice de raideur tangente lors du calcul de la correction des positions.
 === Stratégie de la méthode de Newton Raphson Modifiée  ===
@@ Line 30: / Line 30: @@
 Voici les méthodes disponibles et les options associées.
+== NRStiffnessUpdateCriterion (Méthode Par défaut)==
+<code>
+    mim = metafor.getMechanicalIterationManager()
+    stiffUpdate = NRStiffnessUpdateCriterion(mim)
+    stiffUpdate.setDebug(parameters['StiffUpdateDebug']); Default = False
+    stiffUpdate.setVerbose(parameters['StiffUpdateVerbose']); Default = False
+    mim.setStiffnessUpdateCriterion(stiffUpdate)
+</code>
+== StiffnessUpdateCriterion (Ancienne méthode de mise à jour de la matrice de raideur tangente) ==
+<code>
+    mim = metafor.getMechanicalIterationManager()
+    stiffUpdate = StiffnessUpdateCriterion(mim)
+    stiffUpdate.setNbOfIterationsForStiffUpdate(parameters['NbOfIterationsForStiffUpdate']); Default = 1
+    stiffUpdate.setCpuDependency(parameters['CpuDependency']); Default = False
+    stiffUpdate.setDebug(parameters['StiffUpdateDebug']); Default = False
+    stiffUpdate.setVerbose(parameters['StiffUpdateVerbose']); Default = False
+    mim.setStiffnessUpdateCriterion(stiffUpdate)
+</code>
+== CpuModifiedNRStiffnessUpdateCriterion ==
+<code>
+    mim = metafor.getMechanicalIterationManager()
+    stiffUpdate = CpuModifiedNRStiffnessUpdateCriterion(mim)
+    stiffUpdate.setKeepStepStiffness(parameters['KeepStepStiffness']); Default = True
+    stiffUpdate.setDebug(parameters['StiffUpdateDebug']); Default = False
+    stiffUpdate.setVerbose(parameters['StiffUpdateVerbose']); Default = False
+    mim.setStiffnessUpdateCriterion(stiffUpdate)
+</code>
+== PeriodicModifiedNRStiffnessUpdateCriterion ==
+<code>
+    mim = metafor.getMechanicalIterationManager()
+    stiffUpdate = PeriodicModifiedNRStiffnessUpdateCriterion(mim)
+    stiffUpdate.setMaxNbOfModifiedNRIterations(parameters['MaxNbOfModifiedNRIterations']); Default = 1
+    stiffUpdate.setResidualDecreaseFactor(parameters['ResidualDecreaseFactor']); Default = 1.0
+    stiffUpdate.setKeepStepStiffness(parameters['KeepStepStiffness']); Default = True
+    stiffUpdate.setDebug(parameters['StiffUpdateDebug']); Default = False
+    stiffUpdate.setVerbose(parameters['StiffUpdateVerbose']); Default = False
+    mim.setStiffnessUpdateCriterion(stiffUpdate)
+</code>
+Le paramètre keepStepStiffness permet de conserver la matrice de raideur tangente d'un pas de temps à l'autre, uniquement si le pas de temps est quasi-identique. Sinon, on force la mise à jour de la matrice de raideur tangente.
@@ Line 85: / Line 133: @@
 </code>
+Le paramètre keepAugmentationStiffness permet de conserver la matrice de raideur tangente après une augmentation, uniquement si le premier résidu avant l'augmentation est plus grand que le premier résidu après l'augmentation. Sinon, on force la mise à jour de la matrice de raideur tangente.
+Le paramètre augmentationThreshold permet de spécifier à partir de quel numéro d'augmentation, nous passons d'une méthode du type Newton Raphson à une méthode du type Newton Raphson modifiée. Si cette valeur est prise à -1, nous commençons d'office par une méthode du type Newton Raphson modifiée avant la première augmentation.
+Finalement, l'algorithme est suffisamment malin pour passer à une méthode du type Newton Raphson modifiée si nous ne sommes pas en contact (avant et après impact, après springback), quelque soit la valeur de augmentationThreshold.
 ===== Formulation AIC =====
@@ Line 131: / Line 183: @@
 </code>
- --- //[[gwautelet@ulg.ac.be|gaëtan]] 2017/04/22 22:37//
+ --- //[[gwautelet@ulg.ac.be|gaëtan]] 2017/04/26 22:37//