Télécharger la plaquette Cast3M


fr en

Télécharger ouvfiss.procedur

Retour à la liste

Numérotation des lignes : 

  1. * OUVFISS   PROCEDUR  AM        16/04/04    21:15:01     8883           
  2.  
  3. DEBP OUVFISS TAB1*TABLE IPAS/ENTIER OPT1/MOT;
  4. *
  5. * Procédure OUVFISS Christian La Borderie & M. Matallah
  6. * Cette procédure permet de calculer l'ouverture de fissure à
  7. * Partir de la table de calcul PASAPAS si le modèle utilisé est
  8. * Régularisé par la méthode de Hillerborg
  9. * Les résultats sont dans tab1.ouv
  10. *
  11. * Améliorations à apporter : 
  12. *     utilisation de l'opérateur taille (à valider sur l'ensemble des élémen
  13. *     identification de l'ouverture normale et du glissement
  14. *
  15. *
  16. * Si vous utilisez cette procédure, vous êtes priés de citer :
  17. * M. Matallah, C La Borderie and O. Maurel 
  18. *  'A practical method to estimate crack openings in concrete structures'
  19. *  International Journal for Numerical and Analytical Methods in Geomechanics, 2
  20. *  doi:10.1002/nag.876
  21.  
  22. *
  23. *
  24. *
  25.   SI ('EGA' OPT1  'TR');
  26.     LTRA=VRAI;
  27.     MESSAGE ' Extraction de la trace de l ouverture de fissure';
  28.   SINON;
  29.     LTRA=FAUX;
  30.     LPRIN=FAUX;
  31.     MESSAGE ' Extraction du tenseur d ouverture de fissure';
  32.     SI ('EGA' OPT1 'PRIN');
  33.       LPRIN=VRAI;
  34.       MESSAGE 'Calcul des valeurs principales';
  35.     FINSI;
  36.   FINSI;
  37.    IDIM=VALE DIME;
  38.    SI ('EGA' IDIM 2);
  39.         MOD1='EXTR' TAB1.MODELE 'ELEM' 'TRI3' 'QUA4' 'TRI6' 'QUA8' ;
  40.    SINON; SI ('EGA' IDIM 3);
  41.         MOD1='EXTR' TAB1.MODELE 'ELEM' 
  42.              'CUB8' 'CU20' 'PRI6' 'PR15' 'PYR5' 'PY13' 'TET4' 'TE10';
  43.         SINON;
  44.             ERREUR 'La dimension doit être égale à 2 ou 3';
  45.             QUITTER OUVFISS;
  46.         FINSI;
  47.    FINSI;      
  48.    CHAM1='MANU' 'CHML' MOD1 'SCAL' 1. 'TYPE' 'DEFORMATIONS' 'STRESSES';
  49.    CHAM2='INTG' CHAM1 MOD1 'ELEM' 'SCAL';
  50.    CHAM3=('CHAN' 'STRESSES' MOD1 CHAM2) 'CHAN' 'TYPE' 'DEFORMATIONS';
  51.    CHH=CHAM3**(1./(VALE DIME)) ;
  52.    MAT1='REDU' TAB1 . 'CARACTERISTIQUES' MOD1;
  53.    MAILL1='EXTR' MOD1 'MAILL';
  54.    SI (EXISTE IPAS);
  55.      N1=1;
  56.      MESSAGE 'Pour le pas numéro ' IPAS;
  57.    SINON;
  58.       LECT1='INDEX' (TAB1.TEMPS);
  59.       N1='DIME' LECT1;
  60.       MESSAGE 'Pour tous les pas de temps de la table';
  61.    FINSI;
  62.    TAB1.OUV='TABLE';
  63.    REPETER BOU1 N1; 
  64.      SI (EXISTE IPAS);
  65.        I=IPAS;
  66.      SINON;
  67.         I=LECT1 . &BOU1;
  68.      FINSI;  
  69.      DEP1='REDU' TAB1.DEPLACEMENTS.I MAILL1; 
  70.      EPS1='EPSI' TAB1.DEPLACEMENTS.I MOD1;
  71.      SIG1='REDU' TAB1.CONTRAINTES.I MOD1;
  72.      SIG0I='ELAS' EPS1 MOD1 MAT1;
  73.      SIGAN=SIG0I - SIG1;
  74.      OUF0='ELAS' SIGAN MOD1 MAT1;
  75.      LM1=EXTR OUF0 COMP;
  76.      NCOMP=DIME LM1;
  77.      LM2=MOTS;
  78.      REPETER BOUCOMP NCOMP;
  79.        LM2=INSE LM2  &BOUCOMP 'SCAL';
  80.      FIN BOUCOMP;
  81.      OUF1=OUF0 * CHH LM1 LM2 LM1;
  82.      SI LTRA;
  83.        SI (IDIM EGA 2);
  84.           TROUF1=(('EXCO' OUF1 'EPXX' 'SCAL')
  85.              +('EXCO' OUF1 'EPYY' 'SCAL')) ;
  86.        FINSI;
  87.        SI (IDIM EGA 3);
  88.           TROUF1=(('EXCO' OUF1 'EPXX' 'SCAL')
  89.             +('EXCO' OUF1 'EPYY' 'SCAL') +
  90.              ('EXCO' OUF1 'EPZZ' 'SCAL'));
  91.        FINSI;
  92.        TAB1.OUV . I=TROUF1;
  93.      SINON;
  94.       SI LPRIN;
  95.         OUF2=PRIN OUF1 MOD1;
  96.         TAB1.OUV . I=OUF2;
  97.       SINON;
  98.         TAB1.OUV . I=OUF1;
  99.       FINSI;    
  100.      FINSI;
  101.      I=I+1;
  102.   FIN BOU1;
  103. FINP;
  104.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales