Télécharger la plaquette Cast3M


fr en

Télécharger @coutor1.procedur

Retour à la liste

Numérotation des lignes : 

  1. * @COUTOR1  PROCEDUR  CB215821  15/12/02    21:15:01     8719           
  2. *-------------------------------------------------                              
  3. ******          PROCEDURE @COUTOR1          ******
  4. *-------------------------------------------------                              
  5. *
  6. *
  7. *    CETTE PROCEDURE A ETE MISE GRACIEUSEMENT
  8. *   A DISPOSOTION DE LA COMMUNAUTE  CASTEM2000
  9. *    PAR M. LIBEYRE ( CEA/DSM/DRFC )
  10. *
  11. *     TEL : ( 33 1 ) 42 25 46 03
  12. *
  13. *-------------------------------------------------                              
  14. DEBPROC @COUTOR1 IEL*MAILLAGE
  15.                 CHT*CHPOINT CHN*CHPOINT CHB*CHPOINT ;
  16.  
  17. P2 = IEL POIN INITIAL ; 
  18. P3 = IEL POIN FINAL ;
  19. DIMGEO1 = VALEUR DIME ;
  20. SI (DIMGEO1 > 2);
  21.  NXI2 = EXTR CHN NX P2 ;
  22.  NYI2 = EXTR CHN NY P2 ;  
  23.  NZI2 = EXTR CHN NZ P2 ;
  24.  VN2 = NXI2 NYI2 NZI2 ;  
  25.  TXI2 = EXTR CHT TX P2 ;
  26.  TYI2 = EXTR CHT TY P2 ;  
  27.  TZI2 = EXTR CHT TZ P2 ;
  28.  VT2 = TXI2 TYI2 TZI2 ;  
  29.  BXI2 = EXTR CHB BX P2 ;
  30.  BYI2 = EXTR CHB BY P2 ;  
  31.  BZI2 = EXTR CHB BZ P2 ;
  32.  VB2 = BXI2 BYI2 BZI2 ;  
  33.  NXI3 = EXTR CHN NX P3 ;
  34.  NYI3 = EXTR CHN NY P3 ;  
  35.  NZI3 = EXTR CHN NZ P3 ;
  36.  VN3 = NXI3 NYI3 NZI3 ;  
  37.  TXI3 = EXTR CHT TX P3 ;
  38.  TYI3 = EXTR CHT TY P3 ;  
  39.  TZI3 = EXTR CHT TZ P3 ;
  40.  VT3 = TXI3 TYI3 TZI3 ;
  41.  BXI3 = EXTR CHB BX P3 ;
  42.  BYI3 = EXTR CHB BY P3 ;  
  43.  BZI3 = EXTR CHB BZ P3 ;
  44.  VB3 = BXI3 BYI3 BZI3 ;
  45.  XR = VT3 PSCA VT2 ;  
  46.  YR = VT3 PSCA VN2 ;
  47. * MESS 'XR =' XR   'YR =' YR ;
  48.  ALPHA = ATG YR XR ;
  49. * MESS 'ALPHA =' ALPHA ;
  50.  DS2 = NORM (MOIN P2 P3) ;
  51.  RR = (ALPHA*PI/180.)/DS2 ;
  52.  SI (RR NEG 0.) ;
  53.   R = 1./RR ;
  54.  SINON ;
  55.   R = 1.E99 ;
  56.  FINSI ;
  57. * MESS 'R =' R ;
  58.  XT = VB2 PSCA VB3 ;
  59.  YT = VB2 PSCA VN3 ;
  60. * MESS 'XT =' XT   'YT =' YT ;
  61.  BETA = ATG YT XT ;  
  62. * MESS 'BETA =' BETA ;
  63.  TT = -1*(BETA*PI/180.)/DS2 ;
  64.  SI (TT NEG 0.) ;
  65.   T = 1./TT ;
  66.  SINON ;
  67.   T = 1.E99 ;
  68.  FINSI ;
  69. * MESS 'T =' T ;
  70.  
  71. SINON;  
  72.  NXI2 = EXTR CHN NX P2 ;
  73.  NYI2 = EXTR CHN NY P2 ;  
  74.  VN2 = NXI2 NYI2 ;  
  75.  TXI2 = EXTR CHT TX P2 ;
  76.  TYI2 = EXTR CHT TY P2 ;  
  77.  VT2 = TXI2 TYI2 ;  
  78.  BXI2 = EXTR CHB BX P2 ;
  79.  BYI2 = EXTR CHB BY P2 ;  
  80.  VB2 = BXI2 BYI2 ;  
  81.  NXI3 = EXTR CHN NX P3 ;
  82.  NYI3 = EXTR CHN NY P3 ;  
  83.  VN3 = NXI3 NYI3 ;  
  84.  TXI3 = EXTR CHT TX P3 ;
  85.  TYI3 = EXTR CHT TY P3 ;  
  86.  VT3 = TXI3 TYI3 ;
  87.  BXI3 = EXTR CHB BX P3 ;
  88.  BYI3 = EXTR CHB BY P3 ;  
  89.  VB3 = BXI3 BYI3 ;
  90.  XR = VT3 PSCA VT2 ;  
  91.  YR = VT3 PSCA VN2 ;
  92. * MESS 'XR =' XR   'YR =' YR ;
  93.  ALPHA = ATG YR XR ;
  94. * MESS 'ALPHA =' ALPHA ;
  95.  DS2 = NORM (MOIN P2 P3) ;
  96.  RR = (ALPHA*PI/180.)/DS2 ;
  97.  SI (RR NEG 0.) ;
  98.   R = 1./RR ;
  99.  SINON ;
  100.   R = 1.E99 ;
  101.  FINSI ;
  102. * MESS 'R =' R ;
  103.  BETA = 0.;
  104.  T = 0.;
  105. FINSI;
  106. FINPROC DS2 R T ALPHA BETA ;
  107.  
  108.  

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