From e0218e670bdd3337738cb9f67ee8f6d367159018 Mon Sep 17 00:00:00 2001
From: Jonathan Youett <youett@math.fu-berlin.de>
Date: Mon, 12 Nov 2018 10:07:08 +0100
Subject: [PATCH] Parameter file for the primal-dual solver

---
 resources/pridualcontact.parset | 124 ++++++++++++++++++++++++++++++++
 1 file changed, 124 insertions(+)
 create mode 100644 resources/pridualcontact.parset

diff --git a/resources/pridualcontact.parset b/resources/pridualcontact.parset
new file mode 100644
index 00000000..770be198
--- /dev/null
+++ b/resources/pridualcontact.parset
@@ -0,0 +1,124 @@
+
+path =
+resultPath = /storage/mi/akbib/results/contact/pridual_large_ipopt/
+
+gridName0 = box
+gridName1 = die
+
+minLevel = 0
+
+nCouplings = 1
+
+restart = 0
+loadingStep = 1
+
+loadingStepsVertical = 1 #0 # 5
+loadingStepsSlide = 1 #0
+
+writer = amira
+
+[box]
+  nRefines = 2 #2 #3
+  createGrid = 1
+  type = cuboid
+  lowerCorner = 0 0 0
+  upperCorner = 4 9 3
+  elements = 1 3 1 # 4 8 2
+  tetrahedral = 1
+
+  E = 1
+  nu = 0.3
+
+  localDirichletVertical = 0 0 0
+  localDirichletSlide = 0 0 0
+
+  [box.dirichlet]
+    constantValues = 0 0 0
+    nCriterions = 1
+    [box.dirichlet.criterion0]
+      type = LowerCentre
+      component = 2
+      bound = 0
+
+
+
+[die]
+  nRefines = 1 # 1
+  createGrid = 1
+  type = tubeSegment
+  center = -0.7 4.5 6.2
+  thickness = 0.2
+  length= 5.2
+  innerRadius= 3
+  fromAngle = 3.1415
+  toAngle = 6.283
+  nElementRing =  8 #6 #10 #18 is what I used #25 is too much
+  nElementLength = 3 #4 four is what I used #10
+  closeTube = 0
+  axis = 0
+  tetrahedra = 1
+  parameterizedBoundary = 0
+
+  E = 1000
+  nu = 0.3
+
+  localDirichletVertical = 0 0 -1.4
+  localDirichletSlide = 0 2.1 0
+
+  [die.dirichlet]
+    constantValues = 0 0 0 #-1.4
+    nCriterions = 1
+
+    [die.dirichlet.criterion0]
+      type = UpperCentre
+      component = 2
+      bound = 6.19
+
+[coupling0]
+  nonmortarGrid=box
+  mortarGrid=die
+  overlap=1.3
+  couplingDistance = 5
+  minNormalAngle = 2.2 # -0.1
+  coveredArea = 0.99
+  couplingType = CONTACT
+  backEnd = CONTACT
+
+  [coupling0.die]
+    nCriterions = 2
+
+    [coupling0.die.criterion0]
+      type = LowerCentre
+      component = 2
+      bound = 4.5
+
+    [coupling0.die.criterion1]
+      type = LowerNormal
+      component = 2
+      bound = -1e-1
+
+  [coupling0.box]
+    nCriterions = 1
+    [coupling0.box.criterion0]
+      type = UpperCentre
+      component = 2
+      bound = 2.9999
+
+[primaldual]
+  priDualConstant = 1 #1000
+  priDualFDeps = 1e-6
+  full_linearisation = 1
+  initActiveSetBound = 0.01 # 2 #0.01 #2
+  # verbosity of the filter method
+  verbosity = 2
+  #error tolerance
+  tolerance = 1e-7
+  # maximal iteration number
+  maxIterations =400
+  # multigrid iterations
+  localMaxIterations =1 #0 #3 #10 #500
+  # tolerance for the local filter problems
+  localTolerance = 1e-4
+  # the linear solver to be used within IpOpt, leave empty to use the default
+  linearSolverType = ma97
+ 
-- 
GitLab