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. OPTI EPSI LINEAIRE;
  48.  
  49. TEMPS ;
  50.  
  51. *
  52. *----------- CONSTRUCTION DE LA GEOMETRIE -------------
  53. *
  54. PA=1000. 0.;PB=1000. 100.;PO1=0. 0.;PO2=0. 100.;
  55. CYL = PA D 4 PB;
  56. *
  57. SI (NEG GRAPH 'N');
  58. TRAC 'QUAL' ECLA ELEM CYL;
  59. FINSI;
  60. *
  61. MOD1 = MODE CYL MECANIQUE ELASTIQUE COQ2;
  62.  
  63. *
  64. *--- DECLARATION DE FOURIER NOHARM POUR LES OBJETS ----
  65. *------- QUI SERONT UTILISES POUR PLUSIEURS ----------
  66. *--------------- NUMEROS D'HARMONIQUE -----------------
  67. *
  68. OPTI MODE FOUR NOHARM;
  69. *
  70. *
  71. *-------- CONDITIONS AUX LIMITES SYMETRIQUES ----------
  72. *
  73. SYMB = SYMT CYL DEPL ROTA PA PO1 0.5;
  74. SYMH = SYMT CYL DEPL ROTA PB PO2 0.5;
  75. CDL = SYMB ET SYMH;
  76. *
  77. *----------- MATERIAU ET CARACTERISTIQUES -------------
  78. *
  79. MAT = MATE MOD1 YOUN 20000. NU 0.;
  80. CAR = CARA MOD1 EPAI 1.;
  81. MAT = MAT ET CAR ;
  82. *
  83. *----------- DECLARATION DE FOURIER MODE 0 ------------
  84. *----------- POUR LE CALCUL DES CONTRAINTES -----------
  85. *
  86. OPTI MODE FOUR 0;
  87. MOP = 'MODE' CYL 'CHARGEMENT' 'PRESSION' 'COQ2' ;
  88. MAP = 'PRES' MOP 'PRES' 1. ;
  89. MOD1 = MOD1 ET MOP ; MAT = MAT ;
  90. FP = 'BSIG' MOP MAP ;
  91. *
  92. RIG = RIGI MOD1 MAT;
  93.  
  94. RIGADD = RIGI UT 1.E-4 PA ;
  95.  
  96. AAA = RIG ET CDL ET RIGADD ;
  97. U = RESO AAA FP ;
  98. SIG=SIGMA U MOD1 MAT;
  99. *
  100. *----------- DECLARATION DE FOURIER MODE 2 ------------
  101. *------------ POUR L'ANALYSE DE FLAMBAGE --------------
  102. *
  103. OPTI MODE FOUR 2;
  104.  
  105. PP = MANU CHPO CYL 1 P 1. ;
  106. *
  107. MKSI= KSIGMA MOD1 MAT (SIG * -1.) FLAM;
  108. RIG = RIGI MOD1 MAT ;
  109. *
  110. *------- RECHERCHE DE LA 1ERE FREQUENCE PROPRE --------
  111. *
  112. MODF = VIBR PROC (PROG 0.) (RIG ET CDL) MKSI;
  113. W1 = MODF . MODES . 1 . FREQUENCE;
  114. LAMBDA1=(W1 * PI * 2) ** 2 * (SIGN W1);
  115. ERR1 = 100*(ABS((2E-5/3) - LAMBDA1)/(2E-5/3));
  116.  
  117. MKP = KP PP MOP FLAM ;
  118. MODF = VIBR PROC (PROG 0.) (RIG ET CDL) (MKSI ET MKP);
  119. W2 = MODF . MODES . 1 . FREQUENCE;
  120. LAMBDA2=(W2 * PI * 2) ** 2 * (SIGN W2);
  121. ERR2 = 100*(ABS(5E-6 - LAMBDA2)/5E-6 );
  122.  
  123. SAUT PAGE ;
  124. MESS 'K(SIG) SEUL : ON DOIT TROUVER LAMBDA= 6.66E-6' ;
  125. MESS ' LE CALCUL DONNE LAMBDA=' LAMBDA1 ;
  126. MESS ' SOIT UN ECART DE : ' ERR1 '%' ;
  127.  
  128. MESS 'K(SIG)+K(P) : ON DOIT TROUVER LAMBDA= 5E-6' ;
  129. MESS ' LE CALCUL DONNE LAMBDA=' LAMBDA2 ;
  130. MESS ' SOIT UN ECART DE : ' ERR2 '%' ;
  131.  
  132. TEMPS ;
  133. *------------- CODE DE FONCTIONNEMENT ----------------
  134.  
  135. ERRMAX=MAXI (PROG ERR1 ERR2);
  136.  
  137. SI (ERRMAX < 2e-5 );
  138. ERRE 0;
  139. SINON ;
  140. ERRE 5;
  141. FINSI;
  142.  
  143. FIN;
  144.  
  145.  
  146.  
  147.  
  148.  
  149.  
  150.  
  151.  
  152.  
  153.  
  154.  

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