* fichier : b52c.dgibi
'OPTI' 'ECHO' 0 'DIME' 2 'ELEM' 'QUA8' 'MODE' 'AXIS' ;
GRAPH = 'N' ;
* ajout de option epsilon lineaire car sinon ne converge pas!
OPTION epsilon lineaire;
**
** D3
**
** D4 S1 D2
**
** P1 D1 P2
*
* Tube
X1 = 10. ;
X2 = 11. ;
P1 = X1 0. ;
P2 = X2 0. ;
D1 = 'DROI' 10 P1 P2 ;
S1 = 'TRAN' D1 5 (0. 5.) ;
D2 = 'COTE' 2 S1 ;
D3 = 'COTE' 3 S1 ;
D4 = 'COTE' 4 S1 ;
PP = D2 'POIN' 'FINAL' ;
** THERMIQUE
* Température imposée à l'extérieur et l'intérieur du cylindre
T2 = 'TABLE' ;
T4 = 'TABLE' ;
T2.1 = 1400. ;
T4.1 = 1400. ;
T2.2 = 1400. ;
T4.2 = 1405. ;
T2.3 = 1400. ;
T4.3 = 1410. ;
T2.4 = 1400. ;
T4.4 = 1415. ;
T2.5 = 1400. ;
T4.5 = 1420. ;
T2.6 = 1400. ;
T4.6 = 1425. ;
T2.7 = 1400. ;
T4.7 = 1430. ;
* Facteur multiplicatif des temps
MULT1 = 1.E2 ;
* modèle
MODS1 = 'MODE' S1 'THERMIQUE' 'ISOTROPE' ;
**
* matériau
MATS1 = 'MATE' MODS1 'K' 5. ;
* conductivité
CONDS1 = 'COND' MODS1 MATS1 ;
* Température imposée
* sur la surface externe
BLOQ2 = 'BLOQ' D2 'T' ;
* sur la droite interne
BLOQ4 = 'BLOQ' D4 'T' ;
CHTER = 'TABLE' ;
NT = 7 ;
'REPE' BOULOU NT ;
LEA = (T4.&BOULOU - T2.&BOULOU) / (X1-X2) ;
LEB = T2.&BOULOU - (LEA * X2) ;
CHTER.&BOULOU = 'NOMC' 'T' (LEA * ('COOR' 1 S1) + LEB) ;
'FIN' BOULOU ;
* MECANIQUE
MODME = 'MODE' S1 'MECANIQUE' 'ELASTIQUE' 'FLUAGE' 'CERAMIQUE' ;
* On suppose les coefficient de fluage constant pour l'instant
QFLUGU = -5.3/8.61735E-5 ;
NTFLUGU = 1. ;
NFLUGU = 3.9 ;
AFLUGU = ('EXP' 28.8)/(1.E6**3.9);
* On met le terme dépendant de la température dans une EVOL
YOYO = 1.E9 ;
YOYO = 5.E11 ;
YOYO = 1.E10 ;
NUNU = 0.3 ;
DILA = 9.E-26 ;
RORO = 9000. ;
SIMAX = 1.E6 ;
QFLUGU = -5.3/8.61735E-5 ;
NTFLUGU = 1. ;
NFLUGU = 3.9 ;
AFLUGU = ('EXP' 28.8)/(1.E6**3.9);
***********************************
CHTER3 = 0. ;
'REPE' BOUOU 100 ;
CHTER3 = CHTER3 + 20.;
AVFLUK = AFLUGU * ('EXP' (QFLUGU * (CHTER3**-1.))) ;
'SI' ( CHTER3 '<' 1400.) ;
AVFLUK = 1.E-40 ;
'FINS' ;
'SI' ('EGA' &BOUOU 1) ;
PROT = 'PROG' CHTER3 ;
PROAV = 'PROG' AVFLUK ;
'SINO' ;
PROT = 'INSE' PROT &BOUOU CHTER3 ;
PROAV = 'INSE' PROAV &BOUOU AVFLUK ;
'FINS' ;
'FIN' BOUOU ;
AVFLUKE = 'EVOL' 'MANU' 'T' PROT 'AF1' PROAV;
DILA = 9.E-06 ;
MATME = 'MATE' MODME 'YOUN' YOYO 'NU' NUNU
'RHO' RORO 'ALPH' DILA
'SMAX' SIMAX
'AF1' AVFLUKE 'AF2' NFLUGU 'AF3' NTFLUGU 'TTRA' 0.
'ENDG' 100.01;
MAT2 = 'TAIL' MODME ;
* BLOCAGE DE D1
BLOCO = 'BLOQ' 'UZ' D1 ;
CHAMPI = 'DEPI' BLOCO 0. ;
* BLOCAGE DES POINTS DE D3
BLOCO3 = 'BLOQ' 'UZ' D3 ;
* Table des temps
TATEMPS = 'TABLE' ;
TATEMPS.0 = 0. ;
TATEMPS.1 = 1.*MULT1 ;
TATEMPS.2 = 2.*MULT1 ;
TATEMPS.3 = 3.*MULT1 ;
TATEMPS.4 = 4.*MULT1 ;
TATEMPS.5 = 5.*MULT1 ;
TATEMPS.6 = 6.*MULT1 ;
TATEMPS.7 = 7.*MULT1 ;
TATEMPS.8 = 8.*MULT1 ;
TATEMPS.9 = 9.*MULT1 ;
TATEMPS.10 = 10.*MULT1 ;
TATEMPS.11 = 11.*MULT1 ;
TATEMPS.12 = 12.*MULT1 ;
TATEMPS.13 = 13.*MULT1 ;
TATEMPS.14 = 14.*MULT1 ;
TATEMPS.15 = 15.*MULT1 ;
TATEMPS.16 = 16.*MULT1 ;
TATEMPS.17 = 17.*MULT1 ;
TATEMPS.18 = 18.*MULT1 ;
TATEMPS.19 = 19.*MULT1 ;
* Table d'évolution de la température
TABTER = 'TABLE' ;
TABTER.0 = 'MANU' 'CHPO' S1 1 'T' 0. ;
TABTER.1 = CHTER.1;
TABTER.2 = CHTER.2;
TABTER.3 = CHTER.3;
TABTER.4 = CHTER.4;
TABTER.5 = CHTER.5;
TABTER.6 = CHTER.6;
TABTER.7 = CHTER.7;
TABTER.8 = CHTER.7;
TABTER.9 = CHTER.7;
TABTER.10 = CHTER.7;
TABTER.11 = CHTER.7;
TABTER.12 = CHTER.7;
TABTER.13 = CHTER.7;
TABTER.14 = CHTER.7;
TABTER.15 = CHTER.7;
TABTER.16 = CHTER.7;
TABTER.17 = CHTER.7;
TABTER.18 = CHTER.7;
TABTER.19 = CHTER.7;
* Table du chargement mécanique
FORTOT = CHAMPI ;
CHARMEC = 'TABLE' ;
CHARMEC.0 = FORTOT 'ET' ('DEPI' BLOCO3 (DILA * 5.*('MAXI' TABTER.0)));
CHARMEC.1 = FORTOT 'ET' ('DEPI' BLOCO3 (DILA * 5.*('MAXI' TABTER.1)));
CHARMEC.2 = FORTOT 'ET' ('DEPI' BLOCO3 (DILA * 5.*('MAXI' TABTER.2)));
CHARMEC.3 = FORTOT 'ET' ('DEPI' BLOCO3 (DILA * 5.*('MAXI' TABTER.3)));
CHARMEC.4 = FORTOT 'ET' ('DEPI' BLOCO3 (DILA * 5.*('MAXI' TABTER.4)));
CHARMEC.5 = FORTOT 'ET' ('DEPI' BLOCO3 (DILA * 5.*('MAXI' TABTER.5)));
CHARMEC.6 = FORTOT 'ET' ('DEPI' BLOCO3 (DILA * 5.*('MAXI' TABTER.6)));
CHARMEC.7 = FORTOT 'ET' ('DEPI' BLOCO3 (DILA * 5.*('MAXI' TABTER.7)));
DILAT = (DILA * 5.*('MAXI' TABTER.7));
CHARMEC.8 = FORTOT 'ET' ('DEPI' BLOCO3 (DILAT+2.E-4) );
CHARMEC.9 = FORTOT 'ET' ('DEPI' BLOCO3 (DILAT+4.E-4));
CHARMEC.10 = FORTOT 'ET' ('DEPI' BLOCO3 (DILAT+6.E-4));
CHARMEC.11 = FORTOT 'ET' ('DEPI' BLOCO3 (DILAT+3.E-4));
CHARMEC.12 = FORTOT 'ET' ('DEPI' BLOCO3 DILAT);
CHARMEC.13 = FORTOT 'ET' ('DEPI' BLOCO3 (DILAT-3.E-4));
CHARMEC.14 = FORTOT 'ET' ('DEPI' BLOCO3 (DILAT-6.E-4));
CHARMEC.15 = FORTOT 'ET' ('DEPI' BLOCO3 (DILAT-9.E-4));
CHARMEC.16 = FORTOT 'ET' ('DEPI' BLOCO3 (DILAT-12.E-4));
CHARMEC.17 = FORTOT 'ET' ('DEPI' BLOCO3 (DILAT-15.E-4));
CHARMEC.18 = FORTOT 'ET' ('DEPI' BLOCO3 (DILAT-18.E-4));
CHARMEC.19 = FORTOT 'ET' ('DEPI' BLOCO3 (DILAT-21.E-4));
* Chargement temps force
CHARFOR = 'CHAR' 'MECA' TATEMPS CHARMEC ;
* Chargement temps température
CHARGT = 'CHAR' 'T' TATEMPS TABTER ;
TABA = 'TABLE' ;
TABA.'CHARGEMENT' = CHARGT 'ET' CHARFOR ;
LISTEMPS = 'PROG' 0. 1. 'PAS' 1. 19.;
LISTEMPS = LISTEMPS*MULT1 ;
TABA.'MODELE' = MODME ;
TABA.'CARACTERISTIQUES' = MATME 'ET' MAT2;
TABA.'TEMPS_CALCULES' = LISTEMPS ;
TABA.'TEMPS_SAUVES' = LISTEMPS ;
TABA.'BLOCAGES_MECANIQUES' = BLOCO 'ET' BLOCO3 ;
TABA.'MAXITERATION' = 500 ;
TABA.'PRECISION' = 4.E-4 ;
***********************************
PASAPAS TABA ;
I = 'DIME' TABA.'TEMPS' ;
* On essaie d'extraire les fissures
* on extrait les variables internes
FIFISS = TABA.'VARIABLES_INTERNES'.(I-1) ;
VECI = 'VECT' FIFISS MODME .01 ;
*
'REPE' BOULOU 3 ;
'SI' ('EGA' &BOULOU 1) ;
'SI' ('NEG' GRAPH 'N') ;
'TRAC' S1 VECI ;
'FINS' ;
MA = 'ELEM' 30 S1 ;
'SINO' ;
'SI' ('EGA' &BOULOU 2) ;
MA = 'ELEM' 21 S1 ;
'SINO' ;
MA = 'ELEM' 23 S1 ;
'FINS' ;
'FINS' ;
'REPE' BOILA I ;
* Contraintes
SI1 = 'REDU' TABA.'CONTRAINTES'.(&BOILA-1) MA ;
SILO2 = 'EXTR' SI1 'SMZZ' 1 1 2 ;
'MESS' 'SILO2 = ' SILO2 ;
* Déformations
DE1 = TABA.'DEPLACEMENTS'.(&BOILA-1) ;
EEP = 'REDU' ('EPSI' MODME DE1 ) MA ;
* On retranche la déformation thermique
ETH = 'EPTH' MODME MATME TABTER.(&BOILA-1) ;
EEPTH = 'REDU' ETH MA ;
EEVOU = EEP - EEPTH ;
EZDE2 = 'EXTR' EEVOU 'EPZZ' 1 1 2 ;
TIMER = TABA.'TEMPS'.(&BOILA-1) ;
'MESS' 'EZDE2 ='EZDE2 ;
'SI' ('EGA' &BOILA 1) ;
LOOX = 'PROG' EZDE2 ;
LOOY = 'PROG' SILO2 ;
LOOT = 'PROG' TIMER ;
'SINO' ;
LOOX = 'INSE' LOOX &BOILA EZDE2;
LOOY = 'INSE' LOOY &BOILA SILO2 ;
LOOT = 'INSE' LOOT &BOILA TIMER ;
'FINS' ;
'FIN' BOILA ;
'SI' ('NEG' GRAPH 'N') ;
TSIG = ('EVOL' 'MANU' 'TEMPS' LOOT 'SIGZZ' LOOY) ;
'DESS' TSIG ;
TEPS = ('EVOL' 'MANU' 'TEMPS' LOOT 'EPSZZ' LOOX) ;
'DESS' TEPS ;
TOLO = ('EVOL' 'MANU' 'EPSZZ' LOOX 'SIGZZ' LOOY) ;
'DESS' TOLO ;
'FINS' ;
'FIN' BOULOU ;
'MESS' '***************************************************' ;
'FIN' ;
© Cast3M 2003 - All rights reserved.
Disclaimer