Télécharger intcn1.eso

Retour à la liste

Numérotation des lignes :

  1. C INTCN1 SOURCE PV 09/03/12 21:25:19 6325
  2. SUBROUTINE INTCN1(IPVANO,IPGEOM,IPINTE,IPVAEQ)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. ************************************************************************
  6. *
  7. * I N T C N 1
  8. * -----------
  9. *
  10. * FONCTION:
  11. * ---------
  12. * INTEGRATION NUMERIQUE DANS UN DOMAINE UNIDIMENSIONNEL DU PRODUIT:
  13. * COEF. * TRANSPOSEE( N )
  14. * COEF. : GRANDEUR PHYSIQUE REPRESENTEE PAR UN CHAMELEM
  15. * N : FONCTIONS DE FORME DE L'ELEMENT MAILLANT LE DOMAINE
  16. * CONSIDERE
  17. *
  18. * MODULES UTILISES:
  19. * -----------------
  20. *
  21. -INC CCOPTIO
  22. -INC CCREEL
  23. -INC SMCHAML
  24. -INC SMELEME
  25. -INC SMINTE
  26. -INC SMCOORD
  27. *
  28. * PARAMETRES: (E)=ENTREE (S)=SORTIE (+ = CONTENU DANS UN COMMUN)
  29. * -----------
  30. *
  31. * IPVANO (E) POINTEUR SUR UN SEGMENT MELVAL CONTENANT LES
  32. * VALEURS NODALES DE NOTRE COEFFICIENT
  33. * IPGEOM (E) POINTEUR SUR UN OBJET MAILLAGE ELEMENTAIRE
  34. * DU DOMAINE D'INTEGRATION
  35. * IPINTE (E) POINTEUR SUR UN SEGMENT MINTE CONTENANT LES
  36. * CARACTERISTIQUES D'INTEGRATION DES FACES
  37. * +IDIM (E) VOIR CCOPTIO
  38. * +XZERO (E) VOIR CCREEL
  39. * +XPETIT (E) VOIR CCREEL
  40. * IPVAEQ (S) POINTEUR SUR UN SEGMENT MELVAL CONTENANT LES VALEURS
  41. * NODALES EQUIVALENTES (APRES INTEGRATION NUMERIQUE)
  42. *
  43. * VARIABLES:
  44. * ----------
  45. *
  46. * XE(3,NBPTEL) = CORDONNEES DES ELEMENTS DANS LE REPERE GLOBAL
  47. * FACE = LONGUEUR DE DE L'ELEMENT CONSIDERE
  48. * A = TABLEAU DE TRAVAIL
  49. *
  50. SEGMENT,MMAT1
  51. REAL*8 XE(3,NBPTEL),A(NBPTEL,NBPTEL),SHP(6,NBPTEL)
  52. ENDSEGMENT
  53. *
  54. * CONSTANTES:
  55. * -----------
  56. *
  57. PARAMETER ( O1=1.D0 )
  58. PARAMETER ( O2=2.D0 )
  59. *
  60. * REMARQUES:
  61. * ----------
  62. *
  63. * L'UTILISATION DE CE S-P PRESUPPOSE UN PRE ET POST-TRAITEMENT
  64. * DES SEGMENTS MELVAL PASSES EN TANT QUE PARAMETRES
  65. *
  66. * AUTEUR, DATE DE CREATION:
  67. * -------------------------
  68. *
  69. * DENIS ROBERT,LE 15 AVRIL 1988.
  70. *
  71. * LANGAGE:
  72. * --------
  73. *
  74. * ESOPE + FORTRAN77
  75. *
  76. ************************************************************************
  77. IF (IDIM.EQ.3) THEN
  78. *
  79. * OPTION INDISPONIBLE
  80. *
  81. CALL ERREUR(19)
  82. RETURN
  83. ENDIF
  84. *
  85. * ON RECUPERE LES VALEURS NODALES DU COEFFICIENT
  86. *
  87. MELVA1=IPVANO
  88. SEGACT,MELVA1
  89. NBPTE1=MELVA1.VELCHE(/1)
  90. NEL1=MELVA1.VELCHE(/2)
  91. *
  92. * ON RECUPERE LES CARACTERISTIQUES D'INTEGRATION DES FACES
  93. *
  94. MINTE=IPINTE
  95. SEGACT,MINTE
  96. NBPGAU=POIGAU(/1)
  97. *
  98. * ON RECUPERE UN DES MAILLAGES ELEMENTAIRES DE L'ENVELOPPE
  99. *
  100. MELEME=IPGEOM
  101. SEGACT,MELEME
  102. NBPTEL=NUM(/1)
  103. NEL=NUM(/2)
  104. *
  105. * INITIALISATION DU MELVAL QUI CONTIENDRA LES VALEURS EQUIVALENTES
  106. *
  107. N1PTEL=NBPTEL
  108. N1EL=NEL
  109. N2PTEL=0
  110. N2EL=0
  111. SEGINI,MELVAL
  112. IPVAEQ=MELVAL
  113. SEGINI,MMAT1
  114. *
  115. * BOUCLE SUR LES ELEMENTS
  116. *
  117. DO 10 IEL=1,NEL
  118. *
  119. * ON CHERCHE LES COORDONNEES DES ELEMENTS DANS LE REPERE GLOBAL
  120. *
  121. CALL DOXE(XCOOR,IDIM,NBPTEL,NUM,IEL,XE)
  122. *
  123. * CALCUL DE LA LONGUEUR DE LA FACE CONSIDEREE
  124. *
  125. C DLX=XZERO
  126. C DLY=XZERO
  127. C DO 20 IGAU=1,NBPGAU
  128. C DLONX=XZERO
  129. C DLONY=XZERO
  130. C DO 30 INOE=1,NBPTEL
  131. C DLONX=DLONX+SHPTOT(2,INOE,IGAU)*XE(1,INOE)
  132. C DLONY=DLONY+SHPTOT(2,INOE,IGAU)*XE(2,INOE)
  133. C 30 CONTINUE
  134. C* END DO
  135. CC DLX=DLX+DLONX*POIGAU(IGAU)
  136. CC DLY=DLY+DLONY*POIGAU(IGAU)
  137. C DLX= DLONX
  138. C DLY= DLONY
  139. C 20 CONTINUE
  140. C* END DO
  141. C FACE=SQRT(DLX**2+DLY**2)
  142. *
  143. * BOUCLES SUR LES POINTS DE GAUSS ET SUR LES NOEUDS
  144. *
  145. CALL ZERO(A,NBPTEL,NBPTEL)
  146. DO 40 IGAU=1,NBPGAU
  147. DLX=XZERO
  148. DLY=XZERO
  149. DLONX=XZERO
  150. DLONY=XZERO
  151. DO 30 INOE=1,NBPTEL
  152. DLONX=DLONX+SHPTOT(2,INOE,IGAU)*XE(1,INOE)
  153. DLONY=DLONY+SHPTOT(2,INOE,IGAU)*XE(2,INOE)
  154. 30 CONTINUE
  155. * END DO
  156. C DLX=DLX+DLONX*POIGAU(IGAU)
  157. C DLY=DLY+DLONY*POIGAU(IGAU)
  158. DLX= DLONX
  159. DLY= DLONY
  160.  
  161. * END DO
  162. FACE=SQRT(DLX**2+DLY**2)
  163. AXIS=O1
  164. IF (IFOMOD.EQ.0) THEN
  165. *
  166. * CAS DES ELEMENTS AXISYMETRIQUES
  167. *
  168. DO 41 NP=1,NBPTEL
  169. SHP(1,NP)=SHPTOT(1,NP,IGAU)
  170. SHP(2,NP)=SHPTOT(2,NP,IGAU)
  171. SHP(3,NP)=SHPTOT(3,NP,IGAU)
  172. 41 CONTINUE
  173. * END DO
  174. CALL DISTRR(XE,SHP,NBPTEL,RR)
  175. AXIS=O2*XPI*RR
  176. ENDIF
  177. DO 50 INOE=1,NBPTEL
  178. DO 60 INO2=1,NBPTEL
  179. A(INOE,INO2)=A(INOE,INO2)+SHPTOT(1,INOE,IGAU)*
  180. 1 SHPTOT(1,INO2,IGAU)*POIGAU(IGAU)*AXIS*FACE
  181. 60 CONTINUE
  182. * END DO
  183. 50 CONTINUE
  184. * END DO
  185. 40 CONTINUE
  186. * END DO
  187. IEMIN=MIN(NEL1,IEL)
  188. DO 70 INOE=1,NBPTEL
  189. DO 80 INO2=1,NBPTEL
  190. INMIN=MIN(NBPTE1,INO2)
  191. VELCHE(INOE,IEL)=VELCHE(INOE,IEL)+MELVA1.VELCHE(INMIN,IEMIN)*
  192. 1 A(INO2,INOE)
  193. 80 CONTINUE
  194. * END DO
  195. 70 CONTINUE
  196. * END DO
  197. 10 CONTINUE
  198. * END DO
  199. SEGSUP,MMAT1
  200. SEGDES,MELEME,MELVA1,MINTE
  201. SEGDES,MELVAL
  202. *
  203. END
  204.  
  205.  
  206.  
  207.  
  208.  
  209.  

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