Télécharger comp1_fourier.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : comp1_fourier.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. ************************************************************************
  6. * CYLINDRE COMPOSITE BICOUCHE
  7. * FIBRES ENROULEES -45/+45 AUTOUR DE L'AXE
  8. * PRESSION INTERNE
  9. *
  10. *
  11. * Un cylindre bloqué à sa base en déplacement suivant
  12. * l'axe Z est soumis à une pression interne.
  13. * 2D Fourier COQ2
  14. * Ref : Rapport CEA , DEMT 85-482 , M. Hittinger, 1985
  15. * Remise a plat : BP, 2017-01-23
  16. *
  17. ************************************************************************
  18.  
  19. * SI GRAPH = N PAS DE GRAPHIQUE AFFICHE
  20. * SINON SI GRAPH DIFFERENT DE N TOUS
  21. * LES GRAPHIQUES SONT AFFICHES
  22.  
  23. GRAPH = 'N' ;
  24. * GRAPH = 'O';
  25.  
  26. SAUT PAGE;
  27. SI (NEG GRAPH 'N') ;
  28. OPTI ECHO 1 ;
  29. * OPTI TRAC X ;
  30. OPTI TRAC 'PSC' EPTR 5;
  31. SINO ;
  32. OPTI ECHO 0 ;
  33. FINSI ;
  34.  
  35. TITRE 'CYLINDRE COMPOSITE BICOUCHE SOUS PRESSION INTERNE';
  36. OPTION DIME 2 MODE FOUR 0 ELEM SEG2;
  37.  
  38. DENS 0.1 ;
  39.  
  40. NAX = 1;
  41. NAX = 4;
  42.  
  43. *______________________________________________________
  44. * GEOMETRIE
  45. *______________________________________________________
  46.  
  47. R = 1.05 ; H = 1. ;
  48. PA = R 0.; PB = R H ;
  49. O1 = 0. 0. ; O2 = 0. 1.;
  50.  
  51. CYL = PA DROIT NAX PB ;
  52.  
  53. SI (NEG GRAPH 'N') ;
  54. TRAC CYL 'QUAL' ;
  55. FINSI ;
  56. *______________________________________________________
  57. * DESCRIPTION DU MATERIAU ORTHOTROPE : COUCHE 1
  58. *______________________________________________________
  59.  
  60. MOD1 = MODE CYL MECANIQUE ELASTIQUE ORTHOTROPE COQ2 'CONS' 'couche 1';
  61. MAT1 = MATE MOD1 'DIRE' O2 'INCL' 45
  62. 'YG1' 7.E6 'YG2' 1.3E6 NU12 0.28 G12 5.E5
  63. 'EPAI' 0.05 'EXCE' 0.025 ;
  64.  
  65. * verification graphique de l'orientation :
  66. * 2D fourier coque => 3 vecteurs : V1, V2 et V3
  67. v123 = VLOC MOD1 MAT1;
  68. * mo123 = mots 'V1R' 'V1Z' 'V1T' 'V2R' 'V2Z' 'V2T' 'V3R' 'V3Z' 'V3T';
  69. mo123 = mots 'V1R' 'V1Z' 'V2R' 'V2Z' 'V3R' 'V3Z' ;
  70. ve123 = VECT v123 MOD1 0.07 mo123 (mots 'AZUR' 'OR' 'POUR');
  71. SI (NEG GRAPH 'N');
  72. TITRE 'COMP1 couche 1 : V1(AZUR) V2 (JAUNE) V3(ROUG)';
  73. TRACE ve123 CYL 'CACH';
  74. FINSI ;
  75.  
  76. *______________________________________________________
  77. * DESCRIPTION DU MATERIAU ORTHOTROPE : COUCHE 2
  78. *______________________________________________________
  79. MOD2 = MODE CYL MECANIQUE ELASTIQUE ORTHOTROPE COQ2 'CONS' 'couche 2' ;
  80. MAT2 = MATE MOD2 'DIRE' O2 'INCL' -45
  81. 'YG1' 7.E6 'YG2' 1.3E6 NU12 0.28 G12 5.E5
  82. 'EPAI' 0.05 'EXCE' -0.025 ;
  83.  
  84. * verification graphique de l'orientation :
  85. v123 = VLOC MOD2 MAT2;
  86. ve123 = VECT v123 MOD2 0.07 mo123 (mots 'AZUR' 'BRON' 'POUR');
  87. SI (NEG GRAPH 'N');
  88. TITRE 'COMP1 couche 2 : V1(AZUR) V2 (JAUNE) V3(ROUG)';
  89. TRACE ve123 CYL 'CACH';
  90. FINSI ;
  91.  
  92. * chargement
  93. MOP = 'MODE' CYL 'CHARGEMENT' 'PRESSION' COQ2 ;
  94. * on met -1 pour etre orienté vers les r>0
  95. MAP = 'MATE' MOP 'PRES' -1. ;
  96.  
  97. *______________________________________________________
  98. * CREATION DE LA RIGIDITE
  99. *______________________________________________________
  100.  
  101.  
  102. MODORT = MOD1 ET MOD2 'ET' MOP ;
  103. MATORT = MAT1 ET MAT2 'ET' MAP ;
  104.  
  105. RIG12 = RIGI MODORT MATORT ;
  106.  
  107.  
  108. *______________________________________________________
  109. * CONDITIONS AUX LIMITES
  110. *______________________________________________________
  111.  
  112.  
  113. CDL1 = BLOQ 'UZ' pB ;
  114. CDL2 = BLOQ 'UT' pB ;
  115. CDL = CDL1 ET CDL2 ;
  116.  
  117. *______________________________________________________
  118. * CALCUL
  119. *______________________________________________________
  120. RIGITOT = RIG12 ET CDL ;
  121. FP = 'PRES' MOP MAP ;
  122. *
  123. DEP1 = RESO RIGITOT FP ;
  124.  
  125. DEFO0= DEFO CYL 0.0 DEP1 GRIS;
  126. DEFO1= DEFO CYL DEP1 (EXCO DEP1 'UR');
  127.  
  128. SI (NEG GRAPH 'N') ;
  129. TRAC (DEFO1 et DEFO0) ;
  130. FINSI;
  131.  
  132. * ref = calcul element LC8 tire de [DEMT 85-482]
  133. uref = (0.63615E-5 + 0.64702E-5) /2.;
  134.  
  135. * VALEUR pour le TEST
  136. ur1A = EXTR DEP1 'UR' PA ;
  137. ur1B = EXTR DEP1 'UR' PB ;
  138. MESS ' DEPLACEMENT RADIAL REFERENCE : ' uref;
  139. MESS ' DEPLACEMENT RADIAL EN PA CALCULE :' ur1A ur1B;
  140. RES1= ABS((ur1A - uref) / uref);
  141. MESS 'ECART RELATIF : ' RES1 ;
  142. MESS '************************************************';
  143. TEMPS ;
  144. *______________________________________________________
  145. * CODE BON FONCTIONNEMENT
  146. *______________________________________________________
  147.  
  148. SI (RES1 <EG 5.E-2);
  149. ERRE 0 ;
  150. SINO;
  151. ERRE 5 ;
  152. FINSI ;
  153. FIN;
  154.  
  155.  
  156.  

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