Télécharger la plaquette Cast3M


fr en

Télécharger ella31.eso

Retour à la liste

Numérotation des lignes : 

ella31

  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