Télécharger ella31.eso

Retour à la liste

Numérotation des lignes :

  1. C ELLA31 SOURCE KK2000 14/04/09 21:15:18 8027
  2. SUBROUTINE ELLA31 (R,ZALFAX,ZALFAR,ZALFAF,ZALFAA,
  3. * ZG01,ZG11,ZG02,ZG12,
  4. * ZG03,ZG13,ZG23,ZG33,ZG43,ZZ3,
  5. * ZG04,ZG14)
  6. C
  7. IMPLICIT INTEGER(I-N)
  8. IMPLICIT REAL*8 (A-H,O-Y)
  9. IMPLICIT COMPLEX*16 (Z)
  10. C
  11. C OPERATEUR ELFE LAPLACE ACOU
  12. C
  13. C CALCUL POUR LA POUTRE N› INP DES FONCTIONS DE GREEN CORRESPONDAN-
  14. C TES EN - TC (Num 1)
  15. C - TORSION (Num 2)
  16. C - FLEXION / OY ET / OZ (Num 3) --> MODELE DE TIMOSCHENKO
  17. C - PRESSION ACOUSTIQUE (Num 4)
  18. C AU POINT DE COORDONNEES LOCALES R.
  19. C
  20. C PARAMETRES :
  21. C
  22. C CARACT : TABLEAU DES CARACTERISTIQUE DES POUTRES (10 ,NP)
  23. C R : POINT DE CALCUL DE G ( 0 < R < L)
  24. C
  25. C SORTIES :
  26. C
  27. C ZG'N''M' : VALEUR DE LA DERIVEE 'N'IEME DE LA FONCTION DE GREEN
  28. C NUMERO M.
  29. C
  30. C
  31. C AUTEURS : SAINT-DIZIER ET GORCY
  32. C DATE : 24 JANVIER 1991
  33. C
  34. C -- CALCUL DE GREEN 1 : TRACTION COMPRESSION
  35. C --------------------
  36. ZG01 = CMPLX(.5D0) *EXP(-ZALFAX*R)
  37. ZG11 = CMPLX(-.5D0)*EXP(-ZALFAX*R)
  38. C
  39. C -- CALCUL DE GREEN 2 : TORSION
  40. C -------
  41. ZG02 = CMPLX(.5D0) *EXP(-ZALFAR*R)
  42. ZG12 = CMPLX(-.5D0)*EXP(-ZALFAR*R)
  43. C
  44. C -- CALCUL DE GREEN 3 : FLEXION (/ OY ET OZ)
  45. C --------------------
  46. ZDELTA= ZZ3*ZZ3 - (4.D0,0.D0)
  47. ZDEL = SQRT(ZDELTA)
  48. C
  49. Z1 = (ZZ3+ZDEL)/CMPLX(2.D0)
  50. Z2 = (ZZ3-ZDEL)/CMPLX(2.D0)
  51. ZL1 = SQRT(Z1)
  52. ZL2 = SQRT(Z2)
  53. PR1 = ZL1*ZALFAF
  54. PR2 = ZL2*ZALFAF
  55. C
  56. C -- SELECTION DES PARTIES REELLES POSITIVES
  57. C
  58. IF (PR1.LT.0.D0) ZL1 = -ZL1
  59. IF (PR2.LT.0.D0) ZL2 = -ZL2
  60. C
  61. ZZA = CMPLX(1.D0)/(CMPLX(2.D0)*ZL1*(ZL2*ZL2-ZL1*ZL1))
  62. ZZB = CMPLX(1.D0)/(CMPLX(2.D0)*ZL2*(ZL1*ZL1-ZL2*ZL2))
  63. C
  64. ZE1 = EXP(-ZL1*ZALFAF*R)
  65. ZE2 = EXP(-ZL2*ZALFAF*R)
  66. C
  67. ZG03 = ZZA*ZE1 + ZZB*ZE2
  68. ZG13 = -ZL1*ZZA*ZE1 - ZL2*ZZB*ZE2
  69. ZG23 = ZL1*ZL1*ZZA*ZE1 + ZL2*ZL2*ZZB*ZE2
  70. ZG33 = -ZL1*ZL1*ZL1*ZZA*ZE1 - ZL2*ZL2*ZL2*ZZB*ZE2
  71. ZG43 = ZL1*ZL1*ZL1*ZL1*ZZA*ZE1 + ZL2*ZL2*ZL2*ZL2*ZZB*ZE2
  72. C
  73. C -- CALCUL DE GREEN 4 : PRESSION ACOUSTIQUE
  74. C -------------------
  75. ZG04 = CMPLX(.5D0) *EXP(-ZALFAA*R)
  76. ZG14 = CMPLX(-.5D0)*EXP(-ZALFAA*R)
  77. C
  78. END
  79.  
  80.  
  81.  
  82.  
  83.  

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