Télécharger cou21.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : cou21.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. OPTI ECHO 0 ;
  5. SAUT PAGE ;
  6. MESS ' ';
  7. MESS ' ';
  8. MESS ' TEST AVEC 2 CUBES ET UN JOINT (2D) ';
  9. MESS ' LE CUBE DU HAUT EST SOUMIS A UN DEPLACEMENT IMPOSE. ';
  10. MESS ' L EVOLUTION DU CHARGEMENT EST TELLE QUE LE JOINT ';
  11. MESS ' S OUVRE, PUIS SE REFERME, PUIS SE COMPRIME, ET ENFIN ';
  12. MESS ' S OUVRE DE NOUVEAU ';
  13. MESS ' LE JOINT A AU DEPART UN JEU INITIAL ';
  14. MESS ' ';
  15. MESS ' ';
  16. MESS ' ';
  17. MESS ' B _______________ C ';
  18. MESS ' | | ';
  19. MESS ' | | ';
  20. MESS ' | | ';
  21. MESS ' | | ';
  22. MESS ' | ^ ^ ^ | ';
  23. MESS ' |__|_____|____|_| ';
  24. MESS ' A| |H ';
  25. MESS ' G|_______________|D ';
  26. MESS ' | | ';
  27. MESS ' | | ';
  28. MESS ' | | ';
  29. MESS ' | | ';
  30. MESS ' |_______________| ';
  31. MESS ' K F ';
  32. MESS ' ';
  33. MESS ' JOINT FORME DES NOEUDS GDHA ';
  34. MESS ' CONDITIONS AUX LIMITES : ';
  35. MESS ' - K F BLOQUES EN UX ET EN UY ';
  36. MESS ' - B C BLOQUES EN UX ET EN UY ';
  37. MESS ' CHARGEMENT : ';
  38. MESS ' - DEPLACEMENT IMPOSE SUR AH ';
  39. MESS ' ';
  40. ********************************************************************
  41. *
  42. OPTION ELEM SEG2 MODE PLAN CONT ;
  43. *
  44. * ----------DEFINITION DE LA GEOMETRIE DU CUBE ----------
  45. *
  46. A = 0.00 0.30 ;
  47. B = 0.00 1.30;
  48. C = 4.00 1.30 ;
  49. H = 4.00 0.30 ;
  50. K = 0.00 -1.00 ;
  51. F = 4.00 -1.00 ;
  52. *
  53. G = 0.00 0.00 ;
  54. D = 4.00 0.00 ;
  55. *
  56. * ---------- LIGNE ----------
  57. *
  58. AB = A DROIT 1 B ;
  59. BC = B DROIT 1 C ;
  60. CH = C DROIT 1 H ;
  61. HA = H DROIT 1 A ;
  62. AH = A DROIT 1 H ;
  63. *
  64. DG = D DROIT 1 G ;
  65. DF = D DROIT 1 F ;
  66. FK = F DROIT 1 K ;
  67. KG = K DROIT 1 G ;
  68. *
  69. * ---------- MAILLAGE DU CUBE ----------
  70. *
  71. OPTION ELEM QUA4 MODE PLAN CONT ;
  72. *
  73. CUB1 = DALLER AB BC CH HA PLAN ;
  74. CUB2 = DALLER KG DG DF FK PLAN ;
  75. * quelques definitions pour la verif finale
  76. P4 = NOEUD 4 ;
  77. P1 = NOEUD 1 ;
  78. P8 = NOEUD 8 ;
  79. P7 = NOEUD 7 ;
  80. *
  81. * ----------DEFINITION DE LA GEOMETRIE DU JOINT----------
  82. *
  83. * JOINT FORME DES NOEUDS H A D G
  84. *
  85. * ---------- MAILLAGE ----------
  86. *
  87. OPTION ELEM RAC2 ;
  88. *
  89. JOINT = RACC 0.3001 HA DG ;
  90. JOINT = RACC 0.3001 (INVE HA) ( INVE DG) ;
  91. *
  92. * ---------- MAILLAGE TOTAL ----------
  93. *
  94. CUBE = CUB1 ET CUB2 ;
  95. GEOTOT = CUBE ET JOINT ;
  96. *trac geotot ;
  97. *
  98. * ---------- DEFINITION DES CONDITIONS AUX LIMITES ----------
  99. *
  100. *
  101. *
  102. * CUBE
  103. *
  104. CLBX = BLOQ B UX ;
  105. CLBY = BLOQ B UY ;
  106. CLB = CLBX ET CLBY ;
  107. *
  108. CLCX = BLOQ C UX ;
  109. CLCY = BLOQ C UY ;
  110. CLC = CLCX ET CLCY ;
  111. *
  112. CLKY = BLOQ K UY ;
  113. CLKX = BLOQ K UX ;
  114. CLK = CLKY ET CLKX ;
  115. *
  116. CLFY = BLOQ F UY ;
  117. CLFX = BLOQ F UX ;
  118. CLF = CLFY ET CLFX ;
  119. *
  120. cl5 = bloq ha uy ;
  121. dep5 = depi cl5 1.0 ;
  122. *
  123. CLCUBE = CLB ET CLC ET CLF ET CLK ET cl5 ;
  124. *
  125. * ---------- CONDITIONS AUX LIMITES TOTALES ----------
  126. *
  127. CLTOT = CLCUBE ;
  128. *
  129. * ---------- DEFINITION DU MODELE DES CUBES ----------
  130. *
  131. MODCUB = MODE CUBE 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' ;
  132. MACUB = MATE MODCUB YOUN 100.0 NU 0.000001 ;
  133. *
  134. * ---------- DEFINITION DU MODELE DU JOINT ----------
  135. *
  136. MODJOI = MODE JOINT MECANIQUE ELASTIQUE PLASTIQUE COULOMB JOI2 ;
  137. MAJOI = MATE MODJOI KS 1000.0 KN 100.0 EF 150.0 ECN 0.05
  138. COHE 2.0 FRIC 45.0 ;
  139. *
  140. * ---------- MODELE ET MATERIAU DE TOUTE LA STRUCTURE ----------
  141. *
  142. MODTOT = MODCUB ET MODJOI;
  143. MATTOT = MACUB ET MAJOI ;
  144. *
  145. * ---------- DEFINITION DU CHARGEMENT ----------
  146. *
  147. L51 = PROG 0.0 1.0 2.0 3.0 4.0 ;
  148. L52 = PROG 0.0 1.0 0.0 -0.5 1.0 ;
  149. EV5 = EVOL MANU X L51 Y L52 ;
  150. CAR = CHAR 'DIMP' DEP5 EV5;
  151. *
  152. * ---------- RESOLUTION ----------
  153. *
  154. TAB = TABLE ;
  155. LIS1 = PROG 0.0 PAS 0.5 4.0 ;
  156. DEFH = (MANU CHML MODCUB EIXX 0.0 EIYY 0.0 EIZZ 0.0 GIXY 0.0
  157. STRESSES TYPE 'DEFORMATIONS INELASTIQUES') ET
  158. (MANU CHML MODJOI FRSN 0.0 FRN 0.3
  159. STRESSES TYPE 'DEFORMATIONS INELASTIQUES') ;
  160. TAB.'DEFORMATIONS_INELASTIQUES' = TABLE;
  161. TAB.'DEFORMATIONS_INELASTIQUES' . 0 = DEFH;
  162. TAB.'BLOCAGES_MECANIQUES' = CLTOT;
  163. TAB.'MODELE' = MODTOT;
  164. TAB.'CHARGEMENT' = CAR;
  165. TAB.'TEMPS_CALCULES' = LIS1;
  166. TAB.'CARACTERISTIQUES' = MATTOT;
  167. TAB.'HYPOTHESE_DEFORMATIONS' = 'LINEAIRE' ;
  168. TMASAU=table;
  169. tab . 'MES_SAUVEGARDES'=TMASAU;
  170. TMASAU .'DEFTO'=VRAI;
  171. TMASAU .'DEFIN'=VRAI;
  172. PASAPAS TAB ;
  173. *
  174. * ---------- POST TRAITEMENT ----------
  175. *
  176. *
  177. * -----RECUPERATION DE LA TABLE DES DEPLACEMENTS
  178. *
  179. dtab1=tab.temps;
  180. *
  181. ndime=dime dtab1 ;
  182. *
  183. * ----- EXTRACTION DES DIFFERENTES TABLES
  184. *
  185. * DEPLACEMENTS
  186. *
  187. i=0 ;
  188. REPETER BOU1 (ndime-1) ;
  189. i=i+1 ;
  190. dep1=tab.deplacements.i ;
  191. sig1=tab.contraintes.i ;
  192. var1=tab.variables_internes.i ;
  193. defp=tab.deformations_inelastiques.i ;
  194. eps1 = epsi modjoi dep1 ;
  195. *
  196. * ----- REDUCTION DES CHPOINTS A L'ELEMENT JOINT
  197. *
  198. *
  199. *
  200. *
  201. DPTCU1 = REDU DEP1 CUB1 ;
  202. SIGCU1 = REDU SIG1 CUB1 ;
  203. DPTCU2 = REDU DEP1 CUB2 ;
  204. SIGCU2 = REDU SIG1 CUB2 ;
  205. DPTJOI = REDU DEP1 JOINT ;
  206. SIGJOI = REDU SIG1 JOINT ;
  207. DEFJOI = REDU EPS1 JOINT ;
  208. DEPJOI = REDU DEFP JOINT ;
  209. VARJOI = REDU VAR1 JOINT ;
  210. *
  211. mess ' ' ;
  212. mess ' PAS DE TEMPS ' i ;
  213. mess ' ____________ ' ;
  214. mess ' ' ;
  215. mess ' deplacements du cube 1 ' ;
  216. *list dptcu1 ;
  217. mess ' ' ;
  218. mess ' deplacements du cube 2 ' ;
  219. *list dptcu2 ;
  220. mess ' ' ;
  221. mess ' deplacements du joint ' ;
  222. *list dptjoi ;
  223. mess ' ' ;
  224. mess ' contraintes du cube 1 ' ;
  225. *list sigcu1 ;
  226. mess ' ' ;
  227. mess ' contraintes du cube 2 ' ;
  228. *list sigcu2 ;
  229. mess ' ' ;
  230. mess ' contraintes du joint ' ;
  231. *list sigjoi ;
  232. mess ' ' ;
  233. mess ' def totales' ;
  234. *list defjoi ;
  235. mess ' ' ;
  236. mess ' def inelasti. ' ;
  237. *list depjoi ;
  238. mess ' variables internes ' ;
  239. *list varjoi ;
  240. *
  241. fin bou1 ;
  242. *
  243. DEP11 = EXTR DPTCU1 UY P1 ;
  244. DEP14 = EXTR DPTCU1 UY P4 ;
  245. DEP28 = EXTR DPTCU2 UY P8 ;
  246. DEP27 = EXTR DPTCU2 UY P7 ;
  247. *
  248. SMY1 = EXTR SIGCU1 SMYY 1 1 1 ;
  249. *
  250. SMNJ = EXTR SIGJOI SMN 1 1 1 ;
  251. DEPJ = EXTR DEPJOI FRN 1 1 1 ;
  252. STAJ = EXTR VARJOI STAT 1 1 1 ;
  253.  
  254. *
  255. SI (EGA I 8) ;
  256. *
  257. * DEPLACEMENTS DU CUBE 1
  258. *
  259. DIF11 = DEP11-1.0 ;
  260. DIF14 = DEP14-1.0 ;
  261. T11 = ABS(DIF11) ;
  262. T14 = ABS(DIF14) ;
  263. SI (( T11 <EG 1.E-4 ) ET ( T14 <EG 1.E-4 ));
  264. ERRE 0 ;
  265. SINON ;
  266. ERRE 5 ;
  267. FINSI ;
  268. *
  269. T28 = ABS(DEP28) ;
  270. T27 = ABS(DEP27) ;
  271. SI ( ( T28 <EG 1.E-4 ) ET ( T27 <EG 1.E-4 ));
  272. ERRE 0 ;
  273. SINON ;
  274. ERRE 5 ;
  275. FINSI ;
  276. *
  277. * CONTRAINTES DU CUBE 1
  278. *
  279. DIS1 = SMY1+1.E2 ;
  280. TI1 = ABS(DIS1) ;
  281. SI ( TI1 <EG 1.E-4 ) ;
  282. ERRE 0 ;
  283. SINON ;
  284. ERRE 5 ;
  285. FINSI ;
  286. *
  287. * CONTRAINTES DU JOINT
  288. *
  289. DIJ1 = DEPJ-1.3 ;
  290. DIA1 = STAJ-1.0 ;
  291. TIM1 = ABS(SMNJ) ;
  292. TIJ1 = ABS(DIJ1) ;
  293. TIA1 = ABS(DIA1) ;
  294. SI ( ( TIM1 <EG 1.E-4 ) ET ( TIJ1 <EG 1.E-4 )
  295. ET ( TIA1 <EG 1.E-4 ) ) ;
  296. ERRE 0 ;
  297. SINON ;
  298. ERRE 5 ;
  299. FINSI ;
  300. *
  301. FINSI ;
  302. fin;
  303.  
  304.  
  305.  
  306.  
  307.  
  308.  
  309.  
  310.  
  311.  
  312.  
  313.  
  314.  
  315.  
  316.  
  317.  
  318.  
  319.  
  320.  
  321.  
  322.  
  323.  

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