Télécharger la plaquette Cast3M


fr en

Télécharger shfdt.procedur

Retour à la liste

Numérotation des lignes : 

  1. * SHFDT     PROCEDUR  OF166741  23/02/07    11:58:40     11590          
  2.  
  3. DEBPROC SHFDT TABDON ;
  4. ************************************************************************
  5. *    CALCUL DE FONCTION DE TRANSFERT D'UN CHAMP D'ONDE SH INCLINEE
  6. *    PAR RAPPORT A UN POINT DE REFERENCE
  7. *
  8. *    TABDON : TABLE D'ENTREE AVEC LES INDICES SUIVANTES :
  9.  
  10. *      PDC*POINT: POINT DE CONTROLE A LA SURFACE DU SOL
  11. *      DIR*MOT: DIRECTION DU MOUVEMENT DES PARTICULES DE SOL
  12. *               DIR = UX : ONDE SH SE PROPAGEANT DANS LE PLAN YZ
  13. *               DIR = UY : ONDE SH SE PROPAGEANT DANS LE PLAN XZ
  14. *      ANGLE*FLOTTANT: ANGLE D'INCIDENCE PAR RAPPORT A L'AXE OZ
  15. *      VS*FLOTTANT: VITESSE DE PROPAGATION D'ONDE S
  16. *      PDF*POINT: POINT DE CIBLE QUELCONQUE DANS LE SOL
  17. *      KSI*FLOTTANT: TAUX D'AMORTISSEMENT DU SOL, CONSTANT PAR DEFAUT,
  18. *                    VARIABLE SELON LE MODELE DE RAYLEIGH SI F1 ET F2
  19. *                    EXISTENT
  20. *      F1*FLOTTANT, F2*FLOTTANT: FREQUENCES POUR LE MODELE DE RAYLEIGH
  21. ************************************************************************
  22.  
  23. ANGLE = TABDON.ANGLE ;
  24. VS    = TABDON.VS ;
  25. KSI   = TABDON.KSI ;
  26. LF    = TABDON.LF ;
  27.  
  28. F0 = EXTR LF 1 ;
  29. SI (F0 < 1E-5) ;
  30. I_ENLE = 1 ;
  31. LF = ENLE LF 1;
  32. FINSI ;
  33.  
  34. LOMEG = 2*PI*LF ;
  35. NF    = DIME LF ;
  36. L1    = PROG NF * 1.0 ;
  37. LKSI  = PROG NF * KSI ;
  38.  
  39. SI ((EXIS TABDON F1) ET (EXIS TABDON F2)) ;
  40. OM1 = 2*PI*TABDON.F1 ;
  41. OM2 = 2*PI*TABDON.F2 ;
  42. A0  = 2*KSI*OM1*OM2/(OM1 + OM2) ;
  43. A1  = 2*KSI/(OM1 + OM2) ;
  44. LKSI = (A0*L1 / (2*LOMEG)) + (A1*LOMEG/2.) ;
  45. FINSI;
  46.  
  47. EV1 = EVOL COMP REIM LF L1 (2*LKSI) ;
  48. EV1 = 'RIMP' EV1 ;
  49.  
  50. * Calcul de la racine carre de EV1
  51.  
  52. L2  = (EXTR EV1 ORDO 1)**0.5 ;
  53. L3  = (EXTR EV1 ORDO 2)/2 ;
  54. EV2 = EVOL COMP MOPH LF L2 L3 ;
  55.  
  56. EV2 = 'RIMP' EV2 ;
  57. EVVS = VS*EV2 ;
  58.  
  59. L0     = PROG NF * 0.0 ;
  60. EVI    = EVOL COMP REIM LF L0 L1 ;
  61. EVOMEG = EVOL COMP REIM LF LOMEG L0 ;
  62.  
  63. X0 Y0 Z0 = COOR TABDON.PDC ;
  64. X1 Y1 Z1 = COOR TABDON.PDF ;
  65.  
  66. DX = X1 - X0 ;
  67. DY = Y1 - Y0 ;
  68. DZ = Z1 - Z0 ;
  69.  
  70. SI (EGA TABDON.DIR 'UX') ;
  71.  
  72. EVA = -1*EVI*EVOMEG*DY*(SIN ANGLE)/EVVS ;
  73. EVB = -1*EVI*EVOMEG*DZ*(COS ANGLE)/EVVS ;
  74. EVC =    EVI*EVOMEG*DZ*(COS ANGLE)/EVVS ;
  75.  
  76. SINON;
  77.   SI (EGA TABDON.DIR 'UY') ;
  78.   EVA = -1*EVI*EVOMEG*DX*(SIN ANGLE)/EVVS ;
  79.   EVB = -1*EVI*EVOMEG*DZ*(COS ANGLE)/EVVS ;
  80.   EVC =    EVI*EVOMEG*DZ*(COS ANGLE)/EVVS ;
  81.   FINSI ;
  82. FINSI ;
  83.  
  84. TITR 'EVA1' ;
  85. LMOA = EXP (EXTR EVA ORDO 1) ;
  86. LPHA = 180*(EXTR EVA ORDO 2)/PI ;
  87. EVA1 = EVOL COMP MOPH LF LMOA LPHA ;
  88. EVA1 = 'RIMP' EVA1 ;
  89.  
  90. TITR 'EVB1' ;
  91. LMOB = EXP (EXTR EVB ORDO 1) ;
  92. LPHB = 180*(EXTR EVB ORDO 2)/PI ;
  93. EVB1 = EVOL COMP MOPH LF LMOB LPHB ;
  94. EVB1 = 'RIMP' EVB1 ;
  95.  
  96. TITR 'EVC1' ;
  97. LMOC = EXP (EXTR EVC ORDO 1) ;
  98. LPHC = 180*(EXTR EVC ORDO 2)/PI ;
  99. EVC1 = EVOL COMP MOPH LF LMOC LPHC ;
  100. EVC1 = 'RIMP' EVC1 ;
  101.  
  102. TITR 'FDT ONDE SH' ;
  103. FDTSH = (EVA1 * (EVB1 + EVC1))/2. ;
  104. FDTSH = 'RIMP' FDTSH ;
  105.  
  106. SI (EGA I_ENLE 1) ;
  107. LF    = EXTR FDTSH ABSC ;
  108. LMO   = EXTR FDTSH ORDO 1 ;
  109. LPH   = EXTR FDTSH ORDO 2 ;
  110. LF    = (PROG F0) ET LF ;
  111. LMO   = (PROG (EXTR LMO 1 )) ET LMO ;
  112. LPH   = (PROG (EXTR LPH 1 )) ET LPH ;
  113. FDTSH = EVOL COMP MOPH LF LMO LPH ;
  114. FINSI ;
  115.  
  116. FINPROC FDTSH ;
  117. ************************************************************************
  118.  
  119.  
  120.  

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