Code_Aster ®
Version
8.1
Titrate:
SSNA115 Arrachement of a rigid reinforcement
Date
:
25/11/05
Author (S):
J. LAVERNE Key
:
V6.01.115-A Page:
1/6
Organization (S): EDF-R & D/AMA
Handbook of Validation
V6.01 booklet: Nonlinear statics into axisymmetric
Document: V6.01.115
SSNA115 Arrachement of a rigid reinforcement with
elements with discontinuity
Summary:
This case test has as an aim the numerical study of the wrenching of a rigid reinforcement embedded in a cylinder
hollow. Decoherence is modelled starting from elements with internal discontinuity with a cohesive law CZM_EXP
(see documentation [R7.02.12]) by using modeling AXIS_ELDI. To validate the results us
will support on the analytical solution developed in [bib3]. The interested reader will also be able y
to defer for a thorough study of this case test.
Handbook of Validation
V6.01 booklet: Nonlinear statics into axisymmetric
HT-66/05/005/A
Code_Aster ®
Version
8.1
Titrate:
SSNA115 Arrachement of a rigid reinforcement
Date
:
25/11/05
Author (S):
J. LAVERNE Key
:
V6.01.115-A Page:
2/6
1
Problem of reference
1.1 Geometry
and
loading
Either a hollow roll length L, interior radius RF and external radius R. Soit one
brace rigid circular section of radius RF embedded in its center. One notes I and E them
surfaces interior and external of the hollow roll (see [Figure 1.1-a]). The loading consists with
to apply, at the node of the rigid reinforcement, a displacement I
U.E.
I
Z (U > 0) as well as a displacement
no one on the external edge E.
=
U
With
R
0
I
U = U.E.
Z
L
I
E
E
Z
ur = 0
E
With
R
R F
R
Appear 1.1-a: Schéma of the field and loading
The assumption of an axisymmetric solution is made what enables us to restrict our study has one
rectangular field 2D. Dimensions of the field are as follows
:
RF = 0.5mm, R = 5.5mm, L = 10mm. the loading on the rigid reinforcement will be taken into account
by applying displacement imposed I
U ez on all side I of the field 2D like one
null displacement on the side E to take into account the embedding of the cylinder. Finally one imposes
a radial displacement no one on the faces lower and higher of the field in order to avoid a singularity
dependant on a change of boundary condition at points A and A (see [Figure 1.1-a]). These conditions
in extreme cases will lead to a anti-plane solution (independent of Z) what makes it possible to obtain more
simply an analytical solution.
1.2 Parameters
Material
Values of the Young modulus, the Poisson's ratio, the critical stress and the tenacity of
material are taken in the following way:
- 1
E = 1.5 MPa, = 0, C = 1.1 Mpa, C
G = 0.9 N.mm
(They are of course values “tests” which does not correspond to any material in particular.)
Handbook of Validation
V6.01 booklet: Nonlinear statics into axisymmetric
HT-66/05/005/A
Code_Aster ®
Version
8.1
Titrate:
SSNA115 Arrachement of a rigid reinforcement
Date
:
25/11/05
Author (S):
J. LAVERNE Key
:
V6.01.115-A Page:
3/6
2
Reference solution
The reference solution is an analytical solution draw from [bib3], it even inspired from a study
unidimensional proposed in [bib1] and in a more general way being based on the approach
energetics of the rupture suggested per G.A. Francfort and J.J. Marigo [bib2]. We will not return
in the details of the calculation of this solution, we will present just the analytical value of the answer
total of the structure: displacement imposed U according to the corresponding force F:
()
F L
U F =
+ sign (F)
1
F
-
éq
2-1
2 R µ
F L
2 R
F L
where µ indicates the coefficient of Lamé (µ = E 2 here), the density of energy of cracking (see
documentation [R7.02.12]) and where L = RF ln (R RF) is a length structural feature
decisive for the brutal or progressive evolution of decoherence.
3 Modeling
With
3.1
Characteristics of modeling
Simulation is carried out into axisymmetric. The elements with internal discontinuity allow
to represent the fissure along I. The latter have as a modeling AXIS_ELDI and one
cohesive behavior CZM_EXP. The other elements of the grid are QUAD4 with one
elastic behavior ELAS in modeling AXIS.
3.2
Characteristics of the grid
One carries out a grid structured in quadrangles of the field with 76 meshs in the height and
28 meshs in the radial direction. One lays out a layer of elements with discontinuity interns length
of I using command CREA_MAILLAGE and key word CREA_FISS (see documentation
[U4.23.02]). The orientation of the elements with discontinuity is carried out so that the direction
normal is directed according to - er (the tangential direction is thus according to - ez). The remainder of the field
is divided into linear meshs QUAD4 (see [Figure 3.2-a]).
ez
I
E
E
R
Elements with discontinuity
QUAD4
Appear 3.2-a: Maillage of the field
Handbook of Validation
V6.01 booklet: Nonlinear statics into axisymmetric
HT-66/05/005/A
Code_Aster ®
Version
8.1
Titrate:
SSNA115 Arrachement of a rigid reinforcement
Date
:
25/11/05
Author (S):
J. LAVERNE Key
:
V6.01.115-A Page:
4/6
3.3 Functionalities
tested
Commands
STAT_NON_LINE COMP_INCR
RELATION
CZM_EXP
AFFE_MODELE MODELING AXIS_ELDI
DEFI_MATERIAU RUPT_FRAG
SIGM_C
SAUT_C
CREA_MAILLAGE CREA_FISS
3.4
Sizes tested and results
The tangential constraint T along the fissure (i.e in the elements with discontinuity) corresponds to
opposite of the force F divided by the surface of decoherence: 2 R F L. Moreover while resting on
form density of energy of surface defined in [R7.02.12] and according to [éq 2-1] one deduces
following relation:
()
L
T
= -
+ sign () C
G
U
T
T
ln
T
µ
éq
3.4-1
C
C
The latter will enable us to carry out tests summarized in the table below.
Size tested
Theory
Code_Aster
Difference (%)
Tangential constraint: VI7
7.69747E-01
7.6974726277784E-01
3.41E-05
PG1 of mesh MJ38
Moment: 6.00070E+00
Tangential constraint: VI7
4.34935E-01
4.3493490987033E-01
- 2.07E-05
PG1 of mesh MJ38
Moment: 1.20004E+01
Tangential constraint: VI7
1.28483E-01
1.2848319446210E-01
1.51E-04
PG1 of mesh MJ38
Moment: 1.93334E+01
Displacement DY
1.57674E+00
1.5767415306566E+00
9.71E-05
N5 node
Moment: 1.20004E+01
4
Summary of the results
It is noted that the element with discontinuity allows a good prediction of decoherence, indeed this
last develops in an identical way on all the height of the cylinder moreover the results
numerical are very close to the analytical solution. In addition modeling suggested allows
to correctly reproduce the brutal or progressive evolution of cracking according to
lengths characteristic L of the structure and C
G C of the behavior. This last point is not
highlighted here but the interested reader will be able to refer to [bib3] for more details.
Handbook of Validation
V6.01 booklet: Nonlinear statics into axisymmetric
HT-66/05/005/A
Code_Aster ®
Version
8.1
Titrate:
SSNA115 Arrachement of a rigid reinforcement
Date
:
25/11/05
Author (S):
J. LAVERNE Key
:
V6.01.115-A Page:
5/6
5 Bibliography
[1]
CHARLOTTE Mr., FRANKFURT G.A., MARIGO J.J. and TRUSKINOVSKY L.: Revisiting
brittle fracture have year energy minimization problem: comparison off Griffith and Barenblatt
surface energy models. Proceedings off the Symposium one “Continuous Damage and
Fracture " The dated science library, Elsevier, edited by A. BENALLAL, Paris, pp. 7-18, (2000).
[2]
FRANKFURT G.A. and MARIGO J.J.: Revisiting brittle fracture have year energy minimization
problem. J. Mech. Phys. Solids, 46 (8), pp. 1319-1342 (1998).
[3]
LAVERNE J.: Energy formulation of the rupture by models of cohesive forces:
numerical considerations theoretical and establishments, Thèse de Doctorat of Université
Paris 13, November 2004.
Handbook of Validation
V6.01 booklet: Nonlinear statics into axisymmetric
HT-66/05/005/A
Code_Aster ®
Version
8.1
Titrate:
SSNA115 Arrachement of a rigid reinforcement
Date
:
25/11/05
Author (S):
J. LAVERNE Key
:
V6.01.115-A Page:
6/6
Intentionally white left page.
Handbook of Validation
V6.01 booklet: Nonlinear statics into axisymmetric
HT-66/05/005/A