Télécharger four1.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : four1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Four1.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7.  
  8. *******************************************************
  9. * Test four1.dgibi: Jeux de données *
  10. * --------------------------------- *
  11. *******************************************************
  12.  
  13. * SI GRAPH = N PAS DE GRAPHIQUE AFFICHE
  14. * SINON SI GRAPH DIFFERENT DE N TOUS
  15. * LES GRAPHIQUES SONT AFFICHES
  16.  
  17. GRAPH = 'N' ;
  18.  
  19. SAUT PAGE;
  20. SI (NEG GRAPH 'N') ;
  21. OPTI ECHO 1 ;
  22. OPTI TRAC PSC ;
  23. SINO ;
  24. OPTI ECHO 1 ;
  25. FINSI ;
  26.  
  27. SAUT PAGE;
  28. *******************************************************
  29. * TEST FOUR1
  30. *
  31. * CYLINDRE INFINI SOUS PRESSION EXTERNE(NU=0)
  32. *
  33. * Soit un cylindre infini soumis à une pression externe
  34. *
  35. * Une analyse de flambage permet de déterminer la charge critique
  36. * associée au deuxième mode de Fourier de la structure :
  37. * (u = u°*cos(2*Teta) )
  38. *
  39. * Le calcul est fait deux fois
  40. * 1- Avec les précontraintes dues à la pression externe
  41. * 2- Avec en plus la prise en compte des forces suiveuses (Evolution
  42. * des efforts de pression avec le déplacement de la structure)
  43. *
  44. *******************************************************
  45. TITRE 'CYLINDRE INFINI SOUS PRESSION EXTERNE';
  46. OPTI DIME 2 ELEM SEG2 MODE FOUR 0;
  47.  
  48. TEMPS ;
  49.  
  50. *
  51. *----------- CONSTRUCTION DE LA GEOMETRIE -------------
  52. *
  53. PA=1000. 0.;PB=1000. 100.;PO1=0. 0.;PO2=0. 100.;
  54. CYL = PA D 4 PB;
  55. *
  56. SI (NEG GRAPH 'N');
  57. TRAC 'QUAL' ECLA ELEM CYL;
  58. FINSI;
  59. *
  60. MOD1 = MODE CYL MECANIQUE ELASTIQUE COQ2;
  61.  
  62. *
  63. *--- DECLARATION DE FOURIER NOHARM POUR LES OBJETS ----
  64. *------- QUI SERONT UTILISES POUR PLUSIEURS ----------
  65. *--------------- NUMEROS D'HARMONIQUE -----------------
  66. *
  67. OPTI MODE FOUR NOHARM;
  68. *
  69. *
  70. *-------- CONDITIONS AUX LIMITES SYMETRIQUES ----------
  71. *
  72. SYMB = SYMT CYL DEPL ROTA PA PO1 0.5;
  73. SYMH = SYMT CYL DEPL ROTA PB PO2 0.5;
  74. CDL = SYMB ET SYMH;
  75. *
  76. *----------- MATERIAU ET CARACTERISTIQUES -------------
  77. *
  78. MAT = MATE MOD1 YOUN 20000. NU 0.;
  79. CAR = CARA MOD1 EPAI 1.;
  80. MAT = MAT ET CAR ;
  81. *
  82. *----------- DECLARATION DE FOURIER MODE 0 ------------
  83. *----------- POUR LE CALCUL DES CONTRAINTES -----------
  84. *
  85. OPTI MODE FOUR 0;
  86. MOP = 'MODE' CYL 'CHARGEMENT' 'PRESSION' 'COQ2' ;
  87. MAP = 'PRES' MOP 'PRES' 1. ;
  88. MOD1 = MOD1 ET MOP ; MAT = MAT ;
  89. FP = 'BSIG' MOP MAP ;
  90. *
  91. RIG = RIGI MOD1 MAT;
  92.  
  93. RIGADD = RIGI UT 1.E-4 PA ;
  94.  
  95. AAA = RIG ET CDL ET RIGADD ;
  96. U = RESO AAA FP ;
  97. SIG=SIGMA U MOD1 MAT;
  98. *
  99. *----------- DECLARATION DE FOURIER MODE 2 ------------
  100. *------------ POUR L'ANALYSE DE FLAMBAGE --------------
  101. *
  102. OPTI MODE FOUR 2;
  103.  
  104. PP = MANU CHPO CYL 1 P 1. ;
  105. *
  106. MKSI= KSIGMA MOD1 MAT (SIG * -1.) FLAM;
  107. RIG = RIGI MOD1 MAT ;
  108. *
  109. *------- RECHERCHE DE LA 1ERE FREQUENCE PROPRE --------
  110. *
  111. MODF = VIBR PROC (PROG 0.) (RIG ET CDL) MKSI;
  112. W1 = MODF . MODES . 1 . FREQUENCE;
  113. LAMBDA1=(W1 * PI * 2) ** 2 * (SIGN W1);
  114. ERR1 = 100*(ABS((2E-5/3) - LAMBDA1)/(2E-5/3));
  115.  
  116. MKP = KP PP MOP FLAM ;
  117. MODF = VIBR PROC (PROG 0.) (RIG ET CDL) (MKSI ET MKP);
  118. W2 = MODF . MODES . 1 . FREQUENCE;
  119. LAMBDA2=(W2 * PI * 2) ** 2 * (SIGN W2);
  120. ERR2 = 100*(ABS(5E-6 - LAMBDA2)/5E-6 );
  121.  
  122. SAUT PAGE ;
  123. MESS 'K(SIG) SEUL : ON DOIT TROUVER LAMBDA= 6.66E-6' ;
  124. MESS ' LE CALCUL DONNE LAMBDA=' LAMBDA1 ;
  125. MESS ' SOIT UN ECART DE : ' ERR1 '%' ;
  126.  
  127. MESS 'K(SIG)+K(P) : ON DOIT TROUVER LAMBDA= 5E-6' ;
  128. MESS ' LE CALCUL DONNE LAMBDA=' LAMBDA2 ;
  129. MESS ' SOIT UN ECART DE : ' ERR2 '%' ;
  130.  
  131. TEMPS ;
  132. *------------- CODE DE FONCTIONNEMENT ----------------
  133.  
  134. ERRMAX=MAXI (PROG ERR1 ERR2);
  135.  
  136. SI (ERRMAX < 2e-5 );
  137. ERRE 0;
  138. SINON ;
  139. ERRE 5;
  140. FINSI;
  141.  
  142. FIN;
  143.  
  144.  
  145.  
  146.  
  147.  
  148.  
  149.  
  150.  
  151.  
  152.  
  153.  
  154.  
  155.  

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