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