Télécharger la plaquette Cast3M


fr en

Télécharger mistral_cp.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  mistral_cp.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. *     Test mistral_cp.dgibi: Jeux de donnees          *
  6. *     ---------------------------------------         *
  7. *                                                     *
  8. * repertoire des fichiers "divers"
  9. DIVERS = VENV 'CASTEM_DIVERS';
  10. *
  11.   opti echo 0 ;
  12.   opti dime 2 elem qua8 ;
  13.   opti mode plan cont ;
  14. *
  15. *******************************************************
  16. *                TEST DE VALIDATION                   *
  17. *                ------------------                   *
  18. *  MODELE MISTRAL                                     *
  19. *    ELASTICITE ET PLASTICITE INSTANTANE              *
  20. *                                                     *
  21. *  MAILLAGE:                                          *
  22. *     EPROUVETTE RECTANGULAIRE                        *
  23. *                                                     *
  24. *  CHARGEMENT:                                        *
  25. *     DEPLACEMENT LATERAL IMPOSE MONOTONE CROISSANT   *
  26. *******************************************************
  27. *
  28. * Geometrie
  29. *
  30. * ------------------------------------------
  31.   XX = 1. ; YY = 0.5 ;
  32. * ------------------------------------------
  33. *
  34.    P00 = 0. 0. ; P10 = XX 0.  ;
  35.    P11 = XX YY   ;  P01 = 0. YY  ;
  36. *
  37.    LX0 = droi 1 P00 P10 ;
  38.    LY1 = droi 1 P10 P11 ;
  39.    LX1 = droi 1 P11 P01 ;
  40.    LY0 = droi 1 P01 P00 ;
  41. *
  42. * Maillage
  43. *
  44.    EPROU = dall LX0 LY1 LX1 LY0 ;
  45. *   trac EPROU ;
  46. *
  47. * Modele et materiau
  48. *
  49.    MODL = mode EPROU mecanique elastique orthotrope
  50.                                  viscoplastique mistral ;
  51. * ------------------------------------------------------------
  52.   VEC1 = 1. 1. ;  TETA = 45. ;
  53.   SENSIP1 = -2 ; SENSIP2 = 1 ; ICBASE = 0 ;
  54.   fichier = 'CHAINE' DIVERS '/mimatcp_par' ;
  55.   PDILT    E1 E2 E3 NU12 NU23 NU13 MU12 MU23 MU13
  56.   PNBRE PCOHI PECOU PEDIR PRVCE PECRX PDVDI PCROI PINCR
  57.     = @mistpar fichier SENSIP1 SENSIP2 ;
  58. * ------------------------------------------------------------
  59.    MATER = mate MODL
  60.              'YG1 ' E1 'YG2 ' E2 'YG3 ' E3
  61.              'NU12' NU12 'NU23' NU23 'NU13' NU13
  62.              'G12 ' MU12
  63.              'ALP1' 0.     'ALP2' 0. 'TALP' 0. 'TREF' 300.
  64.  'DILT' PDILT 'NBRE' PNBRE 'COHI' PCOHI 'ACOU' PECOU 'EDIR' PEDIR
  65.  'RVCE' PRVCE 'ECRX' PECRX 'DVDI' PDVDI 'CROI' PCROI 'INCR' PINCR
  66.                     'SIP1' SENSIP1  'SIP2' SENSIP2  'IBAS' ICBASE
  67.              'DIRECTION' VEC1 'INCLINE' TETA ;
  68. *
  69. * Conditions aux limites
  70. *
  71.    CLX0 = bloq ux LX0 ;
  72.    CLX1 = bloq ux LX1 ;
  73.    CLY0 = bloq uy LY0 ;
  74.    CLY1 = bloq uy LY1 ;
  75.    CL = CLX0 et CLX1 et CLY0 et CLY1 ;
  76. *
  77. * Chargement
  78. *
  79. * -------------------------------------------------------
  80.    T1 = 100. ; DT1 = 10. ;
  81.    EPSXY1 = 0.01 ;
  82.    TT0 = 300. ; TT1 = TT0 ;
  83.    PHIT0 = 0. ;
  84. * -------------------------------------------------------
  85.    TEMPS = prog 0. T1 ;
  86.    TEMPSCAL = prog 0. pas DT1 T1 ;
  87. *
  88.    DEPX = depi CLX1 YY ;
  89.    DEPY = depi CLY1 XX ;
  90.    EVEPS = evol manu TEMPS (prog 0. EPSXY1) ;
  91.    CHADE = char dimp (DEPX et DEPY) EVEPS ;
  92. *
  93.    TT = manu chpo EPROU 1 'T' 1. ;
  94.    EVTT = evol manu TEMPS (prog TT0 TT1) ;
  95.    CHTT = char 'T' TT EVTT ;
  96. *
  97.    PHI = manu chpo EPROU 1 'FI' 1. ;
  98.    EVFI = evol manu TEMPS (prog 0. 0.) ;
  99.    CHFI = char 'FI' PHI EVFI ;
  100. *
  101.    CHA = CHADE et CHTT et CHFI ;
  102. *
  103. * Valeurs initiales
  104. *
  105.    VINT0 = zero MODL 'VARINTER' ;
  106.    FIT0 = manu chml MODL 'FIT ' PHIT0 type 'SCALAIRE' 'STRESSES' ;
  107.    VINT0 = VINT0 + FIT0 ;
  108. *
  109. * Calcul
  110. *
  111.    TAB = TABLE ;
  112.    TAB.'VARIABLES_INTERNES' = TABLE ;
  113.    TAB.'BLOCAGES_MECANIQUES' = CL ;
  114.    TAB.'CARACTERISTIQUES' = MATER ;
  115.    TAB.'MODELE' = MODL ;
  116.    TAB.'CHARGEMENT' = CHA ;
  117.    TAB.'VARIABLES_INTERNES' . 0 = VINT0 ;
  118.    TAB.'TEMPS_CALCULES' = TEMPSCAL ;
  119.    TAB.'HYPOTHESE_DEFORMATIONS' = 'LINEAIRE' ;
  120.    PASAPAS TAB ;
  121. *
  122. * Traitement des resultats
  123. *
  124.    SIG = TAB.'CONTRAINTES' ;
  125.    DEP = tab.'DEPLACEMENTS' ;
  126.    VI  = tab.'VARIABLES_INTERNES' ;
  127.    NT = dime TEMPSCAL-1 ;
  128. * ----------------
  129.    ERMAX = 0.001 ;
  130. * ----------------
  131. *
  132.    SIGT = SIG.NT ;
  133.    SIGT_PO = chang chpo MODL SIGT ;
  134.    SIGXX = extr SIGT_PO SMXX P00 ;
  135.    SIGYY = extr SIGT_PO SMYY P00 ;
  136.    SIGZZ = extr SIGT_PO SMZZ P00 ;
  137.    SIGXY = extr SIGT_PO SMXY P00 ;
  138. *
  139.    DEPT = DEP.NT ;
  140.    EPSIT = epsi DEPT MODL 'LINE' ;
  141.    EPSIT_PO = chang chpo MODL EPSIT ;
  142.    EPSXX = extr EPSIT_PO EPXX P00 ;
  143.    EPSYY = extr EPSIT_PO EPYY P00 ;
  144.    EPSZZ = extr EPSIT_PO EPZZ P00 ;
  145.    EPSXY = (extr EPSIT_PO GAXY P00)/2. ;
  146. *
  147.    VIT = VI.NT ;
  148.    VIT_PO = chang chpo MODL VIT ;
  149.    EPSP11 = extr VIT_PO EP01 P00 ;
  150.    EPSP22 = extr VIT_PO EP02 P00 ;
  151.    EPSP33 = extr VIT_PO EP03 P00 ;
  152.    EPSP12 = extr VIT_PO EP04 P00 ;
  153.    EPSP13 = extr VIT_PO EP05 P00 ;
  154.    EPSP23 = extr VIT_PO EP06 P00 ;
  155. *
  156.    mess ;
  157.    mess 'SIGXX SIGYY SIGZZ :    ' SIGXX SIGYY SIGZZ ;
  158.    mess 'SIGXY :                ' SIGXY ;
  159.    mess ;
  160.    mess 'EPSXX EPSYY EPSZZ :    ' EPSXX EPSYY EPSZZ ;
  161.    mess 'EPSXY :                ' EPSXY ;
  162.    mess ;
  163.    mess 'EPSP11 EPSP22 EPSP33 : ' EPSP11 EPSP22 EPSP33 ;
  164.    mess 'EPSP12 EPSP13 EPSP23 : ' EPSP12 EPSP13 EPSP23 ;
  165.    mess ;
  166.  
  167. *   SIGT_XY = exco SIGT 'SMXY' ;
  168. *   trac SIGT_XY MODL ;
  169. *   GAMA_XY = exco EPSIT 'GAXY' ;
  170. *   trac GAMA_XY MODL ;
  171. * ------------------------------------------------------------
  172.    MUXY = 4.E10 ;
  173. * ------------------------------------------------------------
  174.    AA = 400.E6 ; BB = 1.E10 ;
  175.    A0 = AA/(2.*MUXY) ;
  176.    A1 = 1. + (BB/MUXY) ;
  177.    EPSPXY1 = (EPSXY1-A0)/A1 ;
  178.    SIGXY1 = AA + (2.*BB*EPSPXY1) ;
  179.    EPSEXY1 = SIGXY1/(2.*MUXY) ;
  180. * -----------------------------------------------------
  181. *
  182.    ERSIGXY = ABS( (SIGXY/SIGXY1) - 1. ) ;
  183.    EREPSXY = ABS( (EPSXY/EPSXY1) - 1. ) ;
  184.    EREPSPXY = ABS( (EPSP12/EPSPXY1) - 1. ) ;
  185. *
  186.    si (ERSIGXY > ERMAX) ;
  187.     mess ;
  188.     mess 'ABS(erreur relative) sur la contrainte > ' ERMAX ;
  189.     mess  ERSIGXY ;
  190.    finsi ;
  191.    si (EREPSXY > ERMAX) ;
  192.     mess ;
  193.     mess 'ABS(erreur relative) sur la deformation > ' ERMAX ;
  194.     mess EREPSXY ;
  195.    finsi ;
  196.    si (EREPSPXY > ERMAX) ;
  197.     mess ;
  198.     mess 'ABS(erreur relative) sur la deformation plastique > ' ERMAX ;
  199.     mess EREPSPXY ;
  200.     mess ;
  201.    finsi ;
  202.    P_ER = prog ERSIGXY EREPSXY EREPSPXY ;
  203.    ERMA = maxi P_ER ;
  204.    si (ERMA > ERMAX) ;
  205.     ERRE 5 ;
  206.    finsi ;
  207.  
  208.    fin ;
  209.  
  210.  
  211.  
  212.  
  213.  
  214.  
  215.  
  216.  
  217.  
  218.  
  219.  
  220.  
  221.  
  222.  
  223.  
  224.  

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