Télécharger la plaquette Cast3M


fr en

Télécharger stru4.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  stru4.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *           Test Stru4.dgibi: Jeux de données         *
  5. *           ---------------------------------         *
  6. *                                                     *
  7. *******************************************************
  8. *                                                     *
  9. * TEST DE STATIQUE POUR L'ELEMENT TUYAU               *
  10. *                                                     *
  11. * Un tuyau encastre‚ et courbe‚ est soumis a…des      *
  12. * efforts de flexion dans son plan.                   *
  13. *                                                     *
  14. * Ref: Guide VPCS, Test SSLL 06/90                    *
  15. *      Solution analytique; JS Prezemieniecki;        *
  16. *      Theory of matrix structural analysis, New York,*
  17. *      McGraw-Hill, 1968                              *
  18. *                                                     *
  19. * Test effectue‚ le 25/08/92                          *
  20. *                                                     *
  21. *******************************************************
  22.  
  23. GRAPH= 'N';
  24. opti echo 0;
  26. opti elem seg2;
  27.  
  28. ************* GEOMETRIE ********************
  29. P0 = 0. 0. 0.;
  30. PA =  0. 3. 0.;
  31. PB =  3. 0. 0.;
  32. CAB = CERC 20 PA P0 PB;
  33. OEIL = 0. 0. 1000.;
  34.  
  35. SI ( NEG GRAPH 'N');
  36. TRAC OEIL CAB;
  37. FINSI;
  38.  
  39. ************* MODELE ***********************
  40.  
  41. MODL1 = MODE cab 'MECANIQUE' 'ELASTIQUE' 'TUYAU';
  42. MAT1  = MATE MODL1 YOUN 2E11 NU 0.3 EPAI 0.002 RAYO 
  43. 0.01 VECT OEIL;
  44.  
  45. ************* CONDITIONS AUX LIMITES *******
  46.  
  47. BL1 = BLOQU DEPL ROTA PA;
  48.  
  49. ************* RIGIDITE   *******************
  50.  
  51. RG1 = RIGI MODL1 MAT1;
  52. RGTOT = RG1 ET BL1;
  53.  
  54. ************* CHARGEMENT *******************
  55.  
  56. FOR1 = FORC FX 10. FY 5. PB;
  57. MOM1 = MOME MZ 8 PB;
  58. FORT = FOR1 ET MOM1;
  59.  
  60. ************* CALCUL DEPLACEMENT *************
  61.  
  62. DEP = RESOU RGTOT FORT;
  63. UBX = EXTR DEP PB UX;
  64. UBY = EXTR DEP PB UY;
  65. RBZ = EXTR DEP PB RZ;
  66.  
  67.  
  68. ************* CALCUL ANALYTIQUE **************
  69.  
  70. F1 = 10.; F2 = 5.; MZ = 8.; E= 2E11;
  71. I = 4.637E-9; R = 3.;
  72.  
  73. UBX1 = (R**2)/4/E/I*((F1*PI*R)+(2*F2*R)+(4*MZ));
  74. UBY1 = (R**2)/4/E/I*((2*F1*R)+(((3*PI)-(8.))*F2*R)+
  75. (2*(PI-(2.))*MZ));
  76. RBZ1 = R/4/E/I*((4*R*F1)+(2*(PI-(2.))*F2*R)+(2*PI*MZ));
  77.  
  78. ************* COMPARAISON DES RESULTATS *****
  79.  
  80. MESS ' RESULTATS ';
  81. MESS ' --------- ';
  82.  
  84. MESS ' DEPLACEMENT THEORIQUE AU POINT B SUIVANT X;'
  85. UBX1 'M';
  86. MESS ' DEPLACEMENT CALCULE   AU POINT B SUIVANT X;' 
  87. UBX  'M';
  88. R1= 100*(ABS((UBX-UBX1)/UBX1));
  89. MESS '     SOIT UN ECART DE       : '  R1  '%' ;
  91. MESS ' DEPLACEMENT THEORIQUE AU POINT B SUIVANT Y;' 
  92. UBY1 'M';
  93. MESS ' DEPLACEMENT CALCULE   AU POINT B SUIVANT Y;' 
  94. UBY  'M';
  95. R2= 100*(ABS((UBY-UBY1)/UBY1));
  96. MESS '     SOIT UN ECART DE       : '  R2  '%' ;
  98. MESS ' ROTATION THEORIQUE AU POINT B;' RBZ1 'rad';
  99. MESS ' ROTATION CALCULEE  AU POINT B;' RBZ 'rad';
  100. R3= 100*(ABS((RBZ-RBZ1)/RBZ1));
  101. MESS '     SOIT UN ECART DE       : '  R3  '%' ;
  102.  
  103.  
  104. ************* CODE DE FONCTIONNEMENT ********
  105. LR = PROG R1 R2 R3;
  106.  
  107. ERMAX = MAXI LR;
  108.  
  109. SI ( ERMAX < 5);
  110.    ERRE 0;
  111. SINON
  112.    ERRE 5;
  113. FINSI;
  114. FIN;
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  

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