Télécharger joi41.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : joi41.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Joi41.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7. OPTI ECHO 0;
  8. SAUT PAGE;
  9. *
  10. *******************************************************
  11. * *
  12. * TEST JOI41 *
  13. * *
  14. * 2 CUBES SUPERPOSES AVEC UN JOINT JOI4 *
  15. * AU MILIEU ( JOI4 3D ISOTROPE ) *
  16. * *
  17. * Deux CUB8 sont superposes. On place un joint JOI4 *
  18. * entre les 2 cubes. Les noeuds de la surface infe- *
  19. * -rieure du cube du bas sont encastres. Une force *
  20. * de traction verticale (i.e. suivant z) est exercee *
  21. * sur la surface superieure du cube du haut. *
  22. * Deux cas sont etudies : *
  23. * - 1 / On suppose que le joint n'existe pas entre *
  24. * les deux cubes. *
  25. * - 2 / On suppose que le joint existe et a une *
  26. * raideur normale tres grande par rapport a *
  27. * celle de chacun des cubes. *
  28. * Les deplacements UZ des noeuds des deux cubes *
  29. * doivent etre les memes dans les deux cas. *
  30. * Raideur de chacun des cubes : *
  31. * sigma = force / A A = section *
  32. * or force = k * delta k = raideur *
  33. * delta = deplacement *
  34. * de + sigma = E * epsilon E = Young *
  35. * ou epsilon = delta / L L = hauteur du cube *
  36. * ---> k = E * A / L *
  37. * ---> k = 4.2 E11 *
  38. * On choisit pour raideur du joint, la valeur *
  39. * k = 4.2 E15 *
  40. * *
  41. *******************************************************
  42. *
  43. TEMPS ;
  44. *
  45. SAUT PAGE;
  46. MESS ' ' ;
  47. MESS ' _____________________________________________' ;
  48. MESS '! !' ;
  49. MESS '! CAS DE DEUX CUBES SUPERPOSES SANS JOINT !' ;
  50. MESS '! ENTRE EUX !' ;
  51. MESS '!___________________________________________ !' ;
  52. MESS ' ' ;
  53. *
  54. OPTION ELEM QUA4 MODE TRID ;
  55. *
  56. * ----DEFINITION DE LA GEOMETRIE DU CUBE SUPERIEUR-----
  57. *
  58. SA1 = 0.00 0.00 0.00 ;
  59. SB1 = 2.00 0.00 0.00 ;
  60. SC1 = 2.00 2.00 0.00 ;
  61. SD1 = 0.00 2.00 0.00 ;
  62. *
  63. * ---------- MAILLAGE ----------
  64. *
  65. SH1 = SA1 DROIT 1 SB1 ;
  66. SH2 = SB1 DROIT 1 SC1 ;
  67. SH3 = SC1 DROIT 1 SD1 ;
  68. SH4 = SD1 DROIT 1 SA1 ;
  69. *
  70. SL1 = SH1 ET SH2 ET SH3 ET SH4 ;
  71. SS1 = SURF SL1 PLAN ;
  72. *
  73. * ----DEFINITION DE LA GEOMETRIE DU CUBE INFERIEUR-----
  74. *
  75. IA1 = 0.00 0.00 -4.00 ;
  76. IB1 = 2.00 0.00 -4.00 ;
  77. IC1 = 2.00 2.00 -4.00 ;
  78. ID1 = 0.00 2.00 -4.00 ;
  79. *
  80. * ---------- MAILLAGE ----------
  81. *
  82. IH1 = IA1 DROIT 1 IB1 ;
  83. IH2 = IB1 DROIT 1 IC1 ;
  84. IH3 = IC1 DROIT 1 ID1 ;
  85. IH4 = ID1 DROIT 1 IA1 ;
  86. *
  87. IL1 = IH1 ET IH2 ET IH3 ET IH4 ;
  88. IS1 = SURF IL1 PLAN ;
  89. *
  90. OPTION ELEM CUB8 ;
  91. *
  92. VOLS = SS1 VOLU 1 TRAN ( 0.00 0.00 -2.00 ) ;
  93. FACS = FACE 2 VOLS ;
  94. VOLI = IS1 VOLU 1 TRAN ( 0.00 0.00 2.00 ) ;
  95. FACI = FACE 2 VOLI ;
  96. VT1 = VOLS ET VOLI ;
  97. VT1 = ELIM 0.00001 VT1 ;
  98. VTSAN = VT1 ;
  99. *
  100. * ------- DEFINITION DES CONDITIONS AUX LIMITES -------
  101. *
  102. CL11 = BLOQ IC1 UX ;
  103. CL12 = BLOQ IC1 UY ;
  104. CL13 = BLOQ IC1 UZ ;
  105. CL1 = CL11 ET CL12 ET CL13 ;
  106. *
  107. CL21 = BLOQ IH2 UX ;
  108. CL22 = BLOQ IH2 UZ ;
  109. CL2 = CL21 ET CL22 ;
  110. *
  111. CL31 = BLOQ IH3 UY ;
  112. CL32 = BLOQ IH3 UZ ;
  113. CL3 = CL31 ET CL32 ;
  114. *
  115. CL4 = BLOQ IA1 UZ ;
  116. *
  117. CL = CL1 ET CL2 ET CL3 ET CL4 ;
  118. *
  119. * ---------- DEFINITION DU MODELE ----------
  120. *
  121. MOD1 = MODE VT1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' ;
  122. MA1 = MATE MOD1 YOUN 2.1E11 NU 0.3 ;
  123. *
  124. * ---------- MATRICE DE RIGIDITE ----------
  125. *
  126. RI1 = RIGI MOD1 MA1 ;
  127. RI2 = RI1 ET CL ;
  128. *
  129. * ---------- PRESSION UNIFORME ----------
  130. *
  131. CHPO1= PRES 'MASS' MOD1 -100000.00 SS1 ;
  132. *
  133. * ---------- RESOLUTION ----------
  134. *
  135. RE = RESO RI2 CHPO1 ;
  136. *
  137. * ------REDUC. DES DEPLA. AUX DIFFERENTS ELEMENTS------
  138. *
  139. DEPS = REDU RE VOLS ;
  140. MESS ' DEPLACEMENTS DES NOEUDS DU CUBE SUPERIEUR ' ;
  141. LIST DEPS;
  142. *
  143. desusa=deps;
  144. *
  145. DEPI = REDU RE VOLI ;
  146. MESS ' DEPLACEMENTS DES NOEUDS DU CUBE INFERIEUR ' ;
  147. *
  148. deinsa=depi;
  149. *
  150. *
  151. * ---------- CALCUL DES DEFORMATIONS ----------
  152. *
  153. EPS1 = EPSI MOD1 RE ;
  154. *
  155. * ------REDUC. DES DEFOR. AUX DIFFERENTS ELEMENTS------
  156. *
  157. EPSS = REDU EPS1 VOLS ;
  158. MESS ' DEFORMATIONS DU CUBE DU HAUT ' ;
  159. LIST EPSS ;
  160. *
  161. EPSB = REDU EPS1 VOLI ;
  162. MESS ' DEFORMATIONS DU CUBE DU BAS ' ;
  163. LIST EPSB ;
  164. *
  165. * ------- CALCUL DES CONTRAINTES --------
  166. *
  167. SIG1 = SIGMA MOD1 MA1 RE ;
  168. *
  169. * -----REDUC. DES CONTR. AUX DIFFERENTS ELEMENTS-------
  170. *
  171. SIGS = REDU SIG1 VOLS ;
  172. MESS ' CONTRAINTES DU CUBE DU HAUT ' ;
  173. *
  174. SIGI = REDU SIG1 VOLI ;
  175. MESS ' CONTRAINTES DU CUBE DU BAS ' ;
  176. LIST SIGI ;
  177. *
  178. MESS ' ' ;
  179. MESS ' ' ;
  180. MESS' *=============================================*';
  181. *
  182. SAUT PAGE;
  183. MESS ' ' ;
  184. MESS ' ' ;
  185. MESS ' _____________________________________________' ;
  186. MESS '! !' ;
  187. MESS '! CAS DE DEUX CUBES SUPERPOSES AVEC UN !' ;
  188. MESS '! JOINT ENTRE EUX !' ;
  189. MESS '!____________________________________________!' ;
  190. MESS ' ' ;
  191. MESS ' ' ;
  192. *
  193. *
  194. OPTION ELEM QUA4 MODE TRID ;
  195. *
  196. * -----DEFINITION DE LA GEOMETRIE DU CUBE SUPERIEUR----
  197. *
  198. SA1A = 0.00 0.00 0.00 ;
  199. SB1A = 2.00 0.00 0.00 ;
  200. SC1A = 2.00 2.00 0.00 ;
  201. SD1A = 0.00 2.00 0.00 ;
  202. *
  203. * ---------- MAILLAGE ----------
  204. *
  205. SH1A = SA1A DROIT 1 SB1A ;
  206. SH2A = SB1A DROIT 1 SC1A ;
  207. SH3A = SC1A DROIT 1 SD1A ;
  208. SH4A = SD1A DROIT 1 SA1A ;
  209. *
  210. SL1A = SH1A ET SH2A ET SH3A ET SH4A ;
  211. SS1A = SURF SL1A PLAN ;
  212. *
  213. * ----DEFINITION DE LA GEOMETRIE DU CUBE INFERIEUR-----
  214. *
  215. IA1A = 0.00 0.00 -4.00 ;
  216. IB1A = 2.00 0.00 -4.00 ;
  217. IC1A = 2.00 2.00 -4.00 ;
  218. ID1A = 0.00 2.00 -4.00 ;
  219. *
  220. * ---------- MAILLAGE ----------
  221. *
  222. IH1A = IA1A DROIT 1 IB1A ;
  223. IH2A = IB1A DROIT 1 IC1A ;
  224. IH3A = IC1A DROIT 1 ID1A ;
  225. IH4A = ID1A DROIT 1 IA1A ;
  226. *
  227. IL1A = IH1A ET IH2A ET IH3A ET IH4A ;
  228. IS1A = SURF IL1A PLAN ;
  229. *
  230. OPTION ELEM CUB8 ;
  231. *
  232. VOLSA = SS1A VOLU 1 TRAN ( 0.00 0.00 -2.00 ) ;
  233. FACSA = FACE 2 VOLSA ;
  234. VOLIA = IS1A VOLU 1 TRAN ( 0.00 0.00 2.00 ) ;
  235. FACIA = FACE 2 VOLIA ;
  236. VT1A = VOLSA ET VOLIA ;
  237. *
  238. * ------ DEFINITION DES CONDITIONS AUX LIMITES --------
  239. *
  240. CL11A = BLOQ IC1A UX ;
  241. CL12A = BLOQ IC1A UY ;
  242. CL13A = BLOQ IC1A UZ ;
  243. CL1A = CL11A ET CL12A ET CL13A ;
  244. *
  245. CL21A = BLOQ IH2A UX ;
  246. CL22A = BLOQ IH2A UZ ;
  247. CL2A = CL21A ET CL22A ;
  248. *
  249. CL31A = BLOQ IH3A UY ;
  250. CL32A = BLOQ IH3A UZ ;
  251. CL3A = CL31A ET CL32A ;
  252. *
  253. CL4A = BLOQ IA1A UZ ;
  254. *
  255. CLA = CL1A ET CL2A ET CL3A ET CL4A ;
  256. *
  257. * ---- DEFINITION DU MODELE DES DEUX CUBES ----
  258. *
  259. MOD1A = MODE VT1A 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' ;
  260. MA1A = MATE MOD1A YOUN 2.1E11 NU 0.3 ;
  261. *
  262. * - CREATION DE L'ELEMENT JOINT ENTRE LES DEUX CUBES -
  263. *
  264. OPTION ELEM LIA4 ;
  265. VOLJA = LIAISON 0.0001 FACSA FACIA ;
  266. *
  267. VTOTAL=VT1A ET VOLJA ;
  268. *
  269. * ---------- DEFINITION DU MODELE DU JOINT ----------
  270. *
  271. MODJA = MODE VOLJA 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  272. JOI4 ;
  273. MAJA = MATE MODJA KS 4.2E15 KN 4.2E15 ;
  274. *
  275. * ----- MODELE ET MATERIAU DE TOUTE LA STRUCTURE ------
  276. *
  277. MODTA = MOD1A ET MODJA ;
  278. MATTA = MA1A ET MAJA ;
  279. *
  280. * ---------- MATRICE DE RIGIDITE ----------
  281. *
  282. RITA = RIGI MODTA MATTA ;
  283. RITOA = RITA ET CLA ;
  284. *
  285. * ---------- PRESSION UNIFORME ----------
  286. *
  287. CHPO1A= PRES 'MASS' MODTA -100000.00 SS1A ;
  288. *
  289. * ---------- RESOLUTION ----------
  290. *
  291. RETA = RESO RITOA CHPO1A ;
  292. *
  293. * ----- REDUC. DES DEPLA. AUX DIFFERENTS ELEMENTS -----
  294. *
  295. DEPSA = REDU RETA VOLSA ;
  296. MESS ' DEPLACEMENTS DES NOEUDS DU CUBE SUPERIEUR ' ;
  297. LIST DEPSA;
  298. *
  299. DEPIA = REDU RETA VOLIA ;
  300. MESS ' DEPLACEMENTS DES NOEUDS DU CUBE INFERIEUR ' ;
  301. LIST DEPIA ;
  302. *
  303. DEPJ = REDU RETA VOLJA ;
  304. MESS ' DEPLACEMENTS DES NOEUDS DU JOINT ' ;
  305. LIST DEPJ;
  306. * ---------- CALCUL DES DEFORMATIONS ----------
  307. *
  308. EPSTA = EPSI MODTA RETA ;
  309. *
  310. * ---- REDUC. DES DEFOR. AUX DIFFERENTS ELEMENTS ------
  311. *
  312. EPSSA = REDU EPSTA VOLSA ;
  313. MESS ' DEFORMATIONS DU CUBE DU HAUT ' ;
  314. LIST EPSSA ;
  315. *
  316. EPSBA = REDU EPSTA VOLIA ;
  317. MESS ' DEFORMATIONS DU CUBE DU BAS ' ;
  318. LIST EPSBA ;
  319. *
  320. EPSJA = REDU EPSTA VOLJA ;
  321. MESS ' DEFORMATIONS DU JOINT ' ;
  322. LIST EPSJA ;
  323. *
  324. * ---------- CALCUL DES CONTRAINTES ----------
  325. *
  326. SIGTA = SIGMA MODTA MATTA RETA ;
  327. *
  328. * ----- REDUC. DES CONTR. AUX DIFFERENTS ELEMENTS -----
  329. *
  330. SIGSA = REDU SIGTA VOLSA ;
  331. MESS ' CONTRAINTES DU CUBE DU HAUT ' ;
  332. LIST SIGSA ;
  333. *
  334. SIGIA = REDU SIGTA VOLIA ;
  335. MESS ' CONTRAINTES DU CUBE DU BAS ' ;
  336. LIST SIGIA ;
  337. *
  338. SIGJA = REDU SIGTA VOLJA ;
  339. MESS ' CONTRAINTES DU JOINT ' ;
  340. LIST SIGJA ;
  341.  
  342. *
  343. * ---------- CODE DE FONCTIONNEMENT ----------
  344. *
  345. *
  346. desuav=deps;
  347. * difference des deplacements configurations avec/sans
  348. * joint pour les cubes superieurs
  349. desdif=desusa-desuav;
  350. *
  351. deinav=depi;
  352. * difference des deplacements configurations avec/sans
  353. * joint pour les cubes inferieurs
  354. deidif=deinsa-deinav;
  355. *
  356. maxdis = maxi desdif 'ABS' ;
  357. maxdii = maxi deidif 'ABS' ;
  358. *
  359. *
  360. si ((maxdis <eg 1e-9) et (maxdii <eg 1e-9));
  361. erre 0 ;
  362. sino ;
  363. erre 5 ;
  364. finsi ;
  365. *
  366.  
  367. *
  368. FIN;
  369.  
  370.  
  371.  
  372.  
  373.  
  374.  
  375.  
  376.  
  377.  
  378.  
  379.  

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