Sources Esope | Cast3M

bmato
C BMATO     SOURCE    CB215821  16/04/21    21:15:22     8920
      SUBROUTINE BMATO(SX,SY,HXABS,HYABS,HXORD,HYORD,B)
C
C *********************************************************
C ce sous programme calcule la matrice [B] qui lie les
C 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(3,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,3,9)
      DO 20 I=1,9
      B(1,I)=(B11(I)+B12(I))/AIR2
      B(2,I)=(B21(I)+B22(I))/AIR2
      B(3,I)=(B31(I)+B32(I)+B33(I)+B34(I))/AIR2
20    CONTINUE
      RETURN
      END