Télécharger la plaquette Cast3M


fr en

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.       'NON_LOCAL' 'MOYE' 'V_MOYENNE' (MOTS 'EPTI') ;
  35. *
  36.  b1    = 15000.     ; b2       = 0.5     ;
  37.  y0    = 1.e-4      ; ratio_ct = 10.     ;
  38.  ll    = 2.         ; nuc   = 0.         ;
  39.  younc = 1.e4       ; loi   = 1          ; 
  40. *   
  41.  mabet = MATE mobet 'YOUN' younc 'NU' nuc RHO 1. 
  42.                     'Y0' y0 'RATI' ratio_ct 
  43.                     'B1' b1 'B2' b2 'LOI' loi 'LCAR' LL ;
  44. *-------------------------------------------------------------------------------
  45. * Boundary conditions
  46. *-------------------------------------------------------------------------------
  47. * 1) Fixed boundaries
  48.  bcfix = (BLOQUE L13 DEPL) ;
  49. * 2) Prescrite Displacements
  50.  bcdisp = BLOQUE  (L24) UX ;
  51. *-------------------------------------------------------------------------------
  52. * Adjustment of the reference load
  53. *-------------------------------------------------------------------------------
  54.  chdisp1 = MANU CHPO L24 1 'UX' 10*y0 ;
  55.  ffextp  = DEPI bcdisp chdisp1 ;
  56.  evol1  = EVOL MANU Temps (PROG 0. 1.) Depl (PROG 0. 1.) ;
  57.  cha1   = CHAR DIMP evol1 ffextp ;
  58. *-------------------------------------------------------------------------------
  59. * Connectivities (for nonlocal averaging)
  60. *-------------------------------------------------------------------------------
  61.  conn1 = CONN mobet MABET 'NORMAL' 'INTERIEUR' ;
  62. *-------------------------------------------------------------------------------
  63. * Input for pasapas procedure
  64. *-------------------------------------------------------------------------------
  65.  t0 = TABLE                                   ;
  66.  t0.'MODELE'             = mobet              ;
  67.  t0.'CARACTERISTIQUES'   = mabet              ;
  68.  t0.'CHARGEMENT'         = cha1               ;
  69. * il y a un probleme pour le temps 0.1  PV
  70.  t0.'TEMPS_CALCULES'     = PROG 0. 0.05  0.099 PAS 0.05 1.;
  71.  t0.'BLOCAGES_MECANIQUES' =  bcfix  et bcdisp ;
  72.  t0.'MAXITERATION'       = 100                ;
  73.  t0.'NACCE'              = 100                ;
  74.  t0.'MOVA'               = 'D   '             ;
  75.  t0.'CONN'               = conn1              ;
  76.  t0.'FACTMIN'            = .001               ;
  77.  t0.'MAXSTEP'            = 10000              ;
  78.  t0.'PRECISION'          = 1.E-10             ;
  79.  t0.'PRECFLEX'           = 1.E-10             ;  
  80.  t0.'CONVERGENCE_FORCEE' = FAUX               ;
  81.  t0.'ACCELERATION'       = 50                 ;
  82. *
  83.  PASAPAS t0 ;
  84. *
  85. *------------------------------------------------------------------------------
  86. * Evolution of stress.vs.strain and force.vs.displacement
  87. *------------------------------------------------------------------------------
  88.  sig0 = t0.CONTRAINTES        ;
  89.  dep0 = t0.DEPLACEMENTS       ;
  90.  var0 = t0.VARIABLES_INTERNES ;
  91.  for0 = t0.REACTIONS          ; 
  92.  ind0 = INDEX sig0            ;
  93.  nind  = DIME for0            ;
  94. *
  95.  for     = PROG ;
  96.  dep     = PROG ;
  97.  sigma   = PROG ;
  98.  epsilon = PROG ;
  99. *
  100. REPETER indj  nind ;
  101.    i = &indj - 1 ;
  102. *
  103.    sii = EXCO 'SMXX' (sig0.i) 'SCAL' ;
  104.    epi = EXCO 'EPXX' (EPSI mobet (dep0.i)) 'SCAL' ;
  105.    sigp2 = EXTR sii 'SCAL' 1 1 1 ;
  106.    sigma = sigma INSE &indj (sigp2) ;
  107.    epsp2 = EXTR epi 'SCAL' 1 1 1 ;
  108.    epsilon = epsilon INSE &indj (epsp2) ;
  109. *
  110.    fi = EXTR (for0.i) FX p2 ;
  111.    di = EXTR (dep0.i) UX p2 ;
  112.    for      = for    INSE &indj (fi) ;
  113.    dep      = dep   INSE &indj (di) ;
  114. FIN indj ;
  115. *
  116.  SI ((ABS (MAXI ((0.5*sigma) - for))) < 0.5e-5) ;
  117.    ERRE 0 ;
  118.  SINON ;
  119.  mess 'test :' (ABS (MAXI ((0.5*sigma) - for)));
  120.   ERRE 5 ;
  121.  FINSI ;
  122.  
  123. FIN;
  124.  
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  
  133.  
  134.  
  135.  
  136.  
  137.  

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