Esope source files

   1 : C KENTR     SOURCE    PV        09/03/12    21:26:08     6325           
   2 :       SUBROUTINE KENTR                    
   3 : C     
   4 : C***********************************************************************                                                                          
   5 : C                                                                      *
   6 : C                                                                      *
   7 : C      FABRICATION DES MATRICES DE RAIDEUR CENTRIFUGE ET D'EULER       *                                                                           *
   8 : C      (CALCUL DANS LE REPERE TOURNANT                                 *
   9 : C  Syntaxe :                                                           *
  10 : C  --------                                                            *
  11 : C                                                                      *
  12 : C    RIG1 = KCEN MODL1 MAT1 VEC1 (FLAM)                                *
  13 : C                                                                      *
  14 : C  RIG1 : matrice de raideur centrifuge  (TYPE rigidité)               *
  15 : C  MODL1: objet MMODEL, modèle                                         *
  16 : C  MAT1 : objet MCHAML, caractéristiques matériau                      *
  17 : C  VEC1 : objet de type POINT, vecteur vitesse de rotation             * 
  18 : C  FLAM : mot clé si on se sert de la matrice pour faire du flambage   *
  19 : C                                                                      * 
  20 : C            Didier COMBESCURE       Mars 2003                         *
  21 : C                                                                      *
  22 : C***********************************************************************
  23 : C                                                                              
  24 :       IMPLICIT INTEGER(I-N)
  25 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  26 : -INC CCOPTIO                                                                    
  27 : -INC SMCHAML   
  28 : -INC CCHAMP
  29 : c                                                                         
  30 : *                                                                          
  31 :       CHARACTER*4 LPARA1(2), LPARA2(1)
  32 : *  
  33 :       DATA LPARA1/'CENT','EULE'/
  34 :       DATA LPARA2/'FLAM'/
  35 :       NPARA1 = 2                                      
  36 : *                                          
  37 : *********************************************
  38 : *
  39 : *    Matrice de raideur centrifuge et d'Euler
  40 : *
  41 : **********************************************
  42 : C 
  42 :       CALL LIRMOT(LPARA1,NPARA1,INDIC1,0)
  43 : C
  44 :       IF (INDIC1.NE.0) THEN
  45 : C
  46 :            IPCHE1=0
  47 :            IRET=0    
  48 : C
  49 : C Test flambage
  50 : C 
  51 : C                                                                               
  51 :       CALL LIRMOT(LPARA2,1,IFLAM,0)                                             
  52 :       IF (IERR.NE.0) RETURN           
  53 : C
  54 : C  lecture du point vecteur de rotation
  55 : C
  56 :           IF (IFOUR.NE.1) THEN
  56 :             CALL LIROBJ('POINT',IPROTA,1,IRT1)
  57 :             IF (IRT1.EQ.0) RETURN
  58 :           ELSE
  58 :             CALL LIROBJ('POINT',IPROTA,0,IRT1)
  59 :             IF (IRT1.EQ.0) IPROTA =0
  60 :           ENDIF
  61 : C
  62 : C  lecture du modèle
  63 : C
  63 :            CALL LIROBJ('MMODEL',IPMODL,1,IRT2)
  64 :            IF (IRT2.EQ.0) RETURN
  65 : C
  66 : C  lecture du mchaml de caractéristiques matériau
  67 : C
  67 :            CALL LIROBJ('MCHAML',IPCHE1,1,IRT3)
  68 :            IF (IRT3.EQ.0) RETURN
  69 : C
  70 :       IF (INDIC1.EQ.1) THEN
  71 : C
  72 : C    Matrice de raideur centrifuge
  73 : C
  74 : C
  74 :            CALL KCENT1(IPMODL,IPCHE1,IPROTA,IPRIG,IFLAM,IRET)
  75 : C
  76 :       ELSEIF (INDIC1.EQ.2) THEN 
  77 : C
  78 : C
  79 : C    Matrice de raideur d'Euler
  80 : C
  80 :            CALL KEULE1(IPMODL,IPCHE1,IPROTA,IPRIG,IFLAM,IRET)
  81 :       ELSE
  82 :             IRET = 0
  83 :       ENDIF
  84 : C
  85 :            IF (IRET.EQ.1) THEN
  85 :              CALL ECROBJ('RIGIDITE',IPRIG)
  86 :            ENDIF
  87 : C
  88 :          ENDIF
  89 : C
  90 :        RETURN
  91 : C                                                                             
  92 :       END                                                                       
  93 :  
  94 : 
  95 :  
  96 :  
  97 :  
  98 :  
  99 :  
 100 :  
 101 :  
 102 :  
 103 :  
 104 :  
 105 :  
 106 :  
 107 :  
 108 :