Large strain compatibility of reduced integration hexaedral element?

[pierre]

Hi,

I would like to know if the reduced integration (one point) hexahedral of FEAP program is compatible with large strain applications, and if yes, on what theory is it based?

I thank you in advance for your help!
Best regards,

Pierre

[FEAP_Admin]

You can request one-point integration when defining a material but FEAP does not implement any stabilization so you computation will be unstable;
see for example section 6.1 of the manual.

[pierre]

Hi,

Thank you very much for your quick answer!
Best greetings,

Pierre

[Prof. R.L. Taylor]

fld3d1u is a finite deformation element.  Based on paper by Bonet & Bhargava, IJNME, vol 38, 1995, pp 2809--2828.

[pierre]

Hi,

Thank you very much for this hint.

Best greetings,

Pierre

[pierre]

Dear prof. R.L. Taylor,
I could read the very interesting reference you indicate.
Please, i am curious to know how is managed the hourglass parameter k?
Is it a user parameter, set by default to 0, 01 which is the value given in the reference, or do you think a value expressed in terms of material or element geometry can be derived?
I thank you in advance for your helpful comments.
Best greetings,
Pierre

[Prof. R.L. Taylor]

Yes, you can set the parameter using

   DISP UNIF alpha

where alpha is your parameter.  If zero it is reset to 0.001 and the stabilized tangent is proportional to the initial tangent.  See the file fld3d1u.f in /elements/solid3d

[pierre]

Dear Prof. R. L. Taylor,

Thank you very much for your reply.
I fear this element is not part of FEAPpv source files that I could download here (version 3.1).
In /elements/solid3d.f, I sould successfully identify  fld3d1 subroutine, but no fld3d1u.
I guess that the u in fld3d1u stands for under-integraiton, right?
Is there perhaps any possibility to have a copy of this soure file?

Thank you again very much for your help.
Have a very good day.
Best greetings.

Pierre

[FEAP_Admin]

I think this routine is only available in the full version of FEAP.  If you want to do this in FEAPpv you will
need to program it yourself.

[pierre]

Hi,

Thanks for the inputs.
I'm about to look into this.
Thanks again.
Best regards,

Pierre

[Prof. R.L. Taylor]

The element integrates volume by 2x2x2 quadrature but evaluates material at only none point.  You can modify the fld3d1 element for most of it.