Télécharger la plaquette Cast3M


fr en

Télécharger flam1.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  flam1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *           Test Flam1.dgibi: Jeux de données         *
  5. *           ---------------------------------         *
  6. *                                                     *
  7.  
  8. *******************************************************
  9. *  TEST FLAM1
  10. *
  11. *           FLAMBAGE EULERIEN D'UNE POUTRE
  12. *             ENCASTREE A UNE EXTREMITE
  13. *
  14. * Dans cet exemple on se propose d'étudier le flambage
  15. * d'une poutre encastrée à une de ses extrémités et
  16. * soumise à une force parallele à la poutre et
  17. * appliquée à l'extrémité libre de la poutre.
  18. *
  19. * Le flambage est visualisé en fin de calcul et le
  20. * multiplicateur de chargement calculé est comparé à
  21. * celui obtenu par la théorie (41.123).
  22. *******************************************************
  23.  
  24. GRAPH = 'N';
  25.  
  26. SAUT PAGE;
  27. SI (NEG GRAPH 'N') ;
  28.   OPTI ECHO 1   ;
  29.   OPTI TRAC PSC ;
  30. SINO ;
  31.   OPTI ECHO 0   ;
  32. FINSI ;
  33.  
  34. SAUT PAGE;
  35.  
  36. TITRE ' FLAMBAGE D UNE POUTRE ENCASTREE ' ;
  37. OPTI DIME 2 ELEM SEG2  MODE PLAN CONT ;
  38.  
  39.  
  40. TEMPS ;
  41. *
  42. *------------- DEFINITION DU MAILLAGE -----------------
  43. *
  44. DENSITE  1. ;
  45. P1 = 50. 0. ;
  46. P2 = 50. 10. ;
  47. L1 = P1 D P2 ;
  48. *
  49. SI (NEG GRAPH 'N');
  50.   TRAC 'QUAL' L1 ;
  51. FINSI;
  52. *
  53. *--------- CREATION DU MODELE ET DE MMODEL ------------
  54. *
  55. MOD1 = MODE L1 MECANIQUE ELASTIQUE COQ2;
  56.  
  57. *
  58. * -------- DEFINITION DES CARACTERISTIQUES ------------
  59. * ----------- MATERIELLES ET GEOMETRIQUES  ------------
  60. *
  61. MAT1   = MATE MOD1 YOUN   20000.    NU 0.3     ;
  62. CAR1   = CARA MOD1 EPAI 1. ;
  63. MAT1=MAT1 ET CAR1;
  64. *
  65. *--------- CALCUL DES RIGIDITES ELEMENTAIRES ----------
  66. *----------- ET DEFINITIONS DES BLOCAGES --------------
  67. *
  68. RIG1   = RIGI  MOD1 MAT1 ;
  69. BL1    = BLOQUE DEPL ROTA P1 ;
  70. RIG11  = RIG1 ET BL1;
  71. *
  72. *------------- DEFINITION DU CHARGEMENT ---------------
  73. *
  74. FE1 = FORC ( 0 -1. ) P2 ;
  75. *
  76. *------- RESOLUTION ET CALCUL DES CONTRAINTES ---------
  77. *
  78. DEP1 = RESO  RIG11    FE1 ;
  79. SIG1 = SIGM 'LINE' MOD1 MAT1  DEP1 ;
  80. *
  81. *--- CALCUL DE LA MATRICE DES CONTRAINTES INITIALES ---
  82. *
  83. KSI1   = KSIGMA  MOD1 SIG1 CAR1 'FLAM' ;
  84. *
  85. *----------- INITIALISATION DE LA TABLE ---------------
  86. *--------- EN ENTREE DE LA PROCEDURE FLAMBAGE ---------
  87. *
  88. ETAB        = TABLE    ;
  89. ETAB.'CLIM' = BL1 ;
  90. ETAB.'SIG1' = SIG1;
  91. ETAB.'MATE' = MAT1;
  92. ETAB.'LAM1' = 0.001 ;
  93. ETAB.'LAM2' = 100. ;
  94. ETAB.'NMOD' = 1   ;
  95. ETAB.'OBJM' = MOD1;
  96. *
  97. *--------- APPEL A LA PROCEDURE DE FLAMBAGE -----------
  98. *
  99. SAUT PAGE ;
  100. STAB = FLAMBAGE ETAB ;
  101. *
  102. *------------ RECUPERATION DES RESULTATS --------------
  103. *
  104. LA1  =  STAB . 1 . LAMB ;
  105. SAUT PAGE ;
  106. *
  107. *------------ DESSIN DU MODE DE FLAMBAGE --------------
  108. *
  109. MM1 = STAB . 1 . DEPL ;
  110. DE  = DEFO MM1 L1 1. VERT ;
  111. DA  = DEFO MM1 L1 0. NOIR  ;
  112.  
  113. SI (NEG GRAPH 'N');
  114.   TRAC ( DE ET DA ) ;
  115. FINSI;
  116.  
  117. LREF = 41.123;
  118. ERR=100*(ABS(LREF - LA1)/LREF);
  119. *
  120. MESS ' ' ;
  121. MESS ' ' ;
  122. MESS ' MULTIPLICATEUR DU CHARGEMENT ' ;
  123. MESS '                              ' ;
  124. MESS '   SOLUTION ANALYTIQUE  : ' LREF;
  125. MESS '   SOLUTION CALCULEE    : ' LA1 ;
  126. MESS '     SOIT UN ECART DE                : ' ERR '%' ;
  127.  
  128. TEMPS ;
  129.  
  130. *------------ CODE DE BON FONCTIONNEMENT --------------
  131. SI (ERR < 1e-3);
  132. ERRE 0;
  133. SINON;
  134. ERRE 5;
  135. FINSI;
  136.  
  137. FIN;
  138.  
  139.  
  140.  
  141.  
  142.  
  143.  
  144.  
  145.  
  146.  

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