Télécharger coriol.eso

Retour à la liste

Numérotation des lignes :

  1. C CORIOL SOURCE CB215821 16/12/05 21:15:13 9237
  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. IF (IERR.NE.0) RETURN
  49. C
  50. C lecture du mchaml de caractéristiques matériau
  51. C
  52. CALL LIROBJ('MCHAML',IPIN,1,iretou)
  53. IF (IERR.NE.0) RETURN
  54. CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
  55. IF(IR .NE. 1) CALL ERREUR(KER)
  56. IF(IERR .NE. 0) RETURN
  57. c
  58. C lecture du point vecteur de rotation
  59. C
  60. CALL LIROBJ('POINT',IPROTA,0,iretou)
  61. IF (IERR.NE.0) RETURN
  62. C
  63. C lecture du mot clé pour le calcul direct d'impedance (pour Fourier)
  64. C
  65. CALL LIRMOT(LISMO1,1,NUMLI1,0)
  66. IF (IERR.NE.0) RETURN
  67. CALL LIRMOT(LISMO2,1,NUMLI2,0)
  68. IF (IERR.NE.0) RETURN
  69. C
  70. *********************************************
  71. * Matrice de couplage de Coriolis
  72. **********************************************
  73. CALL CORIO1(IPMODL,IPCHE1,IPROTA,NUMLI1,NUMLI2, IPRIG)
  74. C
  75. IF (IPRIG.NE.0 .AND. IERR.EQ.0) THEN
  76. CALL ECROBJ('RIGIDITE',IPRIG)
  77. ENDIF
  78. C
  79. RETURN
  80. END
  81.  
  82.  
  83.  

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