Linear elasticity- boundary conditions


I was trying to run a simple linear elastic problem of Disk under diamteral compression. The disk is model in GMSH with radius 20 with bottom node (0,-20) fully fixed (x and y) and top (node 0,19) x fixed with displacement allowed in -y.

The code is running but not output neither the stress fields and u as exported to para view. Attached Snip of the BCs

If i want to apply a point load at the top how should i modify?

Consider: Implement point source in the variational formulation - #2 by MiroK

Please also note that it is adviced to supply a minimal code example, encapsulated with 3x`. For more info on minimal examples, see: