Sources Esope | Cast3M

elpdi3
C ELPDI3    SOURCE    CHAT      05/01/12    23:36:39     5004
      SUBROUTINE ELPDI3(P0,XN0
     &                  ,A,XN1,XN2,XT1,XT2,XLL
     &                  ,XD,XNU
     &                  ,CS2,CB,ISTAT)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-B,D-H,O-Z)
      IMPLICIT COMPLEX*16(C)
************************************************************************
*
*    PLAQUES EN DYNAMIQUE PAR EQUATION INTEGRALE:
*
*       CALCUL DU TERME DE COIN : 12 COEFFICIENTS A PLACER DANS LES
*       RELATIONS
*
************************************************************************
*
*
          DIMENSION P0(2)
          DIMENSION A (2)
*
          DIMENSION XR (2)
          DIMENSION XN0(2)
          DIMENSION XN1(2)
          DIMENSION XN2(2)
*
          DIMENSION XT1(2)
          DIMENSION XT2(2)
*
          DIMENSION XLL(4)
*
          DIMENSION XG (5 )
          DIMENSION CG (5 )
          DIMENSION CS2 (12)
*
* 1-   CALCUL DE GREEN AU COIN
*
        R = ( (A(1)-  P0(1))**2 + ( A(2) - P0(2)) ** 2 ) ** .5
        XR(1)=(A(1)-  P0(1) )/ R
        XR(2)=(A(2)-  P0(2) )/ R
        IF ( ISTAT .EQ. 1)THEN
          CALL ELPGFS(R,XG)
          DO 10 I=1,5
            CG(I) = XG(I)
 10       CONTINUE
        ELSE
          CALL ELPDGF(R,CB,CG)
        ENDIF
*
 
* 2-   COEFFICIENTS SUR LA PREMIERE LIGNE D'EQUATION
*
         CAUX1 = -1* XD* (1. - XNU ) * CG(1)
         CS2(1) =   CAUX1 / ( ( XLL(3) + XLL(4) ) / 2)
         CS2(2) = - CS2(1)
         CS2(3) = - CAUX1 / ( ( XLL(1) + XLL(2) ) / 2)
         CS2(4) = - CS2(3)
*
*       CALCUL DE MNT EN A
*
         CALL ELPDDR(CG,R,XR,XN1,XT1,XT1,XT1,2,CVA1)
         CALL ELPDDR(CG,R,XR,XN2,XT2,XT1,XT1,2,CVA2)
         CMNT1 = - 1* ( XD* (1. - XNU ) ) * CVA1
         CMNT2 = - 1* ( XD* (1. - XNU ) ) * CVA2
           AX =   XLL(2) / (XLL(2) + XLL(3))
         CS2(5) =  AX * ( - CMNT2 + CMNT1)
         CS2(6) = (1 -  AX) * ( - CMNT2 + CMNT1)
 
*
* 3-   COEFFICIENTS SUR LA DEUXIEME LIGNE
*
         CALL ELPDDR(CG,R,XR,XN0,XN0,XN0,XN0,1,CVA )
         CAUX1 = -1* XD* (1. - XNU ) * (- CVA)
         CS2(7) =   CAUX1 / ( ( XLL(3) + XLL(4) ) / 2)
         CS2(8) = - CS2(7)
         CS2(9) = - CAUX1 / ( ( XLL(1) + XLL(2) ) / 2)
         CS2(10) = - CS2(9)
*
*       CALCUL DE MNT,N0 EN A
*
         CALL ELPDDR(CG,R,XR,XN1,XT1,XN0,XN0,3,CVA1)
         CALL ELPDDR(CG,R,XR,XN2,XT2,XN0,XN0,3,CVA2)
         CMNT1  = - 1* ( XD* (1. - XNU ) ) * (- CVA1)
         CMNT2  = - 1* ( XD* (1. - XNU ) ) * (- CVA2  )
         CS2(11) =  AX        * ( - CMNT2 + CMNT1)
         CS2(12) = (1. -  AX) * ( - CMNT2 + CMNT1)
*
       RETURN
      END