C KENTR     SOURCE    PV        09/03/12    21:26:08     6325           
      SUBROUTINE KENTR                    
C     
C***********************************************************************                                                                          
C                                                                      *
C                                                                      *
C      FABRICATION DES MATRICES DE RAIDEUR CENTRIFUGE ET D'EULER       *                                                                           *
C      (CALCUL DANS LE REPERE TOURNANT                                 *
C  Syntaxe :                                                           *
C  --------                                                            *
C                                                                      *
C    RIG1 = KCEN MODL1 MAT1 VEC1 (FLAM)                                *
C                                                                      *
C  RIG1 : matrice de raideur centrifuge  (TYPE rigidité)               *
C  MODL1: objet MMODEL, modèle                                         *
C  MAT1 : objet MCHAML, caractéristiques matériau                      *
C  VEC1 : objet de type POINT, vecteur vitesse de rotation             * 
C  FLAM : mot clé si on se sert de la matrice pour faire du flambage   *
C                                                                      * 
C            Didier COMBESCURE       Mars 2003                         *
C                                                                      *
C***********************************************************************
C                                                                              
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)                                                 
-INC CCOPTIO                                                                    
-INC SMCHAML   
-INC CCHAMP
c                                                                         
*                                                                          
      CHARACTER*4 LPARA1(2), LPARA2(1)
*  
      DATA LPARA1/'CENT','EULE'/
      DATA LPARA2/'FLAM'/
      NPARA1 = 2                                      
*                                          
*********************************************
*
*    Matrice de raideur centrifuge et d'Euler
*
**********************************************
C 
      CALL LIRMOT(LPARA1,NPARA1,INDIC1,0)
C
      IF (INDIC1.NE.0) THEN
C
           IPCHE1=0
           IRET=0    
C
C Test flambage
C 
C                                                                               
      CALL LIRMOT(LPARA2,1,IFLAM,0)                                             
      IF (IERR.NE.0) RETURN           
C
C  lecture du point vecteur de rotation
C
          IF (IFOUR.NE.1) THEN
            CALL LIROBJ('POINT',IPROTA,1,IRT1)
            IF (IRT1.EQ.0) RETURN
          ELSE
            CALL LIROBJ('POINT',IPROTA,0,IRT1)
            IF (IRT1.EQ.0) IPROTA =0
          ENDIF
C
C  lecture du modèle
C
           CALL LIROBJ('MMODEL',IPMODL,1,IRT2)
           IF (IRT2.EQ.0) RETURN
C
C  lecture du mchaml de caractéristiques matériau
C
           CALL LIROBJ('MCHAML',IPCHE1,1,IRT3)
           IF (IRT3.EQ.0) RETURN
C
      IF (INDIC1.EQ.1) THEN
C
C    Matrice de raideur centrifuge
C
C
           CALL KCENT1(IPMODL,IPCHE1,IPROTA,IPRIG,IFLAM,IRET)
C
      ELSEIF (INDIC1.EQ.2) THEN 
C
C
C    Matrice de raideur d'Euler
C
           CALL KEULE1(IPMODL,IPCHE1,IPROTA,IPRIG,IFLAM,IRET)
      ELSE
            IRET = 0
      ENDIF
C
           IF (IRET.EQ.1) THEN
             CALL ECROBJ('RIGIDITE',IPRIG)
           ENDIF
C
         ENDIF
C
       RETURN
C                                                                             
      END