Esope source files

   1 : C CORIOL    SOURCE    FANDEUR   11/07/19    21:15:52     7042           
   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
  46 :       CALL LIROBJ('MMODEL',IPMODL,1,iretou)
  47 :       IF (IERR.NE.0) RETURN
  48 : C
  49 : C  lecture du mchaml de caractéristiques matériau
  50 : C
  50 :       CALL LIROBJ('MCHAML',IPCHE1,1,iretou)
  51 :       IF (IERR.NE.0) RETURN
  52 : c
  53 : C  lecture du point vecteur de rotation
  54 : C
  54 :       CALL LIROBJ('POINT',IPROTA,0,iretou)
  55 :       IF (IERR.NE.0) RETURN
  56 : C
  57 : C  lecture du mot clé pour le calcul direct d'impedance (pour Fourier)
  58 : C
  58 :       CALL LIRMOT(LISMO1,1,NUMLI1,0)
  59 :       IF (IERR.NE.0) RETURN
  59 :       CALL LIRMOT(LISMO2,1,NUMLI2,0)
  60 :       IF (IERR.NE.0) RETURN
  61 : C
  62 : *********************************************
  63 : *    Matrice de couplage de Coriolis
  64 : **********************************************
  64 :       CALL CORIO1(IPMODL,IPCHE1,IPROTA,NUMLI1,NUMLI2, IPRIG)
  65 : C
  66 :       IF (IPRIG.NE.0 .AND. IERR.EQ.0) THEN
  66 :         CALL ECROBJ('RIGIDITE',IPRIG)
  67 :       ENDIF
  68 : C
  69 :       RETURN
  70 :       END
  71 : 
  72 :