Hyperelasticity#
Author: Jørgen S. Dokken and Garth N. Wells
This section shows how to solve the hyperelasticity problem for deformation of a beam.
We will also show how to create a constant boundary condition for a vector function space.
We start by importing DOLFINx and some additional dependencies.
Then, we create a slender cantilever consisting of hexahedral elements and create the function space V for our unknown.
from dolfinx import log, default_scalar_type
from dolfinx.fem.petsc import NonlinearProblem
from dolfinx.nls.petsc import NewtonSolver
import matplotlib.pyplot as plt
import pyvista
from dolfinx import nls
import numpy as np
import ufl
from mpi4py import MPI
from dolfinx import fem, mesh, plot
L = 20.0
domain = mesh.create_box(MPI.COMM_WORLD, [[0.0, 0.0, 0.0], [L, 1, 1]], [20, 5, 5], mesh.CellType.hexahedron)
V = fem.VectorFunctionSpace(domain, ("Lagrange", 2))
WARNING:py.warnings:/tmp/ipykernel_1859/1146139302.py:14: DeprecationWarning: This method is deprecated. Use FunctionSpace with an element shape argument instead
V = fem.VectorFunctionSpace(domain, ("Lagrange", 2))
We create two python functions for determining the facets to apply boundary conditions to
def left(x):
return np.isclose(x[0], 0)
def right(x):
return np.isclose(x[0], L)
fdim = domain.topology.dim - 1
left_facets = mesh.locate_entities_boundary(domain, fdim, left)
right_facets = mesh.locate_entities_boundary(domain, fdim, right)
Next, we create a marker based on these two functions
# Concatenate and sort the arrays based on facet indices. Left facets marked with 1, right facets with two
marked_facets = np.hstack([left_facets, right_facets])
marked_values = np.hstack([np.full_like(left_facets, 1), np.full_like(right_facets, 2)])
sorted_facets = np.argsort(marked_facets)
facet_tag = mesh.meshtags(domain, fdim, marked_facets[sorted_facets], marked_values[sorted_facets])
We then create a function for supplying the boundary condition on the left side, which is fixed.
u_bc = np.array((0,) * domain.geometry.dim, dtype=default_scalar_type)
To apply the boundary condition, we identity the dofs located on the facets marked by the MeshTag.
left_dofs = fem.locate_dofs_topological(V, facet_tag.dim, facet_tag.find(1))
bcs = [fem.dirichletbc(u_bc, left_dofs, V)]
Next, we define the body force on the reference configuration (B), and nominal (first Piola-Kirchhoff) traction (T).
B = fem.Constant(domain, default_scalar_type((0, 0, 0)))
T = fem.Constant(domain, default_scalar_type((0, 0, 0)))
Define the test and solution functions on the space \(V\)
v = ufl.TestFunction(V)
u = fem.Function(V)
Define kinematic quantities used in the problem
# Spatial dimension
d = len(u)
# Identity tensor
I = ufl.variable(ufl.Identity(d))
# Deformation gradient
F = ufl.variable(I + ufl.grad(u))
# Right Cauchy-Green tensor
C = ufl.variable(F.T * F)
# Invariants of deformation tensors
Ic = ufl.variable(ufl.tr(C))
J = ufl.variable(ufl.det(F))
Define the elasticity model via a stored strain energy density function \(\psi\), and create the expression for the first Piola-Kirchhoff stress:
# Elasticity parameters
E = default_scalar_type(1.0e4)
nu = default_scalar_type(0.3)
mu = fem.Constant(domain, E / (2 * (1 + nu)))
lmbda = fem.Constant(domain, E * nu / ((1 + nu) * (1 - 2 * nu)))
# Stored strain energy density (compressible neo-Hookean model)
psi = (mu / 2) * (Ic - 3) - mu * ufl.ln(J) + (lmbda / 2) * (ufl.ln(J))**2
# Stress
# Hyper-elasticity
P = ufl.diff(psi, F)
Comparison to linear elasticity
To illustrate the difference between linear and hyperelasticity, the following lines can be uncommented to solve the linear elasticity problem.
# P = 2.0 * mu * ufl.sym(ufl.grad(u)) + lmbda * ufl.tr(ufl.sym(ufl.grad(u))) * I
Define the variational form with traction integral over all facets with value 2. We set the quadrature degree for the integrals to 4.
metadata = {"quadrature_degree": 4}
ds = ufl.Measure('ds', domain=domain, subdomain_data=facet_tag, metadata=metadata)
dx = ufl.Measure("dx", domain=domain, metadata=metadata)
# Define form F (we want to find u such that F(u) = 0)
F = ufl.inner(ufl.grad(v), P) * dx - ufl.inner(v, B) * dx - ufl.inner(v, T) * ds(2)
As the varitional form is non-linear and written on residual form, we use the non-linear problem class from DOLFINx to set up required structures to use a Newton solver.
problem = NonlinearProblem(F, u, bcs)
and then create and customize the Newton solver
solver = NewtonSolver(domain.comm, problem)
# Set Newton solver options
solver.atol = 1e-8
solver.rtol = 1e-8
solver.convergence_criterion = "incremental"
We create a function to plot the solution at each time step.
pyvista.start_xvfb()
plotter = pyvista.Plotter()
plotter.open_gif("deformation.gif", fps=3)
topology, cells, geometry = plot.vtk_mesh(u.function_space)
function_grid = pyvista.UnstructuredGrid(topology, cells, geometry)
values = np.zeros((geometry.shape[0], 3))
values[:, :len(u)] = u.x.array.reshape(geometry.shape[0], len(u))
function_grid["u"] = values
function_grid.set_active_vectors("u")
# Warp mesh by deformation
warped = function_grid.warp_by_vector("u", factor=1)
warped.set_active_vectors("u")
# Add mesh to plotter and visualize
actor = plotter.add_mesh(warped, show_edges=True, lighting=False, clim=[0, 10])
# Compute magnitude of displacement to visualize in GIF
Vs = fem.FunctionSpace(domain, ("Lagrange", 2))
magnitude = fem.Function(Vs)
us = fem.Expression(ufl.sqrt(sum([u[i]**2 for i in range(len(u))])), Vs.element.interpolation_points())
magnitude.interpolate(us)
warped["mag"] = magnitude.x.array
Finally, we solve the problem over several time steps, updating the y-component of the traction
log.set_log_level(log.LogLevel.INFO)
tval0 = -1.5
for n in range(1, 10):
T.value[2] = n * tval0
num_its, converged = solver.solve(u)
assert (converged)
u.x.scatter_forward()
print(f"Time step {n}, Number of iterations {num_its}, Load {T.value}")
function_grid["u"][:, :len(u)] = u.x.array.reshape(geometry.shape[0], len(u))
magnitude.interpolate(us)
warped.set_active_scalars("mag")
warped_n = function_grid.warp_by_vector(factor=1)
plotter.update_coordinates(warped_n.points.copy(), render=False)
plotter.update_scalar_bar_range([0, 10])
plotter.update_scalars(magnitude.x.array)
plotter.write_frame()
plotter.close()
2023-11-27 08:04:43.008 ( 5.900s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:43.393 ( 6.284s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.384765, 0.370000, 0.020000 (PETSc Krylov solver)
2023-11-27 08:04:43.540 ( 6.431s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:43.842 ( 6.733s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.301845, 0.280000, 0.020000 (PETSc Krylov solver)
2023-11-27 08:04:43.849 ( 6.740s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 2: r (abs) = 22.2455 (tol = 1e-08) r (rel) = 0.134278(tol = 1e-08)
2023-11-27 08:04:43.989 ( 6.880s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:44.287 ( 7.179s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.298347, 0.280000, 0.020000 (PETSc Krylov solver)
2023-11-27 08:04:44.295 ( 7.186s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 3: r (abs) = 2.43261 (tol = 1e-08) r (rel) = 0.0146837(tol = 1e-08)
2023-11-27 08:04:44.436 ( 7.328s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:44.734 ( 7.625s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.297636, 0.300000, 0.000000 (PETSc Krylov solver)
2023-11-27 08:04:44.742 ( 7.633s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 4: r (abs) = 4.43158 (tol = 1e-08) r (rel) = 0.0267498(tol = 1e-08)
2023-11-27 08:04:44.883 ( 7.774s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:45.182 ( 8.073s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.299429, 0.300000, 0.000000 (PETSc Krylov solver)
2023-11-27 08:04:45.189 ( 8.081s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 5: r (abs) = 0.144189 (tol = 1e-08) r (rel) = 0.000870353(tol = 1e-08)
2023-11-27 08:04:45.331 ( 8.222s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:45.633 ( 8.524s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.302116, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:04:45.640 ( 8.532s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 6: r (abs) = 0.021424 (tol = 1e-08) r (rel) = 0.000129319(tol = 1e-08)
2023-11-27 08:04:45.782 ( 8.673s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:46.081 ( 8.972s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.298337, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:04:46.087 ( 8.978s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 7: r (abs) = 4.80067e-06 (tol = 1e-08) r (rel) = 2.89777e-08(tol = 1e-08)
2023-11-27 08:04:46.227 ( 9.118s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:46.524 ( 9.415s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.297452, 0.290000, 0.000000 (PETSc Krylov solver)
2023-11-27 08:04:46.531 ( 9.422s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 8: r (abs) = 2.69663e-11 (tol = 1e-08) r (rel) = 1.62774e-13(tol = 1e-08)
2023-11-27 08:04:46.531 ( 9.422s) [main ] NewtonSolver.cpp:255 INFO| Newton solver finished in 8 iterations and 8 linear solver iterations.
Time step 1, Number of iterations 8, Load [ 0. 0. -1.5]
2023-11-27 08:04:46.857 ( 9.748s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:47.157 ( 10.048s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.299488, 0.300000, 0.000000 (PETSc Krylov solver)
2023-11-27 08:04:47.307 ( 10.198s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:47.607 ( 10.498s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.300404, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:04:47.614 ( 10.506s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 2: r (abs) = 17.3254 (tol = 1e-08) r (rel) = 0.117842(tol = 1e-08)
2023-11-27 08:04:47.756 ( 10.647s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:48.055 ( 10.947s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.299286, 0.300000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:04:48.063 ( 10.954s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 3: r (abs) = 5.14882 (tol = 1e-08) r (rel) = 0.0350207(tol = 1e-08)
2023-11-27 08:04:48.204 ( 11.095s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:48.502 ( 11.393s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.297947, 0.280000, 0.020000 (PETSc Krylov solver)
2023-11-27 08:04:48.509 ( 11.400s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 4: r (abs) = 7.24003 (tol = 1e-08) r (rel) = 0.0492445(tol = 1e-08)
2023-11-27 08:04:48.648 ( 11.540s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:48.948 ( 11.839s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.299542, 0.290000, 0.020000 (PETSc Krylov solver)
2023-11-27 08:04:48.956 ( 11.847s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 5: r (abs) = 0.777889 (tol = 1e-08) r (rel) = 0.00529096(tol = 1e-08)
2023-11-27 08:04:49.096 ( 11.987s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:49.394 ( 12.286s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.298745, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:04:49.402 ( 12.293s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 6: r (abs) = 1.25525 (tol = 1e-08) r (rel) = 0.00853785(tol = 1e-08)
2023-11-27 08:04:49.542 ( 12.434s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:49.840 ( 12.731s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.297336, 0.290000, 0.000000 (PETSc Krylov solver)
2023-11-27 08:04:49.847 ( 12.738s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 7: r (abs) = 0.00849512 (tol = 1e-08) r (rel) = 5.77813e-05(tol = 1e-08)
2023-11-27 08:04:49.988 ( 12.879s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:50.284 ( 13.175s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.296192, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:04:50.292 ( 13.183s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 8: r (abs) = 0.000192107 (tol = 1e-08) r (rel) = 1.30665e-06(tol = 1e-08)
2023-11-27 08:04:50.432 ( 13.323s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:50.730 ( 13.622s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.298179, 0.300000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:04:50.738 ( 13.629s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 9: r (abs) = 1.70374e-10 (tol = 1e-08) r (rel) = 1.15883e-12(tol = 1e-08)
2023-11-27 08:04:50.738 ( 13.629s) [main ] NewtonSolver.cpp:255 INFO| Newton solver finished in 9 iterations and 9 linear solver iterations.
Time step 2, Number of iterations 9, Load [ 0. 0. -3.]
2023-11-27 08:04:50.952 ( 13.843s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:51.250 ( 14.141s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.297970, 0.280000, 0.030000 (PETSc Krylov solver)
2023-11-27 08:04:51.397 ( 14.288s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:51.694 ( 14.586s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.297618, 0.300000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:04:51.702 ( 14.593s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 2: r (abs) = 10.0011 (tol = 1e-08) r (rel) = 0.0887471(tol = 1e-08)
2023-11-27 08:04:51.840 ( 14.732s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:52.137 ( 15.029s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.296919, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:04:52.144 ( 15.036s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 3: r (abs) = 5.33026 (tol = 1e-08) r (rel) = 0.0472992(tol = 1e-08)
2023-11-27 08:04:52.284 ( 15.176s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:52.583 ( 15.474s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.298380, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:04:52.590 ( 15.481s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 4: r (abs) = 11.9901 (tol = 1e-08) r (rel) = 0.106397(tol = 1e-08)
2023-11-27 08:04:52.730 ( 15.621s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:53.028 ( 15.919s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.297529, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:04:53.035 ( 15.926s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 5: r (abs) = 2.29702 (tol = 1e-08) r (rel) = 0.0203831(tol = 1e-08)
2023-11-27 08:04:53.176 ( 16.067s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:53.474 ( 16.365s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.298293, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:04:53.482 ( 16.373s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 6: r (abs) = 3.90234 (tol = 1e-08) r (rel) = 0.0346282(tol = 1e-08)
2023-11-27 08:04:53.631 ( 16.523s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:53.931 ( 16.822s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.299952, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:04:53.938 ( 16.829s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 7: r (abs) = 0.236535 (tol = 1e-08) r (rel) = 0.00209895(tol = 1e-08)
2023-11-27 08:04:54.078 ( 16.969s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:54.374 ( 17.265s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.296061, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:04:54.382 ( 17.273s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 8: r (abs) = 0.0427142 (tol = 1e-08) r (rel) = 0.000379034(tol = 1e-08)
2023-11-27 08:04:54.523 ( 17.414s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:54.820 ( 17.711s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.296901, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:04:54.828 ( 17.719s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 9: r (abs) = 2.87798e-05 (tol = 1e-08) r (rel) = 2.55384e-07(tol = 1e-08)
2023-11-27 08:04:54.968 ( 17.860s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:55.266 ( 18.158s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.297977, 0.280000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:04:55.274 ( 18.165s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 10: r (abs) = 6.08365e-10 (tol = 1e-08) r (rel) = 5.39846e-12(tol = 1e-08)
2023-11-27 08:04:55.274 ( 18.165s) [main ] NewtonSolver.cpp:255 INFO| Newton solver finished in 10 iterations and 10 linear solver iterations.
2023-11-27 08:04:55.489 ( 18.380s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:55.787 ( 18.678s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.298607, 0.280000, 0.030000 (PETSc Krylov solver)
Time step 3, Number of iterations 10, Load [ 0. 0. -4.5]
2023-11-27 08:04:55.935 ( 18.826s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:56.231 ( 19.122s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.296406, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:04:56.239 ( 19.130s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 2: r (abs) = 5.50693 (tol = 1e-08) r (rel) = 0.0653918(tol = 1e-08)
2023-11-27 08:04:56.378 ( 19.269s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:56.675 ( 19.566s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.296943, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:04:56.682 ( 19.574s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 3: r (abs) = 26.2489 (tol = 1e-08) r (rel) = 0.311692(tol = 1e-08)
2023-11-27 08:04:56.822 ( 19.713s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:57.114 ( 20.005s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.292512, 0.290000, 0.000000 (PETSc Krylov solver)
2023-11-27 08:04:57.122 ( 20.013s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 4: r (abs) = 2.30927 (tol = 1e-08) r (rel) = 0.0274213(tol = 1e-08)
2023-11-27 08:04:57.260 ( 20.152s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:57.555 ( 20.446s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.294071, 0.290000, 0.000000 (PETSc Krylov solver)
2023-11-27 08:04:57.562 ( 20.453s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 5: r (abs) = 14.0562 (tol = 1e-08) r (rel) = 0.16691(tol = 1e-08)
2023-11-27 08:04:57.702 ( 20.593s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:58.000 ( 20.891s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.297415, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:04:58.007 ( 20.898s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 6: r (abs) = 0.222774 (tol = 1e-08) r (rel) = 0.00264532(tol = 1e-08)
2023-11-27 08:04:58.146 ( 21.038s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:58.444 ( 21.335s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.297913, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:04:58.451 ( 21.342s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 7: r (abs) = 0.286671 (tol = 1e-08) r (rel) = 0.00340406(tol = 1e-08)
2023-11-27 08:04:58.592 ( 21.483s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:58.890 ( 21.781s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.297929, 0.290000, 0.000000 (PETSc Krylov solver)
2023-11-27 08:04:58.898 ( 21.789s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 8: r (abs) = 0.000321869 (tol = 1e-08) r (rel) = 3.82203e-06(tol = 1e-08)
2023-11-27 08:04:59.039 ( 21.930s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:59.337 ( 22.228s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.298411, 0.280000, 0.020000 (PETSc Krylov solver)
Time step 4, Number of iterations 9, Load [ 0. 0. -6.]
2023-11-27 08:04:59.345 ( 22.236s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 9: r (abs) = 2.63797e-07 (tol = 1e-08) r (rel) = 3.13244e-09(tol = 1e-08)
2023-11-27 08:04:59.345 ( 22.236s) [main ] NewtonSolver.cpp:255 INFO| Newton solver finished in 9 iterations and 9 linear solver iterations.
2023-11-27 08:04:59.559 ( 22.450s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:04:59.853 ( 22.745s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.294555, 0.270000, 0.020000 (PETSc Krylov solver)
2023-11-27 08:05:00.002 ( 22.893s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:00.297 ( 23.188s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.294291, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:05:00.304 ( 23.195s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 2: r (abs) = 3.19462 (tol = 1e-08) r (rel) = 0.0496479(tol = 1e-08)
2023-11-27 08:05:00.445 ( 23.336s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:00.743 ( 23.634s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.298289, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:05:00.750 ( 23.641s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 3: r (abs) = 7.71429 (tol = 1e-08) r (rel) = 0.119888(tol = 1e-08)
2023-11-27 08:05:00.889 ( 23.780s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:01.184 ( 24.075s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.294400, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:05:01.192 ( 24.083s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 4: r (abs) = 0.850873 (tol = 1e-08) r (rel) = 0.0132235(tol = 1e-08)
2023-11-27 08:05:01.333 ( 24.224s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:01.628 ( 24.519s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.295229, 0.300000, 0.000000 (PETSc Krylov solver)
2023-11-27 08:05:01.635 ( 24.526s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 5: r (abs) = 0.371434 (tol = 1e-08) r (rel) = 0.0057725(tol = 1e-08)
2023-11-27 08:05:01.774 ( 24.665s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:02.073 ( 24.964s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.298723, 0.290000, 0.020000 (PETSc Krylov solver)
2023-11-27 08:05:02.080 ( 24.971s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 6: r (abs) = 0.00215066 (tol = 1e-08) r (rel) = 3.34236e-05(tol = 1e-08)
2023-11-27 08:05:02.220 ( 25.111s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:02.518 ( 25.410s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.298100, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:05:02.526 ( 25.417s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 7: r (abs) = 2.54607e-06 (tol = 1e-08) r (rel) = 3.95687e-08(tol = 1e-08)
2023-11-27 08:05:02.667 ( 25.559s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:02.966 ( 25.857s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.298741, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:05:02.976 ( 25.867s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 8: r (abs) = 4.47166e-13 (tol = 1e-08) r (rel) = 6.94945e-15(tol = 1e-08)
2023-11-27 08:05:02.976 ( 25.868s) [main ] NewtonSolver.cpp:255 INFO| Newton solver finished in 8 iterations and 8 linear solver iterations.
Time step 5, Number of iterations 8, Load [ 0. 0. -7.5]
2023-11-27 08:05:03.188 ( 26.080s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:03.485 ( 26.376s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.296906, 0.260000, 0.030000 (PETSc Krylov solver)
2023-11-27 08:05:03.632 ( 26.523s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:03.930 ( 26.821s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.297361, 0.280000, 0.020000 (PETSc Krylov solver)
2023-11-27 08:05:03.937 ( 26.829s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 2: r (abs) = 2.00649 (tol = 1e-08) r (rel) = 0.0395622(tol = 1e-08)
2023-11-27 08:05:04.079 ( 26.970s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:04.376 ( 27.267s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.296984, 0.300000, 0.000000 (PETSc Krylov solver)
2023-11-27 08:05:04.383 ( 27.274s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 3: r (abs) = 4.60977 (tol = 1e-08) r (rel) = 0.0908914(tol = 1e-08)
2023-11-27 08:05:04.524 ( 27.416s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:04.822 ( 27.713s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.297422, 0.300000, 0.000000 (PETSc Krylov solver)
2023-11-27 08:05:04.830 ( 27.721s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 4: r (abs) = 0.185372 (tol = 1e-08) r (rel) = 0.00365501(tol = 1e-08)
2023-11-27 08:05:04.970 ( 27.861s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:05.265 ( 28.156s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.294954, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:05:05.272 ( 28.163s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 5: r (abs) = 0.024688 (tol = 1e-08) r (rel) = 0.000486777(tol = 1e-08)
2023-11-27 08:05:05.414 ( 28.306s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:05.710 ( 28.601s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.295587, 0.290000, 0.000000 (PETSc Krylov solver)
2023-11-27 08:05:05.717 ( 28.609s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 6: r (abs) = 5.69254e-06 (tol = 1e-08) r (rel) = 1.12241e-07(tol = 1e-08)
2023-11-27 08:05:05.857 ( 28.748s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:06.151 ( 29.042s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.294409, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:05:06.158 ( 29.049s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 7: r (abs) = 2.68147e-11 (tol = 1e-08) r (rel) = 5.2871e-13(tol = 1e-08)
2023-11-27 08:05:06.158 ( 29.049s) [main ] NewtonSolver.cpp:255 INFO| Newton solver finished in 7 iterations and 7 linear solver iterations.
Time step 6, Number of iterations 7, Load [ 0. 0. -9.]
2023-11-27 08:05:06.368 ( 29.259s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:06.663 ( 29.555s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.295333, 0.280000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:05:06.810 ( 29.702s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:07.106 ( 29.997s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.295327, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:05:07.113 ( 30.004s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 2: r (abs) = 1.38506 (tol = 1e-08) r (rel) = 0.0336622(tol = 1e-08)
2023-11-27 08:05:07.254 ( 30.145s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:07.548 ( 30.439s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.293852, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:05:07.555 ( 30.446s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 3: r (abs) = 3.03739 (tol = 1e-08) r (rel) = 0.07382(tol = 1e-08)
2023-11-27 08:05:07.696 ( 30.587s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:07.991 ( 30.882s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.295045, 0.290000, 0.000000 (PETSc Krylov solver)
2023-11-27 08:05:07.998 ( 30.889s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 4: r (abs) = 0.0412386 (tol = 1e-08) r (rel) = 0.00100225(tol = 1e-08)
2023-11-27 08:05:08.137 ( 31.028s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:08.433 ( 31.324s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.295777, 0.280000, 0.020000 (PETSc Krylov solver)
2023-11-27 08:05:08.440 ( 31.331s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 5: r (abs) = 0.00205057 (tol = 1e-08) r (rel) = 4.98364e-05(tol = 1e-08)
2023-11-27 08:05:08.582 ( 31.473s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:08.878 ( 31.769s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.296206, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:05:08.885 ( 31.776s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 6: r (abs) = 1.78866e-08 (tol = 1e-08) r (rel) = 4.34711e-10(tol = 1e-08)
2023-11-27 08:05:08.885 ( 31.776s) [main ] NewtonSolver.cpp:255 INFO| Newton solver finished in 6 iterations and 6 linear solver iterations.
Time step 7, Number of iterations 6, Load [ 0. 0. -10.5]
2023-11-27 08:05:09.095 ( 31.986s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:09.390 ( 32.281s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.294914, 0.280000, 0.000000 (PETSc Krylov solver)
2023-11-27 08:05:09.537 ( 32.429s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:09.831 ( 32.723s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.294083, 0.300000, 0.000000 (PETSc Krylov solver)
2023-11-27 08:05:09.839 ( 32.730s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 2: r (abs) = 1.06336 (tol = 1e-08) r (rel) = 0.031085(tol = 1e-08)
2023-11-27 08:05:09.978 ( 32.869s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:10.272 ( 33.163s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.294070, 0.280000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:05:10.279 ( 33.171s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 3: r (abs) = 2.0477 (tol = 1e-08) r (rel) = 0.0598598(tol = 1e-08)
2023-11-27 08:05:10.419 ( 33.310s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:10.716 ( 33.607s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.296511, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:05:10.723 ( 33.614s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 4: r (abs) = 0.00897719 (tol = 1e-08) r (rel) = 0.000262427(tol = 1e-08)
2023-11-27 08:05:10.863 ( 33.755s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:11.159 ( 34.050s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.295583, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:05:11.166 ( 34.057s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 5: r (abs) = 0.000167422 (tol = 1e-08) r (rel) = 4.89419e-06(tol = 1e-08)
2023-11-27 08:05:11.307 ( 34.198s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:11.604 ( 34.495s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.297155, 0.290000, 0.000000 (PETSc Krylov solver)
Time step 8, Number of iterations 6, Load [ 0. 0. -12.]
2023-11-27 08:05:11.612 ( 34.503s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 6: r (abs) = 3.25048e-11 (tol = 1e-08) r (rel) = 9.50202e-13(tol = 1e-08)
2023-11-27 08:05:11.612 ( 34.503s) [main ] NewtonSolver.cpp:255 INFO| Newton solver finished in 6 iterations and 6 linear solver iterations.
2023-11-27 08:05:11.822 ( 34.713s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:12.117 ( 35.008s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.295178, 0.280000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:05:12.266 ( 35.157s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:12.559 ( 35.451s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.293120, 0.290000, 0.000000 (PETSc Krylov solver)
2023-11-27 08:05:12.567 ( 35.458s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 2: r (abs) = 0.898789 (tol = 1e-08) r (rel) = 0.0309666(tol = 1e-08)
2023-11-27 08:05:12.708 ( 35.599s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:13.005 ( 35.896s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.296976, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:05:13.012 ( 35.903s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 3: r (abs) = 1.38354 (tol = 1e-08) r (rel) = 0.0476679(tol = 1e-08)
2023-11-27 08:05:13.154 ( 36.045s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:13.450 ( 36.341s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.296261, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:05:13.457 ( 36.348s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 4: r (abs) = 0.00185096 (tol = 1e-08) r (rel) = 6.37724e-05(tol = 1e-08)
2023-11-27 08:05:13.597 ( 36.488s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:13.894 ( 36.785s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.296986, 0.290000, 0.010000 (PETSc Krylov solver)
2023-11-27 08:05:13.901 ( 36.792s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 5: r (abs) = 7.87183e-06 (tol = 1e-08) r (rel) = 2.71213e-07(tol = 1e-08)
2023-11-27 08:05:14.041 ( 36.932s) [main ] petsc.cpp:698 INFO| PETSc Krylov solver starting to solve system.
2023-11-27 08:05:14.338 ( 37.229s) [main ] TimeLogger.cpp:28 INFO| Elapsed wall, usr, sys time: 0.296828, 0.290000, 0.000000 (PETSc Krylov solver)
2023-11-27 08:05:14.345 ( 37.236s) [main ] NewtonSolver.cpp:36 INFO| Newton iteration 6: r (abs) = 2.81503e-13 (tol = 1e-08) r (rel) = 9.69881e-15(tol = 1e-08)
2023-11-27 08:05:14.345 ( 37.236s) [main ] NewtonSolver.cpp:255 INFO| Newton solver finished in 6 iterations and 6 linear solver iterations.
Time step 9, Number of iterations 6, Load [ 0. 0. -13.5]
