Télécharger fatigue.dgibi

Retour à la liste

Numérotation des lignes :

  1. opti dime 3 elem cub8 echo 0 ;
  2. *
  3. * operateur FATIGUE
  4. * exemples critere DANG VAN : JL FAYART
  5. *
  6. *
  7. LGRAPH = FAUX ;
  8.  
  9. p_ori = 0. 0. 0. ;
  10. e_x = 1. 0. 0. ; e_y = 0. 1. 0. ; e_z = 0. 0. 1. ;
  11.  
  12. p_A = p_ori plus p_ori ; p_B = p_A plus e_x ;
  13. l1 = d 1 p_A p_B ;
  14. s1 = l1 trans 1 e_y ;
  15. p_C = p_B plus e_y ; p_D = p_C moins e_x ;
  16. p_C = point s1 proc p_C ; p_D = point s1 proc p_D ;
  17. V1 = s1 volu trans 1 e_z ;
  18. p_E p_F p_G p_H = p_A p_B p_C p_D plus e_z ;
  19. p_E = point v1 proc p_E ;
  20. p_F = point v1 proc p_F ;
  21. p_G = point v1 proc p_G ;
  22. p_H = point v1 proc p_H ;
  23.  
  24. mo1 = mode v1 mecanique elastique ;
  25. ca1 = mate mo1 youn 2.e11 nu 0.3 ;
  26.  
  27. bl_A = bloq depla p_A ;
  28. bl_B = bloq uy uz p_B ;
  29. bl_C = bloq uz p_C ;
  30. bl_D = bloq ux uz p_D ;
  31. bl_t = bloq uz (p_E et p_F et p_G et p_H) ;
  32. bl_cis1 = bloq depla(p_A et p_B et p_E et p_F) ;
  33. bl_cis2 = bloq ux (p_C et p_D et p_G et p_H) ;
  34. bl_cis3 = bloq uy uz (p_C et p_D et p_G et p_H) ;
  35.  
  36. ***
  37. * traction
  38. ***
  39. ch_t = depi bl_t 6.5e-4 ;
  40. abs1 = prog 0. 1. 2. 3.;
  41. ev_t = evol manu abs1 (prog 0. -1. 0. 1.) ;
  42. ctrasi = char meca ev_t ch_t ;
  43.  
  44. trasi = table 'PASAPAS' ;
  45. trasi . modele = mo1 ;
  46. trasi . caracteristiques = ca1 ;
  47. trasi . blocages_mecaniques = bl_A et bl_B et bl_C et bl_D et bl_t ;
  48. trasi . chargement = ctrasi ;
  49. trasi . temps_calcules = abs1 ;
  50.  
  51. pasapas trasi ;
  52. *explorer trasi ;
  53. cafatp = manu chml trasi . modele 'AA' -0.33 'BB' 8.e7
  54. type 'CARACTERISTIQUES' stresses ;
  55. chasi1 = char trasi . temps trasi . contraintes ;
  56. chfa1 = FATI mo1 chasi1 cafatp ;
  57. *list chfa1 ;
  58. xtrasi = extr chfa1 dvkp 1 1 5 ;
  59. xtrasi0 = -8.75000E-03 ;
  60. er1 = abs((xtrasi - xtrasi0)/xtrasi0) ;
  61. SI (LGRAPH) ;
  62. trac chfa1 mo1 ;
  63. FINSI ;
  64. chfa1 = FATI mo1 chasi1 cafatp 'DVKP' ;
  65. rits16 = extr chfa1 dvkp 1 1 6 ;
  66. ets16 = extr chfa1 traj 1 1 6 ;
  67. SI (LGRAPH) ;
  68. dess ets16 ;
  69. FINSI ;
  70.  
  71. ***
  72. * traction repetee
  73. ***
  74. abs2 = prog 0. pas 1. 3. ;
  75. ev_r = evol manu abs2 (prog 0. 0.5 1. 0.5) ;
  76. ctrare = char meca ev_r ch_t ;
  77.  
  78. trare = table 'PASAPAS' ;
  79. trare . modele = mo1 ;
  80. trare . caracteristiques = ca1 ;
  81. trare . blocages_mecaniques = bl_A et bl_B et bl_C et bl_D et bl_t ;
  82. trare . chargement = ctrare ;
  83. trare . temps_calcules = abs2 ;
  84.  
  85. pasapas trare ;
  86. *explorer trare ;
  87. chasi2 = char trare . temps trare . contraintes ;
  88. chfa2 = FATI mo1 chasi2 cafatp 'DVKP' ;
  89. xtrare = extr chfa2 dvkp 1 1 5 ;
  90. xtrare0 = -0.41500 ;
  91. er2 = abs((xtrare - xtrare0)/xtrare0) ;
  92. SI (LGRAPH) ;
  93. trac chfa2 mo1 ;
  94. FINSI ;
  95. ritr16 = extr chfa2 dvkp 1 1 6 ;
  96. si (ritr16 > -0.2) ;
  97. etr16 = extr chfa2 traj 1 1 6 ;
  98. finsi ;
  99.  
  100.  
  101. ***
  102. * cisaillement
  103. ***
  104. ch_cis = depi bl_cis2 1.e-3;
  105. ccisa = char meca ev_t ch_cis ;
  106.  
  107. tcisa = table 'PASAPAS' ;
  108. tcisa . modele = mo1 ;
  109. tcisa . caracteristiques = ca1 ;
  110. tcisa . blocages_mecaniques = bl_cis1 et bl_cis2 et bl_cis3 ;
  111. tcisa . chargement = ccisa ;
  112. tcisa . temps_calcules = abs1 ;
  113.  
  114. pasapas tcisa ;
  115. *explorer tcisa ;
  116. chasi3 = char tcisa . temps tcisa . contraintes ;
  117. chfa3 = FATI mo1 chasi3 cafatp 'DVKP' ;
  118. rici16 = extr chfa3 dvkp 1 1 6 ;
  119. eci16 = extr chfa3 traj 1 1 6 ;
  120. xtcisa = extr chfa3 dvkp 1 1 5 ;
  121. xtcisa0 = -3.84615E-02 ;
  122. er3 = abs((xtcisa - xtcisa0)/xtcisa0) ;
  123. SI (LGRAPH) ;
  124. dess eci16 ;
  125. FINSI ;
  126.  
  127. si ((er1 < 0.05 ) et (er2 < 0.05) et (er3 < 0.05)) ;
  128. erre 0 ;
  129. sinon ;
  130. erre 5 ;
  131. finsi ;
  132.  
  133. fin ;
  134.  
  135.  
  136.  
  137.  
  138.  
  139.  
  140.  
  141.  

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