Dear all,
I am working on determining constitutive response of a hyperelastic materials. I can calculate deformation and after doing post processing I can get element stress. However, I want to calculate stress on an arbitrary points and also also want to calculate stress in each nodal points.
here is the sketch of my algorithms:
V=TensorFunctionSpace(mesh, ‘DG’, 0) # Tensor space (fixed)
def stress:
return S . # PK2 stress tensor
PK2 = project(S, V)
s11, s12, s13, s21, s22, s23, s31, s32, s33 = PK2.split(True)
Sxx = s11.vector().get_local() # this is my element stress of S_{11}
Now, I want to determine this Sxx in each nodal points. interpolation is not working. Any suggestion how to make it through expression or interpolation or any other method? Thank you in advance.