Télécharger coriol.eso

Retour à la liste

Numérotation des lignes :

  1. C CORIOL SOURCE CB215821 19/08/01 21:15:33 10279
  2. SUBROUTINE CORIOL
  3.  
  4. C***********************************************************************
  5. C *
  6. C Première option : *
  7. C *
  8. C FABRICATION DE LA MATRICE DE COUPLAGE DE CORIOLIS DANS *
  9. C LE REPERE TOURNANT (ELEM BARR, POUT, TUYAU, COQUE, MASSIF 3D) *
  10. C *
  11. C Syntaxe : *
  12. C -------- *
  13. C *
  14. C RIG1 = CORI MODL1 MAT1 VEC1 *
  15. C *
  16. C RIG1 : matrice d'amortissement construite (TYPE rigidité, SOUS-TYPE *
  17. C amortissement) *
  18. C MODL1: objet MMODEL, modèle *
  19. C MAT1 : objet MCHAML, caractéristiques matériau *
  20. C VEC1 : objet POINT, vecteur rotation *
  21. C *
  22. C *
  23. C Didier COMBESCURE mars 2003 *
  24. C *
  25. C***********************************************************************
  26.  
  27. IMPLICIT INTEGER(I-N)
  28. IMPLICIT REAL*8 (A-H,O-Z)
  29.  
  30. -INC CCOPTIO
  31. c
  32. CHARACTER*4 LISMO1(1),LISMO2(1)
  33. DATA LISMO1 / 'HARM' /
  34. DATA LISMO2 / 'FLAM' /
  35.  
  36. IPMODL = 0
  37. IPCHE1 = 0
  38. IPROTA = 0
  39.  
  40. NUMLI1 = 0
  41. NUMLI2 = 0
  42.  
  43. IPRIG = 0
  44. C
  45. C lecture du modèle
  46. C
  47. CALL LIROBJ('MMODEL ',IPMODL,1,iretou)
  48. CALL ACTOBJ('MMODEL ',IPMODL,1)
  49. IF (IERR.NE.0) RETURN
  50. C
  51. C lecture du mchaml de caractéristiques matériau
  52. C
  53. CALL LIROBJ('MCHAML',IPIN,1,iretou)
  54. CALL ACTOBJ('MCHAML',IPIN,1)
  55. IF (IERR.NE.0) RETURN
  56. CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
  57. IF(IR .NE. 1) CALL ERREUR(KER)
  58. IF(IERR .NE. 0) RETURN
  59. c
  60. C lecture du point vecteur de rotation
  61. C
  62. CALL LIROBJ('POINT',IPROTA,0,iretou)
  63. IF (IERR.NE.0) RETURN
  64. C
  65. C lecture du mot clé pour le calcul direct d'impedance (pour Fourier)
  66. C
  67. CALL LIRMOT(LISMO1,1,NUMLI1,0)
  68. IF (IERR.NE.0) RETURN
  69. CALL LIRMOT(LISMO2,1,NUMLI2,0)
  70. IF (IERR.NE.0) RETURN
  71. C
  72. *********************************************
  73. * Matrice de couplage de Coriolis
  74. **********************************************
  75. CALL CORIO1(IPMODL,IPCHE1,IPROTA,NUMLI1,NUMLI2, IPRIG)
  76. C
  77. IF (IPRIG.NE.0 .AND. IERR.EQ.0) THEN
  78. CALL ECROBJ('RIGIDITE',IPRIG)
  79. ENDIF
  80.  
  81. END
  82.  
  83.  
  84.  
  85.  

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