Télécharger cmct.eso

Retour à la liste

Numérotation des lignes :

  1. C CMCT SOURCE CB215821 19/08/01 21:15:23 10279
  2. SUBROUTINE CMCT
  3. *_______________________________________________________________________
  4. c
  5. c chapeau de l'opérateur cmct
  6. c
  7. c récupération des objets
  8. c quelques tests sur les objets obtenus
  9. c
  10. C ri2 = CMCT mrigid ri1
  11. C mrigid rigidite ( sans dependences)
  12. C ri1 rigidite de dependence
  13.  
  14. *_______________________________________________________________________
  15. IMPLICIT INTEGER(I-N)
  16. IMPLICIT REAL*8(A-H,O-Z)
  17. *
  18. -INC CCOPTIO
  19. -INC SMRIGID
  20. -INC SMELEME
  21.  
  22. *
  23. * lecture d'un objet rigidité
  24. *
  25. CALL LIROBJ('RIGIDITE',IRIG,1,IRET)
  26. IF (IERR .NE.0) THEN
  27. RETURN
  28. ENDIF
  29.  
  30.  
  31. * lecture optionnelle du champ par point
  32. *
  33. CALL LIROBJ('CHPOINT ',ICHP,0,IRET)
  34. IF(IRET .EQ. 1) CALL ACTOBJ('CHPOINT ',ICHP,1)
  35. IF (IERR .NE.0) RETURN
  36.  
  37. if(iret.eq.0) then
  38. C lecture de la matrice a condenser
  39. CALL LIROBJ('RIGIDITE',IRI1,1,IRET)
  40. IF (IERR .NE.0) THEN
  41. RETURN
  42. ENDIF
  43.  
  44. * on verifie que la rigidité ne contient est bien
  45. * de type dependance
  46. mrigid = iri1
  47. SEGACT , MRIGID
  48. DO 1510 I=1,IRIGEL(/2)
  49. C write(6,*) ' irigel(/1) irigel(8,I)',irigel(/1),IRIGEL(8,I)
  50. if(irigel(/1).LT.8.or.IRIGEL(8,I).eq.0) then
  51. CALL ERREUR(888)
  52. SEGDES MRIGID
  53. RETURN
  54. ENDIF
  55. 1510 CONTINUE
  56. segdes mrigid
  57. call depen3(mrigid,ri6)
  58. call scnd2(irig,ri6,irig2)
  59. segsup ri6
  60.  
  61. else
  62. *
  63. * verification que la rigidité ne contient que des relation
  64. * c'est à dire n'est supportée que par des maillage de type22
  65. *
  66. MRIGID = IRIG
  67. SEGACT, MRIGID*NOMOD
  68. DO 100 I=1,IRIGEL(/2)
  69. MELEME = IRIGEL(1,I)
  70. SEGACT , MELEME*NOMOD
  71. IF (ITYPEL.NE.22) THEN
  72. CALL ERREUR(837)
  73. SEGDES MELEME
  74. SEGDES MRIGID
  75. RETURN
  76. ENDIF
  77. 100 CONTINUE
  78. *
  79. CALL CMCT1(ICHP,IRIG,IRIG2)
  80. *
  81. IF ( IERR .NE. 0 ) RETURN
  82.  
  83. endif
  84. *
  85. CALL ECROBJ('RIGIDITE',IRIG2)
  86.  
  87. END
  88.  
  89.  
  90.  

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