CALCUL DE STRUCTURES PRESENTATION DE LA METHODES DES ELEMENTS FINIS

PREMIERS PAS DUN UTILISATEUR ABAQUS

L. PENAZZI

Dcembre 1999

ABAQUS

MANUEL UTILISATEUR

SOMMAIRE

1. LES PRINCIPALES PROCEDURES ET LES FICHIERS......................................1 2. MISE EN DONNEES ...............................................................................................2 2.1 Convention ...............................................................................................2 2.2 Systme dunits.......................................................................................2 2.3 Phases de mises en donnes .....................................................................2 2.4 Mots cls ABAQUS .................................................................................3 2.4.1 Phase 1 : Gnration du maillage (noeuds, lments) ...............3 2.4.2 Phase 2 : Proprits des matriaux (mcanique, physique).......4 2.4.3 Phase 3 : Conditions aux limites et conditions initiales ............5 2.4.4 Phase 4 : Etape de calcul ...........................................................5 3. TYPE D'ELEMENTS................................................................................................9 3.1 Elements 1D solides .................................................................................9 3.1.1 Gomtrie et degrs de libert ...................................................9 3.1.2 Chargement................................................................................9 3.2 Elements 2D solides ..................................................................................10 3.2.1 Gomtrie et degrs de libert ...................................................10 3.2.2 Chargement................................................................................11 3.2.3 Rsultats.....................................................................................11 3.3 Elements 2D axisymtriques....................................................................12 3.3.1 Gomtrie et degrs de libert ...................................................12 3.3.2 Chargement................................................................................13 3.2.3 Rsultats.....................................................................................13 3.4 Elements poutres ......................................................................................14 3.4.1 Gomtrie et degrs de libert ...................................................14 3.4.2 Chargement................................................................................14 3.4.3 Rsultats.....................................................................................14 3.5 Elments coques.......................................................................................15 3.5.1 Gomtrie et degrs de libert ...................................................15 3.5.2 Chargement................................................................................15 3.5.3 Rsultats.....................................................................................16 5. PROCEDURES.......................................................................................................17 5.1 Environnement dexcution .....................................................................17 5.2 Mise en oeuvre des procdures ................................................................17 5.2.1 Pr-processeur...........................................................................17 a. Pr-processeur interactif .....................................................17 b. Sous diteur ........................................................................17 5.2.2 Rsolution.............................................................................................17 5.2.3 Avancement..........................................................................................18 a. Prparation du calcul...........................................................18 b. Arrt li une erreur de mise en donnes...........................18 c. Diagnostic de lerreur .........................................................18 d. Progression du calcul..........................................................19 5.2.5 Conversion postscript dun fichier graphique .......................................19 6. POST-PROCESSEUR GRAPHIQUE ....................................................................20 6.1 Lancement ...............................................................................................20 6.2 Commandes.............................................................................................20

ABAQUS 5.4

page 1

MANUEL UTILISATEUR

1. LES PRINCIPALES PROCEDURES ET LES FICHIERS

ABAQUS PRE

CAO

Editeur

Fichier de donnes XXXinp

ABAQUS

Fichier de procdure XXXlog

Fichier liste de rsultats XXXdat

Fichier binaire de reprise XXXres

Fichier binaire de rsultats XXXfil

Fichier de statistique XXXsta

ABAQUS POST Fichier graphique XXXmpl Fichier de comandes abaqus.jnl

LEGENDE PROCEDURE ABAQUS

ABAQUS PLOT Fichier postcript XXXps

Fichier dentre/sortie

Notes : Les noms de fichiers utiliss dans ABAQUS (XXXinp) sont constitus : - d un nom gnrique, fix par lutilisateur, correspondant par convention la rfrence du problme trait : XXX - dune extension, correspondant la lutilisation du fichier : eg. inp, dat, Chaque procdure de calcul (abaqus, abaqus plot, ) est indpendante.

ABAQUS 5.4

page 2

MANUEL UTILISATEUR

2. MISE EN DONNEES

2.1 Convention
ABAQUS utilise un fichier lecture ASCII pour la mise en donnes (XXXinp). Les diffrentes fonctions de calcul sont dfinies dans un ordre quasi-quelconque (malgr lorganisation de phases de mises en donnes proposes dans les paragraphes suivants) travers des mots cls. Le fichier de donnes comporte 3 types de lignes 1. les lignes de commentaires, prcdes de 2 toiles "**", exemple

** Dfinition des noeuds

2. les lignes de fonctions, comportant un mot cls, prcdes d'une toile "*", exemple

*NODE, ELSET=TOUT
2. les lignes d'arguments, conforment aux prescriptions de la fonction appeles, exemple

1, 1., 0.001, 2.
REMARQUES : Les lignes d'arguments situes entre deux mots cls se rapportent toujours la fonction prcdente. Il n'y a aucune ligne blanche. Elle serait considre comme un argument de la fonction prcdement appele.

2.2 Systme dunits

Comme dans beaucoup de logiciels, il nya pas de systme dunits fix. Cest lutilisateur de se dfinir son propre systme.

Masse M kg g Tonne

Longueur L m mm mm

Temps T s ms s

Force M L T-2 N N N

Pression M L -1 T -2 Pa MPa MPa

Energie M L 2 T-2 J mJ mJ

2.3 Phases de mises en donnes

Pour les applications courantes, on peut distinguer 4 phases dans la mise en donnes : Phase 1 : Gnration du maillage (noeuds, lments) Phase 2 : Proprits des matriaux (mcanique, physique, thermique) Phase 3 : Conditions aux limites Phase 4 : Etape de calcul (chargement, dfinition des rsultats)

ABAQUS 5.4

page 3

MANUEL UTILISATEUR

2.4 Mots cls ABAQUS

Pour les problmes mcaniques et thermiques linaires statiques, on sera couramment amen utiliser les mots-cls donns dans les tableaux suivants.

2.4.1 Phase 1 : Gnration du maillage (noeuds, lments)

Mots-cls NODE Fonction Cration des coordonnes (X, Y, Z) du noeud N Ligne 1 2 etc(2) Arguments/Type(1) *NODE, NSET= N, X, Y, Z I, R, R, R *NGEN, NSET= N1,N2, KN I, I, I *NSET, NSET=, (GENERATE) N1, , Ni I, I, , I *NFILL NS1, NS2, NL, NP I, I, I, I *ELEMENT, TYPE=(3) E, N1, N2, Ni I, I, I, *ELSET, ELSET=, (GENERATE) E1, , Ei I, I, , I *ELGEN E1,NE, KN, KE, NL, LN, LE I, I, I, I, I, I, I, I

NGEN

Gnration des noeuds de N1 N2 avec un pas KN 1 2 Assignation d'un nom un ensemble de noeuds N1, Ni. S'utilise seul ou associ NODE, NGEN, NFILL. Gnration des noeuds de l'ensemble de noeuds NS1 NS2 sur NL lignes avec un pas de NP sur les numros de noeuds. Cration d'un lments E, s'appuyant sur les noeuds N1, N2, (nbre de noeuds dpend du type d'lments, cf. 3) Assignation d'un nom un ensemble d'lments E1, Ei. S'utilise seul ou associ ELEMENT, ELGEN. Gnration partir de l'lment modle E1, de NE lments sur la 1re ligne, avec un pas de KN sur les numros de noeuds, KE sur les numros d'lements, sur NL lignes avec un pas de LN sur les noeuds et LE sur les lments 1 2 etc (2) 1 2 etc (2) 1 2 etc(2) 1 2 etc (2) 1 2 etc (2)

NSET

NFILL

ELEMENT

ELSET

ELGEN

Notes : (1) Type d'arguments : I : entier, R : rel, A : chane de caractres alphanumriques. (2) A rpter sur autant de lignes que ncessaire. (3) Quelques exemples de dsignation dlments sont donnes au paragraphe 3. Pour plus de dtails se reporter la documentation ABAQUS (User Manual I).

ABAQUS 5.4

page 4

MANUEL UTILISATEUR

2.4.2 Phase 2 : Proprits des matriaux (mcanique, physique)

Mots-cls SOLID SECTION BEAM SECTION Fonction Ligne Dfinition de l'paisseur T d'un 1 lment de type solide (2D). 2 Dfinition de la gomtrie de la 1 section d'un lment de type poutre. 2 Dfinition de la gomtrie de la 1 section d'un lment de type coque. 2 Dfinition d'un nom pour les 1 proprits de matriau dfinit ciaprs. Dfinition des proprits E, NU d'un 1 matriau lastique isotrope. 2 Coefficient de dilatation linaire 1 d'un matriau 2 Densit D d'un matriau. 1 2 Coefficient de conduction k d'un 1 matriau. 2 Capacit calorifique Cp d'un 1 matriau. 2 Arguments/Type(1) *SOLID SECTION, ELSET=, MATERIAL= T/R *BEAM SECTION, ELSET, MATERIAL=, SECTION=(RECT, BOX, CIR, PIPE) (3) / R *SHELL SECTION, ELSET=, MATERIAL= T, NI / R, I *MATERIAL, NAME=

SHELL SECTION MATERIAL

ELASTIC EXPANSION DENSITY CONDUCTIVITY SPECIFIC HEAT

*ELASTIC E, NU / R, R *EXPANSION /R *DENSITY D/R *CONDUCTIVITY k/R *SPECIFIC HEAT Cp / R

Notes : (1) : Type d'arguments : I : entier, R : rel, A : chane de caractres alphanumriques. (2) : A rpter sur autant de lignes que ncessaire. (3) : Si SECTION=CIR 2: rayon, SECTION=PIPE 2: rayon intrieur, rayon extrieur SECTION = RECT 2: hauteur, paisseur de la section SECTION = BOX hauteur,largeur de la section, paisseur voile z<0, y>0, z>0, y<0

CI R R Re Ri

PIPE

RECT d h l h

BOX c b l a

ABAQUS 5.4

page 5

MANUEL UTILISATEUR

2.4.3 Phase 3 : Conditions aux limites et conditions initiales

Mots-cls BOUNDARY Fonction Dfinition des conditions aux limites du noeud N ou de l'ensemble de noeuds SN en imposant la valeur V du dgr de libert DL1 DL2 (3) . Fixe une condition initiale de temprature TI au groupe de noeuds N. Ligne 1 2 etc (2) Arguments/Type(1) *BOUNDARY SN ou N, DL1, DL2, V I ou A, I, I, R *INITIAL CONDITIONS, TYPE=TEMPERATURE, TI / I, R

INITIAL CONDITIONS

Notes : (1) Type d'arguments : I : entier, R : rel, A : chane de caractres alphanumriques. (2) A rpter sur autant de lignes que ncessaire. (3) Dgr de libert : 1 = u, 2 = v, 3 = w, 4 = , 5 = , 6 = , 11 = T

2.4.4 Phase 4 : Etape de calcul a. Type de rsolution

Mots-cls STEP END STEP STATIC Fonction Dbut d'une tape de chargement. Fin d'une tape de chargement. Rsolution mcanique quasistatique avec une prcision PTOL sur les forces (dfaut 5.10-3), sur une dure D, avec un pas de temps initial DT, un pas de temps minimal DTm et un pas de temps maximal DTM. Calcul des NF frquences et modes propres avec la valeur de frquence propres maximales VM Rsolution thermique pure statique (STEADY STATE) ou transitoire (par dfaut) avec une augmentation de temprature maximale DELTMX par incrment. Pas de temps DT, priode de calcul T, minimum de pas de temps DTm, maximum de pas de temps DTM. Rsolution thermomcanique couple. Pas de - temps initial DT, priode de calcul T, minimum de pas de temps DTm, maximum de pas de temps DTM. Ligne 1 1 1 2 Arguments/Type(1) *STEP *END STEP *STATIC, (PTOL=) DT, D, DTm, DTM R,R,R,R

FREQUENCY

1 2 1 2

HEAT TRANSFER

FREQUENCY NF, VM I, R *HEAT TRANSFER, (DELTMX=, STEADY STATE) DT, T, DTm, DTM R, R, R, R *COUPLED TEMPERATURE DISPLACEMENT (DELTMX) DTI, T, DTm, DTM R, R, R, R

COUPLED TEMPERATURE DISPLACEMENT

1 2

ABAQUS 5.4

page 6

MANUEL UTILISATEUR

b. Dfinition du chargement
Mots-cls CLOAD Fonction Dfinition d'une charge concentre CC sur le noeud N ou l'ensemble de noeuds SN, suivant la composante du DDL (1, 2 ou 3). Multipli ventuellement par la courbe dvolution dfinie par AMPLITUDE. Si OP=NEW, annulation des CLOAD dfinis dans un STEP prcdent. Dfinition d'une charge rpartie CR sur l lment E ou l'ensemble dlments SE, avec le mots cls de la charge rpartie sur la face DI de llment. Multiplie ventuellement par la courbe dvolution dfinie par AMPLITUDE). Si OP=NEW, annulation des DLOAD dfinis dans un STEP prcdent. Application dun dplacement ou dune rotation impos U0 au noeud N ou l'ensemble de noeuds SN, du degr de libert DDL1 DDL2 (1, 6). Eventuellement couple avec la dfinition damplitude AMPLITUDE. Si OP=NEW, annulation des BOUNDARY dfinies dans un STEP prcdent. Application dune tmprature impose T au noeud N ou l'ensemble de noeuds SN. Ligne 1 2 etc (2) 1 2 etc (2) 1 2 etc (2) 1 *TEMPERATURE(,AMPLITUDE =, OP=) N/SN, T I/A, R *DFLUX (,AMPLITUDE=, OP=) E/SE, MF, FC I/A, A, R *AMPLITUDE, NAME=, DEFINITION= - TABULAR t1, y1, (4 couples/ligne) - ou PERIODIC n, w0 , t0 A1, B1, (2 couples/ligne) *BOUNDARY(,AMPLITUDE=, OP=) N/SN, DDL1, DDL2, U0 I/A, I, I, R *CLOAD(,AMPLITUDE=, OP=) E/SE, DI,CR I/A, A, R Arguments/Type(1) *CLOAD (,AMPLITUDE=, OP=) N/SN, DDL,CC I/A, I, R

DLOAD

BOUNDARY

TEMPERATURE

DFLUX

AMPLITUDE

2 etc (2) Application dun flux de chaleur FC sur l 1 lment E ou l'ensemble dlments SE, avec le mots cls du flux rparti MF, . Multipli 2 ventuellement par la courbe dvolution dfinie etc (2) par AMPLITUDE. Dfinition dune variation damplitude par 1 rapport au temps de manire tabule (TABULAR) ou priodique (PERIODIC : srie trigonomtrique) : 2 n y(t)= (Bi cos(w0 (t+t 0 ))+A i cos(w0 (t+t 0 ))) 2 i=1 3

Notes : (1) Type d'arguments : I : entier, R : rel, A : chane de caractres alphanumriques. (2) A rpter sur autant de lignes que ncessaire. (3) Voir paragraphe de dfinition des lments par des mots-cls de chargement.

ABAQUS 5.4

page 7

MANUEL UTILISATEUR

c. Ecriture des rsultats

Ecriture d'informations dans les fichiers de sortie XXXdat et XXXfil et XXXsta.

Mots-cls EL PRINT

NODE PRINT ENERGY PRINT EL FILE

NODE FILE ENERGY FILE

MONITOR

RESTART

Fonction Dfinition des informations lies aux lments. Ecriture dans le fichier de sortie XXXdat la frquence FREQ (si 0 alors 1 seule criture). Dfinition des informations lies aux noeuds. Ecriture dans le fichier de sortie XXXdat la frquence FREQ (si 0 alors 1 seule criture). Dfinition des informations lies aux nergies. Ecriture dans le fichier de sortie XXXdat la frquence FREQ (si 0 alors 1 seule criture). Dfinition des informations lies aux lments. Ecriture dans le fichier de sortie XXXfil la frquence FREQ (si 0 alors 1 seule criture). Dfinition des informations lies aux noeuds. Ecriture dans le fichier de sortie XXXfil la frquence FREQ (si 0 alors 1 seule criture). Dfinition des informations lies aux nergies. Ecriture dans le fichier de sortie XXXfil la frquence FREQ (si 0 alors 1 seule criture). Ecriture de la valeur dun degr de libert nodale DDL (dplacement : 1,2,3, rotation : 4,5,6, temprature : 11) la frquence F (FREQ) associ au noeud N dans le fichier sta. Ecriture (WRITE) et relecture (READ) de donnes et de rsultats danalyse la frquence F (FREQ) ou partir de lincrement Q (INC) du pas de chargement K (STEP) dans le fichier res. Permet lexploitation graphique de lensemble des rsultats (abaqus post)

Ligne 1 2 etc(2) 1 2 etc(2) 1 2 1 2 etc(2) 1 2 etc(2) 1 2 etc (2) 1

Arguments/Type(1) *EL PRINT (,ELSET=, FREQ=) (3) A, *NODE PRINT (,NSET=, FREQ=) (4) A, *ENERGY PRINT (,FREQ=) (5) *EL FILE (,ELSET=, FREQ=) (3) A, *NODE FILE, NSET=SN (4) A, *ENERGY FILE (5) A, *MONITOR, NODE=N, FREQ=F

DOF=DDL,

*RESTART (,WRITE/READ, STEP=K, INC=Q, FREQ=F)

Notes : (1) Type d'arguments : I : entier, R : rel, A : chane de caractres alphanumriques. (2) A rpter sur autant de lignes que ncessaire. (3) Voir Informations disponibles aux noeuds. (4) Voir rsultats pour chaque type dlment. (5) Voir Informations disponibles pour l'ensemble du problme.

ABAQUS 5.4

page 8

MANUEL UTILISATEUR

Informations disponibles aux noeuds, aux points dintgration ou globalement pour criture dans le fichier XXXdat ou XXXfil.

Variable

NODALE

Nom Vecteur dplacement et rotation Composante de dplacement i (1,2,3) Composante de rotation i (1,2,3) Vecteurs de raction et moments aux appuis Composante de raction aux appuis i (1,2,3) Composante de moment aux appuis i (1,2,3) Temprature Vecteurs de chargement par effort ou moment concentr Composante de chargement par effort i (1,2,3) Composante de chargement par moment i (1,2,3) Toutes les composantes du tenseur des contraintes Composantes ij du tenseur des contraintes i,j (1,2,3) Contraintes principales i (SP1 SP2 SP3) Pression quivalente (- ii/3) Contrainte quivalente de Von Mises ( 3 sijsij/2 ) (1) Toutes les composantes du tenseur des dformations Composantes ij du tenseur des dformations i,j (1,2,3) Dformations principales i (EP1 EP2 EP3) Flux de chaleur i (1, 2) Toutes les densits dnergie Densit dnergie de dformation lastique

ELEMENTAIRE

Argument U Ui URi RF RFi RMi NT CF CFi CMi S Sij SPi PRESS MISES E Eij EPi HFLi ENER SENER
ij - kk ij /3). Il contient la part de

GLOBAL Notes :

(1) : Le tenseur sij est le dviateur du tenseur des contraintes ( sij = cisailement du tenseur des contraintes.

ABAQUS 5.4

page 9

MANUEL UTILISATEUR

3. TYPE D'ELEMENTS
Quelques exemples dlments issus de la librairie dABAQUS.

3.1 Elements 1D solides 3.1.1 Gomtrie et degrs de libert

Dfinition des coordonnes : x, y, z Proprits dlments dfinis par le mots-cls SOLID SECTION (section et proprits physiques).

Type danalyse

Degr(s) de libert nodaux (u, v, w) : 1, 2, 3

Type

Rf.

Schma

MECANIQUE

2 noeuds, linaire

C1D2

3 noeuds, quadratique

C1D3

1
TRANSFERT CHALEUR DE (temprature) : 11 2 noeuds, linaire DC1D2

1
3 noeuds, quadratique DC1D3

1
TEMPERATUREDEPLACEMENT (u, v, w) : 1, 2, 3 (temprature) : 11 3 noeuds, quadratique DC1D3T 2 noeuds, linaire DC1D2T

3.1.2 Chargement

Cas de chargement Charge distribue Flux de chaleur volumique Flux de chaleur surfacique

Mots-cls BX, BY, BZ BF Sn

Description Force de volume /x, /y, /z Flux de chaleur par unit de volume Flux de chaleur par unit de surface sur le cot de llment.

Units FL -3 JL-3T-1 JL-2T-1

ABAQUS 5.4

page 10

MANUEL UTILISATEUR

3.2 Elements 2D solides 3.2.1 Gomtrie et degrs de libert

Dfinition des coordonnes : x, y ou r, Proprits dlments dfinis par le mots-cls SOLID SECTION (paisseur et proprits physiques).

Type danalyse

Degr(s) de libert nodaux (u, v) : 1, 2

Type

Rf.

Schma

MECANIQUE

3 noeuds, linaire

CPE3 /CPS3

3 (3) 1 (1) (2) 2

DEFORMATIONSP LANES / CONTRAINTES PLANES

6 noeuds, quadratique

CPE6 /CPS6

(3) 6 1

(2) 5

(1) 4
(3) 3

(u, v) : 1, 2

4 noeuds, linaire

CPE4 /CPS4

4 (4) 1

(2) (1) 2 (3) 7 3 6 2 (2) (1) 2 (3) 3 (2) (1) 2 (3) 3 (2) (1) 2 (3) 7 3 6 2 (2) (2)

8 noeuds, quadratique

CPE8 /CPS8

4 (4) 8 1

TRANSFERT CHALEUR

DE (temprature) : 11

3 noeuds, linaire

DC2D3

(1) 5 3

(3) 1

4 noeuds, linaire

DC2D4

4 (4) 1

TEMPERATUREDEPLACEMENT

(u, v) : 1, 2 (temprature) : 11

4 noeuds, linaire

CPE4T /CPS4T

4 (4) 1

DEFORMATIONSP LANES / CONTRAINTES PLANES

8 noeuds, quadratique, linaire

dplacement CPE8T temprature /CPS8T

4 (4) 8 1

(1) 5

ABAQUS 5.4

page 11

MANUEL UTILISATEUR

3.2.2 Chargement

Cas de chargement Charge distribue Force de volume Pression Flux de chaleur volumique Flux de chaleur surfacique

Mots-cls BX, BY CENT, GRAV Pn BF Sn

Description Force de volume /x, /y

Force centrifuge ( force de gravit Pression sur la face (n) FL -2 Flux de chaleur par unit de JL-3T-1 volume Flux de chaleur par unit de JL-2T-1 surface sur la face (n) de llment.

Units FL -3 2 ), FL -4,LT -2

3.2.3 Rsultats
Informations en sortie sur l'lment ou aux points d'intgration Contraintes : 11, 22, 12, 33 (non valide en contraintes planes). Dformations : 11, 22, 33, 12, 33 (non valide en dformations planes). Flux de chaleur : HFL1, HFL2.

ABAQUS 5.4

page 12

MANUEL UTILISATEUR

3.3 Elements 2D axisymtriques 3.3.1 Gomtrie et degrs de libert

Dfinition des coordonnes : x, y ou r, Proprits dlments dfinis par le mots-cls SOLID SECTION (paisseur et proprits physiques).

Type danalyse

Degr(s) de libertnodaux (u, v) : 1, 2

Type

Rf.

Schma

MECANIQUE

3 noeuds, linaire

CAX3

3 (3) 1 (1) (2) 2

6 noeuds, quadratique

CAX6

(3) 6 1

(2) 5

(1) 4
(3) 3

(u, v) : 1, 2

4 noeuds, linaire

CAX4

4 (4) 1

(2) (1) 2 (3) 7 3 6 2 (2) (1) 2 (3) 3 (2) (1) 2 (3) 3 (2) (1) 2 (3) 7 3 6 2 (2) (2)

8 noeuds, quadratique

CCAX8

4 (4) 8 1

TRANSFERT CHALEUR

DE (temprature) : 11

3 noeuds, linaire

DCAX3

(1) 5 3

(3) 1

4 noeuds, linaire

DCAX4

4 (4) 1

TEMPERATUREDEPLACEMENT

(u, v) : 1, 2 (temprature) : 11

4 noeuds, linaire

CAX4T

4 (4) 1

8 noeuds, quadratique, linaire

dplacement CAX8T temprature

4 (4) 8 1

(1) 5

ABAQUS 5.4

page 13

MANUEL UTILISATEUR

3.3.2 Chargement

Cas de chargement Charge distribue Force de volume Pression Flux de chaleur surfacique Flux de chaleur volumique

Most-cls BR, BZ CENT, GRAV Pn BF Sn

Description Force de volume /r, /z

Force centrifuge ( force de gravit Pression sur la face (n) FL -2 Flux de chaleur par unit de JL-3T-1 volume Flux de chaleur par unit de JL-2T-1 surface sur la face (n) de llment.

Units FL -3 2 ), FL -4,LT -2

3.2.3 Rsultats
Informations en sortie sur l'lment ou aux points d'intgration Contraintes : 11( rr ), 22 ( zz), 33, 12.

Dformations : 11( rr ), 22 ( zz), 33, 12. Flux de chaleur : HFL1, HFL2.

ABAQUS 5.4

page 14

MANUEL UTILISATEUR

3.4 Elements poutres 3.4.1 Gomtrie et degrs de libert

Dfinition des coordonnes : x, y, z Proprits dlments dfinis par le mots-cls BEAM SECTION (section et proprits physiques).

Type danalyse

Degr(s) de libert nodaux (u, v, w ) : 1, 2, 3 ( , , ) : 4, 5, 6

Type

Rf.

Schma

MECANIQUE

2 noeuds, linaire

B21/B31

PLAN / ESPACE

3 noeuds, quadratique

B22/B32

3.4.2 Chargement

Cas de chargement Charge distribue Force de volume

Mots-cls PX, PY, PZ CENT, GRAV

Description Units Force par unit de longueur FL -1 /x, /y, /z Force centrifuge ( 2 ), FL -4,LT-2 force de gravit

3.4.3 Rsultats
Informations en sortie sur l'lment ou aux points d'intgration. Contraintes : 11, 22 , 12. Dformations : 11, 22 , 12. Forces (SF), moments (SM) : SF1, SF2, SF3, SM1, SM2, SM3.

ABAQUS 5.4

page 15

MANUEL UTILISATEUR

3.5 Elments coques 3.5.1 Gomtrie et degrs de libert

Dfinition des coordonnes : x, y, z Proprits dlments dfinis par le mots-cls SHELL SECTION (section et proprits physiques).

Type danalyse

Degr(s) de libert nodaux (u, v, w ) : 1, 2, 3 ( , , ) : 4, 5, 6

Type

Rf.

Schma

MECANIQUE

3 noeuds, linaire

S3R

3 (3) 1 (1) (3) 3 (2) (1) 2 (2) 2

4 noeuds, linaire

S4R

4 (4) 1

8 noeuds, quadratique

S8R

TRANSFERT CHALEUR

DE (temprature ) : 11, 4 noeuds, quadratique 12, 17

DS4

7 3 4 (3) (4) (2) 6 8 2 1 (1) 5 (3) 3 4 (4) 1 (1) 2 (2)

TEMPERATUREDEPLACEMENT COUPLEE

(u, v, w ) : 1, 2, 3 ( , , ) : 4, 5, 6 (temprature ) : 11, 13, 17

8 noeuds, quadratique

S8RT

7 3 4 (3) (4) (2) 6 8 2 1 (1) 5

3.5.2 Chargement

Cas de chargement Charge distribue Force de volume Pression Flux de chaleur volumique Flux de chaleur surfacique

Mots-cls BX, BY, BZ CENT, GRAV P BF Sn

Description Force par unit de volume /x, /y, /z

Units FL -3 Force centrifuge ( 2 ), force de gravit FL -4,LT-2 Pression positive suivant la normale FL -2 llment. Flux de chaleur par unit de volume JL-3T-1 Flux de chaleur par unit de surface sur la JL-2T-1 face suprieure (1), infrieure (2) de llment .

ABAQUS 5.4

page 16

MANUEL UTILISATEUR

3.5.3 Rsultats
Informations en sortie sur l'lment ou aux points d'intgration. Contraintes : 11, 22 , 12. Dformations : 11, 22 , 12. Forces (SF), moments (SM) : SF1, SF2, SF3, SF4, SF5, SM1, SM2, SM3. h/2 + (SF1, SF2, SF3, SF4, SF5)= ( 11, 22, 12, 13, 23 ) dz , -h/2 h/2 + (SM1, SM2, SM3) = z ( 11, 22, 12) dz . -h/2

ABAQUS 5.4

page 17

MANUEL UTILISATEUR

5. PROCEDURES

5.1 Environnement dexcution

Les logiciels ABAQUS sont installs sur des machines serveurs europe (IBM/RS6000) et italie (SUN/Ultra Sparc 1), accessibles par le rseau EMAC. Pour xcuter une procdure, il faut se connecter sur une des machines, par exemple

$ xon italie xterm -ls

Aprs stre repositionn dans le rpertoire de travail (eg. cd tp/abaqus ), on peut xcuter une des procdures (excution, analyse de rsultats).

5.2 Mise en oeuvre des procdures

Toutes les procdures font appel au nom du problme (job) correspondant la racine des noms de fichiers (XXXinp, XXXdat, ) sans extension. En gnral, ce nom est fix par lutilisateur au moment de la cration du fichier de donnes inp.

5.2.1 Pr-processeur a. Pr-processeur interactif

Actuellement, seulement en local sur opra (SGi/R4000).

$ abaqus pr job = XXX b. Sous diteur

Ouvrir sour diteur (vi, lemacs) le fichier XXX. inp ou le crer. Il nest pas indispensable de se connecter au serveur de licence (italie ou europe) pour faire ldition de fichier. La machine europe ne supporte pas lemacs !

$ lemacs XXX 5.2.2 Rsolution

Le calcul est lanc en entrant la commande

$ abaqus Identifier : XXX

ou

$ abaqus job=XXX User routine file : (return)

Si il existe dj des fichiers de sortie (dans le mme rpertoire), il faut confirmer que l'on veut craser ces fichiers avec la nouvelle excution.

Old job files exist. Overwrite? (y/n):

ABAQUS 5.4

page 18

MANUEL UTILISATEUR

5.2.3 Avancement a. Prparation du calcul

Lors du lancement de lxcution, ABAQUS relie le fichier de donnes (XXXinp) pour crer un programme excutable (transparent pour lutilisateur). Cette phase de prparation du calcul correspond lxcution dun programme abaqus/bin/pre.x.

PID TTY 1263 pts/0 1267 pts/0

TIME CMD 0:00 pout.com 0:00 pre.x

Les rsultats de cette prparation dxcution sont cits dans le fichier XXXdat. En particulier, on trouve dans ce fichier la description des erreurs ventuelles de mise en donnes.

b. Arrt li une erreur de mise en donnes

En cas derreur dtecter lxecution de pre.x, le programme est arrt. La rsolution nest pas engage. Dans le fichier XXXlog, vous constaterez larrt en phase de prparation de calcul.

ABAQUS JOB pout vendredi, 29 fvrier 1999, 10:45:39 WET BEGIN USER INPUT PROCESSING Run /opt/abaqus/bin/pre.x ABAQUS Error: Error during user input processing. vendredi, 29 fvrier 1999, 10:45:41 WET ABAQUS JOB pout COMPLETED
En fin du fichier XXXdat ,vous trouvez linformation

THE PROGRAM HAS DISCOVERED 1 FATAL ERRORS ** EXECUTION IS TERMINATED ** END OF USER INPUT PROCESSING JOB TIME SUMMARY USER TIME = 0.11000 SYSTEM TIME = 6.00000E-02 TOTAL TIME = 0.17000 c. Diagnostic de lerreur
Pour dcouvrir lorigine de lerreur, il suffit de relire le fichier XXXdat et de rechercher le message fatal !

***ERROR: UNKNOWN KEYWORD CARD IMAGE: *TOTO ***NOTE: DUE TO AN INPUT ERROR THE ANALYSIS PREPROCESSOR HAS BEEN UNABLE TO INTERPRET SOME DATA.SUBSEQUENT ERRORS MAY BE CAUSED BY THIS OMISSION *NODE
Revenez au fichier de donnes XXXinp faite la correction qui simpose et relanc le calcul. Patience, cette tape de correction derreur et de test est dsagrables. Mais cest quasi-obligatoire quand on met au point un modle de calcul, surtout comme cest le cas ici, si vous faites directement la mise en donnes travers un diteur !

ABAQUS 5.4

page 19

MANUEL UTILISATEUR

d. Progression du calcul
La progression du calcul (prtraitement + rsolution) peut se consulter soit : - avec les commandes UNIX classiques : + ps (tat des process) + ls (tat des fichiers crs) - avec les fichiers (avec les commandes tail, more, vi ou emacs !) : + XXX dat + XXX log + XXX sta

5.2.5 Post-processeur
Le lancement du post-processeur graphique est ralis avec la commande

$ abaqus post job=XXX

Si aucun nom de fichier nest prcis, il est possible relire le fichier concern une fois la fentre X ouverte avec la commande

$ restart, file=XXX 5.2.4 Conversion postscript dun fichier graphique

1) Les sorties graphiques dabaqus post obtenues par la commande $ hardcopy sont en binaire (fichier XXXmpl). Avant de lancer limpression sur limprimante postcript (N&B et couleur), il suffit de transformer le fichier au format postcript par la procdure

$ abaqus plot job=XXX device=cps

Une suite de questions permettent de prciser les paramtres dimpression (orientation, couleur/N&B, nombre de dessin par page, ordre des dessins). Le fichier de sortie est XXXps. 2) Avant impression, il est possible de vrifier le contenu dun fichier postcript avec lutilitaire (SUN) ghostview

$ ghostview XXX.ps
3) Pour imprimer, ventuellement se connecter sur la machine serveur dimprimante, entrer

$ lpr -s XXX.ps
loption -s vite la recopie du fichier en crant un lien direct entre le buffer et le fichier source imprimer, ce qui rduit le temps de prparation dimpression.

ABAQUS 5.4

page 20

MANUEL UTILISATEUR

6. POST-PROCESSEUR GRAPHIQUE

6.1 Lancement
L'analyse par post-processeur est lance en entrant la commande

$ abaqus post job= nom

Aprs relecture du fichier XXXres ou XXXfil, ouverture d'un cran graphique.

6.2 Commandes
En fonction des informations spcifiques sauvegardes dans le fichier XXXfil, ou pour toutes les informations crites dans XXXres, on peut visualiser les rsultats en utilisant les commandes donnes aux tableaux suivants. Les paramtres sont spars de la commande par une virgule ,, par exemple

$ contour, v=s11
A lexcution dabaqus post un fichier qui contient la liste de toutes les commandes entres est crit. abaqusjnl, qui XXXfil. On peut ainsi se constituer un fichier pour faire des procdures de dpouillement automatique. Tableau 6.1 : Commandes graphiques pour abaqus post. Commande bc display contour Fonction Affichage des conditions aux limites Rsultats dans les lments ou aux noeuds Paramtres (1) variable= (all) Effet Affichage dune condition au limite particulire (toutes) contrainte 11 dplacement u1 tat dforme tat initial affiche la liste des commandes et permet daccder leur fonction. Trac de lvolution dune ou plusieurs variables nodales ou lmentaires en fonction du temps. La liste des variables est donne au tableau 6.2. 15 variables maximum par graphique. Trac dune variable nodale ou lmentaire en fonction dune autre variable nodale ou lmentaire. La liste des variables est donne au tableau 6.2. 15 courbes maximum par graphique.

draw end help history

v=s11 v=u1 Dessin du contour et du displaced maillage undisplaced [dfaut] sortie aide en ligne Trac de courbes en fonction titre x, titre y du temps. var1,ch., lgende, n noeud ou lement var2,ch., lgende,

variable

Trac de courbes.

input

titre x, titre y var1/x,ch., lgende, n noeud ou lement var1/x,ch., lgende, n noeud ou lement var2/x,ch., lgende, n noeud ou lement var2/x,ch., lgende, n noeud ou lement Relecture dun fichier de file=YYY relie et excute en automatique commandes YYY pour le les commandes entres dans post YYY

ABAQUS 5.4

page 21

MANUEL UTILISATEUR

Tableau 6.1 : Commandes graphiques pour abaqus post (suite). Commande load display Fonction Affichage du chragement Informations sur les valeurs particulires d'un lment. Information sur les valeurs particulires d'un noeud. Relecture dun fichier, ventuellement le mme, en fixant le pas et lincrment. Fonction gnrique, permettant de fixer des paramtres et dautoriser des actions. Laction sannule en rptant la commande. (bascule). Paramtres (1) variable= (all) Effet Affichage dun donn (tous ,0)

chargement

report element report node restart

file=XXX,(inc=, step=)

set

el number n number fill hardcopy (2) bc display load display c legend size= c magnification= c eigenmode= v=u v=rf

vector plot

Vecteur de quantits nodales

n des lments n de noeuds remplissage contour copie d'cran visualisation des conditions aux limites visualisation des chargements fixe la taille des carac tres de lgende (dfaut = 0.4) fixe lamplification des dplacements slectionne le mode propre afficher dplacement ractions dappuis

Notes : (1) : Entre parenthse, paramtres optionnels. (2) : Les fichiers de sortie graphiques crs avec la commande "set, hardcopy", sont dans un format binaire (XXXmpl) qu'il faut retraiter pour obtenir un fichier postcript (XXXps). Voir chapitre 5.

Tableau 6.2 : Exemple des variables aux noeuds et aux lments pour abaqus post (suite). Variables Grandeur mcanique Composante de dplacement i (1,2,3) Composante de rotation i (1,2,3) Composante de raction aux appuis i (1,2,3) Composante de moment aux appuis i (1,2,3) Temprature Composante de chargement par effort i (1,2,3) Composante de chargement par moment i (1,2,3) Composantes ij du tenseur des contraintes i,j (1,2,3) Contraintes principales i (SP1 SP2 SP3) Pression quivalente (- ii/3) Contrainte quivalente de Von Mises ( 3 sijsij/2 Composantes ij du tenseur des dformations i,j (1,2,3) Dformations principales i (EP1 EP2 EP3) Flux de chaleur i (1, 2) Toutes les densits dnergie Densit dnergie de dformation lastique Code ui uri rfi rmi nt11 cfi cmi sij spi press mises eij epi HFLi ener sener

NODALE

ELEMENTAIRE

GLOBAL

ABAQUS 5.4

page 22

MANUEL UTILISATEUR

6. INDEX

Elments