This is a read only copy of the old FEniCS QA forum. Please visit the new QA forum to ask questions

Calculation and access of strains at integration points

+3 votes

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.

asked Mar 16, 2017 by titusf FEniCS Novice (570 points)
...