Télécharger la plaquette Cast3M


fr en

Télécharger mistral_axi2.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  mistral_axi2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. *     Test mistral_axi2.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 mode axis elem qua8 ;
  13. *
  14. *******************************************************
  15. *                TEST DE VALIDATION                   *
  16. *                ------------------                   *
  17. *                                                     *
  18. *  MODELE MISTRAL                                     *
  19. *    ELASTICITE ET VISCOPLASTICITE                    *
  20. *                                                     *
  21. *  MAILLAGE:                                          *
  22. *    TUBE CYLINDRIQUE (AXISYMETRIQUE)                 *
  23. *                                                     *
  24. *  CHARGEMENT:                                        *
  25. *    PRESSION IMPOSEE : RAMPE RAPIDE ET PLATEAU       *
  26. *******************************************************
  27. *
  28. * Geometrie
  29. *
  30. * ------------------------------------------
  31.   RI = 49.5E-3 ; RE = 50.5E-3 ; H = 0.010 ;
  32. * ------------------------------------------
  33.    RM = (RI+RE)/2. ;
  34.    RE2MRI2 = (RE**2)-(RI**2) ;
  35.    COEFI = (RI**2)/RE2MRI2 ;
  36.    E = RE-RI ;
  37.    VECR = 1. 0. ;
  38. *
  39.    PIB = RI 0. ;  PEB = RE 0. ;
  40.    PIH = RI H  ;  PEH = RE H  ;
  41. *
  42.    LB = droi 2 PIB PEB ;
  43.    LE = droi 1 PEB PEH ;
  44.    LH = droi 2 PEH PIH ;
  45.    LI = droi 1 PIH PIB ;
  46. *
  47. * Maillage
  48. *
  49.    EPROU = 'DALLER' LB LE LH LI ;
  50. *  'TRAC' EPROU ;
  51. *
  52.    POMB = RM 0. ; POMH = RM H ;
  53.    PMB = EPROU poin proc POMB ;
  54.    PMH = EPROU poin proc POMH ;
  55. *
  56. * Modele et materiau
  57. *
  58.    MODG = mode EPROU mecanique elastique orthotrope
  59.                                viscoplastique mistral ;
  60. *
  61. * ------------------------------------------------------------
  62.   SENSIP1 = 1 ; SENSIP2 = -3 ; ICBASE = 0 ;
  63.   fichier = 'CHAINE' DIVERS '/mimataxi2_par' ;
  64.   PDILT    E1 E2 E3 NU12 NU23 NU13 MU12 MU23 MU13
  65.   PNBRE PCOHI PECOU PEDIR PRVCE PECRX PDVDI PCROI PINCR
  66.     = @mistpar fichier SENSIP1 SENSIP2 ;
  67. * ------------------------------------------------------------
  68.    MAZY = mate MODG 'YG1 ' E1   'YG2 ' E2   'YG3 ' E3
  69.                     'NU12' NU12 'NU23' NU23 'NU13' NU13
  70.                     'G12 ' MU12
  71.                     'ALP1' 0.     'ALP2' 0.      'ALP3' 0. 'TALP' 0. 'TREF' 623.15
  72.  'DILT' PDILT 'NBRE' PNBRE 'COHI' PCOHI 'ACOU' PECOU 'EDIR' PEDIR
  73.  'RVCE' PRVCE 'ECRX' PECRX 'DVDI' PDVDI 'CROI' PCROI 'INCR' PINCR
  74.                     'SIP1' SENSIP1  'SIP2' SENSIP2  'IBAS' ICBASE
  75.                     'DIRECTION' VECR ;
  76. *
  77. * Conditions aux limites
  78. *
  79.    CLB = bloq uz LB ;
  80.    CLH = rela ense uz LH ;
  81.    CL = CLB et CLH ;
  82. *
  83. * Chargement
  84. *
  85.    T1 = 0.001 ; T2 = 2000.001 ;
  86.    DT0 = 0.0001 ;
  87.    DT1 = 5. ;
  88.    TEMPSCAL = prog 0. pas DT0 T1 pas DT1 100.001
  89.                    150.001 200.001 300.001
  90.                    500.001 700.001 1000.001 1400.001 T2 ;
  91.    SIGTT1 = 400.E6 ;
  92.    TT0 = 623.15 ; TT1 = 623.15 ; TT2 = 623.15 ;
  93.    PHIT0 = 2.E25 ;
  94.    PPSI0 = 112.23 ;
  95. *
  96.    TEMPS = prog 0. T1 T2 ;
  97. *
  98.    PRI = E/RI ;
  99.    PRESI = pres mass MODG (PRI) LI ;
  100.    SZPRI = PRI*COEFI ;
  101.    PRHPI = (-1.)*SZPRI ;
  102.    PRESHPI = pres mass MODG (PRHPI) LH ;
  103.    EVPI = evol manu TEMPS (prog 0. SIGTT1 SIGTT1) ;
  104.    CHAPI = char meca (PRESI et PRESHPI) EVPI ;
  105. *
  106.    TT = manu chpo EPROU 1 'T' 1. ;
  107.    EVTT = evol manu TEMPS (prog TT0 TT1 TT2) ;
  108.    CHTT = char 'T' TT EVTT ;
  109. *
  110.    PHI = manu chpo EPROU 1 'FI' 1. ;
  111.    EVFI = evol manu TEMPS (prog 0. 0. 0.) ;
  112.    CHFI = char 'FI' PHI EVFI ;
  113. *
  114.    CHA = CHAPI et CHTT et CHFI ;
  115. *
  116. * Valeurs initiales
  117. *
  118.    VINT0 = zero MODG 'VARINTER' ;
  119.    FIT0 = manu chml MODG 'FIT ' PHIT0 type 'SCALAIRE' 'STRESSES' ;
  120.    VINT0 = VINT0 + FIT0 ;
  121.    PSI0 = manu chml MODG 'PSI ' PPSI0 type 'SCALAIRE' 'STRESSES' ;
  122.    VINT0 = VINT0 + PSI0 ;
  123. *
  124. * Calcul
  125. *
  126.    TAB = TABLE ;
  127.    TAB.'VARIABLES_INTERNES' = TABLE ;
  128.    TAB.'BLOCAGES_MECANIQUES' = CL ;
  129.    TAB.'CARACTERISTIQUES' = MAZY ;
  130.    TAB.'MODELE' = MODG ;
  131.    TAB.'CHARGEMENT' = CHA ;
  132.    TAB.'VARIABLES_INTERNES' . 0 = VINT0 ;
  133.    TAB.'TEMPS_CALCULES' = TEMPSCAL ;
  134.    PASAPAS TAB ;
  135. *
  136. * Traitement des resultats
  137. *
  138.    TEM = TAB.'TEMPS' ;
  139.    SIG = TAB.'CONTRAINTES' ;
  140.    DEP = tab.'DEPLACEMENTS' ;
  141.    VI  = tab.'VARIABLES_INTERNES' ;
  142.    NT = dime TEMPSCAL-1 ;
  143. * ---------------
  144.    ERMAX = 0.01 ;
  145. * ---------------
  146.    T = TEM.NT ;
  147. *
  148.    SIGT = SIG.NT ;
  149.    SIGT_PO = chang chpo MODG SIGT ;
  150.    SIGRR = extr SIGT_PO SMRR PMB ;
  151.    SIGTT = extr SIGT_PO SMTT PMB ;
  152.    SIGZZ = extr SIGT_PO SMZZ PMB ;
  153.    SIGRZ = extr SIGT_PO SMRZ PMB ;
  154. *
  155.    DEPT = DEP.NT ;
  156.    EPSIT = epsi DEPT MODG ;
  157.    EPSIT_PO = chang chpo MODG EPSIT ;
  158.    EPSRR = extr EPSIT_PO EPRR PMB ;
  159.    EPSTT = extr EPSIT_PO EPTT PMB ;
  160.    EPSZZ = extr EPSIT_PO EPZZ PMB ;
  161.    EPSRZ = (extr EPSIT_PO GARZ PMB)/2. ;
  162. *
  163.    VIT = VI.NT ;
  164.    VIT_PO = chang chpo MODG VIT ;
  165.    EPSFRR = extr VIT_PO EP11 PMB ;
  166.    EPSFTT = extr VIT_PO EP12 PMB ;
  167.    EPSFZZ = extr VIT_PO EP13 PMB ;
  168.    EPSFRZ = extr VIT_PO EP15 PMB ;
  169. *
  170.    mess ;
  171.    mess 'T :                          ' T ;
  172.    mess ;
  173.    mess 'SIGRR SIGTT SIGZZ :    ' SIGRR SIGTT SIGZZ ;
  174.    mess 'SIGRZ :                ' SIGRZ ;
  175.    mess 'EPSRR EPSTT EPSZZ :    ' EPSRR EPSTT EPSZZ ;
  176.    mess 'EPSRZ :                ' EPSRZ ;
  177.    mess 'EPSFRR EPSFTT EPSFZZ : ' EPSFRR EPSFTT EPSFZZ ;
  178.    mess 'EPSFRZ               : ' EPSFRZ ;
  179.    mess ;
  180. *
  181. *   SIGT_RR = exco SIGT 'SMRR' ;
  182. *   trac SIGT_RR MODG ;
  183. *   SIGT_TT = exco SIGT 'SMTT' ;
  184. *   trac SIGT_TT MODG ;
  185. *   SIGT_ZZ = exco SIGT 'SMZZ' ;
  186. *   trac SIGT_ZZ MODG ;
  187. *   SIGT_RZ = exco SIGT 'SMRZ' ;
  188. *   trac SIGT_RZ MODG ;
  189. *
  190. * ------------------------------------------------------------
  191.    YOUNGR = 1.1E11 ; YOUNGZ = 0.9E11 ; YOUNGT = 1.E11 ;
  192.    POISRZ = 0.3 ; POISZT = 0.3 ; POISRT = 0.3 ;
  193.    MURZ = 4.E10 ;
  194. * ------------------------------------------------------------
  195.    SIGRR1 = (-1.)*SIGTT1*E/(2.*RM) ;
  196.    SIGZZ1 = SIGTT1*SZPRI ;
  197.    EPSERR1 = (SIGRR1/YOUNGR)-((POISRT/YOUNGR)*SIGTT1)
  198.                             -((POISRZ/YOUNGR)*SIGZZ1) ;
  199.    EPSETT1 = (SIGTT1/YOUNGT)-((POISZT/YOUNGZ)*SIGZZ1)
  200.                             -((POISRT/YOUNGR)*SIGRR1) ;
  201.    EPSEZZ1 = (SIGZZ1/YOUNGZ)-((POISRZ/YOUNGR)*SIGRR1)
  202.                             -((POISZT/YOUNGZ)*SIGTT1) ;
  203. * --------------------------
  204.    EPSFRR1 = -0.00160863 ;
  205.    EPSFTT1 = 0.001695925 ;
  206.    EPSFZZ1 = -0.00008729 ;
  207. * --------------------------
  208.    EPSRR1 = EPSERR1 + EPSFRR1 ;
  209.    EPSTT1 = EPSETT1 + EPSFTT1 ;
  210.    EPSZZ1 = EPSEZZ1 + EPSFZZ1 ;
  211. *
  212.    ERSIGRR = (SIGRR/SIGRR1) - 1. ;
  213.    ERSIGTT = (SIGTT/SIGTT1) - 1. ;
  214.    ERSIGZZ = (SIGZZ/SIGZZ1) - 1. ;
  215.    P_ER = prog ERSIGRR ERSIGTT ERSIGZZ ;
  216.    ERSIGMA = maxi P_ER 'ABS' ;
  217. *
  218.    EREPSRR = (EPSRR/EPSRR1) - 1. ;
  219.    EREPSTT = (EPSTT/EPSTT1) - 1. ;
  220.    EREPSZZ = (EPSZZ/EPSZZ1) - 1. ;
  221.    P_ER = prog EREPSRR EREPSTT EREPSZZ ;
  222.    EREPSMA = maxi P_ER 'ABS' ;
  223. *
  224.    EREPSFRR = (EPSFRR/EPSFRR1) - 1. ;
  225.    EREPSFTT = (EPSFTT/EPSFTT1) - 1. ;
  226.    EREPSFZZ = (EPSFZZ/EPSFZZ1) - 1. ;
  227.    P_ER = prog EREPSFRR EREPSFTT EREPSFZZ ;
  228.    EREPSFMA = maxi P_ER 'ABS' ;
  229. *
  230.    si (ERSIGMA > ERMAX) ;
  231.     mess ;
  232.     mess 'ABS(erreur relative) sur SIGRR, SIGTT ou SIGZZ > ' ERMAX ;
  233.     mess ERSIGRR ERSIGTT ERSIGZZ ;
  234.    finsi ;
  235.    si (EREPSMA > ERMAX) ;
  236.     mess ;
  237.     mess 'ABS(erreur relative) sur EPSRR, EPSTT ou EPSZZ > ' ERMAX ;
  238.     mess EREPSRR EREPSTT EREPSZZ ;
  239.    finsi ;
  240.    si (EREPSFMA > ERMAX) ;
  241.     mess ;
  242.     mess 'ABS(erreur relative) sur EPSFRR, EPSFTT ou EPSFZZ > ' ERMAX ;
  243.     mess EREPSFRR EREPSFTT EREPSFZZ ;
  244.     mess ;
  245.    finsi ;
  246.    P_ER = prog ERSIGMA EREPSMA EREPSFMA ;
  247.    ERMA = maxi P_ER ;
  248.    si (ERMA > ERMAX) ;
  249.     ERRE 5 ;
  250.    finsi ;
  251.  
  252.    fin ;
  253.  
  254.  
  255.  
  256.  
  257.  
  258.  
  259.  
  260.  
  261.  
  262.  
  263.  
  264.  

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