Télécharger la plaquette Cast3M


fr en

Télécharger coriol.eso

Retour à la liste

Numérotation des lignes : 

coriol

  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