By Evalf, and other Nutils contributors
Solves the linear plane strain elasticity problem for an infinite plate with a circular hole under tension. We do this by placing the circle in the origin of a unit square, imposing symmetry conditions on the left and bottom ("sym"), and Dirichlet conditions constraining the displacements to the analytical solution to the right and top ("far").
The script can be run in two modes: by specifying
mode=FCM, the circular
hole is cut out of a regular finite element mesh by means of the Finite Cell
Method; by specifying
mode=NURBS a Non-Uniform Rational BSpline geometry is
created to map a regular domain onto the desired shape. Either mode supports
sub-parameters which can be specified from the command-line by attaching them
in curly braces (e.g.