Télécharger test_sif_3d.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : test_sif_3d.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5.  
  6. ************************************************************************
  7. * Verification & Validation de la procedure SIF
  8. ************************************************************************
  9. * Geometrie : Piece et fissure
  10. * -----------
  11. * 3D, milieu infini (ici cylindre de rayon l1 et hauteur 2*l1) avec
  12. * fissure en forme de disque de rayon l2 (penny-shaped crack)
  13. * soumise a une contrainte sig, dans la direction orthogonale a la fissure
  14. *
  15. * Validation : Test sur le calcul du KI
  16. * -----------
  17. * Valeur theorique KI = 2 * sig * (l2 / pi)**0.5
  18. *
  19. ************************************************************************
  20.  
  21. * Options generales
  22. OPTI 'DIME' 3 'ELEM' 'CU20' ;
  23.  
  24. * Parametres
  25. l1 = 10. ;
  26. l2 = 1. ;
  27. sig = 42.E6 ;
  28.  
  29. * Maillage
  30. * --densites (loin et pres de la fissure)
  31. den1 = l1 / 4. ;
  32. den2 = l2 / 10. ;
  33. * --fissure levre superieure
  34. OPTI 'DENS' den2 ;
  35. p1 = 0. 0. 0. ;
  36. p2 = l2 0. 0. ;
  37. l12 = DROI p1 p2 ;
  38. s1 = l12 ROTA 90 p1 (p1 PLUS (0. 0. 1.)) ;
  39. lff = s1 COTE 2 ;
  40. ELIM s1 1.E-10 ;
  41. s1s = REGE s1 ;
  42. s1i = s1s PLUS (0. 0. 0.) ;
  43. ELIM s1i lff 1.E-10 ;
  44. * --section non fisuree
  45. OPTI 'DENS' den1 ;
  46. p3 = l1 0. 0. ;
  47. p4 = 0. l1 0. ;
  48. l34 = CERC (NBEL lff) p3 p1 p4 ;
  49. s2 = REGL lff l34 ;
  50. * --volumes superieurs et inferieurs
  51. v1 = VOLU 'TRAN' (s1s ET s2) 'DINI' den2 'DFIN' den2 (0. 0. (0.05 * l1 ));
  52. v1 = VOLU 'TRAN' v1 'DINI' den2 'DFIN' den1 (0. 0. (0.95 * l1 ));
  53. v2 = VOLU 'TRAN' (s1i ET s2) 'DINI' den2 'DFIN' den2 (0. 0. (-0.05 * l1));
  54. v2 = VOLU 'TRAN' v2 'DINI' den2 'DFIN' den1 (0. 0. (-0.95 * l1));
  55. vtot = v1 ET v2 ;
  56.  
  57. * Modele et materiau
  58. mo = MODE vtot 'MECANIQUE' ;
  59. ma = MATE mo 'YOUN' 200.E9 'NU' 0.3 ;
  60.  
  61. * Blocages
  62. x y z = COOR vtot ;
  63. px0 = x POIN 'MINI' ;
  64. py0 = y POIN 'MINI' ;
  65. bl1 = BLOQ 'UX' px0 ;
  66. bl2 = BLOQ 'UY' py0 ;
  67. bl3 = BLOQ 'UZ' l34 ;
  68.  
  69. * Chargement
  70. shaut = v1 FACE 2 ;
  71. sbas = v2 FACE 2 ;
  72. f1 = PRES 'MASS' mo shaut (-1. * sig) ;
  73. f2 = PRES 'MASS' mo sbas (-1. * sig) ;
  74. f = f1 ET f2 ;
  75.  
  76. * Resolution
  77. r0 = RIGI mo ma ;
  78. rig = r0 ET bl1 ET bl2 ET bl3 ;
  79. u = RESO rig f ;
  80.  
  81. * Solution de reference
  82. k1ref = 2. * sig * ((l2 / pi) ** 0.5) ;
  83.  
  84. * Solution par SIF
  85. t1 = TABL ;
  86. t1 . 'MAILLAGE' = vtot ;
  87. t1 . 'FRTFISS' = lff ;
  88. t1 . 'LEVRE_1' = s1s ;
  89. t1 . 'MODMIXTE' = VRAI ;
  90. t1 . 'LEVRE_2' = s1i ;
  91. SIF ma u t1 ;
  92. k1sif = t1 . 'K1' . p2 ;
  93. errsif = 100. * (k1sif - k1ref) / k1ref ;
  94.  
  95. * Affichage
  96. OPTI 'ECHO' 0 ;
  97. SAUT 3 'LIGNE' ;
  98. MESS 'BILAN, CALCUL DE KI' ;
  99. MESS 'Solution de reference ' k1ref ;
  100. MESS 'Calcul avec SIF ' k1sif ;
  101. MESS 'Ecart relatif (%) ' errsif ;
  102.  
  103. * Test d'erreur
  104. SI ((ABS errsif) > 1.D0) ;
  105. ERRE 'ERREUR DANS LE CALCUL DU KI' ;
  106. SINON ;
  107. MESS 'CAS TEST PASSE AVEC SUCCES !' ;
  108. FINSI ;
  109. OPTI 'ECHO' 1 ;
  110.  
  111. FIN ;
  112.  
  113.  
  114.  

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