Télécharger la plaquette Cast3M


fr en

Télécharger elpop1.eso

Retour à la liste

Numérotation des lignes : 

elpop1

  1. C ELPOP1    SOURCE    CHAT      05/01/12    23:37:51     5004
  2.       SUBROUTINE ELPOP1(P0,XN0,QI,XN,XT,XD,XNU,XOP)
  3.       IMPLICIT INTEGER(I-N)
  4.       IMPLICIT REAL*8(A-H,O-Z)
  5. ************************************************************************
  6. *
  7. *    PLAQUES EN STATIQUE PAR EQUATION INTEGRALE:
  8. *
  9. *       CALCUL DES OPERATEURS PLAQUE ET DE LEURS DERIVEES
  10. *       POUR UN COUPLE DE PTS P0 , QI
  11. *
  12. ************************************************************************
  13. *
  14. *
  15.           DIMENSION P0(2)
  16.           DIMENSION QI(2)
  17.           DIMENSION XOP(8)
  18. *
  19.           DIMENSION XN (2)
  20.           DIMENSION XN0(2)
  21. *
  22.           DIMENSION XT (2)
  23.           DIMENSION XT0(2)
  24.           DIMENSION XG (5)
  25.           DIMENSION XR (2)
  26. *
  27. * 1- CALCUL DE R ET DU VECTEUR UNITAIRE SUIVANT R
  28. *
  29.          R= ( ( QI(1) - P0(1) ) ** 2 + (QI(2) - P0(2) ) ** 2 ) **.5
  30.          XR(1) = (QI(1) - P0(1)) / R
  31.          XR(2) = (QI(2) - P0(2)) / R
  32. *
  33. *
  34. *
  35. *
  36. * 2- ON CALCULE LES FONCTIONS DE GREEN ET LEURS DERIVEES,R
  37. *
  38.        CALL ELPGFS(R,XG)
  39. *
  40. * 3-   CALCUL DE V
  41. *
  42.        XOP(1) = XG(1)
  43. *
  44. * 4-   CALCUL DE V,N
  45. *
  46.        CALL ELPDER(XG,R,XR,XN,XN,XN,XN,1,AUX1)
  47.        XOP(2) = AUX1
  48. *
  49. * 5-   CALCUL DE MN= -D * (V,NN + NU* V,TT)
  50. *
  51.        CALL ELPDER(XG,R,XR,XN,XN,XN,XN,2,AUX1)
  52.        CALL ELPDER(XG,R,XR,XT,XT,XN,XN,2,AUX2)
  53.        XOP(3) = -1* XD* ( AUX1 +  XNU*AUX2)
  54. *
  55. * 6-   CALCUL DE KN= -D * (V,NNN + (2 - NU)* V,NTT)
  56. *
  57.        CALL ELPDER(XG,R,XR,XN,XN,XN,XN,3,AUX1)
  58.        CALL ELPDER(XG,R,XR,XN,XT,XT,XN,3,AUX2)
  59.        XOP(4) = -1* XD* ( AUX1 +  (2. - XNU) *AUX2)
  60. *****
  61. *      TERMES DERIVES
  62. *****
  63. * 7-   CALCUL DE V,N0
  64. *
  65. *
  66. *
  67. *
  68.        CALL ELPDER(XG,R,XR,XN0,XN,XN,XN,1,AUX1)
  69.        XOP(5) = - AUX1
  70. *
  71. * 8-   CALCUL DE V,NN0
  72. *
  73.        CALL ELPDER(XG,R,XR,XN,XN0,XN,XN,2,AUX1)
  74.        XOP(6) =  - AUX1
  75. *
  76. * 9-   CALCUL DE MN,N0
  77. *
  78.        CALL ELPDER(XG,R,XR,XN,XN,XN0,XN,3,AUX1)
  79.        CALL ELPDER(XG,R,XR,XT,XT,XN0,XN,3,AUX2)
  80.        XOP(7) = - 1* XD* ( - AUX1 -  XNU*AUX2)
  81. *
  82. * 9-   CALCUL DE KN,N0
  83. *
  84.        CALL ELPDER(XG,R,XR,XN,XN,XN,XN0,4,AUX1)
  85.        CALL ELPDER(XG,R,XR,XN,XT,XT,XN0,4,AUX2)
  86.        XOP(8) = -1* XD* (-  AUX1 -  (2. - XNU) *AUX2)
  87.  
  88.        RETURN
  89.       END
  90.  
  91.  

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