Télécharger mcondt.eso

Retour à la liste

Numérotation des lignes :

  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.  
  89.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales