By Evalf, and other Nutils contributors

officialelasticityFCMNURBS

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. FCM{nelems=20,degree=1}).