Update deal.II scripts

perpendicular-flap/solid-dealii/linear_elasticity.prm → perpendicular-flap/solid-dealii/parameters.prm

@@ -14,29 +14,28 @@ subsection Time
 end
 
 subsection Discretization
-  # Time integration scheme
-  # 0 = forward, 1 = backward
-  set theta               = 0.5
-
   # Polynomial degree of the FE system
   set Polynomial degree   = 4
 end
 
 subsection System properties
-  # mu (shear modulus)
-  set mu              = 1538462
+  # Poisson's ratio
+  set Poisson's ratio = 0.3
 
-  # lambda
-  set lambda          = 2307692
+  # Shear modulus
+  set Shear modulus   = 1538462
 
-  # density
+  # Density
   set rho             = 3000
 
-  # body forces x,y,z
+  # Body forces x,y,z
   set body forces     = 0.0,0.0,0.0
 end
 
-subsection Linear solver
+subsection Solver
+  # Structural model to be used: linear or neo-Hookean
+  set Model                     = linear
+
   # Linear solver iterations (multiples of the system matrix size)
   set Max iteration multiplier  = 1
 
@@ -45,6 +44,15 @@ subsection Linear solver
 
   # Linear solver: CG or Direct
   set Solver type               = Direct
+
+  # Number of Newton-Raphson iterations allowed
+  set Max iterations Newton-Raphson = 10
+
+  # Displacement error tolerance
+  set Tolerance displacement        = 1.0e-6
+
+  # Force residual tolerance
+  set Tolerance force               = 1.0e-9
 end
 
 subsection precice configuration

perpendicular-flap/solid-dealii/run.sh

@@ -7,21 +7,5 @@ cd ${0%/*} || exit 1    		    # Run from this directory
 # in another terminal.
 # These scripts present how the two participants would be started manually.
 
-# 1 for true, 0 for false
-nonlinear=0
-if [ "$1" = "-nonlinear" ]; then
-            nonlinear=1
-fi
+./elasticity parameters.prm
 
-linear=0
-if [ "$1" = "-linear" ]; then
-            linear=1
-fi
-
-if [ $linear -eq 1 ]; then
-        ./linear_elasticity linear_elasticity.prm
-elif [ $nonlinear -eq 1 ]; then
-        ./nonlinear_elasticity nonlinear_elasticity.prm
-else
-        echo "No solver type specified. Please specify -linear or -nonlinear as solver type"
-fi