Télécharger kentr.eso

Retour à la liste

Numérotation des lignes :

  1. C KENTR SOURCE CB215821 19/07/30 21:17:05 10273
  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. CALL ACTOBJ('MMODEL ',IPMODL,1)
  67. IF(IERR .NE. 0) RETURN
  68. C
  69. C lecture du mchaml de caractéristiques matériau
  70. C
  71. CALL LIROBJ('MCHAML ',IPIN,1,IRT3)
  72. CALL ACTOBJ('MCHAML ',IPIN,1)
  73. IF(IERR .NE. 0) RETURN
  74. CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
  75. IF(IR .NE. 1) CALL ERREUR(KER)
  76. IF(IERR .NE. 0) RETURN
  77. C
  78. IF (INDIC1.EQ.1) THEN
  79. C
  80. C Matrice de raideur centrifuge
  81. CALL KCENT1(IPMODL,IPCHE1,IPROTA,IPRIG,IFLAM,IRET)
  82. IF(IERR .NE. 0) RETURN
  83. C
  84. ELSEIF (INDIC1.EQ.2) THEN
  85. C
  86. C Matrice de raideur d'Euler
  87. CALL KEULE1(IPMODL,IPCHE1,IPROTA,IPRIG,IFLAM,IRET)
  88. IF(IERR .NE. 0) RETURN
  89.  
  90. ELSE
  91. IRET = 0
  92. ENDIF
  93. C
  94. IF (IRET.EQ.1) CALL ECROBJ('RIGIDITE',IPRIG)
  95. ENDIF
  96.  
  97. END
  98.  
  99.  
  100.  

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