Télécharger rigi.eso

Retour à la liste

Numérotation des lignes :

rigi
  1. C RIGI SOURCE OF166741 24/10/07 21:15:46 12016
  2.  
  3. *---------------------------------------------------------------------*
  4. * __________________________ *
  5. * | | *
  6. * | OPERATEUR RIGIDITE | *
  7. * |________________________| *
  8. * *
  9. *---------------------------------------------------------------------*
  10. * *
  11. * NOUVELLE SYNTAXE : *
  12. * __________________ *
  13. * *
  14. * *
  15. * RIG1 = RIGI MOD1 | CAR1 ; *
  16. * | HOO1 ( CAR1 ) ; *
  17. * *
  18. * ENTREES : *
  19. * *
  20. * MOD1 OBJET MODELE *
  21. * CAR1 CHAMELEM DE CARACTERISTIQUES *
  22. * GEOMETRIQUES ET MATERIELLES *
  23. * *
  24. * EVENTUELLEMENT, *
  25. * *
  26. * HOO1 CHAMELEM REPRESENTANT UNE MATRICE DE HOOKE *
  27. * *
  28. * *
  29. * CHRISTIAN LE BRETON , DECEMBRE 1989 *
  30. * CONVERSION AU NOUVEAU CHAMELEM I. MONNIER AVRIL 90 *
  31. * *
  32. * *
  33. * ANALYSE MODALE : *
  34. * ________________ *
  35. * *
  36. * SYNTAXES POSSIBLES : RIG2 = RIGI M ; *
  37. * RIG3 = RIGI P S ; *
  38. * RIG4 = RIGI M P S ; *
  39. * RIG5 = RIGI B ; *
  40. * RIG6 = RIGI A ; *
  41. * *
  42. * M : OBJET SOLUTION DE SOUS TYPE MODE *
  43. * P : OBJET SOLUTION DE SOUS TYPE SOLU STAT *
  44. * S : OBJET STRUCTURE *
  45. * B : OBJET BASE MODALE *
  46. * A : OBJET ATTACHE *
  47. * *
  48. * RIG2 : RIGIDITE DUE AUX MODES SEULS ( KN ) *
  49. * RIG3 : RIGIDITE DUE AU COUPLAGE DES SOLUTIONS STATIQUES *
  50. * SUR UNE STRUCTURE . *
  51. * RIG4 : RIGIDITE DUE AU COUPLAGE DES SOLUTIONS STATIQUES *
  52. * ET DES MODES . *
  53. * RIG5 : CALCUL AUTOMATIQUE DE RIG2, RIG3 ET RIG4. *
  54. * RIG6 : RIGIDITE DUE AUX LIAISONS MIXTES *
  55. * *
  56. *---------------------------------------------------------------------*
  57.  
  58. SUBROUTINE RIGI
  59.  
  60. IMPLICIT INTEGER(I-N)
  61. IMPLICIT REAL*8(A-H,O-Z)
  62.  
  63. -INC PPARAM
  64. -INC CCOPTIO
  65. -INC SMCOORD
  66.  
  67. character*4 mcle(2)
  68. data mcle/'NOER','DIVE'/
  69. CHARACTER*4 CMOT
  70.  
  71. C On a besoin du MCOORD plus loin
  72. SEGACT,MCOORD
  73.  
  74. NOER = 0
  75. * lecture eventuelle de 'NOER'
  76. call lirmot(mcle,2,NOER,0)
  77. if (ierr.ne.0) returN
  78.  
  79. * -------------------------------------------------*
  80. * ON RECUPERE EVENTUELLEMENT UN OBJET MODELE
  81. * -------------------------------------------------*
  82. IPMODL = 0
  83. CALL LIROBJ('MMODEL ',IPMODL,0,IRET)
  84. if (ierr.ne.0) return
  85. IF (IRET.EQ.0) GO TO 700
  86.  
  87. CALL ACTOBJ('MMODEL ',IPMODL,1)
  88.  
  89. IPCHE1 = 0
  90. IPCHE2 = 0
  91. IPCHA1 = 0
  92. IPCHA2 = 0
  93. *
  94. * ON LIT UN 1-ER MCHAML
  95. *
  96. CALL LIROBJ('MCHAML ',IPCHA1,1,IRET1)
  97. IF (IERR.NE.0) RETURN
  98. CALL ACTOBJ('MCHAML ',IPCHA1,1)
  99. *
  100. * ON LIT EVENTUELLEMENT UN 2EME MCHAML
  101. *
  102. CALL LIROBJ('MCHAML',IPCHA2,0,IRET1)
  103. IF (IERR.NE.0) RETURN
  104. IF (IRET1 .NE. 0)CALL ACTOBJ('MCHAML ',IPCHA2,1)
  105. *
  106. * RANGEMENT DES POINTEURS : 'CARACTERISTIQUES' DANS IPCHE1
  107. * 'MATRICE DE HOOKE' DANS IPCHE2
  108. *
  109. CALL RNGCHA(IPCHA1,IPCHA2,'CARACTERISTIQUES','MATRICE DE HOOKE',
  110. & IPCHE1,IPCHE2)
  111. IF (IERR.NE.0) RETURN
  112.  
  113. IF (IPCHE2.NE.0) THEN
  114. IMAT = 2
  115. ELSE IF (IPCHE1.NE.0) THEN
  116. IMAT = 1
  117. ELSE
  118. IPCHE1=IPCHA1
  119. IMAT = 1
  120. ENDIF
  121.  
  122. * APPEL AU CALCUL DES RIGIDITES ELEMENTAIRES
  123.  
  124. C Les REDUAF sont faits dans RIGI1 sur la partie utile du MMODEL
  125.  
  126. CALL RIGI1(IPMODL,IPCHE1,IPCHE2,IMAT,IPOI6,IRET,noer)
  127.  
  128. IF (noer.EQ.195) then
  129. CALL ECRent(noer)
  130. call soucis(195)
  131. return
  132. endif
  133. IF (IRET.EQ.1) CALL ECROBJ('RIGIDITE',IPOI6)
  134. RETURN
  135. C
  136. C -------------------------------------------------*
  137. C **** RAIDEURS ADDITIONNELLES
  138. C -------------------------------------------------*
  139. 700 CONTINUE
  140. CALL LIRCHA(CMOT,0,IRET)
  141. IF (IERR.NE.0) RETURN
  142. IF (IRET.EQ.0) GOTO 800
  143.  
  144. CALL REFUS
  145. CALL APPUI(1)
  146. RETURN
  147. C
  148. C -------------------------------------------------*
  149. C **** ANALYSE MODALE
  150. C -------------------------------------------------*
  151. 800 CONTINUE
  152. CALL RIGIMO
  153. RETURN
  154.  
  155. END
  156.  
  157.  
  158.  

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