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

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