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.  
  44. -INC PPARAM
  45. -INC CCOPTIO
  46. -INC CCREEL
  47. -INC SMCHAML
  48. -INC SMELEME
  49. -INC SMINTE
  50. -INC SMCOORD
  51.  
  52. PARAMETER (X2Pi=6.283185307179586476925286766559D0)
  53. PARAMETER (X4Pi=12.566370614359172953850573533118D0)
  54.  
  55. DIMENSION XE(3,1)
  56.  
  57. * OPTION INDISPONIBLE
  58. IF (IDIM.NE.1) THEN
  59. CALL ERREUR(19)
  60. RETURN
  61. ENDIF
  62.  
  63. * ON RECUPERE LES CARACTERISTIQUES D'INTEGRATION DES FACES
  64. MINTE=IPINTE
  65. SEGACT,MINTE
  66. NBPGAU=POIGAU(/1)
  67.  
  68. * ON RECUPERE UN DES MAILLAGES ELEMENTAIRES DE L'ENVELOPPE
  69. MELEME=IPGEOM
  70. SEGACT,MELEME
  71. NBPTEL=NUM(/1)
  72. NEL=NUM(/2)
  73.  
  74. C*OF IF ((NBPGAU.NE.1).OR.(NBPTEL.NE.1)) THEN
  75. C*OF WRITE(6,*) 'ERREUR FATALE : INTCN0'
  76. C*OF GOTO 10
  77. C*OF ENDIF
  78.  
  79. * ON RECUPERE LES VALEURS NODALES DU COEFFICIENT
  80. MELVA1=IPVANO
  81. SEGACT,MELVA1
  82. NBPTE1=MELVA1.VELCHE(/1)
  83. NEL1=MELVA1.VELCHE(/2)
  84. C*OF INMIN=MIN(NBPTE1,1)
  85.  
  86. * INITIALISATION DU MELVAL QUI CONTIENDRA LES VALEURS EQUIVALENTES
  87. N1PTEL=NBPTEL
  88. N1EL=NEL
  89. N2PTEL=0
  90. N2EL=0
  91. SEGINI,MELVAL
  92.  
  93. * BOUCLE SUR LES ELEMENTS
  94. DO iElt=1,NEL
  95.  
  96. * COORDONNEES DES NOEUDS DE L'ELEMENT DANS LE REPERE GLOBAL
  97. CALL DOXE(XCOOR,IDIM,NBPTEL,NUM,iElt,XE)
  98.  
  99. * CAS DES ELEMENTS AXISYMETRIQUES et spheriques
  100. IF (IFOUR.GE.12.AND.IFOUR.LE.14) THEN
  101. A11=X2Pi*XE(1,1)
  102. ELSE IF (IFOUR.EQ.15) THEN
  103. RR=XE(1,1)
  104. A11=X4Pi*RR*RR
  105. ELSE
  106. A11=1.
  107. ENDIF
  108. C*OF IEMIN=MIN(NEL1,iElt)
  109. C*OF VELCHE(1,iElt)=MELVA1.VELCHE(INMIN,IEMIN)*A11
  110. VELCHE(1,iElt)=MELVA1.VELCHE(1,iElt)*A11
  111. ENDDO
  112.  
  113. SEGDES,MELVAL,MELVA1
  114. 10 SEGDES,MELEME,MINTE
  115. IPVAEQ=MELVAL
  116.  
  117. RETURN
  118. END
  119.  
  120.  
  121.  
  122.  
  123.  

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