Télécharger t_visk2.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : t_visk2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. graph='N';
  5. saut page ;
  6. *------------------------
  7. * test loi visk2
  8. *------------------------
  9. opti dime 3 elem cub8 ;
  10. p_ori = 0. 0. 0.; e_x = 1. 0. 0. ; e_y = 0. 1. 0. ;
  11. e_z = 0. 0. 1. ;
  12. *
  13. p1 = p_ori ;
  14. l1 = p1 d 1 (p1 plus e_x) ;
  15. s1 = l1 trans 1 e_y ;
  16. v1 = s1 volu trans 1 e_z ;
  17. s2 = face 2 v1 ;
  18. *
  19. misopoin = vrai ;
  20. * ATTENTION le cas faux appelle un peu de travail de l operateur !
  21. *
  22. si faux ;
  23. mo1 = mode v1 mecanique elastique viscoplastique visk2 ;
  24. ma1 = mate mo1 young 2.e11 nu 0.3 rho 7.8e9 alpha 1.5e-5
  25. sigy 200.e6 h (0.1*2.e11) eta (2.*0.1*2.e11) hvis (0.1*2.e11) n 5. ;
  26. finsi ;
  27. si vrai ;
  28. mo1 = mode v1 mecanique elastique viscoplastique visk2 ;
  29. ma1 = mate mo1 young 150.e9 nu 0.3 rho 7.8e9 alpha 1.5e-5
  30. sigy 200.e6 h (0.1*2.e11) eta 0.2e11 hvis (2*0.2e11) n 1 ;
  31. ma2 = mate mo1 young 20.e9 nu 0.3 rho 7.8e9 alpha 1.5e-5
  32. sigy 75.0e6 h (0.05*2.e11) eta 0.2e11 hvis (2*0.2e11) n 1 ;
  33. finsi ;
  34. si faux ;
  35. mo1 = mode v1 mecanique elastique plastique cinematique ;
  36. ma1 = mate mo1 young 150.e9 nu 0.3 rho 7.8e9 alpha 1.5e-5
  37. sigy 200.e6 h (0.1*2.e11) ;
  38. ma2 = mate mo1 young 20.e9 nu 0.3 rho 7.8e9 alpha 1.5e-5
  39. sigy 75.0e6 h (0.05*2.e11) ;
  40. finsi ;
  41. *
  42. chsig1 = zero mo1 contrain ;
  43. chvar1 = zero mo1 varinter ;
  44. *
  45. t_ev = table ;
  46. t_prog = table ;
  47. si misopoin ;
  48. t_prog . 1 = (prog 0. pas 1. 8.)*1. ;
  49. t_prog . 2 = (prog 0. pas 1. 8.)*10. ;
  50. sinon ;
  51. t_prog . 1 = prog 0. 1. 2. 2.3 2.6 3. pas 1. 20. ;
  52. t_prog . 2 = prog 0. 10. 20. 23. 26. 30. pas 1. 40. pas 10. 60. ;
  53. finsi ;
  54. *
  55. *opti impr 'poub' ;
  56. *
  57. si misopoin ;
  58. n_tra = 1 ;
  59. sinon ;
  60. n_tra = 2 ;
  61. finsi ;
  62. repeter b_tra n_tra ;
  63.  
  64. bl1 = (bloq s1 uz) et (bloq p1 depla) et (bloq l1 uy) ;
  65. bl2 = bloq s2 uz ;
  66. chdimp1 = depi bl2 1.e-3 ;
  67. l_temps = t_prog . &b_tra ;
  68. si misopoin ;
  69. l_ordo = prog 0. 1. 2. ((dime l_temps) - 3)*3. ;
  70. l_ordo = l_ordo / 1. ;
  71. sinon ;
  72. l_ordo = prog 0. 1. 2. 2.3 2.6 ((dime l_temps) - 5)* 3. ;
  73. finsi ;
  74. ev_t = evol manu l_temps l_ordo ;
  75. cha1 = char dimp ev_t chdimp1 ;
  76. *
  77. tpas = table ;
  78. tpas . modele = mo1 ;
  79. tpas . caracteristiques = ma1 ;
  80. tpas . chargement = cha1 ;
  81. tpas . blocages_mecaniques = bl1 et bl2 ;
  82. tpas . temps_calcules = prog 0. pas 1. 5. ;
  83. *tpas . temps_calcules = prog 0. 1. 2.;
  84. *
  85. TMASAU=table;
  86. tpas . 'MES_SAUVEGARDES'=TMASAU;
  87. TMASAU .'DEFTO'=VRAI;
  88. TMASAU .'DEFIN'=VRAI;
  89. pasapas tpas ;
  90. *
  91. tpas . caracteristiques = ma2 ;
  92. tpas . temps_calcules = prog 6. pas 1. 8. ;
  93. pasapas tpas ;
  94. *
  95. * post-traitement
  96. o_res = prog ;
  97.  
  98. repeter b_evo (dime tpas . temps) ;
  99. rea1 = reac bl2 tpas . deplacements . (&b_evo-1) ;
  100. si (&b_evo ega 1) ;
  101. res1 = manu chpo 1 p1 fz 0. ;
  102. sinon ;
  103. res1 = resu rea1 ;
  104. finsi ;
  105. pe1 = point 1 (extr res1 mail) ;
  106. o_res = o_res et (prog (extr res1 fz pe1)) ;
  107. fin b_evo ;
  108.  
  109. t_ev . &b_tra = evol manu 'temps(s)' l_temps 'F(1.e3*T)' (o_res/1.e7) ;
  110.  
  111. fin b_tra ;
  112. *
  113. tabsymb = table ;
  114. tabsymb . titre = table ;
  115. tabsymb . 1 = 'MARQ CROI' ;
  116. tabsymb . titre . 1 = 'taux defo 0.001/s' ;
  117. tabsymb . 2 = 'MARQ LOSA' ;
  118. tabsymb . titre . 2 = 'taux defo 0.01/s' ;
  119. si (neg graph 'N') ;
  120. si misopoin ;
  121. dess (t_ev . 1) lege tabsymb
  122. titre ' visk2 / resultante / traction et relaxation' ;
  123. sinon ;
  124. dess (t_ev . 1 et t_ev . 2) lege tabsymb
  125. titre ' visk2 / resultante / traction et relaxation' ;
  126. finsi ;
  127. finsi ;
  128. *
  129. si misopoin ;
  130. si ( abs((extr o_res 6) - 2.41830E+08) > 1.e5) ;
  131. erre 5 ;
  132. sinon ;
  133. erre 0 ;
  134. finsi ;
  135. fin ;
  136.  
  137.  
  138.  
  139.  
  140.  
  141.  
  142.  

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