Télécharger sste2_bcn.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : sste2_bcn.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * TRIAXIAL TEST WITH A NON-HOMOGENEOUS SAMPLE
  5. * TEST: MRS-Lade model
  6. *
  7. *
  8. *------------- OPCIONES GENERALES --------------------------------
  9. *
  10. GRAPH = 'N' ;
  11. OPTION DIME 2 TRAC x ELEM qua4 MODE axis;
  12. * ajout de option epsilon lineaire pour la precision des test!
  13. OPTION epsilon lineaire;
  14. *
  15. precomp = 600.D0;
  16. valdesp = -0.05 * 1.5;
  17. *
  18. npoi = 10;
  19. npoi2 = 10;
  20. *
  21. *------------- CREACION DE LA GEOMETRIA -------------------------
  22. *
  23. p1 = 0. 0.;
  24. p2 = 1. 0.;
  25. p3 = 1. 1.5;
  26. p4 = 0. 1.5;
  27. l12 = d npoi p1 p2;
  28. l23 = d npoi2 p2 p3;
  29. l34 = d npoi p3 p4;
  30. l41 = d npoi2 p4 p1;
  31. malla = daller l12 l23 l34 l41;
  32. mall1 = malla elem 1;
  33. da = 1. / (flot npoi);
  34. db = 1.5 / (flot npoi2);
  35. aa = da 0.;
  36. bb = da db;
  37. cc = 0. db;
  38. la2 = d 1 aa p2;
  39. l4c = d 1 p4 cc;
  40. lcb = d 1 cc bb;
  41. lba = d 1 bb aa;
  42. contaux = la2 et l23 et l34 et l4c et lcb et lba;
  43. mall2 = malla inclus contaux 'BARY';
  44. *
  45. *-------------- MODELO ------------------------------------------
  46. *
  47. ROHOR = 2.5D3;
  48. E_ELAS = 146.D3;
  49. NUHOR = 0.2D0;
  50. PCC = 0.D0;
  51. PAA = 1.D0;
  52. QAA = 1.D0;
  53. EXM = 7.423D-2;
  54. EWW = 0.7D0;
  55. XK1 = 0.2D0;
  56. XK2 = 0.7256D0;
  57. ETB = 1.12D0;
  58. EXV = 1.15D0;
  59. XEP = 7.5D-1;
  60. XNN = -0.5D0;
  61. CCO = 4.D-3;
  62. EXL = 1.0867654D0;
  63. PCO = 5000.D0;
  64. EXR = 1.592D0;
  65. CCA = 0.000159D0;
  66. PHIs = 0.D0;
  67. ALPas = 0.8D0;
  68. MOD1b = modeli MALL2 mecanique elastique plastique mrs_lade;
  69. MAT1b = mater mod1b youn E_ELAS nu NUHOR rho ROHOR pc PCC pa PAA
  70. qa QAA expm EXM e EWW k1 XK1 k2 XK2 etab ETB
  71. expv EXV epsi XEP n XNN ccon CCO expl EXL
  72. pcap PCO expr EXR ccap CCA phi PHIs alp ALPas;
  73. *
  74. etb2 = etb*0.9D0;
  75. pco2 = pco*0.9D0;
  76. MOD1a = modeli MALL1 mecanique elastique plastique mrs_lade;
  77. MAT1a = mater mod1a youn E_ELAS nu NUHOR rho ROHOR pc PCC pa PAA
  78. qa QAA expm EXM e EWW k1 XK1 k2 XK2 etab ETB2
  79. expv EXV epsi XEP n XNN ccon CCO expl EXL
  80. pcap PCO2 expr EXR ccap CCA phi PHIs alp ALPas;
  81. modi = mod1a et mod1b;
  82. mati = mat1a et mat1b;
  83. *
  84. *-------------- PROBLEMA: PRECOMPRESION -------------------------
  85. *
  86. RIGconti = (bloq uz l12) et (bloq ur l41);
  87. valorp = pres mass modi precomp (l23 et l34);
  88. evolu0 = EVOL MANU t (PROG 0. 1.) level (PROG 0. 1.);
  89. evfut0 = CHAR DIMP evolu0 valorp;
  90. maxiter = 10;
  91. ti = TABLE ;
  92. ti.MODELE = MODi;
  93. ti.CARACTERISTIQUES = MATi;
  94. ti.BLOCAGES_MECANIQUES = RIGconti;
  95. ti.CHARGEMENT = evfut0;
  96. ti.ACCELERATION = maxiter;
  97. ti.MAXITERATION = maxiter;
  98. ti.PRECISION = 1.E-8;
  99. ti.K_TANGENT = VRAI;
  100. ti.TEMPS_CALCULES = PROG 0. 1.;
  101. ti.TEMPS_SAUVES = PROG 0. 1.;
  102. PASAPAS ti ;
  103. desppre = ti.deplacements.1;
  104. form desppre;
  105. *
  106. *-------------- PROBLEMA: TRIAXIAL ------------------------------
  107. *
  108. mod1 = mod1a et mod1b;
  109. mat1 = mat1a et mat1b;
  110. RIGcont = (bloq uz l12) et (bloq ur l41);
  111. *
  112. *-------------- MOVIMIENTO IMPUESTO -----------------------------
  113. *
  114. rigsupe = (bloq uz l34);
  115. movyy = (depi rigsupe valdesp);
  116. evolu2 = EVOL MANU t (PROG 0. 1.) level (PROG 0. 1.);
  117. evfut2 = CHAR DIMP evolu2 movyy;
  118. valorp = pres mass mod1 precomp l23;
  119. evolu1 = EVOL MANU t (PROG 0. 1.) level (PROG 1. 1.);
  120. evfut1 = CHAR DIMP evolu1 valorp;
  121. *
  122. *-------------- RESOLUCION --------------------------------------
  123. *
  124. tiempos = prog 0. pas 0.04 0.48 pas 0.02 1.;
  125. maxiter = 15;
  126. t0 = TABLE ;
  127. t0.CONTRAINTES = TABLE ;
  128. t0.CONTRAINTES.0 = ti.CONTRAINTES.1;
  129. t0.MODELE = MOD1;
  130. t0.CARACTERISTIQUES = MAT1;
  131. t0.BLOCAGES_MECANIQUES = RIGcont et rigsupe;
  132. t0.CHARGEMENT = evfut2 et evfut1;
  133. t0.K_TANGENT = VRAI;
  134. t0.ACCELERATION = maxiter;
  135. t0.MAXITERATION = maxiter;
  136. t0.PRECISION = 1.E-10;
  137. t0.SUBSTEPPING = vrai;
  138. t0.NMAXSUBSTEPS = 1000;
  139. t0.NITERINTER_MAX = 8;
  140. t0.PRECISINTER = 1.E-10;
  141. t0.TEMPS_CALCULES = tiempos;
  142. t0.TEMPS_SAUVES = tiempos;
  143. PASAPAS t0 ;
  144. *
  145. *-------------- POSTPROCESO -------------------------------------
  146. *
  147. imax = 38;
  148. fuer = prog 0.;
  149. i = 1;
  150. repeter blocdefi imax;
  151. aux = reac rigsupe t0.deplacements.i;
  152. aux = resu aux;
  153. pbas = extr aux MAIL ;
  154. pbas = POIN 1 pbas ;
  155. aux = extr aux fz pbas;
  156. aux = -1.D0 * aux;
  157. fuer = fuer et (prog aux);
  158. i=i+1;
  159. fin blocdefi;
  160. law = evol manu t (-1.0*(tiempos*valdesp)) f(t) fuer;
  161. SI (NEG GRAPH 'N') ;
  162. dessin law;
  163. FINSI ;
  164. err = (aux - 3.141986E+03)/3.141986E+03 ;
  165. list (aux - 3.141986E+03);
  166. list err;
  167. SI ((ABS err) < 0.5e-5) ;
  168. ERRE 0 ;
  169. SINON ;
  170. ERRE 5 ;
  171. FINSI ;
  172. FIN;
  173.  
  174.  
  175.  
  176.  
  177.  
  178.  

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