Exercice 57 :

1) Détermination la répartition de la contrainte tangentielle :

¿
TS x
τ ( y )=
On a : bI x
1
I x = παβ 3
avec 4
S x=∫ y dA
¿ ¿ ¿
¿
Et A
¿ ¿
dA =a∗dy
L’equation de l’ellipse est :

() ()
x 2 y 2
α
+
β
=1

√β
=± −( ) +1
√()
2
x y y 2
x=±α − +1
α ; β
a
¿
Donc pour une ordonné y on a l’abscisse 2 par conséquent on peut écrire :

a=2 α −
y 2
β
+1
√()
Alors , substituant l’equation () dans () ensuite dans () on aura :

S x=∫ 2αy (−
¿

A
y 2
β
¿
+1)dy
¿ ¿

√()
{ ( )}
β 2 1
y¿
S x=∫ 2αy 1 −
¿ ¿ 2 ¿
dy
y β

( )
¿
y 2 2
1− =u
On pose : β ;
¿ ¿
y dy =−β 2 udu
Pour
Exercice 60 :

1. La valeur de π :
π= pi=3 . 1415926535897932384626433832 D−00
2. L’aire totale de la section :
π
A= a2 =0 .392699081 a 2
8
3. Le centre e gravité de la section :
128 √ 2
y G= a=0 . 54876423 a
105 π
IX
4. Le moment d’inertie 2

45 π +128 2
IX = −( y G ) A
2 1920
4
I X2=0. 022039486 a
5. Détermination la variation de la contrainte
¿
TS x ( y 2 )
τ ( y 2 )= 2

I X b ( y2)
On a : 2

 L’intervalle de variation de l’ordonnée

y2 :

On a :
− y G≤ y 2≤a− y G
yG y2 yG
− ≤ ≤1−
a a a
θ θ y
 Les angles 1 , 2 correspondant à l’ordonnée 2 :
π π
≤θ 1≤
Soit : 4 2
π 3π
≤θ 2≤
2 4
On a :
y 2 + y G=r sin θ1
=(−a cos2 θ1 )sin θ 1
On obtient l’equation suivante :
y2 yG
f (θ1 )=sin θ1 cos 2θ 1 + + =0
a a

Pour
y2 π π
,
[ ]
a donnée dans l’intervalle 4 2 qui s’annulle l’equation suivante
 y2 yG
f (θ1 )=sin θ1 cos 2θ 1 + + =0
a a
On utilise la methode de bipartition “dichotomie “ par voie numerique

π
θ a=
4 ; f (θa )≥0

π
θb =
2 ; f (θb )≤0 avec f (θa ). f (θb )≤0

∃θ1 ∈ [ θ a , θb ]
; f (θ1 )=0

En suite on determine
θ2 = pi−θ1

 Determination la fonction b ( y 2 ) :
On a d’apres la figure :
b( y 2 )
=r cos θ1 =(−a cos 2θ 1 )cos θ 1
2
b( y 2 )
=−2 cos 2 θ1 . cos θ1
Alors : a
¿
Sx ( y2 )
 Determination 2 :
¿
Sx ( y2 ) ¿
 On determine le moment statique 2 numeriquement on discretise la section supperieur A en
des triangles.
a− y G
¿
S x ( y 2 )=
2
∫ y 2 dA
¿

y2
i=NT i=NT
¿ 2 ∑ ΔS x (i)=2 ∑ AT (i). y GT (i )
¿
2
i=1 i=1

Avec NT :le nombre de triangle donnée dans le premier quadrant.

AT (i ) :surface du triangle (i) i=1,2,3,....NT
y GT (i) :coordonnée y G de centre de gravité du triangle (i)
i=1,2,3,....NT

a) Pour un triangle(i) on a :
1
AT (i )= {(x 2 −x 1 )( y 3− y1 )−( y 2 − y 1 )(x 3 −x 1 ) }
2
1
y GT (i)= { y 1 + y 2 + y 3 }
3
b) Position des nœuds dans le repere
Gx 2 y 2 :
Le nombre totale de nœuds est determinée par :
NN =NT +2
Par raison de symetrie par rapport a l’axe 2 on discrétise Gy
Uniquement le premier quadrant ensuite on multipliant le resultat par deux.
On choisi la distribution suivante :
x1 y1 y2
=0 =
a ; a a
π
( )−θ1
2
Δθ= =0
NT

θ=θ1 +(i−2) Δθ
R
=−cos(2θ )
a i=2,3,4....NN
xi
a
=
R
a ( )
cos θ

yi R y
a
=
a ( )
sin θ− G
a
c) numérotation des nœuds de chaque triangle :
les nœuds de chaque triangle (i) ; sont numérotés dans le sens antihoraire par :
i , i +1 ,i +2 (i=1 ,2,3,.... NT )
Finalement :
¿
Sx
( )( )
i=NT
AT ( i ) y GT (i )
a
3
2
=2 ∑ a2 a
i=1

Avec : a 2
=
2 a
−
a {[
AT (i) 1 x (i+1) x (1) y (i+2 ) y (1)
a
−
a
−
a a a
−][
y (i+1 ) y (1) x(i+2 ) x (1)
−
a ][ ][ ]}
{ }
y GT 1 y (1 ) y (i+1 ) y (i+ 2)
i
= + +
Et a 3 a a a

{( ) }
S ¿x
2
3
T a
τ ( y 2 )=
a2
( )
Ix b( y 2 )
2

a 4
a
D’où :

Programme :

Exercice 61 :

1) Détermination l’aire totale de la section :

[∫ ]
c y1 c πx1
h cos
A=∫ dA=∬ dx 1 dy 1 =∫ dy 1 dx 1 =∫ [ y 1 ]0 2c
dx 1
A A −c 0 −c

c c
πx1 πx
=∫ h cos dx 1 =2 h∫ cos 1 dx 1
−c 2c 0 2c (Par raison de symétrie)
=2 h
2c
π
sin
2c ( )
πx 1 c 4 ch
]0 =
π
π 4
sin −0 = ch
2 π { }
4
A= ch=1 . 273239545 ch
Alors π
2) Determination la position du centre de gravité de la section :
x G=0 (par raison de symetrie par rapport à l’axe
oy 1 .
Sx
y G= 1

[∫ ]
y1

[ ]
c c πx
1 2 h cos 2 c1
S x =∫A y 1 dA=∬A y 1 dx 1 dy 1 =∫ y 1 dy 1 dx 1 =∫ y
1
−c 0 −c
2 10
Avec :

c c
1 2 πx 1 πx
= h ∫ cos 2
dx 1 =h ∫ cos 2 1 dx 1
2
2 −c 2c 0 2c (par raison de symétrie).

( πx 1 c 1 2
) { }
c
1 πx 1 1 2 c
=h ∫ 1+cos
2
dx 1 = h x1 + sin = h c
0
2 c 2 π c 0 2
1 2
Sx = h c
Alors : 1 2
1 π π π
y G = h2 c . = h y G= h
Donc : 2 4 ch 8 ; 8

Ix
3) Determination le moment d’inertie 1 :

[ ]
πx 1
y 1=h cos
c 2c

I x =∫A y dA=∬A y dxdy =∫ ∫

2 2 2
1 1 y 1 dy 1 dx 1
1
−c 0
πx 1

{ }
c c
1 3 R cos 1 3 3 πx 1
=∫ y 2 c dx
1 = h (2)∫ cos dx
−c 3 1 0 3 0 2c 1
πx 1 2c 2c
=z x 1= z dx 1 = dz
On pose : 2 c ; π ; π
x 1=0 ; z=0
π
x 1=c ; z= 2
π π
2 2

Alors :
1
2
I x = h3∫ cos 3 z
3 0
2c
π
4
dz= ch3∫ cos 3 zdz
3π 0
( )
π π
2 2
4 3 4
ch ∫ cos z cos zdz= ch ∫ ( 1−sin z ) d ( sin z )
2 3 2
=
3π 0 3 π 0
 4 3
3{
1 3 π2 4 3
0 3π
1 8 3
= ch sin z− sin z = ch 1− = ch
3π 3 3π } { }
8 3
Ix = ch
Alors : 1 9π

Ix
4) Détermination le moment d’inertie 2 :
2
I x 1=I x 2 + y G A

Alors :

128−9 π 3 3
Ix = ch =0 . 08659258 ch
On trouve : 2 144 π

TS ¿x ( y 2 )
τ ( y 2 )= 2

I x2 b ( y 2)
5) Détermination la variation de la contrainte
 Intervalle de variation de l’ordonnée
y2 :

On a :
− y G≤ y 2≤ yG
yG y2 y
− ≤ ≤1− G
h h h
b ( y 2)
 Détermination la fonction :
On a d’après la figure :
π b
y 2 + y G=h cos
4c 2
π y y
f ( b )=cos b− 2 − G =0
4c h h
y 2 yG
cos
π b
2c 2
= +( )
h h ;
π
4c
y y
b=arccos 2 + G
h h ( )
( )
4 ( y2) 4 y y
= arccos 2 + G
c π h h
¿
Sx ( y2 )
 Détermination 2 :
¿
Sx ( y2 )
On determine le moment statique 2 numeriquement on discretise la section
¿
supperieur A
en des triangles.
c− y G
i=NT i=NT
Sx =
¿
2
∫ y 2 dA =2 ∑ Δ S x 2 (i)=2 ∑ AT (i ). y GT (i)
¿ ¿

y2 i=1 i=1

Avec : NT :le nombre de triangle donnée dans le premier quadrant.

AT (i ) :surface du triangle (i) i=1,2,3,....NT
y GT (i) :coordonnée y G de centre de gravité du triangle (i) i=1,2,3,....NT
a) Pour un triangle(i) on a :
1
AT (i )= {(x 2 −x 1 )( y 3− y1 )−( y 2 − y 1 )(x 3 −x 1 ) }
2
1
y GT (i)= { y 1 + y 2 + y 3 }
3
b) Position des nœuds dans le repere
Gx 2 y 2 :
Le nombre totale de nœuds est determinée par :
NN =NT +2
Par raison de symetrie par rapport a l’axe 2 on discrétise Gy
Uniquement le premier quadrant ensuite on multipliant le résultat par deux.

On choisi la distribution suivante :

x1 y1 y2
=0 =
a ; a a

Δx =
b
2 ()
NT

xi b
= −( i−2 ) Δx
h 2
yi
h
=cos
π xi
2 c ( ) i=2,3,4....NN
c) numérotation des nœuds de chaque triangle :
les nœuds de chaque triangle (i) ; sont numérotés dans le sens antihoraire par :
i , i +1 ,i +2 (i=1 ,2,3,.... NT )
Finalement :
¿
Sx
( ATch(i) )( y )
i=NT
T (i)
ch
2
2
=2 ∑ GT
h
i=1

Avec : a 2
=
2 c
−
c {[
AT (i) 1 x (i+1) x (1) y (i+2 ) y (1)
h
−
h
−
h h ][ c
−
c ][
y (i+1 ) y (1) x(i+2 ) x (1)
− ][ ]}
{ }
y GT 1 y (1 ) y (i+1 ) y (i+ 2)
i
= + +
Et c 3 c c c

{( ) }
S¿x
2

T ch2
τ ( y 2 )=
ch
( b(cy ) )
Ix
2 2
3
ch
D’où :
Exercice 62 :
1. L’aire totale de la section :
π
A= a2 =0 .392699081 a 2
8
2. Le centre e gravité de la section :
128 √ 2
x G= a=0 .54876423 a
105 π
IX
3. Le moment d’inertie 2

45 π +128 4
IX = a
2 1920
4
I X2=0. 006964411a
4. Détermination la variation de la contrainte
TS¿x ( y 2 )
τ ( y 2 )= 2

I X2 b ( y 2 )
On a :
Par raison de symetrie par rapport à l’axe
ox 2 on a :
τ ( y 2 )=τ (− y 2 )
, pour cela on s’interesse uniquement pour
y 2 ≥0

 L’intervalle de variation de l’ordonnée

y2 :

On a :
x 1=r cos θ
y 1 =r sin θ
Avec : r =a cos 2 θ sin θ

dy 1
=0
dθ −2 cos 2 θ sin θ+cos 2 θ cos θ=0
;
2 tan 2 θ tan θ=1

2
( 2 tanθ
1−tan 2 θ )
tan θ=1
;
( 4 tan 2 θ
1−tan2 θ
=1
)
1
2 2 2
tan 2 θ=
4 tan θ=1−tan θ ; 5 tan θ=1 ; 5 ;
°
θ=24 . 09484255 =0 . 420534335 rad
1 1 5
cos 2 θ= = =
1+ tan θ 1+ 1 6
2

sin θ=
tan2 θ
2
2
=
1 5 1
1+ tan θ 5 6 6
= ( )( )
5 1 4
cos2 θ=cos 2 θ−sin 2 θ= − =
6 6 6

Alors :
y max =a
4 1
6 √6
4 6
6 6 9
6
( )( )
=a √ =a √ =0. 272165527 a

6 y
− √ ≤ 2≤ √
6
Alors : 9 a 9
y2
-0 . 272165527≤ ≤0 .272165527
a
 Determination les valeurs de
θ1 etθ2 :

On a
0≤θ1 ≤0 . 420534335 rad
 0 . 420534335 rad≤θ 1≤0. 785398163 rad
y 2 =r sin θ1
=(a cos2 θ1 )sinθ 1
=a cos2 θ1 sin θ1
y2
=cos 2 θ1 sin θ 1 θ1 ∈ [ 0 , θm ] θm =0 . 420534335 rad
a ; ;

De même :
y2
a
=cos 2 θ2 sinθ 2
;
[ ]
θ2 ∈ θm ,
π
4 θ =0 . 420534335 rad
; m

y2
Pour a donné dans l’intervalle
0,
√6
[ ]
9 on calcule θ1 dans l’intervalle [ 0 , θm ] et θ2 dans l’intervalle

[ ]
θm ,
π
4 qui annullent respectivement
y2
f (θ1 )= −cos 2 θ1 sin θ 1=0
a
y
g(θ2 )= 2 −cos 2 θ2 sin θ2 =0
a
Pour cela on utilise la méthode de bipartition “dichotomie “ par voie numérique
a) Determination
θ1 :
θ a=0 ; f (θa )≥0
θb =θ m ; f (θb )≤0 avec
f (θa ). f (θb )≤0
∃θ1 ∈ [ θ a , θb ]
pour que f (θ1 )=0
b) Determination
θ2 :
θ a=θ m ;
g(θ a )≤0
π
θb =
4 ; g(θ b )≥0 avec
g(θ a ). g (θb )≤0
∃θ2 ∈ [ θa , θb ]
pour que g(θ2 )=0

 Determination la fonction b ( y 2 ) :
b ( y 2 )=r 1 cos θ 1−r 2 cos θ2
=(a cos2 θ1 )cos θ1 −( a cos2θ 2 )cos θ 2
b( y 2 )
=cos 2 θ1 . cos θ1 −cos 2 θ2 . cos θ 2
Alors : a
S¿x ( y 2 )
 Determination 2 :
A¿ en
¿
Sx ( y2 )
On determine le moment statique 2 numeriquement on discretise la section supperieur
des triangles.
√6 a
9

S x2 ( y 2 )= ∫ y 2 dA
¿ ¿

y2
 i=NT i= NT
¿ 2 ∑ ΔS x (i)= ∑ AT (i). yGT (i)
¿
2
i=1 i =1

Avec NT :le nombre de triangle donnée .

AT (i ) :surface du triangle (i) i=1,2,3,....NT
y GT (i) :coordonnée y G de centre du gravité du triangle (i)
i=1,2,3,....NT

a) Pour un triangle(i) on a :
1
AT (i )= {(x 2 −x 1 )( y 3− y1 )−( y 2 − y 1 )(x 3 −x 1 ) }
2
1
y GT (i)= { y 1 + y 2 + y 3 }
3
b) Position des nœuds dans le repere
Gx 2 y 2 :
Le nombre totale de nœuds NN est determinée par :
NN =NT +2
On choisi la discritisation suivante :
Pour obtenir les positions des nœuds :

x1
=
( b ( y 2)
a ) + cos 2 θ2 cos θ 2−
xG
a 2 a b ( y 2)
(au mlieu de la largeur )
y1 y2
=
; a a
θ −θ
Δθ= 2 1 =0
NT

θ=θ1 +(i−2) Δθ
R
=−cos(2θ )
a i=2,3,4....NN
xi x
a
=
R
a ( )
cos θ− G
a
yi
a
=
R
a ( )
sin θ

c) numérotation des nœuds de chaque triangle :

les nœuds de chaque triangle (i) ; sont numérotés dans le sens antihoraire par :
i , i +1 ,i +2 (i=1 ,2,3,.... NT )
Finalement :
 ¿
Sx
( () y a( i ) )
i= NT
AT ( i )
=∑
2 GT
3 2
a i =1 a

Avec : a 2
=
2 a
−
a {[
AT (i) 1 x (i+1) x (1) y (i+2 ) y (1)
a
−
a
−
a a ][
y (i+1 ) y (1) x(i+2 ) x (1)
−
a
−
a ][ ][ ]}
Et a 3 a {
y GT (i ) 1 y (1) y (i+1 ) y (i+2 )
= +
a
+
a }

{( ) }
S ¿x
2
3
T a
τ ( y 2 )=
a2
( b(ay ) )
Ix
2 2
4
a
D’où :

Programme :
Exercice 63 :
1) Determination l’aire totale de la section :
π
A=∫ dA=∬ dx dy =∫ A y dx=−4∫02 ( a sin3 t ) { 3 a cos 2 t } (−sin t ) dt
A A
π π
=12 a ∫
2 2
0
4 2
sin t cos tdt =12 a 2
∫ ( sin 4 t−sin6 t ) dt
2
0
Connaissant que :
1
sin 4 t= ( cos 4 t−4 cos2 t+3 )
8
1
sin6 t= (−cos6 t+6 cos4 t−15 cos2 t+10 )
Et 32
Alors :

[( )]
π
A=12 a
1 sin 4 t sin 2 t
8 4
−42
2
1 sin 6 t sin 4 t
+3 t − −
32 6
+6
4
−15
sin 2 t
2
+10 t ) ( 0
2

=12 a 2 { ( ) ( )}
3 π 10 π
−
8 2 32 2
3
= πa 2
8
3π 2
A= a
8

2) Determination le moment d’inertie

Ix :
1 3
I x =∫A y dA=∬A y dxdy =∫ [∫ y dy ] dx=∫A y dx
2 2 2
3
π
1 0
= 4 ∫π ( a sin 3 t ) { 3 a cos 2 t } (−sin t ) dt=4 a4 ∫02 sin 9 t cos 2 t sin tdt
3
3 2
π π
=4 a 4
∫ 2
0
10
sin t cos tdt=4 a2 4
∫ {sin10 t−sin12 t } dt
2
0
π
=4 a 4 ∫ 2
0 {63256 −231
1024 } dt=
21
4( ) a
1024 2
π 4
 21 4 4
I x= πa =0 . 128854386 a
512
3) Détermination la variation de la contrainte
TS ¿x ( y 2 )
τ ( y 2 )= 2

x2 2 I b( y )
On a :
Par raison de symetrie par rapport a l’axe ox on a :
τ ( y )=τ (− y ) ; pour cela on s’interesse uniquement pour y≥0
 Intervalle de variation de l’ordonnée
y2 :
π
t=
2 ; x=0
y= a alors
−a≤ y≤a
y
−1≤ ≤1
a
 Détermination la fonction b ( y ) :
Pour b ( y ) donnée dans l’intervalle [ 0,1 ] (‘’moitie supplémentaire’’)
b( y)
=a cos 3 t
2
y=a sin3 t

()
1
y y
=sin3 t sin t= 3
a , a
Connaissant que :

√
2
cost= √1−sin t= 1−
y
a
2
() 3
;
0≤t≤
π
2
Alors :

[ ( )]
2 3
b( y) y 3 2
=a 1−
2 a

[ ( )]
2 3
b( y) y 3 2
=a 1− y
Donc :
2 a 0≤ ≤1
a
¿
Sx ( y2 )
 Détermination 2 :
¿
Sx ( y2 )
On determine le moment statique 2 numeriquement on discretise la section
¿
supperieur A
en des triangles.
a− y i=NT i=NT
S x = ∫ ydA =2 ∑ ΔS (i)=2 ∑ AT (i). y GT (i )
¿ ¿ ¿
x
2
y i=1 i=1
Avec : NT :le nombre de triangle donnée dans le premier quadrant.
 AT (i ) :surface du triangle (i) i=1,2,3,....NT
y GT (i) :coordonnée y G de centre de gravité du triangle (i) i=1,2,3,....NT

d) Pour un triangle(i) on a :
1
AT (i )= {(x 2 −x 1 )( y 3− y1 )−( y 2 − y 1 )(x 3 −x 1 ) }
2
1
y GT (i)= { y 1 + y 2 + y 3 }
3
e) Position des nœuds dans le repère
Gx 2 y 2 :
Le nombre totale de nœuds est determinée par :
NN =NT +2
Par raison de symetrie par rapport a l’axe on discrétise oy
Uniquement le premier quadrant ensuite on multipliant le resultat par deux.

On choisi la distribution suivante :

x1 y1 y
=0 =
a ; a a

Δx =
( a )
b( y )

2 NT

xi
=
( a )
b( y )

− (i −2 ) Δx
a 2

[ ( )]
2 3
yi xi 3 2
= 1−
a a i=2,3,4....NN
f) numérotation des nœuds de chaque triangle :
les nœuds de chaque triangle (i) ; sont numérotés dans le sens antihoraire par :
i , i +1 ,i +2 (i=1 , 2, 3, .. .. NT )
Finalement :
¿
Sx
( )( )
i=NT
AT (i) y G T ( i)
a
3
2
=2 ∑ a
2
a
i=1

Avec : a 2
=
2 a
−
a {[
AT (i) 1 x (i+1) x (1) y (i+2 ) y (1)
a
−
a ][
−
a a ][
y (i+1 ) y (1) x(i+2 ) x (1)
−
a
−
a ][ ]}
{ }
y GT 1 y (1 ) y (i+1 ) y (i+ 2)
i
= + +
Et a 3 a a a
 {( ) }
S¿x
2
3
T a
τ ( y 2 )=
a2 Ix

D’où :
a
2
4 ( b(ay ) )
Exercice 64 :
1. L’aire totale de la section :
A=3 πa 2
2. Le centre e gravité de la section :
5
x G=πa y G= a
; 6
IX
3. Le moment d’inertie 2 :
5 4
I X = πa
2 6
π
I Y = ( 12 π −35 ) a
2 4
2 12
4. Détermination la variation de la contrainte
¿
TS x ( y 2 )
τ ( y 2 )= 2

I X b ( y2)
On a : 2

 L’intervalle de variation de l’ordonnée

y2 :

On a :
− y G≤ y 2≤2 a− y G
yG y2 y
− ≤ ≤2− G
Alors : a a a
5 y 7
− ≤ 2≤
Ou : 6 a 6
 Determination la fonction b ( y 2 ) :
y2
On a pour a donne:
{¿
y2 yG y 5 1 y
cos t=1− − =1− 2 − = − 2
Alors : a a a 6 6 a

t =arccos ( 1 y2
−
6 a )
et :
b( y 2 )
a {
=2 −t +sin t +
xG
a }
=2 {−t +sin t + π }

alors : {¿ 0≤t 1 ≤π
S¿x ( y 2 )
 Determination 2 :
A¿ en
¿
Sx ( y2 )
On determine le moment statique 2 numeriquement on discretise la section supperieur
des triangles.
Par raison de symetrie par rapport a l’axe
Gy 2 on discrétise
 ¿
Sx ( y2 )
Uniquement le premier quadrant ensuite on multipliant le resultat de 2 par deux.

2 a− y G
¿
S x ( y 2 )=
2
∫ y 2 dA
¿

y2
i=NT i=NT
¿ 2 ∑ ΔS x (i)=2 ∑ AT (i). y GT (i )
¿
2
i=1 i=1

Avec NT :le nombre de triangle donnée dans le premier quadrant.

AT (i ) :surface du triangle (i) i=1,2,3,....NT
y GT (i) :coordonnée y G de centre de gravité du triangle (i)
i=1,2,3,....NT

d) Pour un triangle(i) on a :
1
AT (i )= {(x 2 −x 1 )( y 3− y1 )−( y 2 − y 1 )(x 3 −x 1 ) }
2
1
y GT (i)= { y 1 + y 2 + y 3 }
3
e) Position des nœuds dans le repere
Gx 2 y 2 :
Le nombre totale de nœuds est determinée par :
NN =NT +2
On choisi la distribution suivante :
x1 y1 y2
=0 =
a ; a a
t 2=2 π−t 1
t 2 −π
Δt =
NT
¿ i=2,3,4....NN

c) numérotation des nœuds de chaque triangle :

les nœuds de chaque triangle (i) ; sont numérotés dans le sens antihoraire par :
i , i +1 ,i +2 (i=1 , 2, 3, .. .. NT )
Finalement :
¿
Sx
( )( )
i=NT
AT ( i ) y GT (i )
a 3
2
=2 ∑ a2 a
i=1

Avec : a 2
=
2 a
−
a {[
a
−
a ][
AT (i) 1 x (i+1) x (1) y (i+2 ) y (1)
−
a a ][
y (i+1 ) y (1) x(i+2 ) x (1)
−
a
−
a ][ ]}
 { }
y GT 1 y (1 ) y (i+1 ) y (i+ 2)
i
= + +
Et a 3 a a a

{( ) }
S ¿x
2
3
T a
τ ( y 2 )=
a2
( b(ay ) )
Ix
2 2
4
a
D’où :

Exercice 65 :
1. L’aire totale de la section :
A=3 πa 2
2. Le centre e gravité de la section :
5
x G= a y =πa
6 ; G
IX
3. Le moment d’inertie 2 :
π
IX = ( 12 π 2−35 ) a4
2 12
4. Détermination la variation de la contrainte
¿
TS x ( y 2 )
τ ( y 2 )= 2

I X b ( y2)
On a : 2

 L’intervalle de variation de l’ordonnée

y2 :

On a :
− y G≤ y 2≤2 πa− yG
y2
−πa≤ ≤πa
a
y2
−π≤ ≤π
a
 Determination la valeur de
t 2 corespondant a y 2 :
y 2 + y G=a ( t 2 −sin t 2 ) π≤t 2 ≤2 π
;
y2 yG
+ =t 2 −sin t 2 π≤t 2 ≤2 π
a a ;
y2 yG
f ( t 2 )=t 2 −sin t 2− =0 −
a a
y
f ( t 2 )=t 2 −sin t 2− 2 −π=0
a
Pour cela on utilise la methode de dipartition par voie numerique :
t a =π ;
f a ≤0
t b =2 π ;
f b ≥0 avec
f a ( t a ) . f b ( t b ) ≤0
∃t 2 ∈ [ t a , t b ] f ( t 2 )=0
pour que

 Determination la fonction b ( y 2 ) :
On a :
 b ( y 2 )=x1 ( t 2 )
=a ( 1−cos t 2 )

b( y 2 )
=1−cos t 2
a
¿
Sx ( y2 )
 Determination 2 :
S¿x ( y 2 ) ¿
On determine le moment statique 2 numeriquement on discretise la section supperieur A en
des triangles.
Par raison de symetrie par rapport a l’axe
Gy 2 on discrétise
¿
Sx ( y2 )
Uniquement le premier quadrant ensuite on multipliant le resultat de 2 par deux.

πa
S x ( y 2 )=∫ y 2 dA
¿ ¿
2
y2
 i=NT i=NT
¿ ∑ ΔS x (i )= ∑ AT (i ). y GT (i)
¿
2
i=1 i=1

Avec NT :le nombre de triangle donnée dans le premier quadrant.

AT (i ) :surface du triangle (i) i=1,2,3,....NT
y GT (i) :coordonnée y G de centre de gravité du triangle (i)
i=1,2,3,....NT

a) Pour un triangle(i) on a :
1
AT (i )= {(x 2 −x 1 )( y 3− y1 )−( y 2 − y 1 )(x 3 −x 1 ) }
2
1
y GT (i)= { y 1 + y 2 + y 3 }
3
b) Position des nœuds dans le repere
Gx 2 y 2 :
Le nombre totale de nœuds dans la discretisation est determinée par :
NN =NT +2
On choisi la distribution suivante :
x1 5 y1 y2
=−X G =− =
a 6 ; a a
2 π−t 2
Δt =
NT
¿ i=2,3,4....NN

c) numérotation des nœuds de chaque triangle :

les nœuds de chaque triangle (i) ; sont numérotés dans le sens antihoraire par :
i , i +1 ,i +2 (i=1 , 2, 3, .. .. NT )
Finalement :
¿
S x i= NT AT ( i ) y ( i )
a
3
2
=∑
i =1 a
2 ( GT
a )( )
Avec : a 2
=
2 a
−
a {[
AT (i) 1 x (i+1) x (1) y (i+2 ) y (1)
a
−
a
−][ a a ][
y (i+1 ) y (1) x(i+2 ) x (1)
−
a
−
a ][ ]}
{ }
y GT 1 y (1 ) y (i+1 ) y (i+ 2)
i
= + +
Et a 3 a a a
 {( ) }
S ¿x
2
3
T a
τ ( y 2 )=
a2
( b(ay ) )
Ix
2 2
4
a
D’où :
Exercice 66 :
1. L’aire totale de la section :
3π 2
A= a
2
2. Le centre e gravité de la section :
5
x G= a
6
IX
3. Le moment d’inertie 2

21 4
IX = πa
2 32
4. Détermination la variation de la contrainte
TS¿x ( y 2 )
τ ( y 2 )= 2

I X2 b ( y 2 )
On a :
Par raison de symetrie par rapport à l’axe
ox 2 on a :
τ ( y 2 )=τ (− y 2 )
, pour cela on s’interesse uniquement pour
y 2 ≥0

 L’intervalle de variation de l’ordonnée

y2 :

On a :
x 1=r cos θ
y 1 =r sin θ
Avec : r =a ( 1+cosθ )
Alors :
y 1 =a ( 1+cos θ ) sin θ
dy 1
=0
dθ ; −sin θ sin θ+ ( 1+cos θ ) cosθ=0
−sin 2 θ +cos θ+ cos2 θ=0
2 cos 2 θ +cos θ−1=0
θ=π
π π
θ= =60 ° θm =60 °=
3 ; 3
Pour θ=π ;
y 1 =0

θ=
3 ;
1
2 2 ( )
π y =a 1+ 1 √3 = 3 √ 3 a
4
3 3 y 3 3
− √ ≤ 2≤ √
Alors : 4 a 4
y2
-0 . 272165527≤ ≤0 .272165527
a
 Determination les valeurs de
θ1 etθ2 :
 π
0≤θ1 ≤
On a 3
π
≤θ 2≤π
3
y 2 =r sin θ1
=a(1+cosθ 1 )sin θ1
y2
a
=(1+ cos θ1 )sin θ 1
;
[ ]
θ1 ∈ 0 ,
π
3

De même :
y2
a
=(1+ cos θ2 )sinθ 2
;
θ2 ∈
[ ]
π
3
,π

y2
Pour a donné dans l’intervalle
0,
3 √3
[ ]
4 on calcule θ1 dans l’intervalle
0,
π
[ ]
3 et θ2 dans l’intervalle

[ ]
π
3
,π
qui annullent respectivement
y2
f (θ1 )= −(1+ cos θ1 )sin θ 1=0
a
y
g(θ2 )= 2 −(1+cos θ 2 )sin θ2 =0
a
Pour cela on utilise la méthode de bipartition “dichotomie “ par voie numérique
c) Determination
θ1 :
θ a=0 ; f (θa )≥0
π
θb =
3 ; f (θb )≤0 avec
f (θa ). f (θb )≤0
∃θ1 ∈ [ θ a , θb ]
pour que f (θ1 )=0
d) Determination
θ2 :
π
θ a=
3 ; g(θ a )≤0
θb =π ; g(θ b )≥0 avec
g(θ a ). g (θb )≤0
∃θ2 ∈ [ θa , θb ]
Pour que g(θ2 )=0

 Détermination la fonction b ( y 2 ) :
b ( y 2 )=r 1 cos θ 1−r 2 cos θ2
=a(1+cosθ 1 )cosθ1 −a (1+cosθ2 )cosθ2
b( y 2 )
=(1+cos θ1 )cos θ 1−(1+cos θ 2 )cos θ 2
Alors : a
¿
Sx ( y2 )
 Determination 2 :
¿
Sx ( y2 ) ¿
On determine le moment statique 2 numeriquement on discretise la section supperieur A en
des triangles.
 3 √3
4

S x2 ( y 2 )= ∫ y 2 dA
¿ ¿

y2
i=NT i=NT
¿ ∑ ΔS x (i )= ∑ AT (i ). y GT (i)
¿
2
i=1 i=1

i=1,2,3,....NT

Avec NT :le nombre de triangle donnée .

AT (i ) :surface du triangle (i) i=1,2,3,....NT
y GT (i) :coordonnée y G de centre du gravité du triangle (i)
b) Pour un triangle(i) on a :
1
AT (i )= {(x 2 −x 1 )( y 3− y1 )−( y 2 − y 1 )(x 3 −x 1 ) }
2
1
y GT (i)= { y 1 + y 2 + y 3 }
3
c) Position des nœuds dans le repere
Gx 2 y 2 :
Le nombre totale de nœuds NN est determinée par :
NN =NT +2
On choisi la discritisation suivante :
Pour obtenir les positions des nœuds :

x1
=
( b ( y 2)
a ) +( 1+cos θ 2 ) cos θ 2−
xG
a 2 a
y1 y2
=
; a a
θ −θ
Δθ= 2 1
NT

θ=θ1 +(i−2) Δθ
R
=1+cos θ
a i=2,3,4....NN
xi xG
a
= ( )
R
a
cos θ−
a
yi
a
= ( )
R
a
sin θ

d) numérotation des nœuds de chaque triangle :

les nœuds de chaque triangle (i) ; sont numérotés dans le sens antihoraire par :
i , i +1 ,i +2 (i=1 , 2, 3, .. .. NT )
Finalement :
¿
S x i= NT AT ( i ) y ( i )
a
3
2
=∑
i =1 a
2 (
GT
a )( )
Avec : a 2
=
2 a
−
{[
AT (i) 1 x (i+1) x (1) y (i+2 ) y (1)
a a
−
a ][
−
a a ][ a ][
y (i+1 ) y (1) x(i+2 ) x (1)
− −
a ]}
Et a 3 a {
y GT (i ) 1 y (1) y (i+1 ) y (i+2 )
= +
a
+
a }
 {( ) }
S ¿x
2
3
T a
τ ( y 2 )=
a2
( b(ay ) )
Ix
2 2
4
a
D’où :

Programme :

Interpolation Cubic spline

L’interpolation est considérée comme le cœur de l’analyse [Link] l’existence des fonctions qui
sont tabulées et l’expression analytique est .absente sion veut chercher la valeur de la fonction pour
des points qui ne sont pas tabulés , il est necessaire d’interpolation ces points.

Aujourd’hui l’interpolation se fait par ordinateur.

 Interpolation linéaire
 Interpolation quadrature
 Interpolation cubique

Plusieurs auteurs ont developpé des procedures d’interpolations des points tabulé[Link] peut citer à
titre d’information :

 L’interpolation de lagrange.
 L’interpolation par iteration.
 Interpoltion inverse.
 Interpolation de Hermite.
 Interpolation par differences finis.
 Interpolation par cubique spline.

Dans le cas general ,il est possible qu’une séquence d’interpolation polynomiale { n } sur un
p (x)
nombre fixé et fini d’intervalles ne converge pas réguliérement à une fonction lissé,c’est-à-dire ,si on
veut chercher le comportement de ces polynomes entre les points interpolés,on observe dans
certains cas ,les ocillations des polynomes qui quittent volontairement la variation réelle de la fonction
.les polynomes d’interpolation d’interpolation qui possédent un degré grand donnent ces [Link]
chemin pour eviter ce probleme est de considérer une interpolation polynomiale avec un petit degré
sur des petits intervalles de l’intervalles considere . avec cela ,l’oscillation entre ces points n’est pas
significative,ou l’interpolation peut imiter le comportment de la [Link] ,le resultat de la
fonction interpolée à partir des petits intervalles avec des polynomes de degré faibles n’est pas
necessairement encore lisse.

Puisqu’on souhaite d’imiter le comportement de la foction ,il est recommandé que pour les polynomes
d’interpolation soient lisses

Ces fonctions avec un degré maximal sont nommees ‘spline’,parmi celle-ci on cite l’interpolation
‘cubic spline’et voici la description de sa formulation :

[Link] mathematique de l’interpolation’cubic spline’ :

L’interpolation cubique spline(dite par morceau) au point de vue mecanique ressemble à un procédé
qui consiste à joindre une serie des points par une tige elastique (spline).il existe plusieurs types de
splines ,lineaire ,quadratique ,cubique,…la plus utilisée c’est la spline cubique qui consiste à lisser
chaque segment de deux points successifs par une fonction cubique.

Considerons l’intervalle I=[ a , b ] divisée en ( n−1 ) sousintervalles pas necessairement de memes

longueurs

a=a 1 <a 2 <a 3 <a 4 <an−1 < an =b

La fonction spline S ( X ) de degré m est une fonction definie sur l’intervalle I

Par :

1) Considere avec le polynome de degre m sur chacun sous intervalle

I i=[ ai−1 , ai ] i=2,3,4,...,n
2) Possede des derivées continues jusqu’à l’ordre m−1

Les abscisses {
ai }
sont nommées les nœuds du ‘soline’.la fonction ‘spline’ S ( X ) est

l’interpolation des points [

ai , yi ]
si :
S ( xi )= yi i=1,2,3,...n
La forme la plus simple est du premier degré,nommé fonction d’interpolation linéaire qui n’est
pas lisse mais très sfficace si la distance entre les nœuds est [Link] effet,toutes les tables
des valeurs des fonctions qui utilisent l’interpolation linéaire conduire à une approximation
fondamentale des fonctions par ‘’ Spline Linéaire’’.Spline de degré deux peut etre
définie ,mais puisqu’il y a uniquement un degré de libertt=é dans sa défintion ,il y a manque
de symetrie dans leur determination avec relation des points de frontiere de l’[Link] plus
, les fonction resultantes ne sont pas suffisamment [Link] , en generale ,la
meilleure utilisation de’Spline’est celle de ‘Cubic-spline’ qui comporte deux paramètres à
choisir pour donner le comportement aux extrémites de l’intervalle.
Un type des conditions aux limites est tel que :
S'' ( a )=S '' ( b )=0
L’interpolation en ‘Cubic-Spline’qui satisfait ces conditions est nommée’Cubic-Spline
naturelle’
Une deuxieme condition qui generalement donne des meilleurs resultats est :
S' ( a ) =f ' ( a ), S ( b ) =f ( b )
' '

Cela ,évidement demande des donnees sur la derivee aux extrimites .

L’une parmi plusieurs representaions de ‘Spline ‘peut etre donne comme suit :
Soient
hi =a i−a i−1 i=2,3,4,...,n
Cependant, S ( X ) est une fonction cubique ,alors S' ( X ) a une forme quadratique
(2émedegree)et S'' ( X ) a une forme lineaire et [Link],on peut ecrire :
ai −x x−ai −1
S'' ( X ) =M i +Mi
hi hi I
sur i

Pour certaines constantes

M i ,en effet :
 S'' ( a )=M i i=2,3,4,...,n
Integrant l’equation () deux fois ecrivons la fonction lineaire arbitraire sous la forme
indiquee,on obtient :
3 3
( ai−x ) ( x−ai−1 )
S ( X ) =M i −1 +Mi +Ci ( ai −x ) +d i ( x−ai−1 )
6 hi 6 hi sur i
I
Cependant,on souhaite que la fonction ‘Spline’interpole aux nœuds ,on a :
S ( a i−1 ) = y i−1 S ( a i )= y i
,

Ces deux equations detrminent les costantes

Ci et
di ,la fonction S(X) devient :

( ) ( )
3 3
( ai−x ) ( x−ai−1 ) M i−1 hi ai −x M h x−a i−1
S ( X ) =M i−1 +Mi + y i−1− + yi− i i
6 hi 6 hi 6 hi 6 hi
Derivons l’expression () , on obtient :

' ( ai −x ) 2 ( x−ai ) 2 y i − y i−1 M i −M i−1

S ( X ) =−M i−1 + Mi + + hi
2h i 2 hi hi 6 I i=[ ai−1 , ai ]
sur
I i+1 =[ a i , ai+1 ]
Par analogie,sur l’intervalle on a :
2 2
( ai +1 −x ) ( x−ai ) y i +1− yi M i+1 −M i
S' ( X ) =−M i + M i +1 + + h i+1
2 hi +1 2 hi +1 h i+1 6 I i=[ ai−1 , ai ]
sur
On note les valeurs particulieres apres () et ()
hi hi y i − y i−1
S' ( a i )= M i−1 + Mi +
6 3 hi
h h y −y
S' ( a i )=− i+1 M i− i+1 M i+1 + i i−1
3 6 hi

Cependant , la fonction S' ( x ) doit etre continue,ces deux valeurs doivent etre egales , cela
donne ces equations :

hi h i +hi+1 hi+1 y i+1− y i y i − y i−1

M i−1 + M i+ M i+1= −
6 3 6 h i+1 hi i=2,3,4,...,n
On obtient () equations linéaire (à l’exception des deux nœuds des extrémités) dont les inconnus
M M M
sont 1 , 2 , 3 ,…..
M n , d’où deux conditions supplémentaires doivent etre ajoutées pour avoir une
solution unique.

Une fois les valeurs

Mi sontdeterminées ,l’interpolation Spline sera complètement déterminée par
l’equation ().

Pour simplifier l’ecriture de l’equation ()on pose :

y i− yi −1 h I +1 6 ( σ i +1−σ i )
σ i= λ i= di=
hI , h I +h I +1 , μi =1−λ I , h I + h I +1

Apres substitution de () dans () , le système d’equation possede la forme suivante :

μi M i−1 +2 M i +λ i M i+1=d i
Les deux conditions supplementaires possedent la forme :

2 M 1 +λ1 M 2 =d 1

μn M n−1 +2 M n=d n

Et on indique que plusieurs possibilité pour les choix des constantes

λ 1 ,d 1 , μn et d n .sous la forme

matricielle,on aura :

[ ]{ } { }
2 λ1 0 0 0 0 0 M1 d1
μ2 2 λ2 0 0 0 0 M2 d1
0 μ3 2 λ3 0 0 0 M3 d1
⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋅ ⋮ = ⋮
0 0 0 μn−2 2 λ n−2 0 M n−2 d n−2
0 0 0 0 μn−1 2 λn−1 M n−1 d n−1
0 0 0 0 0 μn 2 Mn dn

C’est un système d’equation à matrice tri-diagonales liées.

λ =d1 =μn =d n=0 , àlors

1ere Cas :si on choisit 1
M 1 =M n=0 .ceci donne l’interpolatin ‘Spline
naturelle’

2eme cas :selon la 2emecondition aux limites proposée au paravent :

'
S' ( x )= y '1 ( ) '
, S x = yn

On obtient :

λ 1=μ n=1

d 1=
h2(
6 y 2− y 1
h2
− y1
)
( )
6 y n − y n−1
d n= y n−
hn h2

a
3eme cas :l troisieme possibilité est de choisir i ,i=2,3,4,...,n−1 ,comme des nœuds du [Link]
points des extremités ne sont pas donnés comme des nœuds et d’imposer les conditions :

S' ( a ) = y '1

S' ( a n ) = y 'n

Dans ce cas , l’indice ‘i’ dans l’equation () varie entre 3 et n-2, et les deux conditions additionnelles
sont similaires en forme pour le choix de ().les valeurs de
λ 2 et d 2 peuvent etre rapidement evaluees

en supposant
x=a 1 dans () et posons S ( a ) = y 1 . De meme , en informations additionnelles.
De l’equation () :on sait que :

0< λi <1

0< μi <1 i=2,3,4,...,n−1

Si donc
|λ 1|<2 et |μ n|<2 ,la matrice dans () sera en diagonale [Link] ce cas ,il est possible
de montrer qu’une solution unique de ()existe pour 1 2
d ,d ,.....,d n arbitraire .cependant pour le 1ere cas
2emecas ont est assurée l’existance de l’interpolation Spline.

Similairement dans le 3emecas,la solution toujours existe :

A-2 Resolution d’un système d’equation lineaire àmatrice Tri-diagonale non sysmetrique par la
methode du double balayage de CHOLESKY :

dans le cas des matrices tri diagonales, la methode du pivot est parfois appelée methode du double
balayage de [Link] la donnerons parce que a mise en oeuvreest aisée,et que les matrices tri
diagonale interviennent dans de nombreux problemes,notamment pour interpoler une fonction
donnée sous forme d’un tableau.

Le système d’equations tri giagonale s’ecrit :

b 1 x 1 +c 1 x2 = y 1

a 2 x1 +b2 x 2 +c2 x 3 = y 2

Ka3 x2 +b3 x 3 +c 3 x 4 = y 3

0=M

a i x i−1 +bi x i +c i+1= y i

0=M

a n x n−1 +bn x n = y n

La premier équation de () permet d’exprimer l’inconnue

x2 :

−c 1 y1
x 1= x 2+
b1 b1

Portons cette valeur de x 1 dans la seconde equation de () on obtient :

a2
[ −c 1
b1
x 2+
y1
b1 ] +b 2 x2 +c 2 x 3= y 2

[ −c1 a1
b1 ]
+b2 x 2 =−c 2 x 3 + y 2−a2
y1
b1
 x x
On a ainsi exprimé 2 en fonction de 3 ,et ainsi de suite on exprimera chaque inconnue en fonction
de la [Link] qu’on ait obtenu :

x i−1 =A i−1 x i +Bi−1

Et portons cette valeur dans la ieme equation de (), il vient :

a i ( A i−1 x i + Bi−1 ) + bi x i +c i x i+1 = y i

−ci y i −ai Bi−1

x i= x i+1 +
a i A i−1 +b i ai Ai−1 + bi

On peut écrire donc :

En faisant i=1 dans les equations()et en comparant avec (),on voit que :

{¿

Les equations () et () nous permettent de calculer aisement de proche en proche, etant donnés
A0 =0 et B0 =0 , les couples ( A 1 , B1 ) , ( A 2 , B 2 ) , ( A 3 , B3 ) ,en utilisant un premier balayage de i=1 à n.

Nous allons maintenant effectuer un second balayage , de i=n à 1,qui va nous permettre de calculer
les inconnues .en effet, la derniere équation de () s’écrit :

a n x n−1 +bn x n = y n

x
Ce qui s’écrit en remplaçant n−1 par sa valeur tirée de () ou on fait i=n-1

a n ( A n−1 nn +B n−1 ) +bn x n = y n

D’où :

y n−an Bn−1
x n=
an A n−1 + bn

C’est -a-dire,d’apres() ou i=n

x n =B n

x
On obtenue ainsi la valeur cherchée n . On en deduit :

x n−1 = An−1 x n +Bn−1

Puis
x n−2 ,etc….jusqu’à
x 1 , de proche en proche