Télécharger la plaquette Cast3M


fr en

Télécharger fatigue.dgibi

Retour à la liste

Numérotation des lignes : 

  1. opti dime 3 elem cub8 echo 1 ;
  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.  
  54. * kich : coefficients arbitraires
  55. cafatp = manu chml trasi . modele  'ADVK' -0.33 'BDVK' 8.e7
  56.   'APAP' -0.34 'BPAP' 8.1e7 'ASIN' -0.35 'BSIN' 8.2e7 
  57.   'ACRO' -0.36 'BCRO' 8.3e7 'A_DC' -0.37 'B_DC' 8.4e7   
  58.     type 'CARACTERISTIQUES' stresses ;
  59. chasi1 = char  trasi . temps trasi . contraintes ;
  60. chfa1 = FATI mo1 chasi1 cafatp ;
  61. *list chfa1 ;
  62.  xtrasi = extr chfa1 dvkp 1 1 5 ; 
  63.  xtrasi0 = -8.75000E-03 ;
  64.  er1 = abs((xtrasi - xtrasi0)/xtrasi0) ;
  65. SI (LGRAPH) ;
  66.   trac chfa1 mo1 ;
  67. FINSI ;
  68. chfa1 = FATI mo1 chasi1 cafatp 'DVKP' ;
  69. rits16 = extr chfa1 dvkp 1 1 6 ;
  70. ets16 = extr chfa1 'PTAU' 1 1 6 ;
  71. SI (LGRAPH) ;
  72.  dess ets16 ;
  73. FINSI ;
  74. chfa1 = FATI mo1 chasi1 cafatp 'PAPA' 'SEUI' 'TOUS' ;
  75. rits16 = extr chfa1 papa 1 1 6 ;
  76. ets16 = extr chfa1 'PTAU' 1 1 6 ;
  77. SI (LGRAPH) ;
  78.  dess ets16 ;
  79. FINSI ;
  80. ***
  81. * traction repetee
  82. ***
  83. abs2 = prog 0. pas 1. 3. ; 
  84. ev_r = evol manu abs2 (prog 0. 0.5 1. 0.5) ;
  85. ctrare = char meca ev_r ch_t ;
  86.  
  87. trare = table 'PASAPAS' ;
  88. trare . modele = mo1 ;
  89. trare . caracteristiques = ca1 ;
  90. trare . blocages_mecaniques = bl_A et bl_B et bl_C et bl_D et bl_t ;
  91. trare . chargement = ctrare ;
  92. trare . temps_calcules = abs2 ;
  93.  
  94. pasapas trare ;
  95. *explorer trare ;
  96. chasi2 = char  trare . temps trare . contraintes ;
  97. chfa2 = FATI mo1 chasi2 cafatp 'DVKP' ;
  98.  xtrare = extr chfa2 dvkp 1 1 5 ; 
  99.  xtrare0 = -0.41500 ;
  100.  er2 = abs((xtrare - xtrare0)/xtrare0) ;
  101. SI (LGRAPH) ;
  102.   trac chfa2 mo1 ;
  103. FINSI ;
  104. ritr16 = extr chfa2 dvkp 1 1 6 ;
  105. si (ritr16 > -0.2) ;
  106.   etr16 = extr chfa2 ptau 1 1 6 ;
  107. finsi ;
  108.  
  109.  
  110. ***
  111. * cisaillement
  112. ***
  113. ch_cis = depi bl_cis2 1.e-3; 
  114. ccisa = char meca ev_t ch_cis ;
  115.  
  116. tcisa =  table 'PASAPAS' ;
  117. tcisa . modele = mo1 ;
  118. tcisa . caracteristiques = ca1 ;
  119. tcisa . blocages_mecaniques = bl_cis1 et bl_cis2 et bl_cis3 ;
  120. tcisa . chargement = ccisa ;
  121. tcisa . temps_calcules = abs1 ;
  122.  
  123. pasapas tcisa ;
  124. *explorer tcisa ;
  125. chasi3 = char  tcisa . temps tcisa . contraintes ;
  126. chfa3 = FATI mo1 chasi3 cafatp 'DVKP' ;
  127. rici16 = extr chfa3 dvkp 1 1 6 ;
  128. eci16 = extr chfa3 'PTAU' 1 1 6 ;
  129.  xtcisa = extr chfa3 dvkp 1 1 5 ; 
  130.  xtcisa0 = -3.84615E-02 ;
  131.  er3 = abs((xtcisa - xtcisa0)/xtcisa0) ;
  132. SI (LGRAPH) ;
  133.   dess eci16 ;
  134. FINSI ;
  135.  
  136. si ((er1 < 0.05 ) et (er2 < 0.05) et (er3 < 0.05)) ;
  137.      erre 0 ;
  138. sinon ;
  139.     erre 5 ;
  140. finsi ;
  141. fin ;
  142.  
  143.  
  144.  
  145.  
  146.  
  147.  
  148.  
  149.  
  150.  
  151.  

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