UFR S&T / UEVE Master 2/CS91 – Codes de calcul enseignant :Z-Q Feng
¾ Réaliser un code de calcul des structures en C++ par la méthode des éléments finis (CCMEF)
¾ Valider le code CCMEF avec le code industriel ANSYS
Projet 1
Etude d’un problème de transfert de chaleur par conduction et convection
avec une source de chaleur
Dans un sol, on place des câbles pour chauffer la surface du sol : par exemple dégel de la
piste d'atterrissage d'un aéroport. Les câbles produisent une chaleur de q W/cm. La température de
l'air est de -5 °C et le coefficient de convection est de 0.0034 W/cm2-°C. Les conductivités du sol
sont : kx = ky = 0.018 W/cm-°C.
2 cm
6 cm 4 cm
isolation
Les objectifs de ce projet sont :
- développer un élément triangulaire à 3 nœuds pour la conduction et la convection,
- développer un élément quadrilatéral à 4 nœuds pour la conduction et la convection,
- trouver q telle que la température de la surface du sol soit supérieure à 1 C°,
- étudier l’influence du maillage sur les résultats.
UFR S&T / UEVE Master 2/CS91 – Codes de calcul enseignant :Z-Q Feng
¾ Réaliser un code de calcul des structures en C++ par la méthode des éléments finis (CCMEF)
¾ Valider le code CCMEF avec le code industriel ANSYS
Projet 2
Contraintes dans une plaque percée d'un trou circulaire
Lorsqu'on perce un petit trou circulaire dans une plaque soumise à une contrainte de traction
uniforme σ, une concentration intense de contraintes se représente aux points nn. La théorie de
l'élasticité montre qu'en ces points la contrainte de traction est égale à 3σ.
a a
h
n
σ σ
n
h
Les données du problème sont :
σ = 100 MPa, a = 100 mm, h = 50 mm,
L'épaisseur de la plaque est égale à l'unité. Le rayon du trou est égal à 8 mm.
Matériau élastique: module de Young E = 200000 N/mm2
coefficient de Poisson ν = 0.3
Les objectifs de ce projet sont :
- développer un élément triangulaire à 3 nœuds dans CCMEF
- développer un élément quadrilatéral à 4 nœuds dans CCMEF
- calculer les contraintes dans la plaque avec ces deux types d’élément
- comparer les résultats numériques
Remarque : considérons la symétrie de la structure.
UFR S&T / UEVE Master 2/CS91 – Codes de calcul enseignant :Z-Q Feng
¾ Réaliser un code de calcul des structures en C++ par la méthode des éléments finis (CCMEF)
¾ Valider le code CCMEF avec le code industriel ANSYS
Projet 3
Analyse statique des structures de poutres
I. Poutre en flexion simple
On considère un domaine plan rectangulaire (poutre) de longueur L, de hauteur h et de
largueur b << h, L. Cette poutre est encastrée à une extrémité et soumise à une charge P à l’autre
extrémité (voir figure).
b
h
P
L
Les données du problème sont :
L=200 mm. h=20 mm. b=2mm. P = 100 N
Matériau élastique: module de Young E = 200000 Mpa ; coefficient de Poisson ν = 0.3
II : Poutre courbe
Une poutre courbe de faible épaisseur est articulée sur deux extrémités et supporte une charge
verticale P. Sa section est uniforme dont la ligne moyenne est un demi-cercle.
L P = 200 N
E = 3x105 MPa
L = 100 mm
Section rectangulaire
b a= 1 mm, b=2 mm
a
L L
On demande de :
- développer un élément de poutre 2D à 2 nœuds dans CCMEF
- déterminer les ddl et les contraintes et comparer les résultats numériques et théoriques
- développer un élément triangulaire à 3 nœuds dans CCMEF pour modéliser les poutres en
solide 2D.
PL3
; (II) δ = PL 3π - 1 - 1
3
Solution RDM : (I) δ =
3EI EI 8 2π
UFR S&T / UEVE Master 2/CS91 – Codes de calcul enseignant :Z-Q Feng
¾ Réaliser un code de calcul des structures en C++ par la méthode des éléments finis (CCMEF)
¾ Valider le code CCMEF avec le code industriel ANSYS
Projet 4
Structure de treillis
On demande de :
1. développer un élément de barre 2D à 2 nœuds dans CCMEF
2. développer un élément de barre 3D à 2 nœuds dans CCMEF
3. Calculer les déplacements des nœuds, les forces internes dans les barres et les réactions
d'appuis.
4. Vérifier les résultats de votre code par ceux du code ANSYS
I) Treillis 2D à 15 barres
400 50 kN
80 kN 300 cm
50 kN
200 400
A=10 cm2 , E = 2x107 N/cm2
100
P P=1000 N
III) Treillis 3D
200
Matériau : Aluminum
On demande de
dimensionner les barres
pour que la contrainte
axiale ne depasse pas
150 MPa
200
200
200 mm
UFR S&T / UEVE Master 2/CS91 – Codes de calcul enseignant :Z-Q Feng
¾ Réaliser un code de calcul des structures en C++ par la méthode des éléments finis (CCMEF)
¾ Valider le code CCMEF avec le code industriel ANSYS
Projet 5
Problème de torsion
L’équation régissant un problème de torsion est
d 2φ d 2φ
+ + 2gθ = 0
dx 2 dy 2
avec φ = 0 sur les contours
Considérons un arbre de section carrée, soumis à un couple T.
T 1 cm
1 cm
Données : module de cisaillement g = 8 106 N/cm2 ; angle de torsion θ = 0,01 °/cm.
Les objectifs de ce projet sont :
1. développer un élément triangulaire à 3 nœuds dans CCMEF
2. développer un élément quadrilatéral à 4 nœuds dans CCMEF
3. calculer la valeur de la fonction de contrainte et les contraintes de cisaillement.
4. Calculer le couple T.
5. Modéliser un huitième de la section avec 4 éléments T3.
6. Modéliser un quatrième de la section avec 100 éléments Q4.
Remarque : Par analogie, utiliser l’élément de type thermique (plane55 dans ANSYS).
UFR S&T / UEVE Master 2/CS91 – Codes de calcul enseignant :Z-Q Feng
¾ Réaliser un code de calcul des structures en C++ par la méthode des éléments finis (CCMEF)
¾ Valider le code CCMEF avec le code industriel ANSYS
Projets 6 et 7
Calcul du champ de températures et des contraintes thermomécaniques
Les objectifs de ce projet sont :
¾ développer un élément triangulaire à 3 nœuds et un élément quadrilatéral à 4 nœuds pour la
conduction et un élément linéaire à 2 nœuds pour la convection dans CCMEF pour calculer le
champ de température dans la structure (Projet 6)
¾ développer un élément triangulaire à 3 nœuds et un élément quadrilatéral à 4 nœuds dans
CCMEF pour calculer les contraintes thermique et/ou mécanique dans la structure (Projet 7)
UFR S&T / UEVE Master 2/CS91 – Codes de calcul enseignant :Z-Q Feng
¾ Réaliser un code de calcul des structures en C++ par la méthode des éléments finis (CCMEF)
¾ Valider le code CCMEF avec le code industriel ANSYS
Projet 8 : Contraintes dans une structure 3D
Considérons une structure en T soumise à deux forces concentrées P et F. Les objectifs de ce projet
sont :
1. Développer un élément tédraedrique à 4 nœuds dans CCMEF,
2. Comparer deux solutions I et II en effectuant des calculs par éléments finis,
3. Trouver la contrainte de von Mises maximale.
Solution I
b a F
c
d
d
Solution II
r
b a F
c
d
d
Données du problème :
• Géométrie : a = 120 mm ; b = 300 mm ; d = 250 mm ; c = 240 mm ; r = 40 mm
• Matériau : Module de Young E = 200000 Mpa ; Coefficient de Poisson ν = 0.3
• Sollicitations : P = ± 100 N, F = ± 150 N
• Conditions aux limites : la structure est encastrée sur sa base.
UFR S&T / UEVE Master 2/CS91 – Codes de calcul enseignant :Z-Q Feng
¾ Réaliser un code de calcul des structures en C++ par la méthode des éléments finis (CCMEF)
¾ Valider le code CCMEF avec le code industriel ANSYS
Projet 9
Contraintes dans une structure 3D
Considérons une structure en T soumise à une force P. Les objectifs de ce projet sont :
1. Développer un élément hexaédrique à 8 nœuds dans CCMEF,
2. Comparer deux solutions I et II en effectuant des calculs par éléments finis,
3. Trouver la contrainte de von Mises maximale.
Solution I
b a F
c
d
d
Solution II
r
b a F
c
d
d
Données du problème :
• Géométrie : a = 120 mm ; b = 300 mm ; d = 250 mm ; c = 240 mm ; r = 40 mm
• Matériau : Module de Young E = 200000 Mpa ; Coefficient de Poisson ν = 0.3
• Sollicitations : P = ± 100 N, F = ± 150 N
• Conditions aux limites : la structure est encastrée sur sa base.
UFR S&T / UEVE Master 2/CS91 – Codes de calcul enseignant :Z-Q Feng
¾ Réaliser un code de calcul des structures en C++ par la méthode des éléments finis (CCMEF)
¾ Valider le code CCMEF avec le code industriel ANSYS
Projet 10
Calcul des contraintes dans une charpente
Une charpente d’une section rectangulaire est soumise à une charge repartie p (voir figure).
p = 0.1 MPa
p E = 3x104 MPa
L=4m
c=1m
α
CAS 1 : α = 120 °
a= 30 mm, b=300 mm
b CAS 2 : α = 100 °
a= 30 mm, b=300 mm
a CAS 3 : α = 120 °
c c a= 300 mm, b=30 mm
On demande de :
- développer un élément de poutre 2D (modèle M+N) à 2 nœuds dans CCMEF
- développer un élément triangulaire à 3 nœuds dans CCMEF pour modéliser la charpente
en solide 2D.
- déterminer les DDLs et les contraintes dans la charpente.
- comparer les résultats numériques obtenus par les deux types d’élément au cas 1.
- comparer les résultats numériques obtenus avec l’élément de poutre aux différents cas.
UFR S&T / UEVE Master 2/CS91 – Codes de calcul enseignant :Z-Q Feng
¾ Réaliser un code de calcul des structures en C++ par la méthode des éléments finis (CCMEF)
¾ Valider le code CCMEF avec le code industriel ANSYS
Projet 11
Contraintes dans un barrage
Considérons un barrage dont la section est constante (ABCD avec AB⊥AD, AB//CD et
AB=AD=3CD). Une pression uniforme p exerce sur le côté AD et le barrage est soumis à son poids
propre. La base du barrage est encastrée. Un tunnel de passage (abcd avec ab=cd) est placé dans le
barrage.
D C
d c
A a b B
Les données du problème sont :
p = 50 MPa, AB = 30 m, Aa = 10 m, ab=3 m
Matériau élastique: module de Young E = 20000 N/mm2
coefficient de Poisson ν = 0.3
masse volumique ρ = 2500 kg/m3
Les objectifs de ce projet sont :
- développer un élément triangulaire à 3 nœuds dans CCMEF
- développer un élément quadrilatéral à 4 nœuds dans CCMEF
- calculer les contraintes dans le barrage avec ces deux types d’élément
- comparer les résultats numériques
- refaire le calcul avec 1) Aa=5 m ; 2) Aa=15 m et commenter les résultats
UFR S&T / UEVE Master 2/CS91 – Codes de calcul enseignant :Z-Q Feng
¾ Réaliser un code de calcul des structures en C++ par la méthode des éléments finis (CCMEF)
¾ Valider le code CCMEF avec le code industriel ANSYS
Projet 12
Contraintes dans un anneau
Considérons un anneau soumis à un chargement de traction F.
Les données du problème sont :
Rayon intérieur: 0.05 m; Rayon extérieur: 0.12 m; Epaisseur: 0.005 m
Chargement de traction : F = 4000 N
Matériau élastique: module de Young E = 7 1010 Pa
coefficient de Poisson ν = 0.3
Les objectifs de ce projet sont :
- développer un élément triangulaire à 3 nœuds dans CCMEF
- développer un élément quadrilatéral à 4 nœuds dans CCMEF
- calculer les contraintes dans l’anneau avec ces deux types d’élément
- comparer les résultats numériques
Remarque : considérons la symétrie de la structure.
