C BMATO     SOURCE    CHAT      05/01/12    21:41:12     5004           
      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=DABS(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