Esope source files

   1 : C ELLA31    SOURCE    CHAT      05/01/12    23:34:31     5004           
   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  = .5D0*EXP(-ZALFAX*R)
  37 :         ZG11  =-.5D0*EXP(-ZALFAX*R)
  38 : C
  39 : C -- CALCUL DE GREEN 2 : TORSION
  40 : C                        -------
  41 :         ZG02  = .5D0*EXP(-ZALFAR*R)
  42 :         ZG12  =-.5D0*EXP(-ZALFAR*R)
  43 : C
  44 : C -- CALCUL DE GREEN 3 : FLEXION (/ OY ET OZ)
  45 : C                        --------------------
  46 :         ZDELTA= ZZ3*ZZ3 - CMPLX(4.D0,0.D0)
  47 :         ZDEL  = SQRT(ZDELTA)
  48 : C
  49 :         Z1    = (ZZ3+ZDEL)/2.D0
  50 :         Z2    = (ZZ3-ZDEL)/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 = 1.D0/(2.D0*ZL1*(ZL2*ZL2-ZL1*ZL1))
  62 :         ZZB = 1.D0/(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  = .5D0*EXP(-ZALFAA*R)
  76 :         ZG14  =-.5D0*EXP(-ZALFAA*R)
  77 : C
  78 :         END
  79 :  
  80 :  
  81 :