Télécharger intcn0.eso

Retour à la liste

Numérotation des lignes :

  1. C INTCN0 SOURCE PV 09/03/12 21:25:18 6325
  2.  
  3. C=======================================================================
  4. C= I N T C N 0 =
  5. * ----------- =
  6. *
  7. * FONCTION:
  8. * ---------
  9. * INTEGRATION NUMERIQUE DANS UN DOMAINE UNIDIMENSIONNEL DU PRODUIT:
  10. * COEF. * TRANSPOSEE( N )
  11. * COEF. : GRANDEUR PHYSIQUE REPRESENTEE PAR UN CHAMELEM
  12. * N : FONCTIONS DE FORME DE L'ELEMENT MAILLANT LE DOMAINE
  13. * CONSIDERE
  14. *
  15. * PARAMETRES: (E)=ENTREE (S)=SORTIE (+ = CONTENU DANS UN COMMUN)
  16. * -----------
  17. * IPVANO (E) POINTEUR SUR UN SEGMENT MELVAL CONTENANT LES
  18. * VALEURS NODALES DE NOTRE COEFFICIENT
  19. * IPGEOM (E) POINTEUR SUR UN OBJET MAILLAGE ELEMENTAIRE
  20. * DU DOMAINE D'INTEGRATION
  21. * IPINTE (E) POINTEUR SUR UN SEGMENT MINTE CONTENANT LES
  22. * CARACTERISTIQUES D'INTEGRATION DES FACES
  23. * IPVAEQ (S) POINTEUR SUR UN SEGMENT MELVAL CONTENANT LES VALEURS
  24. * NODALES EQUIVALENTES (APRES INTEGRATION NUMERIQUE)
  25. *
  26. * VARIABLES:
  27. * ----------
  28. * XE(3,NBPTEL) = CORDONNEES DES ELEMENTS DANS LE REPERE GLOBAL
  29. * FACE = LONGUEUR DE DE L'ELEMENT CONSIDERE
  30. * A = TABLEAU DE TRAVAIL
  31. *
  32. * REMARQUES:
  33. * ----------
  34. * L'UTILISATION DE CE S-P PRESUPPOSE UN PRE ET POST-TRAITEMENT
  35. * DES SEGMENTS MELVAL PASSES EN TANT QUE PARAMETRES
  36. C=======================================================================
  37.  
  38. SUBROUTINE INTCN0 (IPVANO,IPGEOM,IPINTE,IPVAEQ)
  39.  
  40. IMPLICIT INTEGER(I-N)
  41. IMPLICIT REAL*8 (A-H,O-Z)
  42.  
  43. -INC CCOPTIO
  44. -INC CCREEL
  45. -INC SMCHAML
  46. -INC SMELEME
  47. -INC SMINTE
  48. -INC SMCOORD
  49.  
  50. PARAMETER (X2Pi=6.283185307179586476925286766559D0)
  51. PARAMETER (X4Pi=12.566370614359172953850573533118D0)
  52.  
  53. DIMENSION XE(3,1)
  54.  
  55. * OPTION INDISPONIBLE
  56. IF (IDIM.NE.1) THEN
  57. CALL ERREUR(19)
  58. RETURN
  59. ENDIF
  60.  
  61. * ON RECUPERE LES CARACTERISTIQUES D'INTEGRATION DES FACES
  62. MINTE=IPINTE
  63. SEGACT,MINTE
  64. NBPGAU=POIGAU(/1)
  65.  
  66. * ON RECUPERE UN DES MAILLAGES ELEMENTAIRES DE L'ENVELOPPE
  67. MELEME=IPGEOM
  68. SEGACT,MELEME
  69. NBPTEL=NUM(/1)
  70. NEL=NUM(/2)
  71.  
  72. C*OF IF ((NBPGAU.NE.1).OR.(NBPTEL.NE.1)) THEN
  73. C*OF WRITE(6,*) 'ERREUR FATALE : INTCN0'
  74. C*OF GOTO 10
  75. C*OF ENDIF
  76.  
  77. * ON RECUPERE LES VALEURS NODALES DU COEFFICIENT
  78. MELVA1=IPVANO
  79. SEGACT,MELVA1
  80. NBPTE1=MELVA1.VELCHE(/1)
  81. NEL1=MELVA1.VELCHE(/2)
  82. C*OF INMIN=MIN(NBPTE1,1)
  83.  
  84. * INITIALISATION DU MELVAL QUI CONTIENDRA LES VALEURS EQUIVALENTES
  85. N1PTEL=NBPTEL
  86. N1EL=NEL
  87. N2PTEL=0
  88. N2EL=0
  89. SEGINI,MELVAL
  90.  
  91. * BOUCLE SUR LES ELEMENTS
  92. DO iElt=1,NEL
  93.  
  94. * COORDONNEES DES NOEUDS DE L'ELEMENT DANS LE REPERE GLOBAL
  95. CALL DOXE(XCOOR,IDIM,NBPTEL,NUM,iElt,XE)
  96.  
  97. * CAS DES ELEMENTS AXISYMETRIQUES et spheriques
  98. IF (IFOUR.GE.12.AND.IFOUR.LE.14) THEN
  99. A11=X2Pi*XE(1,1)
  100. ELSE IF (IFOUR.EQ.15) THEN
  101. RR=XE(1,1)
  102. A11=X4Pi*RR*RR
  103. ELSE
  104. A11=1.
  105. ENDIF
  106. C*OF IEMIN=MIN(NEL1,iElt)
  107. C*OF VELCHE(1,iElt)=MELVA1.VELCHE(INMIN,IEMIN)*A11
  108. VELCHE(1,iElt)=MELVA1.VELCHE(1,iElt)*A11
  109. ENDDO
  110.  
  111. SEGDES,MELVAL,MELVA1
  112. 10 SEGDES,MELEME,MINTE
  113. IPVAEQ=MELVAL
  114.  
  115. RETURN
  116. END
  117.  
  118.  
  119.  
  120.  
  121.  

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