Télécharger grad.eso

Retour à la liste

Numérotation des lignes :

  1. C GRAD SOURCE CB215821 16/12/05 21:39:31 9237
  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.  
  36. -INC CCOPTIO
  37. -INC CCHAMP
  38.  
  39. C 1 - LECTURE DES ARGUMENTS DE L'OPERATEUR
  40. C ==========================================
  41. C 1.1 - Lecture OBLIGATOIRE du modele (IPMODL)
  42. C =====
  43. CALL LIROBJ('MMODEL',IPMODL,1,IRET)
  44. IF (IERR.NE.0) RETURN
  45. C =====
  46. C 1.2 - Syntaxe 1 : lecture d'un CHPOINT de depl. ou T (IPCHP1)
  47. C Transformation du CHPOINT en MCHAML aux noeuds (IPCHE2)
  48. C =====
  49. IRET1=0
  50. CALL LIROBJ('CHPOINT',IPCHP1,0,IRET1)
  51. IF (IERR.NE.0) RETURN
  52. IRET2=0
  53. IF (IRET1.NE.0) THEN
  54. IPCHL1=IPCHP1
  55. CALL CHAME1(0,IPMODL,IPCHP1,' ',IPCHE2,1)
  56. IF (IERR.NE.0) RETURN
  57. C =====
  58. C 1.3 - Syntaxe 3 : lecture d'un MCHAML de depl. ou T (IPCHE2)
  59. C ERREUR si le mode de calcul est "generalise"
  60. C =====
  61. ELSE
  62. IF (IFOUR.EQ.-3.OR.
  63. . (IFOUR.GE.7.AND.IFOUR.LE.11).OR.IFOUR.EQ.14) THEN
  64. MOTERR(1:8)='CHPOINT '
  65. CALL ERREUR(37)
  66. RETURN
  67. ENDIF
  68. IPCHL1=0
  69. CALL LIROBJ('MCHAML',IPIN,1,IRET2)
  70. IF (IERR.NE.0) RETURN
  71. CALL REDUAF(IPIN,IPMODL,IPCHE2,0,IR,KER)
  72. IF (IR.NE.1) CALL ERREUR(KER)
  73. IF(IERR.NE.0) RETURN
  74. ENDIF
  75. C =====
  76. C 1.4 - Lecture FACULTATIVE des caracteristiques geometriques (IPCHE1)
  77. C =====
  78. IRET3=0
  79. CALL LIROBJ('MCHAML',IPCHE1,0,IRET3)
  80. IF (IERR.NE.0) RETURN
  81. IPCH1=0
  82. IF (IRET3.NE.0.AND.IRET2.NE.0) THEN
  83. CALL RNGCHA(IPCHE1,IPCHE2,'CARACTERISTIQUES',
  84. . 'DEPLACEMENTS',IPCHA1,IPCHA2)
  85. IF (IPCHA1.EQ.0.OR.IPCHA2.EQ.0) THEN
  86. CALL ERREUR(805)
  87. RETURN
  88. ENDIF
  89. IPCHE1=IPCHA1
  90. IPCHE2=IPCHA2
  91. CALL REDUAF(IPCHE1,IPMODL,IPCH1,0,IR,KER)
  92. IF (IR.NE.1) CALL ERREUR(KER)
  93. IF(IERR.NE.0) RETURN
  94. IPCHE1 = IPCH1
  95. ENDIF
  96.  
  97. C =====
  98. C 1.5 - Lecture FACULTATIVE du mot 'CONS'
  99. C =====
  100. ISCA=0
  101. CALL LIRMOT('CONS',1,ISCA,0)
  102.  
  103. IF (ISCA.EQ.0) THEN
  104. CALL GRAD1(IPMODL,IPCHE2,IPCHE1,IPCHL1,IRET)
  105. IF(IERR.NE.0) RETURN
  106. ELSE
  107. CALL GRAD2(IPMODL,IPCHE2,IPCHE1,IPCHL1,IRET)
  108. IF(IERR.NE.0) RETURN
  109. ENDIF
  110.  
  111. IF (IRET.EQ.1) THEN
  112. IF (IRET1.NE.0) CALL DTCHAM(IPCHE2)
  113. CALL ECROBJ('MCHAML',IPCHL1)
  114. ENDIF
  115.  
  116. RETURN
  117. END
  118.  
  119.  
  120.  

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