Télécharger nouailhas_b1.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : nouailhas_b1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. **********************************************************************
  5. * *
  6. * fichier : nouailhas_b1.dgibi *
  7. * *
  8. **********************************************************************
  9. * *
  10. * fait par PLG le 03-09-97 *
  11. * *
  12. **********************************************************************
  13. * *
  14. * But : ktest pour la loi viscoplastique NOUAILHAS_B *
  15. * *
  16. **********************************************************************
  17. i_trace = faux ;
  18. *
  19. debproc perso1 plg*TABLE ;
  20. *
  21. li_t = plg.depo.li_tp ;
  22. li_d = plg.depo.li_def ;
  23. *
  24. ind_1 = index (plg.temps) ;
  25. nb_1 = dime ind_1 ;
  26. *
  27. v_t = plg.temps.(ind_1.nb_1) ;
  28. v_d = extr (plg.deformations_inelastiques.(ind_1.nb_1)) EIZZ 1 1 1 ;
  29. *
  30. li_t = li_t et (prog v_t) ;
  31. li_d = li_d et (prog (0.0 - v_d)) ;
  32. ev_d = evol manu ' ' li_t ' ' li_d ;
  33. *
  34. plg.depo.li_tp = li_t ;
  35. plg.depo.li_def = li_d ;
  36. *
  37. dess ((plg.depo.ev_ref) et ev_d) NCLK ;
  38. finproc ;
  39. **********************************************************************
  40. option dime 2 ;
  41. option elem qua8 ;
  42. option mode axis ;
  43. *
  44. a00 = 0.0 0.0 ;
  45. a01 = 1.0 0.0 ;
  46. a02 = 1.0 1.0 ;
  47. a03 = 0.0 1.0 ;
  48. *
  49. la00a01 = droi 1 a00 a01 ;
  50. la01a02 = droi 1 a01 a02 ;
  51. la02a03 = droi 1 a02 a03 ;
  52. la03a00 = droi 1 a03 a00 ;
  53. *
  54. tot = daller la00a01 la01a02 la02a03 la03a00 ;
  55. **********************************************************************
  56. m_model = model tot MECANIQUE ELASTIQUE VISCOPLASTIQUE NOUAILHAS_B
  57. CONS UN ;
  58. *
  59. cst_youn = 260000.0 ;
  60. cst_nu = 0.3 ;
  61. cst_alph = 1.0e-6 ;
  62. cst_rho = 1.0e-8 ;
  63. cst_n = 7.0 ;
  64. cst_kk = 130.0 ;
  65. cst_k0 = 400.0 ;
  66. cst_alf = 0.0 ;
  67. cst_cl1 = 15000.0 ;
  68. cst_dnl1 = 300.0 ;
  69. cst_ptm1 = 0.0 ;
  70. cst_gdm1 = 0.0 ;
  71. cst_cl2 = 6000.0 ;
  72. cst_dnl2 = 100.0 ;
  73. cst_phi = 0.0 ;
  74. cst_b = 0.0 ;
  75. cst_rmax = 150.0 ;
  76. cst_br = 500.0 ;
  77. *
  78. m_mater = mater m_model YOUN cst_youn NU cst_nu ALPH cst_alph
  79. RHO cst_rho N cst_n KK cst_kk K0 cst_k0 ALF cst_alf
  80. CL1 cst_cl1 DNL1 cst_dnl1 PTM1 cst_ptm1 GDM1 cst_gdm1 CL2 cst_cl2
  81. DNL2 cst_dnl2 PHI cst_phi B cst_b RMAX cst_rmax BR cst_br ;
  82. **********************************************************************
  83. m_bloca = bloque UZ la00a01 ;
  84. *
  85. t00temps = 0.0 ;
  86. t01temps = t00temps + 0.01 ;
  87. t02temps = t01temps + 1000.00 ;
  88. *
  89. li_temps = prog t00temps t01temps t02temps ;
  90. *
  91. t00press = 0.0 ;
  92. t01press = 200.0 ;
  93. t02press = 200.0 ;
  94. *
  95. li_press = prog t00press t01press t02press ;
  96. *
  97. ev_press = evol manu 'TEMPS' li_temps 'P' li_press ;
  98. *
  99. m01press = press MASSE m_model 1.0 la02a03 ;
  100. m_charg = charg MECA m01press ev_press ;
  101. **********************************************************************
  102. p01nbpas = 2.0 ;
  103. p02nbpas = 50.0 ;
  104. *
  105. p01pas = (t01temps - t00temps) / p01nbpas ;
  106. p02pas = (t02temps - t01temps) / p02nbpas ;
  107. *
  108. p01lipas = prog (t00temps + p01pas) pas p01pas t01temps ;
  109. p02lipas = prog (t01temps + p02pas) pas p02pas t02temps ;
  110. *
  111. p_lipas = p01lipas et p02lipas ;
  112. **********************************************************************
  113. li_defvg = prog 0.0 0.1066e-3 0.1861e-3 0.3487e-3 0.4539e-3
  114. 0.5944e-3 0.6889e-3 0.7997e-3 0.8792e-3 ;
  115. li_tps = prog 0.0 25.0 50.0 125.0 200.0
  116. 350.0 500.0 750.0 1000.0 ;
  117. ev_ref = evol manu 't' li_tps 'def' li_defvg ;
  118. *dess ev_ref NCLK ;
  119. **********************************************************************
  120. p_table = table ;
  121. p_table.modele = m_model ;
  122. p_table.caracteristiques = m_mater ;
  123. p_table.blocages_mecaniques = m_bloca ;
  124. p_table.chargement = m_charg ;
  125. p_table.temps_calcules = p_lipas ;
  126. *
  127. si i_trace ;
  128. p_table.procedure_perso1 = VRAI ;
  129. p_table.depo = table ;
  130. p_table.depo.ev_ref = ev_ref ;
  131. p_table.depo.li_tp = prog ;
  132. p_table.depo.li_def = prog ;
  133. finsi ;
  134. **********************************************************************
  135. TMASAU=table;
  136. p_table . 'MES_SAUVEGARDES'=TMASAU;
  137. TMASAU .'DEFTO'=VRAI;
  138. TMASAU .'DEFIN'=VRAI;
  139. pasapas p_table ;
  140. **********************************************************************
  141. ind_1 = index (p_table.deformations_inelastiques) ;
  142. nb_1 = dime ind_1 ;
  143. v_d = extr (p_table.deformations_inelastiques.(ind_1.nb_1))
  144. EIZZ 1 1 1 ;
  145. v_d = abs(v_d) ;
  146. v_dif = abs(v_d - 0.8777e-3)/0.8777e-3 ;
  147. si (v_dif > 0.01) ;
  148. erre 5 ;
  149. sinon ;
  150. erre 0 ;
  151. finsi ;
  152. **********************************************************************
  153. fin ;
  154.  
  155.  
  156.  
  157.  
  158.  
  159.  
  160.  
  161.  

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