Val_004c Cook's Membrane

Problem Description

 

Cook’s membrane is a tapered panel made with rubber-like (nearly incompressible) material. One of its sides is clamped, while the opposite side is subjected to a shearing load F = 0.205 MPa for 2D case and F = 1.5 MPa for 3D case. The membrane deforms under a mixture of shear and bending strains. The plastic deformation is based on von-Mises elastoplastic model. Note that, for 2D cases, plane strain condition is assumed.

 

 

Val_004c_01

 

 

Simulation model

 

 

 

Parameter

Value

A (mm)

44

B (mm)

16

C (mm)

48

Young's modulus, E (MPa)

240.565

Poisson's ratio, ν

0.4999

Yield stress, σy (MPa)

2D: 0.45

3D: 5.0

Geometry parameters and material properties

 

 

Mesh Discretisation

 

Fourteen element types have been used in this near-incompressibility test case. The followings illustrate some examples of meshed models that have been used in 2D and 3D cases. In general, higher-order elements and elements with mixed formulation are preferred for Cook’s membrane problem.

 

Val_004c_02

 

Example of meshed models for 2D and 3D cases

 

 

Data File Descriptions

 

The data files for the different element types are in Val_004\cook. The different cases are in different sub-folders according to the dimensions (2D/3D), solver (explicit/implicit) and element type used for each specific case.

 

The basic data includes:

1Geometry_data imports geometry file for 3D cases in *.geo format that has been created using Gmsh and made compatible using ParaGeo pre-processing procedures. In 2D cases, the geometry is created directly through data file without the aforementioned pre-processing.

2Geometry_set groups multiple geometry entities under convenient geometry set names.

3Group_control_data activates geomechanical field for the current simulation group.

4Group_data sets the group name, element type, material name, porous flow type and the associated volume entity.

5Material_data defines the material properties of the model.

6Support_data constrains the displacement freedom on each surface accordingly.

7Global_loads defines the prescribed surface tangential load on the geometry set associated with the shearing load.

8Damping_global_data sets a global damping ratio of 2%.

9Control_data defines the solution algorithm (e.g. 1 for explicit transient dynamic algorithm, 7 for nonlinear implicit algorithm), termination time, etc.

 

Only key data structures are shown below.

 

Group_data

 

Group_data sets the group name, element type, material name, porous flow type and the associated volume entity.

 

Data File

 

 

* Group_data               NUM=1  

! -------------------------------

Group_name            "Volume1"

 Element_type            "HEX8M"

 Material_name     "CookPlastic"

 Porous_flow_type              0

 Volumes  IDM=1

   1

1Group name and material name are defined.

 

2Porous flow type 0 (corresponding to non-porous media) is defined.

 

3Volume set 1 is registered to be associated with the current group data structure.

 

4Element types that are used in this example include

 

Element Type

Element Descriptions

QPM4

2D 4-noded quadrilateral plane strain element

QPM4M

2D 4-noded quadrilateral plane strain mixed element

QPM8

2D 8-noded quadrilateral plane strain element

QPM8M

2D 8-noded quadrilateral plane strain mixed element

TPM3M

2D 3-noded plane strain mixed element

TPM6

2D 6-noded triangular plane strain element

TPM6M

2D 3-noded triangular plane strain mixed mixed element

HEX8_Bbar

3D 8-noded hexahedral Bbar-averaged element

HEX8M

3D 8-noded hexahedral mixed element

HEX20

3D 20-noded hexahedral element

HEX20M

3D 20-noded hexahedral mixed element

TET4M

3D 4-noded tetrahedral mixed element

TET10

3D 10-noded tetrahedral element

TET10M

3D 10-noded tetrahedral mixed element

 

 

 

Material_data

 

Material_data defines the material properties of the model.

 

Data File

 

 

* Material_data               NUM=1

! ---------------------------------

 Material_name       "CookPlastic"

 Grain_density                1000                

 Grain_stiffness             30000

 Elastic_properties          IDM=2              

  /Young's Modulus/        240.565

 /Poisson's ratio/         0.4999

 Plastic_material_type                    7

 Plastic_properties                IDM=1

 /Yield Strength/               5.0

 Hardening_type                            7

 Hardening_properties                IDM=2         JDM=2

  /Effective plastic strain/        0    0.10

  /Yield strength/                        5.0  5.52

 

1Material name is defined as "CylinderMat".

2Grain properties (density = 1000 kg/m3 and stiffness) are defined.

3Elastic properties (Young's modulus = 240.565MPa, Poisson's ratio = 0.4999) are defined.

4Plastic material type 7, corresponding to von-Mises elastoplastic model is defined.

5The yield strength of von-Mises model is defined as 5.0 MPa.

6Hardening properties are defined with prescribed hardening slope that allows visualisation of plastic strain contour at the brink of collapse.

 

 

Global_loads

 

Global_loads defines the prescribed surface normal load on the geometry set associated with the cylinder internal surface.

 

Data File

 

 

* Global_loads

! ---------------------------------

 Surface_load                     IDM=3  JDM=1

 0.0       2.0        0.0  

 Surface_load_type                   3

 Surface_orientation_vector           IDM=3  JDM=1

 /t1/        0        0        1

 Surface_load_geom_set            IDM=1

   "load_surf"

 Surface_load_geom_ass            IDM=1

   1

 

1Tangential shearing load 2.0 MPa is prescribed on the geometry set "load_surf".

 

 

 

Results

 

The load-displacement curve for 2D and 3D cases are plotted separately as the former assumes plane strain conditions. Overall, all tested element types approach asymptotically to a pressure value as the vertical displacement is driven by the shear load. The plateau in the plot is the manifestation of plastic collapse, and therefore the corresponding pressure value is treated as the limit load of the Cook’s membrane.

Val_004c_03

Load-displacement curve for Cook's membrane in plane strain condition.

 

Val_004c_04

Load-displacement curve for Cook's membrane in 3D.

 

 

The following figures show the results of Cook’s membrane deformation. Plastic strain and effective mean stress are plotted at the final time step. As shear load increases, plastic strain develops from compression-dominant region as well as extension-dominant region. These two regions can be identified from the plot of effective mean stress distribution. In 2D case, the distance between the plastic front is closer than the 3D counterpart, indicating that the 2D results are closer to plastic collapse.

 

Val_004c_05

Plastic strain contour on selected element types

 

 

 

Val_004c_06

Effective mean stress contour on selected element types

 

 

 

Val_004c_07

Deviatoric stress contour on selected element types

 

 

References

 

[1] Abaqus Benchmark Manual Version 6.11