Télécharger ellp31.eso

Retour à la liste

Numérotation des lignes :

ellp31
  1. C ELLP31 SOURCE KK2000 14/04/09 21:15:23 8027
  2. SUBROUTINE ELLP31 (R,INP,ZS,ZALFAX,ZALFAR,ZALFAY,ZALFAZ,
  3. * ZG01,ZG11,ZG02,ZG12,
  4. * ZG03,ZG13,ZG23,ZG33,ZG43,ZZ3,
  5. * ZG04,ZG14,ZG24,ZG34,ZG44,ZZ4)
  6. C
  7. IMPLICIT INTEGER(I-N)
  8. IMPLICIT REAL*8 (A-H,O-Y)
  9. IMPLICIT COMPLEX*16 (Z)
  10. C
  11. C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  12. C
  13. C OPERATEUR ELFE LAPLACE POUTRE
  14. C
  15. C CALCUL POUR LA POUTRE N› INP DES FONCTIONS DE GREEN CORRESPONDAN-
  16. C TES EN - TC (N› 1)
  17. C - TORSION (N› 2)
  18. C - FLEXION / OY (N› 3) )
  19. C - FLEXION / OZ (N› 4) ) --> MODELE DE TIMOSHENKO
  20. C AU POINT DE COORDONNEES LOCALES R.
  21. C
  22. C PARAMETRES :
  23. C
  24. C CARACT : TABLEAU DES CARACTERISTIQUE DES POUTRES (12 ,NP)
  25. C R : POINT DE CALCUL DE G ( 0 < R < L)
  26. C INP : NUMERO DE LA POUTRE
  27. C ZS : VALEUR DE S = S0 + I W
  28. C
  29. C SORTIES :
  30. C
  31. C ZG'N''M' : VALEUR DE LA DERIVEE 'N'IEME DELA FONCTION DE GREEN
  32. C NUMERO M.
  33. C
  34. C
  35. C AUTEUR : SAINT-DIZIER
  36. C DATE : 04 JANVIER 1990 (VERSION DU 03 AVRIL 1990)
  37. C
  38. C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  39. C
  40. C
  41. C -- CALCUL DE GREEN 1 : TRACTION COMPRESSION
  42. C --------------------
  43. ZG01 = CMPLX(0.5D0) *EXP(-ZALFAX*R)
  44. ZG11 = CMPLX(-0.5D0)*EXP(-ZALFAX*R)
  45. C
  46. C -- CALCUL DE GREEN 2 : TORSION
  47. C -------
  48. ZG02 = CMPLX(.5D0) *EXP(-ZALFAR*R)
  49. ZG12 = CMPLX(-.5D0)*EXP(-ZALFAR*R)
  50. C
  51. C -- CALCUL DE GREEN 3 : FLEXION / OY
  52. C ------------
  53. ZDELTA= ZZ3*ZZ3 - CMPLX(4.D0)
  54. ZDEL = SQRT(ZDELTA)
  55. C
  56. Z1 = (ZZ3+ZDEL)/CMPLX(2.D0)
  57. Z2 = (ZZ3-ZDEL)/CMPLX(2.D0)
  58. ZL1 = SQRT(Z1)
  59. ZL2 = SQRT(Z2)
  60. PR1 = ZL1*ZALFAY
  61. PR2 = ZL2*ZALFAY
  62. C
  63. C -- SELECTION DES PARTIES REELLES POSITIVES
  64. C
  65. IF (PR1.LT.0.) ZL1 = -ZL1
  66. IF (PR2.LT.0.) ZL2 = -ZL2
  67. C
  68. ZZA = CMPLX(1.D0)/(CMPLX(2.D0)*ZL1*(ZL2*ZL2-ZL1*ZL1))
  69. ZZB = CMPLX(1.D0)/(CMPLX(2.D0)*ZL2*(ZL1*ZL1-ZL2*ZL2))
  70. C
  71. ZE1 = EXP(-ZL1*ZALFAY*R)
  72. ZE2 = EXP(-ZL2*ZALFAY*R)
  73. C
  74. ZG03 = ZZA*ZE1 + ZZB*ZE2
  75. ZG13 = -ZL1*ZZA*ZE1 - ZL2*ZZB*ZE2
  76. ZG23 = ZL1*ZL1*ZZA*ZE1 + ZL2*ZL2*ZZB*ZE2
  77. ZG33 = -ZL1*ZL1*ZL1*ZZA*ZE1 - ZL2*ZL2*ZL2*ZZB*ZE2
  78. ZG43 = ZL1*ZL1*ZL1*ZL1*ZZA*ZE1 + ZL2*ZL2*ZL2*ZL2*ZZB*ZE2
  79. C
  80. C -- CALCUL DE GREEN 4 : FLEXION / OZ
  81. C ------------
  82. ZDELTA= ZZ4*ZZ4 - CMPLX(4.D0)
  83. ZDEL = SQRT(ZDELTA)
  84. C
  85. Z1 = (ZZ4+ZDEL)/CMPLX(2.D0)
  86. Z2 = (ZZ4-ZDEL)/CMPLX(2.D0)
  87. ZL1 = SQRT(Z1)
  88. ZL2 = SQRT(Z2)
  89. PR1 = ZL1*ZALFAZ
  90. PR2 = ZL2*ZALFAZ
  91. C
  92. C -- SELECTION DES PARTIES REELLES POSITIVES
  93. C
  94. IF (PR1.LT.0.) ZL1 = -ZL1
  95. IF (PR2.LT.0.) ZL2 = -ZL2
  96. C
  97. ZZA = CMPLX(1.D0)/(CMPLX(2.D0)*ZL1*(ZL2*ZL2-ZL1*ZL1))
  98. ZZB = CMPLX(1.D0)/(CMPLX(2.D0)*ZL2*(ZL1*ZL1-ZL2*ZL2))
  99. C
  100. ZE1 = EXP(-ZL1*ZALFAZ*R)
  101. ZE2 = EXP(-ZL2*ZALFAZ*R)
  102. C
  103. ZG04 = ZZA*ZE1 + ZZB*ZE2
  104. ZG14 = -ZL1*ZZA*ZE1 - ZL2*ZZB*ZE2
  105. ZG24 = ZL1*ZL1*ZZA*ZE1 + ZL2*ZL2*ZZB*ZE2
  106. ZG34 = -ZL1*ZL1*ZL1*ZZA*ZE1 - ZL2*ZL2*ZL2*ZZB*ZE2
  107. ZG44 = ZL1*ZL1*ZL1*ZL1*ZZA*ZE1 + ZL2*ZL2*ZL2*ZL2*ZZB*ZE2
  108. C
  109. END
  110.  
  111.  
  112.  
  113.  

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