Télécharger kentr.eso

Retour à la liste

Numérotation des lignes :

  1. C KENTR SOURCE CB215821 16/12/05 21:39:53 9237
  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
  43. CALL LIRMOT(LPARA1,NPARA1,INDIC1,0)
  44. C
  45. IF (INDIC1 .NE. 0) THEN
  46. IPCHE1=0
  47. IRET=0
  48. C
  49. C Test flambage
  50. CALL LIRMOT(LPARA2,1,IFLAM,0)
  51. IF (IERR.NE.0) RETURN
  52. C
  53. C lecture du point vecteur de rotation
  54. C
  55. IF (IFOUR.NE.1) THEN
  56. CALL LIROBJ('POINT',IPROTA,1,IRT1)
  57. IF(IERR .NE. 0) RETURN
  58. ELSE
  59. CALL LIROBJ('POINT',IPROTA,0,IRT1)
  60. IF (IRT1.EQ.0) IPROTA =0
  61. ENDIF
  62. C
  63. C lecture du modèle
  64. C
  65. CALL LIROBJ('MMODEL',IPMODL,1,IRT2)
  66. IF(IERR .NE. 0) RETURN
  67. C
  68. C lecture du mchaml de caractéristiques matériau
  69. C
  70. CALL LIROBJ('MCHAML',IPIN,1,IRT3)
  71. IF(IERR .NE. 0) RETURN
  72. CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
  73. IF(IR .NE. 1) CALL ERREUR(KER)
  74. IF(IERR .NE. 0) RETURN
  75. C
  76. IF (INDIC1.EQ.1) THEN
  77. C
  78. C Matrice de raideur centrifuge
  79. CALL KCENT1(IPMODL,IPCHE1,IPROTA,IPRIG,IFLAM,IRET)
  80. IF(IERR .NE. 0) RETURN
  81. C
  82. ELSEIF (INDIC1.EQ.2) THEN
  83. C
  84. C Matrice de raideur d'Euler
  85. CALL KEULE1(IPMODL,IPCHE1,IPROTA,IPRIG,IFLAM,IRET)
  86. IF(IERR .NE. 0) RETURN
  87.  
  88. ELSE
  89. IRET = 0
  90. ENDIF
  91. C
  92. IF (IRET.EQ.1) THEN
  93. CALL ECROBJ('RIGIDITE',IPRIG)
  94. IF(IERR .NE. 0) RETURN
  95. ENDIF
  96. C
  97. ENDIF
  98. C
  99. RETURN
  100. C
  101. END
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  

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