Télécharger relrig.eso

Retour à la liste

Numérotation des lignes :

  1. C RELRIG SOURCE PV 16/11/17 22:01:20 9180
  2. SUBROUTINE RELRIG()
  3. IMPLICIT REAL*8 (A-H,O-Z)
  4. IMPLICIT INTEGER (I-N)
  5. C***********************************************************************
  6. C NOM : RELRIG
  7. C DESCRIPTION : Lecture des données et écriture des résultats
  8. C pour la création des contraintes associées à un
  9. C opérateur linéaire
  10. C
  11. C LANGAGE : ESOPE
  12. C AUTEUR : Stéphane GOUNAND (CEA/DRN/DMT/SEMT/LTMF)
  13. C mél : gounand@semt2.smts.cea.fr
  14. C***********************************************************************
  15. C APPELES :
  16. C APPELES (E/S) :
  17. C APPELE PAR :
  18. C***********************************************************************
  19. C SYNTAXE GIBIANE :
  20. C ENTREES :
  21. C SORTIES :
  22. C***********************************************************************
  23. C VERSION : v1, 26/06/03, version initiale
  24. C HISTORIQUE : v1, 26/06/03, création
  25. C HISTORIQUE :
  26. C HISTORIQUE :
  27. C***********************************************************************
  28. C Prière de PRENDRE LE TEMPS de compléter les commentaires
  29. C en cas de modification de ce sous-programme afin de faciliter
  30. C la maintenance !
  31. C***********************************************************************
  32. -INC CCOPTIO
  33. -INC SMCHPOI
  34. -INC SMRIGID
  35. POINTEUR MLIN.MRIGID
  36. POINTEUR MCON.MRIGID
  37. * Includes persos
  38. CBEGININCLUDE SMMATASS
  39. SEGMENT MATASS
  40. POINTEUR KJPOPA.MLENTI
  41. POINTEUR LINCPA.MLMOTS
  42. POINTEUR MINCPA.MINC
  43. POINTEUR KJPODA.MLENTI
  44. POINTEUR LINCDA.MLMOTS
  45. POINTEUR MINCDA.MINC
  46. POINTEUR PROFMA.PMORS
  47. POINTEUR VALMA.IZA
  48. ENDSEGMENT
  49. CENDINCLUDE SMMATASS
  50. *
  51. INTEGER IMPR,IRET
  52. *
  53. * Executable statements
  54. *
  55. IMPR=0
  56. IF (IMPR.GT.2) WRITE(IOIMP,*) 'Entrée dans relrig'
  57. *
  58. * Lecture de la matrice des contraintes
  59. *
  60. CALL LIROBJ('RIGIDITE',MLIN,0,IRET)
  61. * Sinon, on va essayer de lire un MATRIK
  62. IF (IRET.EQ.0) THEN
  63. CALL RELRI2
  64. ELSE
  65. * Lecture d'un chpoint éventuel précisant les blocages
  66. CALL LIROBJ('CHPOINT',MCHPOI,0,IRET)
  67. *
  68. * Assemblage de la matrice
  69. *
  70. CALL RELR10(MLIN,MATASS,IMPR,IRET)
  71. IF (IRET.NE.0) GOTO 9999
  72. C*
  73. C* Création de la matrice des contraintes avec deux multiplicateurs
  74. C*
  75. C CALL RELR20(MATASS,MCON,IMPR,IRET)
  76. C IF (IRET.NE.0) GOTO 9999
  77. *
  78. * Création de la matrice des contraintes avec un multiplicateur
  79. *
  80. CALL RELR30(MATASS,MCHPOI,MCON,MCHPO2,IMPR,IRET)
  81. IF (IRET.NE.0) GOTO 9999
  82. *
  83. * Ecriture de la matrice des contraintes
  84. *
  85. IF (MCHPO2.NE.0) THEN
  86. CALL ECROBJ('CHPOINT',MCHPO2)
  87. ENDIF
  88. CALL ECROBJ('RIGIDITE',MCON)
  89. ENDIF
  90. *
  91. * Normal termination
  92. *
  93. IRET=0
  94. RETURN
  95. *
  96. * Format handling
  97. *
  98. *
  99. * Error handling
  100. *
  101. 9999 CONTINUE
  102. IRET=1
  103. WRITE(IOIMP,*) 'An error was detected in subroutine relrig'
  104. RETURN
  105. *
  106. * End of subroutine RELRIG
  107. *
  108. END
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  

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