MODELISATION
Pour cette étude, la modélisation de cette structure en 3 Dimensions nous a amené à choisir les
paramètres dimensionnels suivants :
OPTI DIME 3 ELEM CUB8 ;
La directive OPTION sert à préciser des options générales de calcul.
DIME : Mot clé dimension permet d’entrer la dimension de l’espace dans laquelle on veut travailler
Exemple : 1 ; 2 ; 3.
ELEM : mot clé permet de préciser l’élément à fabriquer ;
Exemple : POI1, SEG2, SEG3, TRI3, TRI6, QUA4, QUA8, RAC2, RAC3, CUB8, CU20, PRI6, PR15, PYR5, PY13, TET4,
TE10.
DENS1=0.2 ;
La directive DENSITE sert à définir, par défaut, la taille locale FLOT1 (type FLOTTANT) de la maille
s'appuyant sur les points à construire.
Alternativement, on peut utiliser la directive 'OPTI' 'DENS' FLOT1 ;
➢ Création des points du portique.
A=1.600 ;
B=1.600;
C=2.000;
D=2.000;
XA=A;
XB=B;
XC=C ;
XD=D ;
Densité dens1 ;
P1=0. 0. 0. ;
P2=0. 0. XA ;
P3=XC 0. XB ;
P4=XD 0. 0. ;
➢ Création des lignes.
Pour créer les lignes il faut lier les points entre eux grâce au mot clé droi(te).
L1=P1 DROI 1 P2 ;
L2=P2 DROI 1 P3 ;
L3=P3 DROI 1 P4 ;
L4=P4 DROI 1 P1 ;
➢ Création des surfaces.
SURF1=L1 TRAN (0. 0.20 0.)1 ;
SURF2=L2 TRAN (0. 0.20 0.)1 ;
SURF3=L3 TRAN (0. 0.20 0.)1 ;
SURF4=L4 TRAN (0. 0.20 0.)1 ;
➢ Création des volumes.
VOL1=SURF1 VOLU TRAN (0.20 0. 0.)1 ;
VOL2=SURF2 VOLU TRAN (0. 0. -0.20)1 ;
VOL3=SURF3 VOLU TRAN (-0.20 0. 0.)1 ;
VOL4=SURF4 VOLU TRAN (0. 0. 0.20)1;
ELIM VOL1;
ELIM VOL2;
ELIM VOL3;
ELIM VOL4 ;
PORTIQ1=VOL1 ET VOL2 ET VOL3 ET VOL4 ;
TRAC PORTIQ1 QUAL ;
ELIM PORTIQ1 ;
➢ Caractéristique du béton
****CARACTERISTIQUES BETON*********
OPTI EPSI LINEAIRE ;
at=1.e0 ;
ac=1.2e0;
bc=1.5e+03;
bt=10000.;
epsd0=1.e-04 ;
YG=31.E+9 ;
Mod1=mode PORTIQ1 mécanique élastique isotrope endommagement Mazars ;
mat1=MATE mod1 YOUN YG NU 0.2 RHO 2.3E+03 btra bt ktr0 epsd0 atra at acom ac bcom bc beta
1.06e0;
L'opérateur MATE (MATERIAU) crée un champ de propriétés matérielles et/ou géométriques. Pour
les éléments qui nécessitent des propriétés matérielles et géométriques, on peut soit les introduire
toutes à la fois par MATE, soit introduire les propriétés matérielles par MATE et les propriétés
géométriques par CARA, puis fusionner les deux champs ainsi obtenus par ET.
'YOUN' : module d'Young
'NU ' : coefficient de poisson
'RHO ' : masse volumique
'EPSD0' : Seuil en déformation pour la traction (1.D-04)
'AC’ : Paramètre pour la compression (3000)
'BC’ : Paramètre pour la compression (4.)
'AT’ : Paramètre pour la traction (20000)
'BT’ : Paramètre pour la traction (1.6)
TRAC MOD1 MAT1;
RIG1= RIGI MOD1 MAT1;
MAT1= LUMP MAT1;
MAS1= MASS MOD1 MAT1;
MAS1= LUMP MAS1;
➢ Les conditions aux limites
PZ0= PORTIQ1 COOR 3 POIN MINI ;
CL1=BLOQ PZ0 UX UY UZ ;
TRAC (PZ0 COUL BLEU ET (ARET PORTIQ1)) ;
Calcul
PASAPAS (MECANIQUE) :
Cette procédure permet d'effectuer un calcul non linéaire incrémental
La non linéarité peut provenir, soit du matériau (plasticité), soit des grands déplacements
soit des deux à la fois.
Les résultats sont calculés a des valeurs du paramètre d'évolution (pseudo temps ou temps
réel) définies par l'utilisateur.
Sous l'option MODE FREQ, la procédure résout l'équation dynamique sur la base modale
étendue selon l'approche dite spectrale ou fréquentielle. Implicitement un instant est
interprété comme une fréquence pour l'objet CHARGEMENT sans changement de
terminologie.
La procédure propose un balayage de fréquences par défaut lorsque la liste des
TEMPS_CALCULES n'est pas précisée.
➢ BLOCAGES_MECANIQUES : blocages mécaniques (type RIGIDITE).
➢ CARACTERISTIQUES : Champ de caractéristiques matérielles et éventuellement
géométriques si nécessaire (type MCHAML, sous-type CARACTERISTIQUES)
Ses composantes peuvent être de type :
1) FLOTTANT si la composante est constante sur toute la structure ;
2) MCHAML si la composante dépend uniquement des points de la structure
3) EVOLUTION si la composante varie en fonction d'un seul paramètre.
4) NUAGE si la composante est décrite par une courbe de type EVOLUTION
dépendant d'un seul paramètre.
➢ CHARGEMENT : définition du chargement en fonction du paramètre d'évolution (type
CHARGEME) (crée par l'opérateur CHAR).
Les chargements élémentaires sont obligatoirement nommés :
- la température T
- les concentrations imposées CIMP
- les déplacements imposés DIMP
- les températures imposées TIMP
- Les autres chargements (meca) MECA
- Flux (en consolidation) FLUX
- Les flux de diffusion (diffusion) QCO
- Les flux de chaleur (thermique) Q
- Les températures ext (convection) TECO
- Les températures ext (rayonnement)TERA
- Une déformation imposée DEFI
- Les variables extérieures ?
➢ TEMPS_CALCULES : définition des valeurs du paramètre d'évolution (ou du temps)
pour lesquelles on effectue le calcul (type LISTREEL) (crée par l'opérateur PROG).
Dans le cas où cette donnée est absente, le code passe en ajustement automatique
de pas et considère comme liste des temps à calculer la liste des temps sauves
donnée sous l'indice TEMPS_SAUVES.
➢ MODELE : Objet modèle (type MMODEL) englobant toute la structure excepte celle
correspondant à la RIGIDITE_CONSTANTE.
OPTI ACQU 10 ACQU temps. List ;
ACQU PROGT*LISTREEL 2048 ;
OPTI ACQU 10 ACQU [Link];
ACQU PROGA*LISTREEL 2048 ;
Evol1=evol manu temps (prog 0. 10.) depl (prog 0. 10.) ;
CHX=MANU CHPO PORTIQ1 UX 1. ;
FX=-1. *MAS1*CHX ;
Evol2=evol manu temps (prog 0. 10.) depl (prog 0. 10.);
Evolution temporelle des coordonnées sur base modale
OPTION Evol permet d'obtenir l'évolution temporelle des coefficients de projection d'un
signal instationnaire (par exemple un résultat PASAPAS, DYNAMIC ou EXEC) sur les vecteurs
d'une base modale donnée.
L'opérateur PROG fabrique un objet LREEL1 de type LISTREEL à partir d'un nombre arbitraire
d'objets de type ENTIER ou FLOTTANT.
L'opérateur MANU construit, a partir d'une liste de points ou des points d'un autre maillage,
un maillage forme par ces points, constitue d'éléments dont le type est précise et de la
couleur requise.
F1=PRESSION MASSIF MOD1 (90000) SURF2;
char1=chargement FX EVOL1 MECA;
char2=chargement F1 EVOL2 MECA;
CHAR_TOT=CHAR1 ET CHAR2;
PAS1=0.1;
TAB1=TABLE;
[Link]=MOD1 ;
[Link]=MAT1 ;
[Link]=VRAI ;
[Link]=CHAR_TOT ;
TAB1.TEMPS_CALCULES=PROG 0. PAS PAS1 NPAS 10 ;
➢ La sous-directive PAS permet d'engendrer des nombres régulièrement espaces. Par
ailleurs des options correspondantes a certaines fonctions sont disponibles.
Si le PAS ne divise pas exactement l'intervalle, est retenue la valeur la plus proche le
permettant.
Si le PAS est incohérent, le résultat contiendra la valeur initiale et la valeur finale.
TAB1.BLOCAGES_MECANIQUES=CL1 ;
PASAPAS TAB1 ;
VISUALISATION
Point1=PORTIQ1 POINT PROCH (0. 0. 1.6) ;
DEBP EVDEP Point1*POINT TAB1*TABLE MOT1*MOT ;
N1=DIME [Link];
PROGT= PROG ;
PROGD= PROG;
REPETER BOU1 N1;
IT=&BOU1 - 1;
TPS1=TAB1. TEMPS. IT;
Droit1=EXTR TAB1. DEPLACEMENTS. IT Point1 MOT1 ;
PROGT=INSE PROGT &BOU1 TPS1 ;
PROGD=INSE PROGD &BOU1 Droit1 ;
FIN BOU1 ;
cordo=chaine 'Déplacement ' MOT1 ' (m)';
*on va suivre l’évolution du déplacement avec l’opérateur EVOL
EVOL2=EVOL MANU 'Temps (S)' PROGT CORDO PROGD;
FINP EVOL2;
EV1x=EVDEP Point1 TAB1 UX;
Opti sauv '[Link]';
sauv;
DESS EV1x;
➢ Directive DESSIN : Cette directive permet de tracer une EVOLUTION. Cette évolution
est une éventuelle concaténation de plusieurs sous-évolutions EVOLi.
@cartoon tab1 PORTIQ1 (0. 1.E8 0.1);
La procédure @CARTOON effectue une animation des déformées successives de
GEO1(PORTIQ1) obtenues à partir des déplacements contenus dans [Link]. Dans le
cas ou BLO1 est donne, les vecteurs correspondants aux réactions sont traces. En 3D cette
animation est visualisée suivant le point d'observation ŒIL(0.1 1.E8 0.1). On peut préciser
l'amplitude maximale des déformées avec le coefficient AMPL. Par défaut il est déterminé
automatiquement. Il est possible de désactiver l'option OSCIL qui est prise par défaut, en
utilisant l'option NOSCIL. On obtient alors une simple animation.
ev1=@global tab1 (evol1 ET EVOL2) CL1 FZ;
dess ev1;
ev2=@global tab1 (evol1 ET EVOL2) CL1 FY;
dess ev2;
************SIGMA_xx*************
Ni = dime [Link];
DEP1 = defo [Link]. (0) PORTIQ1 1. ;
SMXX = exco [Link]. (0) 'SMXX' ;
Trac SMXX MOD1 DEP1 (PORTIQ1) ;
Ni = dime [Link];
REPE TER BOUCLE (Ni-1) ;
DEP1=[Link]. (&BOUCLE) ;
DEF1 = defo DEP1 PORTIQ1 1 ;
SMXX = exco [Link]. (&BOUCLE)'SMXX' ;
Trac SMXX MOD1 DEF1 (PORTIQ1);
FIN BOUCLE ;
RESULTAT
COMPARAISON AVEC SAP2000
CAS DE DIAGONAL
Visualisation
Calcul pas à pas
NB : Pour le calcul pas a pas nous avons juste les début et fin car le calcul est beaucoup
Figure 1 Droite du déplacement en fonction du temps
Figure 3 Effort en fonction du déplacement suivant Z
Figure 2 effort en fonction du déplacement suivant z
Figure 4 effort en fonction du déplacement suivant z
CAS DE PORTIQUE AVEC REMPLISSAGE
Les procédures sont les mêmes que le précèdent
Figure 5 déplacement en fonction du temps suivant x
Figure 6 effort en fonction du placement suivant Z
Figure 7 effort en fonction du placement suivant X
SIGMA Suivant les différents
Figure 8 SIGMA_ZZ
Figure 9 SIGMA_YY
Figure 10 SIGMA_XZ
