Numérotation des lignes :

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

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