Esope source files

   1 : C MCONDT    SOURCE    FANDEUR   11/07/19    21:16:37     7042           
   2 : 
   3 : * ******************************************************************* 
   4 : *     CALCUL DE LA MATRICE DE CONDUCTIVITE D'UNE COQUE(COQ2 ET COQ3)* 
   5 : *     PAR UNE METHODE SEMI-ANALYTIQUE                               * 
   6 : ********************************************************************* 
   7 : * 
   8 : * PARAMETRES:  (E)=ENTREE  (S)=SORTIE  (+ = CONTENU DANS UN COMMUN ) 
   9 : *------------ 
  10 : * 
  11 : *      CEL1   (E)  MATRICE BTKB INTERGREE SUR LA SURFACE MOYENNE 
  12 : *      CEL2   (E)  MATRICE HTKH  INTEGREE SUR LA SURFACE MOYENNE 
  13 : *      NBNN   (E)  NOMBRE DE NOEUDS PAR ELEMENT 
  14 : *      EP     (E)  EPAISSEUR DE LA COQUE 
  15 : *      CEL    (S)  MATRICE DE CONDUCTIBILITE TOTALE 
  16 : * 
  17 : * COMMENTAIRE :ON FORME LES SOUS-MATRICES DE LA MATRICE DE CONDUCTIVITE 
  18 : * PAR DES COMBINAISONS LINEAIRES DE LA MATRICE BTKB(CEL1) ET HTKT(CEL2). 
  19 : * LES COEFFICIENTS CORRESPONDENT AUX INTEGRALES  SUR L'EPAISSEUR. 
  20 : * 
  21 : *  AUTEUR,DATE DE CREATION: 
  22 : *  ------------------------ 
  23 : *  PARVIZ DOWLATYARI MAI 1990 
  24 : ************************************************************************ 
  25 : 
  26 :       SUBROUTINE MCONDT(CEL1,CEL2,NBNN,EP,CEL)
  27 : 
  28 :       IMPLICIT INTEGER(I-N)
  29 :       IMPLICIT REAL*8 (A-H,O-Z)
  30 : 
  31 :       DIMENSION CEL1(NBNN,*),CEL2(NBNN,*),CEL(3*NBNN,*)
  32 : 
  33 :       DATA UN,DEUX,TROIS,SEPT/1.D0,2.D0,3.D0,7.D0/
  34 :       DATA HUIT,QUINZE,SEIZE,TRENTE/8.D0,15.D0,16.D0,30.D0/
  35 : 
  36 :       C1  = (DEUX*EP)/QUINZE
  37 :       C2  = SEPT/(TROIS*EP)
  38 :       C3  = EP/QUINZE
  39 :       C4  = -HUIT/(TROIS*EP)
  40 :       C5  = -EP/TRENTE
  41 :       C6  = UN/(TROIS*EP)
  42 :       C7  = (HUIT*EP)/QUINZE
  43 :       C8  = SEIZE/(TROIS*EP)
  44 :       C9  = EP/QUINZE
  45 :       C10 = -HUIT/(TROIS*EP)
  46 :       C11 = (DEUX*EP)/QUINZE
  47 :       C12 = SEPT/(TROIS*EP)
  48 : 
  49 :       NBN2 = 2 * NBNN
  50 : 
  51 :       DO J = 1, NBNN
  52 :         J1 = J + NBNN
  53 :         J2 = J + NBN2
  54 : 
  55 :         DO I = 1, NBNN
  56 :           I1 = I + NBNN
  57 :           I2 = I + NBN2
  58 : 
  59 :           r_z1 = CEL1(I,J)
  60 :           r_z2 = CEL2(I,J)
  61 : 
  62 : *    SOUS-MATRICE K11 = C1*CEL1+C2*CEL2
  63 :           CEL(I ,J ) = C1*r_z1 + C2*r_z2
  64 : *    SOUS-MATRICE K21 = C3*CEL1+C4*CEL2
  65 :           CEL(I1,J ) = C3*r_z1 + C4*r_z2
  66 : *    SOUS-MATRICE K12 = K21
  67 :           CEL(J ,I1) = CEL(I1,J )
  68 : *    SOUS-MATRICE K31 = C5*CEL1+C6*CEL2
  69 :           CEL(I2,J ) = C5*r_z1 + C6*r_z2
  70 : *    SOUS-MATRICE K13 = K31
  71 :           CEL(J ,I2) = CEL(I2,J )
  72 : *    SOUS-MATRICE K22 = C7*CEL1+C8*CEL2
  73 :           CEL(I1,J1) = C7*r_z1 + C8*r_z2
  74 : *    SOUS-MATRICE K32 = C9*CEL1+C10*CEL2
  75 :           CEL(I2,J1) = C9*r_z1 + C10*r_z2
  76 : *    SOUS-MATRICE K23 = K32
  77 :           CEL(J1,I2) = CEL(I2,J1)
  78 : *    SOUS-MATRICE K33 = C11*CEL1+C12*CEL2
  79 :           CEL(I2,J2) = C11*r_z1 + C12*r_z2
  80 : 
  81 :         ENDDO
  82 : 
  83 :       ENDDO
  84 : 
  85 :       RETURN
  86 :       END
  87 : 
  88 :