I am trying to solve the linear elasticity problem. I have some questions:
The assembling of the stiffness-matrix is carried out in the integration points. Then we solve K u = f and get the displacements in the nodes.
Now I want to recover the stress field and therefore need the strain tensor:
epsilon = sym_grad(u)
to get the stress:
sigma = C1trace(epsilon)I+C2*epsilon
Where are the strains and stresses calculated in this case? Are they in nodes or in integration points?
Usually in commercial FE-Software the strains and stresses are calculated in the integration points whithin the elements.
How can I get the values in integration points?
Is there any document where I can read all of this?
Does anybody know how to show results of DG-field in paraview? It only works for DG0. DG1 fields are represented as CG1-field.
Thanks a lot.