Esope source files

   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 : *                                                                               
  37 :        CALL ELPGFS(R,XG)                                                        
  38 : *                                                                               
  39 : * 3-   CALCUL DE V                                                              
  40 : *                                                                               
  41 :        XOP(1) = XG(1)                                                           
  42 : *                                                                               
  43 : * 4-   CALCUL DE V,N                                                            
  44 : *                                                                               
  44 :        CALL ELPDER(XG,R,XR,XN,XN,XN,XN,1,AUX1)                                  
  45 :        XOP(2) = AUX1                                                            
  46 : *                                                                               
  47 : * 5-   CALCUL DE MN= -D * (V,NN + NU* V,TT)                                     
  48 : *                                                                               
  48 :        CALL ELPDER(XG,R,XR,XN,XN,XN,XN,2,AUX1)                                  
  48 :        CALL ELPDER(XG,R,XR,XT,XT,XN,XN,2,AUX2)                                  
  49 :        XOP(3) = -1* XD* ( AUX1 +  XNU*AUX2)                                     
  50 : *                                                                               
  51 : * 6-   CALCUL DE KN= -D * (V,NNN + (2 - NU)* V,NTT)                             
  52 : *                                                                               
  52 :        CALL ELPDER(XG,R,XR,XN,XN,XN,XN,3,AUX1)                                  
  52 :        CALL ELPDER(XG,R,XR,XN,XT,XT,XN,3,AUX2)                                  
  53 :        XOP(4) = -1* XD* ( AUX1 +  (2. - XNU) *AUX2)                             
  54 : *****                                                                           
  55 : *      TERMES DERIVES                                                           
  56 : *****                                                                           
  57 : * 7-   CALCUL DE V,N0                                                           
  58 : *                                                                               
  59 : *                                                                               
  60 : *                                                                               
  61 : *                                                                               
  61 :        CALL ELPDER(XG,R,XR,XN0,XN,XN,XN,1,AUX1)                                 
  62 :        XOP(5) = - AUX1                                                          
  63 : *                                                                               
  64 : * 8-   CALCUL DE V,NN0                                                          
  65 : *                                                                               
  65 :        CALL ELPDER(XG,R,XR,XN,XN0,XN,XN,2,AUX1)                                 
  66 :        XOP(6) =  - AUX1                                                         
  67 : *                                                                               
  68 : * 9-   CALCUL DE MN,N0                                                          
  69 : *                                                                               
  69 :        CALL ELPDER(XG,R,XR,XN,XN,XN0,XN,3,AUX1)                                 
  69 :        CALL ELPDER(XG,R,XR,XT,XT,XN0,XN,3,AUX2)                                 
  70 :        XOP(7) = - 1* XD* ( - AUX1 -  XNU*AUX2)                                  
  71 : *                                                                               
  72 : * 9-   CALCUL DE KN,N0                                                          
  73 : *                                                                               
  73 :        CALL ELPDER(XG,R,XR,XN,XN,XN,XN0,4,AUX1)                                 
  73 :        CALL ELPDER(XG,R,XR,XN,XT,XT,XN0,4,AUX2)                                 
  74 :        XOP(8) = -1* XD* (-  AUX1 -  (2. - XNU) *AUX2)                           
  75 :                                                                                 
  76 :        RETURN                                                                   
  77 :       END                                                                       
  78 :