Télécharger la plaquette Cast3M


fr en

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 1   ;
  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 = 'PRES' 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 ;
  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 = 'BSIG' 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 '************************************************';
  144. TEMPS ;
  145. *______________________________________________________
  146. *                    CODE BON FONCTIONNEMENT
  147. *______________________________________________________
  148.  
  149. SI (RES1 <EG 5.E-2);
  150.    ERRE 0 ;
  151. SINO;
  152.    ERRE 5 ;
  153. FINSI ;
  154. FIN;
  155.  
  156.  
  157.  
  158.  
  159.  

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