Code_Aster ®
Version
4.0
Titrate:
Calculation of load limits by the method of Norton-Hoff-Friaâ
Date:
12/11/97
Author (S):
F. VOLDOIRE, E. SCREWS
Key:
R7.07.01-A
Page:
1/16
Organization (S): EDF/IMA/MN
Handbook of Référence
R7.07 booklet: Analyze limit
R7.07.01 document
Calculation of load limits by the method
of Norton-Hoff-Friaâ
Summary:
The limiting analysis makes it possible to determine the acceptable loadings D `a structure, of geometry fixes given,
composed of a material having a criterion of resistance. One considers the case of loadings made up of
summon of a continued load and of another parameterized by the load factor.
After a recall of the theoretical formulation, one presents the regularized kinematic approach applied to the criterion
of resistance of Von Mises (method of Norton-Hoff-Friaâ) and implemented in Code_Aster. One will be able
to refer to [bib4] for the various possible methods of regularization suggested in the literature. One
expose then the calculation of the solutions of this nonlinear problem and postprocessing providing one
estimate by higher value of the limiting load.
Handbook of Référence
R7.07 booklet: Analyze limit
HI-74/97/023/A
Code_Aster ®
Version
4.0
Titrate:
Calculation of load limits by the method of Norton-Hoff-Friaâ
Date:
12/11/97
Author (S):
F. VOLDOIRE, E. SCREWS
Key:
R7.07.01-A
Page:
2/16
Contents
1 theoretical Formulation of the analysis limits ............................................................................................... 3
1.1 Definition of the load limits ........................................................................................................... 3
1.2 Calculation of the load limits by a kinematic approach ............................................................... 3
1.3 Regularization of the kinematic approach by the method of Norton-Hoff-Friaâ ........................... 5
2 numerical Aspects of the calculation of the load limits ................................................................................. 7
2.1 Establishment in Code_Aster .................................................................................................... 7
2.2 Relation of behavior of Norton-Hoff in STAT_NON_LINE ............................................... 8
2.3 Postprocessing of the calculation of the load limits ................................................................................. 10
3 an example of validation ..................................................................................................................... 12
3.1 Problem of reference .................................................................................................................. 12
3.2 Plane case ........................................................................................................................................ 12
3.2.1 Solution analyzes limit of it ..................................................................................................... 12
3.2.2 Solution analyzes regularized limit of it .................................................................................. 13
3.3 Axisymmetric case ......................................................................................................................... 13
3.3.1 Solution analyzes limit of it ..................................................................................................... 14
3.3.2 Solution analyzes regularized limit of it .................................................................................. 14
3.4 Three-dimensional case ....................................................................................................................... 15
3.4.1 Solution analyzes limit of it ..................................................................................................... 15
3.4.2 Solution analyzes regularized limit of it .................................................................................. 15
4 Bibliography ........................................................................................................................................ 16
Handbook of Référence
R7.07 booklet: Analyze limit
HI-74/97/023/A
Code_Aster ®
Version
4.0
Titrate:
Calculation of load limits by the method of Norton-Hoff-Friaâ
Date:
12/11/97
Author (S):
F. VOLDOIRE, E. SCREWS
Key:
R7.07.01-A
Page:
3/16
1
Theoretical formulation of the limiting analysis
1.1
Definition of the limiting load
One considers a solid occupying a limited field subjected to surface loadings
F + F on the edge and of the loadings of volume F + F on. The loading is distinguished
0
F
0
(F, F), parameterized by positive reality, and the permanent loading (F, F). Conditions of
0
0
Dirichlet homogeneous are applied to the complementary edge U of (an imposed displacement
or an initial anelastic deformation - thermics, plastic… - do not have an effect on the field of
working loads). One can find in [bib5] several other useful properties.
The material constitutive of the solid has a criterion of resistance expressed by a scalar function of
constraints, negative for working stresses. The criterion used for a material of the type
elastoplastic perfect with threshold of von Mises and selected here is:
(
2
2
2
2
G) = J () -
3
D
D
y =
.
-
2
y =
. (-) + (2 - 3) + (1 - 3) - y
2
1
2
D
is the diverter of the tensor of the constraints,
is the threshold of resistance in simple traction (like an elastic limit), possibly
y
variable according to zones' of the solid considered.
being principal constraints of.
I
Being given this criterion of resistance one seeks to calculate the value limits, called limiting load
, for which the structure can support the loadings F + F and F + F.
lim
lim
0
lim
0
Strictly speaking, the value indicates the limit of the bearable loadings, but for
lim
materials obeying Principe of Travail Plastique Maximal, this value is the limit of
supported loadings.
1.2
Calculation of the load limits by a kinematic approach
In design the collapse two approaches are possible: static approach (in variables of constraints)
and kinematic approach (in variables speeds). These approaches provide terminals of
charge limit: undervaluing for the approach static and raising for the kinematic approach. When them
two provide the same result, the limiting load obtained is exact.
The kinematic approach is that used in Code_Aster using finite elements in displacements.
For the loading given (F, F), one defines the space speeds kinematically acceptable and
standardized by:
V 1 = v acceptable, v = 0 on, L
has
U
(v) =
F. v D + F. v ds =
1
F
The power of the permanent loading (F, F) is noted:
(v).
0
0
L0
Handbook of Référence
R7.07 booklet: Analyze limit
HI-74/97/023/A
Code_Aster ®
Version
4.0
Titrate:
Calculation of load limits by the method of Norton-Hoff-Friaâ
Date:
12/11/97
Author (S):
F. VOLDOIRE, E. SCREWS
Key:
R7.07.01-A
Page:
4/16
From the criterion of resistance in constraints G
(), one defines:
· the whole of working stresses by: G
= X, G X 0
X
{() (()) }
()
((
G is convex for the criterion G)
X)
0, if (X) (
G X)
· the indicating function: X =
G (
()) +, if (X)
(
G X)
· the function of support: () = Sup [. - G
(])
R6
Sup in () can be reached only if is selected in
D
(
G, such as: = + µId
X)
/
1 2
(what ensures//D). The optimum corresponds to (
G) =
=
-
0
2. D D
y
3 (. )
D2
()
G
· D
· D
0
D
1
0
Appear 1.2-Error! Argument of unknown switch. : Optimum and graph of the function () in 1D
2
From where the function of support: ((v)) =
.
y
(v).(v) +
Sup (µ.div v). It is observed that
3
µR
function () is not differentiable into 0.
One to date does not treat in Code_Aster possible internal surfaces of discontinuity with the center
solid [bib 4].
The kinematic approach is defined using the convex functional calculus S v, positively homogeneous
E ()
of degree one, for v
1
V
defined on the whole field:
has
S (v) = ((v))D - L (v)
E
0
This functional calculus is the integral on the field of the function of support of convex (
G, calculated
X)
in (v) and is interpreted like the maximum resistant power in the field speeds v (
contribution of resistance of interface on surfaces of discontinuity is supposed to be null). The function
support is positively homogeneous degree 1, and thus the functional calculus (v) also by
consequence.
Handbook of Référence
R7.07 booklet: Analyze limit
HI-74/97/023/A
Code_Aster ®
Version
4.0
Titrate:
Calculation of load limits by the method of Norton-Hoff-Friaâ
Date:
12/11/97
Author (S):
F. VOLDOIRE, E. SCREWS
Key:
R7.07.01-A
Page:
5/16
With the criterion of Von Mises the functional calculus of power S v is:
E ()
2
S (v) =.
(v.)(v) +Sup (q.div v) D - L
E
y
(v)
éq 1.2-1
3
0
Q
R
where it is noted that only fields v pertaining to C = {v V 1a, div v = 0 in}
provide finished values. Fields v must thus check the condition known as of incompressibility
div v = tr (v) = 0. This is why it is necessary to use the incompressible elements for a calculation of
charge limit with the criterion of Von Mises. The processing of the incompressibility is detailed in
[R3.06.05].
The load limits given by the kinematic approach is:
lim
S v
E
lim Inf
=
(V)
()
=
Inf
S v
L v
E
1
v 1
v
V
V
L (v)
=
Sup Inf (() - (() -)
1
has
>0 vV
has
has
L (v) >0
With the optimum one obtains a solution U and the limiting load (not unicity of U but unicity of).
lim
lim
Thus, any loading L (v) + L
(v) with 0 is bearable. Beyond, it
0
lim
lim
problem of balance does not have a solution.
Note:
There are situations where, even if
L (v) is not bearable alone,
0
combination L (v) + L
(v), for, becomes it on a certain interval, and not
1
0
2
only for two parallel loadings.
Note:
The limiting load calculated for a two-dimensional problem, in plane deformations, is
necessarily higher than that obtained for this problem modelled in plane constraints. It
result thus provides one raising. If one wishes to deal with the problem in plane constraints, it is
necessary then to make the kinematic approach on a three-dimensional modeling.
1.3 Regularization of the kinematic approach by the method of
Norton-Hoff-Friaâ
Implementation the numerical of the kinematic approach requires the minimization of the functional calculus
not-differentiable S v. Many techniques of regularization exist [bib4]. Method of
E ()
Norton-Hoff-Friâa is used here [bib2], [bib7]. It rests on precursory work of Casciaro in
1971. It consists in replacing the function of support
() by the function of support regularized and
differentiable NH (). It is adjustable by a parameter of regularization m (1 m2), of which
value limits m +
1 conduit with convergence towards the function of support
():
-
1 m
NH ()
K
=
(())m
m
Handbook of Référence
R7.07 booklet: Analyze limit
HI-74/97/023/A
Code_Aster ®
Version
4.0
Titrate:
Calculation of load limits by the method of Norton-Hoff-Friaâ
Date:
12/11/97
Author (S):
F. VOLDOIRE, E. SCREWS
Key:
R7.07.01-A
Page:
6/16
2
1
In Code_Aster one uses the parameter of regularization N with N =
and K is chosen
y
=
(M-1)
µ
3
to find the incompressible elastic problem when m = 2 is N =1 (2µ being the coefficient of
rigidity of shearing).
One notes the space acceptable speeds adapted to the problem of viscous flow for the law of
Norton-Hoff of command m:
V m 1 = {
m
v L (), and (v)
m
L
(), v = 0 on, L
has
U
(v) = 1}
One defines on this space the regularized functional calculus Sm
E (v):
-
1 m
K
m
Sm (v) =
((v)) D - L
E
(v)
m
0
The problem of Inf minimization
m
is well posed thanks to the properties of
m 1 [(v)]
vV A
Lm spaces () and has a solution um, for which the value of Inf reached is: Mr. One shows
that this problem can be also written in the form of the search of the point-saddle (m, um, pm)
Lagrangian following:
(
With)
m
m
Max Inf
Sup
.((v).(v)) D
.
Q div v D L (v) (L (v))
1
0
V
2
m
Q has
+
-
-
-
R v
L ()
éq 1.3-1
m/2
/2
1
2
-
with: (
With)
m
K - m (
2
1 2
3µ
. In practice the continuation is taken:
3)
m
- m
m
m
=
y
= y
() (3)
m
I
= + -
1 10
2
1,1
1,01
1,001
…
1
1
1
10
100
1000
…
N = M-1
(
With)
m
2µ
…
2
y
3
It is noticed that (
With m) is increasing with m (if E y) and homogeneous with a constraint, and remains
m/2
limited when N +. If one chooses E = y then (
With)
m = 2
y (
. This Lagrangian allows
3)
to impose directly in the operator the condition of incompressibility and standardization on 1 of
power of the loadings. One then builds a decreasing succession of m and the load limits lim
is the limit of this continuation when m +
1 (either N +):
m
m
lim Li
= m
Inf
S
v
S
U
éq 1.3-2
E
E
m
m 1
m1 [
()] lim
=
(()
m 1
v Goes
Handbook of Référence
R7.07 booklet: Analyze limit
HI-74/97/023/A
Code_Aster ®
Version
4.0
Titrate:
Calculation of load limits by the method of Norton-Hoff-Friaâ
Date:
12/11/97
Author (S):
F. VOLDOIRE, E. SCREWS
Key:
R7.07.01-A
Page:
7/16
For the demonstration one will refer to [bib4] and [bib7]. One shows also the following property of
solutions of [éq 1.3-1]. Thus if one ampifie the loading L L
, when L0 = 0, the solutions have
following dependence according to:
U
- 1
-
1 m
m () =
um ()
1
;
pm () =
pm ()
1
D
(U
-
1 m D
m
- m m
m () =
(um () 1);
E
S (um () =
E
S (um ()
1)
One of the essential assets of this method of regularization resides in the properties of fitment
Lm spaces (), which makes it possible to have properties interesting for the continuation of the m, to see it
paragraph 2.3. Thus one shows the following properties [bib7], for a limited field, in
calling =
D and = (A)
m D
:
m
For all 1 m and 1 R S and any function U of V m
has, one a:
1
1
1
1
-
1
-
1
(
With)
m
(U). (U) D
R
((A) m ((U). (U) R)/2D
2
) R
S
S/
S
(U).
m (
(
With)
m (
(U)
)
D
m
)
éq 1.3-3
1
1
m
m
2
1
m
2
y.
(U).(U) D m
y
.
(U).(U)
D
U
éq 1.3-4
3
3
These properties are interesting because applicable if the material is heterogeneous, and one can
to consider the limit of resistance either as measures (density) or like belonging to
the deformation energy.
2
Numerical aspects of the calculation of the limiting load
2.1
Establishment in Code_Aster
To carry out a calculation in Code_Aster in analysis limits with the method of regularization of
Norton-Hoff-Friaâ with the criterion of resistance of Von Mises, it is necessary:
· to define the model 2D (plane or axis) or 3D with the incompressible finite elements,
· to define the characteristics of the materials (Young modulus E y, Poisson's ratio
near to 0 5
, to ensure the operation of the incompressible finite elements [R3.06.03],
limit of resistance y and coefficient of Norton-Hoff N), the coefficient of shearing is
deduced: 2µ = 2nd/3. It should be noted that the limiting load is independent of E and,
· to define the permanent loading and that which is parameterized by,
· to define the standardization of the power of the loading parameterized (key word LIAISON_CHAMNO
command AFFE_CHAR_MECA),
· to carry out a non-linear calculation with the relation of behavior of Norton-Hoff with
order STAT_NON_LINE [U4.32.01],
· post-to treat calculation to obtain the load limits with command POST_ELEM [U4.61.04].
Handbook of Référence
R7.07 booklet: Analyze limit
HI-74/97/023/A
Code_Aster ®
Version
4.0
Titrate:
Calculation of load limits by the method of Norton-Hoff-Friaâ
Date:
12/11/97
Author (S):
F. VOLDOIRE, E. SCREWS
Key:
R7.07.01-A
Page:
8/16
2.2
Relation of behavior of Norton-Hoff in STAT_NON_LINE
The problem is written in variational form in the following way:
For N = 1/(m -)
1 given, to find (m, um, m
p) R ×V
2
× L has () such as:
(
m2
With m).
((um).(um)) (um). v +
v
()
D
p. div D
m
- ml (v) = L (v) vVa
0
q.divu D
= 0 Q
2
m
L ()
L (U
m) = 1
éq 2.2-1
This problem admits a single solution for any N 1 (see [bib4]). For m = 2 or N =1 the problem
is of incompressible linear elasticity type.
One thus obtains an estimate of the limiting load by higher value, the field um giving one
idea of a mode of ruin.
For the processing of the incompressibility, one will refer to the document [R3.06.05]. The treated equation is
in fact:
Q. p
Q. p
q.div U D +
D =, Q
L ()
D
allowing to avoid putting
0
2
, the term
at fault the solvor employed and correspondent with a choice of the Poisson's ratio = 0 4999
.
…
solutions are thus only quasi-incompressible.
The principle of operation of the general algorithm of STAT_NON_LINE is described precisely in
[R5.03.01]. One leads to the following incremental problem:
To find (, U, p) R × V × Q
0
such as:
(U + U
).(v) +
(
B,
v
+) -
D
p
p
(+) L1 (v) = (L + L
0
0) (v)
v V
0
(
B U +,
U Q) = D + D
Q Q
L1 (v) = 1
·
B is a linear operator who contains the boundary conditions homogeneous of Dirichlet,
incompressibility,
·
D describes the data imposed on the solution (boundary conditions of Dirichlet, incompressibility),
·
L is the permanent loading and L the loading controlled by the parameter,
0
1
·
V and
0
Q are spaces of functions discretized on the basis of finite element, and are thus
defined by a vector (U, P) of degrees of freedom.
Handbook of Référence
R7.07 booklet: Analyze limit
HI-74/97/023/A
Code_Aster ®
Version
4.0
Titrate:
Calculation of load limits by the method of Norton-Hoff-Friaâ
Date:
12/11/97
Author (S):
F. VOLDOIRE, E. SCREWS
Key:
R7.07.01-A
Page:
9/16
The tensor of the constraints (U) checks the relation of behavior of Norton-Hoff. The diverter of
constraints associated at the speed of deformation is:
m2
-
1 N
D (U) = ().(D (U.)D (U)) .D (U) D (U) = () N.(D.D) .D
With m
With m
(U)
éq 2.2-2
The problem is solved by the method of Newton, after direct implicit discretization of the relations
of behavior [R5.03.02].
The phase of prediction consists in solving the following system, starting from the current state (U, p) for
to obtain the reiterated first: (U, p.
0
)
0
D
.(U
).
0 (v) D
(
B,
v p
)
0
1
L (v)
L (v)
v V
+
-
=
D
O
0
(U)
(
B U
, Q) = D
0
Q
Q
1
L (U + U
0)
=
1
D
Option RIGI_MECA_TANG provides the tangent operator
.
D (U)
One uses the tangent rigidity, applied to a tensor deviatoric E:
D D
- 2
(U)
(U)
.e = (
R
D
D
With R)
D
D
. (U). (U)
. E + R
- 2
.e
D D
(
)
(
) D
D
(U). (U)
Then one treats the phase of correction, for iteration I:
D
.
(U
- U
+1).(v) D
I
I
+ B (v, p
i+1) - (
+) 1
L (v) = L
+ L (v)
D
(0 0)
(u+ui)
-
[(u+ U
).(v) D
I
+ B (v, p)] v
Vo
B
(U
, Q
i+1
) = D+ D
- B (,
U Q)
Q Q
1
L (u+ U
i+1) = 1
Option RAPH_MECA provides the second member
D
(U
U
+).
. This one is built by call
I
(v)
with the subroutine of Code_Aster of the type NICOMP, which calculates the stress field suitable for
D
law of behavior used. Option FULL_MECA provides in more the tangent operator
. One
D (u+u) I
can only decide updating of this operator with certain iterations I, to avoid one
too frequent expensive assembly.
Handbook of Référence
R7.07 booklet: Analyze limit
HI-74/97/023/A
Code_Aster ®
Version
4.0
Titrate:
Calculation of load limits by the method of Norton-Hoff-Friaâ
Date:
12/11/97
Author (S):
F. VOLDOIRE, E. SCREWS
Key:
R7.07.01-A
Page:
10/16
In our case, the resolution can be made in a nonincremental direct way, but there is interest with
to reactualize the tangent rigidity of time to other to accelerate convergence. The solvor to be used
is: “LDLT”, while having taken care to have specified that standardization is placed at the end of the system
(key word: NUME_LAGR: “APRES” in LIAISON_CHAMNO).
One can also decide to make recoveries, starting from a solution (U, p) obtained previously,
even coming from a resolution for another parameter N (calculations less expensive and better
convergence).
In all the cases, it is advised to start to make a calculation on a coarse grid for
to evaluate the effect of parameter N on Mr.
2.3
Postprocessing of the calculation of the limiting load
Having obtained the solution (m, um, m
p), for a N given, it remains to use the continuation of the m for
to build the approximation of the limiting load. For that one exploits the properties [éq 1.3-2], [éq 1.3-3],
the fact that (
With)
m is increasing and the property resulting from minimization [éq 1.3-1] (see [bib7]).
From these two last, with 1 R S, one deduces that for and ur and custom respective solutions (checking
also the condition of incompressibility and standardization) of [éq 1.3-1] for m = R and m = S:
((Ar)
(
/2
/2
(U). (U) R
)
D) ((
With S)
(
S
(custom).(custom)) D
R
R
)
Associated the property [éq 1.3-2], one draws for 1 R S:
1
1
1
1
-
1
-
1
(
With R)
/2
/2
(U). (U) D
R
((Ar)
(
R
R
S
S
S
(ur).(ur)) D
)
(U).
S
((Ace)
(S (custom)) D
R
R
R
)
éq 2.3-1
~
One notes the terms m of the continuation, which one calculates in practice by postprocessing using um (
external power being unit:
1
1
~
-
1 m
m2
m
m =
U. U
-
m
((A)
m (
/
(
)
m (m)
)
D
L U
éq 2.3-2
) () m
0
~
This continuation m is thus decreasing for N + and it is shown [bib7] that it converges towards lim,
what allows a good control. As one can undervalue (knowing that (
With)
1
2
= y
) the first term
3
of [éq 2.3-1]:
2
~
lim y
(U)
Mr. (U)
m D -
L (U)
m
m
3
0
Handbook of Référence
R7.07 booklet: Analyze limit
HI-74/97/023/A
Code_Aster ®
Version
4.0
Titrate:
Calculation of load limits by the method of Norton-Hoff-Friaâ
Date:
12/11/97
Author (S):
F. VOLDOIRE, E. SCREWS
Key:
R7.07.01-A
Page:
11/16
the continuation is thus calculated!
m decreasing for N + (parameter “CHAR_LIMI_SUP” of the table
exit of POST_ELEM [U4.61.04]) and also converging towards lim,
!
2
lim
m = y (U) Mr. (U) Mandelevium -
L (U)
m
éq 2.3-3
3
0
One judges quality of the approximation of the limiting load by comparison of different
lim
values of!
m which converge towards by excess (out of N +). These terms are calculated by
lim
numerical integration at the points of Gauss of the finite elements.
Another interpretation of the interest which this continuation brings lies in the fact that it exploits
directly the expression of the function of support of convex of resistance, i.e. power
dissipable in the modes of potential ruin, applied to the incompressible and standardized solutions
calculated um.
If the permanent loading (parameter “PUIS_PERMANENTE”) is null: L = 0, one can
0
easily to exploit the stress field (almost statically acceptable) calculated with
solution um and to obtain a value by estimate of the limiting load, which would be him would necessarily be
lower if balance were checked exactly (see [bib4]). One thus calculates the continuation m (parameter
“CHAR_LIMI_ESTIMEE” of the table resulting from POST_ELEM [U4.61.04]), which does not have on the other hand
properties of monotony:
1
3
-
(
With)
m
U.
U
2
=
.
U. U
D. Sup
!
éq 2.3-4
m (
D
D
() ())m
(m) (m)
m
m
m
X
m
y
This maximization (of the function called gauge of convex of resistance) is calculated only with
points of Gauss of the finite elements. Also the value obtained, for each m, lower than!
m [bib4],
can be regarded only as one indication.
Handbook of Référence
R7.07 booklet: Analyze limit
HI-74/97/023/A
Code_Aster ®
Version
4.0
Titrate:
Calculation of load limits by the method of Norton-Hoff-Friaâ
Date:
12/11/97
Author (S):
F. VOLDOIRE, E. SCREWS
Key:
R7.07.01-A
Page:
12/16
3
An example of validation
3.1
Problem of reference
A rectangular plate or a hexahedron or an axisymmetric cylinder are considered. The criterion of
resistance of homogeneous constitutive material checks the criterion of von Mises (with for threshold y).
structure is subjected to pressures on the horizontal edge - F and the vertical edge
- (1 -) F with 0.5. On this very simple problem, an analytical calculation makes it possible to obtain
charge exact limit in the direction of the loading, as well as the estimates produced by
method of regularization. For more details one will refer to [bib4] and [bib5]. This example of
validation corresponds to test SSNV124 [V6.04.124].
The geometry is defined by:
· Interior radius: = m has
1 m,
· External radius: B = 2mm, thickness B - has = m
1 m,
· Height
:
H = 4mm.
Z or y
D
C
H
B
has
R or X
0
With
B
3.2 Case
plan
The structure is subjected to pressures on the edges horizontal: - F and vertical: - 1
(-) F, with:
1/2, and a blocking in Z is exerted. One considers two ways of controlling the loading:
· cas1: the two pressures horizontal and vertical are parameterized by,
· cas2: the horizontal pressure is parameterized by, while the vertical pressure is
constant - (1 -) f0, with F = F
0
0
.
3.2.1 Solution analyzes limit of it
The solution is homogeneous (biaxées constraints:
(
)
xx = F, yy = 1 - F, xy = 0,
plane deformations). One obtains [bib4] the limiting load in these directions of loading, for
criterion of von Mises, in plane deformations, with the threshold y:
2
3 y
cas1: lim. F =
éq 3.2.1-1
.
3
2 - 1
2
3 y 1 -
cas2: lim. F =
+
. F
.
3
0
éq 3.2.1-2
Handbook of Référence
R7.07 booklet: Analyze limit
HI-74/97/023/A
Code_Aster ®
Version
4.0
Titrate:
Calculation of load limits by the method of Norton-Hoff-Friaâ
Date:
12/11/97
Author (S):
F. VOLDOIRE, E. SCREWS
Key:
R7.07.01-A
Page:
13/16
It is checked that if one takes 0 = lim in the cas2, one finds the cas1 then.
3.2.2 Solution analyzes regularized limit of it
The solution is homogeneous. The plane deformations are necessarily of the form:
1 0
0
(U) =
0 - 1 0
;
(U).(U) 2
éq 3.2.2-1
0 0
=
0
By the law of Norton-Hoff, one obtains the deviatoric constraints:
1 0
0
D = () m2 m2
- 2
-
2
1
. 0 - 1 0;
D
() m
2
m
With m
With m
=
3
éq 3.2.2-2
VM
0 0
0
The standardization of the loading leads to:
1
cas1: F =
éq 3.2.2-3
H (B - has) (
2 -)
1
1
cas2: F =
éq 3.2.2-4
H (B - has)
Terms of the continuation!
m of limiting load in these two parameter settings of the loading is then:
2
3 y
cas1: !
.
m F =
m
éq 3.2.2-5
.
3
2 - 1
2
3 y 1 -
cas2: !
.
m F =
+
. F
m
0
éq 3.2.2-6
.
3
Invariance according to m observed here (what is a particular case) results owing to the fact that one is
in an isostatic situation. In the cas1, one can also exploit the continuation of the m:
2
3 y
cas1: . F =
m
éq 3.2.2-7
m
3.
2 - 1
One thus obtains the load limits lim exact when m +
1.
3.3 Case
axisymmetric
In axisymmetric 2D one considers the same geometry, but the solid, on which one imposes a bloquage
axial complete, is only subjected to a pressure on the internal wall: F parameterized par.
Handbook of Référence
R7.07 booklet: Analyze limit
HI-74/97/023/A
Code_Aster ®
Version
4.0
Titrate:
Calculation of load limits by the method of Norton-Hoff-Friaâ
Date:
12/11/97
Author (S):
F. VOLDOIRE, E. SCREWS
Key:
R7.07.01-A
Page:
14/16
3.3.1 Solution analyzes limit of it
One obtains [bib5] the limiting load in this direction of loading, for the criterion of von Mises, in
axisymmetric and null axial deformations, with the threshold y:
2 3
B
lim
. F =
ln
y
éq 3.3.1-1
3
has
3.3.2 Solution analyzes regularized limit of it
The solution is homogeneous. Displacement being only radial, the isochoric deformations are
necessarily of the form:
- 1 0
0
U
R (R) =
;
(U) =
0
0 0
;
(U).(U)
2
éq 3.3.2-1
R
R 2
= R
2
0
0
1
By the law of Norton-Hoff, one obtains the deviatoric constraints:
- 1 0
0
D = () m2 m2 - 2m+2
2
1
2
2
.
. 0 0 0
; D
() m
2
m - m
With m
R
With m
.r
+
2
=
. 3 éq 3.3.2-2
VM
0 0
1
The equilibrium equations axial and radial result in determining the average constraint:
m2
m2 - 2 + 2
2
- m
tr (R) 3 (
With)
m2
.
.r m
=
.
+
3
éq 3.3.2-3
1 - m
where is a constant, which is calculated starting from the boundary condition of null pressure in wall
external. The components of the constraints then are obtained:
- 2m+2
- 2m+2
rr (R) = (B
- R
)
m
- 2
2
- 2m+2
- 2m
With m
+2
zz (R) = (B
- (2 -)
m R
)
()
with: =
éq 3.3.2-4
m
- 1
- 2m+2
- 2m
m
2
- 1 fH
+
(R) =
(B
- (3 - 2)
m R
)
(
) ()
1
The standardization of the loading leads to:
F =
.
H
Terms of the continuation!
m of limiting load for this loading is then:
!
2 3
B
2 3
B
.
m F =
y H
rdr =
ln
m
éq 3.3.2-5
3
has 2
y
R
3
has
The terms of the continuation m of limiting load for this loading are:
- 1
- 2 +2
- 2 +
2 3
B
2
2
2
2
m
m
3 B
-
-
+
- has
. F =
m
m
y
R
rdr.
Max R
m
y
=
éq 3.3.2-6
m
3
has
(
)
] has, [B
3 (
m 1 -)
- 2m
m
has
+2
2 3
B
In m +
1, one finds: +
. F =
ln
1
3
y
, i.e. the same value as!
has
m and lim.
Handbook of Référence
R7.07 booklet: Analyze limit
HI-74/97/023/A
Code_Aster ®
Version
4.0
Titrate:
Calculation of load limits by the method of Norton-Hoff-Friaâ
Date:
12/11/97
Author (S):
F. VOLDOIRE, E. SCREWS
Key:
R7.07.01-A
Page:
15/16
3.4 Case
three-dimensional
In 3D one considers the same geometry, but the solid, unit thickness, is free in the direction
antiplane Z. The solid is subjected to pressures on the walls horizontal: - F and vertical:
- 1
(-) F, with: 1/2. The two pressures horizontal and vertical are parameterized par.
3.4.1 Solution analyzes limit of it
The solution is homogeneous (biaxées constraints:
(
)
xx = F, yy = 1 - F, xy = 0, zz = 0,
deformations). One obtains the limiting load in this direction of loading [bib5], for the criterion
of von Mises, with the threshold y:
y
lim. F =
éq 3.4.1-1
3 2 -
3 + 1
3.4.2 Solution analyzes regularized limit of it
The solution is homogeneous. The isochoric deformations are necessarily of the form:
1 0
0
(U) =
2
0
0
;
(U).(U)
=
(21+ +)
éq 3.4.2-1
0 0 - -
1
By the law of Norton-Hoff, one obtains the deviatoric constraints:
1 0
0
m2
D =
m
2
2
-
.
2
0
0
; D
3 1
with:
With m
=
+ +
=
2 1+ +
VM
()
() ()
0 0 - 1
-
éq 3.4.2-2
2
+
0
0
One deduces from
zz = 0: tr = -
3 (1+). From where constraints: =. 0
+
1
2
0.
0 0 0
3 - 2
The balance of the solid imposes that xx.(1 -) = yy
. One deduces the parameter from it =
.
1 -
3
The standardization of the loading leads to:
1
F =
éq 3.4.2-3
H (B - has) (+ (1 -))
Terms of the continuation!
m of limiting load in this case of loading is thus identical to:
2
3
2 1+ + 2
y
()
!
. F =
.
y
m
=
éq 3.4.2-4
3
+ (1 -)
3 2 -
3 +1
Handbook of Référence
R7.07 booklet: Analyze limit
HI-74/97/023/A
Code_Aster ®
Version
4.0
Titrate:
Calculation of load limits by the method of Norton-Hoff-Friaâ
Date:
12/11/97
Author (S):
F. VOLDOIRE, E. SCREWS
Key:
R7.07.01-A
Page:
16/16
4 Bibliography
[1]
ANGLES J., VOLDOIRE F., Modélisation and calculation of the load limit of a fissured component,
CR-MMN 1522-07, Sept. 96.
[2]
FRIAA A., Loi de Norton-Hoff generalized in plasticity and viscoplasticity, Thèse of doctorate,
1979.
[3]
FRIAA A., FREMOND Mr., Les methods statics and kinematics in design the collapse and in
analyze limit, Journal de Mécanique theoretical and applied, Vol. 11, NO5, 881-905, 1982.
[4]
VOLDOIRE F., Calcul with the rupture and analyze limit of the structures, notes EDF HI-74/93/082.
[5]
VOLDOIRE F., Analyze limits fissured structures and criteria of resistance, notes
EDF/DER HI-74/95/026.
[6]
VISSE E., VOLDOIRE F., Eléments finished incompressible, note EDF/DER HI-75/95/019.
[7]
VOLDOIRE F., Mise in work of the method of regularization of Norton-Hoff-Friaâ for
the analysis limits structures, notes EDF/DER HI-74/97/026.
Handbook of Référence
R7.07 booklet: Analyze limit
HI-74/97/023/A
Outline document