Metafor

ULiege - Aerospace & Mechanical Engineering

User Tools

Site Tools

You are here: metafor.ltas.ulg.ac.be » Commits » 2017 » Commit 2017-04-26
Trace: • Commit 2017-04-27 • Commit 2017-08-22
commit:2017:04_22

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:2017:04_22 [2017/04/25 22:40] – [Lagrangien Augmenté] wauteletcommit: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. 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 === === 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  ===  === 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. 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 54: Line 102:
  
 Voici les méthodes disponibles et les options associées. Voici les méthodes disponibles et les options associées.
 +
 +== AlmCpuModifiedNRStiffnessUpdateCriterion ==
  
 <code> <code>
 +    alm = StandardAugmentedLagrangianManager(metafor)
 +    mim = metafor.getMechanicalIterationManager()
 +    stiffUpdate = AlmCpuModifiedNRStiffnessUpdateCriterion(mim,alm)
 +    stiffUpdate.setAugmentationThreshold(parameters['AugmentationThreshold']); Default = 0
 +    stiffUpdate.setKeepStepStiffness(parameters['KeepStepStiffness']); Default = True
 +    stiffUpdate.setKeepAugmentationStiffness(parameters['KeepAugmentationStiffness']); Default = True
 +    stiffUpdate.setDebug(parameters['StiffUpdateDebug']); Default = False
 +    stiffUpdate.setVerbose(parameters['StiffUpdateVerbose']); Default = False
 +    mim.setStiffnessUpdateCriterion(stiffUpdate)
 +</code>
 +
 +== AlmPeriodicModifiedNRStiffnessUpdateCriterion ==
 +
 +<code>
 +    alm = StandardAugmentedLagrangianManager(metafor)
     mim = metafor.getMechanicalIterationManager()     mim = metafor.getMechanicalIterationManager()
     stiffUpdate = AlmPeriodicModifiedNRStiffnessUpdateCriterion(mim,alm)     stiffUpdate = AlmPeriodicModifiedNRStiffnessUpdateCriterion(mim,alm)
Line 63: Line 128:
     stiffUpdate.setKeepStepStiffness(parameters['KeepStepStiffness']); Default = True     stiffUpdate.setKeepStepStiffness(parameters['KeepStepStiffness']); Default = True
     stiffUpdate.setKeepAugmentationStiffness(parameters['KeepAugmentationStiffness']); Default = True     stiffUpdate.setKeepAugmentationStiffness(parameters['KeepAugmentationStiffness']); Default = True
-    stiffUpdate.setDebug(parameters['StiffUpdateDebug']) +    stiffUpdate.setDebug(parameters['StiffUpdateDebug']); Default = False 
-    stiffUpdate.setVerbose(parameters['StiffUpdateVerbose'])+    stiffUpdate.setVerbose(parameters['StiffUpdateVerbose']); Default = False
     mim.setStiffnessUpdateCriterion(stiffUpdate)     mim.setStiffnessUpdateCriterion(stiffUpdate)
 </code> </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 ===== ===== Formulation AIC =====
Line 112: Line 183:
 </code> </code>
  
- --- //[[gwautelet@ulg.ac.be|gaëtan]] 2017/04/22 22:37//+ --- //[[gwautelet@ulg.ac.be|gaëtan]] 2017/04/26 22:37//
  
  
commit/2017/04_22.1493152800.txt.gz · Last modified: by wautelet
Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki