Télécharger la plaquette Cast3M


fr en

Télécharger nouailhas_a1.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  nouailhas_a1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. **********************************************************************
  5. *                                                                    *
  6. * fichier : nouailhas_a1.dgibi                                       *
  7. *                                                                    *
  8. **********************************************************************
  9. *                                                                    *
  10. * fait par PLG le 03-09-97                                           *
  11. *                                                                    *
  12. **********************************************************************
  13. *                                                                    *
  14. * But : ktest pour la loi viscoplastique NOUAILHAS_A                 *
  15. *                                                                    *
  16. **********************************************************************
  17. i_trace = vrai ;
  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)) EIYY 1 1 1 ;
  29. *
  30. li_t = li_t et (prog v_t) ;
  31. li_d = li_d et (prog 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. si faux ;
  38. dess ((plg.depo.ev_ref) et ev_d) NCLK ;
  39. finsi ;
  40. finproc ;
  41. **********************************************************************
  42. option dime 2 ;
  43. option elem qua8 ;
  44. option mode plan cont ;
  45. *
  46. a00 = 0.0 0.0 ;
  47. a01 = 1.0 0.0 ;
  48. a02 = 1.0 1.0 ;
  49. a03 = 0.0 1.0 ;
  50. *
  51. la00a01 = droi 1 a00 a01 ;
  52. la01a02 = droi 1 a01 a02 ;
  53. la02a03 = droi 1 a02 a03 ;
  54. la03a00 = droi 1 a03 a00 ;
  55. *
  56. tot = daller la00a01 la01a02 la02a03 la03a00 ;
  57. **********************************************************************
  58. m_model = model tot MECANIQUE ELASTIQUE VISCOPLASTIQUE NOUAILHAS_A
  59.  CONS UN ;
  60. *
  61. cst_youn = 260000.0 ;
  62. cst_nu   = 0.3 ;
  63. cst_alph = 1.0e-6 ;
  64. cst_rho  = 1.0e-8 ;
  65. cst_n    = 1.0 ;
  66. cst_kk   = 200.0 ;
  67. cst_k0   = 10000.0 ;
  68. cst_alf  = 0.0 ;
  69. cst_alfr = 0.0 ;
  70. cst_alfk = 0.0 ;
  71. *
  72. cst_cl1  = 7000.0 ;
  73. cst_dnl1 = 0.0 ;
  74. cst_ptm1 = 1.0 ;
  75. cst_gdm1 = 20000.0 ;
  76. *
  77. cst_cl2  = 0.0 ;
  78. cst_dnl2 = 0.0 ;
  79. cst_ptm2 = 0.0 ;
  80. cst_gdm2 = 0.0 ;
  81. *
  82. cst_phi  = 0.0 ;
  83. cst_b    = 0.0 ;
  84. cst_qsta = 0.0 ;
  85. cst_qmax = 1.0 ;
  86. cst_mu   = 0.0 ;
  87. cst_m    = 0.0 ;
  88. cst_gama = 0.0 ;
  89. cst_eta  = 1.0 ;
  90. *
  91. m_mater = mater m_model      YOUN cst_youn NU   cst_nu   ALPH cst_alph
  92.  RHO  cst_rho  N    cst_n    KK   cst_kk   K0   cst_k0   ALF  cst_alf
  93.  CL1  cst_cl1  DNL1 cst_dnl1 PTM1 cst_ptm1 GDM1 cst_gdm1 CL2  cst_cl2
  94.  DNL2 cst_dnl2 PTM2 cst_ptm2 GDM2 cst_gdm2 PHI  cst_phi  B    cst_b  
  95.  QSTA cst_qsta QMAX cst_qmax MU   cst_mu   GAMA cst_gama ALFK cst_alfk
  96.  ETA  cst_eta  ALFR cst_alfr M    cst_m ;
  97. **********************************************************************
  98. m01bloca = bloque UX la03a00 ;
  99. m02bloca = bloque UY la00a01 ;
  100. m_bloca = m01bloca et m02bloca ;
  101. *
  102. t00temps = 0.0 ;
  103. t01temps = t00temps +  0.001 ;
  104. t02temps = t01temps + 10.000 ;
  105. t03temps = t02temps +  0.001 ;
  106. t04temps = t03temps +  5.000 ;
  107. t05temps = t04temps +  0.001 ;
  108. t06temps = t05temps +  5.000 ;
  109. *
  110. li_temps = prog t00temps t01temps t02temps t03temps t04temps t05temps 
  111.  t06temps ;
  112. *
  113. t00press =   0.0 ;
  114. t01press = -500.0 ;
  115. t02press = -500.0 ;
  116. t03press =  -0.1 ;
  117. t04press =  -0.1 ;
  118. t05press = -500.0 ;
  119. t06press = -500.0 ;
  120. *
  121. li_press = prog t00press t01press t02press t03press t04press t05press 
  122.  t06press ;
  123. *
  124. ev_press = evol manu 'TEMPS' li_temps 'P' li_press ;
  125. *
  126. m01press = press MASSE m_model 1.0 la02a03 ;
  127. m_charg = charg MECA m01press ev_press ;
  128. **********************************************************************
  129. p01nbpas =  2.0 ;
  130. p02nbpas = 40.0 ;
  131. p03nbpas =  2.0 ;
  132. p04nbpas =  5.0 ;
  133. p05nbpas =  2.0 ;
  134. p06nbpas = 20.0 ;
  135. *
  136. p01pas   = (t01temps - t00temps) / p01nbpas ;
  137. p02pas   = (t02temps - t01temps) / p02nbpas ;
  138. p03pas   = (t03temps - t02temps) / p03nbpas ;
  139. p04pas   = (t04temps - t03temps) / p04nbpas ;
  140. p05pas   = (t05temps - t04temps) / p05nbpas ;
  141. p06pas   = (t06temps - t05temps) / p06nbpas ;
  142. *
  143. p01lipas = prog (t00temps + p01pas) pas p01pas t01temps ;
  144. p02lipas = prog (t01temps + p02pas) pas p02pas t02temps ;
  145. p03lipas = prog (t02temps + p03pas) pas p03pas t03temps ;
  146. p04lipas = prog (t03temps + p04pas) pas p04pas t04temps ;
  147. p05lipas = prog (t04temps + p05pas) pas p05pas t05temps ;
  148. p06lipas = prog (t05temps + p06pas) pas p06pas t06temps ;
  149. *
  150. p_lipas = p01lipas et p02lipas et p03lipas et p04lipas et p05lipas et
  151.  p06lipas ;
  152. **********************************************************************
  153. li_defvg = prog  0.0     0.02312 0.03720 0.04847 0.05887 0.06899
  154.  0.07903 0.08904 0.09904 0.1090  0.1190  0.1190  0.1306  0.1394 
  155.  0.1467  0.1530  0.1587  0.1642  0.1694  0.1746  0.1797  0.1847 ;
  156. *
  157. li_tps   = prog    0.0     1.001   2.001   3.001   4.001   5.001
  158.    6.0     7.0     8.0     9.0    10.0    15.0    15.5    16.0 
  159.   16.5    17.0    17.5    18.0    18.5    19.0    19.5    20.0 ;
  160. ev_ref  = evol manu 't' li_tps 'def' li_defvg ;
  161. *dess ev_ref NCLK ;
  162. **********************************************************************
  163. p_table = table ;
  164. p_table.modele                    = m_model ;
  165. p_table.caracteristiques          = m_mater ;
  166. p_table.blocages_mecaniques       = m_bloca ;
  167. p_table.chargement                = m_charg ;
  168. p_table.temps_calcules            = p_lipas ;
  169. *
  170. si i_trace ;
  171. p_table.procedure_perso1          = VRAI ;
  172. p_table.depo = table ;
  173. p_table.depo.ev_ref = ev_ref ;
  174. p_table.depo.li_tp  = prog ;
  175. p_table.depo.li_def = prog ;
  176. finsi ;
  177. **********************************************************************
  178. TMASAU=table;
  179. p_table . 'MES_SAUVEGARDES'=TMASAU;
  180. TMASAU .'DEFTO'=VRAI;
  181. TMASAU .'DEFIN'=VRAI;
  182. pasapas p_table ;
  183. **********************************************************************
  184. ind_1 = index (p_table.deformations_inelastiques) ;
  185. nb_1 = dime ind_1 ;
  186. *
  187. v_d = extr (p_table.deformations_inelastiques.(ind_1.nb_1))
  188.  EIYY 1 1 1 ;
  189. *
  190. v_d = abs(v_d) ;
  191. v_diff = abs(v_d - 0.1847) / 0.1847 ;
  192. *
  193. si (v_diff > 0.01) ;
  194.  erre 5 ;
  195. sinon ;
  196.  erre 0 ;
  197. finsi ;
  198. **********************************************************************
  199. fin ;
  200.  
  201.  
  202.  
  203.  
  204.  
  205.  
  206.  
  207.  

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