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?