Télécharger cou31.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : cou31.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. OPTI ECHO 0;
  5. SAUT PAGE;
  6. MESS ' ';
  7. MESS ' ';
  8. MESS ' ';
  9. MESS ' TEST AVEC 2 CUBES ET UN JOINT (3D) ';
  10. MESS ' ';
  11. MESS ' LE CUBE DU HAUT EST SOUMIS A UN DEPLACEMENT IMPOSE ';
  12. MESS ' DE 0.1 DANS LA DIRECTION DE CISAILLEMENT X. ';
  13. MESS ' LE CUBE DU BAS EST BLOQUE EN X, EN Y ET EN Z. ';
  14. MESS ' LE JOINT PLASTIFIE EN CISAILLEMENT ';
  15. MESS ' LA CONTRAINTE DE CISAILLEMENT DANS LA ';
  16. MESS ' DIRECTION X ATTEINT LA COHESION, ALORS QUE LA ';
  17. MESS ' CONTRAINTE DE CISAILLEMENT DANS LA DIRECTION Y ';
  18. MESS ' ET LA CONTRAINTE NORMALE VALENT 0. ';
  19. MESS ' LE DEPLACEMENT DU JOINT DANS LA DIRECTION X EST DE ';
  20. MESS ' 0.1 ';
  21. MESS ' ';
  22. MESS ' ';
  23. MESS ' ';
  24. MESS ' B _______________ C ';
  25. MESS ' | | ';
  26. MESS ' | | ';
  27. MESS ' | CUBE 1 | ';
  28. MESS ' | | ';
  29. MESS ' | | ';
  30. MESS ' |_______________| ';
  31. MESS ' A| JOINT |H ';
  32. MESS ' G|_______________|D ';
  33. MESS ' | | ';
  34. MESS ' | | ';
  35. MESS ' | CUBE 2 | ';
  36. MESS ' | | ';
  37. MESS ' Z |_______________| ';
  38. MESS ' | Y K F ';
  39. MESS ' | / ';
  40. MESS ' | / ';
  41. MESS ' |/_______ X ';
  42. ********************************************************************
  43. *
  44. OPTION ELEM QUA4 MODE TRID ;
  45. *
  46. OEIL = -1.E1 -1.E6 1.E1 ;
  47. *
  48. * ----------DEFINITION DE LA GEOMETRIE DU CUBE SUPERIEUR----------
  49. *
  50. SA1 = 0.00 0.00 0.00 ;
  51. SB1 = 1.00 0.00 0.00 ;
  52. SC1 = 1.00 1.00 0.00 ;
  53. SD1 = 0.00 1.00 0.00 ;
  54. *
  55. * ---------- MAILLAGE ----------
  56. *
  57. SH1 = SA1 DROIT 1 SB1 ;
  58. SH2 = SB1 DROIT 1 SC1 ;
  59. SH3 = SC1 DROIT 1 SD1 ;
  60. SH4 = SD1 DROIT 1 SA1 ;
  61. *
  62. SL1 = SH1 ET SH2 ET SH3 ET SH4 ;
  63. SS1 = SURF SL1 PLAN ;
  64. *
  65. * ----------DEFINITION DE LA GEOMETRIE DU CUBE INFERIEUR----------
  66. *
  67. IA1 = 0.00 0.00 0.00 ;
  68. IB1 = 1.00 0.00 0.00 ;
  69. IC1 = 1.00 1.00 0.00 ;
  70. ID1 = 0.00 1.00 0.00 ;
  71. *
  72. * ---------- MAILLAGE ----------
  73. *
  74. IH1 = IA1 DROIT 1 IB1 ;
  75. IH2 = IB1 DROIT 1 IC1 ;
  76. IH3 = IC1 DROIT 1 ID1 ;
  77. IH4 = ID1 DROIT 1 IA1 ;
  78. *
  79. IL1 = IH1 ET IH2 ET IH3 ET IH4 ;
  80. IS1 = SURF IL1 PLAN ;
  81. *
  82. OPTION ELEM CUB8 ;
  83. *
  84. VOLS = SS1 VOLU 1 TRAN ( 0.00 0.00 1.00 ) ;
  85. VOLI = IS1 VOLU 1 TRAN ( 0.00 0.00 -1.00 ) ;
  86. FACS = FACE 2 VOLS ;
  87. FACI = FACE 2 VOLI ;
  88. VT1 = VOLS ET VOLI ;
  89. * definition pour la verif finale
  90. P3 = NOEUD 3 ;
  91. *
  92. * ----------DEFINITION DE LA GEOMETRIE DU JOINT----------
  93. *
  94. * JOINT FORME DES FACES SS1 ET IS1
  95. *
  96. * ---------- MAILLAGE ----------
  97. *
  98. OPTION ELEM LIA4 ;
  99. *
  100. JOINT = LIAISON 0.001 SS1 IS1 ;
  101. *
  102. *
  103. * ---------- MAILLAGE TOTAL ----------
  104. *
  105. GEOTOT = VT1 ET JOINT ;
  106. *
  107. * ---------- DEFINITION DES CONDITIONS AUX LIMITES ----------
  108. *
  109. CL11 = BLOQ VOLI UX ;
  110. CL12 = BLOQ VOLI UY ;
  111. CL13 = BLOQ VOLI UZ ;
  112. CL1 = CL11 ET CL12 ET CL13 ;
  113. *
  114. CL21 = BLOQ FACI UX ;
  115. CL22 = BLOQ FACI UY ;
  116. CL23 = BLOQ FACI UZ ;
  117. CL2 = CL21 ET CL22 ET CL23 ;
  118. *
  119. CL3 = BLOQ SS1 UX ;
  120. DEP = DEPI CL3 0.1 ;
  121. *
  122. CL = CL1 ET CL3 ;
  123. *
  124. * ---------- DEFINITION DU MODELE DES CUBES ----------
  125. *
  126. MODCUB = MODE VT1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' ;
  127. MACUB = MATE MODCUB YOUN 100.0 NU 0.0000001 ;
  128. *
  129. * ---------- DEFINITION DU MODELE DU JOINT ----------
  130. *
  131. MODJOI = MODE JOINT MECANIQUE ELASTIQUE PLASTIQUE COULOMB JOI4 ;
  132. MAJOI = MATE MODJOI KS 100.0 KN 100.0 EF 150.0 ECN 0.05
  133. COHE 2.0 FRIC 45.0 ;
  134. *
  135. * ---------- MODELE ET MATERIAU DE TOUTE LA STRUCTURE ----------
  136. *
  137. MODTOT = MODCUB ET MODJOI;
  138. MATTOT = MACUB ET MAJOI ;
  139. *
  140. * ---------- DEFINITION DU CHARGEMENT ----------
  141. *
  142. LI1 = PROG 0.0 1.0 2.0 ;
  143. LI2 = PROG 0.0 1.0 1.0 ;
  144. EV = EVOL MANU X LI1 Y LI2 ;
  145. *
  146. CAR1 = CHAR 'DIMP' DEP EV ;
  147. *
  148. * ---------- RESOLUTION ----------
  149. *
  150. TAB = TABLE ;
  151. LIS1 = PROG 0.0 PAS 1.0 2.0 ;
  152. TAB.'BLOCAGES_MECANIQUES' = CL;
  153. TAB.'CARACTERISTIQUES' = MATTOT;
  154. TAB.'CHARGEMENT' = CAR1;
  155. TAB.'TEMPS_CALCULES' = LIS1;
  156. TAB.'MODELE' = MODTOT;
  157.  
  158. TMASAU=table;
  159. tab . 'MES_SAUVEGARDES'=TMASAU;
  160. TMASAU .'DEFTO'=VRAI;
  161. TMASAU .'DEFIN'=VRAI;
  162. PASAPAS TAB ;
  163. *
  164. * ---------- POST TRAITEMENT ----------
  165. *
  166. *
  167. * -----RECUPERATION DE LA TABLE DES DEPLACEMENTS
  168. *
  169. dtab1= tab.temps;
  170. *
  171. ndime= (dime dtab1) - 1 ;
  172. *
  173. * ----- EXTRACTION DES DIFFERENTES TABLES
  174. *
  175. * DEPLACEMENTS
  176. *
  177. i=0 ;
  178. repeter bloc1 ndime ;
  179. i=i+1 ;
  180. dep1=tab.deplacements.i ;
  181. deft=defo geotot dep1 1. ;
  182. ** trac deft ;
  183. ** trac OEIL deft ;
  184. sig1=tab.contraintes.i ;
  185. var1=tab.variables_internes.i ;
  186. defp=tab.deformations_inelastiques.i ;
  187. eps1 = epsi modjoi dep1 ;
  188. epscu = epsi modcub dep1 ;
  189. *
  190. * ----- REDUCTION DES CHPOINTS A L'ELEMENT JOINT
  191. *
  192. *
  193. *
  194. *
  195. DPTCU1 = REDU DEP1 VOLS ;
  196. SIGCU1 = REDU SIG1 VOLS ;
  197. DEFCU1 = REDU EPSCU VOLS ;
  198. DPTCU2 = REDU DEP1 VOLI ;
  199. SIGCU2 = REDU SIG1 VOLI ;
  200. DEFCU2 = REDU EPSCU VOLI ;
  201. DPTJOI = REDU DEP1 JOINT ;
  202. SIGJOI = REDU SIG1 JOINT ;
  203. DEFJOI = REDU EPS1 JOINT ;
  204. DEPJOI = REDU DEFP JOINT ;
  205. VARJOI = REDU VAR1 JOINT ;
  206. *
  207. mess ' ' ;
  208. mess ' TEMPS ' D ;
  209. mess ' ____________ ' ;
  210. mess ' ' ;
  211. mess ' deplacements du cube 1 ' ;
  212. list dptcu1 ;
  213. mess ' ' ;
  214. mess ' deplacements du cube 2 ' ;
  215. list dptcu2 ;
  216. mess ' ' ;
  217. mess ' deplacements du joint ' ;
  218. list dptjoi ;
  219. mess ' ' ;
  220. mess ' contraintes du cube 1 ' ;
  221. list sigcu1 ;
  222. mess ' ' ;
  223. mess ' contraintes du cube 2 ' ;
  224. list sigcu2 ;
  225. mess ' ' ;
  226. mess ' contraintes du joint ' ;
  227. list sigjoi ;
  228. mess ' ' ;
  229. mess ' def totales du joint' ;
  230. list defjoi ;
  231. mess ' ' ;
  232. mess ' def inelasti. du joint ' ;
  233. list depjoi ;
  234. mess ' variables internes du joint' ;
  235. list varjoi ;
  236. mess ' defo. elasti du cube du haut ' ;
  237. list defcu1 ;
  238. mess ' defo. elasti du cube du bas ' ;
  239. list defcu2 ;
  240. *
  241. *
  242. DEP13 = EXTR DPTCU1 UX P3 ;
  243. DIF13 = DEP13-0.1;
  244. XDIF1 = ABS(DIF13);
  245. SI (XDIF1 <EG 1.E-4) ;
  246. ERRE 0 ;
  247. SINON ;
  248. ERRE 5 ;
  249. FINSI ;
  250. *
  251. DEP23 = EXTR DPTCU1 UZ P3 ;
  252. DIF23 = DEP23-0.08;
  253. ZDIF1 = ABS(DIF13);
  254. SI (ZDIF1 <EG 1.E-4) ;
  255. ERRE 0 ;
  256. SINON ;
  257. ERRE 5 ;
  258. FINSI ;
  259. *
  260. SMNJ = EXTR SIGJOI SMN 1 1 1 ;
  261. ASMN = ABS(SMNJ) ;
  262. SI (ASMN <EG 1.E-4) ;
  263. ERRE 0 ;
  264. SINON ;
  265. ERRE 5 ;
  266. FINSI ;
  267. *
  268. SMZC = EXTR SIGCU1 SMZZ 1 1 1 ;
  269. DIFZ = SMZC ;
  270. ASZC = ABS(DIFZ) ;
  271. SI (ASZC <EG 1.E-4) ;
  272. ERRE 0 ;
  273. SINON ;
  274. ERRE 5 ;
  275. FINSI ;
  276. *
  277. fin bloc1 ;
  278. *
  279. FIN;
  280.  
  281.  
  282.  
  283.  
  284.  
  285.  
  286.  
  287.  
  288.  
  289.  
  290.  
  291.  
  292.  
  293.  
  294.  
  295.  
  296.  
  297.  
  298.  
  299.  
  300.  
  301.  
  302.  
  303.  
  304.  
  305.  
  306.  
  307.  
  308.  
  309.  
  310.  
  311.  
  312.  
  313.  
  314.  
  315.  
  316.  
  317.  

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