5. C modeling#

5.1. Characteristics of modeling#

The symmetry of the problem makes it possible to represent only a quarter of the model in 3D: the sphere and the block are represented respectively by the contact surface of the sphere and a quarter of a cylinder, meshed with 3D solid elements CUB8.

A node-mesh contact is defined between the sphere and the block.

An imposed displacement load is applied to the entire surface of the sphere, which is rigidified by kinematic conditions.

_images/10000000000001720000018532F7F696FB8E16D1.png

Boundary condition:

  • Rigid body movements are suppressed by imposing a following link y between node \(E\) belonging to the sphere and node \(S\) belonging to the massif.

Loads:

An imposed displacement is applied to the entire surface representing the sphere (group of nodes \(«\mathrm{SPHSUP}»\)) in the direction \(Y\): Loading from 0 to \(–100.\mathrm{mm}\)

5.2. Characteristics of the mesh#

Number of knots: 6852

Number of meshes and type: 5326 HEXA8, 387 PENTA6 and 183 QUAD4.

5.3. Tested values#

Identification

Movations

Reference

Aster

% tolerance

Reaction \((N)\)

\(d=–20\mathrm{mm}\)

—3.201351E+06

-3.2100211E+06

1

Reaction \((N)\)

\(d=–40\mathrm{mm}\)

—6.402702E+06

-6.1671049E+06

5

Reaction \((N)\)

\(d=–60\mathrm{mm}\)

—9.604053E+06

-9.1689400E+06

5

Reaction \((N)\)

\(d=–80\mathrm{mm}\)

—1.280540E+07

-1.1738899E+07

10

Reaction \((N)\)

\(d=–100\mathrm{mm}\)

—1.600675E+07

-1.4244367E+07

12

5.4. notes#

The results are less accurate than those from 2D models. 3D meshing loses the exact character of the axisymmetric case. In addition, in order to save calculation time and memory space, 3D meshing is less refined than 2D meshing.