Télécharger la plaquette Cast3M


fr en

Télécharger test_cisailnl.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  test_cisailnl.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ********************************
  5. *
  6. * Cas test modele global cisail_nl
  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. Gbar = younbar/(2.*(1+nubar));
  24. *
  25. Sbar = Larg*Long;
  26. Ibar = (1./12.)*(Long**3)*Larg;
  27. *
  28. * Courbe contrainte de cisaillement - déformation
  29. *
  30. Fmax = Sbar*2.1;
  31. UXX1 = prog 0. (0.5*Fmax/(Gbar*Sbar)) (2.*Fmax/(Gbar*Sbar))
  32.        (10.*Fmax/(Gbar*Sbar));
  33. FXX1 = prog 0. (0.5*Fmax/Sbar) (Fmax/Sbar) (1.2*Fmax/Sbar);
  34. DXFXTh = evol manu UXX1 (Sbar*FXX1);
  35. *
  36. DMAXP = 1. - ((extr FXX1 3 )/(Gbar*(extr UXX1 3 )));
  37. DMAXN = DMAXP;
  38. DELAP = (extr UXX1 2);
  39. DELAN = DELAP;
  40. Efiss = (1. - DMAXP)*Gbar;
  41. *
  42. XMONOP = PROG;
  43. YMONOP = PROG;
  44. j0 = 2;
  45. repeter LAB1 ((DIME FXX1) - 2 );
  46.    j0=j0+1;
  47.    YY = (extr FXX1 j0);
  48.    XX = (extr UXX1 j0) - ((extr FXX1 j0)/Efiss);
  49.    XMONOP = INSE XMONOP (j0 - 2) (maxi (prog XX 0.));
  50.    YMONOP = INSE YMONOP (j0 - 2) YY;
  51. fin lab1;
  52. MONOP = EVOL MANU XMONOP (Sbar*YMONOP);
  53. MONON = EVOL MANU XMONOP (Sbar*YMONOP);
  54. *
  55. modbar1 = mode bar1 mecanique elastique plastique cisail_nl timo;
  56. matbar1 = mate modbar1 YOUN younbar NU nubar SECT Sbar
  57.      INRZ Ibar DELP delap DMAP dmaxp DELN delan DMAN dmaxn
  58.      BETA 0.2 ALFA 0.5 TETA 1.  MONP monop MONN monon 
  59.      SECY Sbar RHO 0.;
  60. *
  61. BL0 = BLOQ DEPL ROTA p0;
  62. BL1 = BLOQ ROTA p1;
  63. BL2 = BLOQ UX p1;
  64. DEPP1 = DEPI BL2 1.;
  65. Protemp = prog 0. 10. 20. 101.;
  66. ProHori = PROG 0. 0.001 -0.001 0.;
  67. Procalc = PROG 0. pas 0.25 10.;
  68. CharHori = CHAR DIMP DEPP1 (evol manu Protemp ProHori);
  69. *
  70. TAB1 = TABLE;
  71. TAB1 . 'MOVA'          = 'RIEN';
  72. TAB1 . 'TEMPS_CALCULES' = ProCalc ;
  73. *TAB1 . 'CONVERGENCE_FORCEE ' = 'FAUX';
  74. TAB1 .'MAXITERATION'  = 50;
  75. *TAB1.'PRECISION' = 1.D-5;
  76. TAB1 .'CHARGEMENT' = Charhori;
  77. TAB1 .'DYNAMIQUE' = faux;
  78. TAB1 .'MODELE' = modbar1;
  79. TAB1 .'CARACTERISTIQUES'= matbar1;
  80. TAB1 .'BLOCAGES_MECANIQUES' = BL0 ET BL1 ET BL2;
  81. *
  82. PASAPAS TAB1;
  83. *
  84. *
  85. DEPT = TAB1.DEPLACEMENTS;
  86. CONT = TAB1.CONTRAINTES;
  87. TIME = INDEX DEPT;
  88. DIMT = DIME TIME;
  89. *
  90. DEPX = PROG;
  91. GXYT = prog;
  92. EFFX = prog;
  93. *
  94. EFFXXX = PROG;
  95. REPETER BOU1 (DIMT - 1);
  96.   DEP1 = DEPT.(TIME.(&BOU1 + 1));
  97.   CON1 = CONT.(TIME.(&BOU1 + 1));
  98.   FOR1 = REAC (BL0 et BL1 et BL2) DEP1;
  99.   DEPX = DEPX ET (PROG  (extr DEP1 UX p1));
  100.   EFFX = EFFX ET (PROG (extr FOR1 FX p0));
  101. FIN BOU1;
  102. *
  103. DXFX = EVOL MANU DEPX ((-1.)*EFFX);
  104. *
  105. dreal = extr (ipol (prog 0.0005) DEPX ((-1.)*EFFX)) 1;
  106. dtheo =  extr (ipol (prog 0.0005) 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.  

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