Télécharger la plaquette Cast3M


fr en

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

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