Télécharger TirantLAB2.dgibi

Retour à la liste

Numérotation des lignes :

  1. *
  2. *
  3. * CAS TEST DU 07/09/17 PROVENANCE : TEST
  4. *******************************************************
  5. *******************************************************
  6. * *
  7. * Test Case of Cyclic Loading *
  8. * for liaison_acbe model *
  9. * on an interface element. *
  10. * Pr: AK & CT *
  11. *******************************************************
  12. *******************************************************
  13. OPTI ECHO 0 ;
  14. OPTI DIME 3 ELEM CUB8 ;
  15. *
  16. dia1 = 0.01 ;
  17. sa1 = (pi*dia1*dia1)/4. ;
  18. dx = 0.02;
  19. *
  20. * Malliage
  21. p0 = 0. 0. 0. ;
  22. p1 = dx 0. 0. ;
  23. d1 = droi p0 p1 1 ;
  24. d2 = d1 plus (0. 0. 0.) ;
  25. intf = racc d1 d2 1.e-5 ;
  26. *
  27. *Loi d'Adherence
  28. pul1 = prog 0. 2.e-6 1.E-3 2.E-3 4.E-3 6.E-3 8.E-3
  29. 1.E-2 1.2E-2 1.3E-2 1.4E-2 1.6E-2 0.02 ;
  30. pul2 = prog 0. 2.e6 7.5E6 12.5E6 17.E6 19.5E6 20.E6
  31. 19.5E6 17.5E6 16.E6 15.5E6 15.E6 15.E6 ;
  32. pulop = evol manu (pul1) (pul2) ;
  33. *
  34. * Materiaux et Modeles
  35. mod1 = mode intf mecanique elastique isotrope
  36. plastique liaison_acbe cos2 ;
  37. mat1 = mate mod1 'PULO' pulop 'KN' 1.e15
  38. 'KS' 1.e12 'SECT' sa1 ;
  39. *
  40. * Conditions aux Limites
  41. cl1 = bloq (p0 et p1) 'DEPL' ;
  42. *
  43. * Chargement
  44. cldep = bloq d2 'UX' ;
  45. dep1 = depi cldep (1) ;
  46. *
  47. * Cyclic Loading
  48. prog1 = prog 0. 1. 2. 3. ;
  49. prog2 = prog (0.) (2.E-2) (-5.E-3) (-2.E-2) ;
  50. temps1 = prog 0. pas 0.1 3. ;
  51. *
  52. evo1 = evol manu (prog1) (prog2) ;
  53. cha1 = char dep1 evo1 DIMP ;
  54. cltot = cl1 et cldep ;
  55. *
  56. tab1 = TABLE;
  57. tab1.MODELE = mod1 ;
  58. tab1.CARACTERISTIQUES = mat1 ;
  59. tab1.BLOCAGES_MECANIQUES = cltot;
  60. tab1.CHARGEMENT = cha1;
  61. tab1.TEMPS_CALCULES = temps1 ;
  62. tab1.PLASTIQUE = vrai;
  63. tab1.'PRECISION' = 1.e-5;
  64. *
  65. pasapas tab1 ;
  66. *
  67. tt = (dime(tab1.temps)) - 1 ;
  68. progx = prog ;
  69. i = 1 ;
  70. repeter bou1 tt ;
  71. f1 = reac cldep
  72. TAB1.DEPLACEMENTS.i ;
  73. ff1 = resu f1 ;
  74. FORC = extr ff1 'FX'
  75. (poin (extr f1 mail) initial) ;
  76. progx = progx et (prog forc) ;
  77. i = i + 1 ;
  78. fin bou1 ;
  79.  
  80. fmax = 'MAXI' progx ;
  81. fmin = 'MINI' PROGX ;
  82. FMAX = (FMAX * 1.E-6) /
  83. (PI * DIA1 * DX) ;
  84. FMIN = (FMIN * 1.E-6) /
  85. (PI * DIA1 * DX) ;
  86. errf = abs(FMAX + FMIN) ;
  87. *
  88. si (> errf 1.e-3) ;
  89. erreur (5);
  90. finsi;
  91. *
  92. fin ;
  93.  
  94.  

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