Télécharger la plaquette Cast3M


fr en

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