Télécharger la plaquette Cast3M


fr en

Télécharger rupt3.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  rupt3.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *           Test Rupt3.dgibi: Jeux de données         *
  5. *           ---------------------------------         *
  6. *                                                     *
  7. OPTION echo 1                                         ;
  8. GRAPH = 'N'                                           ;
  9. SAUT PAGE                                             ;
  10. *                                                     ;
  11. ******************************************************;
  12. *                                                     ;
  13. *                                                     ;
  14. *              QUALIFICATION DU CALCUL DE G           ;
  15. *              EN ELASTICITE LINEAIRE SUR             ;
  16. *              UNE PLAQUE A FISSURE INTERNE           ;
  17. *                                                     ;
  18. *                                                     ;
  19. * le calcul est compare a celui obtenu par            ;
  20. * ROOKE et CARTWRIGHT                                 ;
  21. * dans Compendium of Stress Intensity Factors         ;
  22. *                                                     ;
  23. ******************************************************;
  24. *                                                     ;
  25. OPTION DIME 2 ELEM QUA8 MODE PLAN DEFO                ;
  26. *                                                     ;
  27. *----------------- DEFINITION DU MAILLAGE ------------;
  28. *                                                     ;
  29. a = 100. ; b = 200. ; h = 600. ;t = a / 100.          ;
  30. densite  t ; pf = (a 0.)                              ;
  31. c1 = (c ( pf moin (t 0.)) pf ( pf plus (0. t)))
  32.       c pf ( pf plus (t 0.))                          ;
  33. sf = cout pf c1;
  34. r1 = t ; rr1 = t                                      ;
  35. repeter bhomo 7;
  36.    ri  =  r1 + ( 0.3 * r1 )                           ;
  37.    rri = rr1 + ri                                     ;
  38.    dens ri                                            ;
  39.    ci =( c ( pf moin (rri 0.)) pf ( pf plus (0. rri)))  
  40.          c pf ( pf plus (rri 0.))                     ;
  41.    sf = sf et ( cout c1 ci )                          ;
  42.    c1 = ci ; r1 = ri ; rr1 = rri                      ;
  43. fin bhomo                                             ;
  44. dens (a / 3.)                                         ;
  45. p0  = (0. 0.) ; p1 = (b 0.)                           ;
  46. p2  = p0 plus (0. a) ; p3 = p1 plus (0. a)            ;
  47. pi1 = ci poin 1  ; l1  = pi1 d p0                     ;
  48. n = (nbel l1) * -1                                    ;
  49. pi2 = ci poin 4  ; l2  = pi2 d n p2                   ;
  50. pi3 = ci poin 10 ; l3  = pi3 d n p3                   ;
  51. pi4 = ci poin 13 ; l4  = pi4 d   p1                   ;
  52. ci  = inve ci    ;ligh = p2  d   p3                   ;
  53. sc1 = dall l1 (p0 d p2) (inve l2) (ci ELEM comp pi2 pi1)   ;
  54. sc2 = dall l2 ligh (inve l3) (ci ELEM comp pi3 pi2)        ;
  55. sc3 = dall l3 (p3 d p1) (inve l4) (ci ELEM comp pi4 pi3)   ;
  56. sc  = sc1 et sc2 et sc3                               ;
  57. dens (a / 2.)                                         ;
  58. mrest  = ligh tran (0. (h - a)) dini (a / 2.)         ;
  59. lihaut = mrest cote 3                                 ;
  60. phd    = lihaut poin init                             ;
  61. ccp    = sf et sc et mrest                            ;
  62. elim  ccp  0.001                                      ;
  63. cccp   = cont ccp                                     ;
  64. lifis  = cccp ELEM comp pf p0                              ;
  65. libas  = cccp ELEM comp p1 pf                              ;
  66. licot  = ccp  poin droi (0. 0.) (0. 100.) 0.001       ;
  67. licot  = cccp elem appuy stric licot                  ;
  68. SI ( NEG GRAPH 'N' )                                  ;
  69.      TRAC  CCP                                        ;
  70. FINSI                                                 ;
  71. *                                                     ;
  72. *------------------- CREATION DU MODELE --------------;
  73. *                                                     ;
  74. objaf = MODE ccp mecanique elastique isotrope         ;
  75. *                                                     ;
  76. *----------- DEFINITION DES CARACTERISTIQUES ---------;
  77. *-----------   MATERIELLES ET GEOMETRIQUES   ---------;
  78. *                                                     ;
  79. mat = MATE objaf  YOUN 2.e4  NU 0.3                   ;
  80. *                                                     ;
  81. *-------- CALCUL DES RIGIDITES ELEMENTAIRES ----------;
  82. *--------     ET DEFINITION DES BLOCAGES    ----------;
  83. *                                                     ;
  84. rig  = rigi objaf mat                                 ;
  85. cdl1 = bloq Ux licot                                  ;
  86. cdl2 = bloq Uy libas                                  ;
  87. *                                                     ;
  88. *-------------- DEFINITION DU CHARGEMENT -------------;
  89. *                                                     ;
  90. F = pres mass objaf -1. lihaut                        ;
  91. *                                                     ;
  92. *-------- RESOLUTION ET CALCUL DES CONTRAINTES -------;
  93. *                                                     ;
  94. U    = reso (rig et cdl1 et cdl2 ) F                  ;
  95. SIG  = sigm mat objaf U                               ;
  96. *                                                     ;
  97. *------ INITIALISATION DE LA TABLE EN ENTREE ---------;
  98. *-----       DE LA PROCEDURE G_THETA         ---------;
  99. *                                                     ;
  100. SUPTAB = TABLE ;
  101. SUPTAB.'OBJECTIF' = MOT 'J';
  102. SUPTAB.'LEVRE_SUPERIEURE' = lifis;
  103. SUPTAB.'FRONT_FISSURE' = PF ; 
  104. SUPTAB.'MODELE' = objaf;
  105. SUPTAB.'CARACTERISTIQUES' = mat;
  106. SUPTAB.'SOLUTION_RESO' = u;
  107. SUPTAB.'CHARGEMENTS_MECANIQUES' = F;
  108. SUPTAB.'COUCHE' = 4;                                     
  109. *                                                     ;
  110. *--------------- APPEL A LA PROCEDURE G_THETA --------;
  111. *                                                     ;
  112. SAUT PAGE                                             ;
  113. G_THETA suptab                                        ;
  114. *                                                     ;
  115. *----- RECUPERATION DU RESULTAT ET COMPARAISONS ------;
  116. *                                                     ;
  117. G   = suptab.'RESULTATS'                              ;
  118. Kth = 20.8426                                         ;
  119. K = (20000. * G / (1 - (0.3**2)))**0.5                ;
  120. mess '  '                                             ;
  121. mess '  '                                             ;
  122. mess 'taux de restitution energetique  G :' G         ;
  123. mess 'solution calculee             Kcal :' K         ;
  124. mess 'solution analytique           Kthe :' Kth       ;
  125. *                   CODE FONCTIONNEMENT               ;
  126. RESI = abs (( k - Kth ) / Kth )                       ;
  127. SI (RESI <  5E-2)                                     ;
  128.     ERRE  0                                           ;
  129. SINO                                                  ;
  130.     ERRE  5                                           ;
  131. FINSI                                                 ;
  132. fin;
  133.  
  134.  
  135.  
  136.  
  137.  
  138.  
  139.  
  140.  
  141.  
  142.  
  143.  
  144.  
  145.  
  146.  
  147.  
  148.  
  149.  

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