Procedures | Cast3M

   1 : * @RAYO     PROCEDUR  AM        09/06/22    21:15:01     6429           
   2 : 
   3 : 'DEBPROC' @RAYO PF*POINT P1*POINT NBC*ENTIER ;
   4 : * Procedure pour generer un maillage rayonnant autour du 
   5 : * fond de fissure
   6 : &DIME = 'VALEUR' DIME ;
   7 : 'SI' (&DIME 'NEG' 2) ;
   8 :    MESS 'PROCEDURE @RAYO : DIMENSION 2 OBLIGATOIRE' ;
   9 :    SAUT LIGNE;
  10 :    'QUITTER' @RAYO;
  11 : 'FINSI';
  12 : OPTI ELEM QUA8 COUL ROUGE;
  13 : PFP1 = P1 'MOINS' PF ;
  14 : 'SI' ((NBC * 3.14159 / 8.) > 4);
  15 :    MESS 'PROCEDURE @RAYO : ELEMENTS TROP DISTORDUS : DIMINUEZ LE '
  16 :         'NOMBRE DE COURONNES !' ;
  17 :    SAUT LIGNE;
  18 :    'QUITTER' @RAYO;
  19 : 'FINSI';
  20 : *
  21 : FRAC = 1./NBC ;
  22 : EFRAC = FRAC / 100. ;
  23 : DENS  FRAC ;
  24 : PFP01 = FRAC * PFP1 ;
  25 : P01 = PF 'PLUS' PFP01 ;
  26 : L = 'DROIT' (NBC-1) P01 P1 ;
  27 : S1 = L 'ROTA' 5 90.  PF ;
  28 : L11 = S1 'COTE' 1 ;
  29 : L1EXT = S1 'COTE' 2 ;
  30 : L13 = S1 'COTE' 3 ;
  31 : L1INT = S1 'COTE' 4 ;
  32 : PEXTF = L1EXT 'POINT' 'FINAL' ;
  33 : LR = 'DROIT' 1 PF P01 ;
  34 : LLG = L11 'ET' LR ;
  35 : LRG = LR 'TOUR' 90. PF ;
  36 : LT = LRG 'ET' L1INT 'ET' LR ; ELIM LT EFRAC;
  37 : OPTI ELEM TRI6 ;
  38 : VEC1 = PF 'MOIN' (LT 'POIN' 3);
  39 : PT1 = PF 'MOIN' (VEC1 / 2.);
  40 : VEC2 = PF 'MOIN' (LT 'POIN' 4);
  41 : PT2 = PF 'MOIN' (VEC2 / 2.);
  42 : VEC3 = PF 'MOIN' (LT 'POIN' 5);
  43 : PT3 = PF 'MOIN' (VEC3 / 2.);
  44 : VEC4 = PF 'MOIN' (LT 'POIN' 6);
  45 : PT4 = PF 'MOIN' (VEC4 / 2.);
  46 : *
  47 : LT2 = 'CHANGER' 'SEG2' LT;
  48 : EL1 = 'MANU' 'TRI6' PF (LT2 'POIN' 2) (LT2 'POIN' 3) (LT2 'POIN' 4)
  49 :  (LT2 'POIN' 5) PT1 ; 
  50 : EL2 = MANU TRI6 PF PT1 (LT2 'POIN' 5) (LT2 'POIN' 6) (LT2 'POIN' 7) PT2;
  51 : EL3 = MANU TRI6 PF PT2 (LT2 'POIN' 7) (LT2 'POIN' 8) (LT2 'POIN' 9) PT3;
  52 : EL4 = MANU TRI6 PF PT3 (LT2 'POIN' 9) (LT2 'POIN' 10) (LT2 'POIN' 11) 
  53 :   PT4 ; 
  54 : EL5 = MANU TRI6 PF PT4 (LT2 'POIN' 11) (LT2 'POIN' 12) 
  55 :  ((LT2 'ELEM' 12) POINT 'FINAL') (LT2 'POIN' 14) ; 
  56 : S0 = EL1 'ET' EL2 'ET' EL3 'ET' EL4 'ET' EL5;
  57 : S = S1 'ET' S0 ;
  58 : S2 L2EXT SI1 LL2 = S L1EXT S0 LLG SYME 'DROIT' PF PEXTF ;
  59 : P2EXTF = L2EXT POINT 'INITIAL' ;
  60 : 'ELIM'  S  S2 EFRAC ;
  61 : S3 L3EXT SI2 LL3 = S2 L2EXT SI1 LL2 SYME 'DROIT' PF P2EXTF ;
  62 : S3 = S3 COUL JAUNE;
  63 : P3EXTF = L3EXT POINT 'FINAL' ;
  64 : 'ELIM' S2  S3 EFRAC ;
  65 : S4 L4EXT SI3 LLD = S3 L3EXT SI2 LL3 SYME 'DROIT' PF P3EXTF ;
  66 : S4 = S4 COUL JAUNE;
  67 : 'ELIM' S3 S4 EFRAC ;
  68 : SINT = S0 'ET' SI1 'ET' SI2 'ET' SI3;
  69 : ST =  S 'ET' S2 'ET' S3 'ET' S4 ;
  70 : 'FINPROC' ST SINT L1EXT L2EXT L3EXT L4EXT LLD LLG;
  71 : 
  72 :  
  73 :