Télécharger joi_eli.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : joi_eli.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. GRAPH = 'N' ;
  5. SAUT PAGE ;
  6. ********************************************************************
  7. MESS ' ';
  8. MESS ' ';
  9. MESS ' ';
  10. MESS ' ';
  11. MESS ' TEST AVEC UN JOINT (2D) POUR LE MODELE D ANCRAGE ';
  12. MESS ' D ELIGEHAUSEN ';
  13. MESS ' ';
  14. MESS ' LA SURFACE DE HAUT DU JOINT EST SOUMISE A UN DEPL. IMPOSE ';
  15. MESS ' SELON LA DIRECTION TANGENTIELLE ';
  16. MESS ' ';
  17. MESS ' ';
  18. MESS ' ';
  19. MESS ' ';
  20. MESS ' ';
  21. MESS ' ';
  22. MESS ' A _______________ H ';
  23. MESS ' | | ';
  24. MESS ' |_______________| ';
  25. MESS ' G D ';
  26. MESS ' ';
  27. MESS ' ';
  28. MESS ' ';
  29. MESS ' ';
  30. MESS ' ';
  31. MESS ' JOINT FORME DES NOEUDS DG-HA ';
  32. MESS ' ';
  33. MESS ' CONDITIONS AUX LIMITES : ';
  34. MESS ' - DG BLOQUEE EN UX ET EN UY ';
  35. MESS ' ';
  36. MESS ' CHARGEMENT : ';
  37. MESS ' - DEPLACEMENT IMPOSE ET FORCE SUR HA ';
  38. MESS ' ';
  39. MESS ' ';
  40. MESS ' ';
  41. MESS ' DANS CET ESSAI ON COMPARE LA CONTRAINTE TANGENTIELLE ET ';
  42. MESS ' LE DEPLACEMENT NORMAL AUX QUANTITES CORRESPONDANTES ';
  43. MESS ' CALCULEES ANALYTIQUEMENT. ';
  44. MESS ' ';
  45. MESS ' ';
  46. MESS ' ';
  47. ********************************************************************
  48. *
  49. *
  50. OPTION ECHO 1 DIME 2 ;
  51. OPTION ELEM SEG2 MODE PLAN CONT ;
  52. TEMPS ;
  53. *
  54. *------------------------------------------------
  55. * --------- DEFINITION DE LA GEOMETRIE ----------
  56. *------------------------------------------------
  57. *
  58. A = 0.00 0.00 ;
  59. H = 1.00 0.00 ;
  60. *
  61. G = 0.00 0.00 ;
  62. D = 1.00 0.00 ;
  63. *
  64. HA = H DROIT 1 A ;
  65. *
  66. DG = D DROIT 1 G ;
  67. *
  68. *-----------------------------------------------------
  69. * -------------- DEFINITION DU MAILLAGE --------------
  70. *-----------------------------------------------------
  71. *
  72. OPTION ELEM RAC2 ;
  73. *
  74. JOINT = RACC 0.001 (INVE HA ) ( INVE DG ) ;
  75. *
  76. *------------------------------------------------------------------
  77. * -------------- DEFINITION DES CONDITIONS AUX LIMITES ------------
  78. *------------------------------------------------------------------
  79. *
  80. CL1 = BLOQ DG UX UY ;
  81. *
  82. CL5 = BLOQ HA UX ;
  83. *
  84. CLTOT = CL1 ET CL5 ;
  85. *
  86. *-----------------------------------------------------------------
  87. * -------------- DEFINITION DU MODELE E DU MATERIAU --------------
  88. *-----------------------------------------------------------------
  89. *
  90. MODJOI = MODE JOINT 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE'
  91. 'ANCRAGE_ELIGEHAUSEN' JOI2 ;
  92. *
  93. phi0 = 0.010;
  94. *
  95. MAJOI = MATE MODJOI 'KS ' ((30000./(2.*(1+0.25)*0.05))*pi*phi0)
  96. 'KN ' ((30000./0.05)*pi*phi0)
  97. 'S1T' 0.001 'S2T' 0.003 'S3T' 0.010
  98. 'T1T' 15. 'T3T' 5.
  99. 'ALFA' 0.4 'PERI' (pi*phi0);
  100.  
  101. *--------------------------------------------------------------
  102. * ------------------- DEFINITION DU CHARGEMENT ----------------
  103. *--------------------------------------------------------------
  104. *
  105. DEP6 = FORCE FY -10.0 HA ;
  106. DEP5 = DEPI CL5 (-20.0/1000.) ;
  107. *
  108. L56 = PROG 0.0 1.0 1.0 ;
  109. L55 = PROG 0.0 0.0 1.0 ;
  110. L50 = PROG 0.0 1.0 2.0 ;
  111. *
  112. EV5 = EVOL MANU X L50 Y L55 ;
  113. EV6 = EVOL MANU X L50 Y L56 ;
  114. *
  115. CAR5 = CHAR 'DIMP' DEP5 EV5 ;
  116. CAR6 = CHAR 'MECA' DEP6 EV6 ;
  117. CART = CAR5 ET CAR6 ;
  118. *
  119. *-------------------------------------------------------
  120. * --------------------- RESOLUTION ---------------------
  121. *-------------------------------------------------------
  122. *
  123. TAB = TABLE ;
  124. LIS1 = PROG 0.0 1.0 PAS 0.05 2.0 ;
  125. TAB.'BLOCAGES_MECANIQUES' = CLTOT;
  126. TAB.'MODELE' = MODJOI;
  127. TAB.'CHARGEMENT' = CART;
  128. TAB.'TEMPS_CALCULES' = LIS1;
  129. TAB.'CARACTERISTIQUES' = MAJOI;
  130. TAB.'MOVA' = RIEN;
  131. *
  132. TMASAU=table;
  133. tab. 'MES_SAUVEGARDES'=TMASAU;
  134. TMASAU .'DEFTO'=VRAI;
  135. TMASAU .'DEFIN'=VRAI;
  136. PASAPAS TAB ;
  137. *
  138. *--------------------------------------
  139. * ---------- POST TRAITEMENT ----------
  140. *--------------------------------------
  141. *
  142. sig1 = prog ; sig2 = prog ; var1 = prog ; var2 = prog ;
  143. var3 = prog ; var4 = prog ; var5 = prog ; var6 = prog ;
  144. dep1 = prog ; dep2 = prog ; tem0 = prog ; tem1 = prog ;
  145. *
  146. dtab1=index(tab.deplacements) ;
  147. *
  148. ndime=dime dtab1 ;
  149. *
  150. *
  151. i=0 ;
  152. *
  153. REPETER BOU1 ndime ;
  154. *
  155. i=i+1 ;
  156. *
  157. d=dtab1.i ;
  158. *
  159. dep0 = tab.deplacements.d ;
  160. sig0 = tab.contraintes.d ;
  161. var0 = tab.variables_internes.d ;
  162. def0 = tab.deformations_inelastiques.d ;
  163. *
  164. DEPJOI = REDU DEP0 JOINT ;
  165. SIGJOI = REDU SIG0 JOINT ;
  166. VARJOI = REDU VAR0 JOINT ;
  167. DEFJOI = REDU DEF0 JOINT ;
  168. *
  169. * Contrainte de cisaillement
  170. tau1 = extr SIGJOI 'SMSN' 1 1 1 ;
  171. * Contrainte normale
  172. tau2 = extr SIGJOI 'SMN' 1 1 1 ;
  173. *
  174. * Glissement plastique atteint pour un glissement positif
  175. sor1 = extr VARJOI 'EPPP' 1 1 1 ;
  176. * Glissement plastique atteint pour un glissement négatif
  177. sor2 = extr VARJOI 'EPPM' 1 1 1 ;
  178. * Glissement total
  179. sor3 = extr VARJOI 'EPSO' 1 1 1 ;
  180. * Rigidité tangente en glissement
  181. sor4 = extr VARJOI 'KSTA' 1 1 1 ;
  182. * Rigidité tangente en traction/compression
  183. sor5 = extr VARJOI 'KNEL' 1 1 1 ;
  184. *
  185. * Déformations plastiques
  186. pla1 = extr DEFJOI 'FRSN' 1 1 1 ;
  187. pla2 = extr DEFJOI 'FRN' 1 1 1 ;
  188. *
  189. sig1 = inse sig1 i tau1 ;
  190. sig2 = inse sig2 i tau2 ;
  191. *
  192. var1 = inse var1 i sor1 ;
  193. var2 = inse var2 i sor2 ;
  194. var3 = inse var3 i sor3 ;
  195. var4 = inse var4 i sor4 ;
  196. var5 = inse var5 i sor5 ;
  197. *
  198. tem0 = inse tem0 i d ;
  199. *
  200. dep1 = inse dep1 i pla1 ;
  201. dep2 = inse dep2 i pla2 ;
  202. *
  203. fin bou1 ;
  204. *
  205. si ( NEG graph 'N' ) ;
  206. *
  207. gra1 = evol manu U var3 T sig1 ;
  208. gra2 = evol manu t0 tem0 EPPP var1 ;
  209. gra3 = evol manu t0 tem0 EPPM var2 ;
  210. gra4 = evol manu t0 tem0 EPSO var3 ;
  211. gra5 = evol manu t0 tem0 KSTA var4 ;
  212. gra6 = evol manu t0 tem0 KNEL var5 ;
  213. *
  214. dess gra1 ;
  215. dess gra2 ;
  216. dess gra3 ;
  217. dess gra4 ;
  218. dess gra5 ;
  219. dess gra6 ;
  220. *
  221. finsi ;
  222. *
  223. *-------------------------------------------------------------
  224. *------------------ CODE DE FONCTIONNEMENT -------------------
  225. *-------------------------------------------------------------
  226. *
  227. *
  228. * listreel de la contrainte tangentielle de repère en fonction du temps *
  229. *
  230. pro1 = prog
  231. 0.0 0.0 -.47124 -.47124 -.47124 -.42636 -.38148
  232. -.33660 -.29172 -.24684 -.20196 -.15708 -.15708 -.15708 -.15708
  233. -.15708 -.15708 -.15708 -.15708 -.15708 -.15708 -.15708;
  234. *
  235. *
  236. ***** listreel du glissement total en fonction du temps *****
  237. *
  238. pro2 = prog 0.0 0.0 -0.001 -0.002 -0.003 -0.004 -0.005 -0.006 -0.007 -0.008
  239. -0.009 -0.01 -0.011 -0.012 -0.013 -0.014 -0.015 -0.016 -0.017 -0.018
  240. -0.019 -0.020;
  241. *
  242. j1 = 0 ;
  243. j2 = 0 ;
  244. ii0 = 0 ;
  245. num1 = 0. ;
  246. num2 = 0. ;
  247. den1 = 0. ;
  248. den2 = 0. ;
  249. *
  250. repeter bou2 20 ;
  251. *
  252. j1 = j1 + 1 ;
  253. j2 = j2 + 1 ;
  254. ii0 = ii0 + 1 ;
  255. *
  256. ttt1 = extr sig1 j1 ;
  257. vvv1 = extr var3 j2 ;
  258. ttt0 = extr pro1 ii0 ;
  259. vvv0 = extr pro2 ii0 ;
  260. *
  261. num1 = num1 + ((ttt1 - ttt0) * (ttt1 - ttt0)) ;
  262. num2 = num2 + ((vvv1 - vvv0) * (vvv1 - vvv0)) ;
  263. den1 = den1 + (ttt0 * ttt0) ;
  264. den2 = den2 + (vvv0 * vvv0) ;
  265. *
  266. fin bou2 ;
  267. *
  268. err0 = (num1/den1) * 100. ;
  269. err1 = (num2/den2) * 100. ;
  270. *
  271. si ((err0 < 5.) et (err1 < 5.)) ;
  272. *
  273. ERRE 0 ;
  274. *
  275. sinon ;
  276. *
  277. ERRE 5 ;
  278. *
  279. finsi ;
  280. *
  281. fin;
  282.  
  283.  
  284.  
  285.  
  286.  
  287.  
  288.  
  289.  
  290.  
  291.  
  292.  
  293.  
  294.  
  295.  
  296.  

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