I solve a problem and obtain my finite element result for a simple plane stress problem mentioned in this thread .
An example solution for a givn set of material properties looks like the image below (with original mesh)
Now I want to interpolate the solution at a point within the domain, in the un-deformed configuration,which is not necessarily a nodal point in the mesh. I figured I get some result if I try this :
u_fe = solve(A, u_fe.vector(), b) # obtain the solution
pt = (1.95, 1.0) # this point is in the undeformed configuration, and not a nodal point
u_fe_interpolate = u_fe(pt)
But I’m not sure if this is correct. It’ll be really useful to know if this is correct and in case it is not, can I get the result I intend by any other means ?
Thank you