Télécharger la plaquette Cast3M


fr en

Télécharger TirantLAB2.dgibi

Retour à la liste

Numérotation des lignes : 

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

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