Télécharger grad.eso

Retour à la liste

Numérotation des lignes :

  1. C GRAD SOURCE GF238795 18/07/10 21:15:27 9877
  2.  
  3. C=======================================================================
  4. C= G R A D =
  5. C= ------- =
  6. C= =
  7. C= OPERATEUR CAST3M "GRAD" DE CALCUL DU GRADIENT D'UN CHAMP : =
  8. C= ---------------------------------------------------------- =
  9. C= GRA1 = 'GRAD' MODL1 CHP1 ou CHE1 ( CAR1 ) ; =
  10. C= =
  11. C= Cet operateur calcule le gradient d'un champ de temperatures ou =
  12. C= de deplacements. =
  13. C= =
  14. C= ARGUMENTS : =
  15. C= ----------- =
  16. C= MODL1 (MMODEL) Modele associe a la structure etudiee =
  17. C= CHP1 (CHPOINT) | Champ de deplacements (NOEUDS) =
  18. C= CHE1 (MCHAML) | ou de temperatures (ELEMENT) =
  19. C= CAR1 (MCHAML) Caracteristiques geometriques (facultatif) =
  20. C= Sous-type 'CARACTERISTIQUES' =
  21. C= =
  22. C= RESULTAT : =
  23. C= ---------- =
  24. C= GRA1 (MCHAML) Champ de gradients sur la structure =
  25. C= =
  26. C= SUO X.Z., le 16 novemebre 1986. =
  27. C= Christian LE BRETON et Denis ROBERT-MOUGIN, le 31 juillet 1989. =
  28. C= Modifications aux nouvelles normes par I. MONNIER, le 28 mai 1990. =
  29. C=======================================================================
  30.  
  31. SUBROUTINE GRAD
  32.  
  33. IMPLICIT INTEGER(I-N)
  34. IMPLICIT REAL*8 (A-H,O-Z)
  35. INTEGER ZERO
  36.  
  37. -INC CCOPTIO
  38. -INC CCHAMP
  39.  
  40. C 1 - LECTURE DES ARGUMENTS DE L'OPERATEUR
  41. C ==========================================
  42. C 1.1 - Lecture OBLIGATOIRE du modele (IPMODL)
  43. C =====
  44. CALL LIROBJ('MMODEL',IPMODL,1,IRET)
  45. IF (IERR.NE.0) RETURN
  46. C =====
  47. C 1.2 - Syntaxe 1 : lecture d'un CHPOINT de depl. ou T (IPCHP1)
  48. C Transformation du CHPOINT en MCHAML aux noeuds (IPCHE2)
  49. C =====
  50. ZERO=0
  51. IRET1=0
  52. CALL LIROBJ('CHPOINT',IPCHP1,0,IRET1)
  53. IF (IERR.NE.0) RETURN
  54. IRET2=0
  55. IF (IRET1.NE.0) THEN
  56. IPCHL1=IPCHP1
  57. CALL CHAME1(0,IPMODL,IPCHP1,' ',IPCHE2,1)
  58. IF (IERR.NE.0) RETURN
  59. C =====
  60. C 1.3 - Syntaxe 3 : lecture d'un MCHAML de depl. ou T (IPCHE2)
  61. C ERREUR si le mode de calcul est "generalise"
  62. C =====
  63. ELSE
  64. IF (IFOUR.EQ.-3.OR.
  65. . (IFOUR.GE.7.AND.IFOUR.LE.11).OR.IFOUR.EQ.14) THEN
  66. MOTERR(1:8)='CHPOINT '
  67. CALL ERREUR(37)
  68. RETURN
  69. ENDIF
  70. IPCHL1=0
  71. CALL LIROBJ('MCHAML',IPIN,1,IRET2)
  72. IF (IERR.NE.0) RETURN
  73. CALL REDUAF(IPIN,IPMODL,IPCHE2,0,IR,KER)
  74. IF (IR.NE.1) CALL ERREUR(KER)
  75. IF(IERR.NE.0) RETURN
  76. ENDIF
  77. C =====
  78. C 1.4 - Lecture FACULTATIVE des caracteristiques geometriques (IPCHE1)
  79. C =====
  80. IRET3=0
  81. CALL LIROBJ('MCHAML',IPCHE1,0,IRET3)
  82. IF (IERR.NE.0) RETURN
  83. IPCH1=0
  84. IF (IRET3.NE.0.AND.IRET2.NE.0) THEN
  85. CALL RNGCHA(IPCHE1,IPCHE2,'CARACTERISTIQUES',
  86. . 'DEPLACEMENTS',IPCHA1,IPCHA2)
  87. IF (IPCHA1.EQ.0.OR.IPCHA2.EQ.0) THEN
  88. CALL ERREUR(805)
  89. RETURN
  90. ENDIF
  91. IPCHE1=IPCHA1
  92. IPCHE2=IPCHA2
  93. CALL REDUAF(IPCHE1,IPMODL,IPCH1,0,IR,KER)
  94. IF (IR.NE.1) CALL ERREUR(KER)
  95. IF(IERR.NE.0) RETURN
  96. IPCHE1 = IPCH1
  97. ENDIF
  98.  
  99. C =====
  100. C 1.5 - Lecture FACULTATIVE du mot 'CONS'
  101. C =====
  102. ISCA=0
  103. CALL LIRMOT('CONS',1,ISCA,0)
  104.  
  105. IF (ISCA.EQ.0) THEN
  106. CALL GRAD1(IPMODL,ZERO,IPCHE2,IPCHE1,IPCHL1,IRET)
  107. IF(IERR.NE.0) RETURN
  108. ELSE
  109. CALL GRAD2(IPMODL,IPCHE2,IPCHE1,IPCHL1,IRET)
  110. IF(IERR.NE.0) RETURN
  111. ENDIF
  112.  
  113. IF (IRET.EQ.1) THEN
  114. IF (IRET1.NE.0) CALL DTCHAM(IPCHE2)
  115. CALL ECROBJ('MCHAML',IPCHL1)
  116. ENDIF
  117.  
  118. RETURN
  119. END
  120.  
  121.  
  122.  
  123.  

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