Télécharger joi45.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : joi45.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Joi45.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7. OPTI ECHO 0 ;
  8. SAUT PAGE ;
  9. *=====================================================*
  10. * *
  11. * TEST JOI45 *
  12. * *
  13. * ESSAI DE CISAILLEMENT SUR UN JOINT 3 *
  14. * ORTHOTROPE *
  15. * *
  16. * Un joint 3D JOI4 a sa surface inferieure *
  17. * encastree. Sa surface superieure est libre. *
  18. * Un effort de cisaillement est exerce sur le *
  19. * pourtour de sa surface superieure, suivant la *
  20. * direction O-P1. *
  21. * *
  22. * Y ! *
  23. * .!. *
  24. * . ! . *
  25. * . ! . P1 *
  26. * . ! / . *
  27. * . ! / . *
  28. * . ! / . *
  29. * . O ! / . *
  30. * ------------------------------------- X *
  31. * . ! . *
  32. * . ! . *
  33. * . ! . *
  34. * . ! . *
  35. * . ! . *
  36. * . ! . *
  37. * .!. *
  38. * ! *
  39. * *
  40. * Angle (OX,OP1) = 55 degres *
  41. * L axe OP1 definit la premiere direction *
  42. * d orthotropie *
  43. * *
  44. * Calcul de la solution analytique : *
  45. * *
  46. * F *
  47. * delta = ------------- *
  48. * K * A *
  49. * ou *
  50. * *
  51. * delta = deplacement suivant la direction *
  52. * de la force *
  53. * F = Force totale exercee sur la surface *
  54. * superieure, suivant la direction OP1 *
  55. * K = raideur dans la premiere direction *
  56. * d orthotropie *
  57. * A = surface du joint *
  58. * *
  59. * Projection de la solution suivant les axes X et Y : *
  60. * *
  61. * deltax = delta * cos(55) *
  62. * *
  63. * deltay = delta * sin(55) *
  64. *=====================================================*
  65. OPTION ELEM QUA4 MODE TRID ;
  66. *
  67. OO = 00.00 0.00 0.00 ;
  68. ON = 00.00 0.00 1.00 ;
  69. AN = 10.00 ;
  70. *
  71. * -------- DEFINITION DU POINT P1 --------------
  72. * ----- POUR LES DIRECTIONS D'ORTHOTROPIE ------
  73. *
  74. P1 = 0.707106781 0.707106781 0.00 ;
  75. P1 = P1 TOUR AN OO ON ;
  76. *
  77. * ------ DEFINITION DE LA SURFACE BOT DU JOINT --------
  78. *
  79. A1 = 1.414213562 0.000000000 0.00 ;
  80. B1 = 0.000000000 1.414213562 0.00 ;
  81. C1 =-1.414213562 0.000000000 0.00 ;
  82. D1 = 0.000000000 -1.414213562 0.00 ;
  83. *
  84. A1 = A1 TOUR AN OO ON ;
  85. B1 = B1 TOUR AN OO ON ;
  86. C1 = C1 TOUR AN OO ON ;
  87. D1 = D1 TOUR AN OO ON ;
  88. *
  89. * ---------- MAILLAGE ----------
  90. *
  91. H1 = A1 DROIT 1 B1 ;
  92. H2 = B1 DROIT 1 C1 ;
  93. H3 = C1 DROIT 1 D1 ;
  94. H4 = D1 DROIT 1 A1 ;
  95. *
  96. L1 = H1 ET H2 ET H3 ET H4 ;
  97. S1 = SURF L1 PLAN ;
  98. *
  99. * ------- DEFINITION DE LA SURFACE TOP DU JOINT -------
  100. *
  101. IA1 = 1.414213562 0.000000000 0.00 ;
  102. IB1 = 0.000000000 1.414213562 0.00 ;
  103. IC1 =-1.414213562 0.000000000 0.00 ;
  104. ID1 = 0.000000000 -1.414213562 0.00 ;
  105. *
  106. IA1 = IA1 TOUR AN OO ON ;
  107. IB1 = IB1 TOUR AN OO ON ;
  108. IC1 = IC1 TOUR AN OO ON ;
  109. ID1 = ID1 TOUR AN OO ON ;
  110. *
  111. * ---------- MAILLAGE ----------
  112. *
  113. IH1 = IA1 DROIT 1 IB1 ;
  114. IH2 = IB1 DROIT 1 IC1 ;
  115. IH3 = IC1 DROIT 1 ID1 ;
  116. IH4 = ID1 DROIT 1 IA1 ;
  117. *
  118. IL1 = IH1 ET IH2 ET IH3 ET IH4 ;
  119. IS1 = SURF IL1 PLAN ;
  120. *
  121. * ---------- CREATION DU JOINT JOI4 ----------
  122. *
  123. OPTION ELEM LIA4 ;
  124. VOL = LIAISON 0.00001 S1 IS1 ;
  125. *
  126. * -------- DEFINITION DES CONDITIONS LIMITES ----------
  127. *
  128. CL11 = BLOQ A1 UX ;
  129. CL12 = BLOQ A1 UY ;
  130. CL13 = BLOQ A1 UZ ;
  131. CL1 = CL11 ET CL12 ET CL13 ;
  132. *
  133. CL21 = BLOQ B1 UX ;
  134. CL22 = BLOQ B1 UY ;
  135. CL23 = BLOQ B1 UZ ;
  136. CL2 = CL21 ET CL22 ET CL23 ;
  137. *
  138. CL31 = BLOQ C1 UX ;
  139. CL32 = BLOQ C1 UY ;
  140. CL33 = BLOQ C1 UZ ;
  141. CL3 = CL31 ET CL32 ET CL33 ;
  142. *
  143. CL41 = BLOQ D1 UX ;
  144. CL42 = BLOQ D1 UY ;
  145. CL43 = BLOQ D1 UZ ;
  146. CL4 = CL41 ET CL42 ET CL43 ;
  147. *
  148. CL = CL1 ET CL2 ET CL3 ET CL4 ;
  149. *
  150. * ------- DEFINITION DU MODELE DU JOINT ---------
  151. *
  152. MOD1 = MODE VOL 'MECANIQUE' 'ELASTIQUE' 'ORTHOTROPE'
  153. JOI4 ;
  154. MA1 = MATE MOD1 'DIRECTION' P1 'PARALLELE' KS1 4.2E10
  155. KS2 4.2E08
  156. KN 4.2E20 ;
  157. **********************************************
  158. FLOT1 = 0.0 ;
  159. P3 = 3.0 1.0 1.0 ;
  160. ***MA1 = MATE MOD1 'DIRECTION' P1 'INCLINE' FLOT1 P3
  161. *** KS1 4.2E10 KS2 4.2E08 KN 4.2E20 ;
  162. **********************************************
  163. **MA1 = MATE MOD1 'RADIAL' P1 KS1 4.2E10 KS2 4.2E08 KN
  164. * 4.2E20 ;
  165. *
  166. * ---------- MATRICE DE RIGIDITE ---------
  167. *
  168. RI1 = RIGI MOD1 MA1 ;
  169. RI2 = RI1 ET CL ;
  170.  
  171. * ---------- FORCE DE TRACTION ---------
  172. *
  173. FO1 = FORCE ( 57358.00 81915.00 0.00 ) IL1 ;
  174. *
  175. * ---------- RESOLUTION ---------
  176. *
  177. RE = RESO RI2 FO1 ;
  178. SAUT PAGE ;
  179. MESS ' ' ;
  180. MESS ' ' ;
  181. MESS ' ' ;
  182. MESS'!==============================================!';
  183. MESS'! !';
  184. MESS'! TEST JOI45 !';
  185. MESS'! !';
  186. MESS'! ESSAI DE CISAILLEMENT SUR UN JOINT 3D !';
  187. MESS'! ORTHOTROPE !';
  188. MESS'! !';
  189. MESS'! Un joint 3D JOI4 a sa surface inferieure !';
  190. MESS'! encastree. Sa surface superieure est libre. !';
  191. MESS'! Un effort de cisaillement est exerce sur le !';
  192. MESS'! pourtour de sa surface superieure, suivant !';
  193. MESS'! la direction O-P1. !';
  194. MESS'! !';
  195. MESS'! Y ! !';
  196. MESS'! .!. !';
  197. MESS'! . ! . !';
  198. MESS'! . ! . P1 !';
  199. MESS'! . ! / . !';
  200. MESS'! . ! / . !';
  201. MESS'! . ! / . !';
  202. MESS'! . O ! / . !';
  203. MESS'! ------------------------------------- X !';
  204. MESS'! . ! . !';
  205. MESS'! . ! . !';
  206. MESS'! . ! . !';
  207. MESS'! . ! . !';
  208. MESS'! . ! . !';
  209. MESS'! . ! . !';
  210. MESS'! .!. !';
  211. MESS'! ! !';
  212. MESS'! !';
  213. MESS'! Angle (OX,OP1) = 55 degres !';
  214. MESS'! L axe OP1 definit la premiere direction !';
  215. MESS'! d orthotropie !';
  216. MESS'! !';
  217. MESS'! !';
  218. MESS'! Calcul de la solution analytique : !';
  219. MESS'! !';
  220. MESS'! F !';
  221. MESS'! delta = ------------- !';
  222. MESS'! K * A !';
  223. MESS'! ou !';
  224. MESS'! !';
  225. MESS'! delta = deplacement suivant la direction !';
  226. MESS'! de la force !';
  227. MESS'! F = Force totale exercee sur la surface !';
  228. MESS'! superieure,suivant la direction OP1 !';
  229. MESS'! K = raideur dans la premiere direction !';
  230. MESS'! d orthotropie !';
  231. MESS'! A = surface du joint !';
  232. MESS'! !';
  233. MESS'! Projection de la solution suivant !';
  234. MESS'! les axes X et Y : !';
  235. MESS'! !';
  236. MESS'! deltax = delta * cos(55) !';
  237. MESS'! !';
  238. MESS'! deltay = delta * sin(55) !';
  239. MESS'! !';
  240. MESS'! Application numerique : !';
  241. MESS'! !';
  242. MESS'! F = 100000 !';
  243. MESS'! K = 4.2E20 !';
  244. MESS'! A = 4.0 !';
  245. MESS'! !';
  246. MESS'! deltax = 3.416E-7 !';
  247. MESS'! deltay = 4.875E-7 !';
  248. MESS'! !';
  249. MESS'!==============================================!';
  250. MESS ' ' ;
  251. MESS ' ' ;
  252. SAUT PAGE ;
  253. MESS ' Solution Calculee :' ;
  254. MESS ' ' ;
  255. *
  256. SOLAX=3.416E-7 ;
  257. SOLAY=4.875E-7 ;
  258. *
  259. LIST RE ;
  260. *
  261. * ---------- CODE DE FONCTIONNEMENT ----------
  262. *
  263. DEPIAX = EXTR RE UX IA1 ;
  264. RESI1 = ABS( (DEPIAX-SOLAX)/SOLAX ) ;
  265. *
  266. DEPIAY = EXTR RE UY IA1 ;
  267. RESI2 = ABS( (DEPIAY-SOLAY)/SOLAY ) ;
  268. *
  269. SI ( (RESI1 <EG 1E-3 ) ET (RESI2 <EG 1E-3 ) ) ;
  270. ERRE 0 ;
  271. SINO;
  272. ERRE 5 ;
  273. FINSI ;
  274. *
  275. * ---------- CALCUL DES DEFORMATIONS ----------
  276. *
  277. EPS1 = EPSI MOD1 MA1 RE ;
  278. LIST EPS1 ;
  279. *
  280. * ---------- CALCUL DES CONTRAINTES ----------
  281. *
  282. SIG1 = SIGMA MOD1 MA1 RE ;
  283. LIST SIG1 ;
  284. *
  285. * ---------- CALCUL DE BSIGMA ----------
  286. *
  287. BSI1 = BSIGMA MOD1 SIG1 ;
  288. LIST BSI1 ;
  289. *
  290. FIN ;
  291.  
  292.  
  293.  
  294.  
  295.  
  296.  
  297.  
  298.  
  299.  
  300.  
  301.  
  302.  

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