Télécharger mvm_bcn.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : mvm_bcn.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ********************************************************************************
  5. * TEST : Modified Von Mises nonlocal damage model *
  6. * UNITS: *
  7. * Dimensions in mm *
  8. * Loads in N *
  9. * Modulus in N/mm2 (Mpa) *
  10. * *
  11. ********************************************************************************
  12. OPTI DIME 2 ELEM QUA4 MODE PLAN DEFO ;
  13. *-------------------------------------------------------------------------------
  14. * Generation of mesh
  15. *-------------------------------------------------------------------------------
  16. p1 = 0. 0.;
  17. p2 = 1. 0.;
  18. l12 = D 1 p1 p2;
  19. ma1 = l12 TRAN 1 (0. 1);
  20. *
  21. p1 = POIN ma1 PROC (0. 0.) ;
  22. p2 = POIN ma1 PROC (1. 0.) ;
  23. p3 = POIN ma1 PROC (0. 1) ;
  24. p4 = POIN ma1 PROC (1. 1) ;
  25. *
  26. l12 = (CONT ma1) ELEM COMPRIS p1 p2 ;
  27. l13 = (CONT ma1) ELEM COMPRIS p3 p1 ;
  28. l24 = (CONT ma1) ELEM COMPRIS p2 p4 ;
  29. l34 = (CONT ma1) ELEM COMPRIS p4 p3 ;
  30. *-------------------------------------------------------------------------------
  31. * Definition of the model and material parameters
  32. *-------------------------------------------------------------------------------
  33. mobet = MODE ma1 MECANIQUE ELASTIQUE ISOTROPE ENDOMMAGEMENT MVM ;
  34. *
  35. b1 = 15000. ; b2 = 0.5 ;
  36. y0 = 1.e-4 ; ratio_ct = 10. ;
  37. ll = 2. ; nuc = 0. ;
  38. younc = 1.e4 ; loi = 1 ;
  39. *
  40. mabet = MATE mobet 'YOUN' younc 'NU' nuc RHO 1.
  41. 'Y0' y0 'RATI' ratio_ct
  42. 'B1' b1 'B2' b2 'LOI' loi ;
  43. *-------------------------------------------------------------------------------
  44. * Boundary conditions
  45. *-------------------------------------------------------------------------------
  46. * 1) Fixed boundaries
  47. bcfix = (BLOQUE L13 DEPL) ;
  48. * 2) Prescrite Displacements
  49. bcdisp = BLOQUE (L24) UX ;
  50. *-------------------------------------------------------------------------------
  51. * Adjustment of the reference load
  52. *-------------------------------------------------------------------------------
  53. chdisp1 = MANU CHPO L24 1 'UX' 10*y0 ;
  54. ffextp = DEPI bcdisp chdisp1 ;
  55. evol1 = EVOL MANU Temps (PROG 0. 1.) Depl (PROG 0. 1.) ;
  56. cha1 = CHAR DIMP evol1 ffextp ;
  57. *-------------------------------------------------------------------------------
  58. * Connectivities (for nonlocal averaging)
  59. *-------------------------------------------------------------------------------
  60. conn1 = CONN mobet ll 'NORMAL' 'INTERIEUR' ;
  61. *-------------------------------------------------------------------------------
  62. * Input for pasapas procedure
  63. *-------------------------------------------------------------------------------
  64. t0 = TABLE ;
  65. t0.'MODELE' = mobet ;
  66. t0.'CARACTERISTIQUES' = mabet ;
  67. t0.'CHARGEMENT' = cha1 ;
  68. * il y a un probleme pour le temps 0.1 PV
  69. t0.'TEMPS_CALCULES' = PROG 0. 0.05 0.099 PAS 0.05 1.;
  70. t0.'BLOCAGES_MECANIQUES' = bcfix et bcdisp ;
  71. t0.'MAXITERATION' = 100 ;
  72. t0.'NACCE' = 100 ;
  73. t0.'MOVA' = 'D ' ;
  74. t0.'NON_LOCAL' = 'MOT' 'MOYE' ;
  75. t0.'CONN' = conn1 ;
  76. t0.'LVIAM' = (MOTS EPTI) ;
  77. t0.'FACTMIN' = .001 ;
  78. t0.'MAXSTEP' = 10000 ;
  79. t0.'PRECISION' = 1.E-10 ;
  80. t0.'PRECFLEX' = 1.E-10 ;
  81. t0.'CONVERGENCE_FORCEE' = FAUX ;
  82. t0.'ACCELERATION' = 50 ;
  83. *
  84. PASAPAS t0 ;
  85. *
  86. *------------------------------------------------------------------------------
  87. * Evolution of stress.vs.strain and force.vs.displacement
  88. *------------------------------------------------------------------------------
  89. sig0 = t0.CONTRAINTES ;
  90. dep0 = t0.DEPLACEMENTS ;
  91. var0 = t0.VARIABLES_INTERNES ;
  92. for0 = t0.REACTIONS ;
  93. ind0 = INDEX sig0 ;
  94. nind = DIME for0 ;
  95. *
  96. for = PROG ;
  97. dep = PROG ;
  98. sigma = PROG ;
  99. epsilon = PROG ;
  100. *
  101. REPETER indj nind ;
  102. i = &indj - 1 ;
  103. *
  104. sii = EXCO 'SMXX' (sig0.i) 'SCAL' ;
  105. epi = EXCO 'EPXX' (EPSI mobet (dep0.i)) 'SCAL' ;
  106. sigp2 = EXTR sii 'SCAL' 1 1 1 ;
  107. sigma = sigma INSE &indj (sigp2) ;
  108. epsp2 = EXTR epi 'SCAL' 1 1 1 ;
  109. epsilon = epsilon INSE &indj (epsp2) ;
  110. *
  111. fi = EXTR (for0.i) FX p2 ;
  112. di = EXTR (dep0.i) UX p2 ;
  113. for = for INSE &indj (fi) ;
  114. dep = dep INSE &indj (di) ;
  115. FIN indj ;
  116. *
  117. SI ((ABS (MAXI ((0.5*sigma) - for))) < 0.5e-5) ;
  118. ERRE 0 ;
  119. SINON ;
  120. mess 'test :' (ABS (MAXI ((0.5*sigma) - for)));
  121. ERRE 5 ;
  122. FINSI ;
  123.  
  124. FIN;
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  
  133.  
  134.  
  135.  
  136.  

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