Scribd
Importer
92 vues21 pages

Analyse des Contraintes et Déplacements

Le document décrit un modèle aux éléments finis d'une structure en béton armé avec de l'acier soumise à un chargement par déplacement imposé. Le modèle définit la géométrie, les matériaux, les conditions aux limites et effectue une analyse pas à pas pour calculer les contraintes et déplacements à différents pas de temps.

Transféré par

Duclair Hangouang
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats DOCX, PDF, TXT ou lisez en ligne sur Scribd
92 vues21 pages

Analyse des Contraintes et Déplacements

Le document décrit un modèle aux éléments finis d'une structure en béton armé avec de l'acier soumise à un chargement par déplacement imposé. Le modèle définit la géométrie, les matériaux, les conditions aux limites et effectue une analyse pas à pas pour calculer les contraintes et déplacements à différents pas de temps.

Transféré par

Duclair Hangouang
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats DOCX, PDF, TXT ou lisez en ligne sur Scribd

* FICHIER PROG/MEMOIRE/DEA/CAST3M

OPTI DIME 2 ELEM TRI3 MODE PLAN ;

*longueur de la section

a = 0.08 ;

*hauteur de la section

b = 0.08 ;

*rayon de l’acier

r = 0.01 ;

*enrobage

eb = 0.03 ;

*epaisseur

er = 0.01; ;

*****endommagement de MAZARS****

* seuil en déformation pour la traction

e1 = 1E-04 ;

* paramètre pour la compression

e2 = 1.4 ;

* paramètre pour la compression

e3 = 1900 ;

* paramètre pour la traction

e4 = 0.8 ;

* paramètre pour la traction

e5 = 17000 ;

*masse atomique de l’acier

M = 55.8454 ;

*valence du fer

n = 3 ;

Rédigé par HANGOUANG Duclair Page 1


*diametre initiale

d0 = 0.02 ;

*nombre de faraday

F = 96500 ;

*epsilon

esp = 1.20 ;

* humidité relative

HR = 0.62 ;

*temperature

Tp = 313 ;

C3 = -1.64 ;

*rapport E sur C

ec = 0.579 ;

*temps initiale

Ti = 1 ;

*temps d’évolution

t = 20 ;

*longueur de la barre

bb = 1 ;

*densité de l’acier

rhoa = 8000 ;

*densité de la rouille

rhor = 3600 ;

* correction pour le cisaillement

e6 = 1.06 ;

bb =32.13;

a1=-0.29;

Rédigé par HANGOUANG Duclair Page 2


g =-1.64;

enr =3 ;

l=0.5 ;

*les points

CEN1 = 0. 0. ;

P1 = r 0. ;

P2 = (r + eb) 0. ;

P3 = 0. r ;

P4 = r r ;

P5 = (r + eb) r ;

P6 = 0. (r + eb) ;

P7 = r (r + eb) ;

P8 = (r + eb) (r + eb) ;

Pa = (r+er) 0. ;

Pb = 0. (r + er) ;

*nombre d’éléments

N1 = 21 ;

N2 = 4 ;

N3 = 24 ;

N4 = 4 ;

N5 = 21 ;

N6 = 24 ;

N7 = 24 ;

N8 = 4 ;

N9 = 24 ;

Nr = 4 ;

Na = 3 ;

Rédigé par HANGOUANG Duclair Page 3


Nb = 3 ;

*les lignes

L1 = Pa DROIT N1 P2 ;

L2 = P1 DROIT N2 P4 ;

L3 = P4 DROIT N3 P5 ;

L4 = P2 DROIT N4 P5 ;

L5 = Pb DROIT N5 P6 ;

L6 = P4 DROIT N6 P7 ;

L7 = P5 DROIT N7 P8 ;

L8 = P6 DROIT N8 P7 ;

L9 = P7 DROIT N9 P8 ;

arc1 = C Nr P1 CEN1 P3 ;

La = P1 DROIT Na Pa ;

Lb = P3 DROIT Nb Pb ;

*contour

Contot = La ET L1 ET L4 ET L7 ET L9 ET

L8 ET L5 ET Lb ET arc1 ;

*surface

S1 = DALL (La ET L1) L4 L3 L2 COUL BLEU ;

S2 = DALL L3 L7 L9 L6 COUL BLEU ;

S3 = DALL arc1 (L2 ET L6) L8 (Lb ET L5) COUL BLEU ;

STOT = S1 ET S2 ET S3 ;

TRAC STOT;

****définition du modèle et matériaux****

*données matériaux et éléments*


MOD1 = MODE STOT MECANIQUE ELASTIQUE ISOTROPE

ENDOMMAGEMENT MAZARS ;

Rédigé par HANGOUANG Duclair Page 4


MAT1 = MATE MOD1 ‘YOUN’ 3E10 ‘NU’ 0.2 btra e5 ktr0
e1 atra e4 acom e2 bcom e3 beta e6 ;
*conditions aux limites

CL1 = BLOQ L1 UY ;

CL2 = BLOQ L5 UX ;

CL3 = BLOQ arc1 UX UY ;

CL = CL1 ET CL2 ET CL3 ;

*déplacement imposé est MM

*MM = 0.000002 ;

MM = (2.32*10 ** -0*(2.7182818285) ** (0.0312*HR-(4736 / Tp))*t) ;

*MM = ((HR*((1-(ec))**-1.64)*((t-1)**0.71)*t) / (eb)) ;

DEPL1 = DEPI CL3 MM ;

‘TITR’ ‘déplacement impose’ ;

‘TRAC’ (‘VECT’ ((DEPL1 ‘NOMC’ ‘UX’) ET (DEPL1 ‘NOMC’ ‘UY’)) 1.

‘UX’ ‘UY’ ‘ROUG’) (STOT) ;

‘FINS’ ;

*rigidité

RIGI1 = RIGI MOD1 MAT1 ;

RIGT = RIGI1 ET CL ;

*résolution

RES1 = RESOUDRE RIGT DEPL1 ;

*RESULTATS
DXY = EXCO RES1 UX UY ;
DXYABS = ABS DXY ;
DXYMAX = MAXI DXYABS ;
MESS 'DEPLACEMENT EN METRES DEPL1 :'DXYMAX ;
*****TRACE DE LA DEFORMEE ********
DEF0 = DEFORM STOT RES1 0. BLEU ;
DEF1 = DEFORM STOT RES1 1E3 ROUG ;
TITR ‘DEPLACEMENT DU BETON’ ;
TRAC (DEF0 ET DEF1) ;

Rédigé par HANGOUANG Duclair Page 5


******chargement par deplacement impose*******

****TEMPS 1*************

LIST1 = PROG 0. 1. 2. 3. 4. 5. 6. 7. 8. 9. 10.

11. 12. 13. 14. 15. 16. 17. 18. 19. 20. ;

LISF1 = PROG 0. 1. 2. 3. 4. 5. 6. 7. 8. 9. 10.

11. 12. 13. 14. 15. 16. 17. 18. 19. 20. ;

EVF1 = EVOL MANU 'temps' LIST1 'deplacement' LISF1 ;

CHAR1 = CHAR ‘DIMP’ DEPL1 EVF1 ;

SI(NEG GRAPH 'N') ;


DESS (EVF1*10) ;
FINSI ;

@EXCEL1 (EVF1*10) DEFOTEMPhr2 ;


FIN ;

************** TEMPS POUR LES DEPLACEMENTS***************

PDEPAR = 1. ;

PFINALR = 20. ;

PCHARR = PFINALR ;

*CALCUL PASAPAS

TAB1 = TABLE ;
TAB1. ’MOVA' = ‘D’ ;
TAB1. 'CHARGEMENT' = CHAR1 ;
TAB1. ‘MODELE’ = MOD1 ;
TAB1. 'CARACTERISTIQUES' = MAT1 ;
TAB1. 'BLOCAGES_MECANIQUES' = CL ;
TAB1. 'TEMPS_CALCULES' = LIST1 ;
TAB1. 'TEMPS_SAUVES' = LIST1 ;
TAB1. 'ENDOMMAGEMENT' = VRAI ;

TAB1.'NON_LOCAL'= VRAI ;

TAB1. 'GRANDS_DEPLACEMENTS' = VRAI ;


TAB1. 'MAXITERATION' = 5 ;

LC = 0.05 ;
TAB1. ’LCAR' = LC ;

Rédigé par HANGOUANG Duclair Page 6


TAB1. ’CONN' = ‘CONN’ MOD1 LC
‘NORMAL’ ‘INTERIEUR’ ;

TAB1. 'LVIAM’ = MOTS EPTI ;


PASAPAS TAB1 ;
*--- point de chargement
PCHAR = 10. ;
PFINAL = 1. ;
*--- contraintes
SIGMA1 = TAB1.'CONTRAINTES'.PCHAR ;
SIGMAV1 = VMIS SIGMA1 TAB1.'MODELE' ;
TITR 'CONTRAINTES DE VON MISES' ;
TRAC SIGMAV1 TAB1.'MODELE' ;
*--- déplacements
GEOTOT = EXTR TAB1.'MODELE' 'MAIL';
DEP1 = TAB1.'DEPLACEMENTS'.PCHAR;
DEP2 = TAB1.'DEPLACEMENTS'.PFINAL;
DEF0 = DEFO STOT DEP1 0. BLAN;
DEF1 = DEFO STOT DEP1 20. BLAN;
DEF2 = DEFO STOT DEP2 20. ROUG;
TITR 'DEPLACEMENTS EN FIN DE PRESSION';
TRAC (DEF0 ET DEF1) ;
TITR 'DEPLACEMENTS APRES PRESSION';
TRAC (DEF0 ET DEF2) ;
LISOV = PROG 9.0E+3 1.0E+4 1.5E+4 2.0E+4 2.5E+4 3.0E+4
3.5E+4 4.0E+4 4.5E+4 5.0E+4 5.5E+4 6.0E+4 6.5E+4 7.0E+4
7.5E+4 8.0E+4 8.5E+4 9.0E+4 9.5E+4 1.0E+5 1.5E+5 2.0E+5
2.5E+5 3.0E+5 3.5E+5 4.0E+5 4.5E+5 5.0E+5 5.5E+5 6.0E+5
6.5E+5 7.0E+5 7.5E+5 8.0E+5 8.5E+5 9.0E+5 9.5E+5 1.0E+6
1.5E+6 2.0E+6 2.5E+6 3.0E+6 3.5E+6 4.0E+6 4.5E+6 5.0E+6 
5.5E+6 6.0E+6 6.5E+6 7.0E+6 7.5E+6 8.0E+6 8.5E+6 9.0E+6 
9.5E+6 1.0E+7 1.5E+7 2.0E+7 ;
*List LISOV ;
****************ISOVALEURS DES CONTRAINTES****************
OPTI ISOV SURFACE ;
***********************CONTRAINTES***************************
****************CHARGEMENT AUX TEMPS PRECIS***************

SIGMA6= TAB1.'CONTRAINTES'. 6 ;

SIGMA7= TAB1.'CONTRAINTES'. 7 ;

SIGMA8= TAB1.'CONTRAINTES'. 8 ;

Rédigé par HANGOUANG Duclair Page 7


SIGMA9= TAB1.'CONTRAINTES'. 9 ;

****************CONTRAINTE DE VON MISES *******************

SIGMAV6 = VMIS SIGMA6 TAB1.'MODELE' ;

SVABS6 =ABS(sigmav6);

TITR 'contraintes de von mises AU TEMPS 6 ';

TRAC SIGMAV6 TAB1.'MODELE' contot;

TITR ' contraintes de von misesABS AU TEMPS 6 ';

TRAC SVABS6 TAB1.'MODELE' contot;

SIGMAV7 = VMIS SIGMA7 TAB1.'MODELE' ;

SVABS7 =ABS(sigmav7);

TITR ' contraintes de von mises AU TEMPS 7 ';

TRAC SIGMAV7 TAB1.'MODELE' contot;

TITR ' contraintes de von mises ABS AU TEMPS 7 ';

TRAC SVABS7 TAB1.'MODELE' contot ;

SIGMAV8 = VMIS SIGMA8 TAB1.'MODELE' ;

SVABS8 =ABS(sigmav8);

TITR ' contraintes de von mises AU TEMPS 8 ';

TRAC SIGMAV8 TAB1.'MODELE' contot;

TITR ' contraintes de von mises ABS AU TEMPS 8 ';

TRAC SVABS8 TAB1.'MODELE' contot;

SIGMAV9 = VMIS SIGMA9 TAB1.'MODELE' ;

SVABS9 =ABS(sigmav9);

TITR ' contraintes de von mises AU TEMPS 9 ';

TRAC SIGMAV9 TAB1.'MODELE' contot;

TITR ' contraintes de von misesABS AU TEMPS 9 ';

TRAC SVABS9 TAB1.'MODELE' contot;

*************CONTRAINTE DE VON MISES MAX ****************

Rédigé par HANGOUANG Duclair Page 8


SIGMA1= TAB1.'CONTRAINTES'.PCHAR ;

SIGMAV1 = VMIS SIGMA1 TAB1.'MODELE' ;

SVABS =ABS(sigmav1);

sxmax = maxi svabs;

mess 'contrainte maximale de von mises =' sxmax;

TITR 'contraintes de von mises MAX ';

TRAC SIGMAV1 TAB1.'MODELE' contot;

TITR 'contraintes de von mises en ABS MAX ';

TRAC SVABS TAB1.'MODELE' contot;

TITR 'contraintes de von mises + LISOV MAX ';

TRAC LISOV SIGMAV1 TAB1.'MODELE' contot;

TITR 'contraintes de von mises en ABS + LISOV MAX ';

TRAC LISOV SVABS TAB1.'MODELE' contot;

***********************CONTRAINTE XX************************

S_XX6 = EXCO SMXX SIGMA6 ;

S_XXABS6 =ABS(s_xx6);

TITR'CONTRAINTES S_XX AU TEMPS 6 ';

TRAC S_XX6 TAB1.'MODELE' contot;

TITR'CONTRAINTES S_XXABS6 AU TEMPS 6 ';

TRAC S_XXABS6 TAB1.'MODELE' contot;

S_XX7 = EXCO SMXX SIGMA7 ;

S_XXABS7 =ABS(s_xx7);

TITR'CONTRAINTES S_XX AU TEMPS7 ';

TRAC S_XX7 TAB1.'MODELE' contot ;

TITR'CONTRAINTES S_XXABS AU TEMPS 7 ';

TRAC S_XXABS7 TAB1.'MODELE' contot ;

S_XX8 = EXCO SMXX SIGMA8 ;

Rédigé par HANGOUANG Duclair Page 9


S_XXABS8 =ABS(s_xx8);

TITR'CONTRAINTES S_XX AU TEMPS 8 ';

TRAC S_XX8 TAB1.'MODELE' contot ;

TITR'CONTRAINTES S_XXABS AU TEMPS 8 ';

TRAC S_XXABS8 TAB1.'MODELE' contot ;

S_XX9 = EXCO SMXX SIGMA9 ;

S_XXABS9 =ABS(s_xx9);

TITR'CONTRAINTES S_XX AU TEMPS 9 ';

TRAC S_XX9 TAB1.'MODELE' contot ;

TITR'CONTRAINTES S_XXABS AU TEMPS 9 ';

TRAC S_XXABS9 TAB1.'MODELE' contot ;

***********************CONTRAINTE YY************************

S_YY6 = EXCO SMYY SIGMA6 ;

S_YYABS6 =ABS(s_yy6);

TITR'CONTRAINTES S_YY AU TEMPS 6 ';

TRAC S_YY6 TAB1.'MODELE' contot ;

TITR'CONTRAINTES S_YYABS6 AU TEMPS 6 ';

TRAC S_YYABS6 TAB1.'MODELE' contot ;

S_YY7 = EXCO SMYY SIGMA7 ;

S_YYABS7 =ABS(s_yy7);

TITR'CONTRAINTES S_YY AU TEMPS7 ';

TRAC S_YY7 TAB1.'MODELE' contot ;

TITR'CONTRAINTES S_YYABS AU TEMPS 7 ';

TRAC S_YYABS7 TAB1.'MODELE' contot ;

S_YY8 = EXCO SMYY SIGMA8 ;

S_YYABS8 =ABS(s_yy8);

Rédigé par HANGOUANG Duclair Page 10


TITR'CONTRAINTES S_YY AU TEMPS 8 ';

TRAC S_YY8 TAB1.'MODELE' contot ;

TITR'CONTRAINTES S_YYABS AU TEMPS 8 ';

TRAC S_YYABS8 TAB1.'MODELE' contot ;

S_YY9 = EXCO SMYY SIGMA9 ;

S_YYABS9 =ABS(s_yy9);

TITR'CONTRAINTES S_YY AU TEMPS 9 ';

TRAC S_YY9 TAB1.'MODELE' contot ;

TITR'CONTRAINTES S_YYABS AU TEMPS 9 ';

TRAC S_YYABS9 TAB1.'MODELE' contot ;

**********COURBES CONTRAINTE DEPLACEMENT *************

ptt = prog 0 ;

psz = prog 0 ;

ic = 0 ;

repeter mabou 20 ;

ic = ic + 1 ;

dd = tab1 . deplacements . ic;

ptt = ptt et ( prog ( extr dd uy p1 ) ) ;

ptt = ABS(ptt) ;

ss = tab1.contraintes.ic ;

sz = extr ( chan chpo mod1 ss ) smyy p1 ;

psz = psz et ( prog sz ) ;

psz = ABS(psz) ;

fin mabou ;

evsz1 = evol vert manu deplacement ptt contrainte psz ;

dess evsz1 TITRE ‘σ y = f (u y ) au point P1’ ;

Rédigé par HANGOUANG Duclair Page 11


@EXCEL1 EVSZ1 CONTRDEFOhr1;

*EXTRACTION DE L’ENDOMMAGEMENT

dd1 = TAB1.VARIABLES_INTERNES. 1 ;

dd2 = TAB1.VARIABLES_INTERNES. 2 ;

dd3 = TAB1.VARIABLES_INTERNES. 3 ;

dd4 = TAB1.VARIABLES_INTERNES. 4 ;

dd5 = TAB1.VARIABLES_INTERNES. 5 ;

dd6 = TAB1.VARIABLES_INTERNES. 6 ;

dd7 = TAB1.VARIABLES_INTERNES. 7 ;

dd8 = TAB1.VARIABLES_INTERNES. 8 ;

dd9 = TAB1.VARIABLES_INTERNES. 9 ;

dd10 = TAB1.VARIABLES_INTERNES. 10 ;

dd11 = TAB1.VARIABLES_INTERNES. 11 ;

dd12 = TAB1.VARIABLES_INTERNES. 12 ;

dd13 = TAB1.VARIABLES_INTERNES. 13 ;

dd14 = TAB1.VARIABLES_INTERNES. 14 ;

dd15 = TAB1.VARIABLES_INTERNES. 15 ;

dd16 = TAB1.VARIABLES_INTERNES. 16 ;

dd17 = TAB1.VARIABLES_INTERNES. 17 ;

dd18 = TAB1.VARIABLES_INTERNES. 18 ;

dd19 = TAB1.VARIABLES_INTERNES. 19 ;

dd20 = TAB1.VARIABLES_INTERNES. 20 ;

Endo1 = EXCO 'D' DD1 ;

Endo2 = EXCO 'D' DD2 ;

Endo3 = EXCO 'D' DD3 ;

Endo4 = EXCO 'D' DD4 ;

Rédigé par HANGOUANG Duclair Page 12


Endo5 = EXCO 'D' DD5 ;

Endo6 = EXCO 'D' DD6 ;

Endo7 = EXCO 'D' DD7 ;

Endo8 = EXCO 'D' DD8 ;

Endo9 = EXCO 'D' DD9 ;

Endo10 = EXCO 'D' DD10 ;

Endo11 = EXCO 'D' DD11 ;

Endo12 = EXCO 'D' DD12 ;

Endo13 = EXCO 'D' DD13 ;

Endo14 = EXCO 'D' DD14 ;

Endo15 = EXCO 'D' DD15 ;

Endo16 = EXCO 'D' DD16 ;

Endo17 = EXCO 'D' DD17 ;

Endo18 = EXCO 'D' DD18 ;

Endo19 = EXCO 'D' DD19 ;

Endo20 = EXCO 'D' DD20 ;

TRAC ENDO1 TAB1.'MODELE' contot ;

TRAC ENDO2 TAB1.'MODELE' contot ;

TRAC ENDO3 TAB1.'MODELE' contot ;

TRAC ENDO4 TAB1.'MODELE' contot ;

TRAC ENDO5 TAB1.'MODELE' contot ;

TRAC ENDO6 TAB1.'MODELE' contot ;

TRAC ENDO7 TAB1.'MODELE' contot ;

TRAC ENDO8 TAB1.'MODELE' contot ;

TRAC ENDO9 TAB1.'MODELE' contot ;

TRAC ENDO10 TAB1.'MODELE' contot ;

TRAC ENDO11 TAB1.'MODELE' contot ;

Rédigé par HANGOUANG Duclair Page 13


TRAC ENDO12 TAB1.'MODELE' contot ;

TRAC ENDO13 TAB1.'MODELE' contot ;

TRAC ENDO14 TAB1.'MODELE' contot ;

TRAC ENDO15 TAB1.'MODELE' contot ;

TRAC ENDO16 TAB1.'MODELE' contot ;

TRAC ENDO17 TAB1.'MODELE' contot ;

TRAC ENDO18 TAB1.'MODELE' contot ;

TRAC ENDO19 TAB1.'MODELE' contot ;

TRAC ENDO20 TAB1.'MODELE' contot ;

PPP = 300;

D1 = ( extr dd1 'D ' 1 PPP 1 ) ;

D2 = ( extr dd2 'D ' 1 PPP 1 ) ;

D3 = ( extr dd3 'D ' 1 PPP 1 ) ;

D4 = ( extr dd4 'D ' 1 PPP 1 ) ;

D5 = ( extr dd5 'D ' 1 PPP 1 ) ;

D6 = ( extr dd6 'D ' 1 PPP 1 );

D7 = ( extr dd7 'D ' 1 PPP 1 ) ;

D8 = ( extr dd8 'D ' 1 PPP 1 ) ;

D9 = ( extr dd9 'D ' 1 PPP 1 ) ;

D10 = ( extr dd10 'D ' 1 PPP 1 ) ;

D11 = ( extr dd11 'D ' 1 PPP 1 ) ;

D12 = ( extr dd12 'D ' 1 PPP 1 ) ;

D13 = ( extr dd13 'D ' 1 PPP 1 ) ;

D14 = ( extr dd14 'D ' 1 PPP 1 ) ;

D15 = ( extr dd15 'D ' 1 PPP 1 ) ;

D16 = ( extr dd16 'D ' 1 PPP 1 ) ;

D17 = ( extr dd17 'D ' 1 PPP 1 ) ;

Rédigé par HANGOUANG Duclair Page 14


D18 = ( extr dd18 'D ' 1 PPP 1 ) ;

D19 = ( extr dd19 'D ' 1 PPP 1 ) ;

D20 = ( extr dd20 'D ' 1 PPP 1 ) ;

EF_D = PROG D1 D2 D3 D4 D5 D6 D7 D8 D9

D10 D11 D12 D13 D14 D15 D16 D17 D18 D19 D20 ;

*LIST EF_D ;

*EXTRACTION DES CONTRAITES

CONTP1 = TAB1.CONTRAINTES. 1 ;

CONTP2 = TAB1.CONTRAINTES. 2 ;

CONTP3 = TAB1.CONTRAINTES. 3 ;

CONTP4 = TAB1.CONTRAINTES. 4 ;

CONTP5 = TAB1.CONTRAINTES. 5 ;

CONTP6 = TAB1.CONTRAINTES. 6 ;

CONTP7 = TAB1.CONTRAINTES. 7 ;

CONTP8 = TAB1.CONTRAINTES. 8 ;

CONTP9 = TAB1.CONTRAINTES. 9 ;

CONTP10 = TAB1.CONTRAINTES. 10 ;

CONTP11 = TAB1.CONTRAINTES. 11 ;

CONTP12 = TAB1.CONTRAINTES. 12 ;

CONTP13 = TAB1.CONTRAINTES. 13 ;

CONTP14 = TAB1.CONTRAINTES. 14 ;

CONTP15 = TAB1.CONTRAINTES. 15 ;

CONTP16 = TAB1.CONTRAINTES. 16 ;

CONTP17 = TAB1.CONTRAINTES. 17 ;

CONTP18 = TAB1.CONTRAINTES. 18 ;

CONTP19 = TAB1.CONTRAINTES. 19 ;

CONTP20 = TAB1.CONTRAINTES. 20 ;

Rédigé par HANGOUANG Duclair Page 15


S1 = ( extr CONTP1 'SMYY' 1 PPP 1 ) ;

S2 = ( extr CONTP2 'SMYY' 1 PPP 1 ) ;

S3 = ( extr CONTP3 'SMYY' 1 PPP 1 ) ;

S4 = ( extr CONTP4 'SMYY' 1 PPP 1 ) ;

S5 = ( extr CONTP5 'SMYY' 1 PPP 1 ) ;

S6 = ( extr CONTP6 'SMYY' 1 PPP 1 ) ;

S7 = ( extr CONTP7 'SMYY' 1 PPP 1 ) ;

S8 = ( extr CONTP8 'SMYY' 1 PPP 1 ) ;

S9 = ( extr CONTP9 'SMYY' 1 PPP 1 ) ;

S10 = ( extr CONTP10 'SMYY' 1 PPP 1 ) ;

S11 = ( extr CONTP11 'SMYY' 1 PPP 1 ) ;

S12 = ( extr CONTP12 'SMYY' 1 PPP 1 ) ;

S13 = ( extr CONTP13 'SMYY' 1 PPP 1 ) ;

S14 = ( extr CONTP14 'SMYY' 1 PPP 1 ) ;

S15 = ( extr CONTP15 'SMYY' 1 PPP 1 ) ;

S16 = ( extr CONTP16 'SMYY' 1 PPP 1 ) ;

S17 = ( extr CONTP17 'SMYY' 1 PPP 1 ) ;

S18 = ( extr CONTP18 'SMYY' 1 PPP 1 ) ;

S19 = ( extr CONTP19 'SMYY' 1 PPP 1 ) ;

S20 = ( extr CONTP20 'SMYY' 1 PPP 1 ) ;

EF_S = PROG S1 S2 S3 S4 S5 S6 S7 S8 S9 S10

S11 S12 S13 S14 S15 S16 S17 S18 S19 S20;

EF_S = ABS(EF_S);

*EXTRACTION DE DEPLACEMENT

CONTP1 = TAB1.DEPLACEMENTS. 1 ;

CONTP2 = TAB1.DEPLACEMENTS. 2 ;

CONTP3 = TAB1.DEPLACEMENTS. 3 ;

Rédigé par HANGOUANG Duclair Page 16


CONTP4 = TAB1.DEPLACEMENTS. 4 ;

CONTP5 = TAB1.DEPLACEMENTS. 5 ;

CONTP6 = TAB1.DEPLACEMENTS. 6 ;

CONTP7 = TAB1.DEPLACEMENTS. 7 ;

CONTP8 = TAB1.DEPLACEMENTS. 8 ;

CONTP9 = TAB1.DEPLACEMENTS. 9 ;

CONTP10 = TAB1.DEPLACEMENTS. 10 ;

CONTP11 = TAB1.DEPLACEMENTS. 11 ;

CONTP12 = TAB1.DEPLACEMENTS. 12 ;

CONTP13 = TAB1.DEPLACEMENTS. 13 ;

CONTP14 = TAB1.DEPLACEMENTS. 14 ;

CONTP15 = TAB1.DEPLACEMENTS. 15 ;

CONTP16 = TAB1.DEPLACEMENTS. 16 ;

CONTP17 = TAB1.DEPLACEMENTS. 17 ;

CONTP18 = TAB1.DEPLACEMENTS. 18 ;

CONTP19 = TAB1.DEPLACEMENTS. 19 ;

CONTP20 = TAB1.DEPLACEMENTS. 20 ;

S1 = ( extr CONTP1 'UY' P1 ) ;

S2 = ( extr CONTP2 'UY' P1 ) ;

S3 = ( extr CONTP3 'UY' P1 ) ;

S4 = ( extr CONTP4 'UY' P1 ) ;

S5 = ( extr CONTP5 'UY' P1 ) ;

S6 = ( extr CONTP6 'UY' P1 ) ;

S7 = ( extr CONTP7 'UY' P1 ) ;

S8 = ( extr CONTP8 'UY' P1 ) ;

S9 = ( extr CONTP9 'UY' P1 ) ;

S10 = ( extr CONTP10 'UY' P1 ) ;

Rédigé par HANGOUANG Duclair Page 17


S11 = ( extr CONTP11 'UY' P1 ) ;

S12 = ( extr CONTP12 'UY' P1 ) ;

S13 = ( extr CONTP13 'UY' P1 ) ;

S14 = ( extr CONTP14 'UY' P1 ) ;

S15 = ( extr CONTP15 'UY' P1 ) ;

S16 = ( extr CONTP16 'UY' P1 ) ;

S17 = ( extr CONTP17 'UY' P1 ) ;

S18 = ( extr CONTP18 'UY' P1 ) ;

S19 = ( extr CONTP19 'UY' P1 ) ;

S20 = ( extr CONTP20 'UY' P1 ) ;

EF_DE = PROG S1 S2 S3 S4 S5 S6 S7 S8 S9 S10

S11 S12 S13 S14 S15 S16 S17 S18 S19 S20;

S1 = ( extr CONTP1 'UY' P1 ) ;

S2 = ( extr CONTP2 'UY' P1 ) ;

S3 = ( extr CONTP3 'UY' P1 ) ;

S4 = ( extr CONTP4 'UY' P1 ) ;

S5 = ( extr CONTP5 'UY' P1 ) ;

S6 = ( extr CONTP6 'UY' P1 ) ;

S7 = ( extr CONTP7 'UY' P1 ) ;

S8 = ( extr CONTP8 'UY' P1 ) ;

S9 = ( extr CONTP9 'UY' P1 ) ;

S10 = ( extr CONTP10 'UY' P1 ) ;

S11 = ( extr CONTP11 'UY' P1 ) ;

S12 = ( extr CONTP12 'UY' P1 ) ;

S13 = ( extr CONTP13 'UY' P1 ) ;

S14 = ( extr CONTP14 'UY' P1 ) ;

S15 = ( extr CONTP15 'UY' P1 ) ;

Rédigé par HANGOUANG Duclair Page 18


S16 = ( extr CONTP16 'UY' P1 ) ;

S17 = ( extr CONTP17 'UY' P1 ) ;

S18 = ( extr CONTP18 'UY' P1 ) ;

S19 = ( extr CONTP19 'UY' P1 ) ;

S20 = ( extr CONTP20 'UY' P1 ) ;

EF_DER = PROG S1 S2 S3 S4 S5 S6 S7 S8 S9 S10

S11 S12 S13 S14 S15 S16 S17 S18 S19 S20;

evszR = evol vert manu deplacement EF_DER contrainte EF_S ;

dess evszR TITRE ‘σ y = f (u y ) au point P1’ ;

evszP2 = evol vert manu deplacement EF_DE contrainte EF_S ;

dess evszP2 TITRE ‘σ y = f (u y ) au point P1’ ;

ptt = prog 0 ;

psz = prog 0 ;

ic = 0 ;

repeter mabou 19 ;

ic = ic + 1 ;

dd = tab1 . deplacements . ic;

ptt = ptt et ( prog ( extr dd uy p3 ) ) ;

ptt = ABS(ptt) ;

LIST ptt ;

ss = tab1.contraintes.ic ;

sz = extr ( chan chpo mod1 ss ) smyy p3 ;

psz = psz et ( prog sz ) ;

psz = ABS(psz) ;

fin mabou ;

evsz1 = evol vert manu deplacement ptt contrainte psz ;

dess evsz1 TITRE ‘σ y = f (u y ) au point P3’ ;

Rédigé par HANGOUANG Duclair Page 19


evend = evol bleu manu deplacement ptt endommagement EF_D ;

dess evend TITRE ‘endommagement = f( deplacement)’ ;

@EXCEL1 EVEND ENDODEPLhr1 ;

evend1 = evol ROUG manu endommagement EF_D contrainte psz ;

dess evend1 TITRE ‘contrainte = f(endommagement )’ ;

@EXCEL1 EVEND1 CONENDOhr1 ;

TE1 = PROG 1. PAS 1. 20. ;

evend3 = evol BLEU manu TEMPS TE1 endommagement EF_D ;

dess evend3 TITRE ‘endommagement = F( increment )’ ;

@EXCEL1 EVEND3 ENDOTEMPhr1;

FIN ;

Rédigé par HANGOUANG Duclair Page 20


Rédigé par HANGOUANG Duclair Page 21

Vous aimerez peut-être aussi