Télécharger intcn1.eso

Retour à la liste

Numérotation des lignes :

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

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