Télécharger masse.eso

Retour à la liste

Numérotation des lignes :

  1. C MASSE SOURCE BP208322 15/06/22 21:20:30 8543
  2. SUBROUTINE MASSE(ILUMP)
  3. *
  4. *_______________________________________________________________________
  5. *
  6. * operateurs masse et lump
  7. *
  8. *
  9. * ILUMP : =1 indique si il s'agit de l'opérateur LUMP
  10. *
  11. *
  12. *_______________________________________________________________________
  13. IMPLICIT INTEGER(I-N)
  14. IMPLICIT REAL*8(A-H,O-Z)
  15. *
  16. -INC CCOPTIO
  17. -INC SMCHAML
  18. -INC CCHAMP
  19. c
  20. CHARACTER*4 CMOT
  21.  
  22. C Lump d une matrice deja constituée ( formule initiale )
  23. CALL LIROBJ('RIGIDITE',IRIG,0,IR1)
  24. IF (IR1.NE.0) THEN
  25. LMOT=0
  26. CALL LIROBJ('LISTMOTS',LMOT,0,IR2)
  27. IF (IERR.NE.0) RETURN
  28. if (ir2.eq.0.and.ilump.eq.0) then
  29. call ecrobj('RIGIDITE',IRIG)
  30. goto 10
  31. endif
  32.  
  33. *
  34. CALL LUMPIN(IRIG,LMOT,ILUM)
  35. IF (IERR.NE.0) RETURN
  36. *
  37. CALL ECROBJ('RIGIDITE',ILUM)
  38. RETURN
  39. END IF
  40. c
  41. c lecture du modele
  42. c
  43. 10 IPCHE1=0
  44. c
  45. CALL LIROBJ('MMODEL',IPMODL,0,IRT1)
  46. c
  47. IF(IRT1.NE.0) THEN
  48. c
  49. c Cas de la matrice masse sur base element finis
  50. c lecture du 1 ier mchaml
  51. c
  52. CALL LIROBJ('MCHAML',IPCHE1,1,IRT2)
  53. c
  54. c si absent erreur
  55. c
  56. IF (IRT2.EQ.0) THEN
  57. MOTERR(1:16) ='CARACTERISTIQUE'
  58. CALL ERREUR(565)
  59. RETURN
  60. ENDIF
  61. call reduaf(ipche1,ipmodl,ipch,0,ir,ker)
  62. if (ir.ne.1) call erreur(ker)
  63. IF (IERR.NE.0) RETURN
  64. ipche1=ipch
  65. c
  66. c calcul de la masse
  67. c
  68.  
  69. CALL MASSE1 (IPMODL,IPCHE1,IPMASS,IRET,ILUMP)
  70. c
  71. c ecriture du resultat
  72. c
  73. CALL ECROBJ('RIGIDITE',IPMASS)
  74.  
  75. ELSE
  76. IF (ILUMP .EQ. 0) THEN
  77. c
  78. c cas des masses additionnelles
  79. c
  80. CALL LIRCHA(CMOT,0,IRT1)
  81. IF (IRT1.NE.0) THEN
  82. CALL REFUS
  83. CALL APPUI(2)
  84. ELSE
  85. c
  86. c cas de la masse sur base modale
  87. c
  88. CALL MASSMO
  89. ENDIF
  90. ELSE
  91. CALL ERREUR(26)
  92. ENDIF
  93. ENDIF
  94. c
  95. c
  96. c
  97. RETURN
  98. END
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  

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