I obtain wrong displacement values using "ElasticSolver" and I cannot understand why.
The geometry is made up of two identical bodies, one on top of the other:
The dimensions of each body are 0.3(x) x 0.05(y) x 0.05 (z).
The mesh is made up of second order bricks (0.01 x 0.01 x 0.01).
Material:
Youngs modulus = 5e6
Poisson ratio = 0.0 (=> constant cross-section, linear problem)
Forces:
Force 3 = 5e5 (on top assembly face)
Force 3 = -5e5 (on bottom assembly face)
On the common face between the two bodies (similar results by setting the Displacements):
Boundary Condition 1
Target Boundaries(1) = 2
Name = "Ref"
!Displacement 1 = 0.0
!Displacement 2 = 0.0
!Displacement 3 = 0.0
Spring 1 = 1e2
Spring 2 = 1e2
Spring 3 = 1e2
End
The correct displacements on top and bottom faces are +5e-3 and -5e-3 respectively in the z direction.
StressSolver and ElasticSolver stress_zz result: 5e5 everywhere (as expected)
StressSolver top face z displacement (min, mean, max): 4.963e-3, 4.999e-3, 5.034e-3 (nearly correct. Much better results by setting the Displacements to 0 on BC1)
ElasticSolver top face z displacement (min, mean, max): 4.401e-3 4.401e-3 4.401e-3 (WHY???????
![Shocked :shock:](./images/smilies/icon_eek.gif)
The z displacement obtained from the ElasticSolver is 88% the correct displacement everywhere. Why?
![Question :?:](./images/smilies/icon_question.gif)
The other stress and displacement components are virtually zero from both Solvers.
I cannot understand what's wrong.
![Confused :?](./images/smilies/icon_e_confused.gif)
Also, is there a way to output the Elemental Fields? Adding Save Elemental Fields = True in ResultOutputSolver does not make any difference. Aren't stress and strain obtained from the nodal displacement and computed on the elements (Gauss Points) and then nodal stress and strain are obtained from the nodal values?
![Question :?:](./images/smilies/icon_question.gif)
Elmersolver Version: 9.0 (Rev: Release, Compiled: 2024-04-26)
Thanks