1s s+5,l'?

t*

OfC'.IOf fag,-fi.l-I:UlIi\/F-t('i:;'i\|Iirr - /i5TZ

; t

2po3clsJi

UNIVERSITE DE ME
THESE Pourobtenirle gradede DOCTEURDE L'UMVERSITE DB METZ Discipline Mcanique :
Prsente : Par

HamidrezaADIB RAMEZANI
Titre :

de Evaluation la Durede Vie en Fatigue pour des par la Mthode Volumtrique par JointsSouds Point
d'examen publiquementle23 la Septembre devant commission 2@3 Soutenue H. P. LIEIJRADE M. NAIT-ABDELAZIZ G. PLUVINAGE Z. AZAKT J. GILGERT J. Y. BARTMLEMY
CE:TIM(Senlis) Professeur l'Universit de Lille Professeur l'Universit de Metz Professeur l' Universitd' Amiens M. C. l'Ecole Nationale d'Ingnieur de Metz Ingnieur de Recherche,Instilut de Soudure
DE UNIVERSITAIRE METZ BIBLIOT}IEOJE

Rapporteur Rapporteur Examinateur Examinateur Examinateur Invtt

llll illlIilil illl il I lill ffiffiil ililillil ill

022 486098 I

l'| fl rversrtc!

rl ,HL
- JJtetzS

frwtF-R{fue; bfr fu&t kns

AVERTT55EMENT

Cette thse est le fnuit d'un longtrovail opprouv le por jury de soutenonce disponible l'ensemblede lo et commun universitoire aut lorgie. Effe est soumise lo propritintellectuelle l'outeur ou de mme titre que so version popier. Ceci impligue une igotionde citotion de rfrencement obf donslo rdoction , de tousvosdocuments. plagiots,rep?oductions D'outre port, toutes contrefoons, illicitesentronent poursuite pnale. ih, l'outonisoton une Enf jusgu' nouvel de diffusiono t,occorde ordre.
i doniel.michel@scd.univ-metz.fr

CPI orticfes L 122. 4 CPI orticles L 335.2- L 335.10 htt ot/ / www.cf ies.co V2/ leo leo droi.oho coo ml / http://www.culture.gouv.f r/culture/infos-protigues/droits/protection.htm

UNIVERSITEDE}/4.ETZ
THESE Pourobtenirle gradede DOCTEURDE L'UNIVERSITE DE N4ETZ Discipline: Mcanique
Par Prsente :

Hamid rezaADIB RAMEZANI

Titre :

de Evaluation la Durede Vie en Fatigue pour des par la MthodeVolumtrique JointsSouds Point Par
d'examen la 2003 le publiquement 23 Septembre devant commission Soutenue IT P. LIEI.]RADE M. NAIT.ABDE,LAZIZ G. PLUVINAGE Z. AZAHT J. GILGERT J. Y. BARTIIELEMY
CETIM (Senlis) Professeur l'Universit de Lille Professeur l'Universit de Metz Professeur l' Universitd' Amiens M. C. l'Ecole Nationale d'Ingnieur de Metz Ingnieur de Reclrcrche,Institut de Soudure Rapporteur Rapporteur Examinateur Examinateur Emrninateur Invit

|(esume

T\^,

Les joints soudspar point, utiliss dans I'industie automobileet notammenten carrosserie sont gnralement soumis des chargements cycliques. La, caractristique du chargement appliqusur les joints soudspar point conditionnele choix de I'analyse. Dans les deux cas de chargementcyclique ts lev (fatigue oligocyclique) et moyen (endurance dansle mtal de baseet dans limit),I'amorage de la fissure se situe respectivement la zone affecte thermiquement. Les joints soudspar point en cisaillement avec un, tois et cinq points de soudure sont choisis pour obtenir la dure de vie en fatigue basesur I'approchevolumtrique dans les cas de fatigue oligocyclique et fatigue endurance limite. La distanceeffective, la contrainteeffective et le facteur de rduction en fatigue sont calculs par I'approche volumtrique en utilisant la mthodedesElmentsFinis avecdesmodlestridimensionnels. La technique de sous-modelage appliquepour obtenir le champ des conaintes et le est gradient des contraintesau voisinagedu point soud.I-es rsultasobtenus partir de I'approche volumtrique s'avrent en ts bon accord avec ceux obtenus exprimentalementsur des prouvettessoudesdans le cas de fatigue endurancelimite et restent satisfaisantdans le cas de fatigue oligocyclique.

Remerciements

Ce travail a t effectu au Laboratoire de Fiabilit Mcanique, Universit de Metz UFRMIM sous la direction de Monsieur G. PLUVINAGE, directeur du Laboratoire que je tiens remercier pour son soutien et pour son implication dans la direction scientifique de mes recherches. Mes remerciements s'adressent : MessieursZ. AZAKI de I'Universit d'Amiens prsidentdu jury et J. GILGERT de I'Ecole Nationale d'Ingnieur de Metz membre du jury. Messieurs H-P LIEURADE du Cente Technique des Ingnieurs Mcaniques et M. NAITABDELAZIZ de I'Universit de Lille, rapporteurs de cette thse, pollr la pertinence de leurs conseils. Monsieur Jean-YvesBARTHEIMY d'Institut de Soudure,membreinvit du jury. Mademoiselle Jena JEONG, ma collgue pour ses rerurques constructives et sa grande disponibilit. Je ne saurai oublier Monsieur P. JODIN et tous les membres du Laboratoire de Fiabilit Mcaniquepour leur disponibilit et leurs remarques constuctives.

A ma mre

TABLE DESMATIERES
Page Introduction gnrale: ......... I

1.0

Evaluation de la dure de vie en fatigue desjoints soudspar point.

Mthode de prdiction de la dure de vie en fatigue en prsence de concentrationde contraintes.

2.0

gnrale Introduction
Distribution des contraintesen fond d'entaille.

2.r

2.2.1 Distribution des contraintesdansle cas lastique. 2.2.2 Distribution des contraintesen fond d'entaille dans le cas lasto-plastique. 2.2.3 Mthode des Elment Finis dans la cas lasto-plastique. 2.3 2.4 2.5 2.5.r
Concept facteurde concentration contraintes de des Facteur rduction fatigue...... de en Ismthodes d'estimation la durede vie en fatigue. de La mthode I'Eurocode3 pourlesjoints souds. de ..

t4 22 28
45 53 62 64 78 86

2.5.2 Mthodede dformation locale. 2.6 3 3.0 3.1 3.2 3.3

Rsum.

Mthode volumtrique pour la prvision de la dure de vie de structures 87 e n t a i l l e s................ .... .. .. .. .. .. .. .. .. .. .. .. . .. .. .. .. .. .. .. .. gnrale Introduction L'approche volumtrique..... Procdure de prvision de la dure de vie en fatigue par la mthode volumtrique..... Rsum. .. 88 94 99

105

Bibliographie point. Introductiongnrale.

fatigue

des

joints

souds

par

106

4.0 4.1 4.2 4.3 4.4 4.5

5

r07
110 lL2 123

Loi de propagation fissureen fatigue. de par Structure joints souds point. des pour lesjoints souds point. Facteur par d'intensitde contraintes Tenue I'endurance fatiguedespointssouds. en Rsum. par Etude exprimentale joints souds point....... des Matriau.

t32
135 137 138

5.0 5.1

Composition chimiqueet comportement mcanique matriau joints souds 140 du des par point. par Gomtrie joints souds point. des par Procdure opration soudage point... des de Prsentation essais des fatigues.. Rsultats exprimentaux cisaillement...... Rsum. par point. Modlisation lastiqueet elasto-plastique joints souds des Introduction... ... par Modlisationdesjoints souds point Mono Point, Triple Pointset Multiple Points. des joints souds par point en 142 t45 152 157

5.2 5.3 5.4 5.5

5.6 6 6.0 6.1

160 162

r63 r66
168 175 182

par 6 . 1 . 1 Modlisation joints souds point Mono Point. des par 6 . 1 . 2 Modlisation joints souds point Triple Points. des

6 . r . 3 Modlisationdesjoints souds point Multiple Points. par 6.2

Prvision de la vie en fatigue des joints soudspar point par I'approche 188 volumtrique....

6.3 7
Rfrence

Discussion...... Conclusiongnrale..

193 L97 200

bi-axiale par Figure 1-1.Type de sollicitationdesjoinS souds point a) cisaillementb) nsion c) pelage. au de Figure2-l Illustrationde la concentration contraintes bord d'un trou dansuneplaque. autourd'un trou Figure 2-2 Distribution de la contraintecirconfrentielleOr(r,0=/t12) vertical en compression. vertical en action b) Chargement dansune grandetle a) Chargement Figure 2-3 Distribution de contraintesautour d'un hou circulaire dansune grandetle par la y x -thode des ElmentsFinis @gure 2-2a) a) Composante de contrainteb) Composante de contrainte. Figure 24 Distribution de contraintesautour d'un tou elliptique dansune grandetle obtenue par la mthode des Elments Finis @gure 2-2b) a) Composantex de contraintes b) y Composante de contraintes. 6 7 polaires correspondantes. d'une entailleet coordonnes Figure2-5 Schma Figure 2-6 Schma d'une entaille moussedans une plaque soumise une traction au perpendiculaire pland'entaille. des Figwe 2-7 Conditionsaux limites, configuration maillageset distributiondes contraintes 2'6,y comprislI p = 100. dansla Figr.ue pow le modleslectionn de Finis, dessolutions Neuber,Glinka, des avecla mthode Elments Figure2-8 Comparaison pour valeurs rapportI I p 1l = 2a)' du Kujawski,Chenet Paris-Creager diffrentes

Page 3

8 1l

t2

11

t4

r6
18

19

r0
ll

soumise sur lasto-plastiques le ligamentd'une prouvette Figure2-9 Evolutiondescontraintes uneflexion selonWilshaw[23]. sur Figure 2-10 Evolution des contrainteslasto-plastiques le ligament d'une prouvette uneflexion selonBhattacharya soumise [21].

23
25

t2

le lasto-plastiques long du ligamentd'une prouvette 26 Figure 2-11 Evolution des contraintes d'une zoneplastifiede petitetaille selonBhattacharya soumise une flexion pour le problme

l2rl.
l3 d, applique . en d'une entaillecourteet contrainte tractionglobale Figure2-I2 Schma du des Figrue2-13Composants tenseur contraintes. dans les calculs par Figure 2-14 Diffrents comportementsqui s'appliquent gnralement b) multilinairecinmatique c) bilinaire isohope d) ElmentsFinis a) bilinaire cinmatique multilinaire isotrope e) anisotropef) Type Drucker. a) du Figure 2-15 Prsentation comportement cinmatiqueb) anisorope c) Type Drucker dans et lescasbi-dimensionnel ti-dimensionnel. tel isotropeet cinmatique qu'il existedans des matriauxnonFigure 2-16 Comportement lasto-plastiques. linaires et Figure 2-17 Algorithme itratif pour les casnon-linaires lasto-plastiques mthodede non-linairelasto-plastique, Figure 2-18 Algorithme itratif pour un comportement Divergence. b) substitution a) Convergence

26 28 33

t4 15

l6

34

r7
18 19

35 35 38

USTE DESFIGURES

20

Figure 2-19 Algorithme itratif dans la non-linaire lasto-plastique, la mthode de NewtonRaphson modifi.

40
4I

2I

Figure 2-20 Algorithme itratif dans la non-linaire lasto-plastique, mthodede Newtonla Raphson. Figure 2-2I Algoithme itratif, mthode incrmentaleou pas pas avec une itration de Newton-Raphson. Figure 2-22 Algo,rithmegnral itratif dans le cas d'un comportementnonlinaire lastoplastique. Figure 2-23 Distibution des contraintes dansle cas d'un fou dansune plaque,y comprisle chargement gnral. Figure2-24 Distibution des containtesdansle casd'un trou dansune grandeplaquepour un chargementoi = oo , otr = 0 a) Maillageset conditionsaux limites b) Composantes la de
-h

22

42 43
45

23
24

25

48

contrainte dansla directionx c) schma contraintes voisinage trou. des au du 26 Figure 2'25 Comparusonentre la solution analytiqueet le calcul par ElmentsFinis pour le problmed'un petit ou dansune grandeplaqueavecun chargement horizontalde tractio. Figure 2'26 Disnbrtion des contraintes dans le cas d'un trou dans une plaque pour les composantesradiales et circonfrentielles dans une plaque y compris le chargement a) 46

27

49

o ! = o o , o I = o o o ! = o 0 , 6 - - o o .

28

Figure2'27 Distributiondes contraintes dansle casd'une barreen taction D=0.O1 m, d=0.02 m, r=0.01 m, d0 =L MPa a) Schma b) Composant des contraintes c) Gradientde Y contrainte dansla rgioncritique. Figure2-28Schma d'une prouvette entaill avecdfinitiondesparamtres et A. h Figure 2'29 Charry des contrainteset gradientdescontraintes fond d'entaille a) Champdes au contraintesen traction b) Gradient des contraintesen traction c) Champ des conaifts en flexion d) GradientdesconFaintesen flexion. Figure2-30 Schma diagramme Whleret relationde Basquin. du de Figure2-31 Schma diagramme Whlerpour desprouvettes du de lisseset entailles. Figure 2-32Histotre descontraintes a) Amplitude constanteb) Amplitude variable. Figure 2-33 a) Distribution du domainedescontraintesb) Varit descontraintes ordre. par Figure 2-34Diagramne de rsistance fatigue de la mthodede I'Ewocode 3 a) Diagramme en avecuneligne b) Diagramme avecdeuxlignes. Figure 2-35 Schmacumulatif des amplitudesdes conhaintesappliquesa) Speche des contraintes appliqus b) Schma amplitudes contraintes des des dansle diagramme Wohler. e Figure 2-36 Diagramme Normal de I'Eurocode3 avecles limis des domaines contrainte de O, correspondantcinq millions de cycles. Figure2-37 Algorithmedu dimensionnement Eurocode pour un diagramme uneet deux par 3, lignesS-N.

50

29 30

57 60

31 32 33 34 35 36

62 63 65 67 6g

69

JI

7l

38

75

ll

USTE DESFIGURES

39

d'une contrainte de Figure2-38 Schma des mthodes la contraintedansle cas d'existence moyenne. des Figure 2-39 Schma la superposition courbeslastiqueet plastiquepour obtenir la de courbede dformationtotale en fatigue. locale[84]. Figure2-40 Illustrationde la mthode dformation de de Figure2-41 Illustation de I'interprtation la mthode Neuber. de Figure2-42 ltlgonthme gnrald'utilisation de la mthodede Neuber. de Figure243 Illustrationde I'interprtation la mthode Molski-Glinka. de Linaires,Neuberet Molski-Glinka. Figure2-M Comparaison mthodes des Figure3-l hsentationde la distributionschmatique fond d'entaillepour diffrentstypes en d'prouvettes. Figure3-2 Schmade I'intgration voisinage fond d'entaille a) Mthodede la contrainte au du effective b) Mthodede Peterson. du Figure3-3 Schma de la fonctionde poids au voisinage fond d'entailledansla mthode d'intensit champdescontaintes. du des en entailles dansle cas Figwe 3-4 Schma la distribution contraintes fond desdiffrentes de d'un calcullasto-plastique. de Figure 3-5 Schma diagramme du bi-logarithmique la distributionde la contrainte lastoplastiqueet du gradientdescontraintes fond d'entaille. en Figure3-6 Distributon descontraintes fondd'entailleet en fond de fissuredansla zoneIII. en Figure3-7 Algorithmede la mthode volumtrique dansle casde I'existence la courbelisse. de Figure3-8 Algorithmede la mthode volumtrique dansle cassans courbelisse. la Figure4-1 Schma ois modes rupture. des de qn Figure4-2 Schma d'unefissureet lescoordonnes fond d'une fissure. Figure4-3 Schma des dansles mtaux. typiquede la propagation fissures dansles casde fissureslongueset Figure4-4 Schma typiquede la propagation fissures des courtes [84]. par Figure 4-5 Schma la procdure soudage point. de du

76

40

79

4l 42 43 44 45 46

80 81 82 83 84 89

47

9L

48

92

49

95

50

96

5I 52 53 54 55 56 57

98 100 101 lO7 108 110 lll

58 59 60 6l 62

ll2

par Figure4-6 Schma soudage point avecla zonefondue,la zone afft& thermiquement ll3 du et I'interface. par Figure4-7 Schma modle du simpledujoint soud point assimil unepoutre. ll4

Figure4-8 Schma modlesimpledu joint soudpar point avecI'utilisationdes lments 116 du poutre,plaqueet < pin joint >. par Figure4-9 Schma modlepour un joint soud I'utilisationde la thoriedesplaques et du coques [94]. ll7

lll

USTE DESFIGURES

63

par Figure4-10 Schma modlepropos Kan t98l a) Schmagnralb) Modle simplifi du c) Modlebi-dimensionnel. Figure4-11 Schma modlepour un joint soudavecI'utilisationde la theoriedesplaques du et des coquesa) Mthode de la combinaison des lments coqueset solides [85], b) et c) Mthodeavecpoutes et coques 191,941. Figure 4-12 Schma l'quilibre des moments des forcesautourdu joint soudpar point de et

ll7

64

118

65

I2O

teel.
66 Figure 4-13 Effet de I'interfaceentre deux tles mincesdansdesjoints soudspar point en cisaillement. Figure 4-14 Localisationdes fissures dans des diffrentscas a) Fatigue oligocycliqueb) Endurance limite. Figure4-15 Schma dtaillde la fissuredansI'utilisationde la mthode Elments des Finis. Figure4-16 Schma dtaillde la fissureautourdu point soudet I'angle de la fissure partir par de la formulepropose Cooper-Smith tl10l. Figure 4-17 Contraintes radialeset cisaillements circonfrentiels autour du point soudbas selonl'tudede Radaj-Zhang l1]. I Figure 4-18 Cas gnralde I'applicationdes actionsautour du point soudselon SwallemKurath-Lawrence 12]. [1 propose Swallempar Figure4-19 Procdwe d'valuation la valeurde Do selonla mthode de Kurath-Lawrence a) Schmades rsultatsexprimentaux dansle diagrammebiJogarithmique K, - N r b) Utilisationdesrsultats obtenus dansle diagramme logarithmique - (1- R) . K, Figure 4-20 Digrammebi-logarithmique s*-N/(t-R) par fatiguepar I'approchepropose Sheppard. pour valuerla durede vie en I2O

67

l2l

68 69

124 125

70

126

7l

127

72

130

73

133

74

Figure 5-l Ini de comportement cyclique (courbed'amplitudede contrainte-amplitude de dformation). par y Figure5-2 Gomtrie joints souds point en cisaillement comprisMono Point,Triple des Pointset Multiple Poins. par Figure5-3 Aspectmicroscopique d'une prouvette soude point a) Grandissement point du soudX50 b) Grandissement point soud du X500 c) Grandissement point soudX50 d) du X500 F 191. Grandissement point soud du pour la taille du point soud, I'interfaceet Figure5-4 Schma la gomtrie point soud du de de de l'paisseur. poursoudure joints souds point. par applique des Figure5-5 Schma l'lectrode de Figure5-6 Schma paramtres et Z, dansle casdesjoints souds des Z* multi points. par Figure5-7 Schma soudage point, typemulti points. du par point a) sansvariation de l'effort de Figure 5-8 Cycle de soudage joints souds des compression avecvariationde I'effort de compression. b) par Figure5-9 Cycle de soudage desjoints souds point pour les prouvettes slectionnes a)

l4l

75

142

76

143

77

144

78 79 80 81

L45 147 148 149

82

150

lv

USTE DES FIGURES Mono Point b) Triple Pointset Multiple Poinrsu 19l (l Priode= 0.02Sec.). Figure5-10 Schma l'lectode de soudage de pour desjoints souds hmisphrique applique par point. 84 Figure 5-11 Positionnement l'prouvettesoudedans les bridageshydrauliquesde la de < >. machine Schenk Figure5-12 Schma positionnement joints souds point Mono Point, Triple Pointset par du des Multiple Points. 86 Figure 5-13 Schmadu positionnement fissuredesjoints souds par point Mono Point, de Triple Pointset Multiple Poinr. Figure5-14 Propagation fissureen fond d'un point soudavecforce maximale2100 N et de nombre cyclesappliqus800000a) 250000b) 300000. de Figure 5-15 Mode d'endommagement prouvettes par des soudes point en cisaillement a) Mono Point b) Triple Pointset Multiple Points. Figure 5-16 Prsentation I'endommagement joints souds point en cisaillement de par des a) Schma I'endommagement Eprouvettes de b) endommages. 151

152

153

153

87

154

155

156

90
9l

Figure 5-17 Cowbe de Whler pour les prouvettes par soudes points (Mono Points, Triple Pointset Multiple Points). Figure 5-18 Courbede Whler pour les prouvettes par soudes points(Mono Points,Triple Pointset Multple Points)avecI'utilisationdu facteurn . Figure6-l Utilisationdu logicieldesElments Finis par GIII. Figure6-2 Utilisationdu logicielpar la mthode programmation. de Figure6-3 Algorithmede la mthode pour la mthode Elments sousmodelage des Finis. Figwe 64 Algorithme gnral calcul desjoints souds point dansles cas linaireset par du non-linaires. Figure 6-5 Modlisationglobaled'un joint soudpar point Mono Point a) Maillage b) Agrandissement maillage c) Conditionsaux limites et chargement Modle complet e) du d) Agrandissement maillage. du globale Figure6-6 SousModlisation par d'un joint soud point MonoPoint. Figure6-7 Schma sousmodle d'uneprouvette par du soude pointMono Point. Figure 6-8 Prsentationdu dplacementdu modle gnral du joint soud par point en cisaillement Mono Pointpar calcullastique. 158

92 93 94 95

r64 165 t66 t67

96

169

97 98 99

17l r70 172

100

pour le joint soudp point Mono Point (lasto-plastique) t73 Figure 6-9 Champdes contraintes a) Dplacement champdes contraintes b) Champdes contraintesautour du point soud(la et containte maximaledansla directionx et la dformationtotale dansce point). Figure 6-10 Evolution du facteurde concentration contraintes de lastiquesft,,du facteu de concentration contraintes de lasto-plastiques et du facteurde concentration dformations &o de lasto-plastiques dansla direction longitudinaledu joint soudpar point Mono Point ft, ( f est le rapport ene la distancedans la direction x qui commenceau cene du point soudet le 174

101

USTE DESFIGURES

diamtre point soud). du

t02

Figure 6-11 Modlisationglobaled'un joint soudpar point Triple Poins a) Maillage b) Agrandissement maillage c) Conditionsaux limites et chargementd) Modle complete) du Agrandissement maillage. du

176

globaled'un joint soud point Triple Points. par 103 Figure6-12 SousModlisation par 104 Figure6-13 Schma sous-modle soude point Triple Points. du d'une prouvette

t77 t78 t79

105 Figure 6-14 Dplacement modlegnraldu joint soudpar point en cisaillement Triple du
par Poin calcullastque. 106

par pour le joint soud point Triple Points(lasto-plastique) 180 Figure6-15 Champdescontraintes a) Champs des conaintes b) Champ des contraintesautour du point soud (la contrainte maximaledansla direction x et la dformationtotale dansce point). 180

r07 Figure 6-16 Evolution du facteurde concentration contraintes lastique(, du facteurde de

de ko concentration confiainteslasto-plastiques et du facteurde concentration dformations de par lasto-plastiques Ia directionlongitudinale joint soud point Mono Point ( est dans du le rapport entre la distance dans direction x qui cornmenceau centre du point soud et le du diamtre point soud). 108 Figure 6-17 Modlisationglobaled'un joint soudpar point Multiple Points a) Maillageb) Agrandissement maillage c) Conditionsaux limiteset chargementd) Modlecomplet e) du Agrandissement maillage. du globaled'un joint soud point Multple Points. par Figure6-18 SousModlisation

r83

109

r84
185

1 1 0 Figure6-19 Dplacement modlegnral joint soud point en cisaillement par du du Multiple Pointspour le calcullastique. 1 1 1 Figure 6-20 Champ des containtes pour le joint soud par point Multiple Points (lastoplastique) a) Champsdes contraintes b) Champ des contraintesautour du point soud (la contraintemaximaledansla direction x et la dformationtotale dansce point). ll2 Figure 6-21 Evolution du facteurde concenFation contraintes de lastiques(, au facteurde concenfation de contraintes o. lasto-plastiques et du facteurde concentration dformations de par lasto-plastiques la directionlongitudinale joint soud point Mono Point fr, ( 6 est dans du le rapport enhe la distancedans la direction x qui commenceau cenfte du point soudet le du diamtre point soud). 113 Figure 6-22 Schma la directiond'ouverturede la fissureet la containte d'ouvertureen de fatigue. Figure6-23Schma la distribution contraintes de des dansl'paisseur d'un point soud critique. Figure6-24 Disnibutiondescontraintes dansl'paisseur d'un point soud critiquepow obtenir par les facteurs rduction fatiguedesjoins souds point a) Mono Pointb) Triple Pointsc) de en Multiple Points. par Figure6-25Comparaison volumtrique lesrsultats ente lesrsultats obtenus I'approche et par exprimentaux Rsultats joints souds pointMonoPoint considrs a) des cornme rfrence par b) Rsultats points souds point Triple Pointsconsidrs des commerfrence.

r86

186

r 88
189 190

ll4 I 15

116

r92

pour donnantla disfibution des conhaintes analytiques des Tableau 2-1. Rsum mthodes en et mousss fonction de la distance. desentaillesaigues et Tableau2-2 Ditrrents chargements facteur de concenationde contraintespour un petit une plaque. trou dans dans les cas de contraintes des Tableau2-3 Types du changement sectionset chargement planes axisymtriques et [28]. pour l'estimation du facteur de rduction en Tableau 24 Rsumdes relations proposes fatigueftr. pourla relaton entek,, ko et kr. par auteurs proposes diffrents 2-5 Equations Tableau 6 Tableau 3-1 Description des zones de distibution des contraintes dans I'approche volumtrique. pour I'approchevolumtrique. analytiquepropose de Tableau3-2 Prsentation la mthode et en de technique la qualitSolphor acierrephosphor haute 5Tableau 1 Caractristique limi lastique. chimiquede I'aciertudi. Tableau 5-2 Composition pour I'acierSolphor280. Tableau 5-3 Conditiond'essais de 5-4 Proprit I'acierSolphor280. Tableau cyclique. de Tableau Coefficients la loi de Ludwick ou de la loi de comportement 5-5 des points soudspour les prouvettes de Tableau5-6 Rsumde la procdure soudage slectionnes. par pour lesjoints souds point en cisaillement obtenus exprimentaux Tableau Rsultats 5-7 (MonoPoint,Triple Pointset Multiple Points). par absolueet relativepour les joints souds point en cisaillement Tableau5-8 Rsistance (MonoPoint,Triple Pointset Multiple Points).

Page 20

46

51

56 85 99

7 8

lO4 139

9 l0 1l t2 13

l4O I4O l4O I4l 150

L4

I57

15

158

vll

CHAPITRE T

Chapitre 1. Introduction gnrale: Evaluation de la dure par points de vie en fatigue desjoints souds
Le soudagepar points est une techniquelargementutilise dans les conceptions utilisantdes tles mincestellesque les carrosseries d'automobiles, cabinesde les camion,les wagonsde cheminde fer etc. Il est parfoisemploydans desstructures (ex. mtalliques hautetechnologie aviation). de par Le soudage point est largementrpandudans llndustrie automobile.Les tles minceset feuillardsutilises dansles structures assembles cettetechnique par sont de soudage associe celle de la robotisation. est galement Il employdansd'autres industries tellesquel'lectromnager. Il y a 3000pointssouds dansun vhicule normal.Dansles valuations contraintes de desjoints souds point, les chargements considrs par sont commeun cisaillement (TensileShear), tension une bi-axiale(Cross Tension) un pelage(CoalPeel). et La rsistance fatiguedesjoints souds point est contrle la concentration par en par de contrainte bord du point de soudure les criquesde fatiguesont amorces au o sur le ctintrieurdesplats.La concenffation contraintes peuttrevalue la base de sur (c.--d.les contraintes de la valeur maximaledes contraintes, locales dtermines pour un modlede coquetri-dimensionnel la structure). de Gnralement, soudures uneconcentration contrainte les ont de leve la priphrie du point soud. Cette concentration contraintes leveen raison de la forme de est particulire pointset deschargements. des Les contraintes calcules dpendent chargement, diamtreet de la naturedu du du point soud, la naturedu soudage, la distance de de entreles points,de l'paisseur de la gomtrie de l'espace et entredeuxplatsmincessouds. Commeil a t indiqu pcdemment, joints souds point font I'objet de les par nombreuses applications dansles industries automobiles aronautiques. en rsulte et Il desdiffrentscasde chargement correspondent tlpes d'application leurs qui aux et fonctionnalits.

L0

Chapitre l. Introduction gnrale

J
b)

j
c)

par Figure 1-1. Type de sollicitationdesjoints souds point a) cisaillement b) tensionbi-axialec) pelage.

variablesdans le Gnralement, points soudssont soumis des chargements les qui du temps. Ces chargements dpendent temps peuventtre de type impact ou priodique, seposealorsle problme la il de priodique. Dansle casd'un chargement prvisionde la durede la vie en fatigue. priodiques comprenant forceset desmoments des Il estvidentque deschargements par particulirement comportement des le mcanique soudures point. Il conditionnent pratiques. convientdoncde le vrifier pourdesapplications par Selonle chargement, naturedesjoints et les diffrentstypes de joints souds la par point ont t classifis. I-es diffrentstypes de joints souds point sont appels joints en cisaillement, tension et en bi-axiale de pelage @gure l-1). par sont les plus tudis Gnralement, joints souds point de type cisaillement les pour l'valuationde la durede la vie en fatigue. L'influencede la gomtrie des Les sont prisesen considration. soudages soudures la position despoints souds et diffrentselonle nombre par points en cisaillement un comportement mcanique ont (Tractionet Compression). et de points,la localisation lesconditions chargement du

1.0

Clnpitre 1. Introduction gnrale

la Par consquent, prvisionde la due de vie en fatiguede ce type de joint est ts prsentantdiffrents tlpes de importante pour des applicationstechnologiques mcaniquedu point soud en configurationMono Le chargement. comportement priodiquebassur Point, Triple Pointset Multple Points soumis un chargement I'approchevolumtriquesera tudi dans la prsentetude par I'approchede la fatigue. nousnousintressons la prvisionde la durede vie en fatiguedes Danscettetude, joints souds point en cisaillement. par et de du Dansle Chapitre2, le concept facteurde concentration contraintes du facteur de Les mthodes la prvisionde la durede vie en de rduction fatilue sonttudis. est La en fatigue sont prsentes. restrictiond'utilisation des diffrentesmthodes prsente ce chapitre. dans graphique volumtrique prsente.Laprsentation est Dansle Chapitre3, la mthode et I'utilisation de la courbe biJogarithmique sont montres.Une formulation Dans cette formulation, nous amliorede I'approchevolumtriqueest propose. permettons calculerla contrainteeffective,la distanceeffectiveet le facteurde de 3-2. proposes dansle Tableau rduction fatiguepar les quations en Les est Dans le Chapitre4, la bibliographiedesjoints souds prsente. diffrentes mthodes utilisespour I'analyseet la prvisionde la durede vie en fatigue sont dcrites. par point en pour les joints souds la I-e Chapitre5 prsente partie exprimentale et les du mcanique matriau, conditionsd'essais les I-e cisaillement. comportement rsultats sontprsents. pour la prvisionde volumtrique I'approche 6, Dansle Chapitre nousavonsappliqu la durede vie en fatigue.La directionet la contrainted'ouvertureen fatigue sont effectiveet le pour obtenirla contrainte effective,la distance utilises respectivement facteurde rductionen fatigue. critiquesest aussi des au La distributionde contraintes voisinage points de soudure et prsente danslescaslastique lasto-plastique. 7. dansle Chapitre gnrale notretudeestprsente de La conclusion

CHAPITRE 2

Chapitre 2
de 2. lN'fltllrode prdiction de la dure de vie en fatigue de en prsence concentrationde contraintes
2.0 Introduction gnrale
dans de gomtriques conduisent desconcentrations contraintes Les discontinuits dansles de de La de les composants structure. prsence cesconcentrations contraintes On invitableet indsirable. peutobserver est mcaniques un phnomne composants mcaniques' danstous les composants indsirables de les concentrations contraintes Toutefois les meilleures ralisations technologiquesprsentent encore des concentrationsde contraintes dues aux changementsgomtriqueset aux dues I'usinage. de modifications surface et dansla micro structured'inclusionset de dfautsdansles soudures La prsence lors indsirables de de conduitaussi desconcentrations contrainte dansles moulages la ralisation. En outre, I'environnementpeut augmenterla dangerositdes de concentrations conffaintes. de cyclique localis dans une zone de concentration contraintes Un chargement I'amoraged'une fissurede fatigue. L'tudedes localesen fond d'entaille entrane est mentionns, donc d'une contraintesen fond d'entaille,pour les chargements pratiques. pour desapplications essentielle importance Le conceptde fatigue des matriauxa t proposavantle 19mesicle,mais la Des auparavant. avait t dj dcouverte de prsence la concentration contrainte de pour et en de expressions la rsistance fatigue avaient t proposes appliques priodiques.On a Elles s'appliquent des chargements rsoudreces problmes. observle premierproblmede fatigue dansles axesde wagonsde chemin de fer rpt. soumis un chargement que par Rankinelll a tI'un despremiers tudierlesruptures fatigueet suggrer et dansleursconceptions privilgierunerelationde continuit devraient les ingnieurs d'viterles concentrations D'autrepart,il a propos technologiques. les constructions utiliss par I'industrie mcaniques des contraintessvresdans les composants

ferroviaire. Des accidentssrieux d'avion [2] ont prouv que plus de 50Zo des incidentsavaientt produitsen raison de I'amorage d'une fissurede fatigue dans deszonesde concentration levede contraintes, exemple,les trous de boulons, par lescrous les filets d'attache. et Durant la priode 1927-1981,des tudes statistiques indiquent qu'aucune tude complte n'a t ralise sur les phnomnes fatigue dans les composants de mcaniques queles trous,les filets et lesjoints souds. tels Is exigences accrues scuritdansles avionset dansles vhiculesde transport de imposentmaintenant d'analyserI'amorageet la propagation fissuresde fatigue des (concept tolrance dommage). de au

2.L Distribution descontraintesen fond d'entaille

<< sont gomtriques appels entailles>> Dansles tudesde fatigue,les discontinuits . la par et canctises trois paramtres, longueura, le rayonp et I'angle d'entarllety les les tous les dfauts, imperfections, Cettedfinition permetde dcriresimplement et de variations gomtrie lesjoints souds. Des contraintes localement une lvationde contraintes. Ces entaillesconduisent et par sont dues aussi aux excentricits rapport I'axe de chargement excessives danslesjoints souds. et installations de dfauts de proviennent mauvaises dansdes tudesde la Cettedfinition du conceptd'entailleest utlise couramment l'tude de I'influence des paramtresde I'entaille est fatigue. Par consquent, prs des zonescritiques' pour connatreles distributionsdes contraintes ncessaire pour descassimplesavecun comportement L'tudede la distributiondescontraintes analytiquement, exemplele casd'un petit du Pil lastique matriaupeuttreralise (Figure2-l).L'ttde du problmementionn trou dansune tle de grandedimension et de de fournit une bonneidede phnomne concentration contraintes de ci-dessus des la redistribution chargements. autourd'un petit trou dansune tle Dansla Figure 2-l,ladistribution de contraintes Pour obtenirla distributiondesefforts autourdu de grandedimensionestprsente. trou dans une plaque, on se base sur I'approcheclassiquede la rsolution des et la d'quilibre, relationcontrainte-dplacement I-es d'lasticit. quations problmes dansle casde l'lasticitlinaire. l'quationde compatibilitdoiventtre satisfaites ci-dessous: sontcrites d'quilibre polaires, quations les En coordonnes

au de Figure2-1. Illustration de la concentration contraintes bord d'un trou dansuneplaque.

Dansla Figure 2-1, a, r , 0 et os reprsentent rayondu trou, les composantes le des coordonnes polaireset la contrainte globalerespectivement. En I'absence forcesgravitationnelles, equations de (2-1) et (2-2) sont satisfaites les partir de la fonction des contraintes, dont drivent les contraintes radiales, circonfrentielles scissions. et

=o +*!u!:0 +o,-o,+F, r r r

(2-r)
(2-2)

!u{**+*2'," +4 =o
r0 r r

En I'absence forcesmassiques et ), les relations(2-l) et (2-2) dewonttre de (F, satisfaites par la fonction QQ,O) exprime partir des contraintes radiales, circonfrentielles de scissions. et Elles sontprsentes ci-dessous : -r
f =---r-J

ra | a2

rr'r"Ae'

(2-3)

2 o, =j ctr-

(2-4)

=- a (ra\ ",' arl;ao )

coordonnespolaires par:
=-

(2-s)

La relation dplacement-dformationdans le cas mentionn peut tre reprsenteen

u
dr

(2-6)
u

t'=;*|

lAy

(2-7) (2-8)

v lu v T"=t*; ae-;

o, u , v , ,, ee et y sontrespectivement dplacements la directionx et y, les dans les dformations radiales, circonfrentielles de cisaillement. et Finalement l'quation compatibilit de s'crit:

| &*+e-*?b-!L=!'r,, * r2 r,, ' dr' r' 0' r r r r r r0

e 9

Q-9)

2.1

Chapitre2. Mthodedeprvisiondeladuredevieen fafieue

polaires,l'quationde compatibilitavecintroductionde la fonction En coordonnes de contraintes pour ce problme d'lasticit bi-dimensionnelle classique et par est sans axisymtrique forcesgravitationnelles reprsente la relationsuivante: (2-10*) et convenable l'applicationdesconditions une En considrant fonctionde contraintes la aux limites appropries, solution analytiqueexige des conditionsde symtrie peut bi-dimensionnelle trecrit: en Le ci-dessous. casdu problme lasticit dcrites

= -5). + za] o, o,l(, (,. +).",

(2-r0)

= .5)-. ", io,l(' (' T)*"'1t

,,,=-Io,('-T.+)sin2o

(2-rr)

(2-12)

La distributionde la conffainteautourd'un petit trou dansune gtandetle soumise dansla Figure2-2. Pourr=Q,0=90",la bi-axialest prsente un chargement diminuer. maximaleo, atteintla valeur 3o, Pourensuite circonfrentielle contrainte La contrainte radiale augmentepour atteindre une valeur maximale et ensuite zro. diminuerjusqu' autourd'un petit trou, dansune grandetle, peut tre La distributiondes contraintes diffrents.Cela conduit des rsultats modifie par applicationde chargements et analytiquesdans le cas axisymtrique lastiquelinaire bi-dimensionnel. Des autres cas sont montrs dans la Figure 2-2. Ces chargementsde traction et de maximale conduisent une valeur de la contraintecirconfrentielle compression (2or du lieu 4o, ) de de de Is phnomnes concentration contrainteen fond d'entaille dpendent la gomtried'entailleet du type de chargement [3-6]. L'influence de ces paramtres attnuentcet effet de Les sera analysepour le dimensionnement. concepteurs du par de concentration contrainte unemodificationde la gomtrie composant.

10

a)

b)
Figure Disuibution la contrainte 2-2 de circonfrentielle ar(r,0=1T12)autour d'un troudansune grande a) Chargement tle verticalentraction b) Chargement vertical compression. en L'analyse de la distribution des contraintes n'est pas chose une facile pour toutes les gomtries.L'utilisation de rsultatsanalytiquespour des distributions des contraintes en fond d'entaille est limite quelquesmodlesclassiques. L'utilisation des rsultats obtenus par la mthode des Elments Finis est pratique et efficace pour analyser les distributions des contraintespar superposition pour toutes les gomtriescomplexes.

. Exemple
Nousavonspris une gomtrie pour un petit trou dansunegrande plaque(la solution analytique pour ce casexiste) partir d'un codedu calcul (ANSYS).Les Figures2-3 2-4 reprsentent distribution des contraintes la autour d'un petit trou dans une grande tle soumise diffrentschargements la mthode Elments par des Finis. Les Figures 2-3 et 2-4, prsententles distributions des contraintessuivant les directionsx et y pour le casvoqu cettedistributiona t calcule la mthode par ; desElments Finis.

ll

2.1

Chapitre2. Mtlade deprvisionde Ia durede vie en fatieue

Dans les Figure 2-3a et 2-4a,la distribution des contraintes est diffrente. Elle est obtenue pour une mme gomtrie mais avec un chargement diffrent. La comparaisonde rsultatsobtenusmontre que la contrainte maximale selon la direction gale 2o, et 4or. x est respectivement

a)

b)
tle par la mthode autourd'un trou circulairedansune grande Figure2-3 Distributionde conhaintes y x Finis @gure 2-2a) a) Composanede contrainte b) Composante de contrainte. desElments

12

2.1

Chapitre 2. Mthode de prvision de la dure de vie en .fatieue

a)

b)
Figure2-4. Distributionde contraintes autourd'un trou elliptiquedansune grandetle obtenue la par mthodedes ElmentsFinis (Figure 2-2b) a) Composante de contraintes Composante de x y b) contraintes.

13

2.2.1.Distribution descontraintesdansle caslastique

La connaissance la distribution des contraintes fond d'entaille est ncessaire de en pour dterminerle comportement mcaniquedes structures.Des mthodesanalytiquesrpour obtenircette distributiondes contraintes fond d'entaille,ont t proposes plusieurs en par auteurs et de Linaure [7-8, 10-11,14,16]. Creager Parisont utilisla solution la Mcanique de la Rupture(MLR) aprsI'avoir corriget7-81.En coordonnes polaires,dansla Figure 25, ils proposent la reprsentation la contrainte que de d'ouverture s'critcommesuit:

o,=#*(*)*'(i).h*'(
panir du fond d'entaille. Le facteur d'intensit de contraintesK7 , pour une fissure de longueur a , s'cit comme:

(2-13)

o p est le rayon d'entaille, K7 est le facteur d'intensit de contraintesen mode de rupture I et(r,0) sont les coordonnes polaires partir d'une origine dcale la distancepl2 de

K, = FoorJ
contrainte globale. La contrainte d'ouverture maximale en fond d'entaille vaut:

(2-14)

o F" et os sont respectivementle facteur de correction gomtriquenon dimensionnel et la

of* = o, +Zf"or.rll
"\p
En d'autrestermes,le facteur de concentrationde contrainte s'crit comme:

(2-rs)
(2-r6)

k ,= L + 2 F , ^ l z lp

ln

Figure2-5.Schma d'uneentailleet lescoordonnes polaires correspondantes.

La connaissance la distributiondescontraintes la pointede I'entaillepermetd'valuerle de facteur d'intensit contraintes de qui d'entaille caractrise partiede la distribution[9-l l]. une L'quation Paris-Creager donnela distributiondes contraintes pour des entaillesprofondes

t4

2.2.1

Chnpitre 2. Mthode de prvision de Ia dure de vie en .fatipue

quandI'angle d'entaille est gal zro @gure 2.5). D'autres expressions uniquement pour cettedistribution.En considrantla solution suggres ont analytiques t galement du champ des contraintesde Neuber pour les entailles ll2], la contrainte d'ouverture par du au s'exprime voisinage fondd'entaille [13] : o . . = o**
OU' x

\!p+ax

(2-r7)

et d.a* sont respectivementla distance mesure panir du fond d'entaille et la.valeur

maximaledeoo. Chen et Pan [17] ont proposune quation lgrementdiffrente: lo oo=o**1ffi

(2-18)

En comparant la relation (2-18) avec la solution exacte obtenuepar Elments Finis, pour une entaille aigue, l'quation de Chen-Pan donne une valeur de 5 L5Vomoins leve. La valeur de I'erreur devient plus importante si I'entaille devient plus mousse. est Glinka et Newport [4] ont conclu que l'quation de Creager-Paris raisonnablementexacte pour les entailles aigues. Pourtant de, prsentedes valeurs sous estimes pour les entailles moussest151. Glinka et Newport ont suggr une distribution des contraintes d'entaille diffrente pour les entailles mousses et les entailles aigues sous la forme d'une approximation polynomiale. Is deux expressionsde la contrainte d'ouverture se prsentent ci-dessous:

r o, =n,orlr z szo(L).*r(d)" - ono, .t *t[;)' (;)'

k,<4'5

]
k, > 4'5

- rrr(;) o, =r,,o- o.zzs(LJ" . t rr( -, *r(;)'f rlt ;)"

globale, la distancedans la direction x et le rayon d'entaille.

(2-re)
o, ft, ,o s , x et p sont respectivementle facteur de concentrationdes contraintes,la contrainte

pour plusieursconfigurationsd'entaille[4]. (2-19) sont desrelations exactes Les expressions pour n'est pas satisfaisante Kujawski [11] a montr que I'approximation dcrite prcdemment une entaille avec une concentrationde contraintefaible k,33. Il a donc proposune mthode pour amliorer la relation de Creger-Parispar I'application d'un facteur de correction /:

++ o,= rylO * rLf'5(1rir"f

r5

(2-20)

2.2.1

Chapitre 2. Mthode de prvision de Ia dure de vie en fatigue

(2-2r)
Dans des manuels,le facteur de correction gomtriquenon dimensionnel, pour diffrent type de fissure, peut tre obtenu partir d'une analyse de la contrainte. Il faut remarquer que le facteur de concentration des contraintes est relatif la contrainte dans la section moyenne danscesmanuels. La formule (2-L6) a t tudrepour des gomtriesd'entailles diffrentes par la mthode des ElmentsFinis. Les relations analytiques dcrites prcdemment sont des mthodes approximatives pour connatrela distribution des contraintesau voisinage de I'entaille ou de la fissure. La mthode des Elments Finis peut tre aussi applique pour vrifier cette distribution de contrainte. Dans la Figure 2-6,1agomtrietudieest psente.

q
Figure 2-6. Schma d'une entaillemousse dansune plaquesoumise une tractionperpendiculaire plan au d'entaille.

. Exemple
Nous avonspris la gomtrie pour uneentaillemousse suivante dansune plaquesoumise une traction perpendiculaire plan d'entaille par la mthodedes ElmentsFinis. I-es au pour hypothses lesmodles d'Elments Finissontdonnes ci-dessous: o Modulede Young
o e o Coefficientde Poisson

E =2L0 GPa
v =0.30

Hypothsed'un comportementlastiquelinaire Contrainte globale 6e =l.o MPa

t6

2.2.1

Chnpitre 2. Mthode de orvision de Ia dure de vie en .fatieue

Hypothse de contraintesplanes Symtrie double pour le modle desElments Finis

Longueurde plaque Largeurde Plaque . de paisseur la plaque

L =3.20m W =2.40m B =0.001m

Dans la Figure 2-7, nous avons prsent les conditions aux limites, la configuration du maillage et la distribution de contrainte pour le rapport I I p =1gg . On peut constater dans la Figure 2-8 que les mthodes de Glinka, Kujawski et Chen sont pour connatre la concentration de contrainte. La mthode de Cregar-Paris exacteset adaptes est plus exacte pour des concentrationsde contraintestrs leves,ce qui se produit dans des Le entailles aigues ou des fissures mousses. choix de la mthode dpend donc des valeurs du facteur de concentration de contraintes.

a)
EWEE Fs.rcrpldt pF1 av =l Drffiao a! * =aDo w z-ffitaB Er2 il-1 'Dr9?.993 rE! r33S.8ll .a =50,I84 z-lwtl EF3 .DrrFl5.3{A .rr 439-O5t . =a2.t9 z-tulttR

b)

t7

2.2.1

Chanitre 2- Mthode de nrvision de Ia dure de vie en tbtipue

c)
Figure 2-7. Conditions aux limites, configuration des maillages et distribution des contraintes pour le modle slectionndans la Figure 2-6, y compris I I p = 100.

La mthodedu Glinka donnedes valeursprochesdes rsultatsnumriques, mais avec un domainede validit rduit. Elle est moins exactequandon s'loignedu fond d'entaille.I-es mthodes Chenet de Kujawski et Neubersous-estiment distributiondescontraintes. de la La pourlesentailles mthode plusconvenable la aigues la mthode Creager. est de a)n
l0 8 --El,ntnis --O-Nerber --A-Cilhka -jX- Kujavtski --X-Chen -0-Geager

q.
b

18

2.2.1

Chapitre 2. Mthode de prvision de Ia dure de vie en -fatipue

b) ro
t4 T2
o
>l

l0 8 6 4 2 0

o
b

-{-Elerrentfinis --O-Neuber -*-Glinka -iK-Kujawski --X-Chen --O-Geager

o,

o.

Ilp47

c) zo
18 16 l4 *12

-{-O----iK--X--O-

Elementfinis Nerer Clinka Kqiaunki Chen Geapr

Iro
b8
6,

6 4 2 0

llFts d) zt
22 20 l8 16

-{-Elenentfinis --O- Neuber ---Glinka -*-Itujaunki -*-Clrcn --O-Geager

>r2 bto
8 6 4 2 0 0 llp=tgg

dr4

xlp

de Finis, dessolutions Neuber,Glinka,Kujawski,Chen des avecla mthode lments Figure2-8. Comparaison pour valeurs rapportI I p fl = 2a). du et Paris-Creager diffrentes

19

2.2.1

de de Chnpitre2. Mthode orvision la durede vie en .fatisue

pour l'valuation de la distribution des contraintes sont relatives aux Les mthodesprsentes Dans les cas rels, des cas lastiques.Il faut remarquer que ce sont des solutions approchees. effets de la plasticit apparaissenten fond d'entaille. C'est d au comportement mcanique la des matriaux et au fait que la contrainte en fond d'entaille dpasse limite d'lasticit. I-e Tableau 2-1, rsume les diffrentes mthodes dcrites et la distribution de la contrainte d'ouverture en fonction de la distance.

pour aigues des donnant distribution contraintes desentailles la analytiques des Tableau Rsum mthodes 2-1.
et mousss fonction de la distance. en METTIODES Timoshenko LE MATHEMATIOI.]ES POI.JR COMPOSANTDE LA CONTRAINTE D'OWERTURE EXPRESSIONS

(lesr) tl8l
Neuber

o =%. +(,. .i(r-;)o [' ;)"

(re61) tlel
CreagerParis(1967)

= oyy^ffi "^

o_.,.

tlsl
Chen-Pan

-*(*)."'(T) .'-(f""g) = *ffi"o"(tt')

E 8x \ p+
\s4 .s
k, > 4'5

(1e78) t17l
GlinkaNewport (1e87) tl4l

= o_^. a*,

oo

= ho,lr-n3o(!-).r *{z)" -oror .'*'[;)'] (;)' -o,r'[r)'] =r,o,[r-o r"(o) -' rr(;)., zs(-r)"

o,

(leel) url

Kujawski

L.o.z
p

L, o.z
p
Xu KeweinJiawen

m = constante

(re92) l20l

o, = or + oil + onri I et IITraction,III Compression (le9s) [21]

BhatacharKumar

6 o,=o**|tr, o(x<f,r*,-t>
(w:o) M= < x <(w - a) o,,= ((w: ")-"), \ " I [ 2 2 I
M = Moment I = Moment Inertie

20

2.2.1
METHODES

Clwpitre2. Mthodede orvisionde Ia durede vie en fartgue

DPRESSIONS MATTTEMATIQIJES POUR LE COMPOSANT DE LA CONTRAINTE D'OIMRTIIRE

(rees) l22l

R. X. Xu

oo

=,,ry(r.

)r['. [')'],

ooo5ft, a=0 425+

f-

= facteur decorrectionde gomtrtepourles prouvettessoumises unetraction

2l

2.2.2.Distribution des contraintes en fond d'entaille dans le cas lasto-plastique

La distribution des contraintes en fond d'entaille change selon la gomtrie et le mode de chargement. Si les valeurs des contraintes sont suprieures la limite d'lasticit, une relaxation plastique intervient en fond d'entaille. Dans les matriaux mtalliques, on observeun comportementde nature lasto-plastique.En une zone en plastifie et des contraintes lasto-plastiquesdoivent tre pris en consquence, considration Dans le cas des matriaux fragiles, tel que I'acier haute rsistance,le verre, la cramique et ne la fonte la zone plastifie et les contrainteslasto-plastiques serontpas considres. L'objectif principal de cette tude est concentr sur les aciers. Ces aciers sont gnralement considrscomme ductiles en raison de leur composition chimique. L'tude de la zone plastifie dans le cas d'application de chargementspriodiques est trs importante pour valuer la dure de vie en fatigue. De multiples travaux ont t effectus dansdesconditionslasto-plastiques 121,23,241. pour tudierla distributiondescontraintes Les mthodes suivantes sont donnes titre d'exemple et montrent l'volution de la en distribution lasto-plastique fond d' entaille. Wilshaw l23l a propos une relation poru cette distribution dans un matriau lastique parfaitement plastique. L'volution de la distribution lasto-plastiqueest divise en quatre zonesdistinctes. I-es relations proposeset les limites de ces zones distinctives sont exclusivement obtenues pour des entailles sollicites en flexion.

Zone I: distribution logarithmique des contraintes:

o,(x)= o,[r.*(t.;)}

0( x( r,

Q-22)

et la oo" et p sontrespectivement limite d'lasticit le rayond'entaille.

: o ZoneII: la contrainte resteconstante d'ouverture

o,(x)=o,[r*,n[t+)}rr 3 xl

ro

Q-23)

22

2.2'2

Clnpitre 2. Mthode de prvision de Ia dure de vie en fatisue

En raison de la continuit des contraintes pour .r = ru I'interface des zones I et II, la

taille de lazone plastifie s'crit :

,"=plel;) -rl .(
) o 14 et p sont respectivementI'angle et le rayon d'entaille. o Zonefr:

( ("-w\

(2-24)

(2-2s)

o o, , ko et L sont respectivementla contrainte globale, le facteur de concentration de contraintes lasto-plastiques et la distance partir d'un fond d'entaille qui

correspond la fibre neutre o la contrainte s'annule. Le facteur de concentration de contraintes lasto-plastiques s'exprime selon la mthode de Wilshaw [23] par :

o,=[,+h(r .i)]

(2-26)

Figure 2-9' Evolution descontraintes lasto-plastiques le ligamentd'une prouvette sur soumise une flexion selonWilshaw 231

23

2.2.2

Clwpitre 2. Mthode de prvision de la dure de vie en fartpue

o ZoneIV: Partieen compression la flexion lastique. due Dansla Figure2-9,|a contrainte quandon s'loignedu fond d'entaille,et atteintla augmente valeur maximale pour.r= ry. Ensuite, elle demeure constanter= r dans I'intervalle rr 3 x < 7. Aprs,les contraintes puis diminuentplastiquement lastiquement. L'quation (225) estla mmequecelle de Neuberqui confirmele comportement lastique dansla zoneIII. En raisondu chargement flexion, on observe zoneen compression. en une La distinctiondeszonesen compression tractionestdonne I'axe neutre(x = rn). La par et compression devientmaximale pour x = D. Bhattacharya lasto-plastiqr:es en [21] a proposune mthodepour valuerles contraintes fond d'entaille. Cette mthode est utilise pour une entaille soumise la flexion. La distributiondescontraintes s'exprimepar :
6r=O,+oil+Ou+ON+Ov

(2-27)

o or , a n , a n1, ory sont des contraintesen traction. ov estune contrainte en compression. La formule donnant la distribution de ceffe contrainte s'crit @gure 2-10):

o,=o,lr*r"(r.;)} 0 3 x 3 r ,
6tt =

(2-28) (2-2e) (2-30) (2-3r) (2-32)

",lr*-il'
o,o,fft,

rrlxlrp

oat =

rO1a1l

\ =ll oru M ( w - a -')' z M(w ) o'=71- a -'I z

0(x< (w-a)12 (w-a)|2< x<(w-a)

o o, , P, V , k,,or, M, I,'W, a, r, et ro sont respectivementlimite d'lasticit, rayon la le d'entaille,l'angled'entaille,le facteurde concentration contraintes de lastiques, contrainte la globale,le momentde flexion , le momentd'inertie,la longueur, profondeurd'entallle,la la qui distance correspond lazone plastifie la distance et correspondantla fin de la zoneII.

2.2.2

Chapitre2. Mthodedeprvision de la durede vie en fatisue

=Or+Oil+Ou+ON+Ov

N ct
!>
I

^
cl

F
d'une prouvette Figure2-10.Evolutiondescontraintes elasto-plastiques le ligament sur soumise une flexion selonBhattacharya [21].

pour un petit coulement prsente. est Dansla Figure2-Ll,l'volution descontraintes Les (2-33)et (2-37)sontla solution gnrale la mthode Bhattacharya. de de quations on = o, + oil + orrr+ oN Q-33)

o, dr, tr,o,osont respectivement contraintes traction. o'est une contrainte des en en : compres sion.I-esquatreexpressions expriment s'

o, = o,lr**(t.;)],
o, = k,o,

o1x1 ry

(2-34)

p p+4x

-r) rr3x=t=LOQ','
(w-a)t2 o<x< x @- a)t2s <(w a)

(2-3s)
(2-36)

or,=+(T-,), =+(ry -,), orv

(2-37)

la o o, , p , V , k,,or, M, f, W, a, r, et ? sont respectivement limite d'lasticit,le rayon d'entaille, l'angle d'entaille, le facteur de concentration de contraintes lastiques , la contrainte globale, le moment de flexion, le moment d'inertie, la longueur, la profondeur d'entaille, la distance qui correspond la zone plastifie et la distance correspondant la fin de la zone II.

25

2.2.2

Clupitre2. Mthode prvision Ia duredevie en fatigue de de Figure 2-11 y compris les

L'volution des contraintesglobales est illustre dans la contraintes tractionet en compression. en

oo=Q*ot*onlotv

N l>
I

^
ctl I f>

Figure 2-11. Evolution des contraintes lasto-plastiques long du ligamentd'une prouvette le soumise une flexion pour le problme plastifie petitetaille selonBhattacharya d'une zone de [21]. Comme il a t dit prcdemment, la relation entre la contrainte locale et la dformation en

fond entailleest essentielle. Dans I'approche locale,la dtermination la relation entrela de contrainteglobaleet la dformation fond d'entailleest ncessaire en pour la prvisionde la durede vie en fatigue. Pour les entaillescourtes,une autrerelation est propose M. par Zheng-E. Niemi et Shi-Puls124-261. gomtrie cesentaillescourtesest prsente La de dans la Figwe2-12.

+ +

(x
+

Q-+
6-

+5

-o ^

+ +

x
Figure2-12.Schma d'uneentaillecourteet la contraine bactionglobale en applique, d,

26

2.2.2

de de de Clnoitre2. Mthode prvision ln dure vieen fartoue

uniaxialor. Si la et L'entaille courte a une forme semi-elliptique soumise un chargement profondeurd'entaille a est petite par rapport la largeur W, la contraintelastque au base sur en directement utilisant une mthode voisinage fond d'entaillepeut tre obtenue du la thoriede l'lasticit.Quandla contrainteatteintla limite d'lasticit,la valeurmaximale descontraintes tractiono** s'exprimepar l24l: en

o*,=o)t+2|!l

r -\o'l

"L \p) I|

(2-37)

la o o,** ,or, p et d sont respectivement contraintemaximale,la contrainteglobale,le d'entaille. rayonet la profondeur Dans l'quation (2-37), la deuxime partie peut tre considrecomme le facteur de lastique. concentration contrainte de

k.=r+219l ' lp)

\0.5

t2-38)

pour en lastiques gnral de L'quation(2-38) donnele facteurde concentration contraintes la Dansle caslastique, toutesles entaillescourtesy comprispour la forme semi-elliptique. contraintede traction selon la direction y dans la rgion prs du fond entaille peut tre par obtenueapproximativement la solution de Neuber.Cette solution est prochede celle pour le casd'unezoneplastifiede petitetaille parBhattacharyhal2ll. obtenue prs En conclusion, l'volution descontrainteslasto-plastiques du fond d'entailleprsente : les caractristiques suivantes o La contraintemaximale,dansI'analyselasto-plastique, n'est pas en fond d'entaille mais une certaine distance. o par L'volution des contraintesdans la zone plastifieepeut tre reprsente une fonctionlogarithmique. o Pour une zoneplastifiede petitetaille, la contrainte est ma,rimale plus faible, et est lastiques. de par donne le facteurde concentration contraintes o Dansle cas gnral maximalen'est pas la d'une zoneplastifieimportante, contrainte dansunezoneo elle estconstante. en seulement localise un point,maisrepartie o Pour les fissuresprofondes, plastique limit peut tre aussi I'analysede l'coulement utilise [25]. pourune approche locale. et commeapproximatives utilises Cesmthodes sontconsidres

27

2.2.3.MthodedesElmentsFinis dansle caslasto-plastique

L'utilisation de la mthodedes ElmentsFinis n'est pas limite aux tudesdes matriaux lastiques linaires. Elle est aussifrquemment pour utilisedansles tudes lasto-plastiques par analyser exemple,la zone plastiqueau voisinage fond d'entaille dansdes matriaux du d'acqurirpralablement connaissances ductiles.Cela demande les thoriques fondamentales utilises dansla mthode Elments Finis [27]. des L'expressionde la relation entre les contrainteset les dformationss'crit corru,rre cidessous :

{o}= [a{e'}

(2-3e)

o, {o}, [O]"t{e"'} sont respectivementles vecteursdes contraintes,la matrice de la rigidit

et le vecteurdesdformations lastiques. Les expressions vecteurs contraintes de dformations : des de et sont {o}={o{e}={e*

6o
,

oz
eu

t,
,

tyz T*I'
eyz t*}'

(240) (241)

par Les expressions contraintes des sontprsentes le vecteurou le tenseur contraintes des de la Figure 2-13. Les tenseurs symtriques contraintes dformationssont donnspar des et (242): l'quation

lo* r, '*1 lt* lol=lr, 6o ,,, l,le)=| r,o L"- r* o. ) leo

o o

t"l
e.)

'*1
(2-42)

r,

):", (,""
,6u O ro.----.,"' i v

*?" "l
/ o" a o_

o,
/oo
I

Ortn-t_

Figure2-13.Composants tenseur containtes. du des

2.2.3

de Chnpitre2. Mthode prvisionde la durede vie en fafiguz

En gnrale,les dformations s'crivent:

{e}= {e"'}+{e'o} o, {e}, {r'}

(2-43)

des les .t {e'o} sont respectivement vecteurdes dformations, dformations

thermiques. lastiques desdformations et sont: thermiques I-esvecteurs dformations des {t*}= LT {ao Z=1-Tsa, a, qa o o o}' (2-44)

qn , dyyet ausont le changement temprature, coefficientde le de

x dansles directions , y et z respectivement. dilatationthermique peut tre dduite des et L'expressionde la relation entre les dformations les contraintes (2-39)et (2-43): quations

{e}= {e'}+ {r" }= [af'{o}+ {e"}

peuttrecritecommeci-dessous : La matricede raideurdansle caslastique
1 -vry -ve

(2-4s)

E,
-vyx Ey -vu

E,
I Ey -v4

E,
-Vyz Ey 1

[of'=

E,

E,

E"

0
I Gry

0 L
Gv"

0 0 1
Go

les o, E,,Eret E. sont respectivement modulesde Young dans les directionsx, y et z. x, dans de sontles modules cisaillement les directions, y et z.vo,vn,vo sont G,y,Gyr,G*. de les coefficients Poisson. de Les caractristiques symtrie de la matrice des tenseursde contraintessimplifient la relation(2-46)selon: b --v r, . vo - vu , E , Er ' E , E r' vr, =vo Ey E, ( 2 47 )

lastique isotropeconduit du L'utilisation de l'quation,(2-47) et I'hypothse comportement la relationsuivante:

2.2.3

de Clnpitre 2. Mthode prvisionde la durede vie en farteue

E r = E r = E, , vr,y =v' y ,xt = v y , ' 'ry

(248)

(2-49). linaireisoffope,quation Cecicondui l a ll o i de Hooke pourun matriau lastique uit r de E - l- v E -l v E 0 0 0

-- v : l

-v0
T E -v0

e# o ea , r, exz

E -- v t E 0 0 0

1 T-

0 'l
0 0 ll"0 I G 0

Io
E 0l

o o:|,z
0 0ll"
GJ

+LT

(2-4e)

0 0

0llz + 1l 0
(2-s0)

le le o, E, v,G et a sontrespectivement modulede Young, le coefficientde Poisson, modulede cisaillement le coefficientde dilatationthermique. et L'utilisation du principedestravauxvirtuelss'crit:

&l =N
gnral, peut crire : on &J = &Jr* M, N =6Vr+Nr+N,

o, U et V sont respectivementl'nergie de dformation et les travaux extrieurs.Dans le cas

(2-sr)
(2-s2) (2-s3)

tu, = !{at[r}av'
v

(2-39)et (2-43)dansl'quation(2-53)donne: L'applicationdesquations

N,=!<{a'Y [o]t"]-{a'f lol{e'^\av

v

(2-s4)

La relationentreles dformations dplacements nuds et des s'exprime:

{"}=[r[u]

<2-ss)

le des la o, {e}, [n]et{r} sontrespectvement vecteur dformations, matricede dformationLa dplacement et nodaux. base lesfonctionsd'interpolation le vecteurdesdplacements sur substitution l'quation(2-55)dansl'quation(2-54)donnela relation(2-56): de

N, ={aY trltrln > - {'Y <ftafl {aI Bffin[nlat

vv

(2-57)

: &Jrest l'nergielastique la fondationlastique due

30

2.2.3

Chnnitre 2. Mthode de prvision de la dure de vie en fatioue

i

tuz=I{n"y {"1^,
Al

(2-s8)

o, 6wn, o et A, sont respectivement dplacement perpendiculaire la surface,la le contrainteet la surface correspondante. dplacement perpendiculaire la conffainte Le et s'crivent dans cas: ce

{.,}=hr,J[r]
{o}=k{r^}

(2-58.)

(2-se)

(2-58-)et (2-59) dans l'quation(2-58), nous avonsla relation En mettantles quations suivante :

tu, ={6"Y [iv, ]'[iv,]dA){ G

Al

(2-60)

de [N,] sontla rigidit de la fondationqui s'exprimepar la force sur les surfaces la "t fondationet la matricede fonctiond'interpolation pour la fondation. A prsent, nousavonsI'expression l'nergiepotentielle: de

o, k

N,=-{tYrrJtnJ' S {"} tiv1an

v

(2-6r) (2-62)
(2-63)

N, ={a,f [[N,]{r}ae,
Ao

N,={arY {ri!

o, tl,, Nret dV,sontrespectivement nergies potentielles les bases les effetsd'inertie, sur la pression applique la surface , lesforcesnodales sur A" extrieures. quantits{f}, p et Les {ryI correspondent respectivement vecteur la pression, la densit aux vecteurs au de et des

qui forcesextrieures s'appliquent nuds. aux (2-57),(2-58)et (2-6I) (2-63)donneunerelation: L'utilisation desquations

+{af {,Y( tN.] JtaItpfrlav1{a}-{aY @l[nlav 4JIN"f d,e;{}= nf

-{aYpllrtf Wl7fu}*{alrjhr" f {pIa,+{af {rf}

L'quation(2-64)va tresimplifiepar l'liminationdu vecteur{}r:

(2-64)

* rlr"* lrI \t"t- {r:}=[u"H"]{ry}*{,y } I

OU,

(2-6s)

31

2.2.3

de d.e de Chnpitre2. Mthode orvision la d.ure vie en fafieue

lr,l= nT[ona)av
il, lx l=r J[iv,J'lr,rl
AI

(2-66) (2-67) (2-68)

lr:^l= lol'"}u [lnf v

lu"l=olllv| Wlav v
u=VIul
.. "-,

(2-6e)
(2-70)

{r} lr:')=![x,Y ae,

AP

(2-7r)

linaireisotropen'existe pas toujours.D'autres comportements lastique I-e comportement (Figure2-14).Dansla Figure2-14,6 diffrents diffrentssont envisager pour desmatriaux les respectivement Les sontprsents. Figures2-l4a et 2-l4b prsentent rhologique modles Ces cinmatiques. comportements et bilinairescinmatiques multilinaires comportements Finis. Is Figures2utilisspour les acierspar les logicielsdesElments ont t largement En isotropes. et isotropes multilinaires bilinaires I4c et 2-L4dmontrentles comportements dansles Figures2-14"et et anisotropes de Druckersontprsents dernier,les comportements dansla Figure 2isotropede type Druckerestprsent cinmatique 2-l4f.l comportement l5c. utilispour les solset les btons. de le Gnralement, comportement typeDruckerestsouvent est et cinmatique isotrope psentdanslaFigate2-16. du I-e concept comportement Finis, la rsolutiondu systme des Dans les calculsnonlinaires,par la mthode Elments dans la non-linaire est trs importante. Dans un problme non-linaire apparaissent physiques deuxtypes: de formulationdesproblmes de physiques indpendants {z} dansun modlelinaire, sontsupposs 1 . Lesparamtres des ... de telsque le modulede Young,lescoefficients conductivit etc.Ils deviennent non et le fonctionsde{z}dans casde la plasticit descoulements newtoniens. dans apparaissent les du par Destermesnon-linaires rapportaux inconnues problme quations aux drives partielles, mme lorsque les proprits physiques sont on dplacements, a: avecgrands en de indpendantes {r} par exemple lasticit

32

2.2.3

Chanitre 2. Mthode de prvision de la dure de vie en fatipue

u 1fav\' t'=E.zla)

(2-7r)

o t, , uet v sontla dformation dansla directionx, le dplacement dansla direction x et le dplacement dansla directiony respectivement. La mthodedes ElmentsFinis conduit une formulationdiscrtise problmesnondes qui linaires peuts'crire:

ffr=6v+N =talr(tr(r[t )-{r})=o

pour {a}' * o

(2-72)

L'quation(2-72)peut s'criresousla formesimplifie(2-73).

e)
Figure 2-14. Ditrrentscomportements s'appliquent qui gnralernent dansles calculspar Elmen Finis a) bilinaire cinmatique b) multilinaire cinmatique c) bilinaire isotrope d) multilinaire isotrope e) anisotrope f) Type Drucker.

33

2.2.3

Chapitre 2. Mthode de prvision de Ia dure de vie en fatigue

a)

2D

Ot= 02= 03

2D
b)

Ot= Oz= 03

c)

3D

2D

Figure 2-15 Prsentation comportementa) cinmatique b) anisotrope c) lype Drucker dans les cas bidu dimensionnels tri-dimensionnels. et

34

2.2.3

Chapitre 2. Mthode de prvision de Ia dure de vie en .fartsue

Surfacecoulementsubsquente

;e coulementsub subsquente

initiale Surfacecoulement

rface coulementinitiale rent

\ //
I I t

I I

a)
Figure 2-16 Comportement isotrope et cinmatique plastiques.

b) tel qu'il existe dans des matriaux nonlinaires lasto-

lx(u)I{"}={r} ou {n(,)}= [x(")l{"}=o {r}le dplacement, vecteurde la forcenodaleet le vecteurdersidu.

Dans certains cas en plasticit, la forme de l'quation (2-73) seraincrmentale:

(2-73)

la o tr(r)] , {u}, {r'} et {n(r)} sont respectivement matricede rigidit, le vecteur de

lr(,)h,)=br]
rena {n(z)} nul. que exacte proche possible zro.Lasolution de rsidu {n(r)} aussi itrative@gure 2-L7). La recherche la solution {r} sefait demanire de

(2-74)

non-linaire(2-73), il faut chercherun vecteur {z} qui rend le Pour rsoudrece systme

l. Changer I'algorithme 2. Changerla solution initiale

et Figure2-17 Algorithmeitratifpour lescasnon-linaires lasto-plastiques.

35

2.25

Chapite 2. Mthode de prvision de la dure de vie en fatioue

d'quationslinaires chaque conduit rsoudreun systme La majorit des algorithmes : doit tenir comptedeplusieursfacteurs de itration.Le choix d'un algorithme rsolution ou ou 1. I-e type de nonlinarit : localise non,prpondrante non solutions d'uneou plusieurs 2. L'existence d'une solutionapproche de d'unemthode construction 3. La disponibilit dsires et 4. La prcision la rapiditdesconvergences 5. Le risquede divergence de valablepourtous les cas; la stratgie il En pratique, n'existepasde mthodegnrale et donne en faisantappel de par doit s'adapter exprience une classe problmes rsolution : de des unecombinaison trois mthodes basesuivantes . Mthodede substitution

o Mthodede Newton-Raphson o Mthodeincrmentale

Mthode de Substitution
une consiste construire suitede solution.trt), {ut} ... {r'} calcules partir Cettemthode
de {z-t} en rsolvant le systmelinaire :
r .1

[,*(r-)k' ]={r} i =r,2,3...

le en Ce qui peuts'criresousformeincrmentale introduisant resiau{n'}:

(2-7s)

{r}{*(u"')}= [r(,'")][u"'] {n'}=

["(r'-'l[or-']= {n'}
{r'}= {r-'}*{or'}

(2-76)

les Le vecteur{u-t} etantconnu,il est possiblede construire matriceslmentaittt[ft(r"t)], pour obtenir [t(r"tl, puis de les assembler Lu'parI'unedesmthodes. est L'algorithmecorrespond^nt(2'76) le suivant: linaire(2-76) en le et enfin rsoudre systme

nulle. n"ntuellement l. Calculerunesolutionapproche{uo}, lmentairetff} . des 2. Construire{f} p- assemblage vecteurs

36

3 . i =1,2,...(pourchaque itration) 4. Pourchaque lment: a. Extraire valeurs dplacemenrs a" b-r} les des {":'I
b. Calculerla matricede rigidit lmentaire, [t(r"tl

c. Calculer rsidu le lmentaire {f }- VrX":,} t)= d. Assembler comme dans problme un linaire :

[r] oans [r] {'} a-s {n'}

5. Rsoudre commedansun problme linairer [f]{Ar'}= {n,} 6. construire nouvelle la esrimation la solurion {u'}= {u-t}* "'} de :

(2-77) (2-78)

7. calculer norme o" {or,} * ll.ll ae{n'} la llrll

8. Appliquer testdeconvergence un urilisant oullz{l llzll
Le facteurde sur-relaxattonautilis dansle stade6 permetsouvent d'amliorerla vitesse de convergence. valeuroptimaledpend problme; elle estdtermine I'exprience Sa du par en calcul numrique. Dansles problmes plasticitpar exemple, de elle se situeentre 1.7 et 1.9. La mthode sans sur-relaxation correspond a=1. Pourle testde convergence I'itrationi, il estncessaire d'utiliser soit la normemrimale :

uaxltu,l' llmll= ou uy,tn,l' llnll=

soit la normedesmoindres cars :

(2-7e)

y"11=,[ffi]o, p1=ffin1
En pratique,desnonnesrelatives sontfrquemment utilises :
| . tj

(2-80)

llnll=M@pl " " t l"tl

s'crit:

(2-8r)

Si z, est trs petit par rapport la valeurmoyennedes termes de{u} ou gale zro,u,

11,;1=.7]'{-T ''

,lb'Y {",}

(2-82)

Le processus itratif estarrtlorsque:

37

2.2.3

Clapitre 2. Mthode de prvision de la dure de vie en .fatisue

pour la normale(2-82): e = 0.05 Parexemple, graphiquement Dansun problme unevariable, I'algorithmede substitutionestprsent sur la Figure(2-18)de manire souligner possibilits divergence. les de K( " ') x ( "' )

ll,ll. "

(2-83)

a)
K(u)-u = P R o= ,K(r') ,((r')

rc("')

b)
Figure 2-18 Algorithme itratif pour un comportement non-linairelasto-plastique, mthodede substitution a) b) Convergence Divergence.

38

2.2.3

Clnpitre 2. Mthode de prvision de Ia dure de vie en fatisue

Reformulons l'quation (2-76) en dcomposant[f] constante d'une matrice [r<,,] fonctionde z : et

en la somme d'une maffice K]

(t",* [r, 0- I[o", {n, I ]= }

lx,], lx*], br')

(2-84)

par les sont respectivement matricesconstruites assemblage des "t {n'} matrices lmentaireslinaires [t,]. Cette maffice est fonction de u le dplacement, pour les dplacements du sidu pour l'valuationde I'exactitudedescalculs. incrment et Dans I'algorithme de l'quation (2-76). il faut assembler dcomposer[f] et a chaque

itration,ce qui est trs coteux.En ngligeafi lK,,t] Oansl'quation(2-84),1'quaton suivante s'crit:

[r,]{nu'}= {n'}
{"'}={u'' }* {or'}

(2-8s)
(2-86)

La matrice [f,] peut tre dcompose seulefois ; A chaqueitrationil suffit de calculer une {n'} puisd'valuer{ou'} partiro" [r,] djdcompose.

Mthode de Newton-Raphson modifie

l. Calculerunesolutionapproch"{ro}, nulle. ventuellement 2. Construire{f} p* assemblage vecteurs des lmentairerff}. 3. Construire[ ] p* assemblage matrices des lmentaires linairerlfr,] . 4. Dcomposer[K,]. 5. i =I,2,... (pourchaque itration) 6. Extraire{ri'} a" {u''}no* chaque lment. 7. Calculerle rsidu {^'} p- assemblage rsidus des lmentaires :

{'}={f }- [ril,-]
8. Rsoudre o. [&ilar']= {*'}u partir [r,] dcompose.

e-87)

g. Calculer =I ) fu'I={r-'}*fur'} t"nsupposantar

10. calculerllnll. 11. deconvergence Test utilisantllzll .

39

2.23

Chapitre 2. Mthode de prvision de Ia dure de vie en fafioue

dansle cas d'une variablesur la Figure 2-I9. graphiquement Cet algorithmeest prsent modifi,la matric" [K,] est assemble de Dansl'algorithmede la mthode Newton-Raphson et il une et dcompose seulefois, alorsquedansI'algorithmede substitution, faut assembler modifi est souvent dcomposet [K] chaqueitration.L'algorithme de Newton-Raphson faiblement non-linaire.Par contre, pour les employ dans le cas d'un comportement est de problmes fortementnonlinaires,la mthode Newton-Raphson plus souventutilise quela mthode substitution. de
x, = x,(uo)

K(u)-u= P

ltt

tl2

u?

u*t

modifi. de la Figure2-19 Algorithmeitratifdansla nonlinairelasto-plastique, mthode Newton-Raphson

pour les problmesfortement non-linairesconvergeen La mthodeNewton-Raphson gnr plus rapidement. al

La mthodede Nervton-Raphson
u'-t de la solution qu' I'itrationi-1, nousavonsobtenuune approximation Supposons telle quele sidu ne soit pasnul.

o {r}k(u"'))= [r(u-')[r'"]*
: ci-dessous ui A I'itrationi, uneapproximation de la solutionpeuttrepropose

(2-88)

e {*(r')}={^(u-'+ az')}=

(2-8e)

de ce L'algorithmeestobtenuen dveloppant rsiduen sriede Taylor au voisinage ui-t :

40

2-2'3

Clnpitre 2. Mthode de prvision de la dure de vie en fatieue

* {o(,'2,)}= )}.[u#],*,' ... 0 t^0{-,}* =

D'o, en ngligeant termes les d'ordre suprieur 1 :

(2-e0)

-[#].-,,... )} '' = ]* {^(,"'

(2-er)

ou:

,]= [x,p-,1f,t {*(,"')}

{"'}={r'" }* {or'}

(2-e2)

L'expression de la matrice tangente lX,@-t)f s'obtienten drivantl'expression (2-73)du rsidu :

-[#] -[#] rr,(,)t== .rol.[# e]]

Dans le cas o F est indpendantde u :

\z-gz)

(r)]= )l [r<, tr(,*|W rrl

L'expression (2-94)s'crit sousla formeindicielle, ci-dessous :

(2-e4)

(r,),, t<,,+1.*u, = --u \--t t,j

| U,

(2-es)

L'algorithme correspondant l'quatron(2-92) ressemble I'algorithme de la mthodede substitution, cependant estremplac parftl. [f] 2-2O, dans casd'uneseule le variable. Il estprsent graphiquement la Figure sur

Figure 2-20 Ngorithme itratif dansla nonJinairelasto-plastique, mthodede Newton-Raphson. la

4l

globale [f,] La matricetangente tangentes ]. lmentaires ft

"r,

des par en obtenue pratique assemblage matrices

Mthode incrmentale (Pas Pas) Selon prcdentes. le itratives dans initialejoueun rle important lesmthodes La solution non vers ou peuvent diverger converger unesolution les choixde cettesolution, mthodes
acceptable. La mthode incrmentaleconsiste remplacerla rsoluton de

t,{ro}={r} [r(r)J[r]=
de parla rsolutionsuccessive

(2-e6)

= ["0,)k,]),,{F,}
o ["(rr)],{rr} {r}= 4{r"}

(2-e7)
sont la matrice de rigidit au j eme pas, le vecteur

", s'exprime Dansla relation(2-97), .,1., global et le vecteurdesforcesappliques. dplacement : ci-dessous corrlme

tr,=Ar,4,...,trn

(.2-98)

La solutioninitiale utilisepour calculer u,est la solution zr-,obtenue l'tape prcdente. Chaquetape constitue un problme non-linairequi se rsout avec une ou plusieurs modifie. ou itrationsde la mthode Newton-Raphson de la mthodede Newton-Raphson utilise une itrationde Newton-Raphson chaquetapeet s'crit La mthodeincrmentale pourun niveaudonn de sollicitations2, :

K,-. K,(u')

ou Figure 2-21 Algorithme itratif, mthodeincrmentale pas pasavecune itrationde Newbn-Raphson.

42

2.2'3

Clwpitre 2. Mthode de prvision de Ia dure de vie en fatigue

-[.(r,-, )J[r,-, {o(r,-, tr,-,{Fo} )}= ]

{",}={r,-,}*{our}

[",(r,-, "}= {n(, @,1,-,\r,} Xo ,-,)I* -

(z-ee)

Cet algorithme prsent est graphiquement unevariablesur la Figtne2-21. pour La mthodeincrmentale utilisant plusieursitrationsde Newton-Raphson s'crit pourun niveaudonnde sollicitations,1, :

+ -,a,,_,\ro} [*,("':l[o"; {*(,;")} (s", ]=

[,',]={,;'}. {o,;} i = 2,3,...
pouri = 1, on utilisedirectement (2-99).

(2-100)

Toutesles mthodes prcdentes peuvent ramener un seulalgorithme se qui, un niveaude sollicitationdonn,estschmatis ainsi :

Choix de la solution approche Correction: Amliorationde la solution .Calculdu r*idu {n'} .Rsolution ae [r" [lu']=

{n' }

. {r*'}={u'}+at$u'}

Test de oonvergenoe

Figure 2-22 Algorithmegnralitratif dansle casd'un comportement non-linairelasto-plastique.

du Quelle que soit la mthodeutilise,I'expression rsidu {R} restela mmecar elle est caractristique l'quation rsoudre. de Par contre I'expression la matrice [f"] varie de d'une mthode I'autreet influencela vitesse convergence de : [f.]= [f] [f,]= [f,] pourla mthode substitution de po* la mthode Newron-Raphson de modifie

43

[rK.]= [f,]

de po"t la mthode Newton-Raphson

itratives les ae[r"], nouspouvonsretrouver mthodes d'autresexpressions En choisissant, d'quationslinaires.I-e facteurde sur-relaxationa les utilisespour rsoudre systmes Le permet souventd'amliorer la vitessede convergence. facteur de sur-relaxationest trs efficace dans les cas de grands problmes qui sont rsoudre gnralement interne la mthode Il frquemment. peut tretrouvpar essaiet effeur avecun algorithme Finis. desElments dpenddu problmetrait.Il est parfoisncessaire I_echoix d'une desoptionsprcdentes ou d,utiliserdespasde sollicitationtrspetitspour desraisonsde convergence parceque la commeen plasticit.Il peut alors suffire formationdu problmeest de natureincrmentale pas rduit d'itrationprr p*. Les mthodes pasavecutilisationdesmthodes d,un nombre dans le calcul en modifie sont souventappliques et Newton-Raphson Newton-Raphson Finis. par plasticit les logicielsd'Elments

2.3 Conceptde facteur de concentrationdescontraintes

La distributiondescontraintes voisinaged'une entaillejoue un rle trs important.Il est au notamment ncessaire d'valuer le maximum de cette distributionen fond d'entaille. Le facteurde concentration contraintes propos des est pourvaluer maximum. ce Dan le domainelastique Peterson dfinit le facteurde concentration contraintes [4] de de deuxfaons: 1. Le rapportde la contrainte maximale fond d'entailleet la contrainte nominale. 2. Le rapportde la contrainte maximale fond d'entailleet la contrainte globale. La formegnrale facteurde concentration contraintepeuts'crire: du de , ^, - o-* *
o,

(z_t0t)

o k,, o** et os sont respectivementle facteur de concentrationde contraintes,la contrainte lastiquemaximale fond d'entaille et la contrainte globale appliqueloin du fond d,entaille. Il reprsentele facteur d'amplification de la contrainte globale applique. Dans le cas d,un petit trou dans une plaque, la concentration de contrainte dpend aussi du chargement appliqu (Figure 2-23).

o;

o:

o:

o;
Figure2-23 Distributiondescontraintes dansle casd'un hou dansuneplaquey comprisle chargement gnral.

45

23

Chapitre 2. Mthode de prvision de Ia dure de vie en fatisue

Dans la Figare 2-23, o*,,

o!, oi, k, et f(o!,oi)

la sont respectivement contrainte

maximale,la contrainteglobale horizontale,la contrainteglobale verticale, le facteur de et de concentration contraintes la fonctiondu chargement. pour les de les I Tableau2-2 rassemble valeursdu facteurde concentration contraintes cas diffrents de chargement:
pow un petit trou dansuneplaque. de et Tableau2-2Ditr&entschargements facteurde concentration contraintes TYPEDU CHARGEMENT

o!
oo oo oo

oi
0
oo -oo

r("X,oI)
oo oo oo

k, 3

o*,=k,f\o!,oi)
o^ =3oo
or* =2oo o** -- 4Jo

2 3

2
4

I-es valeurs du facteur de concentration des contraintes sont obtenues par les calculs analyiques. Par exemple, le calcul analytique pour le type de chargement I du Tableau 2-2 a t dj donn par les quations (2-10) (2-L2).Is quationsmentionnessont crites ci-

: dessous

= -5).(,. ",,,] il(' +-+), i

= 5)-.#)"o' zef +l('. [' 3
t,, --1[ r -5": *41sin2o
oo 2[ r* r' )
- / ^ a. ^ z\

(2-ro2) (2-ro3)
o

(2-r04)

La contrainte circonfrentielleest la plus levee et le rapport entre les contraintes et circonfrentielles globales sera utilis pour dfinir le facteur de concentrationde peuts'crire: maximale circonfrentielle La contraintes. contrainte

(2-l0s)
appliqu dpenddonc du type du chargement de Le facteur de concentration contraintes (Tableau2-2).

46

Les facteursde concentration contraintes de peuventt obtenuspar descalculsanalyiques, destravauxexprimentaux calculpar Elments ou Finis.

. Exemple
Nous avons pris la gomtrie suivante pour calculer le champ de contraintes par le code du calcul Elments Finis. Dans la Figure 2-24, le calcul par Elments Finis et la sclution analytique sont comparspour le cas d'un petit trou dansune grande plaque. Is conditions du calcul par Elments Finis pour le problme mentionn sont donnescidessous: o o o o . o o o o Module de young Coefficient de Poisson

E=270 GPa
v =0.30

Hypothsed'un comportementmcaniquelastiquelinaire Contrainte globale o! =t MPa et oJ = 0 MPa

Hypothsede contraintesplanes Symtrie double pour le modle ElmentsFinis Longueurde plaque Largeurde plaque Rayon du trou L= 0.800 m

W=0.600 m e0.050 m

Epaisseur la plaque 8=0.001m de

a)

47

I@

m.l @.1 -l
ffi.0 m '.ml q r-,qt{la G.r.lq

.st&

.9655

l.lta

:.{t

!.14

b)

c)

Figure 2-24 Distribution des contraintes dans le cas d'un Eou dans une grande plaque pour un chargement 6i = ao
h,^

, oi =U

a) Maillages et conditions aux limites

b) Composantesde la conFainte dans la direction

x c) schmades contraintes au voisinage du trou.

Dans la Figure 2-24,|a distribution de contrainte oo est prsente. L'utilisation de contrainte globale o! =t MPa permet de reprsenterla distribution des facteurs de concentration de contraintes. La comparaisonentre le calcul par Elments Finis et la solution analyique est montre dans laFigure2-25.
3.50 3.00 2.50
3oo

-+ -# + *-

o0/o0 SolutionClassique orlo0 SolutionClassique o0/o0 Solutiond'eles finis orlo0 Solutiond'elerems nis

o 2.fi) b 1.50 1.00 0.50

30o

0.m
0.00 1.00

1.50 2.w fla

2.50

3.00

Figure2-25 Comparaison enhela solutionanalytique le calculpar Elments et Finis pour le problme d'un petit trou dansunegrande plaqueavecun chargement horizontal traction[28]. de

48

2.3

Chapitre 2. Mthode de orvision de Ia dure de vie en .farteue

o o, ,oeetoo sont respectivementla contrainte radiale, la contrainte circonfrentielle et la contrainte globale.

2,50

+ -0-

des otd Solution Elnpnsfinis finis des orlo0Solution Elnrcnts

0,00 0,50 a)
6,00 5,00 4,00 3,00
b

r,00 r,sO

250 3,00 3,50 4,00

+o0/o0 .orlo0

Solution Elnrens des finis Elrrensfinis Solutiondes

2, 1,00 0,00 -1,00

0,00

0,s0 1,00 l,so

?E

z,so 3,00 3,50 4,00

b)
radiales et dansle casd'un trou dansune plaquepour les composantes Figure2-26 Distributiondescontraintes circonfrentielles dans une plaque y compris le chargement a) o! = Oo O", = oo b)

o ! = o o , o ',= - o o .

des aussidela gomtrie prouvettes. dpend de Le facteurde concentration contraintes

49

Dans la Figure 2-27, une bane en traction est considreet un huitime de la gomtrie de la bare est prsent.

a)
Ftrt soilfilot

,mPrI tB .I lDE.l /,D^!DED rt lv6l 3tt5.0 DID('.I?-03 $.-.0112 S(.1'6e9

.7t797 .353t63 'ttg66:l .rtlo u barE aa rlloa ulrl L clrEga! Dou

-,0tl2a

.9@irr5

b)
rtD^L l0t0Tr0r stlD.l
gE flE'l .l

L " ('Tr tgllt.0

DD(.,r Is r-,01 tl!( '1,39

,/E(I,DID

31.-.0r1? {3

LG qhcCt

.l?106 utrl Dost &

.93566? bqE aD llGloD

c)
dansle cas d'une barreen action D=0.04 m, d=0.02 m, 10'01 m, Figure 2-2? Distributiondes contraintes Oo =l MPa de dansla rgioncritique' c) y a) Schmab) Composantesdescontraintes Gradient contrainte

50

2.3

Chaoitre2. Mthode prvisionde la durede vie en Tatieue de

Le facteur de concentration des contraintes change avec la variation des rapports D/d et r ld,. L'volution du facteur de concentration de contraintes en fonction de la gomtrie est

gnralement prsente dansdesformulaires. Dansle Tableau2-3,les typesde changement sectionsontprsents. changements de Ces de gomtrie provoquent modifications facteurde concentration contrainte. des du de
Tableau 2-3 Types de changement des sections et chargement dans les cas de contraintes planes et axisymtriques [28].

TYPE DU CHANGEMENT DES SECTIONS

CHARGEMENT

Contrainte nominale

Eshelle

Cas

Jht "
f

Traction
D
P

P 2bt

h
ir;

Flexion

m
P bt

3M

Traction
D b
M P

Flexion B

b't

6M

Tension

P 2bt

3Mh Flexion
iti

2t(c3 -h3)

51

2.3

Chanitre 2. Mthode de prvision de la dure de vie en fatipue

Tension
h

Flexion
d=2b
P T

4M
-

h
h

Cisaillement

r.23V
lb'

Torsion

2M,
7b3

+ ln:,Jbta-s.8 A+ ->c lr-Jht"-4.7 =2.4+c ff - Jht"

+ k, -J.{
/

*c7 -f3zl

1214567

JilE'

'* J5

Abaquede Neuberpour calculerle facteurde concentratron conrarns. de

I-escasprsents dansle Tableau sontutilissdansle diagramme Neuberpour calculer 2-3 de les facteurs concentrations contraintes de de de t28]. Il est ncessaire calculerles paramtres paramtre o dansle Tableau2-3.1-e, Jbt. et la valeurdansla 5mecolonne ""tJW (la valeur de <<n> dpenddu type de chargement prsent dans le Tableau2-3) donnentle JW point < A >. Une horizontale passant par.(A> coupeI'axe verticalen <<B>. Point << > et le Le B paramtrecalcul JW o donnent le facteur de concentrationdes contraintesavec

I'utilisation la 4mecolonne. exemple, de pour n=3 et Jbtt =5.8 avec Jhla =4.7 Par pour <<e et Jht"=2.4 > gale 3.4. pout<<f la valeurdu facteur concentration contraintes >>, de de est

2.4 Facteur de rduction en fatigue

La limite de rsistance fatigue oD,Npour les prouvettes entailles moindreque celle est en d'une prouvettelisse oo,r. I-e rapport de la limite de rsistanceen fatigue pour les prouvettes lisseset entailles avecla mmesectioncritiquedfinit le facteurde rductionen fatiguek,

k, =%'

(2-to6)

ao-u

les en o or., et oD,N sont respectivement rsistances fatigue d'une prouvettelisse et entaille. I facteurde rductionen fatigueest moindreque le facteurde concentration contraintes de (I3kr <k,). La diffrence plus grande lastique entre k, et k,sera dansle casdesentailles aigues.La sensibilit I'entaille est gnralement mesurepar I'indice de sensibilit I'entaille(q).n est dfini commela valeur du rapportdu facteurde rduction fatigue en

k, diminud'une unit sur le facteurthoriquede concentration contraintes,k,galement de

diminu d'une unit :
e =J'

k, -l

k,-l

(2-r07)

La connaissance cet indice q qr:jmcaniques matriau de dpenddescaractristiques du et en particulier de la charge la rupture,doit permettreen principe d'estimer ft, ; ce qui conduit la connaissance la limite d'endurance de sousentaille sansavoir raliserdes essais.La dfinition dans l'quation (2-106) est prcisment utilise dans le cas d'une contrainte normale(tractionet flexion). Commele facteurde concentration contraintes, de la dfinition gnrale crit: s' k, =%, " oo,* k', = b , k,=o^o o, k,=*

(2-108)

To,r

t,

(2-r0e)

o o,,,"* et os To,s,o,u, T^et

sont la contraintenormale maximaleet la contrainteglobale normale. les tssont respectivement limites d'endurance d'prouvette lisse et

entaille dansle casdu cisaillement, cisaillement maximalet du cisaillement global. du

53

2.4

Clnpitre 2. Mthode de prvisionde la dure de vie en fatieue

Une autre dfinition du facteur de rductionen fatigue est le rapport de I'amplitude de contrainte d'uneprouvette lisseet cellede la contrainte d'uneprouvette entaillee. os(Nn) ft- = ' o*(N) (2-llo)

k,=* '

tr(N *)

e-nD

o o, (N* ) , d" (N^ ) sont respectivementles amplitudes de conffainte d'une prouvette lisse et entaille conduisant au mme nombre de cycles rupture. zr(N*) et r"(N^) sont les

amplitudes de scission d'une prouvette lisse et entaille conduisant au mme nombre de cycles rupture . Le facteur de rduction en fatigue est un paramtreexprimental obtenu partir d'essais de fatigue. Ce facteur prend une valeur comprise entre I'unit et le facteur de concentration de contrainteslastique( ; il dpenddu nombre de cycles rupture. l<kr <k, Des relations ont t proposespar les diffrents auteurspour relier k, et k, ; on cite: Q_lL2)

Peterson[29] Q=-,
'l

I
qP ,

ar=ar(ou)

(2-r13)

lT-

a, estun coefficient qui dpendde la contrainte ultime, et p le rayon d'entaille.

ap=l -

(zto\'''
\ou )

e-rr4)

ou est la contrainte ultime en MPa et dp en mm.

Neuber [30]
Q= -, I QN= a*(ou) 12-LIS)

l+-f'3tlp

div un coefficientqui dpend la contrainte de ultime. o Heywood[3]-32]

54

2.4

Clnpitre 2. Mtlnde de orevisionde Ia darede vie en latieuc

+=+-,

k'

aH=or(or,4,,Ar)

(2-r16)

lr.rli)

ultime, de un o a" , 4et A, sontrespectivement coefficientqui dpend la contrainte la sectionnominaleet la sectionglobale. stieler-siebel[33]

kt=L. k, t+,tSo Z'

Sn= S(ao.r)

(2-Lt7)

un o Sn et Zsont respectivement coefficientqui dpendde la limite d'lasticit llel=yf*# au O.2Vo gradientdescontraintes fond d'entai etle Kuhn-Hardraht [34]
, eK = a*(ou)

(2-rr8)

o a* , p, a et ou sont respectivementun coefficient qui dpend de la contrainte ultime,le rayon d'entaille, I'angle d'ouverture de I'entaille et la contrainte ultime.

BuchSwitech t35-361

l=+l'-#l'
simultanment. Topper-ElHaddad[37]

po=po(A'h)

(2-t19)

du o A et h sont descoefficientsqui dpendent matriauet du type de l'prouvette I-e respectivement. paramtre Po est une constantequi dpend de A et de h

(2-r20)
o ao, Fo et 1osont respectivementla profondeur de I'entaille, le facteur de gomtrie et une constantequi dpenddu matriau.

Selon Topper-ElHaddad,le facteurde sensibilit relationsuivante:

peut tre dfini par la

2.4

Chapitre 2. Mtlnde de nrvision de Ia dure de vie en fatipue

p It 4Io n=Il"

Pour Pour

p 34Io p>4Io

(2-r20)

Les quationsproposes Topper-El Haddadtrouvent leur fondementdans une par approche type mcanique de linairede la rupture. Le Tableau24 prsente rsumdesdiffrentes un proposes cesauteurs. quations par Ces relations donnent estimation facteurde rduction fatiguek, . une du en

Tableau Rsum relations pour I'estimation facteurde rduction fatigue ft, 2-4 proposes des du en PARAMETRE DU MATERIAU
o=-1

Person [29]

'

k, -l

=-

k,-l

t+b p
-I

a, = ar(ou)

^ k,-l
kr-

k,-l

,,* lA 1P
1
dn =aH(ru,4,er)

Heywoodl3l-321

r-4
kr1

tr-W
Kuhn-Hardraht [34]

So = S(40., )

l*'^E p ft-(t)l
L=1[,- z'rh f
k, AL P+po)

a* = a*(ou)

Po = Po(A,h)

Io = Io(o,ou)

I-es quations proposes Peterson, par Neuber,Heywood,Stieler et Topper comportent un seul paramtre<< matriau>. Ce paramtrematriaua t suppos constantpour tout type d'prouvettes de chargement. et

2.4

Chapitre 2. Mthode de orvision de Ia dure de vie en fafisue

de le Selon cette hypothse facteur de rductionen fatgue et le facteur de concentration Cette hypothsen'est pas toujours contraintessont gaux pour une entaille mousse. Le par L'quationpropose Buch- Switechutilise deux paramtres. rapportentrele correcte. une facteur de rduction en fatigue et le facteur de concentrationde contraintesreprsente autredfinitiondu facteurde sensibilit I'entailleq . Cerapports'crit commesuit :

f,= f @,'+,n)

(2-r2r)

des le o p, A et h sont respectivement rayon d'entaille et les constantes matriaux.I-es mais pas constantes et h ne dpendent du rayon d'entaille ni de la taille des prouvettes, A la Dans la mthodede Buch-Switech, du dpendent type du matriauet du chargement. en des contrainteo, = Aoo correspond I'amorage fissures fatigue.Dansla Figure2-28,un schma la mthode donn; de est

oi

=6,

h et des entaille dfinition paramtres A. avec prouvette Figure 2-28Schma d'une

la Dans la Figure 2-28,W, p , o*, , on,an7 et od sont respectivement largeurde la plaque, le rayon d'entaille, la contrainte maximale, la contrainte nominale, I'endurance en fatgue pour une prouvetteentaille et I'enduranceen fatigue pour une prouvettelisse .

57

2.4

Clapite

2. Mthode de prvision de la dure d,evie en fartpae

On considre que A et h sont des constanteset oy la limite d'lasticit. I gradient des contraintesdansce cas peut tre dfini par l'quation suivante : do __or** - o, dy h Q-122)

Selonla distribution thorique contraintes, gradient contraintes des le des s'crit: do _Co** dvp

(2-r23)

o, p , a**et C sont respectivement rayon de I'entaille, la contraintemrimaleet le le coefficient de proportionnalitobtenu par la solution lastiquedans diffrents cas de chargement.
La considration quations(2-122)et (2-t23) conduit l'quation (2-124) : des do _ oro* - oy =9o*-

dy

;=

o^=m

oy

(2-r24)

l'-v)
(2-rzs)

La substitution o, = Aoodans (2-124) de l'quation donne relation la suivante : Ao, 6r**=ffi

l'-v)
L'quation(2-L25)peuttredivisepar la contrainte nominale En raisonde l'galit de la o". contrainte nominaleon ave& limite d'endurance fatiguepour les prouvettes la en entailles (2-126) critesous forme: o*,la relation est la o** _

on

Aoo

(2-126)

"^(, i)

Selonles dfinitionsdu facteurde concentration contraintes du facteurde rductionen de et fatigue,la relation(2-126)estrcritesousla forme suivante :

(2-127)

OU'

9=1h- "ol=(p,A,h) k, A[ P)

(2-r28)

2.4

Chnpitre 2. Mthode de prvision de la dure d.evie en fatipue

(2-123), inffoduitune on des infiniespourle gradient contraintes Afin d'viterdesvaleurs p0 et petite d'entaille dutypeduchargement. du longueur qui dpend rayon kr ---r r ( k, A(

ch \

p+ po)

fz-Lzg)

a,rial et d'une entaille au centred'une plaque Dans le cas de I'utilisation d'un chargement (2-L29)s'critdenouveau: (Figure2-28),1'quation

kr r (. z.rh \ -'-.=- 6.3h * =-l 1- p+ po l, po=po(h,A)(3- A)

k, A[ t

(2-130)

h h La relation proposepar Buch-Switechcontient deux paramtres et A. Le paramtre A au dpenddu gradientdes contraintes voisinagede I'entaille. Le paramtre dpenddu et cyclique. Dansl'quation(2-l3O),le gradientdescontraintes le coefficientd'crouissage simultanment. sontconsidrs coefficientdu matriau I-e gradientdescontraintesZ est la valeurde la pentede la tangentede la distributiondes maximale[38]. contraintes fond d'entaillerapporte la valeurde la contrainte en
T = llIfl-

I do(x\ .. r{ o'"* dx

(2-13r)

7 s'expimeenmm-I. pour et Dansla Figure 2-29,1'volution champdescontraintes du gradientdescontraintes du I-e de sont prsents. une contraintenominale on et le facteurde concentration contraintes Le dansles casde tractionet de flexion estaussiprsent. gradientdes champdescontraintes montren fond d'entaille. contraintes exclusivement est

I I I I "---.---.$'"I I I I

a)

59

2.4

Chapitre2. Mthodedeprvision de la durede vie en futieue

b)

c)

o*,
Figwe 2-29 Champdes contraintaset gradientdes containesau fond d'entaille a) Champdes contraintesen haction b) Gradientdescontraintesen traction c) Champdescontraintes flexion d) Gradientdescontaintes en en flexion.

60

2.4

Chnnitre 2. Mthode de nrvision de Ia dure de vie en fatipue

Il estncessaire calculerle gradientdescontraintes voisinage entaillesen sebasant au des de prsente sur la procdure dansla Figve2-29. Dansla Figure2-28,d, D, F, M, on et o** sontrcspectivement petit diamtre I'arbre,le le de grand diamtre de I'arbre, la force, le moment, la contraintenominale et la contrainte maximale. par Les contraintes : maximales nominales et s'expriment lesrelationssuivantes

or**=kr6n=fr
4
O- =-

(2-r32)

'

!a,
4

(2-r33)

Le facteurde concentration contraintes dansl'quation(2-132) s'cit dansle cas de la de traction@gure 2-29a et 2-29b):

ke

l'ffi'

ko=

(2-134)

1l;

o a et p sontrespectivement profondeur la d'entailleet le rayond'entailledansla Figare228. pour les L'utilisation du conceptde gradientdes contraintes se fait pas exclusivement ne calculslastiques. le serautilis pour Dans le cas lasto-plastique, gradientdes contraintes trouverle volumed'laboration processus fatigue. de du

2.5 Les mthodes d'estimationde la dure de vie en fatigue

La rsistance fatiguedesmatriaux traditionnellement en prsente le diagramme est par de Whler. Ce diagramme s'appellegnralement diagramme S-N. I diagramme V/hler, de est dcritpar deux paramtres. premierest I'amplitudedescontraintes Le ( appliques o) et le deuxime le nombrede cycles rupture (N^). I diagramme Whler est prsent est de sousforme semiJogarithmique bi-logarithmique. relationde BasquindonnantI'allure ou La de ce diagramme donne : par est ao = o'r(2N*)' o o, , d, (2-135)

N* et sont respectivement I'amplitude de la contrainte,le coefficient de

Basquin, nombrede cycleset I'exposant la relationde Basquin. le de Dansle diagramme Whler,trois zones de sonttraditionnellement distingues:
1000

1,0E{-03 l,0E+04

1,0E+05

1,0E+06

l,0E+07

1,0E+08

NR
Figure2-30Schma diagramme Whleret la relationde Basquin. du de

La Figure2-30 montrele schma d'un diagramme Whlerclassique. de Les zones prsentes dansla Figure2-30 sonthabituellement dtermines le diagramrne pour lOa et 106 dans S-N cycles. En raisonde l'tendue nombrede cyclestudi,le diagramme Whlerestprsent du de sous la formelogarithmique.

62

2.5

de 2. de Chapitre Mthode prvision la durede vieen fatieue

La prvisionde la durede vie en fatigue est faite partir du diagrammede Whler. I-e I-es courbesde Whler facteurde rductionen fatigue est dfini partir de ce diagramme. d'prouvettes de lisseset entailles sontncessaires le problme la prvisionde la dure dans de vie en fatigue.Dans la Figure 2-3I, un schmasimpledes rsultatsen fatigue pour des prouvettes entailles lissesestpsent et
1000

a
I

l.t
R

bs

fs\

+Courbe Lisse * Courbentaille

^-l-

a-r

d
a

,(C
J(.

br

100 l,0E{{4

1,0Er-05

NR

1'oE{6

l,0E{7

pourdes Figure 2-31Schma diagramme Whler du de prouvettes et entailles. lisses Dans la ligwe 2-31, o*, oa,,, o'fi,o'r, bs et bn sont respectivementI'amplitude des

contraintes appliques dans le cas d'une prouvette lisse, I'amplitude des contraintes appliques dans le cas d'une prouvette entaille, les coefficients de Basquin dans les cas d'une prouvettelisse et entaille et les exposantsde la relation de Basquin dans les cas d'une prouvette lisse et entaille. Le facteur de rduction en fatigue se dfinit par:
nf

, , ---o ^
'oM

(2-136)

La mthoded'valuationde la durede vie en fatigue partir du diagramme Whler est de gnralement utilise dans les codes de construction.Par exemple, I'Eurocode3 utilise pour estimerla durede vie enfatigue[39]. l'tenduedescontraintes

63

2.5.1La mthodede ltBurocode3pour lesjoints souds

Le chapitre9 de I'Eurocode3prsenteune mthodegnralede calcul de fatigue des Cetteprocdure de en soudes acier soumis desvariationsreptes contraintes. structures n'est valablequepour les aciersde construction. prsente le chapitre3 de I'Eurocode3 dans sontutiliss: suivants les par propose I'Eurocode3, paramtres Dansla mthode des 1. L'tendue contraintes. moyennes. 2. I-escontraintes rsiduelles. 3. Les contraintes nominales. 4. I-escontraintes locales. des 5. I-esconcentrations contraintes 6. I-eseffetsde tailles. pour la prvisionde la fatigue paramtres sontpaspris en comptepar I'Eurocode3 ne Certains : desjoints souds (f 1. La frquence S4OHz). de 2. Les dfauts soudage. des 3. Le comportement matriaux. chimiquede I'acier. 4. La composition mtallurgique. 5. L'homognit de 6. Le procs soudage. soudage. aprs thermique 7. Le traitement apssoudage. de 8. Le traitement la surface 9. L'humiditet I'irradiation. par et prenden compteles effetsde la gomtrie de la fabrication de La mthode I'Eurocode3 de I'utilisationdesfacteurs scurit. de le Gnralement, phnomne fatigue tire son origine de I'applicationd'un chargement soit constante soit considre amplitude est variable. L'histoire descontraintes normalement avec amplitudeconstanteou variable est amplitudevariable.L'histoire des contraintes la dansla Figure2-32,o o^,, cf,;o, om, oa et Ao sontrespectivement conainte montre

&

2.5.1

Clnoitre 2. Mthode de orvision de Ia dure de vie en fatisue

maximale,la contrainteminimale, la contraintemoyenne,I'amplitude de contrainteet la diffrence entrela contrainte maximaleet la contrainte minimale.

no g 1oo
0
U) at

(/)

-100

-200

a)
500 400 300 200

Tinp

(t)

g 100
0 -100

.t) (rt

-m
-300

b)

Tin

Figure2-32Histoire contraintes Amplitude des a) constante Amplitude b) variable.

Dans le cas d'un chargement cyclique avec amplitude constante,la contrainte moyenne est I'amplitude quivalente. Mais dans le cas de I'histoire des contraintes avec amplitude variable, il est ncessairede trouver une dfinition partir de l'histoire des contraintes.

65

2.5.1

de 2. Chanitre Mthode orvisionde Ia durede vie en .fatigue

qui ressemble la mthode rservoir >> L'Eurocode3 utilise gnralementla mthode du << << Rainflow" [40]. Ces deux mthodesne donnent pas exactementle mme rsultat, mais leur dfinition de I'endommagementen fatigue est la mme. donnent cependantle mme rsultat en Rservoir >> Rainflow > et de << Is mthodes de << fatigue dansle cas d'une histoire des contraintestrs longues. La mthode la plus cofitmune pour prsenter une histoire des contraintes variables est d'utiliser I'histogramme des contraintes. Dans cet histogramme, les contraintes sont squentiellement slectionnes @gure 2-33). Ce concept est l'origine de la mthodePalmgren-Miner l4l-421que I'Eurocode3 utilise dans le cas des amplitudes de contraintesvariables. L'histogramme des contraintes classespar entendue peut tre approchepar la distribution de Weibull [43] deux paramtres. le Dans la Figure 2-33, n,, Lo,, no et ft sont respectivement nombre de cycles rupture pour une tendue de contrainte\,o,, l'tendue de contrainte, le premier paramtre de la distribution de \Meibull et le deuximeparamtrede la distribution de Weilbull.

La distribution de Weibull s'exprime:

-:

L,o ( , - - l Los@)\o Log(nr) oo I )

(2-137)

L,o

a)

66

2.5.1 Lo aoo

Chapite 2. Mthode de prvision de Ia dure de vie en fatipue

Lo
a6o

Log(n,)

LoS(n)

b) Figure2-33 a) Distribution domaine contraintesb) Varit contraintes ordre. du par des des

o Aoo , no, k sont respectivementI'amplitude maximale de contrainte et les coefficients de la mthodeWeibull. L'Eurocode3 est une mthode essentiellement base sur des contraintes. Il utilise

gnralementle diagramme Whler comme la recommandationECCS @uropeanConvention for Construction Steelwork Fatigue Recommandaton). La courbe de rfrence normale de la mthode d'Eurocode3 est montre dans laFigne 2-34. Cette courbe peut tre considre comme constitue d'une ligne ou deux lignes dans un di agrammebi-logarithmique. LogA,o
La courbe de rfrencede la rsistance fatigue en

Lo, La limite de la rsistance fatigue en

Loo

a)

Nc=2.08+06 Np=5.0E+06

Log N

67

2.5.1

Chaoitre

2. Mthnde dc nr,visinn c In rhre

p vie pn fntiotto

La courte derfrence la rsistance fatigue de en

La limite de la rsistance fatgue en

7nnefr,
N c = 2 . 0 8 + M / V o= 5 . 0 E + 0 6

b)

Nr =10.0E+06 Log N

Figure Z-34Diagramme de rsistanceen fatigue de la mthode de I'Eurocode 3 a) Diagramme avec une ligne b) Diagramme avec deux lignes.

L'utilisation d'une ou de deuxlignesdansI'Eurocode3 de la responsabilit I'utilisateur. est de DanslaFigure 2-34a,deuxzonesdiffrentes existent. Dansla zonel,la valeurdepenteest gale 3. Normalement, cettezoneestgnralement utilisedansla mthode I'Eurocode3. de La zoneII estla zonedersistance fatigue. en Dansla Figure2-34b,trois zonesspcifiques peuttredistingues. type de diagramme Ce de I'Eurocode3est plus compliquque celui prsent dansla Figure 2-34a. Le domainede rsistance fatigueestdivis enzoneI et zoneII. en Dansla zoneI, la pentea la mmevaleurque danslazone I de la Figure2-34a.Dansla zone II, Figure2-34b,la penteestmodifiepourtrouverun dimensionnement conomique. plus I problme fatiguedansla zonem n'estpasle plusimportant. de L'Eurocode3proposeune seulerelation utilisabledanstous les cas (Dans les diagrammes avecune ligneet deuxlignes): TtLo,ou<LI* Tu e-r3g)

o, yt, Tu, Lo,ouet Lo* sont respectivement coefficientde scuritdansle domaine le quivalent contraintes, coefficientde scurit la rsistance fatigue,l'tendue des le sur en des contraintes quivalentes l'tendue rsistances fatigue. et des en

68

la Le coefficient de scuritT1 augmente valeurde I'amplitude@uivalentedescontraintes diminuer I danslesjoints souds. coefficientde scurity* fadrt de pour desraisons scurit en la rsistance fatiguetel qu'il estprvupar I'Eurocode3. pour chaque type de dansI'Eurocode3 sontdiffrentes de Les valeursdescoefficients scurit utiliss. et soudage d'aciers

a)

Lon

Lot Lo,

b)

N.

N2

N,

ivr N

a) appliques Spectredes containtes appliques des cumulatif desamplitudes contraintes Figure2-35 Schma de dansle diagramme Whler' des des b) Schma amplitudes conftaintes

69

2.5.1

Chnpite 2. Mthode de prvision de la dure de vie en .fatisue

n'est pasconstante. applique joints souds aux Dansles casrels,I'amplitudedescontraintes La mthode de Palmgren-Miner,proposepour ce type de problmes,est base sur par de I'endommagement cumulatif.Cetteloi est obtenue I'hypothse la sommationlinaire pour de domaine contrainte. desendommagements chaque de s'exprimepar: mathmatique la mthode Palmgren-Miner de L'expression

n, D=\ +L+b+...+nk-=i Nr N2 N3 N- frN'

(2-r3e)

o ni et N, sont respectivement le nombre de cycles avec amplitude constante des contraintesL,o,et le nombre totale de cycles d'endommagementsous amplitude constantedes contrainteso, . Un dimensionnementacceptableavec la mthodede Palmgren-Miner conduit la condition :

D<l

(2-r40)

quivalentes Loru est utilis gnralement dansle I-e conceptdu domainedes contraintes pour les joints souds. concept Ce dansI'Eurocode3 ECCS.C'est une mthodede rfrence par propose Palmgren-Miner. dfinition du domainedes La est plus simpleque la mthode de contraintes quivalentes dansun diagramme Whler avecune seuleligne s'crit corlme suit :
I f

=laLo"nu

li",oor1'
)l
|

(2-l4l)

lI,, t '

le des quivalentes o, Ao,o,, ni, Loi, k et m sontrespectivement domaine contraintes Ao,, les domaines appliques, nombrede cyclesavecamplitude le constante contraintes des constante la pentedu diagramme et correspondants n,,le nombrede cyclesavecamplitude (m deWhlercorrespondante I'Eurocode3 = 3 ou m=5). dansun diagramme Whler deux L'expressiondu domainequivalentdescontraintes de une dans lignesdansI'Eurocode3 estplus complexe. L'Eurocode3 propose relationgnrale (2-t4l)).Selon le diagramme I'Eurocode I'amplinrdedecontrainte tousles cas(quation de 3, dans correspondantecinq millions de cyclesestreprsente la Figure2-36 par o, qui est Des valeurs suffrieures ou une limite pour les domainesde contraintedu diagramme. Loi infrieures L,oo sontprsentes appeles respectivement et Lo,. et

70

2.5.1

de Chnpitre2. Mthode prvisionde la durede vie en fafieue

^ lLo, L6=too,

pour N < N, pour N>N,

(2-r42)

est Palmgren-Miner la sommation basesur la mthode La dfinition de I'endommagement normal de I'Eurocode3, desrapportsentre n, etN,. Dansle casd'utilisation du diagramme (2 : l'quation - 142)s'critci-dessous
ivslvo N>ilo

o=*L*i3
N,
Figure2-36.

(2-r43)

HN,

normal sontprsentes dansla au I-es quations Basquin[44] correspondantes diagramme de

LoSL

' ou =(tdr.rl(lo) Lo= no'r.r(N)-t't 7r7

Lo,

t L,o=to'r,r(Nl''tot,u =(t4,r)P(oo)

+
Loi

Nc =2.08+06 ND =5.08+06

LoS N

Figure 2-36 Diagramme Normal de I'Eurocode3 avec les limies des domaines de containte LOo correspondantcinq millionsde cycles.

prsentes dans la Figure2-36 dansl'quation (2-143),donne La substitutiondes quations (2-144): l'quation

NS/VD lv>flo

' D=iL*iL=i N, N, # ?,

ft

* , *i,1, =.o ' a\,o,3 utolt'

ft

a- ( L ,o ", r ) t,b =( Lo /s ) t

( z - l4 4)

2.5.1

Clnnitre 2. Mthode de orvision de la dure de vie en fatisue

L'quation(2-L44)s'crit :

ft ,=+l fry+fn,rcll blf e F. 'l

lu.J
Dans la Figure 2-36,|e point d'intersectiondeslignes donnela relation ci-dessous : No=a\oo-3 =bLoos et encore : al

,2-t4s
(2-146)

b= Lor,

e-r47)

(2-147)et (2-146)dansI'quatton (2-145)conduit la relation: L'utilisation desquations

p =-9-.
NoLo"o'

o g=fn,nol +fn,d,"r(*)'

(2-148)

I calculpar la mthode I'Eurocode3 peuttreralisarbitrairement partir de l'quation de (2-143) (2-148). ou (2-143)et (2-148)donnentle calcul de I'endommagement. est ncessaire Les quations Il de calculer les paramtres Lo"ou et Lo* (2-138). I'Eurocode par I'equation 3 I domaine quivalent contraintes doncdfini comme: des est D= N <l
N "nu

dans le cas de I'application de la mthodede

(2-r4e)
et le nombre
j=l

o N et N,*sont respectivement nombre le totalede cycles 1f,l=fn,*fn,)

i=l

quivalent cyclesd'endommagement fatigue. aux en Le nombrede cyclesd'endommagement fatigue N,n, pour m=3 s'crit: en
kt

'ru--

l,n,+ln,
-A' j=r

"Lo-,;" (2-144)et (2-150) donne: L'utilisation desquations

KT

a\,o,j,

(2-rs0)

ln,+ln, o= , E =-ft= *f!r, a\o; aAoi bLo

fr o=,

(2-rsr) oo a=(L,o'r.r)r,b=(ao'r.)'
,

72

2.5.1

Chapitre 2. Mthode de prvision de Ia dure de vie en fafigue

Le domainequivalent des contraintesL,o*, sera changen raison de I'utilisation de m=3: Lo"n, + Lo"qus

(2-rs2)

(2-L46),(2-148)et (2-152) La relation (2-153)est donc critepar I'utilisation desquations (2-15 1). dans l'quation

= Lolo,,s

kl

(2-rs3)

ln,+ln, j=r i=l

La relation(2-153)peuttresimplifiecomme:

Lo"q*z=l o , -

[^]*

I In,*I', I I i=r L t=l

o g=fn,toi+fn,t"r(*)'

(2-rs4)

Un autre paramtre ncessairepour le calcul est o*. fatigue Lo"nu.

Il lui correspond la rsistance en

par Lo"q,,t et L'o*,, Dans lecas de m=3, les paramtresAo,o, et L,o* sontreprsents respectivement. o*., s'exprime par :

o"r(+)u lo*,r(lr) = Loo(N = rc'r,r--!4't4o*!+Lo*o,r= o)I (2-154) (2-t55)s'crit: et entre l'quation Le rapport

(2-rss)

(2-ts6)

L'utilisation

des quations (2-140), (2-148)et

(2-L56)donne I'expression de

dans I'endommagement le casd'un diagramme deuxlignesde I'Eurocode3.

73

2.5.1

Clwpitre2. Mthodedeprvisiondeladuredevieen

fatisue

tl

13

n =l1or"''| =-9[ o*,, )

A,o'oN o

31

(z-rsl)

peuttreutilisepowm = 5 . Dansce cas,l'quation(2-150)s'crit: Cetteprocdure

kl

D=

:\ - - i=t ., i=t a\,o"], o\o-

r'

l'.,+ln,

(2-15g)

L'expression du conceptde domaine quivalent descontraints donne les quationssuivantes :

ln,+ln, o= =- 5! =f,+.22=*, oa=(Lo'r,)',6=1\,o'1,)s uia b\,o"l, bLoj" N r\,oi' fr aLo,'

(2-Lse)
ktkt

p-i=t

ln,+l n,

i=t

ba,o"j,

Q -

N,L,o3o

=od ,p-

b.:

ln ,+l n, i=r
i=j

b\,o"j,

ba,o-j

= + --l - equ. Q , L o z , 'o . =

f =t n,* ij=t r,

(2-r60)
En simplifiant:

'li
(2-16L)

=l Lo,qu.s, noo'? |

lI",*In, l I i=r \,=l

oo"f&)i '[ryJ
(2-162\

(2-155) Lo*,, peuts'crire comme l'quation :

na*,r(iv)i= Loo(N = tdr.r-g;-'-'+Lo*,r= o)i

(2-161) (2-162) (2-163) L'utilisationdesquations et donne l'quation :

z r5

D=lao"o*t I=
[ o*r,, )

L,o'oN o

=t

(2-163)

par: I-e paramtre dansle casgnral D s'exprime

74

2.5.1

Chapitre 2. Mthode de prvision de Ia dure de vie en fartgue

,=J
t
[ff)'=''

N<N. [fri)'=',
N>N,

(2-r64)

Dansla Figure2-37,l'algorithme gnral dimensionnement I'Eurocode3 prsent. par du est Les facteurs scurit de sontissusde I'Eurocode3.

L' utilisationdu diagramme normalavecuneligne

L' utilisationdu diagramme normalavecdeuxlignes

_l r=r

L;''1

[",oal' I Lo*,r=o"r(+)t

(o\' . ao,qu,=(.;,J
L,o*n,r= L.o

Lo"qu,s = Lana.s =

( o\' tNJ
( Lorl

roo,*<L9,.;

Figure 2-37 Algorithmedu dimensionnement Eurocode3 par pour un diagramme une et deux lignesS-N.

Les paramtresAo"n,,^ L,o*,, peuventtrecalculs partir de I'utilisationdu diagramme et (2-154),(2-155),(2-16l) et (2-162). normalavecuneet deuxligneset partir desquations La mthodede I'Eurocode3 une mthode prvisionde la durede vie en fatiguebase est de sur les mthodes contraintes. de Cette mthode considre I'effet des contraintes moyennes pffiLo"qu.^. Il existe des diffrentesmthodes dans le cas de la prsence contraintes de movennes.

75

2.5-l

Clnpilre2. Mthodedeprvisiondeladuredevieen

fatipe

Dansla Figure2-38,les mthodes Goodman, de Gerberet Soderberg sont prsentes dansle diagramme ao-o^ [49]. I-esquations Goodman, de Gerberet Soderberg s'crivent:

o"=or^(t 3)
(r.r)

(2-r6s)

oo= onl t-l % | | -(.

\o") )

(2-r66)

o"=o,^(l-%l -\
o,,)

(2_r67)

Les quations(2-165), (2-166) et (2-167) sont respecrivement relations de Goodman, les Gerberet Soderberg.
oo, o^, oyu, o^ et or" sont respectvement I'amplitude des contraintes, la contrainte

moyenne, la rsistanceen fatigue coffespondant l'amplitude des contraintes dans le cas o^ =O, la rsistance tractionet la limit d'lasticit. en

oa

[*)']

Goodman 'r"=rr(r-fr)

ow

o,,

om

Figure 2-38 Schma mthodes la conhaintedansle casd'existenced'une containte moyenne. de.s de

La comparaison entreles mthodes Goodman, de Gerberet Soderberg montrequela mthode Soderberg la plus conservative. est Ces mthodes considrent contraintemoyennedans la

76

2.5.1

Chnpitre2. Mthodedeprvisiondeladuredevieen fafipue

localepeuventaussi de I'utilisation d'une approche contrainte.Des mthodes dformation en treutilises.

77

2.5.2Mthodesde dformationlocale
qui les de Les mthodes dformationlocaleutilisent gnralement dformations existentau dansle cas de la fatigue gnralement proposes sont des voisinage entailles.Cesmthodes au tre aussiappliques casde la localepeuvent de Les oligocyclique. mthodes dformation fatigue grandnombrede cycles. La et lastiques plastiques. composante totaleest la sommedesdformations La dformation et lastique en fonctiondu nombrede cycles: s'exprime partir de la contrainte lastique

ryo +, '2
Ac

oo=o't(zN)b

(2-r68)

E, ao, o'r, N, et b sontrespectivementl'amplitudedeladformationlastiq

en le le modulede Young, I'amplitudede la contrainte, coefficientde rsistance fatigue,le en de nombrede cycles ruptureet I'exposant la rsistance fatigue. plastique la dformation s'exprimepar la relationde Coffin 145-481: de La composante
r

- e',(ZN ,)"
J' r'

(2-t6g)

a ",
2.

J-

l'amplitude de la dformationplastique,le 'r,N J ret c sont respectivement

de coefficientde la ductilit en fatigue,le nombrede cycles ruptureet I'exposant la ductilit en fatigue. entre-0.50 des L'exposantde la ductilit en fatigue c prendgnralement valeurscomprises et -0.70[49]. et est simple descourbeslastiques plastiques montr.La Dans la Figure 2-39, un schma de superposition la courbe lastiqueet de la courbe plastique dans un diagrammebilogarithmiquedonne la courbe totale. Cette sommations'exprimesous forme de loi de Manson-Coffin:

(2-r70)
utilises[49] : sont suivantes gnralement ductiles,les approximations Pourles matriaux o't = ot,| = r par proposes Monow [50] : Les coefficients et c sontdfinispar lesquations Q-llI)

78

2.5.2

Clmpitre 2. Mthode de prvision de la dure de vie en fatigue

-1
c=-

I+5n' -n' . l+5n'

(2-r72)

b=n'c=

(2-r73)

Log@")
s'l

o', E 1 Elastique

e o = e a e + e a , p=212Nr)b

+e!r(2N)"

To"tale

1""

,,,=+-%=LpnSu
2N 1,i

Plastique

=+= e'r(zN ,,.0 )"

Ing(2N 1)

Figure 2-39 Schma la superposition courbeslastiqueet plastiquepour obtenir la courbede dformation de des totaleen fatigue.

o dr , J et z' sont respectivement contraintede rupture,la dformation ruptureet la I' exposant croui d' ssage dformati cyclique. en on L'utilisation de l'quation de Manson-Coffin est la mthode locale la plus gnrale. L'utilisation de l'quation de Manson-CoffinentraneI'utilisation d'une auffe hlpothse dansla mthodede dformation locale.Cette hypothse l'galit de la durede vie en est fatigued'une prouvette lisse et d'une prouvette entaillesoumise la mmeamplitudede dformationavecle mmematriau. Dansla Figure2-40, un schma cettehypothse de est prsent. La mthodedes ElmentsFinis permet le calcul des contraintes des dformations et au voisinagede I'entaille et est basesur le comportement matriaux.Les mthodes des de

79

2.5.2

Chapitre 2. Mthode de prvisionde Ia dure de vie en Tatigue

dformation localeproposent relationsbases le comportement matriaux fond des sur des en d'entaillesonttrsrapides miseen uvre,ce qui en fait leur intrt. de

H tl
N,,

_q

N,

lv

Nc, = Ncz -- N,

Figure2-40 Illustration la mthode dformation de de locale[84].

o Lamthode de linarisation La mthodede linarisationest basesur l'galit du facteur de concentration de contraintes du facteurde concentration dformations. relationpropos": et de La dans la mthode linarisation s'exprimepar : de k,=k"

(2-r74)

o ( et fr, sontrespectivement facteurde concentration contraintes le facteur le de et deconcentration dformations. de (2-17 s'exprime: L' quation 4) *"& k,=kr+e"=k,,

(2-r7s)

o eo et s sontrespectivement I'amplitudede dformation la dformation globale. et I-es rsultatsobtenuspar la mthodede linarisation sont correctsdansle cas de la '.' dformationplane. La mthodede Neuber [9] La mthode Neubercommed'autresmthodes dformation de de localeproposeune qui relationapproximative voisinage entailles s'crit : au des
kok, = ft2

(2-r76)

2.5.2

Clnpitre 2. Mthode de prvision de Ia dure de vie en .fatieue

le o k,, k, et ko sont respectivement facteur de concentrationd" ,ontruinr", et lasto-plastiques le facteur de le lastiques, facteurde concentration dformations de I-e de concentrationde contrainteslasto-plastiques. facteur de concentration est contrainte lasto-plastiques doncdfini: ko = 3 o,
fJ

(2_177)

en I'amplitudede contraintelasto-plastique fond o o, et os sont respectivement globale. d'entailleet la contrainte pour le facteurde concentration dformations de L'utilisation desrelationsproposes de lasto-plastiques le facteur de concentration contraintesen fatigue dans la et (2-177): (2-L76)donne l'quation relation

ii=k?

ou ooo=klo,e,

(2-r78)

L'quation (2-178) contient deux variables inconnues.Pour la rsoudre,il est ncessaire d'utiliser une autre relation qui est la relation entre l'amplitude des et contraintes celle desdformations:

= * =?.(?)' o o," o,p

(2-r7e)

I'amplitude I'amplitudedesdformations, o eo, oo, E , K' , fl' sontrespectivement des contraintes,le module de Young, le coefficient d'crouissagecyclique et cyclique. l' exposant crouissage d'
a

(sY
\K')

k,a,

irl,oi1

ii
l,A(",o;)....'
i--:/ i ./i tl

lr. v. \ 1., :J.--^ = -l/cd-3 lk.c |'= - k: of , '\. - 2' 2\' i**
ti

li./

iI

1..

iJ/,dFc=;ooeo iz L'hypothsede la mthodede Neuber

i../i

s-*\s*ro

ic(s;,o)
s Krg ee o

de de Figure241 Illustrationde I'interprtation la mthode Neuber.

8l

2.5.2

Clnpitre 2. Mthode de prvision de la dure de vie en latigue

Dansla Figure 24L,l'interprtation de la mthode Neuberestprsente.L'galit de dessurfacesSoa", et S-,ro donnela relationdeNeuber. Neuberconsidreque l'nergie lastiquedansle cas lastique gale l'nergie est dansle casrel,pourles mmes valeurs amplitudes contraintes dformations. des de et

L'algorithmede la mthode Neuberdansle casgnral prsent de est danslaEigne2-42.

ooo=klore,

LE

A -i-'l;
Ero'r relrrbrc

Hfits -"n\-A

ILE
oarEa

d. eo=J(2Nr)" +er(2Nr)
E

Figure2-42 Ngonthmegnral d'utilisationde la mthode Neuber. de

par Les rsultats obtenus la mthode Neubersontgnralement de bonsdansle casde planes. contraintes

La mthodede Molski-Glinka l5l-521 La mthode Molski-Glinkacommed'autresmthodes dformation de localepropose uneapproximation voisinage entailles s'critcommesuit : au des et

(k+s) =4 *-"" (".\*'

2E 2E l+ ti\K' )

(2-180)

82

2.5.2

de Chapitre2. Mthode prvisionde la dured.evie en .fatigue

le o fr,, as, E, oo, n' et K' sont respectivement facteur de concentration de contraintes, la contrainte globale, le module de Young, I'amplitude des contraintes, I'exposant d'crouissagecyclique et le coefficient d'crouissagecyclique. L'quation (2-lS0) contient deux variables inconnues. Pour la rsoudre il faut aussi faire appelle la relation entre I'amplitude des contraintes et des dformations (2-

r79\.
Dans la Figure 2-43,l'interprtationde la mthodede Molski-Glinka est prsente. L'galitdessurfacesSoo,r,et 56o* donnela relationde Molski-Glinka.

oo = Eeo

-oo

E k,o
oa

r(

4)7

K)

B(k,er,k,or)

s- o= or\*,r = orr, r) )n7 " I@,

L'hydothsede la mthodede Molski-Glinka

,ooroL,*ro

de de Figure2-43 Illustrationde I'interprtation la mthode Molski-Glinka.

Dansla Figure2-43,lessurface S*,ro et S6*o peuts'crire:

s*,ro -

1^

,kiorer

( 2 - 18 1)

="To"r+.+(+)' soo,o ,")I-'yo"

s =n!,?. @po. .#(+)+ = .hQal' )* # v:fr # o,o,o +W)+ =

(2-182) 83

2.5.2

Clnpitre 2. Mthode de prvision de la dure de vie en fatisue

(2-181), (2-182) e, =| I-esquations et

I

conduisent l'quation suivante :

oI' _(k,o)2 *_(\7 2E l+ n'\K' ) 2E

oi estuneamplitude contrainte de arbitraire.

(2-r83)

L'quation(2-183) l'quation la mthode Molski-Glinka. est de de Dansl'quaton (2-183),

L'algorithmede la mthode Molski-Glinkaestdonndansla Figure2-42 danslaquelleon de remplace l'quationde Neuberpar cellede Molski-Glinka. I-esrsultats par obtenus la mthode Molski-Glinkadonneune estimation de situeentreles mthodes linarisation de Neuber.Gnralement, mthodes linarisationet de de et les de Neuberconduisent desvaleurslimites.Dansla Figure244,les mthodes leurs quations et (k,er,k,or) . sontprsentes. calculs dcroissentpartirdu pointB de la courbe Is lastique Les courbesdes mthodesde linarisation,Neuber et Molski-Glinka coupent la courbe reprsentative comportementdes matriaux respectivement du aux pointsl(k,er,o't), ni'G,oi'G1 et F(e:'*,oI'*). La mthodede Molski-Glinkadonne gnralement des plus que estimations exactes lesmthodes linarisation Neuber. de et par diffrents auteurssont prsentes. Dans le Tableau2-5, les relationsproposes Is mthodes Linarisation, de Neuberet Molski-Glinka sont utilisesen I'absenced'effet de contrainte moyenne.
oa
oo= Eo

(p)'

qeo=(ore,

4
2E
C;.=I,

Ktc

o'"

e;"

Figrue24

Comparaison mthodes des linaires,Neuberet Molski-Glinka.

84

2.5.2

de de Chapitre2. Mthode prvision la durede vie en fatisue

La prise en compte des effets de contrainte moyenne dans la mthode de dformation locale est trs discutable et difficile.

par pour ft, proposes diffrents entre , ko et k" . auteurs la relation Tableau Equations 2-5

AUTEUR Neuber[19] Makhutov[53]

RELATION PROPOSEE kl = kok,

*:=n,,(T)'1, (#)]
l+n'

Molski-Glinka

t5r-szl k?

kok,

#.h(9'l #'"'(*o)'
t=T*

rl

l+.".(?)r
Morozov-Pluvinage

Ls4l

peuttrepris en comptepourestimerla moyennes Dansce cas,I'effet d'une desdformations par : durede vie. Les relationssuivantes propos{es diffrentsauteurs sont o La mthode Morrow [55] : de

t,=T(zN)b

o" -o

+;(2N)"

(2-rB4)

o Lamthodede Smith-Watson-Topper (SWT) [56] : (2N eao*xE = (o'r)2 ,)'u * o're'rE(2N , o*, = on + oa ,)u*" (2-ls5)

globalemaximale,la contrainte globale la o oor*, onin et o. sontrespectivement contrainte minimaleet la contrainte moyenne. prsentes Les mthodes dans les quations(2-184) et (2-185) sont normalement utilises dansle casde la fatigueoligocyclique.

85

2.6Rsum
et de Dans ce chapitre,le facteur de concentration contraintes le facteur de rductionen autourd'un par I sontprsents. champde contraintes fatigueproposs diffrents auteurs et plaqueest vrifi par I'approcheanalytique la simulationnumrique. trou dansune grande Le calcul par ElmentsFinis est ralisavecle logiciel ANSYS. On a trouv le champde prsent 2-2. dansla Figtlr:e pourdiffrentscasdechargement contraintes dansles la respectivement distributiondes contraintes Les Figures 2-3 et 2-4 reprsentent directions et y. x Nous avons comparle champ des contraintesobtenuspar I'approche analytiqueet la par Finis danslesFigures2-24,2-25et2-26. Le champdescontraintes mthode Elments des des la mthode ElmentsFinis pour une barreest montrdansla Figure 2-27. Cettefigure pour diffrentesgomtries. Nous avonscalculle prsente concentration contraintes de la pourcetteprouvette. de facteurde concentraton contraintes est dansle et dansle caslastique lasto-plastique considre La distributiondescontraintes proposes notrersultatobtenupar la et formulations les chapitre On a compar diffrentes 2. Finis dansla Figure2-7 et2-8. mthode Elments des est en La distributiondescontraintes fond d'entailledansle caslasto-plastque tudiepar les de Dansle casnumrique, techniques calculpar et les approches analytiques numriques. Elments Finis sontprsentes. par dansla section est I facteurde rductionen fatigueprsent diffrentsauteurs considr dansles localesont prsentes 3 de 2.4.La mthode Errrocode et les mthodes dformation sections suivantes. Kuhn Hardraht et Topper-El Haddad Les mthodesPeterson,Neuber, Stieler-Siebel, par prcisent la connaissance rayon d'entaille.Dans le cas desjoints souds point, la du dfinitionclairepour rayon d'entailleet la courbelissen'existentpas.Il est doncvidentque mentionnes pasenvisageable. n'est I'utilisationdesmthodes dansles structures acier. en communs La mthode Eurocode est limite auxjoints souds 3 par ne Lescoefficients scurit sontpasdfinisdansle casdejoints souds point. de

86

CHAPITRE 3

87

Chapitre 3
3. Mthode volumtrique pour la prvision de la dure de vie de structuresentailles
3.0Introduction gnrale
La mthode volumtrique une approche est semilocalequi utilise le conceptde I'intensitdu (stress champdes contraintes field intensityapproach, SFI) pour pvoir la durede vie de structure entaillet57-581. DansI'approche I'intensitdu champdescontraintes, rupture de la par fatigue est ralisepar I'accumulationde I'endommagement dans la zone plastique. L'endommagement dans la zone plastifieau voisinagedu fond d'entaille dpendde la distributiondescontraintes lasto-plastiques du gradientde contraintes et simultanment. Sa valeurcritique pour la rupturepar fatiguedansla zoneplastifieest donc caractristique. La contrainteeffective du champ des contraintesest calculepar I'intensit du champ des contraintes.

o!,,=| !f",>oe)dv

( 3- 1)

o oS, ,V ,f(or) ,0(7) et ) sontrespectivement contrainte propospr la mthode la d'intensit du champ des contrainte,le volume de la zone plastique, la fonction des contraintes, fonctionpoidset le volumed'intgration. la La contrainte critiquedpend matriau de soncomportement du et dansle casd'une grande concentration contrainte. de Dansla Figure 3-1, un schma simplede la mthode d'intensit du champdescontraintes prsent est dansdescasgnraux gomtrie de chargement. de et La mthoded'intensitdu champdes contraintes donc une mthodesemi locale.Cette est mthodencessite connatrele champ des contraintes voisinagedu fond d'entaille. de au Ensuite,le volumede la zone plastifieest considr pour calculerI'intgralede l'quation (3-l). Cettemthode ncessite dfinir deux fonctions,la fonction descontraintes la fonction de et poids. 88

3.0

Chapitre 3. Mthade volumtrtaue pour la prvision de Ia dure de vie dc structures entailles

H lli
l.

,q
I

, !

T-- j i
\--

..'1-..-."

La surfaced'entaille

Figure 3-l

Prsentation de la distribution schmatique des contraintes en fond d'entaille pour diffrents tlpes

d'prouvettes.

est Le facteurde rductionen fatigue par la mthoded'intensitdu champdes contraintes prsent:

o,=

!r<t)oe)dv

(3-2)

le o , et f (t)sont respectivement facteurde rductionen fatigue et la fonction des globale. par contraintes normalises la contrainte du avec de Le volume d'laboration processus fatigueau voisinage fond d'entaillechange du planeset de dformations Dansle cas de contraintes la gomtrie le type du chargement. et planes,le volume reprsente zone plastifiedans le cas bi-dimensionnel.Dans d'autres la cas,une valeur prcisedu volume de la zone plastifieen fond d'entaille ne peut pas tre en prouvettes fond d'entaille avec obtenue.Dans la Figure 3-1, les volumesde quelques des gomtries Selonles difficultslies la configuration volumes, diffrentes sontmontrs. d'intensitdu champdescontraintes. il est ncessaire simplifier I'intgralede la mthode de et le respectivement cas bi-dimensionnel uniLes quations(3-3) et (34) reprsentent dimensionnel:

89

3.0

Chapitre3. Mthode pour la prvisionde la durede vie de structures volumtrtaue entailles

o,=*!rrQe)da k,=!tlrou>oa>a*

(3-3) (3-4)

o A, D ,l ,L,0(7) et Q@) sontrespectivement valeurde la surface la zoneplastique, la de la surface d'intgration, longueurde lazone plastique, longueur la la d'intgration,la fonction poidsdansle casbi-dimensionnel la fonctionpoidsdansle casuni-dimensionnel. et Dansle casuni dimensionnel partirde la dfinitronf () =+ 4) s'critcomme dessous ci :
Itc
6,

et Q@)=l, la relation(3-

, o"tr lo'<aa* ooo"r 1' =7 ot=l=f,l !""<.>*

(3-s)

o frr , or, or(x) et o,r sontrespectivement facteurde rduction fatigue,la contrainte le en globale, contrainte la d'ouverture fatigueet la contrainte en effectiveen fatigue. L'quation (3-5) reprsentela relation du facteur de rduction en fatigue base sur I'hlpothsede la mthode la contrainte de effective34,19,59-621. La contraintemoyenneeffective dans l'quation (3-5) peut s'exprimer par la mthode Peterson [61]: (3-6) Ainsi, selonla mthode Peterson, facteurderduction fatigues'exprimepar: de le en
Kf =-

'

o *( x')

'os

(3-7)
sont respectivement distance du fond d'entaille et la contrainte la

o .r'et or(x')

correspondantela distancex.. Dansla Figure (3-2),le schma mthodes la contrainte des de effectiveet de la mthode de Peterson prsent. est Dansla mthode d'intensitdu champdescontraintes, fonctiondes Ia contraintesf (o,i) dpenddu tenseurde contraintes [63]. La dfinition de la fonction des contraintes dansla mthode d'intensit champdescontraintes du s'crit 164-651:

3.0

Chapitre 3. Mthode volumtrtaue pour Ia prvision de la dure de vie de structures entailles oD@

=r"Soorr'r*-=trfurUr* o* o,.txo

a)

b)
du Figure 3-2 Schma de I'intgrationau voisinage fond d'entaille a) Mthode de la contrainte effective b) Mthodede Peterson. T @ = rl o1 + o', - oxoy +3rI = o'"0

(3-8)

est La fonction poids dans la mthoded'intensitdu champdes contraintes gnralement : dfiniecommeci-dessous Q (i )= L - c r ( 1 + si n (9 ))

(3-e)

de la o /(i) , c , r , et 0 sontrespectivement fonctionpoids,la constante la fonctionpoids polaires. et lescomposantes Cettefonction appliquni de la gomtrie. La fonction poids ne dpendni du chargement des dpenddu comportement matriaux.La fonction poids permet une discriminationdes points les plus importantsdans le calcul de I'intgral. Cette fonction est soumise des particulires : conditions 91

3.0

Chapitre 3. Mthode volumtriaue pour Ia prvision de la dure de vie d.estructures entailles

1. o30(v)<r
2. Si G=Vo(O)=0=/(0)=1 3. SiG=io(r)=O+Q(7)=l :. s'exprime
I 6' zqv t -l t 0?,0)= lr(l + sind) L on'o) f ,,

(3-10) (3-11) (3-r2)

la d'intensit champ contraintes, relation(3-9) du des Selonla dfinitionde la mthode

(3-13)

o o'eqn,6p*k, r , et d sontrespectivement contrainte la quivalente dansl'lment<< >>, i la polaires. contrainte maximale les composantes et Dansla Figure3-3,Ie diagramme schmatique la fonctionpoidsdansle i melment de est prsent. maillagescorrespondent un grand trou dans une plaque pour le cas d'une Is planes. hypothse contraintes de

lil= 0,o ** ,e r-, Figure 3-3 Schma de la fonction de poids au voisinagedu fond d'entaille dansla mthoded'intensitdu champdesconffaintes.

Dans la Figure 3-3,la contraintemaximaleet la dformation maximalese trouventen fond d'entaille.La fonction poids dpend donc de la valeurmoyenne danschaquelment. est Il d'intgrerla contrainteutilise dansle conceptd'intensitdes contraintes donc ncessaire sousformediscte:

=+|"'*r(r)s, on
92

(3-14)

3.0

Chaoitre 3. Mthode volumtriaue pour la prvision de Ia dure de vie de structures entailles

o o", , 5 , o'a", QQ), ^S, n sontrespectivement contrainte la mthoded'intensit et la de du champdescontraintes, surfacede I'intgration,la contrainte la quivalente dansl'lment << >>, fonction poids dansl'lment<< >>, surface i la i la dansl'lment << >> le nombredes i et maillages. La procdure mentionne peut s'crire dans la mthode d'intensit du champ des dformations[63]. Les quations(3-8), (3-13) et (3-14) sont transformes sous la forme suivante: (3-15)

(3-16)

(3-17)
La mthode volumtrique utilise le concept de la mthode d'intensit du champ des contraintes avec diffrentes fonctions poids dans le cas des calculs lasto-plastiques.La mthode volumtrique et la mthode d'intensit du champ des contraintes sont des mthodes semi locales. La moyenne des contraintes dans le volume d'laboration est la contrainte effective. Le volume d'laboration est dfini soit par la zone plastifie (M.I.F), soit par le minimum du gradient des contraintes(mthodevolumtrique).

93

3.L Ltapprochevolumtrique
Dansle casde I'existence du d'unezoneplastifieau voisinage fond d'entaille,la distribution prsente caractre descontraintes lasto-plastique. Dansles cassimples,cettedistribution un il est donnepar une solutionanalyique; dans les cas plus compliqus, est ncessaire en d'obtenir le champdes contraintes lasto-plastiques fond d'entaillepar I'utilisation de la par mthode Elments Finis. Les rsultats obtenus les calculsaux ElmentsFinis et les des que la contrainte lasto-plastique maximalen'est pasen fond approches analytiques montrent Au d'entaille mais une certainedistancede celui-ci 124-261. del de ce maximum, la distribution est souvent de type puissanceet gouvernepar le facteur d'intensit de contraintes d'entailleet la distanceX"u estdfiniede faonparticulire. La mthodevolumtrique effectiveet le facteurde rductionen fatigue calculela contrainte par les quations suivantes [66-67]:
X'

o"tr=!' ct

lo-<r)(t- xl)drc, o x"n d ,o' "^

I "^*",
x"

'v

oo(x)

d-r

,. -o" "'- o^-

to,@)(r-xz)drc 0

(3-1e)

o o, , on, X"u, or(x)

et Z sont respectivementla contrainte effective, la contrainte

nominale, la distance effective (selon cette mthode), la contrainte d'ouverture en fatigue, le gradient de contraintes. La comparaison entre l'quation (3-4) et (3-19) proposes par la mthode d'intensit du champ des contrainteset I'approche volumtrique montre la diffrence entre les dfinitions de la fonction des contrainteset la fonction poids [68-70] :

=o'(x) f (o,,) 'oo * Q@)=r- 2, o l=

I do-(x) on@) dr

(3-20) (3-2r)

la Dansle casde I'existence contraintes cisaillements, relation(3-19)s'exprimesousla des de

forme suivante : x.

k,=4=+^ k,(x)(rxt)dx ' o tnX"a

1 dt-(x) 7 'v/ = - 4 co@) dx

(3-22) (3-23)

94

3.1

pour Chapitre.3. Mthode volumtrique laorvision Ia dure vied.e de de structures mtailles

o T"r, tn, X"r, ro(x) et Z sorftrespectivement cisaillement le effectif, le cisaillement nominal,la distance effective,le cisaillement fatigueet le gradient cisaillements. en des En gnrale, contrainte la principale.Cettecontrainte d'ouverture fatigueestunecontrainte en est perpendiculaire la directionde la propagation fissures. conceptde la contrainte des Le d'ouvertureen fatigue et la direction d'intgrationsont les paramtres plus importants les dansle calcul de la durede vie en fatiguepar la mthode volumtrique. Dansla Figure 34, la contrainted'ouvertureen fatigue et la directiond'intgrationpour diffrentes prouvettes sont schmatiquement prsentes. L'valuation des contraintes au voisinage du fond d'entaille prsente toujoursun maximum dfini prsdu fond d'entaille.La taille de la zone plastifieedpend type de chargement de la gomtrie prouvettes desentailles. du et des et

tl g
% lasto-plastique.

ffi

Figure 3-4 Schma la distributiondes containtesen fond desdiffrentes de entaillesdansle cas d'un calcul

I-e rayon d'entaille changela morphologiede la zone plastifie au voisinagedu fond d'entaille.Plus le rayonaugmente, la zoneplastifie plus augmente vice-versa. calculde et I-e la distributiondescontraintes fond d'entailledpend type d'prouvettes de la valeur en du et du rayond'entaille.

95

3.1

pour Ia orvision la durede vie de structures entailles de volumtrtaue Chapitre3. Mthode

La fonction poids dpend du gradient des contrainteset de la distanceen fond d'entaille. La fonction poids contient les deux premires parties de I'expansion de Taylor. Elle est soumise aux conditionssuivantes:

1. 0<Q@)<r 2. Q(0)=r
3. Si7=0=Q@)=l I calcul de la distance effective est un lment essentiel de la mthode volumtrique et consiste dterminer la position du minimum du gradient des contraintesdans un diagramme biJogarithmique.

o,60

ros(o o(")) = tog(c)- a Log(x) ----o-o

O ,1 O; q
o
a

o
2

A"

-o.40f:
- : l

E
h

^ v
H

o -0,90;.
H

. r-J

c)
\., t

-r,40v

,E

- 1,90

o,1

&nax

x (nlr

xen

&r1,O

et lasto-plastique du de bi-logarithmique la distributionde la contrainte du Figure 3-5 Schma diagramme gradient contraintes fond d'entaille. en des

sont prsents et lasto-plastique le gradientde contraintes Dansla Figure3-5, la contrainte Trois zonesparticiilires biJogarithmiquede la mthodevolumtrique. dans le diagramme : peuvent tredistingues c jusqu' la contrainte en lasto-plastique fond d'entaille augmente Zone I : La contrainte maximaleo,**.

96

3.1

pour la prvisionde Ia durede vie de structures Ch.apitre Mthodevolumtrtaue 3. entailles

o o

Zone fI : La contrainte lasto-plastiquediminue jusqu' la distance effective X

"u

Zone III : L'volution de la contrainte lasto-plastiquea un comportement linaire dans le diagramme logarithmique. Ce comportements'exprime par une loi en puissance[68]:

Log (o r(x)) = In SG) - a Log(x) La relation(3-24)s'crit: o*(x) =a "x'

(-

(3-24)

(3-2s)

o orr(.x) , C et a sont respectivementla contrainte d'ouverture en fatigue, un coefficient du matriau et I'exposant de la pseudo singularit des contraintes.

soir(3-25),
o ,t'(x\ '

K
(2m)"

'

Kp=C(2n)o

(3-26)

Ko estle facteurd'intensit d'entaille 166,68,7I-741. Les coefficients C et q dpendent du comportement du matriau, de la l'entaille et du type du chargement: C = C(E,v,6y,n',K', p,W,a, R) q = a(E,v,oy,fl', K', p,W,a, R) o E, v, 6y, fl', K', p,

gomtriede

(3-27) (3-28)

W, a et R sont respectivement module de young, le le

coefficient de Poisson,la limite d'lasticit,I'exposantd'crouissage cyclique, le coefficient d'crouissage cyclique, le rayon d'entaille, la taille du ligament,la profondeurd'entaille et le rapport de contrainte ( R = o-o o** ).

Les paramtres mentionns sont utiliss dans la mthode des Elments Finis. Les paramtres , v , oy, n' et K' sont intgrs dans la loi de comportement du matriau. E D'autres paramtres (p, W , a , R) caractrisentla gomtrie et le type du chargement

respectivement. Is coefficients C et a sont exclusivement obtenus par calcul par la mthode des Elments Finis. Dans la Figure 3-6,le schmade la zone III est compar la relation conventionnelle la de mcaniquelinaire de la rupture. La distribution des conraintes lasto-plastiques est prsente dans I'intervalle entre la distanceeffective X,u et la distancede la fin de la zone In, xn .

97

3.1

pour la prvisionde Ia durede vie de structures entailles volumtrtaue Chnpitre3. Mthode

I'os@)

Mecaniquelinairede la rupture

7,tal=ffiru<el
=oo(x,) ou(x,)
o"fr
oi

K..,.,,

ffi

= K ,u@ ,y

- yy\--' (2n)" "l<.t=* a<0.5 Kp=C(2n)"

Mcaniquede

X"u1x1xo

0.5

xi

x" I-og(x)

et dans III. des en Figure Distribution contraintes fondd'entaille enfonddefissure la zone 3-6

La distance xnde situe la fin de la zone linaire. I point (*,,o,) est I'intersection des courbesde la mcaniquelinaire de la rupture et la mcaniquede rupture d'entaille. La procdurementionnepeut tre utilise en remplacantla contrainte par la densit d'nergie de dformations(W = [ode) 75-76l.La relation (3-25) s'exprimepar:

Lw = c{
xun

(3-zs)

o I4/, C" et dw sont la densit d'nergie de dformations et les coefficients de la distribution linaire dans le diagrammebiJogarithmique (17 - x).

: des Les coefficients C* et dy, dpendent matriauxet du type du chargement Cw = Cw(E,v ro,n', K', p,W,a,R) = c4o an(E,v,o,n',K',p,W ,a,R)

(3-30) (3-3 1)

le o E, v, oy, fl', K', p, W, a et R sont respectivement modulede Young,le I'exposant d'cruissage cyclique,le coefficient la coefficientde Poisson, limite d'lasticit, d'crouissage cyclique, le rayon de I'entaille, le ligament, la profondeurd'entaille et le rapportde contraintes.

98

3.2 Procdure de prvision de la dure de vie en fatigue par Ia mthodevolumtrique

La mthodevolumtriquencessite connaissance la contrainteeffective et la distance la de effectiveet permetle calcul du facteurde rductionen fatigue.I facteur de rductionen fatigueest utilis dansla procdure prvisionde la durede la durede vie en fatiguepar de la mthodevolumtrique.La procduretraditionnellede prvision de vie en fatigue par I'approchevolumtrique I'utilisation de la mthode graphique. prsentation est graphique La de la mthode volumtrique montre est dansla Figure(3-5). L'volutiondescontraintes fond d'entailleestdivise trois zones. Tableau en en Au 3-1, ces zonessont individuellement prsentes la zone Itr peut tre prsente et ave,c une relation analytique.
Tableau Description zones distribution contraintes 3-l des de des dansI'approche volumtrique.
NO

DESCRIPTION
La conhainteaugmente la valeurma:cimale dtermine,0 ( X ( X.",* et est La contraintediminue et la valeureffective estdtermine,x-* La contraintea un comportement linaire, X "t 3 x 3 xn S* < X "f

EQUATION

I
II

ru

ln g(oo@))= LosG) - a Lo 4,

La mthodegraphique I'approche de volumtrique donnela distance effectivecorrespondant au gradientde la zone linairerespectivement. L'utilisation de la distance effectiveet de la relation(3-19)conduitau facteurde rduction fatigue. en Le facteurde rductionen fatiguecalculpar la mthode volumtrique s'applique la courbe de rsistance fatigue d'une prouvette en lisse dansle diagramme Whler et permetde de retrouverla courbede rsistance fatigued'une prouvette en entaille. Dans la Figure (3-7),I'algorithmed'valuationde la durede vie en fatigueest prsent. La gomtrie,le chargement le comportement matriauxsont introduits dansle calcul et des lasto-plastique la mthodedes ElmentsFinis. I-e comportement matriauxen par des prsence chargements de par cycliques reprsent la relationde Ramberg-Osgood est :

,.=+.(?)'
I

Q-32)

99

3.2

Chnpitre 3. Mthod.e volumtrique pour Ia prvision de la dure de vie de structures entailles

o E , K' et n' sont respectivement module de Young, le coefficient d'crouissage le cycliqueet I' exposant d'crouissage cyclique. L'quation(3-32)reprsente comportement desmatriaux le rel dansle cas de chargement cyclique.

Gomtrie Chargement et

t desMatriaux

Calcul par Elments Finis

Courbeentaille

Figure 3-7 Atgorithme de la mthode volumtrique dans le cas de I'existence de la courbe lisse.

Le facteurde rductionen fatigue,relation (3-33), permet partir de la courbede V/hler d'uneprouvette lissede retrouver celled'uneprouvette entaille.
.oo kro = ] ' N

(3-33)

o k/ , os et dN sont respectivement facteur de rductionen fatigue, la rsistance le en fatigued'une prouvette lisseet la rsistance fatigued'uneprouvette en entaille. L'algorithmementionn dansla Figure 3-7 ne peut tre utilis que si la courbede V{hler d'uneprouvette lisseestdisponible. exemple, Par pour lesjoints souds point I'obtention par de cette courbe n'est pas possible. La mthode volumtriquepropose alors un autre pour la due de vie en fatigue. algorithme Dansla Figure3-8, lalgorithme de la mthode volumtrique prsent. courbelisseest est Ia par obtenue I'utilisation desrsultats exprimentaux d'prouvettes entailles avecles mmes gomtriques. courbelisseobtenue I'utilisation de I'algorithmede la caractristiques La par Figure 3-8, peut tre utilise dansla Figure 3-7 pour des prouvettes entaillesde mme gomtrie maisavecdiffrents rayonsd'entaille.

3.2

Chapite 3. Mthode volamtrique pour la prvision de Ia dure de vie de structures entailles

cat semiJocale, cettemthodeutilise la contrainte La mthode volumtrique une approche est effective. sur effectivecalcule commeunemoyenne la distance

Gomtrieet Chargement

Courbeentaille

Finis Calcul par Elments

Courbelisse virtuelle

Utilisation de I' algorithme dansla Figure 3-7

volumtrique dansle cassansla courbelisse. Figure3-8 Algorithmede la mthode

exprimentaux serontcompars I-es rsultatsobtenuspar la courbeentailleet les rsultats pourexaminer et I'exactitude la fiabilit de cetteapproche. Cettefaonde procder bi-logarithmique. La mthode volumtrique utilise unereprsentation d'utiliser une formulation plus simple. Pour n'est pas facile en pratique.Il est ncessaire en de surmonter difficults,le lissage la courbeest utilis.La distributiondescontraintes ces fond d'entaille doncdfiniepar l'quationQ4D: est

o,(x) --Zo,r'
i=d)

(3-34)

I-e gradientdescontraintess'exprimepar :
y(y\ - L '!\"' = oo@) doyy(x) =ro,*'-' g dx i o'*'

(3-3s)

et la fonctionpoidss'crit:
g

= Q @) l -

xlia,xt-t '-:

Zo'*'
i=0

3.2

Chapitre 3. Mthode volurntrtaue pour la prvision de la dure de vie de structures entailles

La distance effective peut tre obtenuepar I'expansionde Taylor au voisinagedu fond

: de d' entaille. Elle correspond au point minimumd'volutiondu gradient contraintes

az@) o _
dx

(3-37)

L'utilisation de la relation(3-35)dansla relation(3-37)donne:

' -r^. d Z - ; :f(a,i'x-' o''

- o,fut-')
-

dx

Ei
(,,'')'

-0

(3-38)

fo'' i=0

L'utilisation de I'expansionde Taylor avec quatretermeset la relation (3-37) donnentla pourz = 3 , lesrelationsde (3-34) (3-38) s'crivent: effective.Parexemple, distance
or(x)= ao* qx+ a"x' + q*t --, \
)L\4) -

(3-3e)
.'i' (ar+2arx+3arxz)z (ao + qx + arx2 + arx3)2

a" + 2arx +3arx2 ao+ arx+ arx' + arx3 2ar+6arx ao+ a1x+ arxz + q*t _

(3-40) (3-4r)

df(x)_ dJc

L'expansion de Taylor avec quatre parties au voisinagedu fond d'entaille s'exprime:

dz@)=L* drc%%

(2a'-1 3 %-W- - Qq o o. : ) q " oo'** a;

x2+

(3-42)
x3+ O(xa)

(3-42) donnela distance effective partir de L'utilisation de l'quation(3-37)et I'expression I'expansion Taylor en fond d'entaille: de

x"f=ww
volumtrique s'crit: en Le facteurde rduction fatiguepar la mthode

(343)

"r=-

, -1"-?*'"-?*'*)
Zln L02

(344)

La relation (344)

est une approximationdu facteur de rduction en fatigue base sur

de a,, I'approchevolumtrique. Dans les quations (3-39) (34),les coefficients i=1..n

3.2

Chaoitre 3. Mthode volumtrioue pour la prvision de Ia dure de vie de structures entailles

dpendentde la distribution des contraintesen fond d'entaille. Cette formulation de la mthode volumtrique conduitaux avantages suivants: 1. Il est difficile d'obtenir la distributiondes contraintes le gradientde contraintes et partir du diagramme biJogarithmique, dansl'applicationde la mthodegraphique de I' approche volumtrique. 2. La formulationanalytique pluspratiqueet efficace. est 3. La formulation analytiquen'est pas restreinte la distribution des contraintes. La procdure peut tre utilis pour la distributiondes densits mentionne d'nergiede dformations la distributiondesdformations voisinage fond d'entaille. et au du or(x) + IV(x) oo@) + eyy(x)

(34s)
(3-46)

L'expressionprsente dansl'quation (344) est exactepour la distributiondes contraintes avec = 3 . l-e Tableau3-2 reprsente n respectivement distribution,la distance la effectiveet le facteurde rduction fatigue. en Il est ncessaire faire le lissagele plus exactpossible. de L'ordre du dveloppement obtenu dpendde la gomtrie, type du chargement du comportement matriauxutiliss. du et des Au Tableau3-2, les valeur obtenues pour les paramtres X* et k, I'ont t pour I'ordre quatre. exemple distance Par la effectiveet le facteurderduction fatiguedansle casn=2 en

peuvent crits: tre = X,tr !+O@) z4z

(3-47)

kr=T+o(4)

,.-l^-?'')

f^t

(3-48)

A partir de I'ordre cinq dansle Tableau3-2,ladistanceeffectiveresteconstante quel que soit I'ordre polynomial. prsente peutpastreutilisedansle casde I'existence La procdure ne d'un maximumdans la distribution. La distribution des contraintes lasto-plastiques prsentetoujours ce maximum.Il est videntque dansd'autrescas,la formulationprsente devraitpas tre ne utilise. L'expression mathmatique cetterestrictions'exprime: de

103

3.2

Chnpitre 3. Mthode volurntriaue pour Ia prvision de la dure de vie de structures entailles

lo@^) o*'=to(".,.) x**) + x"o.-o ,r*=o = a e i ,*=u

,x,**t0 -

X*uO

(3-4e)

propose pour I'approche volumtrique. Tableau analytique 3-2.Prsentation la mthode de

EQUATION

r@)
rG) =

X
i-l *"f =-;

kr

o.,(r) =

3 ia ; r L
i=0 '

TL

L ltl;X j=0 '

?.

a:x

i=0

kJ

l^?.,]
dn

o rr(x')=

L arx' ' i=0

z@=+

3 .tA:X i - r L
L a;x' i=0

,. = ^"tr

-ar t4o! n"?"1 2r2-E- ,*1 2

kf =

l"-?-'"-?r'"]

o,-.")=,\o,,t

z\xr=-

i-l . L Ul;X ' i=0

s L a;r i=0

.*=;-#.
-.
""f = h2-,*:

-g"?"3, eoo"rol +
tsa5, kI=

?ral 1",-?'?u-?*2u
on

onn(x) = .L^a;x t=u

i=0 zG) =
:i

, ,,rt-t '
L AiX i=0 '

- at t4ol . -n'?"1+ aaoaras,

ooF) = 2 atx' '

i=0

6;

zG) =

i=0 ' 6; E atx' ' i=0

f, ,,,'-t

-. = ^"fr

- ar t4ol

-n'?"1++aoaral

2"2-E'

"4

-? -2 -? -? - :ul V, -b -h -b -b i',T
on

La mme procdurepeut tre utilise en remplaantl'volution des contrainteslastoplastiques l' volutiondesdformations par lasto-plastiques.

104

3.3Rsum
pour la prvisionde la durede vie en fatigueest tudiedansce La mthode volumtrique chapitre.L'approchevolumtriqueest une mthodesemi-locale. facteur de rductionen I-e fatigue est calcul partir de la valeur moyenneou effectivede contraintes dansune zone dfinie par la distanceeffective. Cette mthodencessite connaissance la distance la de effective,la fonctionpoidset la distributiondescontrainteslasto-plastiques voisinage au du plus pour la prvisionde la durede vie en fond d'entailleet est gnralement performante fatigue[70]. La mthodevolumtrique calcule gnralement contrainte la effectiveen fond d'entaille en moyennant queles valeurscalcules pour dansun volume,mais destudes ralises montrent le facteur de rductionen fatigue partir d'un volume,d'une surfaceou d'une ligne sont approximativement mmes[67]. les gnralement un logarithmique.Dans cette La mthode volumtriqueest prsente sur prsentation, qui nousavonspropos nouvelleformulaton analyiquepour cetteapproche une nouspermetde calculerla distance effective,la contrainte effectiveet le facteurde rduction en fatigue partir du lissage rsultats par des obtenus Elments Finis. La distribution des contraintes lasto-plastiques fond d'entaille dpenddu module de en Young, du coefficient de Poisson,de la limite d'lasticit,de I'exposant d'crouissage cyclique, du coefficient d'crouissage cyclique, du rayon d'entaille, de la gomtriede l'prouvette du rapportde contrainte. et L'obtentiond'uneformulationgnrale analytique est donc difficile par calcul. I-a mthodedes ElmentsFinis nous permet d'obtenir cette distribution. La courbe aprs lissage de cette volution des contrainteslasto-plastiques conduit uneexpression analytique. Cetteexpression utilisedansla nouvelleformulation est ..1 La formulationanalytique la mthode par de volumtrique obtenue I'expansionde Taylor est au voisinagedes entailles.La valeur de la distanceeffective est toujours trs petite et I'utilisationde I'expansion Taylorestdoncvalide. de propose pour I'approche volumtrique.

105

CHAPITRE 4

Chapitre 4
4. Bibliographie sur la fatigue desjoints soudspar point
4.0 Introduction gnrale
Les joints souds par point sont souvent utiliss dans I'industrie automobile. Le par tudipar la mcanique linaire comportement joints souds point est gnralement des par de la ruptureassocie la propagation fatigue.L'tude desjoints souds point par la en de la du mthode la mcanique de linairelastique la rupturedemande connaissance concept defacteurs d'intensitde contraintes. DansI'analysede Lesfacteurs sont d'intensitde contraintes dfinispour l'tudedesfissures. la distribution des contraintesautour d'un dfaut, trois modes de rupture peuvent tre distingus les modes II et III. : I, relatif desdeuxsurfaces la fissurequi setraduitpar de I-e Mode I correspond un cartement une ouvertureangulaire. Mode II correspond un dplacement I-e dans le mmeplan des dansune directionperpendiculaire front de la au deuxsurfaces la fissure,en sensoppos, de dan fissure.I Mode III correspond un dplacement le mmeplan des deux surfaces la de au Dans fissure, sensoppos, en maiscettefois dansunedirectionparallle front de la fissure. la Figure4-1, les trois modes rupturesontprsents. de

modc lll

Figurezt-l Schma trois modesde rupture. des

r07

4.0

Chnpite 4. Bibliopraphie sur la fatieue desjoints soudspar point

globalepour dterminer conditions rupture.En revanche, de les Griffith a utilisuneapproche en locale partir du champdescontraintes fond d'une Irwin l77l a dvelopp une approche fissure.Dans ce critre, la rupture se produit lorsquele facteur d'intensit de contraintes atteint une valeur critique. Invin propose de dcrire la distribution des contraintesau voisinage fond d'une fissurepar la relation: du

ou=hFti@)

(4-r)

I facteur d'intensit de contraintes reprsenteI'amplitude de la distribution des contraintesen mode de rupture I. I-es paramtresr et 0 dfinis dans la Figure 4-2 polaires la zonetudie fond de fissure. en de correspondent coordonnes aux

en Figure4-2 Schma d'unefissureet lescoordonnes fond d'une fissure.

infinies la pointe de la fissure pour r + 0 . kwin montreque les contraintes deviennent due d'un matriaucontinu et Ceci constitueune singularitmathmatique aux hypothses d'une fissure infiniment mince. Le coefficient angulaireFr(O) et le facteur d'intensitde L'expansiondu facteur contraintes dterminent K, I'intensitdu champlocal descontraintes. d'intensit contraintes de s'crit:

*,= r"(#)",o K, , ,"1;

( .\

(4-2)

le , os et a sont respectivement facteur d'intensit de contraintes,

gomtrique, contrainte globaleet la longueurde la fissure. facteurde correction la 108

4.0

Chaoitre 4. Bibliographie sur la fatisue desjoints soudsoar point

pour une gnrale facteurd'intensitde contraintes du L'quation(4-2) est une expression plaqueayant une largeurW, qui contientune fissurea soumise une tension o, avecla
/\

correction gomtri F"l +l . de " "\w/

La mcanique lastiquelinaire de la rupture peut tre utilise pour la propagationdes priodiques. des L'tude de la propagation fissures dansle casd'applicationde chargements cycliques est ralise partir de lois de fissures par I'application des chargements propagation.

109

4.L Loi de propagationde fissureen fatigue

La propagation fissuressousun chargement fatigueest d'une importanceessentielle des de dansl'tude de la fatigue.Diffrentesmthodes sont proposes les auteurs[78-81]. En par pratique, propagation fissures dfiniepar unerelationentrela vitesse propagation la des est de et le chargement. Dans la Figure 4-3, un schmagnralde propagationdes fissuresest montr.
d4

dN

(mml cycle)

= l"',. ".

M,n

LK WPaJ;) Figrue 4-3 Schma typique la propagation fissures desmtaux. de des pour Dans la Figure 4-3, trois zonespeuvent tre distingues: o Zone I est trs proche de la rgion d'une rgion seuil.

o ZoneII estla zoneintermdiaire. vitesse fissuration La de bi-logarithmique ffi ""

'do Iinairedans diagramm"(fr un "

^rK) .

o zone III, la vitesse fissuration jusqu' la rupture. de augmente rapidement L'expression mathmatique ceszoness'exprimepar : de

[o +N l]r t*1. = d
L*

LK < LK* LK,h<LK<K"

LK>K"

(4-3)

Paris une pourlazone cetterelation II. t78l a propos relation s'crit:

il0

4.1

Chapitre 4. Biblioeraphie sur Ia fatieue desjoints soudsoar point

dN

tu = (LK\^ C

(4-4)

o c, N , C, AK et rn sont respectivement profondeurde la fissure,le nombre de cycles, la une constante,I'amplitude du facteur d'intensit de contrainteset une constante. La relation de Paris est uniquement valable dans la zone II. Cette relation est modifie, pour prendreen considrationla zonem [82], par la relation de Forman qui s'crit:

dN

tu =

c (LK)^
(l - R)K" - AK

(4-s)

L'expressionde fissuration(4-5) est utilise gnralement dans les zones tr et ltr. Cette expressionne donne pas une bonne valuationde la vitesse.Selon une hypothsedu dplacement d'ouverture rsiduelle au voisinage de la fissure, Donahue I83l propose I'expression suivante :
daA dN 8",

(tx' - LK:,)

(4-6)

o E, a" et A sont respectivementle module de Young, la limite d'lasticit et une constante qui dpendde la courbe. La propagation des fissures courtes n'obit pas la mcanique de la rupture [84]. Dans la Figure 4-4, les types de fissures longues et courtes sont prsents.

Figure4-4 Schma typiquede la propagation fissures des dansle casdesfissures longues coutes [84]. et

o a, N, AK, m, Leoet n' sontrespectivement longueur la fissure, nombredes la de le cycles,le facteurd'intensitde contraintes, qui I'exposant dpend matriau, dformation du la plastique I'exposant et d'crouissage cyclique.
lll

4.2 Structure desjoints souds par point

Les joints soudspar point sont gnralement utiliss pour assemblerdes tles minces. La mthode du soudagedesjoints soudsest diffrente par rapport d'autres types de soudages communs.Dans le cas d'applicationdu soudage point, l'lectroden'est pas consommable. par Par passaged'un courant lectriqueentre deux lectrodes, soudagepar point est ralis le selon trois procduresdiffrentes : 1. Accostage 2. Soudage 3. Forgeage Dans le premier stade,les lectrodes sont rapproches I'une de I'autre. Les deux tles minces sont situes entre les lectrodes.I-e soudagepar point est ralis par un appareil manuel ou par machines automatiques. Dans I'industrie de fabrication des automobiles,les robots qui quipent les appareils de soudagespar point rduisent le travail. Le soudageest fait aprs la premire phase. Pendant le passagedu courant, une force est applique pour garder les tles minces dans les conditions prvueset obtenir une bonne qualit du soudage.Au cours du forgeage, le courant est arrt et la pression sur les lectrodes diminue. Ce cas amliore la qualit du soudage.Dans la Figure 4-5, les diffrents stades de soudagepar point sont prsents.

ge Accostage

Soudage

Forgeage Forl

Figure 4-5 Schmade la procdrue du soudagepar point.

r12

4.2

Chapitre 4. Bibliosraohie sur Ia fatisue desjoints soudspar point

La procdure soudage trs importante pour une rsistance du est optimaledesjoints souds par point.Dansla Figure4-6, un schma point soud de sazonefondueestprsent. du et
ZAT

ZONEFONDUE Figure Schma soudage pointavec zone 4-6 par du la fondue, zoneaffentn la thermiquement et I'interface.

La zone fondue, la zone affecte thermiquement et le mtal de base sont reprsents respectivement.L'utilisation du soudagepar point donne des zones fondues rsistantesqui sont gnralement deux fois plus rsistantes que le mtal de base.Dans la modlisationdes joints soudspar point, la soudure est considrecomme une partie rigide. Cette hypothse simplifie I'analyse structurelledesjoints soudspar point. Les joints soudsen cisaillement (Tensile Shear) seront considrsdans cette tude. Ce type de joints soudsest frquemment utilis pour attacher des tles minces soumises un chargement axial. L'valuation de comportementdesjoints soudsen cisaillement est considreselon diffrentes hypothses. La plus simple est I'utilisation du conceptde bande [85]. Dans cette mthode,les effets de cisaillementdes bords de la bande slectionne sont ngligs.Le joint soud par point est considrcomme une poutre (Figure 4-7). Selon la thorie des poutre en flexion, I'application de la force rial (F) dplacementdans le modle propos: d o w - ' " d ' z w --U * *d*r * (4-7) conduit au

@ \E

(4-8)

k, w et t sont respectivement force axiale,le modulede Young, le moment la

ll3

4.2

nar sur Chapite 4. Biblioeraohie lafatipue desioints souds point

d'inertie (I =t3 ll2

:urrrit largeur), le paramtre canctnstique, le dplacementdu joint de

soudet l'paisseurde la tle mince. Is conditions aux limites s'expriment par:

=0,+= x=o:w ' dx,' -+ EI

x = I : w = g . dd x = O 't

Ml,=o o Mo=

(4-e)
(4-10)

(4-9) et (4-10) reprsentent conditionsaux limites pour le modledu les Les expressions (4-7), (4-9) et (4-10)donnela par La soudage point en cisaillement. solutiondesquations rotation(x = 0 ) qui s'crit: ( 4- 1 1) (4-12)
La contrainte maximal o** est localise au voisinage du point soud. Cette contrainte est

responsablede la fissuration des joints soudspar point. La contrainte maximale peut tre

par : obtenue la relationsuivante F(. , 6'**=;l*-9 d ) 4t) ok, (4-13)

maximale, force,l'paisseur la la de F , t , d et / sontrespectivement contrainte la

par du tle mince,le diamtredu soudage point et demi-longueur joint soudpar poirrten cisaillement.

4,

...'tO*,

par Figure4-7 Schma modlesimpledejoint soud point assimil unepoutre. du

tt4

4.2

par 4. Ch.apitre Bibliopraohiesur laatieue desjoints souds point

Dans la Figure 4-7 , la localisation de la contrainte maximale est donnepour .r = d l2 . L'expression (4-13) montre que la contrainte maximale est un peu moins que quatre fois la globale applique(ou=l' contrainte uo bt"t ,o"=L). Il est ncessaire mentionnerque de

dans la mthode de bande, une bande de largeur unit est slectionne. Le, facteur de correction f s'crit:

= o!**+(, #)= +o,(r-L)''=('-*), =f.(+o,) o'**

= 1^ 1-9d r6I

(4-r4)
(4-15)

La contrainte maximale du joint soud par point est approximativement gale 4 o r . L e facteur f est voisin de I'unit. La contrainte maximale du joint soudpar point en cisaillement peut tre obtenu selon Radaj [85] :

(4-16)

-u,'"[4.=r) o.",[,1.=, ) ok=

Ry, Mo, M, moment de flexion pour la localisation de contrainte maximale.

(4-r7)

la I, F , t et M sont respectivement force verticale,le moment droite, le

moment x = 0, la demi longueurde poutre,la force axiale,l'paisseurde la tle mince et le

La valeur de la contrainte maximale est donc dfinie par les quations (4-14) et (4-17). La

(4-18)

La mthode la bandeest doncbase un simplemodled'valuation comportement de sur du par plus exactsont aussit proposs mcanique joints souds point. D'autresmodles des par diffrentsauteurs[87-93].Par exemple, dansla Figure4-8, le modlede combinaison d'lments solideset d'lments rigidesest montr.SelonI'existence d'une granderotation

115

4.2

Chapitre 4. Bibliosraphie sur lafatioue desjoints soudspar ooint

dansle modle. rigidessontutilises dansle point soud, rotuleset despoutres des

Pin Joint

Figure 4-8 Schma du modle simple du joint soud par point avec I'utilisation d'lments poutre de plaque et de < pin joint >.

Ls modles prsentssont tous des modles bi-dimensionnels. L'valuation de la distributiondescontraintes voisinage joint souda besoinde modlesplus dtaills. au du pour obtenir une L'utilisation des lments solidesest ncessaire coquesou d'lments possible. distributiondescontraintes plus exacte la propos Salviniet al. estprsent par Dansla Figure4-9,|e modle thorique Dans [94-97]. ce modle,le point soudest considr commerigide. La thoriedesplaques coques et est appliquepour obtenir la distributiondes contraintes autourdu point soud.Une solution qui analyiqueestpropose donnela fonctionde contrainte commeci-dessous : f t I lA: I .= 0) QQ,0) = ro?) + 40 * l fr(r) + lL r 0 lsind + l rr'Q) + 5 r 0 lcosg + ) (r) sin(m +
LzJLzJ^-_z@_D)

m=2

Lfl?)"os(mo)

(4-19): L'utilisationdesconditions limitesde ce problme aux simplifiel'quation Q?,0) = Lrgsine +(nrrt + crr-t + Drrlogr)cosl (4-20)

t16

4.2

Chanitre 4- Bibliopranhie sur Ia fatipue des ioints souds nar noint

4, B, C, Dr sont des constantesqui dpendentdes conditions aux limites du joint soud par point.

(4-20)danslesquations (2-3) (2-5)donne L'utilisation l'quation ( de l'quation 4 - 2 r ) :

o, =l(4+ D,)r-t +2Brr-2Crr4lcos7

o, =168r,+2Crr-3 Drr4lcosl +

(4-2r) (4-22)
(4-23)

-2crr-3+ er-rJsind r,, =p,Brr

Figure 4-9 Schmadu modle pour un joint soudpar I'utilisation de la thorie des plaques et coques [94].

Is dformations et les dplacements peuvent tre calculs partir des relations obtenuesde la loi de Hooke.

Figure4-10 Schma modlepropos Kan [98] a) Schma par du gnralb) Modlesimplifi c) Modlebidimensionnel.

tt7

4.2

par Chapitre4. Biblioeraohiesur la.fatiguedesjoints souds ooint

I modle, montr dans la Figure 4-9, est un modle tri-dimensionnel.Dans la Figure 4-11, les modles tri-dimensionnels sont proposs par I'utilisation des lments coques. Les mthodes prsentesdonnent le comportement global mcanique des joints souds. Il est fortement ncessaire d'utiliser les mthodes tri-dimensionnelles avec des lments

compltementsolides. Cette utilisation permet d'obtenir une distribution des contraintesplus exacteautour du point soud.
Joint SoudPar Point cisaillement Elmentcooue

b)

lil
c)

Figure 4-11 Schma du modle pour un joint soud avec I'utilisation de la thorie des plaques et des coquesa) Mthode de la combinaison des lmentscoqueset solides [85], b) et c) Mthode avec poutres et coques 191,941.

A partir de l'quilibre statique, relationsapproximatives proposes des sont [99]. Danstoutes les expressions pour les joints souds point en cisaillement, contrainte obtenues par la en tractionet la contrainte flexion sont importantes. contrainte en La radiales'exprimepar :

1r8

4.2

Chapitre 4. Biblioeraphie sur hfafieue desjoints soudspar point

o, = o F,cos(p)+ o rr sin() - 6 r, - a u,sin(p) t, a wcos(p) avec : oo= F-

(4-24)

(4-2s)
(4-26)

F. orr=

Fr 10

(4-27)

4 >0
oMa

(4-28)

or, Fr, F y, Fr, M r,M y, t , oFr, orr, oFr, ou, ,aut et A

(4-2e)
les sont respectivement

du de composantes la force, les composantes moment, l'paisseurdes tles minces, les sur bases I'applicationdes forceset des moments composantes tenseurdes contraintes du du angulaire point soud. et correspondant la coordonne I-es autourdu point soud. en La contrainte tractionet la contrainte flexion sontcombines en et contraintes cisaillementont t tudispar Radaj forces,les momentset les diffrentes des et est des [99].Dansla Figure4-l2,le diagramme moments forces montret les analyses des et sur sontbases l'quilibrestatique la compatibilit dformations forceset desmoments par I-es relationsproposes Radaj [99] sont relativementexacteset n'est pas considre. conservatifs. donnent rsultats des : pour lesjoints souds divissen trois catgories proposs sont Les modles approximatifs analytiques 1. Les modles ; bi-dimensionnels numriques 2. Ismodles ; tri-dimensionnels numriques 3. I-esmodles [89-90,100-107]. ne bi-dimensionnels considrent colnmeles modles approximatifs Les modles analytiques que le comportementgnral des joints par point en cisaillement.Les modles tripar globalet local desjoints souds point. Ils ne le dimensionnels considrent comportement globale jointssouds. des sontpaslimits I'analyse

tt9

4-2

Chaoitre 4. Bibliographie sur Ia fatisue desjoints soudspar point

rs@$Dy
Il, Sj:)Fr -r.*
ssotll$ /*llB'
qfl

4.rRn

'J

il4l"* H: t3J"* [:t

rt6r

F*

|iln-

bcces
ogtt

.,Jt

;1.:*
l11i o19r

rn

I Ir
\ \ \
b./
\ -n/

r'il/

gl {r

Figure 4-12Schma l'quilibre moments desforces de des et autour joint soud point[99]. du par Il est ncessaired'envisager I'effet de I'interface entre deux tles minces dans I'analyse des joints soudspar point en cisaillement.La Figure 4-13 reprsente I'interface desjoints souds par points en cisaillement.

yt =,t' v+o)=orIr"T)
Figure 4-13 Effet de I'interface entre deux tles minces dans desjoints soudspar point en cisaillement.

Dansla Figure4-13, 6 estI'interface entredeuxtlesminces. L'paisseur tlespource des type de joints soudsest trs petite.Le rapport entre la valeur de I'interfaceet l'paisseur n'est lui gnralement petit. L'excentricitdes joints soudslie I'existencede pas I'interface provoque un moment supplmentaire. valeur de ce moment ne peut tre La nglige dansles casdesjoints souds pointencisaillement Mono Point. par et 120

4.2

Chaoitre 4. Bibliopraohie sur la fatipue des ioints soudsnar noint

Le comportement mcanique des joints soudspar point dpenddu type de chargement appliqu de la valeurde la forceapplique et [90].

ZONEAFFECTE TIIER.MIQI,JEMENT

METALDEBASE

a)

ZONEAFFECTE THERMIQIJEMENT

METALDEBASE

b)
Figure4-l4Incalisation desfissures dansdesdiffrents cas a) Fatigue oligocycliqueb) Endurance limite

Dansla Figure 4-l4,la localisation la fissureau voisinage point soud prsente. de du est La fissureestrespectivement localise dansle mtalde baseet en7-one AffecteThermiquement (ZAT) dans le cas de la fatigue oligocycliqueet dans le cas de I'endurancelimite. I phnomne la fatigueoligocyclique peuttrerencontr de dansle casd'applicationdesforces globales ou contraintes leves. Dansce cas,la fissurc au voisinage point soudcritique du s'estformedansle mtalde base. endurance En limite,la forceapplique plus raible.La est fissure au voisinage du point soud se forme dans la ZAT. Des travaux exprimentaux

tzr

4.2

Chaoitre 4. Biblioeraphie sur lafatigue desioints soudspar point

montrentque la zonefondueprsente hautelimite d'lasticit. valeur obtenue une I-a dansla zonefondueest gnralement deuxfois plus levequela limite d'lasticitdansle mtalde base 105]. [100-102, par Selonla procdure soudage point (forgeage) du montredansla Figure 4-5, une limite d'lasticitplus ou moinsimportante obtenue zonefondue.La zonefonduedansle cas est en par du forgeage compresse les lectrodes. est Selonle changement phasedansle procd de du soudage, cela conduit une zone fondue homogne rsistante. point soudest et Le considr commerigide en raisonde cettersistance leve autourde la zonefondue.Cette hypothse diminuele tempsde calculnumrique. Elle reprsente assez bien le comportement par mcanique joints souds point.En raisonde cettersistance des lev,e, zoneplastique une ne seforme pasdansla zonefondue.Selonla localisation la formationdesfissures de autour du point soud,l'analyseest diffrenteet diviseen deux approches:propagation avec la mcanique la ruptureet amorage fatigueavecunemthode de en approprie.

r22

4.3 Facteur d'intensit de contraintespour les joints soudspar point

dansle cas des fissuresse linairede la ruptureest applique de La mthode la mcanique par des d'intensitde contraintes joints souds propageant I-esfacteurs dansle mtalde base. Les par point en cisaillement ttudis diffrentsauteurs ont U08-1131. solutionsproposes numrique : analytique par I'approche soit parcesauteurs tobtenues par I'approche ont soit o Chang-Muki [108]:

K,=cqlrrr(n)""]
*,,=, r . o.tt(!)"') [T[r

(4-30)

(4-31)

en le K, Kn , , d et t sontrespectivement facteurd'intensitde contraintes mode de I et II, le cisaillementnominal,le diamtredu point soudet l'paisseur la tle mince. o L. P. Pook[109]:

*,={Gf'*-(n)""]
K= ns s(4)"' a o aLlo, * o.+s ]

(4-32)

(4-33)

en le K, Kr, F , d et t sontrespectivement facteurd'intensitde contraintes mode et de du I et II, la forceaxiale,le diamtre point soud l'paisseur la tle mince. o J. F. Cooper-R. SmithUl0l A. (4-34).

Kr=

(4-3s)

W4
123

4.3

Chapitre 4. Bibliooraphie sur laTatipue des ioints soudsoar point

Km=

,,[L 16, -la\ I

/-\

(4-36)

z"1

K,, K,r, K,, , v , p, u' , v' , r4i sont respectivement facteur d'intensit de le contraintes modeI, II et III, le coefficientde Poisson, modulede cisaillement en le et les dplacements les directions y et z. dans x,

Figure4-15 Schma dtaillde la fissuredansI'utilisationde la mthode Elments des Finis.

Les dplacements , { et vi u'

sont calculs partir des dplacement nodaux.Is

relationssuivantes expriment dplacements ces :

- u,)"ot - @,- r,) rinrpf u,- -l(u,

(4-37) (4-38)

-'') ,' =lu'

2
w=

(u,- u,)sin - (wi- *,)"ore) e

(4-3e)

tti, ttj, v., vj, wi, wjet a sont respectivement dplacements les nodauxdansles directions y, z, et I'angleautourdu point soud. x, Dansla Figure4-16,le schma dtaill fissures des autour point soud prsent. du est Selonl'tudede Cooper-Smith, I'anglede bifurcationde la fissures'exprime par:

r24

4.3

Chapitre 4. Biblioeraphie sur Ia fatieue desjoints soudspar point

(4-40)

Figure 4-16 Schma dtaill de la fissure autour du point soud et I'angle de la fissure partir de la formule par Cooper-Smith[110]. propose

Les calculs obtenus par la mthode des Elments Finis sont prsentssous forme non
3 ( s\.t

(-r dimensionnelle

t^Kt[; - I l- zA\.

Cooper et Smith ont utilis des mthodes numriques pour obtenir les facteurs d' intensit de contraintes. o Zhan9-Radaj[111]

K, =\o,,-o**at,-ouNt

J?,

\/-

(4-4r)
(4-42)

Kn=I@,,-o,,Vt

--f{r,, -r,,)J7 Ko,

(4-43)

la oui, oti , ouo oto, ui , qi et t sontrespectivement contrainte I'intrieurde la tle , la suprieure, contrainte I'intrieurde la tle infrieure, contrainte I'extrieurde la la tle suprieure, contrainte I'extrieurde la tle infrieure,le cisaillement la le I'intrieurede la tle suprieure, cisaillement I'intrieur de la tle infrieureet l'paisseur la tlemince. de

r25

4'3

Chaoitre 4. Biblioeraphie sur lafatisue desjoints soudspar point

Zhangpropose corlme solution:

,

- (t-y1t l+(d, * ofi - o,,o,,)t (ouo or)' -2(o,,ouo ouioro ouooti o,,o,)] * + + * L'\"ai
488

(4-44) J , E , v , t sontrespectivement I'intgraleJ, le modulede Young, le coefficientde Poisson l'paisseur la tlemince. et de Dansla Figure4-17 schma (4-41),(4-42),(4-43) utilispourobtenir quations les ,le et (4-M) estmontr.
TOLE SUPERIEURE a

ZONE FONDUE
\l

__
f

-/
I a

i
KI'Kil'KilI

TLE INFERIEURE

Figure 4-17 Contraintes radialeset cisaillements circonfrentiels autourdu point soudbasselonl'tude de ZhangRadaj[111].

Les quationsde (4-41) (4-44) sont obtenues approcheanalytique.Zhanget par Radaj ont nglig I'effet de I'interface et de la zone affect thermiquement. Is relationsmentionnes restreintes joints souds point en cisaillement. sont par aux Zheng$rll K, = {3F' Zrrdlt

(4-4s)
(4-46)

Ktr=#

126

4'3

Chapitre 4. Bibliographie sur laTatieue des ioints soudspar noint

(447)

(4-48)
F , d , t, E , v, J sont respectivement force axiale, le diamtre du point soud, la l'paisseurde la tle mince, le module de Young, le coefficientde Poissonet la valeur de l'intgraleJ. o Swellam-Kurath-Lawrencelll2-1I4, 116l

Swellam, Kurath et Lawrence ont propos un paramtre de dimensionnementK,. Lawrence [15] relation: a tudi le facteur de concentration de contraintes et propose la

k,=l+ a

(4-4e)

k,, a , t, et p sont respectivementle facteur de concentration de contraintes, une qui constante dpenddu type de joints souds, l'paisseurde la tle mince et le rayon d'entaille au voisinagedu point soud.La constanted est gale 0.43 pour les joints par souds point en cisaillement. La relation (4-49) est difficile utiliser, car en pratique le rayon d'entaille en fond du point soudn'est pas prcisment dfini. L'utilisation de I'expressionprsente dans l'quation (4-49) est en pratiqueimpossible.

'< ..illlliiiiiiiiiiiiiiiiiiiiiiiiiiiiiiliiil|i

Figure 4-18 Cas gnral I'applicationdes actionsautour du point soudselon Swallem-Kurathde Lawrence[112].

127

4.3

Chapitre 4. Bibliosraphie sur laTatieue desjoints soudsoar point

Dans la Figure 4-18, le schma gnraldes forces et moments autour du point soud est prsent.I paramtre K, est calcul selon la mcanique linaire de la rupture et la combinaison des modes I et [I. Les ffois facteurs d'intensit de contraintessont donc dfinis:

(4-s0)

K^o"n=h
Krhror=h

(4-51)

(4-s2)

sontrespectivementlacomposanteverticalede P, Q, M, r, K*rot, K^.nn et K,h"o, la raction du point soud, la composantehorizontale de la raction du point soud,le moment du point soud, le rayon du point soud, le facteur d'intensit de contraintes bas sur I'effet axial, le facteur d'intensit de contraintes bas sur I'effet de moment dans le point soud et le facteur d'intensit de contraintes bas sur I'effet du cisaillement dans le point soud. Selon la formulation de Swallem-Kurath-Lawrence [12], les facteursd'intensit de contraintes s' expriment par:

K, K, = K*at* K^o^"n--14-5E"(4-5\ =+*+ 2 nJ m ' 2r ' Jt r

(4-53)

Krr=Krh"o,-S

K,,=#

(4-s4)

[ facteur d'intensit de contraintesquivalent est dfini par cette mthodecomme cidessous:

K"q =

Kl +B,Kl,

(4-s5)

K"q, K,,

K,, et , sont respectivementle facteur d'intensit de contraintes

quivalent, le facteur d'intensit de contraintes en mode I, le facteur d'intensit de contraintes en mode II et une constantequi dpenddu matriau. Le facteur d'intensit de contraintes quivalent (4-55) ne contient pas les effets de la o** gomtrie et le rapport de contrainte (R ). I-e facteur de correction gomtrique
O.^

r28

4.3

par jointssouds point 4. Chapitre Biblioeraphie bfatieuedes sur doit donc tre dfini. Ce facteur de correction gomtrique considreles effets par: mentionns G9Wt4 Wtz +4rs 13

(4-s6)

W , t et r sont respectivementla largeur du joint soud, l'paisseur de la tle mince et le rayon du point soud. Selon de Swallem-Kurath-Lawrence,le paramtre K, maximum s'crit :

K Le paramtre estdfini par : K,

(4-s7)

-.'* K,=f@)c
I
K,, KT* , Ki" , Ktr , G , f*

.,lu-o'*

(4-58)

R, bo sontrespectivement paramtre Swallemle de

Kurath-Lawrence, valeurmaximaledu paramtre Swallem-Kurath-Lawrence, de la le facteur d'intensit de contraintesmaximal en mode I, le facteur d'intensit de contraintes maximal en mode II, le facteurde correctiongomtrique, constante une qui dpenddu matriau,le rapportde contrainteentre la contraintemaximaleet la contrainte minimaleapplique un coefficientqui peuttreobtenupar essaiet erreur et dansle diagramme K,- N t . I-e coefficient bo est obtenupar lissage la courbereprsentant paramtreK, en de le fonction du nombrede cycleN, . Cettecourbeest obtenue exprimentalement. Dans un diagramme on biJogarithmique, obtientla relation:

x, =e(tt ,Y
diagramme logarithmique - N t . K, bi

(4-5e)

qui A, h sont descoefficients peuvent tre calculspar lissagede la courbedansle

Dans la Figure 4-l9,la procdure d'valuationdu coefficient Doest psente.Il est ncessaire biJogarithmique K,- N t pour diffrentes d'obtenirle diagramme valeurs du rapport de contrainte(R-d'io ). Is valeurs obtenues partir de rsultats o*, exprimentaux peuvent tre reprsentes dans le diagramme bi-

r29

4.3

Chaoitre 4. Bibliosraphie sur la fatipue desjoints soudsoar ooint

. du logarithmique - (1- R) . I lissagede la courbedonnela valeurbo L'expression K, : diagramme -(1-R) estdonne K, I'os(K,) = Ing(K) + boIng(r- R) ou

<o o

(4-60)

LoS(K,)

\ = gtN ....., K, In d4 A)+eh I-o ) = Inde A)+ehlr,glNr) ......tkK,

f,= ft= ft= ft=

.t

hK, = I'od&'q)+\h I-odN r)

& = h dN ..., K, Inde A)+& I^o r)

a)

bs(
LoS(K,)

f)

Los(K,) - I-os(K ) + boLog(1- R)

*,K,
*rK
*tK, *oK

1-& b)

l- R,

t-&

1-&

Ins Q- R)

Figure 4-19 Procdured'valuationde la valeur de bo selon la mthodeproposede Swallem-KurathLawrence a) Schma rsultats bi-logarithmique - N, K, des exprimentaux dansle diagramme desrsultats obtenus dansle diagramme logarithmique - (1- R) . r(, b) Utilisation

1 30

4.3

Chapitre 4. Bibliooraphie sur la fatieue desjoints soudspar point

dpend donc des rsultats La formulation propose par Swallem-Kurath-Lawrence analytique Une approche d'intensitde contraintes. et exprimentaux desvaleursdesfacteurs estutilisepour obtenirla formulationde cettemthode. . Sheppard Ning Pan-Sherry [106]

K, = o^J

e)-o.o* +o.ooo8f 3r.-:(e\' I t a^\t

) , , ,o_);,

!) 3)- o.r4 s( g\- o.oorrg( +o.os23fu)'l

t) d.\f ) o^\t ) \t/ |

(4-62) en le K, K rr, 6^, o6, a et t sontrespectivement facteurd'intensitde contraintes modede rupture I, le facteur d'intensit de contraintesen mode de rupture II, la contraintede de de de membrane, contrainte flexion,la profondeur la fissureet l'paisseur la tle mince. la par de et de La contrainte membrane la contrainte flexion s'expriment :
F o^=a

(4-63) (4-64)

6M ou=w

la F , M et p, sontrespectivement force axiale,le momentdu point soudet le coefficient du de gomtrie joint soud. par sont et de I-esrelations Ning-Pan SherrySheppard obtenues descalculsnumriques.

131

4.4 Tenue ltenduranceen fatigue despoints souds

L'tude de la propagation fatiguedespointssouds trestreinte la mcanique en a linaire de la rupture.La formulationde I'analyseen terme de mcanique linaire de la ruptureest plus avance que I'analyse en tenue d'enduranceen fatigue. L'analyse de la mcanique linairede la rupturedesjoints souds point n'est plus valide dansles casde I'existence par du phnomne d'endurance limit. Selonla valeur du chargement, lieu d'amorage la le de fissure autour du point soud et le mcanismede I'endommagement sont diffrents. L'amorage fissures localisdanslazone affecte des est (ZAT). thermiquement . .ll7-1181

Sheppard

a propos que le domaine de la contrainte structurelleau

voisinage point soud du dpend la vie en fatiguedesjoints souds point. C'estde par -direque les fissuresqui se propagent autourdespoints souds dansla directionde l'paisseursont ngligeables. contraintestructurelleest calcule en condition La lastique.La contraintestructurellelastiqueest la combinaisondes contraintes de flexion, contraintede membraneet contraintesaffectespar la force axiale. I contrainte structurelle mentionne s'crit ci-dessous :
coltraiDt mbrane ooDtraiDte enflexion contraine p8r h fqce axiah

AJ'=

:LO

to)

6Ltr['

fw

+-

AP t'

(4-6s)

LP , LQ, 1, ), LIil' et W sont respectivement composante la verticale de la raction(Figure4-18),la composante horizontale la raction,l'paisseur la tle de de mince,la largeureffective ,=4, soud. Les valeursdesforceset desmoments par sont obtenues le calcul desElments Finis avec des lments poutreset des lments plaques(plate elements). valeurde la La contraintestructurelle propose Sheppard par dpend l'paisseur la tle mince. de de Danscetteapproche, durede vie de la propagation fissures gale la dure la des est de vie totaleen fatiguedesjoints souds point.Une relationestobtenue lissage par par de la courbedansun diagramme (N.* bi-logarithmique -ft, t (4-66) 3' le momentdu point soudet la largeurdu joint

#=c(^s".)*
t32

4.4

par sur souds point Chapite4. Bibliosraphie hfatieuedes.ioints la N/, R, C, S*, et D sont respectivement durede vie totale y comprisla propagation le des fissures,le rapport de contraintes, coefficientqui dpendde la courbe de lissage,la contraintestructurellemaximale et un autre coefficient qui dpend aussi de la courbe de lissage dans le diagramme bi-logarithmique S** N. l-R

I-es Dans la Figure 4-20, un schmadu diagrammebiJogarithmiqueest prsent. coefficientsC et b sontdfinis partirde cettecourbe.
Izg(S*.)

Iog(AS*=rr <r- I-c (#) ) > ,

br

Figure 4-20 Digramme bilogarithmique ^Sr., -+

l-l(

pour valuer la dure de vie en fatigue par

propose Sheppard. par I'approche

Nin Pan-Sheppard de en [05] a utilis I'approche I'endurance fatiguepour lesjoints soudpar I'application du domainede la dformationmorimale de l'quaticn de Morrow [50] et de la relationdeRamberg-Osgood:

,.=$eNy
I

+e'r(2N)"

(4-67)

,"=+.(?)'

(4-68)

I'amplitude eo, dt, om, N, t'r, b, c, oo, E, K' et n'sont respectivement de dformation,le coefficient de la rsistance fatigue, la contraintemoyenne,le en nombrede cycle la rupture,le coefficientde la ductilit en fatigue,I'exposant la de rsistance fatigue,I'exposant la ductiliten fatigue,I'amplitudede contrainte, en de le

133

4.4

par Chapitre4. Bibliooraphiesur hfatieue desjoints souds ooint

module de Young, le coefficient d'crouissage cyclique et I'exposantd'crouissage cyclique. Le domaine de la dformation est obtenu par calcul aux Elments Finis. Ia dformation quivalente est utilis dans cette valuation. L'application de I'endurance fatigue danscettetudemontreque ce conceptpeut tre aussiutilis en dansles conditions fatigueoligocyliqueet d'endurance de limite. Henrensontl00-1021 a tudi les joints soudspar point avec les approches la de fatigue.Danscettetudela mthodeSWT est appliqu pour valuerla durede vie en fatigue desjoints souds. zoneaffecte La thermiquement considr dansles est tudes numriques. Selon les rsultats obtenus par comparaisondes travaux exprimentaux numriques, et I'effet des contraintes rsiduelles limin dansles est premierscyclesd'applicationde la force. La mesuredes contraintes siduellespar les mthodes exprimentalesmontre l'limination des contraintes rsiduelles introduitspar I'oprationdu soudage.

134

4.5Rsum
dansce de sur succincte la fatigueet la mcanique la ruptureest prsente Une bibliographie par chapitre.Les diffrentsmodlespour I'analysedesjoints souds point en cisaillement d'utilisationdu conceptde poutre. la La sontprsents. mthode plus simpleest la mthode par les on Dans cette approche, considre joints souds point en cisaillementcomme une peuvent calculs partir de cettehypothse. tre poutre.Le dplacement la contrainte et permettantde simplifier le problme utilisent diffrentes hypothses D'autres mthodes par I'utlisation deslments par des d'analyse joints souds point en cisaillement, exemple, detype poutreet coque. La mthode de Kan t98l utilise un modle bi-dimensionneldans lequel I'effet de et des dansl'paisseur tles infrieureet suprieure I'interface des concentration contraintes des dansl'paisseur tles des L'amorage fissures les entres deuxtlesmincessontngliges. de dansla mthode Salvini [94-97]. minces n'est pasconsidr Cettemthode par sur proposes Radaj[99] sontbases l'quilibre statique. Les formulations que cette mthodedonne C'est--dire pas la compatibilitdes dplacements. ne considre par Les statiques. conditionsimposes de uniquement rsultats la rsolutiondesquations les danscettetude. lesconditions limitesne sontpasappliques aux est La modlisationnumriquetri-dimensionnelle la mthodela plus exacte et la plus convenablepour considrerles effets des contraintesautour du point de soudureet la dansla directionde l'paisseur distributiondescontraintes Dans les sections4-3 et 44, l'ttrde des joints soudspar point basesur la mcanique linairede la ruptureet la fatigueestprsente. par l'tude des joints souds du de Selonla localisation la fissureau voisinage point soud, en linairede la ruptureou par I'endurance point en cisaillement traitepar la mcanique est fatigue. Dansles casde de de Dansla Figure 4-t4,la localisation I'amorage la fissureestprsente. respectivement dansle limite, la fissureapparat la fatigueoligocycliqueet de I'endurance (ZAT). La localisationde cettefissureest thermiquement mtal de baseet la zone affecte par trs importantedansl'tude desjoints souds point. Si la fissureest localisedansle par est I'amorage ngligeable rapport la est mtalde base(la force applique importante), linairede la rupturepeuttreutilise.Si propagation fissures. Dansce cas,la mcanique des (la est thermiquement force applique faible), la la fissureest localise dansla zoneaffecte

135

4.5

Chapie 4. Bibliograohie sur lafatigue desjoints soudspar point

durede vie en amorageest beaucoup plus importante que celle en propagation.Dans ce cas, la prvision de la dure de vie en fatigue partir de la courbe de Whler peut tre applique. Dansla section 4.4,les approches bases la fatiguesontmontres. sur Shappard a [117-118] que la propagation fissuresest ngligeable. considr des Elle a proposque l'amplitude de contrainte maximaledfiniedansl'quation(4-65)estcellereprsente le diagramme dans de Whler ^,S-,. Nin Pan [05] et Henrenson [100-102]ont utilis I'approchede dformation locale patir desmthodes Morrow [50-55] et S\ryT[56]. L'utilisaton de cesmthodes de est restreinte. Les coefficientsde rductionen fatigue au voisinagedu fond d'entaille sont pour obtenir la durede vie en fatigue.L'tude desjoints souds point en ncessaires par cisaillement n'est actuellement assez pas tabliepour une estimationexactede leur durede vie en fatigue.

136

CHAPITRE 5

t37

Chapitre 5
par point 5. Etude exprimentale joints souds des
5.0 Matriau
Le, matiau utilis est I'acier Solphor280. Cet acier est classdansla catgoriedes aciers galvaniss double face. I-eur capacit protgerI'acier de la corrosion fait des tles par galvanises matriaude choix pour tous les secteurs I'industrie.[-es tles galvanises un de sont utilisesdans de nombreuses applicationsindustriellestouchantdes domainesaussi varisque I'automobile, btiment, travauxpublics,l'lectromnager... Les tles le les etc. d'acier galvanissont obtenues un procdqui combineen une seule tape,le recuit par continude recristallisation la banded'acier lamine froid et sonimmersiondansun bain de de zinc ou d'alliage de zinc fondu. La banded'acier dfile sansarrt dans les diffrentes zonesde traitementet la sortie du bain, un systme d'essorage lames d'air ou d'azote permetde rglerla quantitde revtement dpose lesdeuxfaces. sur Les revtements base de zinc protgentI'acier contre la corrosion.Ils jouent non seulement rle de barrire le physiqueen isolant I'acier de I'environnement, mais ils scrtent aussi des oxydes qui recouvrent protgent zones I'acier a tmis nu (tranches et les o cisailles, rayures, etc.). Bien que cela paraisse paradoxal, est ncessaire protgerles tles galvanises un il de par traitementsupplmentaire. effet, lorsqu'il est dansune atmosphre En confine,le zinc ne forme plus d'oxyde protecteur mais dveloppe contraireun compos pulvrulentappel au "rouille blanche".L'espaceentredeux spiresd'une bobineo entredeux tles d'un paquet sont autant d'endroitso la "rouille blanche"peut sedvelopper les conditionsde stockage si ne sontpasadaptes. I-es tles galvanisessont donc livres chromatesetlou huiles afin de limiter temporairement risques de "rouille blanche". Un revtementorganiquemince peut les galement appliqu surface, garantissant excellente tre protection en une temporaire contrela corrosion.Les tles galvanises, demandes expressment aucuntraitementde surface, sans sontliwessousla seuleresponsabilit demandeur. du Le Solphor280 est un acierrephospho avecune hautelimite lasticit. Cet acierpeut tre utilis pour I'opration du soudage. Pour viter le phnomne rouille blanche,il est de 1 38

5.0

Chapitre 5. Etude exprtmentale des ioints souds par point

ncessaire prendre toutes les prcautionspour liminer les risques d'humidit et de pendant le stockagedes tles galvanises. Tableau 5-1 reprsente I-e les condensation gnrales caractristiques techniques I'acier Solphor220,Solphor260et Solphor280. de par d'adapteret de recalerrgulirement En soudage rsistance lectrique, est ncessaire il par peut tre amliore ragrage en I'intensitdu courant.La durede vie deslectrodes et utilisant une qualit approprie(alliage Cu-Cr-Zr par exemple).En soudage I'arc, il est ncessaire rduirela vitesse afin de etlou d'utiliser desmtauxd'apportappropris d'viter la formationde porosits duesaux vapeurs zinc. de
T"bt""" 5-1 Crt*trrtq
Eseide lredion m Qualo Eolphor R. {MPE} R{MPs} FosSbilll rvlemenl de el
fp

A%ru mini

mini

n mini 0,19 0,18 0,18

$lsndsrd

Galllex

220 2S tro

22f?69 34ry40 260/310 380/440 ?801340 38tr440

't,7
1.5

r,5

9c ruEssai senstra{ers eur prowette $t- Q,2 enGakallia f90 20 x 80. p - 2 nhen GalrrEllie . = f a i s a b l eN = N o n f a i s a b l e , = nousconsullgr.

Le Solphor 280 a aussi une bonne capacitpour supporter des oprationsde formageet d'assemblage comme d'autres aciersutiliss dans I'industrie automobile,le btiment,les et habituellement travauxpublics,l'lectromnager. techniques formage d'assemblage Is de pour la tle galvanise. Toutefois, utilisespour les tles d'acier non revtu,restentadaptes plus le revtement pais,plus I'aptitude la dformation au soudage et diminue.Il est est que doncncessaire la quantitde revtement demande compatible soit non seulement avec le niveaude protectioncontrela corrosionmais aussiavecles exigences formageet de de soudage.

r39

5.1 Composition chimique et comportement mcanique du matriau desjoints soudspar point

par La composition chimiquedu Solphor280 utilis pour les joints souds point est donne par le Tableau 5-2.

par sontdtermines desessais tractionralise de I-esproprits mcaniques matriau du sur sontprsentes le Tableau dans unemachine type INSTRON.I-esconditions du d'essais 5-3.
I'ableau5-3 Condluon d'essaisDour l'acler

TYPE DE LA MACHINE Frequence(Hz) ('C) Temprature

INSTRON

30
Ambiante

mcaniques, caractristiques les mcaniques prsentes Tableau sont au 5-4. Selonles essais
zEO.

ableau5-4

de I'acier

o, (MPa) Limite d'lasticit 299

ou (MPa) Rsistance ultime

AVo
Allongement

43r

32.2

ultime ou sont obtenues partir I-es valeursde la limite d'lasticito, et de la rsistance desessais tractionsimple. de cycliquesont aussiimportants des cycliques de la loi de comportement ou Les rsultats essais (courbed'crouissage cyclique).Le trac de la courbecyclique permet de dterminerles cycliquesstables matriau. loi de comportement du La cycliquedu matriau caractristiques simultanment dansle et les coefficientsde la loi de Ludwick (oo = K'ei.o) sont prsents Tableau et la Figure5-1. 5-5

IN

5.1

Chapitre 5. Etude exprimentale des ioints souds par point

Tableau Coefficients la loi de Ludwick ou de la loi de 5-5 de ,

rque.

K' (MPa)
576

0.1001

^s

c) .E

()
o. o c) E

c-=-

oo

"

205000 \576)

Tt

,(o"')orool |

E=205 GPa K'=576 MPa n'= 0.1001 v =0.3

Amplitude de dformation ea Figure 5-1 Loi de comportement cyclique (courbede amplitudede conhainte- amplitudede dformation).

La courbede comportement cycliqueestobtenue qui avecuneprouvette a la mmepaisseur que les prouvettes joints souds point ralisavecdestlestrs minces.Il est donc par des ncessaire d'utiliser la mmeconditionde gomtrie avecle mmematriau. Pour viter le problme du flambage de tle mince, deux renforts ont t appliqus.I-es valeurs de I'exposantd'crouissage cyclique et du coefficientd'crouissage cyclique sont obtenues partir de la loi de comportement cyclique(Tableau 5-5).

L4l

5.2 Gomtrie joints souds des par point

Pourla prvisionde la durede vie en fatigue,desprouvettes dejoints souds point avec par diffrentes gomtries utilises. joints souds point sontclassifis sont Is par commesuit : 1. Jointssoud point en cisaillement par Mono Point 2. Jointssoud point en cisaillement par Triple Points 3. Jointssoud point en cisaillement par Multiple points Selonl'utilisationdesjoints souds dansI'industrie automobile I'assemblage pices et des de carrosserie vhicules, tlesmincesavecdiffrentes des des paisseurs peuventtre utilises. En pratique,l'paisseur gnralement est gale.L'analysedesjoints souds point avec par diffrentes paisseurset diffrents comportementsmcaniquesest plus complique. Gnralement, tles mincesavecles mmespaisseurs les sont utilisespour fabriquerles assemblages. Dansla Figure5-2,|a gomtrie joints souds point en cisaillement comprisMono des par y Point,Triple Pointset Multiple Pointsestprsente ll07l.

Figure 5-2 Gomtrie joints souds des par point en cisaillementy comprisMono Point, Triple points et Multiple Points.

I-es tudesdesjoints souds point en cisaillements par pour la prvisionde la vie en fatigue sont restreintesaux joints soudspar point en cisaillementMono Point. L'tude de

t42

5.2

par des Chapitre5. Etudeexprimentale joints souds point

l'influence du nombre de points soudspour la dure de vie en fatigue

i-pote. Pour "rt cette raison, des joints soudsMono Point, Triple Point et Multiple Points sont utiliss pour valuer la dure de vie en fatigue avec diffrents gomtrieset nombres de points souds. est dans la Figure 5-3. L'aspect microscopiquedesprouvettes slectionnes prsent

c)

d)

par Figure5-3 Aspectmicroscopique du X50 b) d'une prouvette soude point a) Grandissement point soud du du X50 d) Grandissement poini soud Grandissement point soudX500 c) Grandissement point soud du x50o [l19].

Dans la Figure 5-3, un grossissementdu point soud pour les cts gauche et droit est prsent. La zone fondue, la zone affecte thermiquement et le mtal de base sont respectivementlocaliss. I-e centre est la zone fondue et la ZAT et le mtal de base peuvent tre distingus autour de la zone fondue. L'interface entre les tles minces est visible. Une observation macroscopique donne la valeur moyenne de la dimension de I'interface. Cette

143

5.2

Chapitre5. Etudeexorimentale ioints souds point par des

valeur sera considredans le calcul par la mthode des Elments Finis. Il est ncessairede savoir que la valeur de I'interface changeles valeurs de la contrainte maximale. Cette prise en compte de I'interface conduit une prvision de la dure de vie en fatigue moindre qu'en son absence. Dans la Figure 5-4, les dimensions du point soud mesures et utilises sont indiques.
ZAT

ZONEFONDUE

Figure 5-4 Schmade la gomtrie du point soudpour la taille du point soud, de I'interface et de l'paisseur.

Dansla Figure5-4,|a dimension l'interface de mesure partir d'un examen microscopique estobtenue. SelonI'utilisationdesprouvettes soudes, deuxplaques avecla mmepaisseur de tle sont utilisespour appliquerla force. I.es tles mentionnes sont montres dansla Figure5-2.

t4

par 5.3 Procduredesoprationsde soudage point

Les joints souds par point sont utiliss selon une procdure spcifique de soudage. L'lectrode utilise dans I'opration de soudage est gnralement consommable. Cette lectrode possde une composition chimique particulire. Cette composition particulire permet d'avoir desjoints soudsavec un haut comportementmcanique.L'lectrode utilise pour lesjoints soudspar point est en cuivre. Dans la Figure 5-5, le schmade l'lectrode des joints soudspar point est prsent.

Electrodesuprieure cuiwe en
Tle

Forcede compression Pointde soudure

Tle infrieure Forcede compression

Elctrode infrieureen cuivre

Figure 5-5 Schmade l'lectrode applique pour souduredes joints soudspar point.

localemententre deux lectrodes Dans la Figure 5-5, les pices soudersont serres en par qui est alliagede cuivre. L'ensemblepices| lectrodes travers un courantde soudage par provoqueune lvationde temprature effet Jouleet la fusion localisedesdeux pices dansla zonede positionnement deuxlectrodes. des la Lorsqu'un courant lectriqueparcourtun conducteur, quantit de chaleur dgage est proportionnelle carrde I'intensitdu courantet au temps.L'expression mathmatique au de
ces paramtres s'exprime par: Wr = R*xI*2 xt*

(s-1)

145

5.3

Chapitre 5. Etude exprtmentaledesjoints soudspar point

Wt , R*, I* et tw sont respectivement l'nergie de soudage, rsistance la lectriquedu conducteur, I'intensitdu courant lectrique le tempsde soudage. et La rsistance d'un conducteur, temprature ambiante est: 1. proportionnelle salongueur 2. inversement proportionnelle sa section 3. fonctiondu matriau constitutifde l'me du conducteur L'expression la rsistance de lectrique d'un conducteur s'critdonc:
D vf

&,=Y ,\, "

Rw , Pw, I* et matriauconstitutif de l'me, la longueur et la section.

G-2)
sont respectivementla rsistancelectrique du conducteur,la fonction du

Il est ncessairede mentionner que diffrents types de mise en pression des lectrodes existent.Ils peuventtre classifis : l. Commandemanuelle

2. Commandepneumatique 3. Commandehydraulique 4. Commandemcanique Chaque commande a exclusivement son utilisation particulire. La commande hydraulique peut tre faite par pressiond'eau ou d'huile. Dans I'utilisation desjoints soudspar point, le diamtre du point est un paramtreimportant. Il est fonction des paramtressuivants: l. de I'intensitdu courant

2. du temps de passage courant du 3. de la section de contact des lectrodes 4. de I'effort appliqu aux lectrodes I diamtre de contact de l'lectrode doit tre environ deux fois l'paisseur de la pice souderplus 3 mm. L'expression de cefte restriction s'crit: d" = 2t +3 (5-3)

d" et t sont respectivementle diamtre de contact entre l'lectrode et les tles (en mm) et l'paisseurde la tle applique(en mm). La distance entre les points de soudureet le bord de la pice est donc dfinie par une norme. La distance entre la rive de la pice et le point de soudure doit tre environ deux fois

r46

5.3

Chapitre5. Etudeexprtmentale joints souds point des oar

l'paisseur du point de soudure plus 4 mm pour viter un talement du point de soudure sous l'effort de la pression hydraulique du noyau en fusion. L'expression de cette rgle s'exprime par: Ln=2t+4

(5-4)

L^ et t sont respectivementla distance entre la rive de la pice et le point de soudure (en mm) et l'paisseurde la tle mince (en mm). La distance minimale entre deux points de soudure est donc dfinie par une rest:iction particulire. La distance entre deux points soudsdoit tre environ trois fois le diamtre du point de soudurepour viter un shuntagedu courant de soudage.Si la distanceentre les points est infrieure trois fois le diamtre du point, il est ncessaired'augmenter I'intensit pour compensercet effet de shuntage. LP>_3d

(5-s)

L, et d sont respectivementla distance entre les deux points de soudure et le diamtre de point de soudure. Dans la Figure 5-6, la distanceminimale entre deux points de soudureet la distanceminimale entre la rive de la pice et le point soudpar point est aussiprsente.

Figure5-6 Schma paramtres et Z" dansle casdesjoints souds des Z^ multi points.

Dans le cas du soudagemulti points, I'excutionsimultane plusieurspoints sans de dplacement pices permetd'amliorerla productivit des surtoutlorsquele nombrede points raliserestimportant.I-e courantparcourt successivement chaque coupled'lectrodes. Dansla Figure5-'7,leschma soudage point,le typedu multi pointsestprsent. du par Ce type de soudageest frquemmentutilis pour souderles tles minces dans I'industrie automobile pour les ca:rosseries vhicules. des

147

5.3

Chapitre 5. Ende exprimentaledes ioints soudspar ooint

La procdure soudage point estdiviseen deuxmthodes par du : 1. Soudage sansvariationde I'effort de compression 2. Soudage avecvariationde I'effort de compression

Figure 5-7 Schmadu soudagepar point, type multi points.

Le soudagesans variation de I'effort est le plus frquemment utilis. Le soudageavec variationde I'effort est utilis dansdes applications particulires. Dans la Figure 5-8, les cyclesde soudage pour les deuxcassontprsents.

z
o fI{

()
i:

Courbede la Force

= O

F_

I
Temps(Sec.)
Forgeage Tempsnort

a)

148

5.3

Chapitre5. Etudeexprtmentale iointssouds ooint des par

z
C)

!?

tr sq U

Courbede la Force

F-

'emps (Sec.)
Tempsmort

b)
Figure 5-8 Cycle de soudage des joints soudspar point a) sansvariation de I'effort de compression b) avec variation de I'effort de compression.

Dansla Figure 5-9, les procdures soudage joints souds de des Mono Point, Triple Point et Multiple Point sontrespectivement prsents.

t2w 1000

- -

Force de compression Courantdesoudage

z
(t) !)

15g
o bo 'o

.9 800 U)

E 600
(.)

rc3
() rJ
q)

; 400 (J
tr.i

!)

2W 0 0 10

O ()

0
RI (h

ru o
C

E
JO

C) 0

20

a)

3040 Priode

60

70

149

5.3

-Force de compression - - Courant de soudage

^ 2500

3
g

.9 zooo (t)
(

M
c) o ct) c)

o c)

F ls00

a.

; looo
sl

F
r\

s00

l0 b)

30 40 Priode

Figure5-9 Cyclede soudage joints souds point pour les prouvettes des par slectionnes Mono Point b) a) Triple Pointset Multiple PointsI 19] (1 Priode= 0.02Sec.). Le cycle de soudage sans variation de l'effort

de compression est utilis pour toutes les

prouvettes soudespar point en cisaillementy compris Mono Point, Triple Points et MultiplesPoints.La procdure soudage considre trois tapes de est de : 1. Accostage 2. Soudage 3. Forgeage Dansla Figure5-9a,la procdure pour lesjoints souds point Mono Point estmontre. par La procdure pour les joints souds points Triples Pointset Multiple Points est identique. par Cetteprocdure prsente la Figure5-9b. est dans
Tableau Rsum la procdure soudage pointssouds 5-6 de de des pourlesprouvettes slectionns.

TYPE
Mono Point Multiple Poin Triple Points

ACCOSTAGE (Priode)

SOI'DAGE (Priode)

FORGEAGE COURANT (Priode) (KA)

FORCE (KI.I)

ELECTRODE

40

9 9 9

9 9 9

7.6 7.6 7.6

1. 08

Hnisphrique /6

40 40

2.52 2.52

Hmisphrique /6
Hmisphrique/6

1 50

5.3

Chapitre 5. Ende ewrimentale des ioints soudspar point

I Tableau reprsente 5-6 brivement procdure soudage pointsde soudure la de des dansles cas Mono Point, Triple Points et Multiple Points.L'lectrodeutilise dans I'oprationde soudage joints souds point est une lectrode par des hmisphrique [119]. I diamtrede cettelectrode de 6 mm. Dans la Figure 5-10, l'lectrode de soudage est hmisphrique et sondiamtre sontschmatiquement prsents.
Tle mince

Electrodes Hmisphriques Cuiwe en Figure 5-10 Schma l'lectrode soudage de de hmisphrique pour jointssouds point. applique des par

151

5.4Prsentation essais fatigues des de

Le principe de I'essai consiste appliquerune force axiale sur une prouvettesoude par point.I-e type de chargement obtenuest sinusoidal amplitude constante. machine La utilise estmontre dans Figure5-11. la

g
r1) \c)

(n
0
k It)

F)
q) .F )!)

Eprouvette
(l)

o
.F

Figure 5-11 Positionnementde l'prouvette soudedans les bridages hydrauliques de la machine << Schenk >.

Les essaissont effectussur machine<< > Schenk du Laboratoirede Fiabilit Mcanique. Dans la Figure 5-11, les vrins hydrauliques la machine<Schenk> appliquentle de chargement. mors infrieur est fixe et le mors suprieurest libre. I schmades Le prouvettes soudespar points est prsentdans la Figure 5-12. Dans cette figure, le positionnementdes prouvetteset la condition d'application de la force axiale sont schmatiquement prsents. I-e phnomne fatiguesediviseen deuxparties: de l. Amorage fissures des 2. Propagation fissures des De petitesfissures peuvent qui treprovoques le procd soudage par de s'amorcement dans la premire partie.I-espetitesfissures peuventsetrouvergnralement autourdu point soud

152

5.4

Chanitre 5- Etude exnrimentale des ioints souds nar noint

car une grande concentration des contraintes est observe au voisinage du point soud infrieur.
Force ,e Force

a a a

Bridage Fix
Figure 5-12 Sch

BridageFix

par ts souds point Mono Point, Triple Points et Mul tiple Foints.

Dansla Figure ns Fi fa(rteur un facteur trs

fissures prsente. nombre de joints souds est I-e est fissurati on.

'

{-

.a

Brid Fi) re Figure 5-13 Sch Multip e Points. tiple Poi desjoins par point Mono Point,
le Points et

15 3

5.4

par Chapitre5. Etudeexprtmentale joints souds point des

[s fissures sont frquemment perpendiculaires la direction d'application de la force. Des fissure d'amorage peuvent tre trouves sur le ct de l'application de la force. Le type de fissuresrencontresest montr dansla Figure 5-14 ainsi que la micrographie correspondante.

I I

N 4* = 2100 N^ = 800000 cycl

a - O.5l mm N = 250000 cyc

a)

I I I

4n* = 2100N N^ = 800000cycles a - 0.61mm Ni = 300000cycles

b)
\--t/

Figure 5-14 Propagation fissuredansun point soudavecforce maximale2100 N et nombrede cycles de appliqus 800000 a) 250000b) 300000.

154

5-4

Chaoitre 5. Etude exprtmentaledesjoints soudspar point

Pourcetteforce,un nombrede cycles la rupturegal 800000a pu treestim.La longueur de la fissureaprs a:rt.250000et 300000cyclesestrespectvement 0.51 mm et 0.61 de mm. La fissureprsente dansla Figure5-14sepropage versla surface suprieure la tle. de L'influence du nombrede points souds importante pour la prvisionde la durede vie est en fatigue. Toutefois, le mode d'endommagement changele comportementgnralde la fissurationdes prouvettes soudes. mode de la fissurationdes joints soudspar point Le Triple Pointset Multiple Pointsrestele mme.I modede fissurationdesjoints souds mono pointsestquant lui diffrent.Dansla Figure5-15,les modes fissurationdesjoints souds de par point slectionns prsents. sont

a)

b)
Figure5-15 Mode d'endommagement prouvettes des par soudes point en cisaillement a) Mono Point b) Triple Pointset Multiple Points.

155

5.4

Chapitre 5. Etude exprimentale des ioints souds par point

I-ocalisationde I'endornnagement d'une prouvettesoude Mono Point

Localisationde I'endommagenpntd'une

Localisationde I'endornrnagernt d'une prouvette soude Multiple Points

a) :.
FFr 5o
50tr

E> lo
Eo) E> t)t
q)

.F
trt t () () r
E v)

E 'o a
G) l) C) o

a tr c)

E
!)

q) E o. J \t)

{)

., r tE Bc
E

b) des Figure 5-16 Prsentation I'endommagement joints soudspar point en cisaillementa) Schmade de I'endommagement Eprouvettes b) endommages.

La frquencede la force sinusoidaleet le rapport de contraintessont I =30H2 et R='=*

F

&*

=0.10 respectivement. pour toutesles prouvettes Cetteconditionest conserve

(MonoPoints,Triple Pointset MultiplesPoints). par soudes point slectionnes

r56

5.5 Rsultat exprimentaux des joints souds par point en cisaillement

dans la Les joints soudspar point sont utiliss avec les conditionsd'essaismentionnes par des dans section 5-3. Dansla Figure5-17,les rsultats joints souds point sontprsents un diagramme Log(F**)- LoS(N^).

10 0 0 0

ocD

E Mono point A Triple points O Multiple points

z
rII

tr
I 000 1,0E+05

Nombre de cycles(Nn)

1,0E+06

1,0E+07

Figure 5-17 Courbe Whler pour les prouvettes des joints soudspar points (Mono Points, Triple Points et Multiple Points).

I-escoefficients la loi de Basquin4"* = F;(N *)u pourla force axialesontdonns de dansle Tableau5-7 sousforme de rsistance fatigueavecle facteur de corrlation en correspondant par R2 auxjoints souds point en cisaillement, Mono Point,Triple Pointset Multiple Points.
Tableau 5-7 Rsultats exprimentaux obtenus pour les joints souds par point en cisaillement (Mono Point, Triple Points et Multiple Points).

I'ype
Mono Potnt , fl=l

Rsistanceabsolue (KN) F-u*=34.I 90(N*;-o'zo+o

Corrlation R2=0.8998 R2=0.96+t R2=0.9696

Inple Polnts , m=3 Multiple Points,m=5

lo+z F."*=300.424N*;-o' tr F-*=3 33.394(N*;-o'r t

157

5.5

Chaoitre 5. Etude exprimentaledesjoints soudsoar ooint

La rsistance prsente dansla Figure 5-17 est la rsistance globaleou la rsistance absolue desjoints souds pointen cisaillement. par Dansla Figure5-18,la rsistance la fatiguedes de joints souds point pour Mono Point, Triple Pointset Multiple Pointsest respectivement par divisepar 1 et 3. Mais, lesjoints souds point Multiple Pointssontplus faiblesque les par
joints soudsTriple Points en valeur relative.

10000
El Mono Point ln, n=l A Triple Points /n, n=3 o Multiple Points ln, n=3

e
x
fr

3-o

tr

1000 1,0E+05

1,0E+06

1,0E+07

Nombre de cycles (Nn)

Figure5-18 Courbede Whlerpour les prouvettes joints souds point (MonoPoints,Triple Pointset des par Multiple Points)avecI'utilisationdu facteur . n

Tableau Rsistance 5-8 absolue relativepour lesjoints souds point en cisaillement et par (Mono Point,Triple Pointset Multiple Points). Type Mono Potnt , ffi=l

Keslstanceabsolue(l(N ) F.ar=34. I 90(Np-o'zole Fna'=300.424(Nn)-o''*'

Kesrstance retatlve (I(N ) Fma*=34.90(N*;-o'zo+e I

Corrlation

R2=0.8998

lnple Pornts , m=3

F."*= I 00. 14I (Nn)-o'3M2 R2=0.9641

MultiplePoints, m=5

F.*=333.394(N*;-o'rt:t

F^ =66.678(N6-o'rtrt

R2=0.9696

I-e Tableau5-8 prsente respectivement rsistance la absolue,relative et le facteurde la corrlation.la rsistance absolueest le rsultat obtenu par les essaisde fatigue, et la rsistance relative est le rapport de la rsistance absolueet le nombre de points souds appliqus dansles joints souds point en cisaillement. par L'expression mathmatique la de rsistance relativeestdoncdfiniepar:

r58

5.5

Chapitre 5. Etude emrimentale desjoints soudsoar point

^F;

AFL,. = ml

uo*poiot

,"f;

-34.19(b!)4'2ge

=34.19(N)4.?Ms

KN

(s-6)

et pour lesjoints souds point Triple Pointset Multiple Points,les expressions par suivantes s'crivent: ^4** =
*F;

o4* m3
nIou"

tnon"o.o ,^F;

-300.424(!'l*)u't*'=100.141(\;-o.ro+z KN
-333.394(N)4'3r3t =66.67g(Np)4.3r3r KN

(s-7)

*F;,

m5 oF**, m et N* sontrespectivement rsistance la relative,la rsistance absolue, le

""oo."oto ,"r;

(5-g)

nombrede pointssouds le nombrede cycles. et La comparaison rsistances des relativesmontreque lesjoints souds point Triple Points par ont unersistance relativesimilaireauxjoints souds point Multiple Points.Malgr,le fait par que lesjoints souds point Multiple Pointsont deuxpointssouds plus, leur rsistance par en relative est moins que la rsistance relativedesjoints souds point Triple Points.Il est par donc videntque deux points souds plus danslesjoints souds point Multiple Points en par ne conduisent unersistance pas supplmentaire. laFigure 5-I7,la comparaison Dans des rsultats montreque lesjoints souds point Triple Pointset Multiple Pointsont la mme par rsistance fatigue. en Dansla Figure5-18,la comparaison indiquequela rsistance joints des par souds point Triple Pointset Multiple Pointsest un peu moins de trois fois plus grande que les joints soudspar point Mono Points. Par exemple, pour le nombre de cycle, N n = lE + 06 , la rsistance absolue est gale respectivement oF;* = 2.O76 , KN ^Fo,"*=4.492 et oF**=4.409KN pour KM Mono Point,Triple Pointset MultiplePoints ll20). Le seul avantagedes joints soudspar point Multiple points est une meilleure rsistance fatigueoligocyclique. en

159

5.6Rsum
L'tude exprimentale joints soudspar point en cisaillementest prsente des dans ce chapitre.Le matriauutilis pour les joints soudspar point en cisaillementest l'acier Solphor280. Ce type d'acier estclass dansla catgorie aciers hautelimite d'lasticit. des Il est galvanis doubleface.La composition sur chimiqueet le comportement mcanique de cetaciersontprsents lesTableaux dans 5-2,5-3,5-4et 5-5 et la Figure5-1. La gomtrie joints souds point en cisaillement montre des par est dansla Figure 5-2.[-es joints soudspar point Mono Point, Triple Points et Multiple Points sont respectivement prsents. Dans le cas desjoints souds point Triple Points,la distanceentreles points par souds gale 24 mm. Dans le cas desjoints soudspar point Multiple Points, cette est distance rduite 12 mm. La distance est entreles points souds trs importante est dansla prvisionde la durede vie en fatigue.I-es rsultats obtenuspar les travauxexprimentaux montrentquelesjoints souds point Triple Pointset Multiple Pointsont desdures vie par de en fatiguesimilaires.Ce rsultatmontrequela distance entredespointsde soudure change la rsistance la fatigue.La distance entrelespointsde soudures normalise, exemple par est les distances et L, sontutilises L* danslesquations et 5-5. 5-4 Dans la Section 5.3, la procduredes oprationsde soudagepar point est prsente. L'lectrode utilise dans I'opration de soudagepar point n'est pas consommable. La procdure soudage compose trois stades: Accostage2) Soudage3) Forgeage. de est de 1) Dansle Tableau5-6, la procdure soudage pointssouds point,Mono Point,Triple de des par Pointset Multiple Pointsestrsume. Dans la section5.4, les essais fatiguesont prsents.I-esjoints souds point Mono de par Point,Triple Pointset Multiple Pointssontsollicitssur unemachine type << de Schenk>. L'influence du nombre de points soudsest observepar la comparaison des modes d'endommagement joints souds point Mono Point, Triple Pointset Multiple Points. par des Selonles observations exprimentales, joints souds point Triple Points et Multiple les par Pointsont les mmesmodesd'endommagement. joints souds point Mono Point ont Is par un mode d'endommagement diffrent @gure 5-15 et Figure 5-16). Dans le cas de I'endommagement joints souds des Mono Point, il y a une granderotation pour les tles mincessuprieure infrieure. et Pourlesjoints souds Triple Pointset Multiple Points, il y a une lgrerotationet lesjoints souds sontcoups deuxmorceaux. en Dansla section5.5, les rsultats exprimentaux joints souds point en cisaillement des par sontprsents. Selonles

160

5.5

Chaoitre 5. Etude exprimentale desjoints souds oar point

par les exprimentaux, joints souds point Triple Pointset Multiple Pointsmontrent rsultats que des (Tableau 5-7).Il estvident la rsistance joints souds absolues rsistances lesmmes en multi-points(Triple Pointset Multiple Points)est plus grande fatiguegrceau nombrede par relativedesjoints souds point en points de soudure. Dans le Tableau5-8, la rsistance par I-es est cisaillement prsente. joints souds point Triple Points sont les plus rsistants. dans le cas des joints soudsMono Point, la rsistance Selon le mode d'endommagement par relative est faible par rapport celle desjoints souds point Triple Points et Multiple par relativedesjoints souds point Multiple Pointsestplusfaible Points.De plus la rsistance relativede ce type dejoints souds par quelesjoints souds point Triple Points.La rsistance et entreles points de soudure le nombre donc de la distance par point en cisaillement dpend de pointssouds. les En pratique,les joints soudspar point Triple Points et Multiple Points reprsentent supplmentaires deux points souds En mcaniques. d'autrestermes, mmescomportements I'effet du nombre en pas de n'engendrent un accroissement la rsistance fatigue.En revanche, plus lgrement grandeque celle desjoints absolue conduit une rsistance de pointssouds Dansle casd'endurance par souds point Triple Pointsdansle casde la fatgue oligocyclique. limite, les joints soudspar point Multiple Points montrentune dure de vie en fatigue moinsleve. sensiblement

l6l

CHAPITRE 6

t62

Chapitre 6
6. Modlisation lastique et lasto-plastiquedes joints souds points par
6.0 Introduction
Pourdterminer distributions contraintes, logiciel ANSYS a t utilis.ANSYS est les des le un code gnral calcul bassur la mthodedesElments de Finis conu pour I'analysedes structures. Il permetde traiter des problmes d'lasticitlinaire,des problmes non linaires(lastovisco-plastique), des problmesdynamiquestransitoireset stationnaires, des problmes thermiques transitoires stationnaires, problmes et des mcaniques fluide (CFD) et des du problmes acoustiques. Poureffectuer calcul,ce logiciel passe les trois stades un par suivants : o Stade : dfinitionde la gomtrie modle: I du 1. Introduction lments des pour convenables le typed'analyse considr 2. Dfinition desdonnes caractrisant modle: le o Iscoefficients matriaux des utilisspour le modle; o Constructionde la gomtrie dans les cas bi-dimensionnelet tri(Is surfaces, volumeset la combinaison surfaces dimensionnel les des et volumes); o Ismaillages modleconsidr du ; o Vrificationdestaillesde maillages obtenues. o Stade : solutiondu problme II dfini l. Applicationdesconditions limites ; aux 2. Applicationdeschargements ; 3. Arangementdestypesd'analyse ; 4. Vrificationcomplte modleavantla solution; du 5. Solutiondu problme o Assemblage matrices rigidit et de masse la structure des de de complte

r63

6.0

Chapitre 6. Modlisation Iastiaueet lasto-plastique desjoints souls par point

o o o

Application des conditions limites Application des chargementsexternes Rsolution du systme d'quation linaire d'quilibre ou des systmes non linaires.

Stade III : analyse et post-traitement des rsultats qui peuvent tre, soit des quantits locales comme les dplacements,les contraintes,les dformations. Soit des quantits globales, telles que l'nergie de dformation.

L'analyse des joints soudspar point a t ralisepar ce programme. Il permet de faire un paramtragepour les diffrents calculs par Elments Finis. L'analyse des Elments Finis par ANSYS peut tre ralisepar trois mthodes[L2L-I22]: l. Interactive (Graphical User Interface)

2. Programmation par le fichier d'analyse 3. Mixage des deux mthodesmentionnes Dans la mthode interactive ou GUI, les donnes peuvent tre entres par les fentres prpares par ce logiciel Elments Finis. Ces fentres sont plus simples pour les petits

calculs, mais elles ne sont pas adaptespour les calculs compliqus et rptitifs. A titre d'exemple,une fentreest prsente dansla Figure 6-1.

ld(r||wahdrfficsmrd(ffi4

.lFl

lvnr

lro*l

lqr{

Figure6-1 Utilisationdu logicield'Elments Finis par GUI.

r64

6'0

Chapitre 6. Modlisation lastiqueet Lasto-plastiaue ioints soudspar rtoint des

Dans la mthode de la programmation pour le fichier d'entre, les paramtres du comportement des matriaux, la gomtrie et la rsolution sont facilement contrls par la mthode de programmation. Un autre avantagede cette application est la mise en vidence dansles calculscompliqusqui se font en plusieurspas,en utilisant par exemple,la technique << sub-modeling >, << superlment > et le calcul de la mcaniquelinaire de la rupture. I-e Mixage des deux mthodes mentionnes(GUI et Programmation) est gnralementplus efficace dans les calculs par Elments Finis. Dans la Figure 6-2 qui prsenteune utilisation gnralede la programmation, le fichier TEXT est frquemmentutilise.

@*i.*o.n. - S
* esshJ((J t lod (:) lqr i Efg V

U) ,P1 t 2)e-*a- ,gr*,,W*.

. " l=l l:-J i"iert.'lt i i'

i*l ,#':

rxr X ift fil

fr wtLi+ B
lfo

rnlt:i adcr

&

'

$ffi*;-

C ?3 hiilo.dt Htri)ffi

ri,
:r o iSDb

srt* 'f u
# i5 mt #

* i #i3

Figure Utilisation logiciel la mthode programmation. 6-2 du par de

Dans la modlisation des joints souds, la mthode mixage est utlise. Selon les caractristiques particulires des joints souds par point dans le cas de I'analyse tridimensionnelle,la mthodedes sous modelages<< sub-modeling> est donc utilise tl23l.I-a longueur et la largeur des joints soudssont trs grandespar rapport leurs paisseurs. est Il vident qu'il est alors trs difficile d'obtenir de bons maillages dans l'paisseurdes joints souds.Il est ncessaire d'avoir des maillagesraffins dansl'paisseurdes tles minces.Pour obtenir cette exigence, le technique du sous modelage est donc choisie pour les calculs des joints soudspar point. En plus cette mthodediminue le temps du calcul.

r65

6.L Modlisation desjoints soudspar point Mono Point, Triple Points et Multiple Points
Le mixage de la programmation et GUI est utilis pour le calcul des Elments Finis. Le paramtragede la gomtrie, du chargementet les critres de la solution sont appliqus. En raison de la prsence d'paisseur trs mince dans les joints souds par point, un sous modelageest donc ralis.L'algorithme de cette analyseest prsentdans la Figure 6-3.

tI
|

Modle Global
| SousModle

Figure 6-3 Algorithme de la mthode sous modelagepour la mthode des Elments Finis.

Les calculs par la mthode des Elments Finis pour les joints soudspar point sont donc raliss partir du mixage des deux mthodesmentionnes,et I'algorithme de cette analyse pour tous les joints soudsest prsentdans la Figure 6-4. Dans cette figure, les paramtres des matriaux et des gomtriessont donns.Selon le type de calcul utilis (lastiqueou lasto-plastique),la solution du problme avec les paramtresdu contrle de la procdureest

r66

6.1

Chaoitre 6. Modlisation lastiaue et lasto-plastiauedes ioints soudsoar point

donc ralise. Dans la prparation du calcul, ces paramtres sont arangs. Ce pas est trs important dans le cas des calculs non linaires. mthode calcul nonlinaire a dj t La de mentionne dansla section2.2.3.

Insertiondesparamtres matriaux

Insertiondesparamtres la gomtrie modle du

Calcullastique
lasto-plastiq

Prparation calcul du lasto-plastique

Solutiondes systmes linaires

Solution des systmes

Vrificationdesrsultats
par Figure6-4 Algorithmegnral calculdesjoints souds point danslescaslinaires non-linaires. du et

par Dansla modlisation joints souds point mentionne des dansla Figure64,la solution (Figure6-3). dessystmes linaires non-linaires et utilise le sousmodelage

t67

6.1.1. ModIisationdesjoints souds point Mono Point par

La modlisation des joints soudspar point Mono Point est ralise dans deux cas, lastique et lasto-plastique. Dans la Figure 6-5, des maillages, I'agrandissement maillage et les du conditions aux limites sont prsents.

b)

168

6.1.1

Chapitre 6. Modlisation lastique et lasto-plastiquedes ioints soudspar point

c)

d)

e) globale par Figure6-5 Modlisation d'unjoint soud pointMonoPoint a) Maillageb) Agrandissement du e) maillage Conditions limites chargement Modle c) aux et d) complet Agrandissement dumaillage.

r69

6.1.I

Chapitre 6. Modlisation lastique et lasto-plastiaue des ioints soud.s p,oint par

Dans la Figure 6-5, la modlisation globale d'un joint soud par point Mono Point est prsente.La Figure 6-5c reprsenteles conditions aux limites et le chargement.La symtrie est applique pour diminuer le temps de calcul par Elments Finis. La force est applique du ct droit du modle et le ct gaucheest fix par bridage (Figure 5-12). Dans la Figure 6-6, le sous-modle prsent. est

Vu Haut Vu Devant

Vu Isomtrique

b=48 mm, t=1.2 mm, d=6 mm, ft=J mm, Gap=9.13otFigure6-6 SousModlisation globaled'un joint soud point Mono Point. par

Dans la Figure 6-6, le sous-modle pour une prouvette soude point avecla gomtrie par mentionne dansla Figure5-2, estprsent. Dansla Figure6-6, b, t, d, R et Gap sont respectivement largeurdu joint soud point, l'paisseur la tle mince,le diamtre la par de du point soud, rayondu point soud I'interfaceentrelestlessuprieure infrieure. le et et I-esjoints souds point sontanalyss calcullastique calcullasto-plastique par par et pour valuer champdescontraintes la distribution contraintes le et des autour joints souds des par

170

6.LI

Chapitre 6. Modlisation lastiaue et lasto-plastiauedes ioints soudsnar ooint

point. L'utilisation de la technique sous-modelagedoit tre vrifie le modle global et le sous-modledans les lignes de coupure des frontires (cut boundary lines) ll24l. Dans la Figure 6-7,|e sousmodleest schmatiquement prsent.

Sousmodle

----- sourooer. Figure6-7 Schma sousmodle du d'uneprouvette par soude pointMono Point.

La distributioncritiquedescontraintes autourdu point soud doncconsidre est pour valuer lesjoints souds point. Dansle casde I'existence par d'un seulpoint soud, point estle ce point critique. La distribution des contraintespour les joints soudsest gnralement reprsente l'volution descontraintes par dansla directionlongitudinaledesjoints souds. Les rsultats permettent de dfinir les facteurs de concentrationde contrainte et de dformation. facteurs Ces sontdfinispar:
tlt

,- -- o'* -

o,

(6-l) (6-2) (6-3)

oTo, c, et , sont respectivement facteur de concentration le de

No--

,. - oio o, ioes
kr, oi,

k' = -

k,, k*

contraintes,le facteur de concentrationde contrainteslasto-plastiques, facteur de le dformations lasto-plastiques, contraintelastiquedans la direction x ou la direction la t7l

6.1.I

Chanitre 6- ModIlsation Ia.rtiaue et Ia.to-nlnstiaue des ioints soul nnr noint

longitudinale du joints soud, la contrainte lasto-plastiquedans la direction x, la contrainte globale et la dformation globale. La contrainte globale et la dformation globale sont donc dfinies par les relations suivantes:

o"=!
" o" t, =

bxt

rc4)
(6-5) et E sont respectivementI'amplitude de la force appliquee, la largeur du joint

Fo, b , t

soud,l'paisseurde la tle mince et le module de Young. Dans la Figure 6-8, le dplacement et les valeurs des dplacements sont respectivement prsents.

a)
toD t 50tnT10r sllP.l 5ll8 .I 1I G.l /D(IDED

l\l.

ust[

R5?5.0 DE(

,1t,|712 .20?365 .622tSA

.211

1.4s2

1 . 4 56

b) Figure Prsentation dplacement modle gnral joint soud pointencisaillement par 6-8 du du du MonoPoint
par calcul lastique.

172

6,1.1

Chapitre 6. Modlisation lastiaue et Iasto-plastiouedes joints soudspar point

a)

!POX..OO34?a SX-{s{, t08

b)

Figure6-9 Champdescontraintes pour le joint soudpar point Mono Point (lasto-plastique) Champdes a) contraintesb) Champdescontraintes autourdu point soud(la contrainte maximaledansla directionx et la dformation totaleen ce point).

L'volution des contraintes fonction de la distanceest frquemment en utilise dans ces courbesU071. Cetteprsentation rsultats des considre tous les cas d'analysed'Elments Finis (lastique et lasto-plastique). L'volution des contraintes mentionnesa un comportement antisymtrique rapportau centredu point soud. par Mais, dansles conditions

173

6.1.1

Chapitre 6. Modlisati.onlastiaue et lasto-plastiaue des ioints soudspar point

antisymtriques. Selondiffrentesconditions rellesces volutionsne sontpas parfaitement n'est pasparfait. aux limites desdeuxctsdujoint soud, phnomne ce par mentionns ralise descalculslastiques Dansla Figure6-10, l'volution desfacteurs et lasto-plastiques prsente. est

22 20 18 I6 T4 L2 10 8 6 4 2 0 0

,k6,k -o- kt --E- ko --- ke

0,25

0,5

0,75

Figure 6-10 Evolution du facteur de concentration de containtes lastiques(, au facteur de concentration de contraintes lasto-plastiques ko et du facteur de concentration de dformations lasto-plastiques dans la direction longitudinale du joint soud par point Mono Point ( f est le rapport entre la distance dans la

direction x qui commence au centre du point soudet le diamtre du point soud).

que les contraintes lastiques, lasto-plastiques les et Dans la Figure 6-10, on remarque dformations augmentent voisinagedu fond d'entaille.Elles diminuentavecla distance au dans la direction longitudinale du joint soud par point. Ce modle d'volution est frquemment obtenu. montrentI'existence contraintes Les valeursdesfacteursde concentration contraintes des des autourdu point de soudure. pratique, En I'amorage la propagation fissures et des trsleves commencent cet endroit. Dans l'approche volumtrique,on trace le diagrammebilogarithmique l'volution de la contrainte de d'ouverture fatigueen fonctionde la distance. en de Pour tracer cette courbe, il est ncessaire dterminerla direction d'ouverture et la contrainte ouverture fatiguecorrespondante. d' en

174

6.1.2Modlisationdesjoints souds point Triple Points par

La modlisation joints souds point Triple Pointsestralise par des dansdeux cas,lastique et lasto-plastique. la Figure6-11,les maillages, Dans I'agrandissement maillageet des du conditions limites sontprsents. aux

b)

r75

6t2

Choni.trp

6.

Mnd.,lj:ation

,lnstinrtc p.t .lasto-nln.sti.nu.e. d.es ioi.n.t. nud.s nnr nninl

c)

d)

e) Figure6-l l Modlisation globale joint soud pointTriplePoints a) Maillageb) Agrandissement par d'un du
maillage c) Conditions aux limites et chargement d) Modle complet e) Agrandissementde maillage.

176

6.I .2

Chapitre 6. ModIisation lastiaue et lasto-plastiquedes ioints soudspar point

Dans la Figure 6-11, la modlisation globale d'un joint soud par point Triple Points est prsente. La Figure 6-llc reprsente des conditions aux limites et le chargement. La

symtrie est applique pour diminuer le temps de calcul par Elments Finis. La force est appliquedu ct droit du modle et le ct gaucheest fix par bridage (Figure 5-12). Dans la pour le point soudcritique. Figure 6-12,le sous-modle prsent est

Joint soudpar point cri

Vu Haut

Vu Isomtrique b=48 mm, 1=l.l mm, d=6 mm, R=3 mm, Gap=Q.l3mm Figure 6-12 SousModlisation globale d'un joint soudpar point Triple Points.

pour par Dansla Figure6-12,le sous-modle uneprouvette soude point avecla gomtrie Dans Figure 6-12, b , t, d , R et Gap sont mentionne dans la Figure 5-2 est prsent. par respectivement largeurdujoint soud point,l'paisseur la tle mince,le diamtre la de du point soud, rayondu point soud I'interface le entreles tlessuprieures infrieures. et et

177

6.1.2

Chapitre 6. Modlisation lastioue et lasto-plastiaue des ioints soudspar ooint

L'utilisation de la technique sous-modelage doit tre vrifie pour le modle global et le sous modle ainsi que les lignes de coupure des frontires (cut boundary lines). Dans la Figure 613,Ie sous-modleest schmatiquement prsent.

Figure 6-13 Schmadu sous-modled'une prouvettesoudepar point Triple Points.

La distribution des contraintes autour du point soud pour le point critique est donc considrepour valuer les joints soudspar point. Dans la Figure 6-14, le dplacement et les valeurs des dplacements sont respectivement prsents.

a)

t78

6.L2

Chaoitre 6. Mod.lisationlastioue et lasto-plastiauedes ioints soudspar point

DD( .I. I

b)

Figure6-14 Dplacement modle gnral joint soudpar point en cisaillement du du Triple Poins par calcul lastique.

a)

t79

6.1.2

Chapitre 6. Modlisation lastiaue et lasto-plastioue des ioints souds oar point

b)
Figure6-15 Champdescontraintes pour le joint soud point Triple Points(lasto-plastique) Champs par a) des contraintes b) Champ des contraintesautour du point soud(la contraintemaximale dans la direction x et la dformationtotaleen ce point).

Dans la Figure 6-16, l'volution des facteursmentionnsest donc ralisepar calculs lastiques lasto-plastiques. et

22 20 18 t6 t4 L2 10 8 6 4 2 0 0

,k6,k

0,25

0,5

0,75

E
Figure 6-16 Evolutiondu facteurde concentration contraintes de lastiques(, du facteurde concentration de contrainteslasto-plastiques et du facteur de concentation de dformations ko lasto-plastiques dans la direction longitudinaledu joint soudpar point Mono Point ( f est le rapport ente la distance dansla direction x qui commence centredu point soudet le diamtredu point soud). au

r80

6.1.2

ChaPitre 6. Modlisation lastique et lasto-plastique des ioints soudspar point

La Figure 6-16 montre que les contraintes lastiques, lasto-plastiques les dformations et augmentent voisinage fond d'entaille.Elles diminuentavecla distance au du dansla direction longitudinale joint soud point. du par

l8l

par Modlisationdesjoints souds point Multiple Points 6.1..3

par point Multiple Points est ralisedans deux cas, La modlisation joints souds des lastique lasto-plastique. la Figure6-17,lesmaillages, et Dans I'agrandissement maillage du et desconditionsaux limites sontprsents.

a)

b)

r82

6.L3

Chaoitre 6. Modlisation lastique et lasto-plastiauedes ioints soudspar point

c)

d)

e)
Figure 6-17 Modlisation globale d'un joint soud par point Multiple Points a) Maillage b) Agrandissementdu maillage c) Conditions aux limites et chargementd) Modle complet e) Agrandissementdu maillage.

183

6.1.3

Chapitre 6. ModIisation lastique et lasto-plastique des ioints soudspar point

Dansla Figure6-17, la modlisation globaled'un joint soud point Multiple Pointsest par prsente. Figure6-l7c reprsente conditionsauxlimites et le chargement. symtrie La les La pour diminuerle tempsde calcul par Elments est applique Finis. La force est applique du ctdroit du modleet le ctgauche fix par bridage(Figure 5-I2). Dansla Figure6-18, est le sous-modle prsent pour le point soud est critique.

par Joint soud point critique

A.1

Vu Isomtrique

-b=48 mm,t=l.l mm,d=6uun" R=3mm,Gap=g.13 Figure6-18 SousModlisation globaled'un joint soud point Multiple Points. par

Dansla Figure 6-18, le sous-modle pour une prouvette par soude point avecla gomtrie mentionne la Figure5-2estprsent. dans La distribution contraintes des autourdu point soud point critiqueest donc considre par pourvaluer joints souds point. les par Dans la Figure 6-19, le dplacement Ies valeursdesdplacements respectivement et sont prsents.

r84

6.1.3

Chaoitre 6. Modlisation lastioue et lasto-plastiaue des ioints souds Dar Doint

IODII

SOLI'TI(II

AN

SIIE .I TIIE.I ,/DOTDED (aYcl u$r! RSIS.0 DID( sID(

Multiple Pointspour le Figure6-19 Dplacement modlegnral joint soudpar point en cisaillement du du calcullastique.

Dans la Figure 6-20, la distribution des contraintes autour du point de soudurecritique mentionn dansla Figure6-18 estprsente.

E iLo.tt

rtta.}s :tt.ttl 2t2-1t a!.t:t -tlt-a,S -ta.iail.. oal-llla tlrr..|l

a)

r85

6.1.3

Chqpitre 6. Modlisation lnstiaue et lasto-plastique des ioints souds oar point

b)
Figure 6-20 Champde contraintespour le joint soudpar point Multple Points (lasto-plastique)a) Champdes contraintesb) Champdes contraintes autourdu point soud(la contrainte maximaledansla directionx et la dformation totaleen ce point).

Dansla Figure6-2L,l'volutiondes facteurs mentionns donc alisepour les calculs est lastiques lasto-plastiques. et

22 20 18

k7,k6 ,kg

r6
l4

--o- kt --tr- ko --- ke

t2 10
8 6 4 2 0 0 0,25
0r5
Figure6-21 Evolutiondu factew de concentationde contraintes lastiques du facteurde concentration k,, de contrainteslasto-plastiques et du facteur de concentration dformations fto. de lasto-plastiques dans la directionlongitudinaledu joint soudpar point Mono Point *, (( estle rapport entre la distancedansla

0,75

directionx qui commence centre point soud le diamtre point soud). au du et du

186

6.1.3

Chapitre 6. Modlisation lastiaue et lasto-plastiquedesjoints soudspar point

lastiques, lasto-plastiques les et Dans la Figure 6-21, on remarqueque les contraintes augmentent voisinagedu fond d'entaille.Elles diminuentavec la distance au dformations par du dansla directionlongitudinale joint soud point. par Il est doncvidentque tousles typesdejoints souds point en cisaillement Mono Point, autour du point de soudure Triple Points, Multiple Points) ont le mme comportement maximalesetrouveen fond du point soud. critique.La contrainte

187

6.2 Prvtsionde la vie en fatigue desjoints soudspar point par I'approchevolumtrique

Dans I'approchevolumtrique,la prvision de la dure de vie en fatigue dpendde la distribution des contraintes lasto-plastiques dans la direction d'ouverture.Dans les joints souds par point, ces paramtres doivent tre bien distingus. Des observations exprimentales montrentque la fissures'amorceen fond du point souddansla directionde l'paisseur. partir de cesobservations, contrainte A la d'ouverture sadirectionsontd6finies et et prsentes laFigure 6-22. dans

Directiond'ouverture la fissure de

Contrainte d'ouverture fatigue en

Figure 6'22 Schmade la direction d'ouverture de la fisswe et la contrainte d'ouverture en fatigue.

La comparaison entre la Figure 6-22 et les Figures6-6, 6-12 et 6-18 montre comnientla formulationde I'approche volumtrique peuttreadapte sousla forme:

k, = * lo*(z) (L- z)dz ' "' Z"or

' tk)= ky l-.doryk) o*\z) os, o*(z), ctz

ttq

6-4)

(6-5)

Z"t et Z sont respectivement le facteur de rduction en fatigue, la

contrainte globale, la contrainte d'ouverture x dans la direction d'ouverture de la fissure, la distanceeffective et le gradient de contrainte relatif. La distribution des contraintesd'ouverture pour la direction d'ouverture est prsentedans la Figure 6-23. L'volution des contraintes lasto-plastiqtresd'un point de soudure critique dans les joints soudspar point est prsente. dbut de cette volution, la contrainte lasto-plastiquea la Au jusqu' la contraintemaximale. Aprs ce maximum, la valeuro,. Cette contrainteaugmente contrainte diminue comme le montre laFigve 6-23.

188

6.2

Chnpitre 6. ModIisation lastiaue et lasto-plastiouedes ioints soudspar point

Compression

ZoneFondue

Figure6-23Schma la distribution contraintes de des dansl'paisseur d'un point soud critique.

Lescourbes utilises pourI'approche volumtrique prsentes la Figure6-24IlO7l. sont dans

0,50 0,25
500 A

(z*o*)

(z"tt,o"t)

0,00 -0,25

x00 N

o
300

x x

-0,50

II
o

ilI
-o,75

x
200

-1,00

+- .0,05 0,1
Zm

-----__

-1,25

a)

0,15
Zeff

zw)

1 89

6.2

Chapitre 6. Modlisation lastique et lasto-plastiaue des ioints souds par point

600

2,O0

500 { ^400 N x
cl

(4o,o- ) ---

1,00

(4n,o.n )
0,00

x
yl \l

a
\l

-1,00 -rd rrr ---\ -2,00

300

II

}-i

iI I I I I I I

200

-3,00

0,05
Z1

0,1

0,15
Zeff

b)

o,2 (mm) z
1,50 1,00

(z-,oo' )

(2"n,6"t)

0,50 0,00

-0'tg
I IIt III
-1,00 -1,50 -2,00 -2,50

0 ,0 5 c)

0,10

0,15
Zeff

zm

z (mm)

0, 20

Figure6-24Distribution contraintes l'paisseur pointsoud des pourobtenir facteurs dans d'un critique les de rduction fatigue jointssouds point a)MonoPoint b) TriplePoints c) Multiple par en des Points. Selon les volutions prsentes dans la Figure 6-24,les valeurs de la contrainte maximale, de la contrainte effective et de la distance effective peuvent tre dfinies. La mthode volumtrique utilise gn&alement le facteur de rduction en fatigue ainsi que la courbe de

190

6.2

Chapitre 6. Modlisation lastique et Iasto-plastiouedes ioints soudspar point

rsistanceen fatigue sur des prouvetteslisses pour obtenir la courbe de rsistanceen fatigue d'une prouvette entaille. La procdure de la prvision de la dure de vie en fatigue a dj t mentionnedans les Figures (3-7) et (3-8). Cette courbe peut tre utilise pour obtenir les courbes entailles pour le mme type d'prouvettes, mais avec des caractristiques gomtriquesd' entailles diffrentes. Cetteprocduren'est pas conseilledansle cas de I'existencedes courbeslisses.Pourtant,la ralisation d'prouvetteslisses pour tous les types d'prouvettesn'est pas facile ou pas par faisable.Dans le cas de I'utilisation desjoints souds points,la dfinition d'une courbede rfrence en fatigue n'est pas possible. La procdure prsente dans la Figure (3-8) est applique pour obtenir la prvision de la dure de vie des joints souds par point en cisaillement. La prvision obtenue par la mthode volumtrique est prsentedans la Figure 6-25. Dans la Figure 6-25,les rsultats de l'approche volumtrique et les rsultats exprimentaux sont compars.Dans la Figure 6-25a,la courbe pour les Mono Points est utilise pour obtenir la courbe lisse virtuelle. Les rsultats obtenuspour les joints soudsTriple Points et Multiple Points prsententun bon accord avec les rsultatsexprimentaux. Pour les joints soudspar point Mono Points, les rsultats obtenus par I'approche volumtrique et les rsultats exprimentaux sont identiques. Avec cette mthode, dans la Figure 6-25b, les rsultats obtenuspar les essais,pour desjoints soudspar point Triple Points, sont considrscornme rfrence.

I 0000

z
frr
cr,

I 000
O A o -2-3s M o n o p o i n t - R s u l t a te x p r i m e n t a u x Triplepoints-Rsultats xprimentaux e Multiplepoints-Rsultats xprimentaux e v l - M o n o p o i n t - C o u r b ed e I ' a p p r o c h e o l u m t r i q u e T r i p l ep o i n t s . C o u r b e e I ' a p p r o c h e o l u m t r i q u e d v d v M u l t i p l ep o i n t s - C o u r b e e I ' a p p r o c h e o l u m t r i q u e

100 1,0E+05
a)

N ombr"t Nn d"ttofcles,
191

l,0E+07

6.2

Chapitre 6. Modlisation lastioue et lasto-plastioue desjoints souds oar point

10 0 0 0

z
rt
1000 1 tr
A M o n o p o i n t - R s u l t a t se x p e r i m e n t a u x Trip le p o ints -Rs ultats experimentaux M ultip le p o ints -R sultats exp erimentaux 1 - M o n o p o i n t - C o u r b e d e I ' a p p r o c h ev o l u m t r i q u e 2 -Trip le p o ints -C o urb e d e I'ap pro che vo lumtriq ue M u l t i p l e p o i n t s - C o u r b e d e I ' a p p r o c h ev o l u m t r i q u e

o
-3-

100 1,0E+05

1,0E+06

1,0E+07

Nombre de cycles,Np
b)

Figure6-25 Comparaison entreles rsultats par obtenus I'approche volumtrique les rsultats et exprimentaux a) Rsultats joints souds point Mono Pointsconsidrs des par corlme rfrence b) Rsultats joints souds des par point Triple Pointsconsidrs cornme rfrence.

Dansce cas,les rsultats exprimentaux ceuxde I'approche et volumtrique sontidentiques. I-es rsultatsobtenusdans la Figure 6-25 montrent que les prvisionspar la mthode volumtrique donnede bonsrsultats pour lesjoints souds point en cisaillement. par

192

6.3 Discussions
Dansl'tudedesjoints souds point, deuxstades par peuvent tre distingus I'amorage la et propagation. Dansla Figure5-13,cesstades sontindiqus. L'expression nombrede cycles du la rupturepeuts'crirecommesuit: N*=N,+N, (6-6)

N^, N, et N, sontrespectivement cycles la rupture,les cyclespour I'amorage les les et cyclespour la propagation fissures. localisation I'amorage fissureschange des La de des en position avec la valeur du chargement pour les joints soudspar point en acier [90]. Autrementdit, plus la valeurdu chargement augmente, plus la fissurese dplacede la zone affecte (ZAT) au mtalde base. thermiquement Dansle casd'existence d'un chargement trs leven fatigueoligocyclique, nombrede cycle d'amorage ngligeable. le est Dansle casde I'endurancelimite, le nombre de cycles de propagationest petit vis--vis de celui de I'amorage. conditions Ces sontprsentes la formesuivante: sous

Fatigue oligocycliq*=l+ + o +lv*=lrr)

(6-7)

Endurance limite = (+ -+0 +n"= n

Il est donc ncessaire d'utiliser une mthode convenable pour chaquecas.En pratique,les joints souds par point travaillentdans le cas de I'endurance limit (105-10? cycles). Commepar exemple, dansle casdestlesmincespour les carrosseries d'automobiles. Cette hypothse t applique SherrySheppard a par |17-1181, Nin Pan [105] et H-F. Henrenson joints souds pointencisaillement Il que par [100-102]. estdoncvident l'tudedes conduit s'intresser I'amoragedesfissuresau voisinagedu point soudcar en gnral,les joints soudssont utiliss pour assembler deux tles minceset donc la force appliqu n'est pas grande. L'analyse des joints soudspar point avec une approche la fatigue peut tre faite par de diffrentesmthodes. L'approchevolumtriqueet les mthodes< Point Chaud)>peuvent potentiellement utilisespour I'estimationde la durede vie en fatigue.I-es mthodes tre << Point Chaud> utilisentgnralement contrainte la maximaleau voisinage fond d'entaille du constitupar le point soud.I comportement fatiguene dpendpas uniquement la en de valeur maximalede la contrainte fond d'entaille.Il est donc ncessaire considrer en de la

t93

6.3

Chapitre 6. Modlisation lastiaue et lasto-plastioue des ioints souds par point

taille de la zoneplastifieet le gradientdes contraintes La t66-671. mthodevolumtrique calculela valeurde la contrainte effectiveet de la distance effectivecorrespondante point au minimal du gradient contrainte de d'ouverture. L'intgrattonpeuttrefaite dansun volumeen fond d'entaille.Mais,l'tude de N. Kadi [67] montrequele calculde la contrainte effectiveet du facteurde rductionen fatiguedansun volume,surunesurface sur uneligne conduisent et sensiblement la mme valeur. Il est donc suffisant d'utiliser la formulation basesur I'intgrationsur une ligne (quations (3-18) et (3-19)).L'approchevolumtrique les et mthodes Neuberet Molski-Glinkaont djtcompares dansdesdiffrentscasde durede vie en fatigue I70).La comparaison avecles rsultats exprimentaux montreque I'approche volumtrique conduit uneestimation convenable la durede vie en fatigue. de Il est ncessaire d'avoir une courbede rfrenceen fatigue pour I'approchevolumtrique (Figure3-7). Dansl'tude desjoints souds point en cisaillement, courbede rfrence par la en fatigue est pratiquementimpossible obtenir. Par consquent, une autre procdure prsente dansla Figure3-8 estutilisepour estimerla durede vie en fatigue. La courbede rfrence fatgue est obtenue partir desrsultatsexprimentaux du facteurde rduction en et en fatigue desjoints soudsmono point. Cette courbeest utilise pour pourvoir appliquer I'approchevolumtrique dansles casdejoints souds fiple pointset multiple points (Figure 6-25a).Dans la Figure 6-25b, la courbe de rfrence fatigue est obtenue partir des en rsultats exprimentaux du facteurde rductionen fatigue patir desjoints souds et Triple Points.La comparaison pour les trois typesdejoints souds entreles rsultats exprimentaux par point en cisaillement ceuxobtenus I'approche et par volumtrique donnedoncune bonne estimation la durede vie en fatigue. de Les valeursdes facteursde rductionsen fatigue sont obtenues moyen de calculs aux au ElmentsFinis. Le facteur de rduction en fatigue est calcul partir de la contrainte d'ouverturede fatigue dans la directiond'ouverturedes fissures.I-es Figures6-22 et 6-23 reprsentent l'volution de la contraintelasto-plastique d'ouverturede la fissure dans la directionde I'intgration.En raisonde l'paisseur mincedestlesutilises trs danslesjoints par souds point en cisaillement, technique sous-modelage utilise. la du est L'approchevolumtrique calculela contrainte effectivedansle volumedtermin le point par du gradientde la contrainterelative minimum. Cette approche ralise partir d'une est intgration numrique. consquent, nombrede maillesdansla directionde l'intgration Par le doit tresuffisantpour obtenirunebonne valuation contraintes des lasto-plastiques.

t94

6.3

Chaoitre 6. Modlisation lasti.oueet lasto-olastiaue des ioints sou.dsoar ooint

prenanten compteI'interface entre les L'utilisation d'un modlecomplet tri-dimensionnel plus exacts.L'influence de I'interfaceest tles infrieureet suprieure donnedes rsultats trs importantepour les joints soudspar point mono point [24]. L'interface n'est pas prcdentes considre dansles tudes celle-ci : [90-97,100-104] Dans le calcul tri-dimensionnel, symtrieest utilisepour tous les types de joints souds la par point en cisaillement. Des conditionsaux limites similaires la situationrelle sont par appliques. Selonles tudes exprimentales lesjoints souds points,le comportement sur (ZAT) disparatds le dbut du chargement particulierde la zone affectethermiquement par cyclique 105-106]. dure vie en fatigue La de estime cettesimplification un est [100-101, peu moindreque la durede vie relle.L'tude exprimentale ZAT estgnralement de faite <ABI> (Automated par la mthode prcdentes Ball ldentation)ll25-L261. Selondestudes par la des des Finis a t [00-101, 105-106], modlisation joints souds la mthode Elments ralisesansconsidration effets de ZAT. Cette approchesimplifie le problmed'un des par calcullasto-plastique compliqu joints souds point en cisaillement. des L'exactitudedes par rsultats obtenus la mthode volumtrique par les rsultats et exprimentaux montreque I'effet deZAT estngligeable dansle casd'endurance limit. Dans l'tude desjoints souds, modede dplacement diffrentpar rapportaux joints le est soudspar point. Selon le nombre de points soudsdans les joints soudspar point en cisaillement, moded'endommagement le change. pour lesjoints souds point Mono Point et par La comparaison par entreles rsultats obtenus point Triple PointsmonEeque lesjoints souds multi pointsont unersistance fatigueplus en que leve celle desjoints souds monopoint. L'tude desjoints souds multi pointsmontre quela gomtrie, localisation, le diamtre point soudsontdesparamtres la et du importants pour la durede vie en fatigue.Selonles rsultats exprimentaux obtenus, rsistances les en fatiguedesjoints souds Triple Pointset Multiple Pointssontapproximativement gales. Ce rsultatvient de I'effet de la gomtrie joints souds. que des C'est--dire la distance entre les pointssouds aussiimportante la disposition que est dansla rsistance fatigue.Dansle en casdesessais fatiguesur desprouvettes de par soudes point Triple Points,la distance entre les points souds plus grandeque dansles joints souds est Multi Points.L'influence de la distance changeen fonction de la sollicitationmcanique joints souds tensionou en des en compressionU241. Le, comportementdes joints souds par point Triple points en compressionest meilleur celui des joints soudsMultiple Points [1241. Les tudes exprimentales montrentque le mtalfondu a une limite d'lasticit une contrainte et ultime

195

6.3

Chapitre 6. Modlisation lastiaue et lasto-Dlastiaue des ioints souds par point

plus leveque celle du mtal de baseU00-1011. des Selonle comportement joints souds par pour I'analysedesjoints souds point en par point, diffrentes mthodes sont proposes de par de cisaillement, exemplela mthode Bandeet la mthode Taylor [85, 9U. Multi Points montre que ce joint soudest plus L'tude exprimentale joints souds des dansle casde I'endurance rsistant dansle casde la fatigueoligocycliqueet moinsrsistant par est limite par rapportaux joints souds point Triple Points.Ce phnomne attribuaux pour lesjoints souds Multi Le deuxpointssouds supplmentaires. nombrede pointssouds, la Pointsinfluencedonclgrement rsistance la fatigue.

196

CHAPITRE 7

Chapitre 7
7. Conclusiongnrale
7.0 Conclusiongnrale
par et en Desjoints souds point, utilissdansI'industrieautomobile notamment carrosserie du cycliques. La caractistique chargement sont gnralement soumis des chargements par le appliqu lesjoints souds point conditionne choix de I'analyse. sur cycliquetrs lev (fatigueoligocyclique)et dansle cas d'un Dansle cas d'un chargement chargementcyclique moyen (endurancelimit), I'amorage de la fissure se situe (ZAT). thermiquement respectivement le mtaldebaseet dansla zone affecte dans de limite, la Selon que nous sommesen prsence fatigue oligocycliqueou d'endurance peuventrespectivement observes. propagation I'amorage fissures tre L'analysedes et des joints souds point est doncdiviseselonune approche type mcanique par linairede la de ruptureou de fatigue. par limite En pratique,les joints souds point travaillentdansle domainede I'endurance -107 cycles). (105 pour l'tudedesjoints souds, tproposes diffrentsauteurs. par Diffrentes mthodes, ont proposes introduisent hypothses des simplificatrices. Lesmthodes pour les joints souds par L'analyse tri-dimensionnelle I'approchela plus convenable est point. Dans ce cas,les conditionsaux limites et la gomtrie joints souds du point des et soud sontprisesencomptedansle calcul. par L'effet de I'interfaceest importantdans le cas desjoints souds point en cisaillement mono point. Ce paramtrc augmente I'excentricitde I'applicationde la force et donc, la par Les valeurde la contrainte mrimaleprsdu point soud. joints souds point Mono Point sontdoncplus sensibles la gomtrie I'interface. de des la Pourrendrre comptede I'amorage fissures des dansl'paisseur tlesminces, conffainte effective,la distance effectiveet le facteurde rduction en fatiguesont calculsdanscette rgion. (ZAT), pour les joints soudspar point en L'effet de la zone affectethermiquement dansle cas de I'endurance limite, disparatau dbutdu chargement cyclique. cisaillement 198

7.0

Chapitre 7. Conclusion

Dans le cas de I'endurance limit, le nombrede cyclespour la disparitiondes contraintes rsiduelles, ngligeable rapportau norhbre cycles la rupture. est par de L'obtention d'une courbe de rfrenceen fatigue pour les joints soudspar point est pratiquement impossible.I-es courbesde rfrence pour I'approchevolumtriquesont donc obtenues partir des rsultatsexprimentaux pour les joints soudsmono point et hiple points.La comparaison avecles rsultats exffrimentaux et I'approchede fatigueconfirmela validitde I'approche volumtrique pour I'estmation de la durede vie en fatigue.

r99

REFERENCES

200

tI 1 12 I t3 I I4 | t5 I t6 I 17 I I8 ] t9 I t l0 I t 1l I I 12 I t 13 I I 14 1 t 15 I t 16 I

W.J.M.Rankine,Minutesof theProceedingof thelnstituteof CivilEngineers,London,2,lO5,(1843). G. S. Campbell, Lahey,International R. Journal Fatigue, No. 1,25, (1985). of 6, H. Neuber, Verlag,FRG, (1958). Kerbspannungslehre, Springer R. E. Peterson, Stress concentration factors, JohnWiley & Sons,Inc. (1974). G. N. Savin,V. I Tulchii, Handbook thestress on concentration, High SchoolPublisher House,(1976). S. J. Hardy,N. H. Malik, International Journal Fatigue,14,No; 3, 147,(1992). of H. Tada,P. C. Paris,G. Irwin, The stress analysisof crackshandbook, Corporation,1973, @el Research second edition,ParisProductions, Inc., St.[.ouis,(1985). M. Creager, C. Paris,International P. Journal Fracture of Mechanics, No. 4, U7,(1967\. 3, D. P. Rooke,F. I. Baratta, J Cartwright, D. Engineering Fracture Mechanics, 14,3g7,(1981). C. S. Shin,R. A. Smith,Proceeding International of Symposium Microstructure mechanical and behaviorof (1985). Materials, Xian, China, D. Kujawski, Estimationof stressintensity Fatigue and Fracture of factor for small cracksat notches, (1991). Engineering Materials,14, No. 10,953-965, H. Neuber, Kerbspannungslehre, Verlag,FRG, (1958). Springer V. Weiss,ASME Paper (1962). No. 62-WA-270, AmericanSocietyof Mechanical Engineers, G. Glinka, A. Newport, Universalfeatures of elastic notch tip stress field, InternationalJournalof Fatigue, (1987). 9, No. 3, 143-150, M. Creager, P.C. Paris,Elasticfield equations blunt crackswith reference stess corrosioncracking, to for International (1967). Journal Fracture of Mechanics, 3, 17-255, No. S. Usami,CurrentResearch FatigueCracks, on The Society MaterialSciences, T. Tanaka, Jono, of Ed. M. K. Koami,Kyoto,(1985). C.Chen,H.Pan, Collectionof papersonFractureof (Ed.C.Chen),MetallurgylndustryPress, Metals, Beijing,197,(1978). S. Thimonshenko, N. Goodier,Theoryof Elasticity, J. MacGraw-Hill, New York, (1951). H. Neuber,Theory of stressconcentration shear-strained prtsmaticbodieswith arbitrary nonlinear for stress-strain Journal AppliedMechanics, law, of 28,pp.254, (1961). Xu Kewein, He Jiawen, Prediction on notched fatigue limits for crack initiation and propagation, Engineering Fracture Mechanics, I (3), pp 504-5 (L992). 4 10, S. Bhattacharya,A. N. Kumat, Rotationnlfactor using bending momentapproach under elnstic-plastic; Situation Notch 3PB Geomery,Engineering I: (195). Fracture Mechanics, 50(4),pp 495-505, R.X. Xu, J. C. Thompson, H. Topper,Practical stress T. eryressions stress concentration regians, for FatigueandFracture Engineering of Materials pp (1995). andSEuctures, 18(718), 885-895, T. R. Wilshaw,C. A. Rau,A. S. Tetelman, generalmodelto predict the el.astic-plastic A stress distribution andfracture strengthof notchedbars in plane strain bending,Engineering FractureMechanics,Vol. l, 1912ll, (1968).

l,l7 I

t18l tle l 1.20 |

I2r l 122 l
123 |

1.24 1

M.7eng, E. Niemi, Analysis the stress of concentration facbr for a slnllow notch by the slip-linefield rneod,International Journal fatigue,l9(3), pp. l9l-194, (1997). of S. Q. Shi, M. P. Puls,A simple methodof estimatingthe maximumnormal stress and plastic zonesizeat a shallownotch,IntenntionalJournalofpressurevesselandpiping,64(1),pp.67-71,(L995).

l.zs I

126 I I27 I 1.28I

S. Q. Shi,M. P. Puls,S. Sagat, 208,pp. 243,(1994). Journalof NuclearMaterials, ANSYS Inc. TheoryReference manual. A. C. Ugural,S. K. Fenster, AdvancedStrength and AppliedElasticity,AmericanElsevierPublishing Inc.,

(le7s).
129 1 t 30 ] I 3l I .32 I t 33 I 1,34I t 35 I I 36 I 1,37I t 38 I t 39 I t 40 I t 41 I L42 I | 43 I R. E. Peterson, ^S/ress concentratianfactors, JohnWhiley andSons, New York, (1974). H. Neuber, Theoryof Notchstresses, Edwards fW. Co., Ann Arbor,Michigan,(1946). R. B. Heywood,Designingbyphoto-elasrrclry, London,Chapman Hall, (1952). and R. B. Heywood,Designing against and fatigze, London,Chapman Hall, (1962). E. Siebel, Steiler,Non uniform stress M. distributionduing fatigueloading,VDI-2,97, (1955). P. Khun, H. F. Hardraht,An Engineering methodfor estimating notch sizeeffectinfatigue tests steel, on NACA, Technical Note,No. 2805,(1952). A. Buch, Analytical approachto size and notch size fficts in fatigue of aircraft material specimens, MaterialScience Engineering, and Vol. 15,(1974),pp. 75-85. A. Buch,Fatiguestrength calculation,USA, MaterialScience Surveys, 6, (1988),pp. 465. No. T. H. Topper,M. H. El Haddad, predictionof notches Fatiguestrength basedon fracture thresholds, First (1981), International pp777-797. Conference, Stockholm, Yol.2,EMAS, Warley, A. Brand, R. Suleitin, Calcul des pieces ln fatigue, mthode gradient, PublicationSenlis-France de CETIM, 1980. Eurocode Designof SelStructures; 3, EI.tVl993-l-1: Part l.I, Generalrulesand rulesfor building,CEN,

(1ee3).
European Conventionfor Construction Steelworks:Recomrnendations the Fatigue Design of Steel for Stuctures.ECCSPublication (1985). 43, A. Palmgreem Ingenieure, (lg%), p. 339. 58, , BertchriftdesVereines M. A. Miner, Journalof AppliedMechanics, (I954),p. A-159. 12, W. Weibull, A statisticaldistribution (1951), function of wide applicability,Journalof Applied Mechanics, pp.437452. 144 I t 45 I 146 I O. H Basquin, Proc.,ASTM, 10,PartII, (1910), P.625 R. W. Smith,M. H. Hirschberg, S. Manson, S. NASA TN D-1574,NASA, April (1963). S. S. Manson, M. H. Hirschberg,Fatigue: An interdisciplinary Approach,SyracuseUniversiry Press, p. Syracuse, Y., (1964), 133. N. 147 I t 48 I 1,49I t 50 I 151 I | 52 I L. H. Coffin,Jr.,Trans., ASME,76,(1954),p.931. J. F. Tavernelli, F. Coffrn,Jr.,Trans., L. p.438. ASM,5l, (1959), R. W. Hertzberg, Defortnation and FractureMechanicsof Engineering Materials,John Willey and Sons, usA (1976). Jo DeanMorrow, Internalftiction, Dampingand Cyclic Plasticity, ASTM STP378,(1965),p.45. G. Molski, G. Glinka, A methodof elastic-plasticstress strain calculationat the notch roots, Mateial (1981); 93-100. Science Engineering 50: G. Glinka, Energy density approach to calculation of inelastic stressstrain near notchesand cracks, -508. ( Engineering FractureMechanics,I 985); 22:485

202

t 53 I t 54 I t 5S 1 t 56 I

N. A. Makhutav, The deformation fracture criteria and the strength calculationof the machinepart, (I Maschinostroyeine,98I ), pp.272. Integral, by and of Study stress strain concentationcoefficient path htdependent E. Morozov,G. Pluvinage, Issue(1996):53'U. Special of Problems Strength, energyas a citerion for fatigue Fracture,Trans. C.E. Feltner,J.D Morrow, Micro-plasticstrain hysteresis (1961). 60-MET-2 A.S.M.E., K. N. Smith, P. Watson,and T. H. Topper,A stress-strain function for thefatigue of metals,Journalof -778. pp. Vol. 5, No. 4, (1970), 767 Materials, (1992);45: 105. Technology Science W. Yao, Composite Joumalof Fatigue,(1993); 15: intensity approachfor predictingfatigue life,Internanonal W. Yao, Stress 243. Sin. Astronaut. (1988),9, A424. H. Nie, F. Wu, Acta Aeronaut. M. Hoffmann, T. Seeger,ProceedingInternationalConferenceon Fatigue of Engineering Materials and l,ondon,Vol. 1, ( I 986),pp.293-3M. Structures, (1959), pp'293-306' (Ed.G.Sines), McGraw-Hill, MetalFatigue R. E. Peterson, (1974). OxfordUniversityPress, N. E. Frost,K. J. Marsh,L. P. Pook,MetalFatigue, Yao Weixing,Xia Kiaquan,Gu Hi, On thefatigue notchfactor, k 1 ,InternationalJournalof Fatigue,Vol. (1995). 17,No. 4,pp.245-251, to intensity.approach Da-KangWang,Ming Li, Wei-Xing Yao,local stress-strainfreld Shang, De-Guang fatigue life prediction under random cyclic loading, InternationalJournal of fatigue, 23(10), (2001), pp. 903-910. of Weixing Yao, Bin Ye, Lichun Zheng,A vertficationof the assurnption antfatigue design,International 23(3),(2ffi1), pp.27| -277. of Journal Fatigue, (2000). de Universit Metz,France, de G. Qylafqu,Thse Doctorat, (2001). de de Universit Metz,France, N. Kadi, Thse Doctorat, proposal fatigue life Applicationof newrnodel N. G. Qylafku,Z. Azaui, Kadi, M. Gjonaj,G. Pluvinage, for 53-760. InternationalJournalof Fatigue,(1999);21:7 prediction on notches and key-seats, to subjected Fatigueof specimens Z. G. Qylafku, N. Kadi, J Dobranski, Azari,M. Gjonaj,G. Pluvinage, ofFatigue,(2001);23(8):689-701. pressure, Journal International loadingrole ofhydrostatics cornbined of H. Adib, G. Pluvinage,Theoreticaland numericalaspects the volumetricapproachfor fatigue life (2003);25(l\:67-76. of Journal Fatigue, components,International in estimation notched intensity The useof therctch stress T. factor T. Boukharouba, Tamine,L. Niu, C. Chehimi,G. Pluvinage, ( FractureMechanics, 1995);152(3): 503'512' as a fatigue crack initiation paramete4Engineering of G. Pluvinage,Application of notch effect in high cycle fatigue, Proceeding the 96 International Vol. 3, Fatigue of metallic and nonConferenceof on Fracture,ICF9, Advancesin FractureResearch, (1997).p. 1239-L25O. Pergamon metallic materialsand structures. G. Pluvinage,Application of notch Fracture Mechanicsto fracture enunatingfrom stressconcentrators, (1997),p. 2I3Sciences, Engineering of Engineering Congress Computational in Advances Computational

157 1 t 58 I t 59 I t 60 I t 61 I 16Z 1 t 63 I

164 I

t 65 I t 66 I 167 I t 68 I t 69 I t 70 I t 7l I

1,72|

t73 |

2r8.

203

l'14 I

H. Adib, A. Gowhari-Anaraki, Fatiguelift estimation based notchfracture mechanics, International on 5 and 9 Annul Mechanical Engineering Conference.Iranian Society of Mechanical Engineering, ISME2O0l;4: 163-169, Qffit1.

t 75 I

S. Bentachfine,G. Pluvinage,G. Gilgert, Z. Azari, D. Bouami, Notch effect in low cycle fatigue, International (1999); Journal Fatigue, of 2l(l):421430. H. Adib' Application of strain energybasedvolumetricapproach fatigue life prediction, Internal Report for oflran UniversityofScienceandTechnology, (1999)Qn Persian). G. R. hwin, Analysis of stresses and strains near the end of a crack traversing a plate, Journalof Applied -364. Mechanics, ol. 24, ( 1957),pp. 361 Y P. Paris,F. Erdogan,Critical analysisof crackpropagationlnws,Trans.ASME, J. BasicEng.,85, pp. 528s34,(1e63). Xiulin Zhang,M. A. Hirt, Fatiguecrackpropagationin sfeels, Engineering Fracture Mechanics, (1983), 18, pp.965-973. C. Bathias,J. P. Bailon (eds.),/-a fatigue des matriauxet des structures,L.es presses I'universitde de ( Montral,Montral,Canada, 198 1). M. G. Yan, The fatigue crack propagation mechanismof aeronautical structural materials and their engineering applicationin Progress Material Science (Editedby C.X.Chi),pp. 113-158, in Scientific Publishing House,Beijing, China,(1986). R. G. Forman,V. E. Kearney,R. M. Engle, Numerical analysis crack propagation in cyclic loaded of structures, Trans.ASME, J. BasicEng.,89, (1967),pp.459464. R' J. Donahue, M. Clark,P. Atanmo,R. Kumble,A. J. McEvily, Crackopeningdisplacement H. and the rate offatigue crackgrowth,International Journal Fracture of (1972). Mechanics, pp.209-215, 8, C. S. Shin, Fatigue crack growth from stress concentratian and fatigue life prediction in notched components,Handbook of fatigue crack propagation in metallic structures,Ed. A. Carpinteri, (1994) ElsevierScience V., pp. 613-652. B. D' Radaj, Z. Zheng, Geometrically nonlinear behaviorof spot welded joints in tensileand com.pressive shearloading,Engineering Fracture Mechanics, 5l(2), (l9g1), pp. 2Bl-294. D. Radaj, S. Zhang, W. Mohrmawt, load stress parametersat the spot welds of various specimcns, Engineering FractureMechanics, 37(5),(1990),pp. 933-951. D. Radaj, Design and analysis fatigue resistant of weldedstructures,Abington Publishing,Cambridge

176 I

177 I t 78 I

I79 I

t80 I t 8l I

I 82 I t 83 I t 84 I

t 85 I t 86 I t 87 I t 88 I t 89 I I 90 I

(leeo).
H. S. Reemsnydet, Modelingof thefatigue resistance single of spot-welded steelsheet,IlW-Doc.In-ggs-gz andXIII- I 469-92, (1992). M. Fujimoto,N. Mori, S. Sakuma,IIWDoc. III-721-82,(lgB2). T. Satoh,H' Abe, K. Nishikawa, Morita, A studyof threedimensionql M. elastic-plasticstress analysis of spot weldedioint, Prcce'ding of the 4 InternationalConference fatigue and fatigue threshold, 15-20 on July,(1990)in HonoluluHawaii. t 9l I Y. Zhang,D. Taylor, Fatigue life prediction of spot weldedcomponents, Proceeding the 7 international of fatiguecongress, ol. 2, Fatigue99, 127 - 1280, 1999). Y ( 5

204

le2|

PingshDong, Xin Sun, Innmiao Zhang,F. Lu, A generalizedspot weld modelfor finite elemcntanalysisof Papers, Doc. No. I 999-0 -3191, ( 1999). auto-body structures, I SAE Technical Y. Nakahara, A. Masao Fujimoto, M. Takahashi, Kawamoto,Nobon Tomioka, Method of fatigue lift (2000). estimationforspot-welded Doc. No. 2000-01-0779, structures,SAE TechnicalPapers, Pietro Salvini, Francesco Vivio, Vincenzo Vullo, A spot weWfinite element for structural modeling, InternatonalJournalof Fatigue,22(8), (2W); 645-656. P. Salvini, F. Vivo, V. Vullo, Use of an appropriatespot weld element fatigue life preventionof for joints, European (Spain), (1999). complex AutomotiveCongress Barccelona 99, P. Salvini,F. Vivio, V. Vullo, La resistenzg afatica di structurein acci.ano innossidiable saladate punti : a sintesidelleanalisi teoriche, numeriche sperimentali, e PartsI, II, m. ATA (1998);51(6-7):338-346;15(89); 438449; I 5(I 0):53G43.

t e 3I l e 4l t e sl t e 6l

l.e7 l t e 8l t e el
tl00l

joints, P. Salvini,E. Sccardeccia, Demofonti,A procedure fatigue life prediction of spot weWed G. for (1997);20(8):1117-1128. Fatigue Fracture Engineering and of MaterialStructures, Y. Kan,Fatigueresistance spotwelds-an analyticalstudy,Met.Eng.Quart(1976); 16(4):26-36. of D. Radaj, Theory of forces and stresses spot welded overlapjoints, Archive of Applied Mechanics67 in (1996) v2,22-34. H-F. Henryson,F. Abdulwahab, L. Josefson, Fermr,Residualstresses resistance B. M. in spot welds: Finite elementsimulations,X-Ray measurements influenceon fatigue beluvior, VTT-SYMP-182, and (1998), 179-190. Espoo-Finland, PP.

tr0l l

H-F Henryson, Abdulwahab, Josefson, Feremer, F. BL M. Residualstresles resistance in spot welds: finite elementsimul.ation,X-Ray measurementland influence on Fatigue behaviaur, Welding in the World (1999); 43(1):55-63. H-F Hemyson, B. L. Josefson,Fatigue crack initiation at spot welds, Fracture Mechanics or Strain life approach,Proceeding ofthe 7 international fatiguecongress, Fatigue99,(1999);2:L251-1256 X. Deng,W. Chen,G. Shi, Threedimensional analysisof the mechanical behaviourof spot finite element (2W);35:17 -39. welds, FitttteElementAnalysisandDesign, W. Chen,W. Deng, Performance shell elements modellingspot weldedjol'nts,Finite ElementAnalysis of in (200O);35:41-57. andDesign, N. Pan S. Sheppard, Spot welds life prediction with cyclic strain range, InternationalJoumal of Fatigue, (20U2) 24(5):5 L9-528. ; N. Pan, S. Sheppard, ,Stressintensity faclors in the spot welds, Engineering Fracture Mechanics, (2003);70(5): r-684. 67 H. Adib, J. Gilgert, G. Pluvinage,Fatigue life prediction durationfor spot weldsby volwnetic approach, -rO. InternationalJoumalof Fatigue,(2(rf,/l1;26117.9, D. J. Chang,R. Muki, Stress distrtbution in a lap joint undcr tensionshear,IntenrationalJournalof Solid (197a); Structures, l0(5):503-5 17. L. P. Pook,Approximatestress intensity factorsfor spot and similnr wel, NationalEngineeringLaboratory UK, (1975);NEL rc,port 588. J.F Cooper,R. A. Smith, Initial fatigue crack growth at spot welds, proceding of internationalconference on fatigueof engineering materials stuctures and ,Yol.2, (1986),pp.283-288.

tl02 l t 1 0 3l tl04l

t l o sl
l 1106 ll07 l t 1 0 8l tl09l lll0l

205

t r 1 1l ltrz l l1l3 l til4 l tlls l t 1 1 6l

S. Zhang,Stress intensities the spots,International at Journal ofFracture,(1997\;88:167-185. MH. Swellam,P. Kurath, FV. Lawrence,Electric-Potentinl-drop studiesof fatigue crack development in tensileshearspotwer.ds, ASTM STP(1991);ll22:383-401. MH' Swallam,FV. Lawrence,A fatigue designparameterlor spot welds,FractureControl ProgramReport No. 157,Universityof Illinois at Urbana-Champaign, (1gl). MH. Swallam, PhD Thesis, Universityof Illinois at Urbana-Champaign, (1991) F. V. Lawrence,P. C. Wang, H. T. Corten, An empirical methodfor estimatingthe fatigue resistanceof tensileshearspotwelds,SAIE paper830035, (1983). H. Kang, M' E. Barkey, Y. Ie, Evaluation of multiaxial spot weld fatigue parametersfor proportional loading,lnternational Journal Fatigue, of 22(8);(200O):691-702. SD. Sheppard,Estimation of fatigue propagation life in resistancespot welds, Advances in fatigue philadelphia, predictive techni4ues,2d Vol., ASTM STp 1211, (1993):169-1g5. SD. Sheppard,FurTher refinementof a methodologyfor fafigue life estimation in resistancespot weld connections,Advances in fatigue predictive techniques,3d Vol., ASTM STp 1292, philadelphia (1996):265-282. H. Adib, Effect of the geometryin tensile shear spot wetdsfor fatigue tife prediction,Internal report of Laboratoire FiabilitMcanique, July 2000). de (28 H. Adib, J. Gilgert, G. Pluvinage,Tensile-Shear spot welded joints fatigue life estimationbasedon hot volumetric approach,New Trendsin FatigueandFracture NT2F2Congress, IIHamammet-Tunesia, 12-13 May, (2003). ANSYS Inc. GettingStarted manual. ANSYS Inc. BasicAnalysisProcedure Guidemanual. ANSYS Inc. AdvanceAnalysisTechnical Guidemanual. H' Adib, J. Jeong, G. Pluvinage, Three dimensionalfinite element analysisof hot spot welded joints in tensile and compressiveloadings, Journal of "strength of Material", Accepted for Publication, (August 2OO3). F. M. Haggag,G. E. C. Bell, Measurement Yield strengthand flow properties in spot welds and their of HA7'sat varoius strain Rates,Presented the ASM, 3th InternationalConference the Trendsin Welding at on Research, June,1992Gatinburg,Tennessee, l-5 USA, (1992). F' M. Haggag,In situ Measurements Yield Strength andflowproperties Using a Novel ABI system, of Presented the ASTM Symposiumon small specimentest techniques their applicationsto the nuclear at and reactorvessel Thermalannealing plant life extension, and New Orleans, Lousiana, January1992,published at ASTM STP,(1992).

[rr7) tlls l

tile l [r20 )

trzr I
t r 2 2I lt23l I r 2 4l

I r z sl

t r 2 6l

206