Télécharger joi44.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : joi44.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Joi44.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7. OPTI ECHO 0 ;
  8. SAUT PAGE ;
  9. TEMPS ;
  10. *=====================================================*
  11. * *
  12. * TEST JOI44 *
  13. * *
  14. * CALCUL DES CONTRAINTES *
  15. * THERMIQUES SUR UN JOINT 3D *
  16. * *
  17. * Un joint 3D JOI4 a sa surface inferieure *
  18. * encastree. Sa surface superieure est libre. *
  19. * Un champ de temperatures est donne, de type *
  20. * CHAMELEM, de sous type temperatures avec le *
  21. * MCHAML s'appuyant d'une part sur un objet de *
  22. * type MAILLAGE, d'autre part sur un objet de *
  23. * type MMODEL. *
  24. * On cherche a calculer le champ des contraintes *
  25. * thermiques induit par ce champ de temperatures.*
  26. * *
  27. * Solution analytique : *
  28. * *
  29. * SIGMA = K * ALPN * DELTA(T) *
  30. * ou *
  31. * K = raideur dans la direction perpendiculaire *
  32. * au plan du joint *
  33. * ALPN = coefficient de dilatation thermique *
  34. * suivant la direction n, perpendiculaire *
  35. * au plan du joint *
  36. * ( ALPN est homogene a une longueur *
  37. * divisee par des degres, la longueur *
  38. * representant l'epaisseur du joint) *
  39. * DELTA(T) = CHAMALEM de temperatures *
  40. * SIGMA = 4.2E20 * 2.5 * 1000.0 *
  41. * = 1.05E24 *
  42. * *
  43. *=====================================================*
  44. OPTION ELEM QUA4 MODE TRID ;
  45. *
  46. SOLNUL = 0.0 ;
  47. SOLANA = 1.05E24 ;
  48. *
  49. * ------- DEFINITION DE LA SURFACE TOP DU JOINT -------
  50. *
  51. A1 = 0.00 0.00 0.00 ;
  52. B1 = 2.00 0.00 0.00 ;
  53. C1 = 2.00 2.00 0.00 ;
  54. D1 = 0.00 2.00 0.00 ;
  55. *
  56. * ---------- MAILLAGE ----------
  57. *
  58. O1 = 0.00 0.00 10.00 ;
  59. H1 = A1 DROIT 1 B1 ;
  60. H2 = B1 DROIT 1 C1 ;
  61. H3 = C1 DROIT 1 D1 ;
  62. H4 = D1 DROIT 1 A1 ;
  63. *
  64. L1 = H1 ET H2 ET H3 ET H4 ;
  65. S1 = SURF L1 PLAN ;
  66. *
  67. * ------- DEFINITION DE LA SURFACE BOT DU JOINT -------
  68. *
  69. IA1 = 0.00 0.00 0.00 ;
  70. IB1 = 2.00 0.00 0.00 ;
  71. IC1 = 2.00 2.00 0.00 ;
  72. ID1 = 0.00 2.00 0.00 ;
  73. *
  74. * ---------- MAILLAGE ----------
  75. *
  76. IH1 = IA1 DROIT 1 IB1 ;
  77. IH2 = IB1 DROIT 1 IC1 ;
  78. IH3 = IC1 DROIT 1 ID1 ;
  79. IH4 = ID1 DROIT 1 IA1 ;
  80. *
  81. IL1 = IH1 ET IH2 ET IH3 ET IH4 ;
  82. IS1 = SURF IL1 PLAN ;
  83. *
  84. * ---------- CREATION DU JOINT JOI4 ----------
  85. *
  86. OPTION ELEM LIA4 ;
  87. VOL = LIAISON 0.00001 S1 IS1 ;
  88. *
  89. * -------- DEFINITION DES CONDITIONS LIMITES ----------
  90. *
  91. CL11 = BLOQ IA1 UX ;
  92. CL12 = BLOQ IA1 UY ;
  93. CL13 = BLOQ IA1 UZ ;
  94. CL1 = CL11 ET CL12 ET CL13 ;
  95. *
  96. CL21 = BLOQ IB1 UX ;
  97. CL22 = BLOQ IB1 UY ;
  98. CL23 = BLOQ IB1 UZ ;
  99. CL2 = CL21 ET CL22 ET CL23 ;
  100. *
  101. CL31 = BLOQ IC1 UX ;
  102. CL32 = BLOQ IC1 UY ;
  103. CL33 = BLOQ IC1 UZ ;
  104. CL3 = CL31 ET CL32 ET CL33 ;
  105. *
  106. CL41 = BLOQ ID1 UX ;
  107. CL42 = BLOQ ID1 UY ;
  108. CL43 = BLOQ ID1 UZ ;
  109. CL4 = CL41 ET CL42 ET CL43 ;
  110. *
  111. CL = CL1 ET CL2 ET CL3 ET CL4 ;
  112. *
  113. * --------- DEFINITION DU MODELE DU JOINT ----------
  114. *
  115. P1 = 0.707 0.707 0.707 ;
  116. MOD1 = MODE VOL 'MECANIQUE' 'ELASTIQUE' 'ORTHOTROPE'
  117. JOI4 ;
  118. MA1 = MATE MOD1 'DIRECTION' P1 'PERPENDICULAIRE'
  119. KS1 4.2E10 KS2 3.0E5 KN 4.2E20 ALPN 2.5;
  120. *
  121. * ---------- MATRICE DE RIGIDITE ----------
  122. *
  123. RI1 = RIGI MOD1 MA1 ;
  124. RI2 = RI1 ET CL ;
  125. *
  126. * ---------- CONTRAINTES THERMIQUES ----------
  127. *
  128. * CHAMELEM S'APPUYANT SUR UN OBJET MODELE
  129. *
  130. C1T =MANU 'CHML' MOD1 'T' 1000.0 'TYPE' 'TEMPERATURES';
  131. SI1 = THET MOD1 MA1 C1T ;
  132. *
  133. MESS ' CHAMELEM S APPUYANT SUR UN OBJET MODELE ' ;
  134. MESS ' ' ;
  135. MESS ' ____________________________________________ ' ;
  136. MESS '! !' ;
  137. MESS '! JOINT 3D JOI4 !' ;
  138. MESS '! !' ;
  139. MESS '! - Surface inferieure du joint encastree. !' ;
  140. MESS '! - Surface superieure du joint libre . !' ;
  141. MESS '! - Chamelem de temperatures donne DELTA(T) !' ;
  142. MESS '! s appuyant sur un objet de type modele !' ;
  143. MESS '! !' ;
  144. MESS '! ==> !' ;
  145. MESS '! Champ des contraintes thermiques = Sigma !' ;
  146. MESS '! !' ;
  147. MESS '! Sigma = raideur * coeff. dilat. thermi. !' ;
  148. MESS '! * DELTA(T) !' ;
  149. MESS '!____________________________________________!' ;
  150. MESS ' ' ;
  151. MESS ' ' ;
  152. MESS ' Solution Analytique :' ;
  153. MESS ' ' ;
  154. MESS ' SMS1 =' SOLNUL ;
  155. MESS ' SMS2 =' SOLNUL ;
  156. MESS ' SMN =' SOLANA ;
  157. MESS ' ' ;
  158. MESS ' ' ;
  159. MESS ' ' ;
  160. MESS ' Solution Calculee :' ;
  161. LIST SI1 ;
  162. *
  163. * ---------- CODE DE FONCTIONNEMENT ----------
  164. *
  165. MAXCO1 = MAXI SI1 'ABS' ;
  166. RESI1 = ABS( (MAXCO1-SOLANA)/SOLANA ) ;
  167. *
  168. * CHAMELEM S'APPUYANT SUR UN OBJET MAILLAGE
  169. *
  170. C2T = MANU 'CHML' VOL 'T' 1000.0 'TYPE' 'TEMPERATURES';
  171. SI2 = THET MOD1 MA1 C2T ;
  172. MESS ' ' ;
  173. MESS ' CHAMELEM S APPUYANT SUR UN OBJET MAILLAGE ' ;
  174. MESS ' ' ;
  175. MESS ' ____________________________________________ ' ;
  176. MESS '! !' ;
  177. MESS '! JOINT 3D JOI4 !' ;
  178. MESS '! !' ;
  179. MESS '! - Surface inferieure du joint encastree. !' ;
  180. MESS '! - Surface superieure du joint libre . !' ;
  181. MESS '! - Chamelem de temperatures donne DELTA(T) !' ;
  182. MESS '! s appuyant sur un objet de type maillage !' ;
  183. MESS '! !' ;
  184. MESS '! ==> !' ;
  185. MESS '! Champ des contraintes thermiques = Sigma !' ;
  186. MESS '! !' ;
  187. MESS '! Sigma = raideur * coeff. dilat. thermi. !' ;
  188. MESS '! * DELTA(T) !' ;
  189. MESS '!____________________________________________!' ;
  190. MESS ' ' ;
  191. MESS ' ' ;
  192. MESS ' Solution Analytique :' ;
  193. MESS ' ' ;
  194. MESS ' SMS1 =' SOLNUL ;
  195. MESS ' SMS2 =' SOLNUL ;
  196. MESS ' SMN =' SOLANA ;
  197. MESS ' ' ;
  198. MESS ' ' ;
  199. MESS ' ' ;
  200. MESS ' Solution Calculee :' ;
  201. MESS ' ' ;
  202. LIST SI2 ;
  203. *
  204. MAXCO2 = MAXI SI2 'ABS' ;
  205. RESI2 = ABS( (MAXCO2-SOLANA)/SOLANA ) ;
  206. *
  207. SI ( (RESI1 <EG 1E-4) ET (RESI2 <EG 1E-4) );
  208. ERRE 0 ;
  209. SINO;
  210. ERRE 5 ;
  211. FINSI ;
  212. *
  213. FIN ;
  214.  
  215.  
  216.  
  217.  
  218.  
  219.  
  220.  
  221.  
  222.  
  223.  

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