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.  
  31. -INC PPARAM
  32. -INC CCOPTIO
  33. c
  34. CHARACTER*4 LISMO1(1),LISMO2(1)
  35. DATA LISMO1 / 'HARM' /
  36. DATA LISMO2 / 'FLAM' /
  37.  
  38. IPMODL = 0
  39. IPCHE1 = 0
  40. IPROTA = 0
  41.  
  42. NUMLI1 = 0
  43. NUMLI2 = 0
  44.  
  45. IPRIG = 0
  46. C
  47. C lecture du modèle
  48. C
  49. CALL LIROBJ('MMODEL ',IPMODL,1,iretou)
  50. CALL ACTOBJ('MMODEL ',IPMODL,1)
  51. IF (IERR.NE.0) RETURN
  52. C
  53. C lecture du mchaml de caractéristiques matériau
  54. C
  55. CALL LIROBJ('MCHAML',IPIN,1,iretou)
  56. CALL ACTOBJ('MCHAML',IPIN,1)
  57. IF (IERR.NE.0) RETURN
  58. CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
  59. IF(IR .NE. 1) CALL ERREUR(KER)
  60. IF(IERR .NE. 0) RETURN
  61. c
  62. C lecture du point vecteur de rotation
  63. C
  64. CALL LIROBJ('POINT',IPROTA,0,iretou)
  65. IF (IERR.NE.0) RETURN
  66. C
  67. C lecture du mot clé pour le calcul direct d'impedance (pour Fourier)
  68. C
  69. CALL LIRMOT(LISMO1,1,NUMLI1,0)
  70. IF (IERR.NE.0) RETURN
  71. CALL LIRMOT(LISMO2,1,NUMLI2,0)
  72. IF (IERR.NE.0) RETURN
  73. C
  74. *********************************************
  75. * Matrice de couplage de Coriolis
  76. **********************************************
  77. CALL CORIO1(IPMODL,IPCHE1,IPROTA,NUMLI1,NUMLI2, IPRIG)
  78. C
  79. IF (IPRIG.NE.0 .AND. IERR.EQ.0) THEN
  80. CALL ECROBJ('RIGIDITE',IPRIG)
  81. ENDIF
  82.  
  83. END
  84.  
  85.  
  86.  
  87.  

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