NF04

1
NF04 Pratique des EF avec ANSYS
Benjamin Hagège, Mabrouk Ben Tahar Département
Génie Mécanique
NF04 UTC
2
Plan du cours (2 CM 3 TD)
0) Présentation générale dANSYS
I) TD1 plaque en conduction thermique
II) TD2 poutre en traction et flexion
III) TD3 barre en thermoélectricité
BONUS Macro VBA dimport sous Powerpoint
NF04 UTC
3
0) Présentation générale dANSYS
NF04 UTC
4
ANSYS, code n1 mondial
NF04 UTC
5
ANSYS, code n1 mondial
Entreprises en France -ALSTOM tout ce qui est
vibratoire -AREVA solution unique mondiale, du
coup EdF aussi pour l'EPR -EADS-ASTRIUM impact
sur les satellites -SAFRAN montée en puissance
de Workbench, perte de vitesse d'Abaqus -Certains
centres ST-GOBAIN impact sur les
parebrises -ARCELOR pour LS-DYNA (distribué par
ANSYS et inclu dans le code) -VINCI, BOUYGUES
CONSTRUCTION notamment pour le projet à la
Défense de la nouvelle tour de 260m de
haut -VEOLIA -SUEZ -DCN -DGA -ALKAN
(missiles) -THALES -SDBM (Monaco) -ITER, CERN et
bientôt le CEA -Renault Sport F1
NF04 UTC
6
ANSYS, code n1 mondial
NF04 UTC
7
Travailler avec ANSYS
Ouverture dANSYS CLASSIC (Mechanical
APDL) -Ouvrir le launcher -Définir sa licence
(Ansys Academic Introductory) et son répertoire
de travail
Règles de base -Réserver sur le HD un espace,
permanent, pour la sauvegarde des fichiers
(NF04-TD1) et un espace, temporaire, pour le
calcul (0Calculs) -8 caractères pour un nom, pas
despaces, pas daccents -ANSYS ne fait pas de
différences entre majuscules et miniscules Ee !
Méthodologie de mise en oeuvre -Cliquer dans le
Graphic User Interface sur une fonction -Ouvrir
le Session Editor pour construire son script APDL
(dans un fichier .dat)
Structurer son job -Paramètres définition de la
session -Mise en donnée EF type dEF, real
dEF, matériaux, CAO, maillage -Analyse CL (T,
U, ROT, V), chargements (Q, F, M, I), type de
calcul, sorties -Post-traitement le plus
difficile car dépend du problème ! Des
iso-contours, des champs de tenseur et des bilans
textuels
NF04 UTC
8
Travailler avec ANSYS
Mauvaise approche des EF -Vouloir cliquer
interactivement pour tout faire et ignorer le
prompt (mode du clickngo) -Refuser dadmettre
que faire des éléments finis, cest avant tout
coder un job Bonne approche des EF -Coder ses
jobs à lancienne avec un authentique
préprocesseur, à savoir un éditeur de texte
performant (EditPlus, UltraEdit, PSPad) -Ne se
servir du GUI que pour définir ses
vues Apprendre à scripter lAPDL -Le .log
contient tout ce qui exécuté (clics ou commandes
tapées dans le prompt) -Pour apprendre cliquer
pour effectuer laction désirée dans le GUI,
ouvrir le .log (ou le session editor), copier la
partie APDL correspondante et la stocker dans ses
tablettes pour la réutiliser -Se servir de laide
en tapant help, nom2commande
NF04 UTC
9
Organisation du GUI 1
Modules dANSYS CLASSIC pour NF04 -/PREP7
pré-processeur, /SOLU CLChargements, analyse,
/POST1 post-processeur
NF04 UTC
10
Organisation du GUI 2
Output Window et Pan/Zoom/Rotate -Output Window
contrôle des commandes effectuées et affichage
des commentaires et résultats textuels (bilans)
NF04 UTC
11
Help
Aide dANSYS sur la commande ou lélément
xxxx -Dans le prompt taper help,xxxx -Laide
comporte de nombreux exemples de références
Verification Manual
Exemples de vérification -VM1 Statically
Indeterminate Reaction Force Analysis (S.
Timoshenko, Strength of Material, Part I,
Elementary Theory and Problems, 3rd Edition, D.
Van Nostrand Co., Inc., New York, NY, 1955, pg.
26, problem 10) -help,1
NF04 UTC
12
vm1.dat en détail
/COM,ANSYS MEDIA REL. 11.0 (10/27/2006) REF.
VERIF. MANUAL REL. 11.0 /VERIFY,VM1 /PREP7
/TITLE, VM1, STATICALLY INDETERMINATE REACTION
FORCE ANALYSIS ET,1,LINK1 ! TYPE DELEMENT
NUMEROTE 1 R,1,1 ! CROSS SECTIONAL AREA
(ARBITRARY) 1 MP,EX,1,30E6 ! MATERIAU DE
MODULE DYOUNG 30E6 (psi) N,1 ! NŒUD 1 EN
(0,0) N,2,,4 ! NŒUD 2 EN (0,4) N,3,,7 ! NŒUD 3
EN (0,7) N,4,,10 ! NŒUD 4 EN (0,10) E,1,2 !
DEFINE ELEMENTS (ENTRE LES NŒUDS 1 ET
2) EGEN,3,1,1 ! GENERER 3 ELEMENTS A PARTIR DE
LELEMENT 1 PAR INC. DE 1 D,1,ALL,,,4,3 !
BOUNDARY CONDITIONS AND LOADING
(DDISPLACEMENTqinc) F,2,FY,-500 ! APPLIQUER
UNE FORCE AU NŒUD 2 SELON Y DE -500 (lb)
(FFORCEFext) F,3,FY,-1000 ! APPLIQUER UNE
FORCE AU NŒUD 3 SELON Y DE -1000 (lb) /SOLU
ANTYPE,STATIC ! STATIC ANALYSIS OUTPR,BASIC,1
! SORTIR LES GRANDEURS BASIQUES OUTPR,NLOAD,1 !
SORTIR LES CHARGEMENTS NODAUX SOLVE ! LANCER LE
CALCUL /POST1 NSEL,S,LOC,Y,10 ! SELECTION SUR
LE PLAN Y10 DES NOEUDS FSUM ! SOMMER LES
EFFORTS ET MOMENTS SELON X, Y, Z GET,REAC_1,FSUM,
,ITEM,FY ! ATTRAPER LA VALEUR DE LA SOMMATION
PRECEDENTE NSEL,S,LOC,Y,0 FSUM
GET,REAC_2,FSUM,,ITEM,FY
NF04 UTC
13
vm1.dat en détail - amélioration
Post-traitement -le post-traitement répond
uniquement à la question posée quelles forces
de réactions pour ce problème hyperstatique
? -Obtenir les forces de réactions dans un
fichier rapport.post /POST1 NSEL,S,LOC,Y,10
FSUM GET,REAC_1,FSUM,,ITEM,FY NSEL,S,LOC,Y,0
FSUM GET,REAC_2,FSUM,,ITEM,FY
/OUTPUT,rapport,post /COM !!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! /COM !!!! CALCUL
DES REACTIONS !!!!!! /COM !!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!! /COM REACTION HAUT
(lb) REAC_1 /COM REACTION BAS (lb)
REAC_2 /OUTPUT
NF04 UTC
14
Calcul en pratique
2 modes dutilisation -Session intéractive mise
en donnée, vérification, petit calcul. On
travaille par copier/coller du .dat vers le
prompt dANSYS -Batch pas dinterface
graphique, gros calculs validés à la suite. On
exécute un programme préparé
Comment batcher un job ANSYS ? -2 fichiers sont
nécessaires sous Windows !EXEC.bat et
Liste.cmd -!EXEC.bat fichier à exécuter qui va
appeler le .cmd. Il contient _at_echo off set
ANSYS130_PRODUCTaa_t_a (anciennement aa_t_i) set
ANS_CONSECYES call Z\0Calculs\Liste.cmd -
Liste.cmd fichier qui contient la liste des
jobs et lemplacement de lexécutable ANSYS CD
Z\0Calculs\ "C\Program Files\Ansys
Inc\v130\ANSYS\bin\Intel\ansys130" -b -i
0NF04_ANSYS-1.dat -o 0NF04_ANSYS-1.out "C\Program
Files\Ansys Inc\v130\ANSYS\bin\Intel\ansys130"
-b -i 0NF04_ANSYS-2.dat -o 0NF04_ANSYS-2.out
NF04 UTC
15
Calcul en pratique
Savoir sorganiser -Existence de 3 répertoires de
travail sur Z/ -Répertoire de
sauvegarde PERMANENT pour sauvegarde des
jobs -Répertoire de travail TEMPORAIRE pour
édition des jobs ou run intéractif -Répertoire de
calcul TEMPORAIRE pour exécution des jobs en
batch
NF04 UTC
16
I) Exemple de TD plaque en conduction thermique
NF04 UTC
17
Physique du problème
Thermique -Type de transferts conduction,
convection (libre/forcée), rayonnement
? -Grandeurs à analyser T (C), q (W/m²), F
(W) gt graphiques iso-contours, champs de
vecteur gt feuille de calcul bilan des flux
NF04 UTC
18
Job paramètres 1
!! SYSTEME mm, C, mW !

TA44.51151 !TEMP ANALYTIQUE POINT
A TB36.40567 !TEMP ANALYTIQUE POINT
B !!PROBLEME a1000. !DIMENSION X
PLAQUE b500. !DIMENSION y PLAQUE k1. !COND
UCTIVITE THERMIQUE BETON !! !! !!MAILLAGE !
!! !! !!POST NOMJOBNF04_ANSYS-1 !NOM DU
JOB (CHAINE DE CARACTERES) ! SFACTOR10
!FACTEUR DE ZOOM POUR LES VECTEURS
Remarque en violet et avec ! sont indiqués
les endroits où il faudra modifier le code en TD
afin de réaliser lensemble des calculs et des
post-traitements automatiquement.
NF04 UTC
19
Job mise en donnée 2
!
/NOPR !PAS DE
COMMENTAIRES DANS L'OUTPUT WINDOW /TITLE,NOMJOB
!TITRE SUR LA FENETRE D'AFFICHAGE /FILENAME,NOMJ
OB !NOM DES FICHIERS GENERES ! ! ! ! ! ! !

/PREP7 !ENTREE
DANS LE PREPROCESSEUR ! ET,1,PLANE55 !TYPE
D'EF AVEC L'ORDRE, ICI BILINEAIRE ! !
MP,KXX,1,k !MATERIAU ! RECTNG,0,a/2.,0,b, !C
AO ! ! ! !MAILLAGE MSHAPE,0,2D !MeshSH
APE, ici QUADRANGLE ! MSHKEY,1 !MAILLAGE
REGLE AESIZE,ALL,50, !AreaElementSIZE50mm SUR
TOUTES LES SURFACES ! AMESH,ALL !AreaMESH SUR
TOUTES LES SURFACES SELECTIONNEES
NF04 UTC
20
Job CL, chargement et résolution 3
!
/SOL !ENTREE
DANS LE SOLVEUR ! ANTYPE,0 !ANALYSE REGIME
PERMANENT 0STATIC ! nsel,,loc,X !SELEC BORD
X- nsel,a,loc,Y,b !RAJOUTER BORD
Y nplot !AFFICHER LES NOEUDS D,ALL,TEMP,0. !
CL SUR TOUS LES NOEUDS SELECTIONNE T0C !
nsel,,loc,Y !SELEC BORD Y- nplot D,ALL,TEMP,10
0. ! allsel !TOUT SELECTIONNER solve !LANCE
R LE CALCUL SUR LA SELECTION ACTIVE
NF04 UTC
21
Job post-traitement graphique 4
/POST1 !ENTREE DANS LE POSTPROCESSEUR !!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!! !Inverser la
video! /RGB,INDEX,100,100,100, 0 /RGB,INDEX,
80, 80, 80,13 /RGB,INDEX, 60, 60, 60,14
/RGB,INDEX, 0, 0, 0,15 !!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!! /SHOW,PNG,,0 !DEBUT CAPTURES
D'ECRAN /pbc,RFORCE,,1 !AFFICHER LES FORCES DE
REACTION Eplot !AFFICHER LES ELEMENTS DE LA
SELEC. ACTIVE /pbc,RFORCE,,0 !NE PAS AFFICHER
LES FORCES DE REACTION /CVAL,ALL,20,30,36.40567,40
,44.51151,60,80,100 !CODE DE COULEUR plns,temp
!CHAMP DE TEMPERATURE /CVAL !RESET DU
CODE DE COULEUR ples,tf,sum !FLUX DE
CONDUCTION ELEM ISOCONTOURS plns,tf,sum
!FLUX DE CONDUCTION NODAL ISOCONTOURS ! /VSCAL
E,1,SFACTOR,0 !FACTEUR DE ZOOM DES
VECTEURS plvect,tf,,,,VECT,ELEM,ON !FLUX DE
CONDUCTION CHAMP DE TENSEUR /SHOW,CLOSE !FIN
CAPTURES D'ECRAN ! ! /SYS,MKDIR
1RESULTATS !CREER LE REPERTOIRE
1RESULTATS /SYS,MOVE .PNG 1RESULTATS/ !Y
DEPLACER TOUTES LES CAPTURES DECRAN
NF04 UTC
22
Job post-traitement des bilans 5.1
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !! DONNES
MODELE !! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ALLSEL
GET, NN, NODE, 0, COUNT !ATTRAPER LE NOMBRE DE
NOEUDS GET, NE, ELEM, 0, COUNT !ATTRAPER LE
NOMBRE D'ELEMENTS ! ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!! !! TEMPERATURES CIBLES
!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Lab
NANODE(a/2,b/2,0.) !ATTRAPER LE LABEL DU NOEUD
EN A LabNBNODE(a/4,b/2,0.) !ATTRAPER LE LABEL
DU NOEUD EN B GET, TAMEF, NODE, LabNA,
TEMP !ATTRAPER LA TEMPERATURE DU NOEUD EN A GET,
TBMEF, NODE, LabNB, TEMP !ATTRAPER LA TEMPERATURE
DU NOEUD EN B
NF04 UTC
23
Job post-traitement des bilans 5.2
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !! BILAN DES FLUX
!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !! nsel,,loc,X
!SELECTION DES NOEUDS DU PLAN
X0 FSUM !SOMMATION DES FLUX NODAUX GET, CXM,
FSUM, 0, ITEM, HEAT !ATTRAPER LA SOMME DES FLUX
NODAUX PRECEDENTE !! nsel,,loc,Y,b !SELECTION
DES NOEUDS DU PLAN Yb nsel,u,loc,X !UNSELECT
DES NOEUDS DU PLAN X0 FSUM GET, CYP, FSUM, 0,
ITEM, HEAT !! nsel,,loc,Y nsel,u,loc,X FSUM GET,
CYM, FSUM, 0, ITEM, HEAT !! nsel,all !SELECTION
DE TOUS LES NOEUDS FSUM GET, PHITOT, FSUM, 0,
ITEM, HEAT !FLUX TOTAL DE REACTION
NF04 UTC
24
Job sortie des bilans dans un fichier 6
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !! SORTIE FICHIER
!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!! /OUTPUT,rapport_
NF04-1,post,,APPEND !DIRIGER LA SORTIE TEXTE VERS
LE FICHIER .post /COM /COM JOBNOMJOB /COM
!!!!!!!!!!!!!!!!!!! /COM !! MODELE !! /COM
!!!!!!!!!!!!!!!!!!! /COM NB NOEUDSNN /COM NB
ELEMNE !!!! ! /COM !!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!! /COM !! TEMPERATURES !! /COM
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! /COM ERATAMEF
! /COM ERBTBMEF ! /COM !!!!!!!!!! /COM
!! FLUX !! /COM !!!!!!!!!! /COM FLUX YCYP mW
! /COM FLUX X-CXM mW ! /COM FLUX
Y-CYM mW ! /COM /OUTPUT !REDIRIGER LA
SORTIE TEXTE VERS LOUTPUT WINDOW
NF04 UTC
25
Champ de température
NF04 UTC
26
Forces de réaction sur un maillage raffiné
NF04 UTC
27
Densité de flux de conduction thermique
NF04 UTC
28
Densité de flux de conduction thermique très
raffiné
Problème physique ?
NF04 UTC
29
TD 1 cylindre en conduction
FZ
FRadial?
FZ-
NF04 UTC
30
Erreurs à ne pas faire !
Savoir présenter des résultats physiques -Courbes
discrètes -Chiffres significatifs
Savoir rendre un rapport -Introduction
conclusions -Livrer dans le délai imparti suivant
la procédure
NF04 UTC
31
II) TD2 poutre en flexion
NF04 UTC
32
Physique du problème
Mécanique des structures minces -Milieux
curvilignes efforts modélisés par des torseurs
(MQ01) -Grandeurs à analyser U et ROT (mm et
), sij (MPa), RF et RM (N et Nmm) gt
graphiques iso-contours et diagrammes de
cohésion gt feuille de calcul bilan des
efforts, flèche max, s max
NF04 UTC
33
Physique du problème
3 cas -CAS 2  encastré-libre effort tranchant
FZ à lextrémité libre -CAS 2  encastré-libre
effort tranchant FY à lextrémité libre -CAS 4 
appuyé à ses deux extrémités effort tranchant
FZ au centre
NF04 UTC
34
Orientation dune section délément 188 aligné
avec Z
ZlocXlocYloc
Dans le plan (X,Y) si Xloc nest pas colinéaire à
Y
Zloc
Yloc
Xloc
Vue dune section par SECPLOT sens contraire
de Xloc
NF04 UTC
35
Théorie cinématique dune poutre droite épaisse
?
NF04 UTC
36
Torseur de cohésion ?
NF04 UTC
37
Cinématique et torseur de cohésion dune section
dans le plan
HPP  zoom x100  au moins
NF04 UTC
38
Job paramètres 1
!! SYSTEME mm, N !
!!
PROBLEME L5000. !LONGUEUR POUTRE EN
mm FEXT10000. !N (1T ici) !!GEOMETRIE DE LA
SECTION W1100. !LARGEUR SEMELLES W2W1 W3200
. !HAUTEUR AME t13. !EPAISSEUR t2t1 t3t
1 !!MATERIAU Young190000. !MPa NU0.3
!!MAILLAGE NBELEM2 !NOMBRE ELEMENTS PAIR
POUR POINT C !!POST NOMJOB'NF04_ANSYS-2' CAS'4
!
NF04 UTC
39
Job mise en donnée 2
!
!TITRE SUR LA
FENETRE D'AFFICHAGE /TITLE,NOMJOB_CASCAS /FI
LENAME,NOMJOB_CASCAS !NOM DES FICHIERS
GENERES ! ! ! ! ! ! !

/PREP7 !ENTREE DANS LE
PREPROCESSEUR ! !!POUTRE !TYPE D'EF AVEC
L'ORDRE, ICI BILINEAIRE ET,1,188 KEYOPT,1,3,2 !FO
RMULATION ISOPARAM(0) OU SUPERPARAM
(2) KEYOPT,1,7,1 !SMAX AUX POINTS
D'INTEGRATION ! ! MP,EX,1,YOUNG !MATERIAU
MP,NUXY,1,nu !MATERIAU ! !
K,1,0.0.,0 !CAO KEYPOINT 1 K,2,L,0.,0. !CAO
KEYPOINT 2 L,1,2 !CAO LIGNE DE 1-gt2 lplot
NF04 UTC
40
Job définition des sections et maillage 3
!MAILLAGE ! !ASSIGNER LES SECTIONS AUX
KEYPOINTS 1 ET 2 SECTYPE,1,BEAM,I ! SECDATA,W1
,W2,W3,t1,t2,t3 ! !! SECTYPE,10,TAPER !TAPERE
D BEAM SECDATA,1,0.,0.,0. !STARTING LOCATION OF
TAPERED BEAM SECDATA,1,L.,0.,0. !ENDING
LOCATION OF TAPERED BEAM !! !! !! TYPE,1 MAT,1
SECNUM,10 !ACTIVER TYPE 1, MAT 1 ET SECTION
10 LESIZE, 1, L/NBELEM, !L/NBELEMTAILLE DE
MAILLE lmesh,all !MAILLER TOUTES LES LIGNES
SELECTIONNEES !!!!!!!!!!!!
NF04 UTC
41
Job CL, chargement et résolution 4
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !Inverser la
video! /RGB,INDEX,100,100,100, 0 /RGB,INDEX,
80, 80, 80,13 /RGB,INDEX, 60, 60, 60,14
/RGB,INDEX, 0, 0, 0,15 !!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!! !Definir une vue
sympa (à partir du sessions editor) /USER, 1
/VIEW, 1, 0.565040617968 , -0.611926071170
, 0.553421704913 /ANG, 1, -61.6518891670
/REPLO /AUTO,1 /REP,FAST !!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! /eshape,1,1 !AFFICHER LA POUTRE AVEC SA
GEOMETRIE 3D /SHOW,PNG,,0 !DEBUT CAPTURES
D'ECRAN secplot,1 !AFFICHER LA SECTION
1 eplot !AFFICHER LES ELEMENTS /SHOW,CLOSE !F
IN CAPTURES D'ECRAN
NF04 UTC
42
Job CL, chargement et résolution 5
!
/SOL !ENTREE
DANS LE SOLVEUR ! ANTYPE,0 !ANALYSE REGIME
PERMANENT 0STATIC ! !
! CLCHARGEMENT !
NSEL,,LOC,X !CL ! D,ALL,UX ! D,ALL,UY
! D,ALL,UZ ! D,ALL,ROTX ! !! NSE
L,,LOC,X,L !CL ! D,ALL,ROTZ ! D,ALL,UZ
! !! NSEL,,LOC,X,L/2 !CHARGEMENT FORCE
CONCENTREE ! F,ALL,FZ,FEXT ! !
! ! ! allsel !TOUT
SELECTIONNER solve !LANCER LE CALCUL SUR LA
SELECTION
NF04 UTC
43
Job affichage des CL et chargements 6
!
/POST1
!ENTREE DANS LE POSTPROCESSEUR ! ! ! ! /eshape
,1,1 !AFFICHER LA POUTRE AVEC SA GEOMETRIE
3D ! ! ! /PBC,U,,1 !AFFICHER
CLCHARGEMENTS /PBC,ROT,,1 /PBC,FORCE,,1 /PBC,MOME
NT,,1
NF04 UTC
44
Job générer les diagrammes de cohésion 7
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!
TORSEUR DE COHESION !! !!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!! !!ECRITURE DES TABLES aux noeuds
I et J DE CHACUN DES EF !N ETABLE,FX_I,SMISC,1 ETA
BLE,FX_J,SMISC,14 !FY ETABLE,FY_I,SMISC,6 ETABLE,F
Y_J,SMISC,19 !FZ ETABLE,FZ_I,SMISC,5 ETABLE,FZ_J,S
MISC,18 !M ETABLE,MX_I,SMISC,4 ETABLE,MX_J,SMISC,1
7 !MFY ETABLE,MY_I,SMISC,2 ETABLE,MY_J,SMISC,15 !M
FZ ETABLE,MZ_I,SMISC,3 ETABLE,MZ_J,SMISC,16
-Utilisation des tables définition des SMISC et
NMISC dans la librairie déléments (help,188
ici). -Afficher une table PRETAB,NOM_TABLE
NF04 UTC
45
Job post-traitement des diagrammes de cohésion
8
!!ISO-CONTOUR DES COMPOSANTES DU TORSEUR DE
COHESION EN REPERE LOCAL DE LA POUTRE !!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!! !Definir une vue sympa (à partir du
sessions editor) /VIEW,1,,,1 /ANG,1
/REP,FAST !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! /SHOW,PNG,,0 !DEB
UT CAPTURES D'ECRAN /TITLE,NOMJOB_CASCAS_N P
LLS,FX_I,FX_J,1. !1.FACTEUR DE
ZOOM /TITLE,NOMJOB_CASCAS_FY PLLS,FY_I,FY_J,1.
/TITLE,NOMJOB_CASCAS_FZ PLLS,FZ_I,FZ_J,1. /TI
TLE,NOMJOB_CASCAS_M PLLS,MX_I,MX_J,1. /TITLE,
NOMJOB_CASCAS_MFY PLLS,MY_I,MY_J,1. /TITLE,NOM
JOB_CASCAS_MFZ PLLS,MZ_I,MZ_J,1. /SHOW,CLOSE
!FIN CAPTURES D'ECRAN /TITLE,NOMJOB_CASCAS
NF04 UTC
46
Job post-traitement graphique 9
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!! !Definir une vue sympa (à
partir du sessions editor) /USER, 1 /VIEW,
1, 0.565040617968 , -0.611926071170 ,
0.553421704913 /ANG, 1, -61.6518891670
/REPLO /AUTO,1 /REP,FAST !!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! !!!!!!!!!!!!!!!!!!!!!!!!!!!!! !! ISO-CONTOURS
!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!! /GRAPHICS,POWER
!MODE POWERGRAPHICS POUR SECTION
3D /SHOW,PNG,,0 !DEBUT CAPTURES
D'ECRAN ples,s,eqv !VON MISES AUX
ELEMENTS ples,s,x !SXX AUX ELEMENTS plns,u,sum
!NORME DE U AUX NOEUDS /SHOW,CLOSE !FIN
CAPTURES D'ECRAN
NF04 UTC
47
Job calcul des erreurs sur la flèche et la
contrainte 10
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !! DONNES MODELE
!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ALLSEL GET,
NN, NODE, 0, COUNT !ATTRAPER LE NOMBRE DE
NOEUDS GET, NE, ELEM, 0, COUNT !ATTRAPER LE
NOMBRE D'ELEMENTS ! ! ! !!!!!!!!!!!!!!!!!!!!!!!
!! !! FLECHE MAX !! !!!!!!!!!!!!!!!!!!!!!!!!! LabN
BNODE(L/2,0.,0.) !ATTRAPER LE LABEL DU NOEUD EN
B ! GET, FLEMEF, NODE, LabNB, UZ !ATTRAPER
LA FLECHE NUMERIQUE DU NOEUD C ! !! GET, IYY,
SECP, 1, PROP, IYY !INERTIE DE LA SECTION
1 ! RIGIYOUNGIYY ! FLEANA(FEXTL
LL)/(48RIGI) !FLECHE ANALYTIQUE ! ! ! !!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !! CONTRAINTE
XX MAX !! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
GET,SIGMEF,SECR,ALL,S,X,MAX !PLUS GRANDE
CONTRAINTE XX DU MODELE !! SIGANA(FEXTL/4)/IYYW
3/2. !SXX ANALYTIQUE !
NF04 UTC
48
Job sorties dans un fichier 11
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !! SORTIE FICHIER
!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!! /NOPR /OUTPUT,ra
pport_NF04-2,post,,APPEND /COM JOBNOMJOB_CASCA
S /COM !!!!!!!!!!!!!!!!!!! /COM !! MODELE
!! /COM !!!!!!!!!!!!!!!!!!! /COM NB
NOEUDSNN /COM NB ELEMNE /COM
!!!!!!!!!!!!!!!!!!!!!!!!!!! /COM !! FLECHE MAX
!! /COM !!!!!!!!!!!!!!!!!!!!!!!!!!! /COM FLECHE
MAX FLEMEFmm /COM ERFLECHE(FLEMEF-FLEANA)/FL
EANA100 /COM !!!!!!!!!!!!!!!!!!!!!!!!! /COM !!
SIGMA MAX !! /COM !!!!!!!!!!!!!!!!!!!!!!!!! /CO
M SIGMAXX MAX SIGMEF MPa /COM
ERSIGMA(SIGMEF-SIGANA)/SIGANA100 /COM
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! /COM !!
BILAN DES EFFORTS !! /COM !!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!! prrs !SORTIR LES REACTIONS
NODALES CALCULEES
NF04 UTC
49
Job sorties dans un fichier 12
/COM /COM !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! /COM !!
DEPLACEMENTS !! /COM !!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!! PRNSOL,DOF !SORTIR LES INCONNUES NODALES
CALCULEES /COM /COM
/COM

/OUTPUT ! ! ! ! ! !!!!!
!!!!!!!!!!!!!!!!!!!!! !! RANGER LES RESULTATS
!! !!!!!!!!!!!!!!!!!!!!!!!!!! /SYS,MKDIR
1RESULTATS /SYS,MOVE .PNG 1RESULTATS/ SAVE
NF04 UTC
50
Torseur de cohésion cas 4
NF04 UTC
51
Différents cas pour le TD 2
Inertie dune section ? Elancement dune poutre ?
NF04 UTC
52
Inerties principales de la section en I
élancement l
NF04 UTC
53
Résolution analytique flowchart
Notion de déformation faire apparaître le
déplacement recherché
Diagramme de cohésion équilibre dune section
(MQ01)
Loi de comportement équation différentielle du
problème
CL et chargements intégration de la solution
NF04 UTC
54
Résolution analytique cas 3 flexion PURE
NF04 UTC
55
Résolution analytique cas 3 flexion PURE
Flèche max ?
Contrainte max ?
NF04 UTC
56
Résolution analytique cas 4 flexion SIMPLE
NF04 UTC
57
Résolution analytique cas 4 flexion SIMPLE
NF04 UTC
58
Von Misès cas 4
CL supplémentaires pour assurer lisostatisme
NF04 UTC
59
Contrainte normale cas 4
NF04 UTC
60
Flèche cas 4
Formulation superparamétrique
NF04 UTC
61
Contrainte de cisaillement transverse cas 4
Cisaillement transverse gt section réduite pour
les efforts tranchants
NF04 UTC
62
Flèche totale dune poutre en flexion simple
NF04 UTC
63
III) TD3 barre en thermoélectricité
NF04 UTC
64
Physique du problème
Aluminium
Thermoélectricité des circuits 3D -Milieux 3D
le courant électrique dissipe de la chaleur par
effet joule qui doit séchapper de la barre (CL
Dirichlet ou Cauchy) -Grandeurs à analyser E et
i (V/m et A/m²), T et F (C et W) gt
graphiques iso-contours, champs de vecteur gt
feuille de calcul bilan des flux thermique et
électrique
NF04 UTC
65
Système dunité
Thermique -Lmm -TC -kSIWC/mmWC/mm -L
es flux sont donc en mW
Electricité -Lmm (choix thermique) -PUI
mW (choix thermique) -UmV et IA gt
PmVAmW OK -RU/ImV/AmOhm gt
rmOhm.mm
Exemple de conversion dunité r30E-9
Ohm.m 30E-9 E3mOhm.E3mm 30E-3 mOhm.mm
NF04 UTC
66
Job paramètres 1
!! SYSTEME DUNITES ! THERMIQUE mm, C, mW gt
condmW/mm/CW/m/CSI ! ELECTRICITE mm, mV, A
gt PUImVAmW, RU/ImV/AmOhm ! ! !

!!CL VE100. !POTENTIEL
ENTREE EN mV VS0. !POTENTIEL SORTIE EN
mV TE100. !TEMPERATURE ENTREE EN
C TS100. !TEMPERATURE SORTIE EN
C !!GEOMETRIE L1000. !LONGUEUR BARRE EN
mm a100. !SECTION axb ba !!MATERIAU
ALUMINIUM RES30E-3 !RESISTIVITE EN
mOhm.mm COND231 !CONDUCTIVITE EN
mWC/mm !!MAILLAGE TM50 !TAILLE DE MAILLE en
mm !!POST NOMJOB'NF04_ANSYS-3' CAS'1'
NF04 UTC
67
Job mise en donnée 2
!
!TITRE SUR LA
FENETRE D'AFFICHAGE /TITLE,NOMJOB_CASCAS /FI
LENAME,NOMJOB_CASCAS !NOM DES FICHIERS
GENERES ! ! ! !
/P
REP7 !ENTREE DANS LE PREPROCESSEUR !
!!BARRE 3D ET,1,69 !H8 THERMOELECTRIQUE R,1 !
! MP,KXX,1,COND !MATERIAU
MP,RSVX,1,RES !MATERIAU ! !
BLOCK,0,L,-a/2,a/2,-b/2,b/2 !CAO ! ! !
!MAILLAGE AESIZE, ALL, TM !TAILLE DE MAILLE
TM SUR TOUTES LES SURFACES TYPE,1 MAT,1
REAL,1 !ACTIVER TYPE 1, MAT 1 ET REAL
1 VMESH,all !MAILLER VOLUME
NF04 UTC
68
Job CL, chargement et résolution 3
!
/SOL !ENTREE
DANS LE SOLVEUR ! ANTYPE,0 !ANALYSE REGIME
PERMANENT 0STATIC ! tunif,20 tref,20 toffst,273.
15 !
! CLCHARGEMENT ELECTRIQUE !
NSEL,,LOC,X !CL D,A
LL,VOLT,VE NSEL,,LOC,X,L !CL D,ALL,VOLT,VS !
!
CLCHARGEMENT THERMIQE !
NSEL,,LOC,X !CL D,ALL,TE
MP,TE NSEL,,LOC,X,L !CL D,ALL,TEMP,TS !
! ! allsel !TOUT
SELECTIONNER solve !LANCER LE CALCUL LA
SELECTION
NF04 UTC
69
Job afficher CL post-traitement graphique 4

/POST1 !ENTREE DANS LE POSTPROCESSEUR ! !
! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!
CLCHARGEMENTS !! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!! /SHOW,PNG,,0 !DEBUT CAPTURES
D'ECRAN /PBC,VOLT,,1 /PBC,TEMP,,1 eplot /PBC,VOLT,
,0 /PBC,TEMP,,0 /SHOW,CLOSE !FIN CAPTURES
D'ECRAN ! ! ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!! !! CHAMPS DE VECTEURS !! !!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!! /SHOW,PNG,,0 !DEBUT
CAPTURES D'ECRAN PLVECT,EF, , ,
,VECT,ELEM,ON,0 !champ electrique PLVECT,JS, , ,
,VECT,ELEM,ON,0 !champ deffet joule PLVECT,TF,
, , ,VECT,ELEM,ON,0 !flux de conduction thermique
NF04 UTC
70
Job post-traitement graphique 5
!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !! ISO-CONTOURS
!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!! plns,TEMP !TEM
PERATURES NODALES plns,VOLT !TENSIONS
NODALES /SHOW,CLOSE !FIN CAPTURES D'ECRAN
NF04 UTC
71
Job bilans 6
!!!!!!!!!!!!!!!!!!!!!! !! BILANS
!! !!!!!!!!!!!!!!!!!!!!!! ! !!EFFET
JOULE etable,DENS_JOULE,jheat !TABLE DE DENSITE
DE CHALEUR DE JOULE etable,VOL,volu !TABLE DES
VOLUMES ELEMENTAIRES smult,FLUX_JOULE,DENS_JOULE,V
OL !TABLE DES FLUX DE CHALEUR JOULE
ELEMENTAIRES ssum !SOMMER LES FLUX
ELEMENTAIRES get,JOULE,SSUM,FLUX_JOULE !ATTRAPER
LA VALEUR TOTALE DE L'EFFET JOULE ! !!REACTIONS
AUX CL NSEL,,LOC,X !SELECTIONNER NOEUDS
ENTREE FSUM !SOMMER LES REACTIONS get,FLUX_E,F
SUM,,ITEM,heat !ATTRAPER LA VALEUR TOTALE DU
FLUX DE CONDUCTION get,INT_E,FSUM,,ITEM,amps !AT
TRAPER LA VALEUR TOTALE DE L'INTENSITE !! NSEL,,LO
C,X,L !SELECTIONNER NOEUDS SORTIE FSUM !SOMME
R LES REACTIONS get,FLUX_S,FSUM,,ITEM,heat !ATTR
APER LA VALEUR TOTALE DU FLUX DE
CONDUCTION get,INT_S,FSUM,,ITEM,amps !ATTRAPER
LA VALEUR TOTALE DE L'INTENSITE ! !!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!! !! DONNES MODELE
!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ALLSEL GET,
NN, NODE, 0, COUNT !ATTRAPER LE NOMBRE DE
NOEUDS GET, NE, ELEM, 0, COUNT !ATTRAPER LE
NOMBRE D'ELEMENTS
NF04 UTC
72
Job sorties dans un fichier 7
!
!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!! !! SORTIE FICHIER !! !!!!!!!!!!!!!!!!!
!!!!!!!!!!!!! /NOPR /OUTPUT,rapport_NF04-3,post,,A
PPEND /COM JOBNOMJOB_CASCAS /COM
!!!!!!!!!!!!!!!!!! /COM !! MODELE !! /COM
!!!!!!!!!!!!!!!!!! /COM NB NOEUDSNN /COM NB
ELEMNE /COM
/COM TE
(C) TE /COM TS (C) TS /COM VE (mV)
VE /COM VS (mV) VS /COM !!!!!!!!!!!!!!!!!!!!!
!!!!! /COM !! BILAN !! /COM
!!!!!!!!!!!!!!!!!!!!!!!!!! prrs !SORTIR LES
REACTIONS NODALES CALCULEES
NF04 UTC
73
Job sorties dans un fichier 8
/COM !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! /COM !!
BALANCE EF !! /COM !!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!! /COM
/COM FLUX
THERMIQUE ENTREE (kW) FLUX_E1E-6 /COM FLUX
THERMIQUE SORTIE (kW) FLUX_S1E-6 /COM
DISSIPATION EFFET JOULE (kW) JOULE1E-6 /COM
BALANCE (mW) FLUX_EFLUX_S-JOULE /COM
INTENSITE ENTREE (kA) INT_E1E-3 /COM
INTENSITE SORTIE (kA) INT_S1E-3 /COM BALANCE
(A) INT_EINT_S /COM

/COM RESISTANCE TOTALE EF (mOhm)
(VS-VE)/INT_S /COM !!!!!!!!!!!!!!!!!!!!!!!!!!!!
/COM !! EF vs ANA !! /COM
!!!!!!!!!!!!!!!!!!!!!!!!!!!! /COM RESISTANCE
TOTALE ANALYTIQUE (mOhm) /COM CHAMP ELECTRIQUE
ANALYTIQUE (V/m) /COM EFFET JOULE ANALYTIQUE
(kW) /COM
/COM

/OUTPUT !!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !! RANGER LES
RESULTATS !! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!! /SYS,MKDIR 1RESULTATS /SYS,MOVE .PNG
1RESULTATS/ SAVE
NF04 UTC
74
Température cas 1
CL pour assurer lisostatisme
NF04 UTC
75
Potentiel cas 1
CL pour assurer lisostatisme
NF04 UTC
76
Champ électrique cas 1
IRE IRS 0
Rcircuit(VS-VE)/IRS
NF04 UTC
77
Densité de flux de conduction thermique cas 1
FjouleU²/R
FRE FRS Fjoule
NF04 UTC