Remove Neumann bc from fenics-example

Author: srosenbu

benchmarks/linear-elasticity-plate-with-hole/data/output_003125_p0_000000.vtu

@@ -5,6 +5,7 @@
 import numpy as np
 import ufl
 from dolfinx.fem.petsc import LinearProblem
+from petsc4py.PETSc import ScalarType
 from mpi4py import MPI
 from pint import UnitRegistry
 
@@ -180,8 +181,8 @@ def as_tensor(v):
 a = df.fem.form(ufl.inner(sigma(u_), eps(v_)) * dx)
 
 
-f = df.fem.form(ufl.inner(ufl.dot(stress_function, ufl.FacetNormal(mesh)), u_) * ufl.ds)
-#f = df.fem.form(ufl.inner(ufl.Constant(mesh, ((0.0),(0.0))), u_) * ufl.ds)
+# set rhs to zero
+f = df.fem.form(ufl.inner(df.fem.Constant(mesh, np.array([0.0,0.0])), u_) * ufl.ds)
 
 bc_right = df.fem.dirichletbc(u_prescribed, dofs_right)
 bc_top = df.fem.dirichletbc(u_prescribed, dofs_top)

benchmarks/linear-elasticity-plate-with-hole/data/output_003125_p0_000001.vtu

@@ -4,5 +4,5 @@
     "input": "data/input_003125.json",
     "mesh": "data/mesh_003125.msh",
     "output": "data/output_003125.vtk",
-    "max_mises_stress": 299783353.3539022
+    "max_mises_stress": 299783353.3384597
 }

benchmarks/linear-elasticity-plate-with-hole/data/output_0125_p0_000000.vtu

@@ -4,5 +4,5 @@
     "input": "data/input_0125.json",
     "mesh": "data/mesh_0125.msh",
     "output": "data/output_0125.vtk",
-    "max_mises_stress": 300129622.7217106
+    "max_mises_stress": 300129622.7217104
 }

benchmarks/linear-elasticity-plate-with-hole/data/output_0125_p0_000001.vtu

@@ -4,5 +4,5 @@
     "input": "data/input_025.json",
     "mesh": "data/mesh_025.msh",
     "output": "data/output_025.vtk",
-    "max_mises_stress": 299791507.5586348
+    "max_mises_stress": 299791507.5586332
 }