Télécharger la plaquette Cast3M


fr en

Télécharger rupt12.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  rupt12.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. OPTION ECHO 1 ;
  5.  complet = faux;
  6. * complet = vrai;
  7. *
  8. *************************************************
  9. *  pour calcul complet mettre complet à : vrai;
  10. ***************************************************
  11. **************************************************
  12. *                                                I
  13. *    VALIDATION DE LA PROCEDURE                  I
  14. *        G_THETA EN DYNAMIQUE.                   I
  15. *                                                I
  16. *     PLAQUE EN TRACTION PURE AVEC               I
  17. *          CHARGEMENT 'f(t)'                     I
  18. *                                                I
  19. *----------------------------------------------- I
  20. * GEOMETRIE :                                    I
  21. *   Longueur de la plaque   :  40  mm            I
  22. *   Largeur de la plaque    :  104  mm           I 
  23. *   Longueur de la fissure  :  24 mm             I
  24. *                                                I
  25. * MATERIAU :                                     I
  26. *   Module d'Young          : 7.561 daN mm^(-2)  I
  27. *   Coefficient de poisson  : 0,286              I
  28. *   Masse volumique         : 2.45E-6 kg/mm^(-3) I
  29. *                                                I
  30. **************************************************
  31. OPTI echo 1 ;
  32. OPTION DIME 2 ELEM TRI6 MODE PLAN DEFORMATION ;
  33. graph='N';
  34. *---------------------------------------*
  35. * CONSTRUCTION DU MAILLAGE (RAYONANT)   *
  36. *---------------------------------------*
  37.  A =  0. 0.    ;  B = 52. 0.   ;
  38.  C = 52. 40.   ;  D = 0. 40.   ;
  39.  P1 = 0. 20.   ;  P2 = 0. 20.  ;
  40.  P3 = 12.  20. ;  P5 = 52. 20. ;
  41.  P4 = 6.  20. ;
  42. *
  43. SS COUR1 l1 l2 l3 l4 LVD LVG = @rayo P3 P4 6 ;
  44.   P4 = L1 POINT INITIAL ;
  45.   P6 = L4 POIN 'INITIAL' ;
  46.   P7 = L3 POIN 'INITIAL' ;
  47.   P7 = L2 POIN 'INITIAL' ;
  48. *1ER CONTOUR
  49.   AB    = DROI  A  B  'DINI' 5. 'DFINI' 5. ;
  50.   BP5   = DROI  B  P5 'DINI' 5. 'DFINI' 2. ;
  51.   P5P7  = DROI  P5 P7 'DINI' 2. 'DFINI' 1. ;
  52.   P4P2  = DROI  P4 P2 'DINI' 2. 'DFINI' 1. ;
  53.   P2A   = DROI  P2 A  'DINI' 1. 'DFINI' 5. ;
  54. *2EM CONTOUR
  55.   P1D   = DROI  P1 D  'DINI' 1. 'DFINI' 5. ;
  56.   DC    = DROI  D  C  'DINI' 5. 'DFINI' 5. ;
  57.   CP5   = DROI  C  P5 'DINI' 5. 'DFINI' 2. ;
  58.   P6P1  = DROI  P6 P1 'DINI' 1. 'DFINI' 1. ;
  59. *-----------------------------------------------
  60.   CT2 = P4P2 ET P2A ET AB ET BP5 ET P5P7 ET L2 ET L1 ;
  61.   CT3 = P6P1 ET  P1D ET DC ET CP5 ET P5P7 ET L3 ET L4 ;
  62.   MAIL1 = SURF CT2 PLANE ;
  63.   MAIL2 = SURF CT3 PLANE ;
  64.   MAIL = MAIL1 ET MAIL2 ET SS ;
  65. SI ( ega graph 'O');
  66.   TRAC MAIL ;
  67. finsi;
  68. *      LEVRES SUP ET INF DE LA FISSURE
  69.    LSUP =  LVD ;
  70.    LINF =  LVG ;
  71. *---------------------------------------*
  72. * CONSTRUCTION DU MODEL
  73. *---------------------------------------*
  74.      YOU1 = 7.561E7;
  75.      NU1  = 0.286 ;
  76.      RH1  = 2.45E-6;
  77. MO=MODE MAIL MECANIQUE ELASTIQUE ISOTROPE;
  78. MAT = MATE MO YOUN YOU1 NU NU1 RHO RH1;
  79. *---------------------------------------*
  80. * RIGIDITE 
  81. *---------------------------------------*
  82.      RIG1 = RIGIDITE MO MAT;
  83. *---------------------------------------*
  84. * CONDITIONS AUX LIMITES ;
  85. *---------------------------------------*
  86.   cl1 = bloque UX P1D ;
  87.   cl2 = bloque UX P2A ;
  88.   cl4 = bloque UX CP5 ; 
  89.   cl5 = bloque UX BP5 ;
  90.   CL3 = BLOQUE UY  P5 ;
  91.   CLT = CL1 ET CL2 ET CL3 ET CL4 ET CL5 ;
  92. *---------------------------------------*
  93. *                                       *
  94. *    Procedure PASAPAS DYNAMIQUE        *
  95. *                                       *
  96. *---------------------------------------*
  97. * Matrices de rigidite et de masse
  98. RIGT = RIG1 ET CLT;
  99. MAST = MASS MO MAT ;
  100. *---------------------------------------*
  101. *     Instant de calcul
  102. *
  103.   TI = 0. ;
  104. si complet;
  105.  TF = 20E-6; NP1 = 40 ;
  106. sinon;
  107. tf = ( 20e-6 / 40 ) * 6 ; np1= 6;
  108.  finsi;
  109.   DT  = (TF - TI ) / NP1 ;
  110. *---------------------------------------*
  111.   TITRE ' CHARGEMENT : IMPACT DE 0.4E3 ';
  112. *---------------------------------------*
  113.    F0 =  0.4E6  ;
  114. L_TEMPS  = PROG 0.0 'PAS' DT 'NPAS' (NP1+1);
  115. L_TEMP2  = PROG 0.0 'PAS' DT 'NPAS' NP1 ;
  116. F_TEMPS  = PROG  ( NP1 + 2 )*F0 ;
  117. *
  118. EVO_FOR = EVOL MANU 'temps' L_TEMPS 'f(t)' F_TEMPS ;
  119.  CH1_FORC = PRES 'MASS' MO -1. AB ;
  120.  CH3_FORC = PRES 'MASS' MO -1. DC ;
  121.  CHP_FORC = CH1_FORC ET CH3_FORC ;
  122.  CHARG_T  = CHARG 'MECA' CHP_FORC  EVO_FOR;
  123. *
  124.   FR_COUP = 0.25 /DT  ;
  125.   TAB_DYN = TABLE ;
  126.   TAB_DYN.DYNAMIQUE = VRAI;
  127.   TAB_DYN.MODELE = MO;
  128.   TAB_DYN.CARACTERISTIQUES = MAT;
  129.   TAB_DYN.BLOCAGES_MECANIQUES = CLT;
  130.   TAB_DYN.CHARGEMENT = CHARG_T;
  131.   TAB_DYN.TEMPS_CALCULES = L_TEMP2;
  132.   PASAPAS TAB_DYN;
  133.  
  134. **************** tracés  *******************
  135. SI ( ega graph 'O');
  136.   trac (vect CHP_FORC 'FORC' 'BLEU') MAIL;
  137.   dess EVO_FOR;
  138.   isosig2 = prog 0. PAS 0.25E6 5E6;
  139.   ntemp2 = dime L_TEMP2;
  140.   itemp2 = -1;
  141.   repe BDEFO ntemp2;
  142.     itemp2 = itemp2 + 1;
  143.     temp2 = extr L_TEMP2 (itemp2 + 1);
  144.     u2t = TAB_DYN . DEPLACEMENTS . itemp2;
  145.     sig2t = EXCO 'SMYY' (TAB_DYN . CONTRAINTES . itemp2);
  146.     trac sig2t MO (DEFO u2t MAIL 10) isosig2
  147.         'TITR' (chai 'u(t=' temp2 's') 'NCLK';
  148.     xpause = 0.;
  149.     repe bpause 1000;
  150.        xpause = xpause + 1.;
  151.     fin bpause;
  152.   fin BDEFO;
  153. finsi;
  154.  
  155. *opti donn 5 ;
  156. *=======================================*
  157. *
  158. * CALCUL DE G PAR LA METHODE G_THETA
  159. *
  160. *=======================================*
  161.    SUPTAB = TABL ;
  162.    SUPTAB.'OBJECTIF' = MOT 'J_DYNA';
  163.    SUPTAB.'LEVRE_SUPERIEURE' = LSUP;
  164.    SUPTAB.'LEVRE_INFERIEURE' = LINF;
  165.    SUPTAB.'FRONT_FISSURE' = P3 ; 
  166.    SUPTAB.'SOLUTION_PASAPAS' = TAB_DYN;
  167.    SUPTAB.'COUCHE' = 3;                 
  168. *----------------------------------*
  169. * Appel DE la procedure G_THETA
  170. *----------------------------------*
  171.    G_THETA SUPTAB;
  172. *------------------------------------*
  173. * Evolution de G en fonction du temps.
  174. * -----------------------------------*
  175. SI ( ega graph 'O');
  176.    DESS SUPTAB . 'EVOLUTION_RESULTATS' ;
  177. finsi;
  178.  
  179. *------------------------------*
  180. *     Calcul de K
  181. *------------------------------*
  182. ab1 = EXTR (SUPTAB.EVOLUTION_RESULTATS) ABSC;
  183. ORD1 = EXTR (SUPTAB.EVOLUTION_RESULTATS) ORDO;
  184. ORD2 = ((ORD1 * YOU1) /(1 - (0.286**2)))**0.5 ; 
  185. K1 = EVOL MANU 'TEMPS' AB1 'K1(t)' ORD2;
  186. SI ( ega graph 'O');
  187.   DESS K1 ;
  188. finsi;
  189.  
  190.  
  191. **************************************
  192. ****  SOLUTION ANALYTIQUE ************
  193. **************************************
  194. *anciennes valeurs de reference (formulation Attigui)
  195. *li2 = prog 82275 2701582 4420771  5525946 ;
  196. *nouvelles valeurs de reference (formulation B Prabel)
  197. KREF = prog 74056. 3.12957E+06  4.34363E+06  5.64263E+06;
  198.  
  199. ***RECUPERATION DU RESULTAT ET COMPARAISONS***
  200. ERR1MAX = prog;
  201. si complet; nfois = 4; 
  202. sinon; nfois = 1;
  203. finsi;
  204. I = 1; 
  205. REPETER BOU1 nfois;
  206.   I = I + 5;
  207.   J = &BOU1;
  208.   K1REF = extr KREF J ;
  209.   K1CAL = extr ord2 I;
  210.   ERR1 = abs ((K1CAL - K1REF)/K1REF);
  211.   MESS (CHAI 'I=' I ' K1REF=' KEXP ' K1CAL='K_CAL 
  212.              '  Erreur relative =' ERR1); 
  213.   ERR1MAX = ERR1MAX et  ERR1;
  214. FIN BOU1; 
  215. *LIST ERR1MAX;
  216.  
  217. ERR1 = maxi ERR1MAX;
  218. *SI (ERR1 <  5.E-2);
  219. SI (ERR1 <  1.E-2);
  220.       ERRE 0;
  221. SINO;
  222.       ERRE 5;
  223. FINSI;
  224. FIN;
  225.  
  226.  
  227.  
  228.  
  229.  
  230.  
  231.  
  232.  
  233.  
  234.  
  235.  
  236.  
  237.  
  238.  
  239.  
  240.  
  241.  
  242.  
  243.  

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