Strains in shell and birck elements
Strains in shell and birck elements
(OP)
Hi to everyone,
I am an aerospace stress engineer and I need support in the following question:
I modeled a beam fixed-free, two times, in hypermesh:
First, with shell elements and the second one with brick elements.
Then, I applied the same vertical load at each free side of beams and I ran the model with a sol101 in Nastran.
Now, I'm looking the output with hyperview and I see somethings doesn't feel right:
I have, in the same points, the same stress values in beam direction (and in the other directions) for both the beams, and these values are also in according to the values that I find using the classical beam theory.
The problem is that I find at the same points different strains between the two beams (shell and brick).
So I did this check:
i applied the stress-strain relationship (sigma =E*eps) to find the stress in beam direction by the strain in this direction:
the stress calculated by the solid beam strain (in beam direction) are in according to the stress values for classical theory and in according with the stress result by the FEM.
But, obliviosly, the shell beam strain output and stress strain relationship give me a wrong values.
Why does this happen?
Thanks in advance.
I am an aerospace stress engineer and I need support in the following question:
I modeled a beam fixed-free, two times, in hypermesh:
First, with shell elements and the second one with brick elements.
Then, I applied the same vertical load at each free side of beams and I ran the model with a sol101 in Nastran.
Now, I'm looking the output with hyperview and I see somethings doesn't feel right:
I have, in the same points, the same stress values in beam direction (and in the other directions) for both the beams, and these values are also in according to the values that I find using the classical beam theory.
The problem is that I find at the same points different strains between the two beams (shell and brick).
So I did this check:
i applied the stress-strain relationship (sigma =E*eps) to find the stress in beam direction by the strain in this direction:
the stress calculated by the solid beam strain (in beam direction) are in according to the stress values for classical theory and in according with the stress result by the FEM.
But, obliviosly, the shell beam strain output and stress strain relationship give me a wrong values.
Why does this happen?
Thanks in advance.
RE: Strains in shell and birck elements
RE: Strains in shell and birck elements
it does seem "odd" to model something (a slab ?) as a shell and also as a solid. The shell collapses the thru thickness freedom (that the 3D model models); a shell has not analytical thickness, the 3D model has actual thickness. The fact at the stresses are consistent says that the shell element code is correctly accounting for things not modelled.
I imagine your stress/strain relationship is "wrong". In the first place there are probably 3D strains. In the second place I'd apply (or at least read up) the element coding ... how does the element convert strain to stress, both shell and 3D. Remember these elements typically assume a strain function (linear of a 4 node shell, parabolic for a 8 node shell) and stress is calculated from strain. And not necessarily stress = E*strain ... often they hide "easter eggs" there, factors accounting for other effects.
"Hoffen wir mal, dass alles gut geht !"
General Paulus, Nov 1942, outside Stalingrad after the launch of Operation Uranus.
RE: Strains in shell and birck elements
"Hoffen wir mal, dass alles gut geht !"
General Paulus, Nov 1942, outside Stalingrad after the launch of Operation Uranus.
RE: Strains in shell and birck elements
If you want us to help more, please share the results.
RE: Strains in shell and birck elements
I have just founded the solution; I was looking only at membrane strain and not at bending strain component (as FEA way He advised to do).
RE: Strains in shell and birck elements
Elements in solid mechanics use displacement and rotation interpolants in elements. In post-processing, the definitions of stress resultants (bending, twisting, shear, normal forces) - which are functions of derivatives of displacement and sometimes also rotation - for the applied dimension reduction model (bar, beam, plate, membrane, shell) are used to evaluate the stress resultants in each element and at nodes. Finally, normal stresses and shear stresses are evaluated by utilizing their definitions (functions of bending, shear etc.). Stress invariants (e.g., von Mises stress) are calculated in the final steps before display of results.
In the 3D model, no dimension reduction is performed and thus 3D kinematical (displacement-strain) and constitutive (strain-stress) relations are directly applied to give normal stresses, shear stresses and invariants.
Mixed models deviate from this, but they are rarely applied in commercial FEA software.