Version

Code_Aster default
Titre : SSNP173 - Contact entre deux sphères concentriques Date : 21/04/2017 Page : 1/7
Responsable : DROUET Guillaume Clé : V6.03.173 Révision :
7c9e8dfd9b46

SSNP173 - Contact enters two concentric spheres

Summary:

This problem corresponds to a quasi-static analysis of a problem of mechanics with contact without friction
whose analytical solution is known. One is interested particularly here in integration of the terms of contact by
the method Mortar Room Average Contact (LAKE) on a curved surface. The grids are incompatible.

Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2017 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
 Version
Code_Aster default
Titre : SSNP173 - Contact entre deux sphères concentriques Date : 21/04/2017 Page : 2/7
Responsable : DROUET Guillaume Clé : V6.03.173 Révision :
7c9e8dfd9b46

1 Problem of reference
1.1 Geometry
One is consideredE sphere digs external ray 30 mm and of internal ray 20 mm , included in a hollow
sphere of external ray 40 mm and of internal ray equal to the external ray of the sphere included.

20

40

30

1.2 Properties of material

The material is elastic isotropic whose properties are:
• E=2000 MPa
•  =0.3

Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2017 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
 Version
Code_Aster default
Titre : SSNP173 - Contact entre deux sphères concentriques Date : 21/04/2017 Page : 3/7
Responsable : DROUET Guillaume Clé : V6.03.173 Révision :
7c9e8dfd9b46

1.3 Boundary conditions and loadings

Conditions of symmetry DX=0 on surface AD and DY = 0 on surface BC are forced in order
to model only one quarter of the model.

X
B C

A pressure P=25 MPa is imposed on surface AB , inside the first sphere. The two spheres are
maintained in contact with the interface.

Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2017 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
 Version
Code_Aster default
Titre : SSNP173 - Contact entre deux sphères concentriques Date : 21/04/2017 Page : 4/7
Responsable : DROUET Guillaume Clé : V6.03.173 Révision :
7c9e8dfd9b46

2 Reference solution
2.1 Method of calculating
There exists a known analytical solution with this problem, if one considers R 1 the interior ray and
R 2 the external ray, then the radial constraint expressed in polar coordinates is written:
3 3 3
R R −r
σ rr (r)= 3 1 3⋅ 2 3 ⋅P .
R 2− R1 r

Thus one finds well σ rr (R 1)=P and σ rr (R 2)=0 .

2.2 Sizes and results of reference

One tests the contact pressure on the interfaces, in r=30 mm , on both sides of discontinuity. With
R 1=20 mm and R 2=40 mm , one has then:

σ rr (r)= 4.894179894 MPa .

2.3 Uncertainty on the solution

Analytical solution.

Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2017 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
 Version
Code_Aster default
Titre : SSNP173 - Contact entre deux sphères concentriques Date : 21/04/2017 Page : 5/7
Responsable : DROUET Guillaume Clé : V6.03.173 Révision :
7c9e8dfd9b46

3 Modeling A
3.1 Characteristics of modeling
A modeling is used AXIS.

3.2 Characteristics of the grid

The grid contains 851 elements of the type TRIA6. Surfaces of contact are not in conformity.

3.3 Sizes tested and results

It is tested the contact pressure on the entirety of surfaces of main contact and slave.

Identification Type of reference Value of reference Precision

LAGS_C ‘ANALYTICAL’ 4.894179894 1%

Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2017 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
 Version
Code_Aster default
Titre : SSNP173 - Contact entre deux sphères concentriques Date : 21/04/2017 Page : 6/7
Responsable : DROUET Guillaume Clé : V6.03.173 Révision :
7c9e8dfd9b46

4 Modeling B
4.1 Characteristics of modeling
A modeling is used AXIS.

4.2 Characteristics of the grid

The grid contains 1096 elements of the type QUAD8. Surfaces of contact are not in conformity.

4.3 Sizes tested and results

It is tested the contact pressure on the entirety of surfaces of main contact and slave.

Identification Type of reference Value of reference Precision

LAGS_C ‘ANALYTICAL’ 4.894179894 1%

Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2017 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
 Version
Code_Aster default
Titre : SSNP173 - Contact entre deux sphères concentriques Date : 21/04/2017 Page : 7/7
Responsable : DROUET Guillaume Clé : V6.03.173 Révision :
7c9e8dfd9b46

5 Summary of the results

One seeks on this example very simple to test a novel method of integration of the terms of contact
based on the subdivision by subelements. This method is activated when one chooses ALGO_CONT='
LAC' for a zone of contact in formulation CONTINUOUS. One thus validates in 2D the method for the
quadratic elements, in the case of curved surfaces.

Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2017 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)