Sources Esope | Cast3M

bgrmat
C BGRMAT    SOURCE    CB215821  16/04/21    21:15:21     8920
      SUBROUTINE BGRMAT(SX,SY,HXABS,HYABS,HXORD,HYORD,B)
C
C *********************************************************
C ce sous programme calcule la matrice [BGR] qui
C lie les  courbures aux d.d.l. de l element
C                 -------
C *********************************************************
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
      DIMENSION SX(3),SY(3),HXABS(9),HYABS(9),HXORD(9),HYORD(9),B(4,9)
      DIMENSION B11(9),B12(9),B21(9),B22(9),B31(9),B32(9),B33(9),B34(9)
      DO 10 I=1,9
        B11(I)=SY(2)*HXABS(I)
        B12(I)=SY(3)*HXORD(I)
        B21(I)=-SX(2)*HYABS(I)
        B22(I)=-SX(3)*HYORD(I)
        B31(I)=-SX(2)*HXABS(I)
        B32(I)=-SX(3)*HXORD(I)
        B33(I)=SY(2)*HYABS(I)
        B34(I)=SY(3)*HYORD(I)
10    CONTINUE
C .....................................................
      AIR2=ABS(SX(2)*SY(3)-SX(3)*SY(2))
      CALL ZERO(B,4,9)
      DO 20 I=1,9
         B(1,I)=(B11(I)+B12(I))/AIR2
         B(4,I)=(B21(I)+B22(I))/AIR2
         B(2,I)=(B31(I)+B32(I))/AIR2
         B(3,I)=(B33(I)+B34(I))/AIR2
20    CONTINUE
      RETURN
      END