Télécharger la plaquette Cast3M


fr en

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 ;
  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 '%' ;
  128.  
  129. MESS 'K(SIG)+K(P)  : ON DOIT TROUVER LAMBDA=  5E-6'     ;
  130. MESS '               LE CALCUL DONNE LAMBDA=' LAMBDA2   ;
  131. MESS '     SOIT UN ECART DE                : ' ERR2 '%' ;
  133.  
  134. TEMPS ;
  135. *-------------  CODE DE FONCTIONNEMENT ----------------
  136.  
  137. ERRMAX=MAXI (PROG ERR1 ERR2);
  138.  
  139. SI (ERRMAX < 2e-5 );
  140.   ERRE 0;
  141. SINON ;
  142.   ERRE 5;
  143. FINSI;
  144.  
  145. FIN;
  146.  
  147.  
  148.  
  149.  
  150.  
  151.  
  152.  
  153.  
  154.  
  155.  
  156.  
  157.  
  158.  

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