Esope source files

   1 : C CUCALB    SOURCE    CHAT      07/10/22    21:15:37     5921           
   2 :       SUBROUTINE CUCALB(BKSIP,XNOE,B,AJAC)
   3 :       IMPLICIT REAL*8(A-H,O-Z)
   4 :       IMPLICIT INTEGER (I-N)
   5 :  
   6 : * 
   7 : *     ------------------------------------------------------------------ 
   8 : *      CALCUL DE LA MATRICE B AU POINT D INTEGRATION DONNNE 
   9 : * 
  10 : *                                                    H. BUNG     06-98 
  11 : *     ------------------------------------------------------------------ 
  12 : * 
  13 : *  ENTREES : 
  14 : *     XNOE(24) : COORDONNEES DES NOEUDS 
  15 : *     IPINT      : NO DU POINT D INTEGRATION 
  16 : *     MOT        : NOM DE L ELEMENT FINI 
  17 : *  SORTIE 
  18 : *     B(NBLIB,NBN) : MATRICE B 
  19 : *     AJAC         : JACOBIEN 
  20 : * 
  21 : *
  22 :       DIMENSION XNOE(24),BKSIP(3,8),B(3,8)
  23 :  
  24 : *      REAL *8 XNOE(24),BKSIP(3,8) 
  25 : *      REAL *8 B(3,8),AJAC 
  26 : * 
  27 : *----   VARIABLES LOCALES 
  28 : * 
  29 : *      REAL *8 DJ(3,3),UJ(3,3) 
  30 : *      INTEGER  LRET,I,J,K,NBN 
  31 :        DIMENSION DJ(3,3),UJ(3,3)
  32 : * 
  33 : *      REAL *8  ZER, UN 
  34 : * 
  35 :       NBN=8 
  36 :       ZER=0D0 
  37 :       UN=1.D0 
  38 : * 
  39 : *---   DJ = BKSIP * TRANSPOSE(XNOE) 
  40 : * 
  41 :       DO I=1,3 
  42 :         DO J=1,3 
  43 :         DJ(I,J) = 0. 
  44 :         END DO 
  45 :       END DO 
  46 :       DO I=1,3 
  47 :         DO J=1,3 
  48 :            DO K=1,NBN 
  49 : C             DJ(J,I)=DJ(J,I)+BKSIP(J,K)*XNOE(I,K) 
  50 :              DJ(J,I)=DJ(J,I)+BKSIP(J,K)*XNOE((K-1)*3+I) 
  51 :            END DO 
  52 :         END DO 
  53 :       END DO 
  54 : * 
  55 : *-----   UJ(J,I)  MATRICE INVERSE DE DJ(J,I) 
  56 : * 
  57 :        
  57 :       CALL INV33(DJ,UJ,AJAC,LRET) 
  58 : C TEST SI ELEMENT TROP DEFORME: CROISEMENT 
  59 :       AJAC=ABS(AJAC) 
  60 : * 
  61 : *-----   MATRICE ( B ) = UJ * BKSIP 
  62 : * 
  63 :       DO I=1,3 
  64 :          DO J=1,NBN 
  65 :          B(I,J)=0. 
  66 :          END DO 
  67 :       END DO 
  68 :       DO K=1,3 
  69 :         DO J=1,3 
  70 :           DO I=1,NBN 
  71 :             B(J,I)=B(J,I)+UJ(J,K)*BKSIP(K,I) 
  72 :           END DO 
  73 :         END DO 
  74 :       END DO 
  75 : * 
  76 :       END 
  77 :  
  78 :