République Algérienne Démocratique et Populaire
Ministère de l’Enseignement Supérieur et de la
Recherche Scientifique
ـــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــ
ــــــــــــ
Université Amar Telidji- Laghouat
Faculté de Technologie
Département de mécanique
Spécialité : Maintenance Industrielle
ـــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــ
ــــــــــــــــــــــــــــــــــــ
ــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــــ
ــــــــــــــــ
TP 1
Poutre Encastrée à une Extrémité
Réalisé par : HAMMANI MOHAMMED ELHABIB & GAFSI ABDELHAKIME
Instructeur: DR, Abdelkader Djilal Hammou
But du Tp :
Cette étude a pour objectif d’explorer les vibrations d'une poutre en encastrée à une
extrémité afin d'identifier ses fréquences propres. En combinant des calculs analytiques
et des simulations numériques réalisées avec le logiciel ABAQUS, nous analysons le
comportement dynamique de la structure. Cela permet de comparer les résultats
théoriques avec ceux issus des
simulations.
I. Partie Théorique
Détermination des premières fréquences propres :
U(x, t) = (A1 cos(βnx) + A2 sin(βnx) + A3 cosh(βnx)
+ A4 sinh(βnx)) (𝐵1𝑠𝑖𝑛𝝎𝑡 + 𝐵2𝑐𝑜𝑠𝝎𝑡)
La partie spatiale : 𝐴1, 𝐴2, 𝐴3, 𝐴4 est trouvée en utilisant les
La partie temporelle∶ (𝐵1𝑠𝑖𝑛𝝎𝑡 + 𝐵2𝑐𝑜𝑠𝝎𝑡) est trouvée en utilisant
conditions aux limites.
les conditions initiales.
Nous travaillons uniquement sur la recherche des formes de mode, nous
pouvons donc ignorer La partie temporelle.
Conditions aux limites :
x=0 (condition encastré):
𝑈(0, 𝑡) = 0
𝐴1 𝑐𝑜𝑠(0) + 𝐴2 𝑠𝑖𝑛(0) + 𝐴3 𝑐𝑜𝑠ℎ(0) + 𝐴4 𝑠𝑖𝑛ℎ(0) = 0
Depuis 𝑐𝑜𝑠(0) = 1, 𝑠𝑖𝑛(0) = 0, 𝑐𝑜𝑠ℎ(0) = 1, 𝑠𝑖𝑛ℎ(0) = 0,
Cela se simplifie en : A1+A3 =0
𝑈′(0, 𝑡) = 0∶
𝑈′(0) = −A1βnsin(0) + A2βncos(0) + A3βnsinh(0)
+ A4βncosh(0) = 0
Depuis 𝑐𝑜𝑠(0) = 1, 𝑠𝑖𝑛(0) = 0, 𝑐𝑜𝑠ℎ(0) = 1, 𝑠𝑖𝑛ℎ(0) = 0,
Cela se simplifie en : A2+A4=0
𝑈(𝐿, 𝑡) = 0∶
x=L (condition libre):
−𝐴1𝛽𝑛2𝑐𝑜𝑠(𝛽𝑛𝐿) − 𝐴2𝛽𝑛2𝑠𝑖𝑛(𝛽𝑛𝐿) +
𝐴3𝛽𝑛2𝑐𝑜𝑠ℎ(𝛽𝑛𝐿)
+ 𝐴4𝛽𝑛2𝑠𝑖𝑛ℎ(𝛽𝑛𝐿) = 0
𝑈′(, 𝑡) = 0∶
𝐴1𝑠𝑖𝑛(𝛽𝑛𝐿) − 𝐴2𝑐𝑜𝑠(𝛽𝑛𝐿) + 𝐴3𝑠𝑖𝑛ℎ(𝛽𝑛𝐿)
+ 𝐴4𝑐𝑜𝑠ℎ(𝛽𝑛𝐿) = 0
À partir des conditions aux limites à x=0 ,nous avons :
A1+A3 =0 & A2+A4=0
Cela implique :
A3=-A1 & A4=-A2
On Remplacez A3=-A1 & A4=-A2 dans les conditions aux limites à
(x=L) :
Pour 𝑈(𝐿, 𝑡) = 0∶
−𝐴1𝛽𝑛2𝑐𝑜𝑠(𝛽𝑛𝐿) − 𝐴2𝛽𝑛2𝑠𝑖𝑛(𝛽𝑛𝐿) +
𝐴3𝛽𝑛2𝑐𝑜𝑠ℎ(𝛽𝑛𝐿)
+ 𝐴4𝛽𝑛2𝑠𝑖𝑛ℎ(𝛽𝑛𝐿) = 0
𝐴1(𝑐𝑜𝑠(𝛽𝑛𝐿) + 𝑐𝑜𝑠ℎ(𝛽𝑛𝐿)) + 𝐴2(𝑠𝑖𝑛(𝛽𝑛𝐿) + 𝑠𝑖𝑛ℎ(𝛽𝑛𝐿)) = 0
Donc:
(1)
Pour 𝑈′(𝐿, 𝑡) = 0∶
𝐴1𝑠𝑖𝑛(𝛽𝑛𝐿) − 𝐴2𝑐𝑜𝑠(𝛽𝑛𝐿) + 𝐴3𝑠𝑖𝑛ℎ(𝛽𝑛𝐿)
+ 𝐴4𝑐𝑜𝑠ℎ(𝛽𝑛𝐿) = 0
Donc:
𝐴1(𝑠𝑖𝑛(𝛽𝑛𝐿) − 𝑠𝑖𝑛ℎ(𝛽𝑛𝐿)) − 𝐴2(𝑐𝑜𝑠(𝛽𝑛𝐿) + 𝑐𝑜𝑠ℎ(𝛽𝑛𝐿)) = 0
(2)
Les équations (1) et (2) peuvent être combinées pour donner
𝑐𝑜𝑠(𝛽𝑛𝐿)−𝑐𝑜𝑠ℎ(𝛽𝑛𝐿)
𝐴2 = −A1
(𝑠𝑖𝑛(𝛽𝑛𝐿)+𝑠𝑖𝑛ℎ(𝛽𝑛𝐿)
(3)
En insérant l'équation (3) dans (1) ou (2), on obtient l'équation
de fréquence pour une poutre en porte-à-faux,
1 + 𝑐𝑜𝑠(𝛽𝐿)𝑐𝑜𝑠ℎ(𝛽𝐿) = 0
Cette équation est transcendante et doit être résolue
numériquement. Les solutions sont :
βL = 1,875; 4,694; 7,855; 10,996; …
Calculs détaillés :
Section rectangulaire : b=0,10 m, h=0,10 m ,L=1 m
A = b⋅h
Étape 1 : Calcul de la surface de la section transversale (A)
A = 0,10 m⋅0,10 m = 0,01 m²
Substituons les valeurs :
Étape 2 : Calcul du moment d'inertie (Iy)
Iy = b⋅h³/12
Substituons les valeurs :
Iy = 0,10⋅(0,10)³/12 = 8,33×10⁻⁶ m⁴
Étape 3 : Calcul de la densité linéaire (μ)
μ = ρ⋅A
Substituons les valeurs :
μ = 7850 kg/m³⋅0,01 m² = 78,5 kg/m
Étape 4 : Calcul de la rigidité en flexion (EI)
EI = E⋅Iy
Substituons les valeurs :
EI = 210×10⁹ Pa⋅8,33×10⁻⁶ m⁴ = 1,7493×10⁶ Nm²
Étape 5 : Calcul des fréquences naturelles (fn)
La formule des fréquences naturelles pour une poutre en porte-à-faux est
𝛚𝑖 = 2𝜋𝑓𝑖 donc: 𝑓𝑖 = 𝛚𝑖/2𝜋
:
où βL est une constante dépendant du mode.
Mode 1:
(βL)² = (1,875)² = 3,5156
ω₁ = ((β₁L)² / L²) * (EI / μ)^(1/2)
ω₁ = ((1.875)² / (1 m)²) * (1.75 x 10⁶ Nm² / 78.5 kg/m)^(1/2)
ω₁ = (3.5156 / 1) * (22305.088)^(1/2) ω₁ = 3.5156 * 149.35
ω₁ ≈ 524.1 rad/s
1er Fréquence naturelle :
f₁ = w₁/2π = 524.1 /2π = 83,41 Hz
Mode 2 :
ω2 = ((βL)² / L²) * (EI / μ)^(1/2)
ω2 = (22.024 / (1 m)²) * (1.75 x 10⁶ Nm² / 78.5 kg/m)^(1/2)
ω2 = 22.024 * (22305.088)^(1/2)
ω2 = 22.024 * 149.35 ω₁ ≈ 3289.4 rad/s
2éme Fréquence naturelle :
f₂ = w₂/2π = 3289.4 /2π = 523,52 Hz
II. Partie Pratique Simulation (ABAQUS) :
Longueur (L) 1m
Largeur (b) 0.10 m
Hauteur (h) 0.10 m
Module de Young (E) 210 Gpa
Densité de l'acier (𝜌) 7850 kg/m³
Coefficient de Poisson 0,3
Mode : 01 Frequence = 83.283 Hz
Mode : 02 Frequence= 499.45 HZ
Mode Frequence Theorique Frequence Pratique Difference relative
Mode 1 0.21%.
83,59 83.41 Hz
Mode 2 523,52 Hz 499.45 Hz 4.70%
L'analyse : met en évidence une excellente corrélation pour le Mode 1, avec un
écart relatif très faible de 0.37%..Ces résultats confirment la validité du modèle
théorique pour cette fréquence.
pour le Mode 2, un écart relatif plus important de 4.62% est observé. Cet écart
pourrait être attribué à des effets non linéaires, tels que des déformations non
proportionnelles à la sollicitation, des interactions non linéaires entre les différents
degrés de liberté, ou des approximations du modèle théorique qui deviennent
moins pertinentes pour les modes de vibration supérieurs
Conclusion :
Les résultats de cette étude ouvrent des perspectives intéressantes pour
l'amélioration des modèles numériques utilisés en analyse vibratoire. En identifiant
les limites d'ABAQUS pour les modes supérieurs, cette étude souligne la nécessité
de développer des modèles plus sophistiqués capables de prendre en compte les
effets non linéaires et les interactions complexes entre les différents composants
d'une structure. De plus, l'intégration de données expérimentales dans les modèles
numériques pourrait permettre d'améliorer encore la précision des prédictions.