Télécharger rigi.eso

Retour à la liste

Numérotation des lignes :

rigi
  1. C RIGI SOURCE JK148537 26/10/19 21:15:14 12674
  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. CALL ACTOBJ('MMODEL ',IPMODL,1)
  87. CALL MODOPE (IPMODL,IPMOD5,IR1,1)
  88. IF (IERR.NE.0) RETURN
  89. IPMODL = IPMOD5
  90.  
  91.  
  92. IPCHE1 = 0
  93. IPCHE2 = 0
  94. IPCHA1 = 0
  95. IPCHA2 = 0
  96. *
  97. * ON LIT UN 1-ER MCHAML
  98. *
  99. CALL LIROBJ('MCHAML ',IPCHA1,1,IRET1)
  100. IF (IERR.NE.0) RETURN
  101. CALL ACTOBJ('MCHAML ',IPCHA1,1)
  102. *
  103. * ON LIT EVENTUELLEMENT UN 2EME MCHAML
  104. *
  105. CALL LIROBJ('MCHAML',IPCHA2,0,IRET1)
  106. IF (IERR.NE.0) RETURN
  107. IF (IRET1 .NE. 0)CALL ACTOBJ('MCHAML ',IPCHA2,1)
  108. *
  109. * RANGEMENT DES POINTEURS : 'CARACTERISTIQUES' DANS IPCHE1
  110. * 'MATRICE DE HOOKE' DANS IPCHE2
  111. *
  112. CALL RNGCHA(IPCHA1,IPCHA2,'CARACTERISTIQUES','MATRICE DE HOOKE',
  113. & IPCHE1,IPCHE2)
  114. IF (IERR.NE.0) RETURN
  115.  
  116. IF (IPCHE2.NE.0) THEN
  117. IMAT = 2
  118. ELSE IF (IPCHE1.NE.0) THEN
  119. IMAT = 1
  120. ELSE
  121. IPCHE1=IPCHA1
  122. IMAT = 1
  123. ENDIF
  124.  
  125. * APPEL AU CALCUL DES RIGIDITES ELEMENTAIRES
  126.  
  127. C Les REDUAF sont faits dans RIGI1 sur la partie utile du MMODEL
  128.  
  129. CALL RIGI1(IPMODL,IPCHE1,IPCHE2,IMAT,IPOI6,IRET,noer)
  130.  
  131. IF (noer.EQ.195) then
  132. CALL ECRent(noer)
  133. call soucis(195)
  134. return
  135. endif
  136. IF (IRET.EQ.1) CALL ECROBJ('RIGIDITE',IPOI6)
  137. RETURN
  138. C
  139. C -------------------------------------------------*
  140. C **** RAIDEURS ADDITIONNELLES
  141. C -------------------------------------------------*
  142. 700 CONTINUE
  143. CALL LIRCHA(CMOT,0,IRET)
  144. IF (IERR.NE.0) RETURN
  145. IF (IRET.EQ.0) GOTO 800
  146.  
  147. CALL REFUS
  148. CALL APPUI(1)
  149. RETURN
  150. C
  151. C -------------------------------------------------*
  152. C **** ANALYSE MODALE
  153. C -------------------------------------------------*
  154. 800 CONTINUE
  155. CALL RIGIMO
  156. RETURN
  157.  
  158. END
  159.  
  160.  
  161.  
  162.  

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