Cours MEF2
Cours MEF2
Cours MEF2
Paul
Sabatier
Toulouse III
M.E.F.
Michel SUDRE
Sept 2013
M.E.F.
M SUDRE
Rappels
1 Mthode de Ritz
L
(1)
(2)
v(x) = a. sin (
x )
L
W =
Fi . u i
Si on admet que les forces Fi drivent dun potentiel U, alors le rsultat prcdent
peut sexprimer comme une condition dextmum:
Fi = - U
ui
Fi . ui = - U do (W+U) = 0
U est loppos de , travail des forces Fi calcul dans le dplacement virtuel (en supposant les forces Fi constantes).
W- est lnergie potentielle totale (EPT). On montre que, pour un corps en quilibre
stable, cet extrmum de lnergie potentielle totale est un minimum absolu.
On nonce le principe du minimum de lnergie potentielle totale ainsi:
Parmi tous les champs cinmatiquement admissibles,
celui qui rend minimale lnergie potentielle totale
correspond la solution.
M.E.F.
M SUDRE
La meilleure approximation est celle qui rend extrmale lnergie potentielle totale.
Elle est donc obtenue par les N relations:
(EPT)
=0
i=1..N
ai
La mthode peut tre illustre par le problme de flexion suivant:
A
2L
2
2
Mz2 dx
1.
1 . ( EIz d v ) dx
=
2
dx
2 0 EIz
2 EIz 0
-F.v(L) -F.v(2L)
M.E.F.
M SUDRE
(W- ) = 0
a1= - 11FL
8EIz
(W- ) = 0
a2=
a1
a2
F
4EIz
3
a1.(2L)2 + a2.(2L)3 = -
7FL3
2EIz
2.1 Principe
La Mthode des Elments Finis consiste dcouper la structure en lments de forme simple et choisir une approximation du dplacement sur chaque subdivision.
Cest une mthode de Ritz par morceaux qui sadapte aux gomtries les plus complexes.
lment
noeud
Les subdivisions sont les lments et les connexions entre lments sont les noeuds.
Il est utile que les inconnues soient des paramtres physiques. Cest pourquoi on choisit les composantes de dplacements des noeuds (dplacements nodaux).
M.E.F.
M SUDRE
Par exemple: (
I, I ,
I)
(
I
I, I ,
I) sont
.
.
.
ui
.
.
.
Lnergie potentielle totale (W- ) est calcule en fonction des n variables ui.
Il existe 3 grandes familles dlments. Les lments uni-dimensionnels, bi-dimensionnels, tri-dimensionnels.
Lnergie lastique W est calcule en fonction des dplacements nodaux:
En 1D, grce aux relations de la thorie des poutres.
Un modle uni-dimensionnel sappuie donc sur les hypothses restrictives de la rdm.
En 2D, en utilisant, selon les cas, les thories des membranes, des plaques
et des coques.
Un modle bi-dimensionnel sappuie sur les hypothses restrictives de ces diffrentes
thories.
En 3D, en exploitant les relations de llasticit sans aucune hypothse restrictive.
Le travail des efforts extrieurs sexprime simplement en fonction des dplacements
des noeuds auxquels sont appliqus ces efforts.
La solution recherche correspond un minimum de lnergie potentielle totale.
Elle est obtenue par les conditions:
(EPT)
=0
ui
i=1..n
[K].{U} = {F}
dplacements nodaux
efforts appliqus
M.E.F.
M SUDRE
3 Exemple dapplication
exemple
On considre un assemblage de 4 ressorts de raideur k lis au bti
et soumis 4 efforts dintensit F.
F
k
F
k
u1
u2
u3
u4
u5
uJ
uJ uI
L
Donc lnergie lastique de llment (i) vaut:
= { 1, 1 }
2
1
1
W i = --- k = ---
2
2
uI uJ
uI
uJ
k k
uI
uJ
3.2 assemblage
Lnergie lastique totale sobtient par addition des contributions de chaque lments.
1
W = W 1 + W 2 + W 3 + W 4 = ---
2
u1 u2 u3 u4 u5
k k 0 0
k 2k k 0
0 k 2k k
0 0 k 2k
0 0 0 k
u1
u2
u3
u4
u5
M.E.F.
M SUDRE
EPT = W
u1 u2 u3 u4 u5
X1
F
F
F
F
3.4 minimisation
La minimisation de lnergie potentielle totale conduit lcriture du systme:
k k 0 0
k 2k k 0
0 k 2k k
0 0 k 2k
0 0 0 k
u1
u2
u3
u4
u5
X1
F
= F
Ce systme ne peut tre rsolu car il nexiste pas de solution unique tant quun dplacement densemble de la structure reste possible. La matrice [K] obtenue par assemblage est singulire.
0 0 0
2k k 0
k 2k k
0 k 2k
0 0 k
u1
u2
u3
u4
u5
= F
3.6 rsolution
u1
u2
u3
u4
u5
F
- = k 7
10
M.E.F.
M SUDRE
F
-X1 = k k 0 0 0
0
4
7
9
10
= 4F
On trouve:
N1 = k 1 = 4F
= k ( u J u I ) = { k, k }
N 2 = k 2 = 3F
uI
uJ
N3 = k 3 = 2F
N4 = k 4 = F
M.E.F.
M SUDRE
Nastran
1 Prsentation
NASTRAN (NAsa STRuctural ANalysis) est un logiciel conu par la NASA dans les annes 60. La version utilise ici est la version MSC/NASTRAN dveloppe et maintenue
par MacNeal-Schwendler Corporation.
Des ensembles pr-programms en langage DMAP (Direct Matrix Abstraction Programming) sont disponibles pour effectuer les tches les plus courantes:
sol 101
Analyse Statique Linaire
sol 103
Frquences naturelles et modes de vibration
sol 105
Flambement
sol 106
Analyse Statique Non Linaire
sol 200
Optimisation
.....
Le langage DMAP offre lusager la possibilit dcrire ses propres programmes pour
des besoins particuliers.
Plusieurs logiciels compatibles avec MSC/NASTRAN peuvent tre utiliss pour gnrer
les modles. Les plus connus sont PATRAN, MSC/NASTRAN for Windows, I-DEAS.
Bulk Data
M.E.F.
M SUDRE
3 Paramtres standardiss
SID Set identification number
Grid point
Degree of freedom
EID
PID
1
2
3
4
5
6
Tx
Ty
Tz
Rx
Ry
Rz
4 Elments de NASTRAN
C....
P....
Dimension 1
CELAS
PELAS
CROD
PROD
CBAR
PBAR
Dimension 2
CSHEAR
PSHEAR
Dimension 3
CTETRA CHEXA
PSOLID
Element
nu
CROD
Traction
Torsion
CBAR
Traction
Flexion
Torsion
Tranchant
Membrane et Flexion
Cisaillement
Transverse
Cisaillement
CSHEAR
CTETRA CHEXA
Toutes sollicitations
10
M.E.F.
M SUDRE
5 Commandes de NASTRAN
5.1 Gnralits
En format standard, une carte NASTRAN est constitue de 10 champs de 8 caractres.
Le premier champ contient le nom de la commande (GRID, CBAR ..). Les champs de
2 9 contiennent les paramtres associs cette commande.
Le 10me champ ne contient pas de donne. Il est rserv la marque de continuation.
En effet, certaines commandes ncessitent plus de 8 champs de donnes. (voir cidessous la commande CBAR )
En format libre, les champs sont spars par des virgules. Ils ne doivent pas dpasser
8 caractres.
Un entier ne peut pas contenir de point dcimal. Un rel doit contenir un point dcimal.
Un rel peut tre crit de diffrentes manires:
7.0
.7E1
0.7+1
70.-1
GRID
dfinit la position dun noeud, les directions de ses dplacements et ses SPC.
ID
CP
X1
X2
X3
cylindrique
x
r
sphrique
cartsien
Repre
ID:
CP:
X1,X2,X3:
CD:
PS:
11
CD
PS
M.E.F.
M SUDRE
EID
CROD
PID
G1
G2
EID:
PID:
G1,G2:
PID
PROD
MID
NSM
PID:
MID:
A:
J:
C:
NSM:
CBAR
EID
PID
GA
GB
X1
X2
X3
+A1
PA
PB
W1A
W2A
W3A
W1B
W2B
EID:
PID:
GA,GB:
X1,X2,X3:
PA,PB:
W1A,W2A,W3A,W1B,W2B,W3B:
12
+A1
W3B
M.E.F.
M SUDRE
autre configuration:
CBAR
EID
PID
GA
GB
G0
PA
PB
W1A
W2A
W3A
W1B
W2B
W3B
G0:
PBAR
PID
MID
I1
I2
NSM
C1
C2
D1
D2
E1
E2
F1
F2
K1
K2
I12
G3
G4
G5
G6
PID:
MID:
A:
I1,I2:
J:
NSM:
C1,C2,D1,D2,E1,E2,F1,F2:
K1,K2:
I12:
SPC1
SID
G1
G2
+A1
G7
G8
G9
....
SID:
C:
Gi:
13
+A1
M.E.F.
M SUDRE
autre configuration:
SID
SPC1
G1
THRU
G2
Pour que la carte SPC1 soit prise en compte, il faut placer la commande SPC=SID
dans la partie Case Control Deck du fichier de donnes.
MPC
MPC
SID
G1
C1
A1
+A1
G3
C3
A3
....
G2
C2
Ai.ui = 0 .
A2
+A1
SID:
Gi:
Ci:
Ai:
Pour que la carte MPC soit prise en compte, il faut placer la commande MPC=SID
dans la partie Case Control Deck du fichier de donnes.
SID
FORCE
CID
N1
N2
N3
SID:
G:
CID:
F:
N1,N2,N3:
Pour que la carte FORCE soit prise en compte, il faut placer la commande LOAD=SID
dans la partie Case Control Deck du fichier de donnes.
14
M.E.F.
M SUDRE
SID
MOMENT
CID
N1
N2
N3
SID:
G:
CID:
M:
N1,N2,N3:
Pour que la carte MOMENT soit prise en compte, il faut placer la commande
LOAD=SID dans la partie Case Control Deck du fichier de donnes.
LOAD
SID
S1
+A1
S4
L4
....
L1
S2
L2
S3
L3
+A1
SID:
S:
Si:
Li:
Pour que la carte LOAD soit prise en compte, il faut placer la commande LOAD=SID
dans la partie Case Control Deck du fichier de donnes.
MID
MAT1
NU
RHO
MID:
E:
G:
15
TREF
GE
M.E.F.
M SUDRE
NU:
RHO:
A:
TREF:
GE:
RBE2
EID
GN
CM
GM1
GM2
EID:
GN:
CM:
GMi:
16
GM3
....
M.E.F.
M SUDRE
Exercice 1:
Une poutre de longueur 0.5m et de section A= 8.E-5 m2 est soumise un effort de
tension F= 4000 N.
1 question:
Ecrire le fichier *.bdf.
y
.5m
4000 N
x
2 question:
Effectuer le calcul avec Nastran et lire lallongement dans le fichier *.f06.
3 question:
Vrifier lallongement par un calcul analytique.
17
M.E.F.
M SUDRE
Exercice 2:
Voici lexemple suivant dont le fichier *.bdf est fourni. Llment 1-2 a pour section A1
et llment 2-3 a pour section A2. Laxe local 2 coincide avec la direction y.
1 question: Complter les cartes "Load".
y
1000 N
1m
1m
500 Nm
x
1
section A2:
2
F
1
E
0.05
0.05
1
E
0.1
0.05
$ Exemple 2
SOL 101
CEND
$
DISPLACEMENT(SORT1,REAL)=ALL
STRESS(SORT1,REAL,VONMISES,BILIN)=ALL
SPC = 20
LOAD =
$
BEGIN BULK
PBAR,1,1,5.000E-3,1.042E-6,4.167E-6,2.861E-6,0.0
,2.500E-2,5.000E-2,-2.50E-2,5.000E-2,-2.50E-2,-5.00E-2,2.500E-2,-5.00E-2
,8.338E-1,8.333E-1,0.0
CBAR,1,1,1,2,0.,1.,0.
PBAR,2,1,2.500E-3,5.208E-7,5.208E-7,8.791E-7,0.0
,2.500E-2,2.500E-2,-2.50E-2,2.500E-2,-2.50E-2,-2.50E-2,2.500E-2,-2.50E-2
,8.334E-1,8.334E-1,0.0
CBAR,2,2,2,3,0.,1.,0.
MAT1,1,2.1+11,,.3
GRID,1,,0.,0.,0.
GRID,2,,1.,0.,0.
GRID,3,,2.,0., 0.
LOAD,
SPC1,20,123456,1
FORCE,
MOMENT,
ENDDATA
18
M.E.F.
M SUDRE
19
M.E.F.
M SUDRE
PID
MID1
+A1
Z1
Z2
MID4
MID2
RI
RS
MID3
NSM
+A1
MID1
Fxy
MID2
y
x
Fy
My
y
x
Mxy
Fxy
Mx
Fx
Mxy
Qy
Qx
I=
Ta
Ta
Tb
Tb
T = Ta + Tb
MID3
I=
T3
12
RS=5/6
RI=1
ha
hb
Ta 3 Tb3
2
2
+ Ta ha + Tbhb
+
12 12
20
RS=TS/T
M.E.F.
M SUDRE
CQUAD4
EID
+A1
PID
G1
G2
G3
G4
T1
T2
T3
T4
THETA
(MCID) ZOFFS
+A1
EID:
PID:
Gi:
THETA:
MCID:
ZOFFS:
Ti:
Exercice 3:
Soit une plaque encastre de dimensions 100 mm x 100 mm, modlise par un lment
QUAD4, charge par un effort de membrane de 10 000 N et par un effort transverse
de 100 N.
z
4
100 N
5.
10 000 N
100 N
100.
100.
10 000 N
matriau: E=70 000 MPa / nu=.3
1 question:
Ecrire le fichier *.bdf en demandant le calcul des dplacements, contraintes et forces.
21
M.E.F.
M SUDRE
2 question:
Retrouver par le calcul une partie des rsultats fournis par le fichier *.f06:
Dplacements:
POINT
ID.
TYPE
T1
D I S P L A C E M E N T
T2
T3
V E C T O R
R1
R2
R3
Contraintes:
ELEMENT
ID.
1
S T R E S S E S
I N
Q U A D R I L A T E R A L
E L E M E N T S
FIBER
STRESSES IN ELEMENT COORD SYSTEM
PRINCIPAL STRESSES (ZERO SHEAR)
DISTANCE
NORMAL-X
NORMAL-Y
SHEAR-XY
ANGLE
MAJOR
MINOR
-2.5E-3
2.5E-3
Forces:
F O R C E S
ELEMENT
ID.
I N
FX
Q U A D R I L A T E R A L
- MEMBRANE FORCES FY
FXY
22
E L E M E N T S
- BENDING MOMENTS MX
MY
MXY
M.E.F.
M SUDRE
1 Objectif.
Soit une plaque troue de dimensions 100 mm x 200 mm, dpaisseur 5 mm.
Elle est encastre sur les bords infrieur et suprieur et soumise une pression de matage p() lintrieur du trou.
200
zone2
zone1
p=200.sin() N/mm2
R2
R1
long
fonction: sinr()
100
larg
Lobjectif est:
-de paramtrer la gomtrie et le maillage grce au langage pcl,
-de crer une fonction pour appliquer la pression de matage,
-de tester les possibilits de visualisation des rsultats.
Le problme tant symtrique, il suffit de considrer la moiti de la plaque. La zone1
sera maille en "isomesh", la zone2 en "paver". Des lments rod de section 0.01
mm2 seront installs en bord de trou.
Les paramtres dimensionnels sont:
-larg = 50 mm
-long = 200 mm
-R1 = 5 mm
-R2 = 10 mm
Les paramtres de maillage sont:
-nR = 20 (nombre dlments en bord de trou)
-ny = 5 (nombre dlments entre les 2 cercles)
-elem = 3 mm (taille des lments en zone 2)
23
M.E.F.
M SUDRE
Pour paramtrer le problme, il faut crer le fichier p3prolog.pcl dans lequel sont
dclares les variables:
global real larg=50.,long=200.,R1=5.,R2=10.
global real elem=3.
global integer nR=20,ny=5
Enregistrer ce fichier dans le rpertoire dans lequel dmarre Patran.
2 Paramtrage de la gomtrie
Lancer lenregistrement dune session appele "geo".
Crer le point de coordonnes [0. , -R1 , 0.] puis le cercle de rayon R1.
Avec la mme mthode, crer le point [0. , -R2 , 0.] et le cercle de rayon R2.
Puis crer les points 5 8:
24
M.E.F.
M SUDRE
25
M.E.F.
M SUDRE
fichier "geo.ses":
STRING asm_create_grid_xyz_created_ids[VIRTUAL]
STRING sgm_sweep_curve_rev_created_ids[VIRTUAL]
STRING asm_line_2point_created_ids[VIRTUAL]
STRING sgm_surface_2curve_created_ids[VIRTUAL]
STRING sgm_surface_trimmed__created_id[VIRTUAL]
STRING sgm_create_curve_ch_created_ids[VIRTUAL]
INTEGER sgm_create_curve_cha_segment_id
asm_const_grid_xyz( "1", "[0 `-R1` 0]", "Coord 0", @
asm_create_grid_xyz_created_ids )
asm_const_grid_xyz( "2", "[0 `-R2` 0]", "Coord 0", @
asm_create_grid_xyz_created_ids )
sgm_const_curve_revolve( "1", "Coord 0.3", 180., 0., "Coord 0", "Point 1", @
sgm_sweep_curve_rev_created_ids )
sgm_const_curve_revolve( "2", "Coord 0.3", 180., 0., "Coord 0", "Point 2", @
sgm_sweep_curve_rev_created_ids )
asm_const_grid_xyz( "5", "[0 `-long/2` 0]", "Coord 0", @
asm_create_grid_xyz_created_ids )
asm_const_grid_xyz( "6", "[0 `long/2` 0]", "Coord 0", @
asm_create_grid_xyz_created_ids )
asm_const_grid_xyz( "7", "[`larg` `-long/2` 0]", "Coord 0", @
asm_create_grid_xyz_created_ids )
asm_const_grid_xyz( "8", "[`larg` `long/2` 0]", "Coord 0", @
asm_create_grid_xyz_created_ids )
point_label( TRUE )
asm_const_line_2point( "3", "Point 3", "Point 4", 0, "", 50., 1, @
asm_line_2point_created_ids )
asm_const_line_2point( "4", "Point 4", "Point 6", 0, "", 50., 1, @
asm_line_2point_created_ids )
asm_const_line_2point( "5", "Point 6", "Point 8", 0, "", 50., 1, @
asm_line_2point_created_ids )
asm_const_line_2point( "6", "Point 8", "Point 7", 0, "", 50., 1, @
asm_line_2point_created_ids )
asm_const_line_2point( "7", "Point 7", "Point 5", 0, "", 50., 1, @
asm_line_2point_created_ids )
asm_const_line_2point( "8", "Point 5", "Point 2", 0, "", 50., 1, @
asm_line_2point_created_ids )
asm_const_line_2point( "9", "Point 2", "Point 1", 0, "", 50., 1, @
asm_line_2point_created_ids )
sgm_const_surface_2curve( "1", "Curve 1", "Curve 2", @
sgm_surface_2curve_created_ids )
surface_label( TRUE )
trim_surface_label( TRUE )
sgm_create_curve_chain_v1( "10", "Curve 2 8:4:-1", TRUE, @
sgm_create_curve_cha_segment_id, sgm_create_curve_ch_created_ids )
sgm_create_surface_trimmed_v1( "2", "Curve 10", "", "", TRUE, TRUE, TRUE, @
TRUE, sgm_surface_trimmed__created_id )
ga_display_lines_set( "general", 2 )
26
M.E.F.
M SUDRE
3 Paramtrage du maillage
Lancer lenregistrement dune session appele "mesh".
Dfinir un Mesh Seed sur les arcs de cercle (nR) et sur le segment entre les 2 cercles (ny).
27
M.E.F.
M SUDRE
28
M.E.F.
M SUDRE
29
M.E.F.
M SUDRE
contrainte (MPa)
()
30
M.E.F.
M SUDRE
Model validity
Checks
1 Etude modale structure libre.
La structure tant libre, on excute une sol103 et on vrifie la prsence de 6 modes
rigides (frequences infrieures ou gales 1.0E-04 Hz)
5 Calcul de la masse.
Dans le BULK, ajouter la carte:
PARAM,GRDPNT,1234
numro du noeud
o sont calculs
les moments dinertie
31
M.E.F.
M SUDRE
6 Mesh Check.
Cette vrification concerne le maillage.
Loption Show Only Free Edges permet de visualiser si les condensations (quivalence)
sont bien ralises.
avant "equivalence"
aprs "equivalence"
32
M.E.F.
M SUDRE
7 Geometry Check.
Cette vrification concerne la forme des lments Quad:
Aspect Ratio:
Warp Angle:
a
b
h
a
un message si
b
>5
a
un message si
h
> 5%
a
Taper:
Skew Angle:
3
A1
A2
A3
un message si <30
Aire
33
A4
M.E.F.
M SUDRE
8 Vrification du rsultat.
Aprs rsolution du systme:
[K].{U} = {F}
le terme: {F} = {F} - [K].{U}
par la formule:
{U}t.{F}
{U}t.{F}
34
M.E.F.
M SUDRE
1 Mthode de Newton
Soit par exemple calculer:
f ( x ) dx
1
f(x)
-1
Ide: considrons 4 points P1, P2, P3, P4 uniformment distribus. La fonction f(x)
prend respectivement en ces points les valeurs f1 , f2, f3, f4.
f(x)
f2
f4
f3
f1
P3
P2
P1
P4
a0 + a1 x1 + a2 x2 + a3 x3
1
et de calculer
1
1
( a + a x + a x + a x ) dx
1
35
2
= 2 a 0 + --- a 2
3
M.E.F.
M SUDRE
2 Mthode de Gauss
Ide: considrons 2 points P1, P2 distribus symtriquement, leur position tant repre par .
La fonction f(x) prend la valeur f1 en P1 et la valeur f2 en P2 .
Les 2 points P1 et P2 sont nomms points de GAUSS.
f(x)
f1
f2
P2
P1
a0 + a1 x1 + a2 x2 + a3 x3
f1 + f2
Il vient:
Or
= 2 ( a0 + a2 ) .
2
( a + a x + a x + a x ) dx
1
1
= ------- pour obtenir:
3
1
= 2 a 0 + --- a 2
1
1
- + f -------
f ( x ) dx f 1 + f 2 = f -----
3
3
1
36
M.E.F.
M SUDRE
Postprocessing
1 Prsentation dun trac Fringe.
Cest un trac pour lequel une zone de couleur reprsente un niveau de valeur du
rsultat reprsent.
Slectionner le subcase.
Slectionner le rsultat tracer.
Current Viewport: par dfaut, tous les lments du modle qui sont affichs lcran
sont concerns.
sinon:
Elements: slectionner des lments.
Groups: slectionner un ou des groupes contenant des lments.
Materials: slectionner des lments par le matriau associ.
Properties: slectionner des lments par la proprit associe.
Element Types: slectionner des lments par leur type.
Free Faces: par dfaut, seules les faces libres sont concernes.
Faces: toutes les faces des lments slectionns sont concernes.
Free Edges: seuls les bords libres sont concerns.
Edges: tous les bords sont concerns.
4 Options daffichage
37
M.E.F.
5 Options de trac
M SUDRE
Coordinate Transformation:
Les composantes de vecteurs (dplacements) ou de tenseurs (contraintes)
afficher peuvent tre exprimes dans diffrents repres.
Loption par defaut As Is correspond aux repres locaux des lments (comme dans
le fichier f06). Loption CID permet de choisir nimporte quel repre cr par lutilisateur.
Scale Factor:
Permet dappliquer un coefficient multiplicateur aux rsultats affichs.
Filter Values:
Permet de filtrer les valeurs infrieures ou suprieures un certain seuil ou
bien comprises dans un intervalle donn.
Averaging Domain:
Cette option concerne les rsultats "Elments" qui doivent tre affichs aux
noeuds. Plusieurs valeurs se trouvent attaches un mme noeud et il faut choisir
une option pour moyenner ces quantits.
Par dfaut toutes les valeurs affectes au noeud sont moyennes. Mais il est possible
de se limiter uniquement aux lments qui partagent le mme matriau, la mme
proprit ou le mme type.
Averaging Method:
Lorsque des rsultats ne sont pas directement les composantes de vecteurs
(dplacements) ou de tenseurs (contraintes) mais proviennent dun calcul partir de
ces composantes, le rsultat peut varier si:
-on calcule dabord puis on moyenne ensuite,
-on moyenne dabord puis on clalcule ensuite.
Lutilisateur doit dcider de la formule qui lui semble la plus juste.
Extrapolation Method:
Les rsultats qui doivent tre affects aux noeuds sont parfois connus aux
points dintgration. Il existe diffrentes mthodes pour extrapoler ces valeurs des
points de Gauss vers les noeuds.
-Shape Function: la fonction de forme utilise pour lintgration sert
effectuer lextrapolation aux noeuds,
-Average: la moyenne aux points de Gauss est calcule et affecte
chaque noeud,
-Centroid: le rsultat calcul au centre de llment est affect chaque
noeud.
6 Options danimation
M.E.F.
M SUDRE
projet PATRAN-NASTRAN
Cette voilure a fait lobjet dun essai dynamique suivant Z (balayage en frquence):
Z
Y
X
3 modes ont t obtenus pour des frquences de 14.75 Hz, 66.75 Hz et 93.25 Hz.
M.E.F.
M SUDRE
root section
Main spar
Rib
1500.
Secondary spar
175.
Matriaux:
40
250.
250.
rib7
200.
rib6
175.
rib5
150.
rib4
175.
rib3
150.
rib2
150.
rib1
root section
150.
M.E.F.
M SUDRE
Gomtrie:
Rib
123.8
45.33
39.70
52.50
58.32
86.30
R10.0
172.6
198.8
2.4
2.4
Rib + Connectors
8.
6.
18.
6.
24.
8.
20.
20.
20.
Spar
Connector
Proprits:
41
20.
M.E.F.
M SUDRE
Corner connectors
30.
axe du rivet
noeuds dpendants
noeud indpendant
42
M.E.F.
M SUDRE
Les nervures (rib1 rib7) seront dupliques partir de la nervure rib0 par la commande Group/Transform/Translate.
43
250
250
rib7
200
rib6
175
rib5
150
rib4
175
rib3
150
rib2
150
rib1
rib0
150
M.E.F.
M SUDRE
root section
150
1
150
2
150
175
4
150
c
se
d
on
y
ar
ar
sp
ain
b
we
sp
b
we
r
a
175
200
250
250
44
M.E.F.
M SUDRE
4 Modlisation de la peau.
Dfinir le profil de la voilure (profil.bdf) partir dlments rods reliant les noeuds du
contour de la rib0 puis gnrer la peau par la commande Sweep.
Respecter la concordance des noeuds entre longerons et peau suivant x.
La peau et des longerons seront lis par des RBE2 (pour simuler le collage).
RBE2
RBE2
45
M.E.F.
M SUDRE
Tip.
noeud matre
centre de gravit du Tip
[1520. -175. 0]
170g
root section
30.
60.
30.
7 Calcul.
Caractriser les premiers modes et comparer les rsultats obtenus avec les valeurs exprimentales. Excuter des sol108 et sol111. Conclusions.
46
M.E.F.
M SUDRE
Annexes
47