Difference between revisions of "Contrib:KeesWouters/solids/mooney-rivlin"

Revision as of 14:52, 30 September 2011

Material behaviour Mooney-Rivlin

Simple block of Mooney-Rivlin material

september 2011 - Salome 6.3 - Code Aster VERSION DE DEVELOPPEMENT 11.00.10 - Ubuntu 11.04

This is a simple example of a block of Mooney-Rivlin material. As usual, the input and partly the output files can be found at the end of this contribution.
Errors are all mine.
Creative criticism welcome (of course I rather have positive than negative remarks).
Have fun.

Geometry

This is really annoying. A block of 5x12x43 [mm]. My only excuse is that it does show the behaviour of the material.
Four groups are defined on the block: two areas, top and bottom areas, and two node groups: centre node of the bottom plane and two nodes on the extreme positon of the y-axis. The are used to define boundary conditions and pressure (load).

 :

Loads and boundary conditions
A pressure of 6 [MPa] is applied to the top area Atop.
The boundary conditions are applied to the bottom area: all z displacements of the bottom area are restricted.
On node Nfixx the displacement in y direction is restricted.
On nodes Nfixy (two nodes on the extremes of the x axis) the displacement in x direction is restricted.

Material properties

The material properties are defined by Mooney-Rivlin.
This material behaviour is defined by a number of parameters:

• C01 = 2.3456;
• C10 = 0.709;
• C20 = 0.0;
• NU = 0.499
• K = (6*(C10+C01))/(3*(1-2*NU))

and in Code Aster called by the hyper elastic material module:

RUBBER1=DEFI_MATERIAU(ELAS_HYPER=_F(C10=C10,
                                  C01=C01,
                                  C20=C20,
                                  K=K,
                                  RHO=1000.0),);