Télécharger la plaquette Cast3M


fr en

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