Change generalized alpha

alpha_m and alpha_f must be 0 to obtain the Newmark Beta scheme. We currently cannot support alpha_f != 0 with preCICE as long as Waveform Iteration is not implemented.

Author: BenjaminRodenberg

perpendicular-flap/solid-fenics/perp-flap.py

@@ -79,8 +79,8 @@ def neumann_boundary(x, on_boundary):
 bc = DirichletBC(V, Constant((0, 0)), fixed_boundary)
 
 # alpha method parameters
-alpha_m = Constant(0.2)
-alpha_f = Constant(0.4)
+alpha_m = Constant(0)
+alpha_f = Constant(0)
 gamma = Constant(0.5 + alpha_f - alpha_m)
 beta = Constant((gamma + 0.5)**2 / 4.)