Télécharger test_infill.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : test_infill.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ********************************
  5. *
  6. * Cas test modele global infill_uni
  7. *
  8. * D. COMBESCURE - EMSI 1999
  9. *
  10. ********************************
  11. opti dime 2 elem seg2;
  12. *
  13. p0 = 0. 0.;
  14. p1 = 0. 1.;
  15. bar1 = p0 d 1 p1;
  16. *
  17. Long = 0.25;
  18. Larg = 0.25;
  19. *
  20. younbar = 30000.;
  21. *
  22. nubar = 0.;
  23. *
  24. Sbar = Larg*Long;
  25. *
  26. * Courbe contrainte axiale - déformation
  27. *
  28. Fmax = Sbar*2.1;
  29. UXX1 = prog 0. (0.5*Fmax/(younbar*Sbar)) (2.*Fmax/(younbar*Sbar))
  30. (10.*Fmax/(younbar*Sbar)) (20.*Fmax/(younbar*Sbar));
  31. FXX1 = prog 0. (0.5*Fmax/Sbar) (Fmax/Sbar) (1.2*Fmax/Sbar)
  32. (0.6*Fmax/Sbar);
  33. DXFXTh = evol manu UXX1 (Sbar*FXX1);
  34. *
  35. DMAXP = 1. - ((extr FXX1 3 )/(younbar*(extr UXX1 3 )));
  36. DELAP = (extr UXX1 2);
  37. Efiss = (1. - DMAXP)*younbar;
  38. *
  39. XMONOP = PROG;
  40. YMONOP = PROG;
  41. j0 = 2;
  42. repeter LAB1 ((DIME FXX1) - 2 );
  43. j0=j0+1;
  44. YY = (extr FXX1 j0);
  45. XX = (extr UXX1 j0) - ((extr FXX1 j0)/Efiss);
  46. XMONOP = INSE XMONOP (j0 - 2) (maxi (prog XX 0.));
  47. YMONOP = INSE YMONOP (j0 - 2) YY;
  48. fin lab1;
  49. MONOP = EVOL MANU XMONOP (SBAR*YMONOP);
  50. *
  51. modbar1 = mode bar1 mecanique elastique plastique infill_uni
  52. barr;
  53. matbar1 = mate modbar1 YOUN younbar NU nubar SECT Sbar
  54. DELA delap DMAX dmaxp GAMM 0. GAMP 0.
  55. BETA 0.2 ALFA 0.5 TETA 1. MONO monop RHO 0.;
  56. *
  57. BL0 = BLOQ DEPL p0;
  58. BL1 = BLOQ UX p1;
  59. BL2 = BLOQ UY p1;
  60. DEPP1 = DEPI BL2 (-1.);
  61. Protemp = prog 0. 10. 20. 101.;
  62. ProHori = PROG 0. 0.002 -0.002 0.;
  63. Procalc = PROG 0. pas 0.25 10.;
  64. CharHori = CHAR DIMP DEPP1 (evol manu Protemp ProHori);
  65. *
  66. TAB1 = TABLE;
  67. TAB1 . 'MOVA' = 'RIEN';
  68. TAB1 . 'TEMPS_CALCULES' = ProCalc ;
  69. *TAB1 . 'CONVERGENCE_FORCEE ' = 'FAUX';
  70. TAB1 .'MAXITERATION' = 50;
  71. *TAB1.'PRECISION' = 1.D-5;
  72. TAB1 .'CHARGEMENT' = Charhori;
  73. TAB1 .'DYNAMIQUE' = faux;
  74. TAB1 .'MODELE' = modbar1;
  75. TAB1 .'CARACTERISTIQUES'= matbar1;
  76. TAB1 .'BLOCAGES_MECANIQUES' = BL0 ET BL1 ET BL2;
  77. TAB1 .'HYPOTHESE_DEFORMATIONS' = 'LINEAIRE' ;
  78. *
  79. PASAPAS TAB1;
  80. *
  81. *
  82. DEPT = TAB1.DEPLACEMENTS;
  83. CONT = TAB1.CONTRAINTES;
  84. TIME = INDEX DEPT;
  85. DIMT = DIME TIME;
  86. *
  87. TIMP = PROG;
  88. DEPX = PROG;
  89. GXYT = prog;
  90. EFFX = prog;
  91. *
  92. EFFXXX = PROG;
  93. REPETER BOU1 (DIMT - 1);
  94. DEP1 = DEPT.(TIME.(&BOU1 + 1));
  95. CON1 = CONT.(TIME.(&BOU1 + 1));
  96. FOR1 = REAC (BL0) DEP1;
  97. DEPX = DEPX ET (PROG (extr DEP1 UY p1));
  98. EFFX = EFFX ET (PROG (extr FOR1 FY p0));
  99. TIMP = TIMP ET (PROG (TIME.(&BOU1 + 1)));
  100. FIN BOU1;
  101. *
  102. DXFX = EVOL MANU ((-1.)*DEPX) (EFFX);
  103. TIFX = EVOL MANU TIMP (EFFX);
  104. *
  105. dreal = extr (ipol (prog -0.001) DEPX EFFX) 1;
  106. dtheo = extr (ipol (prog 0.001) UXX1 (Sbar*FXX1)) 1;
  107. prec = abs (dreal - dtheo);
  108. *
  109. xlim0 = 1.D-7;
  110. SI (prec > xlim0);
  111. 'ERRE' 5;
  112. SINON;
  113. 'ERRE' 0;
  114. FINSI;
  115. *
  116. FIN;
  117. *
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  

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