Télécharger rupt16-weib.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : rupt16-weib.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *----------------------------------------------------------------------
  5. *----------------------------------------------------------------
  6. * *
  7. * Test du critère de Weibull pour un cylindre *
  8. * en traction modelisé en axisymétrique *
  9. * sigu constante ou évolution constante *
  10. * *
  11. *----------------------------------------------------------------
  12. *----------------------------------------------------------------
  13. * *
  14. opti echo 1 ;
  15. * paramètres pour calculer la contrainte de Weibull *
  16. m = 22. ; v0 = 0.000125 ; ic = 1. ; n = 2. ;
  17. *sigu = cste *
  18. sigu1 = 2530. ;
  19. *sigu = evol_cste *
  20. lte1 = prog -1000. 1000. ; lsigu = prog 2530. 2530. ;
  21. sigu2 = evol manu 'temper' lte1 'sigu' lsigu ;
  22. * *
  23. *****************************************************************
  24. * *
  25. * MODELE AXISYMETRIQUE *
  26. * *
  27. *****************************************************************
  28. *
  29. *****************************************************************
  30. * Maillage *
  31. *****************************************************************
  32. opti elem qua8 dime 2 mode axis ;
  33. * *
  34. * densites utilisees *
  35. l1 = 0.25 ;d1 = 50e-3 ; d2 = 120e-3 ; d11 = d1/10. ;
  36. * *
  37. * coordonnees des points *
  38. po = 0. 0. ; pb = l1 0. ; pc = 0. d1 ; pd = 0. l1 ;
  39. pa = d1 0. ;
  40. * *
  41. * definition des lignes *
  42. ligoa = droi po pa 1 ;
  43. ligab = droi pa pb 'DINI' 0.08 'DFIN' 0.12 ;
  44. ligob = ligoa et ligab ;
  45. surf1 = ligob tran 1 pc ;
  46. *
  47. lig1 = cote 3 surf1 ;
  48. * *
  49. * création de la surface *
  50. surf3 = lig1 tran 'DINI' 0.08 'DFIN' 0.12 (0. (l1-d1)) ;
  51. surf2 = surf1 et surf3 ;
  52. * *
  53. * extraction des lignes *
  54. lib = cote 1 surf2 ;
  55. lic = cote 3 surf2 ;
  56. * *
  57. elim surf2 1e-3 ;
  58. * *
  59. *****************************************************************
  60. * création du modéle *
  61. *****************************************************************
  62. mod1 = mode surf2 mecanique elastique plastique ;
  63. * *
  64. *****************************************************************
  65. * définition du matériau *
  66. *****************************************************************
  67. * *
  68. * courbe de traction *
  69. peps = prog 0. 0.0026 0.0027 0.0030 0.0035 0.0045 0.0065 0.0085 0.0105 0.0125 0.0176 0.0227 0.0329 ;
  70. psig = prog 0. 503. 503. 503. 503. 503. 503. 503. 503. 516. 538. 557. 594. ;
  71. tracevt = evol manu 'eps' peps 'sig' psig ;
  72. * *
  73. mat1 = mate mod1 youn (503./0.0026) nu 0.3 tract tracevt alpha 1.e-5 ;
  74. * *
  75. *****************************************************************
  76. * rigidité *
  77. *****************************************************************
  78. cdl0 = bloq uz lib ;
  79. cdl2 = bloq uz lic ;
  80. dep1 = depi cdl2 4e-3 ;
  81. forcevt = evol manu abs (prog 0. 1.) ord (prog 0. 1.) ;
  82. char1 = char 'DIMP' dep1 forcevt ;
  83. * *
  84. tabtps = table ;
  85. tabtps. 0 = 0. ;
  86. tabtps. 1 = 1. ;
  87.  
  88. tabT = table ;
  89. tabT. 0 = manu 'CHPO' surf2 1 'T' 100. ;
  90. tabT. 1 = manu 'CHPO' surf2 1 'T' -100. ;
  91.  
  92. char2 = CHAR 'T' tabtps tabT ;
  93.  
  94. *****************************************************************
  95. * definition de la table d'entree *
  96. *****************************************************************
  97. tab1 = table ;
  98. tab1.caracteristiques = mat1 ;
  99. tab1.modele = mod1 ;
  100. tab1.blocages_mecaniques = cdl0 et cdl2 ;
  101. tab1.chargement = char1 et char2 ;
  102. tab1.temps_calcules = (prog 0. pas 0.5 1.) ;
  103. tab1.talpha_reference = 20. ;
  104. * *
  105. *****************************************************************
  106. * procedure de calcul *
  107. *****************************************************************
  108. TMASAU=table;
  109. tab1 . 'MES_SAUVEGARDES'=TMASAU;
  110. TMASAU .'DEFTO'=VRAI;
  111. TMASAU .'DEFIN'=VRAI;
  112. pasapas tab1 ;
  113. * *
  114. *****************************************************************
  115. * Calcul de la contrainte de Weibull axi *
  116. *****************************************************************
  117. *sigu est une constante
  118. * *
  119. * zone du critere *
  120. surfcrit = elem surf2 1 ;
  121. zone1 = redu mod1 surfcrit;
  122. * *
  123. enttab1 = table ;
  124. enttab1.'XMULT' = 2. ;
  125. enttab1.'OBJMO' = zone1 ;
  126. enttab1.'WEIBULL' = vrai ;
  127. enttab1.'TNONL' = tab1 ;
  128. enttab1.'M' = m ; enttab1.'V0' = v0 ; enttab1.'SIGU' = sigu1 ;
  129. enttab1.'IC' = ic ; enttab1.'N' = n ;
  130. enttab1.'TEMPER'= 0 ;
  131. * *
  132. sortab1 = critloc enttab1;
  133. * *
  134. sigw1 = (sortab1.'SIGW').1. ;
  135. * *
  136. *****************************************************************
  137. * Calcul de la contrainte de Weibull axi *
  138. *****************************************************************
  139. *sigu est une evolution constante
  140. * *
  141. * zone du critere *
  142. surfcrit = elem surf2 1 ;
  143. zone1 = redu mod1 surfcrit;
  144. * *
  145. enttab2 = table ;
  146. enttab2.'XMULT' = 2. ;
  147. enttab2.'OBJMO' = zone1 ;
  148. enttab2.'WEIBULL' = vrai ;
  149. enttab2.'TNONL' = tab1 ;
  150. enttab2.'M' = m ; enttab2.'V0' = v0 ; enttab2.'SIGU' = sigu2 ;
  151. enttab2.'IC' = ic ; enttab2.'N' = n ;
  152. enttab2.'TEMPER'= 1 ;
  153. enttab2.'TEREF'= 150. ;
  154. * *
  155. sortab2 = critloc enttab2;
  156. * *
  157. sigw2 = (sortab2.'SIGW').1. ;
  158. * *
  159. *****************************************************************
  160. * Test d'erreur *
  161. *****************************************************************
  162. ecart = ( abs(sigw1-sigw2) ) / sigw1 ;
  163. mess 'ecart relatif : ' (ecart*100.) '%' ;
  164. si (ecart < 0.1 ) ;
  165. erre 0 ;
  166. sinon ;
  167. erre 5 ;
  168. finsi ;
  169. * *
  170. fin ;
  171.  
  172.  
  173.  
  174.  
  175.  
  176.  

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