Télécharger testICQ4.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : testICQ4.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test TestICQ4.dgibi: Jeux de données *
  5. * ------------------------------------ *
  6. * *
  7. *
  8. * FICHIER GIBIANE POUR TESTER LES ELEMENTS
  9. * INCOMPRESSIBLES ICQ4
  10. *
  11. * Maillage d'une plaque de forme carrée avec les
  12. * éléments ICQ4.
  13. *
  14. * Cas déformation plane.
  15. *
  16. * On soumet la plaque à un champ de contrainte
  17. * hydrostatique (contribution: alpha=1) et à un champ
  18. * de contrainte déviatorique (contribution: beta=0.001)
  19. *
  20. * Les résultats théoriques sont déterminés dans le cas
  21. * incompressible (coefficient de Poisson proche de 0.5)
  22. *
  23. * cf rapport DMT sur ces elements
  24. *
  25. *
  26. *-----MAILLAGE-----------
  27. *
  28. opti dime 2 mode plan defo elem qua4 ;
  29. *
  30. p1 = 0. 0. ;
  31. p2 = 100. 0. ;
  32. p3 = 100. 100. ;
  33. p4 = 0. 100. ;
  34. *
  35. li1 = d 10 p1 p2 ;
  36. li2 = d 10 p2 p3 ;
  37. li3 = d 10 p3 p4 ;
  38. li4 = d 10 p4 p1 ;
  39. *
  40. sur1 = dall li1 li2 li3 li4 'PLAN' ;
  41. *
  42. *
  43. *-----MODELE MECANIQUE------------
  44. *
  45. * Valeur du coefficient de Poisson
  46. *
  47. nu1 = 0.499 ;
  48. *
  49. mod1 = modele sur1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  50. icq4 ;
  51. mat1 = mate mod1 'YOUNG' 2.e11 'NU' nu1 ;
  52. rigi1 = rigi mod1 mat1 ;
  53. *
  54. *
  55. *-----CONDITIONS AUX LIMITES------------
  56. *
  57. bl1 = bloq li1 'UY' ;
  58. bl2 = bloq li4 'UY' ;
  59. bl3 = bloq p1 'UX' ;
  60. *
  61. *
  62. *----------------CHARGEMENT-----------------
  63. *
  64. * Coefficients de contribution de chargement
  65. * (alpha-->hydrostatique)
  66. * (beta-->déviatorique)
  67. *
  68. alpha = 1. ;
  69. beta = 1.e-3 ;
  70. *
  71. *-----Chargement déviatorique------------
  72. *
  73. *-----sur li1-----
  74. *
  75. chpx1 = coor 1 li1 ;
  76. chpresy1 = pres 'MASS' mod1 (chpx1 * (-2.)*1.e5) ;
  77. chpres1 = beta * chpresy1 ;
  78. vecfor1 = vect chpres1 1.E-7 'FX' 'FY' rouge ;
  79. *
  80. *-----sur li2-----
  81. *
  82. chpy2 = coor 2 li2 ;
  83. chpresx2 = pres 'MASS' mod1 (200.*1.e5) li2 ;
  84. chpresy2 = pres 'MASS' mod1 (chpy2 * (-2.)*1.e5) ;
  85. chforx2 = chpresx2 ;
  86. chfory2 = exco 'FX' chpresy2 'FY';
  87. chpres2 = beta * ( chforx2 et chfory2 ) ;
  88. vecfor2 = vect chpres2 'FX' 'FY' 1.E-7 rouge ;
  89. *
  90. *-----sur li3-----
  91. *
  92. chpx3 = coor 1 li3 ;
  93. chpresx3 = pres 'MASS' mod1 (-200.*1.e5) li3 ;
  94. chpresy3 = pres 'MASS' mod1 (chpx3 * (-2.)*1.e5) ;
  95. chforx3 = exco 'FY' chpresx3 'FX';
  96. chfory3 = chpresy3 ;
  97. chpres3 = beta * ( chforx3 et chfory3 ) ;
  98. vecfor3 = vect chpres3 'FX' 'FY' 1.E-7 rouge ;
  99. *
  100. *-----sur li4-------
  101. *
  102. chpy4 = coor 2 li4 ;
  103. chpresy4 = pres 'MASS' mod1 (chpy4 * (-2.)*1.e5) ;
  104. chfory4 = exco 'FX' chpresy4 'FY';
  105. chpres4 = beta * chfory4 ;
  106. vecfor4 = vect chpres4 1.E-7 'FX' 'FY' rouge ;
  107. *
  108. *
  109. *-----Chargement hydrostatique------------
  110. *
  111. *-----sur li1-----
  112. *
  113. presy1 = pres 'MASS' mod1 (1.e5) li1 ;
  114. pres1 = alpha * presy1 ;
  115. vfor1 = vect pres1 1.E-5 'FX' 'FY' vert ;
  116. *
  117. *-----sur li2-----
  118. *
  119. presx2 = pres 'MASS' mod1 (1.e5) li2 ;
  120. pres2 = alpha * presx2 ;
  121. vfor2 = vect pres2 'FX' 'FY' 1.E-5 vert ;
  122. *
  123. *-----sur li3-----
  124. *
  125. presy3 = pres 'MASS' mod1 (1.e5) li3 ;
  126. pres3 = alpha * presy3 ;
  127. vfor3 = vect pres3 'FX' 'FY' 1.E-5 vert ;
  128. *
  129. *-----sur li4-------
  130. *
  131. presx4 = pres 'MASS' mod1 (1.e5) li4 ;
  132. pres4 = alpha * presx4 ;
  133. vfor4 = vect pres4 1.E-5 'FX' 'FY' vert ;
  134. *
  135. *
  136. *---------CALCUL DES DEPLACEMENTS --------------------
  137. *
  138. *
  139. * Solution analytique au point P3: uyan.
  140. *
  141. x3 = coor 1 p3 ;
  142. y3 = coor 2 p3 ;
  143. mu = 1.e11 / ( 1 + nu1 ) ;
  144. uyan = ( 2 * x3 * y3 ) / ( 2 * mu ) * 1.e5 * beta ;
  145. *
  146. * Solution Castem au point P3: uycas.
  147. *
  148. deptot = reso (rigi1 et bl1 et bl2 et bl3)
  149. (chpres2 et chpres3 et chpres4 et chpres1
  150. et pres2 et pres3 et pres4 et pres1) ;
  151. uycas = extr deptot 'UY' p3 ;
  152. *
  153. * Calcul de l'erreur sur le déplacement.
  154. *
  155. erry = abs (( uyan - uycas ) / ( uyan )) * 100 ;
  156. *
  157. * Test de précision sur le calcul du déplacement uy
  158. * (<1%).
  159. *
  160. si ( erry >eg 1. ) ;
  161. erreur 5 ;
  162. finsi ;
  163. *
  164. *
  165. *-----------CALCUL CONTRAINTES-------------------------
  166. *
  167. *
  168. * Solution analytique: pression = alpha*1.E5.
  169. *
  170. * Solution Castem : pression maximale = maxpres1.
  171. *
  172. sig1 = sigma mod1 deptot mat1 ;
  173. sxx = exco 'SMXX' sig1 'P' ;
  174. syy = exco 'SMYY' sig1 'P' ;
  175. szz = exco 'SMZZ' sig1 'P' ;
  176. sp = (sxx + syy + szz) / 3. ;
  177. maxpres1 = maxi sp ;
  178. *
  179. * Calcul de l'erreur sur la pression.
  180. *
  181. errpress = abs ((maxpres1 + (alpha * 1.E5)) /
  182. (alpha * 1.E5)) * 100. ;
  183. *
  184. * Test de précision sur le calcul de la pression (<1%).
  185. *
  186. si ( errpress >eg 1. ) ;
  187. erreur 5 ;
  188. finsi ;
  189. *
  190. *
  191. fin;
  192.  
  193.  
  194.  
  195.  
  196.  
  197.  
  198.  

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