Chap 3 - Methodes de Resolution

Download as pdf or txt
Download as pdf or txt
You are on page 1of 53

INBTP

COURS DE MECANIQUE LINEAIRE DE


LA RUPTURE

CHAPITRE 3:
METHODES DE RESOLUTION EN
MECANIQUE DE LA RUPTURE
Rostand MOUTOU PITTI
HDR

RMP RUPTURE - Février 2023


PLAN DU CHAPITRE

• Introduction

1. Quelques éprouvettes de fissuration

2. Méthodes expérimentales

3. Méthodes numériques par éléments finis

4. Notions avancées de rupture par fatigue

• Conclusions

RMP RUPTURE - Février 2023


Plusieurs éprouvettes de fissuration ont été utilisées pour déterminer le taux de re
QUELQUES ÉPROUVETTES DE FISSURATION
d’énergie critique et la ténacité des matériaux en mode I. On peut citer par exemple :

L’éprouvette SENT (Single Edge Notch Tensile) généralement utilisée dans la m


ténacité. A cause de sa taille, elle ne présente pas une plage de propagation observable

Figure 14. Eprouvette SENT [9]


Eprouvette SENT

L’éprouvette Cantilever à inertie variable, contrairement à l’éprouvette SENT, elle


une plage de propagation stable et sa géométrie est adaptée aux essais de fissuration
mixte. 3
RMP RUPTURE - Février 2023
Eprouvette SENT 4
RMP RUPTURE - Février 2023
mixte.

Figure 15. Eprouvette Cantilever à inertie variable [10]


Eprouvette CANTILEVER

L’éprouvette DCB (Double Cantilever Beam), généralement utilisé


du taux de restitution d’énergie, ce type d’éprouvette présente 5un in
RMP RUPTURE - Février 2023
CHAPITRE III. STRATEGIE EXPERIMENTALE
G (N/mm)
2.5

2.0 Initiation de fissure


(a) : éprouvette DCB

R
L 1.5

1.0
Stabilité Rupture par instabilité
0.5

0.0
R

30 40 50 60 70 80 90 100
Longueur de fissure (mm)
T

(b) : Taux de restitution d’énergie

Eprouvette DOUBLE
Figure III.2.1. Stabilité de l’éprouvette
CANTILEVER BEAM DCB [DUB 97] [DUB 02]

Cependant, cette éprouvette présente une fissuration instable6 indui


RMP RUPTURE - Février 2023
Eprouvette DCB 7
RMP RUPTURE - Février 2023
In the past several specimen types CTSR for(Compact-Tension-Shear-Rotation)-specimen
2D-mixed-mode and loading device
problems
CHAPITRE III. STRATEGIE have been proposed [2, In
EXPERIMENTALE
theAmong
15]. last few years
others,an the
advanced AFM-specimen,
CTS-specimen together with its loading device [16, 17] has Tension Shear Rotation)
so-called CTSR (Compact
proven its applicability. Therefore -specimen (Figure 6) for cycling loading was
only the CTS-specimen
developed [22]. The appropriate loading device
will be discussed in the following.
is shown in Figure 7. This new specimen and
The loading device (Figure 4) allows applying pure
corresponding loading device enables any
Mode I-, pure Mode II- as well as almost every of
combination 2D-mixed-
mixed-mode loading including
mode-loading combination to the CTS-specimen
pure Mode I-, pure by using
Mode II- and pure Mode III-
just a uniaxial tension testing machine. ForBecause
loading. the purpose
of the of
holes, which are installed
varying the mixed-mode-loading only the loading
circularly angleatα a bolt circle around the
in 15°-steps
has to be changed. specimen, it is possible to generate any ratio of
Depending on the mixed-mode-portion
Mode I tothe crack
Mode growsIII (Figure 7).
II/Mode
By rotating the
into a new direction, i. e. the crack kinks under a specific so-called turret inside the Figure 6. CTSR-specimen
angle φ0 (see Eq. 7). Some fracturedloading device thewith
CTS-specimens ratioanof Mode II- or Mode III-
Figure 4. Loading device and
load is set (Figure 8).
increasing KII/KI-portion are shown in [18]. the adjustment of
the loading angle α

EXPERIMENTAL INVESTIGATIONS ON 3D-MIXED-MODE-CRACKS

For experimental investigations of three-dimensional crack growth and fatigue for


spatial-mixed-mode loading some types of specimens are available [2, 3, 19]. But in
fact only view of these specimens covers the full range of all basic fracture modes or all
combinations thereof.
Due to their importance for the researches of mixed-mode-loaded crack problems, in
this paper the AFM-specimen and the CTSR-specimen with their corresponding loading
devices (cf. Figs. 5 and 7) are described in more detail.

AFM (All-Fracture-Mode)-specimen and loading device


The AFM-specimen, developed by Richard, enables the
Figure III.2.2. Eprouvette
Eprouvette CTS
investigationCTS [VAL 89] [CAU
bois
of spatial-mixed-mode 87] in arbitrary
problems
combination by using a simple uniaxial testing machine [20].
Some experiments under static load Figure[21] had device:
proven the
7. Loading
applicability of this loading device. Due to its high weight and
AFM
Adjustment
specimen
of the ratio of Mode I to Mode II/Mode III

III.2.2. Eprouvettehigh
à deformation only low test frequency for fatigue is
propagation stabilisée
possible. Consequently this
enverymode mixte 8
RMP fact leads -to Février
RUPTURE high
2023 duration of
DCB CTS
Nous ne pouvons pas afficher
l’image.

2MCG
Eprouvette MMCG

9
RMP RUPTURE - Février 2023
Eprouvette 2MCG
(Mixed Mode Crack Growth) Talon supérieur

Bras en PVC

Talon inférieur

R Eprouvette bois

L
Trou de sollicitation Congé de raccordement

FII + I (β=45°) Trous de fixation


des bras
FI (β=0°)
FII (β=90°) mm
Entaille initiale
20

FI (β=0°)

FII (β=90°)

FII + I (β=45°)
10
Eprouvette
RMP RUPTURE MMCG:
- Février 2023 conception
A proposed mixed-mode fracture specimen for wood under creep loadings 201

Fig. 8 Geometric definition


ure specimen for wood under creep loadings 201

(c) FII + I ( β = 60° )

FI ( β = 0° )
FII ( β = 90° )
initial crack length

tightness t

FI ( β = 0° )
10mm

Fig. 9 Reducing of the crack growth path tightness


FII ( β = 90° )
PVC Arcan
(a)
FII + I ( β = 60° )
G 1 and G 2 are energy release rate part for open and II + IF ( β = 60° )
shear modes, s
respectively. G 1 and G 2 designate sepa-s Fig. 10 MMCG specimen and its fixation
(b) I F ( β = 0° )
ration work rates for each fractureF mode. ( β = 60° )
κ is a phe- F ( β = 90°)
II
II + I
nomenological coupling parameter. If we considerβ =G s initial crack length
3.3 Normalize energy release rate
F ( 0° ) 1
and G s2 as intrinsic material parameters, expressions
I
F ( β =tightness
II 90° ) t

(1)–(4) can be rewritten as follow initial crack length


I F ( β = 0° ) The optimized design of the geometry is based on the
10mm
tightness t
respect of crack growth stability, Eqs. (3) or (18), in
Fig. 9 Reducing of the crack growth path tightness
Phase of crack growth F ( β =initiation:
0° ) f <1 (16) an exploitable zone onF the II ( β = 90crack
°) growth path. In this
Eprouvette MMCG: conception
I PVC Arcan
mm
case, evolutions of the energy release rate or threshold
F ( β = 60° )
II + I
Phase
growth path tightness
G 1 andof Gcrack propagation:
2 are energy release rate fpart ≥ for
1 open and (17)F ( β = 90°)function can be interpreted as relative. In this case, let
G s1 and G s2 designate sepa-
II
shear modes, respectively. Fig. 10 MMCG specimen and its fixation
us consider the normalize energy release rate G̃ 11 such as
PVC Arcan

ration work rates for each fracture mode. κ is a phe- F ( β = 60° )


∂RMPf RUPTURE
II + I
lease rateCrack
part forgrowth
open andcoupling
nomenological stability: f = 1Ifand
parameter. we consider <G0s1 3.3- Normalize
(18) G 1rate
Février 2023energy release G2
12
Eprouvette MMCG: conception
Soliman, Moutou Pitti et al. CRACK PATH 2015 RMP RUPTURE - Février 2023
Eprouvette MMCG: zone de stabilité
1
(
G J / m2 ) 2
(
G J / m2 )
4,0E-01 8,0E-03

b = 90 °
3,5E-01 b = 0° 7,0E-03

3,0E-01 6,0E-03

2,5E-01 5,0E-03

2,0E-01 4,0E-03 Zone de


Zone de stabilité
1,5E-01 stabilité 3,0E-03

1,0E-01 2,0E-03

5,0E-02 1,0E-03
a(m) a(m)
0,0E+00 0,0E+00
0,04 0,05 0,06 0,07 0,08 0,09 0,10 0,04 0,05 0,06 0,07 0,08 0,09 0,10 0,11 0,12 0,13 0,14
Zones de stabilité
1
(
G J / m2 ) 2
(
G J / m2 )
5,0E-02 7,0E-03
b = 15 ° b = 75 °
4,5E-02
6,0E-03
4,0E-02
b = 30 ° b = 60 °
3,5E-02 5,0E-03

3,0E-02 Zone de
4,0E-03
Zone de b = 45° stabilité b = 45°
2,5E-02
stabilité 3,0E-03
2,0E-02
1,5E-02
b = 60 ° 2,0E-03 b = 30 °
1,0E-02
1,0E-03
5,0E-03 b = 75 ° b = 15 °
0,0E+00 0,0E+00
13
0,04 0,05 0,06 0,07 0,08 0,09 0,10 0,04 0,05 0,06 0,07 0,08 0,09 0,10 0,11 0,12 0,13 0,14
a(m) RMP RUPTURE - Février 2023 a(m)
MEHODES EXPERIMENTALES AVEC CONTACT

Méthode énergétique

W (a) + W (a + Δa)
G=−
bΔa
14
RMP RUPTURE - Février 2023
Méthode de la complaisance

Dd
C=×
DP

2
P æ ¶C ö
GC = R ×ç ÷
2 × b è ¶a ø d
15
RMP RUPTURE - Février 2023
dessous résume les principales méthodes optiques développées.
METHODES EXPERIMENTALES PAR MESURE DE
CHAMPS SANS CONTACT
Mesure de Mesure de Mesure de Mesure de Mesure de
déplacements déformations déplacements déplacements et contraintes
surfaciques surfaciques hors plan : déformations
tridimensionnelles
Suivi de Suivi de Moiré Corrélation Photoélasticimétrie
marqueurs marqueurs d’ombre d’mages 3D bidimensionnelle
Corrélation Corrélation Moiré de
d’image d’image projection Photoélasticimétrie
tridimensionnelle
Grille Grille Interférométrie

Granularité
laser

Méthode
Méthodes de mesure
optiques desdechamps
de mesure champs
sans contact
I.3.1. Méthode de suivi de marqueurs

Cette méthode consiste à suivre n marqueurs, situés en surface de la pièce étudiée,


16 au
RMP RUPTURE - Février 2023
Suivi dedes
donne la déformation dans la direction marqueurs
taches. L’analyse de n taches permet d’obte
champ de déformations.

Figure 18. Méthode de suivi de marqueurs [11]

I.3.2. Méthode de corrélation d’images 17


RMP RUPTURE - Février 2023
Suivi de marqueurs: résultats expérimentaux

Capteur LVDT
F (N)
Presse 900
FR
800

700 a = 40 mm
600
FC
500
Mode I Caméra
400
CCD
300

Eprouvette 200
2MCG 100
d (mm)
0
0,0 0,2 0,4 0,6 0,8 1,0 1,2 1,4 1,6 1,8 2,0 2,2 2,4 2,6
Dispositif expérimental en mode I (β=0°, hêtre)
Courbe force-déplacement en mode I (β=0°)

F (N)
FR
3000
a = 40 mm
2500 FC

Mode 2000

Mixte 1500

1000

500

d (mm)
0
0,0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0

Courbe force-déplacement en mode mixte (β=45°, hêtre) Courbe force-déplacement en mode mixte (β=45°, hêtre) 18
RMP RUPTURE - Février 2023
Suivi de marqueurs: résultats expérimentaux
F (N)
1600

1400 FR
1200

1000

800 F (N)
600 FR
FC 1000
400
b = 15 °
200 800
d (mm)
0
0 0,5 1 1,5 2 2,5 3 600
F (N)
1800 400
FR
1600

1400 200
FC b = 60 °

1200
0
1000 0 0,2 0,4 0,6 0,8 1 1,2 1,4
800 FC
600

400

200
d (mm) b = 30 °
0
0 0,5 1 1,5 2 2,5
19
RMP RUPTURE - Février 2023
Corrélation d’images

F
DIC:

20
RMP RUPTURE - Février 2023
Méthode de la grille

F
GM:

21
RMP RUPTURE - Février 2023
q MMCG et méthode de la grille

Part du mode 1

MMCG avec grille

Grille et camera Dispositif

- Transfert de la grille selon le protocole [1]


- p : 400 mm
- 2 tailles de ZOI (pour les essais structures)
- u et ε obtenus selon une méthode [2]
- Sensicam camera 1040x1376 pixels, 12bit,
cooled
[1] Piro J.-L., Grédiac M., Exp. Tech, 28(4):23-26, 2004
Lachhab et al. COST FP0904, Lisbon 2013 RMP RUPTURE[2]- C.Février
Badulescu
2023 et al. Meas. Sci. Technol, 2009
Maillage MMCG Filtre de Sobel
Evolution de a vs. temps

1. Méthode de la complaisance

2. Méthode numérique Mθ

3. Méthode analytique ou énergétique

23
Lachhab et al. COST FP0904, Lisbon 2013 RMP RUPTURE - Février 2023 Mode 1
Fissuration Aluminium 2074 - Castellium

Schéma
Images réelles

— Aluminium 2074, Entaille 2,5mm


— Limite élastique 372 MPa
— Module d’Young 72 Gpa
— Coefficient de Poisson 0,305
24
RMP RUPTURE - Février 2023
Machine d’essai de traction 200 KN

Deux éclairage LED

Eprouvette avec grille

Camera PCO 2000


- Dynamique 14 bits
- Temps d’exposition 5 millisecondes
- Fréquence d’acquisition : 1 image/0,5 s

25
RMP RUPTURE - Février 2023
Cartes de déformation

Données de l’essai aluminium avec entaille 6mm


Figure 15: Courbe force-temps 32
ERISATION DE LA FISSURATION PAR MESURES DE CHAMPS CINEMATIQUES
32 / 41: ALUMINIUM
Synthèse des paramètres de l’essai
I Grille de 10trs/mm (ancienne)
I Durée de l’essai : 22minutes
I 5 forces critiques

U Force(N ) a(mm) b C G
0 0 0 4 0 0
0.0812 67570 0 4 1.2.10-6 0
0.13 1930 0,681 4 6.9.10-5 47.22
0.25 1260 0,056 4 1.9.10-3 697.15
0.38 140 0,561 4 2.74.10-3 118.74
Table 5: synthèse global de l’essai 26
RMP RUPTURE - Février 2023
La méthode consiste à enregistrer sur un même film photographique de haute densité les
Granularité
taches de granularité avant et après laser
sollicitation de la pièce. L’analyse de ce cliché par
diffraction d’un faisceau laser donne une série de franges parallèles dont le pas et l’inclinaison
permettent d’accéder aux composantes dans le plan du déplacement.

(a) (b)

Figure 20. Méthode de granularité laser ; dispositif d’enregistrement (a) et dispositif de


dépouillement (b) [11]
27
RMP RUPTURE - Février 2023
I.3.4. Méthode de la grille
METHODES NUMERIQUES PAR MEF
q Formalismes analytiques
Couronne d’intégration C0 x2
x2
Contour d’intégration
Domaine intérieur Ci

x1 x1
dl
Domaine extérieur Ce dS
Domaine fissuré Ω

C0 : continu de (1,0) à (0,0)

(a) (b) (c)

28
Moutou Pitti et al. IJSS 2009 RMP RUPTURE - Février 2023
q Formalismes analytiques
C0 x2
x2
Chemin d’intégration
Ci

x1 x1
dl
Ce dS
Domaine fissuré Ω

Intégrales de Rice et Gθ

Intégrales M et Mθ (modes mixtes)

29
Moutou Pitti et al. IJSS 2009 RMP RUPTURE - Février 2023
Intégrale T en fissure statique

x2

Δa
fissure
x1

Intégrale A en propagation

A1: Terme classique A2 : Effet de la variation de ΔT

A3 : effets de pression sur les lèvres de la fissure A4 : effet de la propagation de la fissure

30
Moutou Pitti et al. IJF 2010 RMP RUPTURE - Février 2023
Interprétation physique
Champs virtuels:
Taux de restitution K Iu × K Iv u
K II v
× K II - élastique: champ d’Irwin
d’énergie virtuel A= +
C1 C2 - orthotrope: forme de Shi

Séparation des modes mixtes


Mode I Mode II

K Iu × K Iv K II
u v
× K II K Iu × K Iv K II
u v
× K II
A= + A= +
C1 C2 C1 C2

1 0 0 1

K Iu = A × C1 u
K II = A × C2

A K Iu × K Iu A
u
K II u
× K II
G1 = G2 =
C1 C2
31
Moutou Pitti et al. Springer 2015 RMP RUPTURE - Février 2023
q Applications : milieux isotropes, Mode I pure (Int M)

2b

c2c1c3
2h

c4 c
2a c6 5
c7

Matériau Chargement
E = 20000 daN/mm2 ν = 0,3 σ = 1 daN/mm2

Géométrie
2b = 400 mm 2h = 1200 mm 2a = 200 mm

Calcul en déformations planes

32
RMP RUPTURE - Février 2023
q Applications : milieux isotropes, Mode mixte (Int M)

2b

2a
2h

30°
c2c1c3
c4 c
c6 5
c7

Calcul en contraintes planes

33
RMP RUPTURE - Février 2023
q Applications : milieux isotropes, Axisymetrie (Int M)

368 Fig. 4 Crack growth process position in axisymmetric structure [19] R. Moutou Pitti et

(a) z (b) z (c) β = 90 °

A7 β = 45 °
β
Γ1 circular crack a

€ n1 A7
A1 Γ2
A2 B1 Stee
d = 100 mm
x2 = z β = 0°
Wood specimen
B2 n r A1
dθ B1 2

x1 = r V
€ Crack tip
β


€ B7

Fig. 2 Integration domain Fig. 5 The modified Compact Tension shear specimen

Axisymétrique : coordonnées (a) ;


domaine
Applying the Gauss–Ostrogradski d’intégration
theorem to the second (b)
terms; of Eq. (10) leads to
4 Numerical applications
! ! • Eprouvette CTS (c)
∂F 4.1 Numerical specimen and elastic components
δL = κk δadtdV
∂ xk In the axisymmetric configuration, the circular crack evolution and the integratio
t V
! ! " # [19]. The crack !tip is
"oriented in the radial direction
# and the continuum map θ s
∂F • ∂F • In order to compute •
∂ Fthe energy release •
∂ F rate in the axisymmetric material,
+ u i,k + vi,k κk n jShear
δadtdS − Fig. 5, initially
specimen, + by vRichard
u i,kproposed δadV
i,k κk[20] and used by Moutou P
∂u i, j ∂vi, j ∂u ∂v
i,t This specimen
material has been considered. i,t is used in computational
34finite elem
t ∂V mixed-mode V in the case of axisymmetric configurations. The initial square
ratios
Moutou Pitti et al. Acta
Meccanica 2011 RMP RUPTURE - Février 2023 
q Applications : milieux isotropes, Axisymetrie (Int M)

Ri t

c2c1c3
a a c4 c
c6 5
c7

P
Matériau Chargement
E = 2 1011 Pa ν = 0,3 P = 1 106 N

Géométrie
Ri = 1 m t = 0,1 m a = 0,05 m

35
RMP RUPTURE - Février 2023
q Applications : milieux isotropes, pression sur les lèvres

Propriétés matériaux
E = 20000 daN/mm2 ν = 0,3

Paramètres géométriques
2b = 400 mm 2h = 1200 mm
2a = 200 mm

Force appliquée
σ = 1 daN/mm2
Géométrie et chargement Maillage MEF Déformée

Couronnes d’intégration

c2c1c3
c4 c5
c6
c7

Comparaison Taux de restitution FIC

36
RMP RUPTURE - Février 2023 Moutou Pitti et al. CFRAC 2015, Springer 2015
q Application: milieux isotropes, intégrale A en mode I

(a) – Plaque rectangulaire soumise à un champ thermique (a) – Géométrie, (b) – Maillage EF, (c)
– Maillage avec champ thermique, (d) – Etat de déformation

Temperature field distribution


!! −! !
!! = 2!! !! ∈ , !! = 100!°!
! 2 2 37
Riahi, Moutou Pitti et al. MTDM 2014 RMP RUPTURE - Février 2023
(a) – Indépendance du domaine d’intégration en GI ; (b) – KI numérique et analytique

38
Riahi, Moutou Pitti et al. MTDM 2014 RMP RUPTURE - Février 2023
q Application: milieux isotropes, intégrale A et pression sur les lèvres

(a) Distribution du champ de température;


(b) indépendance du domaine d’intégration KI

39
Riahi, Moutou Pitti et al. MTDM 2014 RMP RUPTURE - Février 2023
(R, T, L), pour le bois.
q Formalisme généralisé : milieu orthotrope
L
R

Coupe transversale

T Coupe radiale

Coupe tangentielle

Figure I.1.3. Orthotropie locale : deux plans de symétrie matérielle.


Contrainte plane Variation de température
x2
Thèse de Doctorat de Rostand MOUTOU PITTI
8
E2 Soutenue le 15/10/2008

E1 x1
Hyp1 : γ = f(E1,v12,α1)

A1: Classical term A2 : temperature variation effect 40


Riahi, Moutou Pitti et al. MTDM 2014 RMP RUPTURE - Février 2023
q Application : milieu orthotrope, modes mixtes
Arcan
Maillage rayonnant
et champ θ

Wood

Eprouvette CTS Maillage

Indépendance du domaine pour ΔT=0 41


Riahi, Moutou Pitti et al. MTDM 2014 RMP RUPTURE - Février 2023
Indépendance du domaine pour ΔT= +10°C

Indépendance du domaine pour ΔT= -10°C 42


Riahi, Moutou Pitti et al. MTDM 2014 RMP RUPTURE - Février 2023
(a)l’humidité
q Application : Effet de (b) (c) (

(a) (b)
Mesh of wood MMCG specimen
(b)
specimen
(c) (d)

(a) (b) (c) (d)


In mode II

43
Riahi, Moutou Pitti et al. MTDM 2014 RMP RUPTURE - Février 2023
(a) (b)
q Application : Effet de l’humidité
Numerical mesh of wood specimen under
Thermo-visco-hydro-Mechanical (TVHM)
fields at 12% of internal moisture content

(a) (b) (c)


ct ion
e
x dir
Ex

(a) (b) (c)

Ey
yd
i re
ctio
n

Real wood specimen

44
Riahi, Moutou Pitti et al. MTDM 2014 RMP RUPTURE - Février 2023
Evolution of GI versus crack a for β=45°

45
RMP RUPTURE - Février 2023
Evolution of GII versus crack a for β=45°

46
RMP RUPTURE - Février 2023
q Aciers à haute limite d’élasticité

Evolution of GII versus crack a for β=45°


47
RMP RUPTURE - Février 2023
110

100

90

80

70
load (kN)

60

50

40

30

20

10

0
0 1 2 3 4 5 6
displacement (mm)

48
RMP RUPTURE - Février 2023
q Deformations

x
49
RMP RUPTURE - Février 2023
εxx
EXPERIMENT FE SIMULATION
-3 -3
x 10 x 10
4 4
40
100 3 3
35
200
2 30 2
300
1 25 1
400
500 0 20 0

600 -1 15 -1
700
-2 10 -2
800
-3 5 -3
900
-4 0 -4
200 400 600 800 1000 1200 0 10 20 30 40

x
50
RMP RUPTURE - Février 2023
εyy
EXPERIMENT FE SIMULATION
-3 -3
x 10 x 10
15 15
40
100
35
200
30
300 10
10
400 25

500 20
600 15
5 5
700
10
800
5
900
0 0 0
200 400 600 800 1000 1200 0 10 20 30 40

51
RMP RUPTURE - Février 2023
εxy
EXPERIMENT FE SIMULATION
-3 -3
x 10 x 10
5 5
40
100
35
200
300 30

400 25

y(mm)
500 0 0
20
600
15
700
10
800
900 5
-5 0
200 400 600 800 1000 1200 -5
0 10 20 30 40
x(mm)

52
RMP RUPTURE - Février 2023
UFR – MR2

COURS DE MECANIQUE LINEAIRE DE


LA RUPTURE

CHAPITRE 3:
METHODES DE RESOLUTION EN
MECANIQUE DE LA RUPTURE
Rostand MOUTOU PITTI
MCF - HDR

RMP RUPTURE - Février 2023

You might also like