Télécharger grandksi.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : grandksi.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. complet = faux;
  5. *
  6. *************************************************
  7. * pour calcul complet mettre complet à : vrai;
  8. ***************************************************
  9. opti echo 0;
  10. ********************************************************
  11. * *
  12. * calcul fleche d'une plaque sous son propre poids *
  13. * *
  14. * la plaque est horizontale : on joue sur inclinaison *
  15. * *
  16. * de G et de la Force de Traction *
  17. * *
  18. * TRAITE par METHODE des GRANDS DEPLACEMEMTS *
  19. * *
  20. * *
  21. * *
  22. * *
  23. ********************************************************
  24. *
  25. *
  26. * geometrie
  27. *
  28. *
  29. opti dime 3 elem qua4 ;
  30. *
  31. inclin = 0. ;
  32. *
  33. teta = 90. ;
  34. *
  35. temp = 30. ;
  36. mess ' effort traction en Newton 0 si pas d'' efforts sinon 15 ' ;
  37. eft = 15;
  38. *obte eft ;
  39. *
  40. epu = 0. ;
  41. ru=19520.;
  42. epc=7E-04 ;
  43. rc=1400.;
  44. *
  45. rcf=2.*epu*ru/(2*epu+epc)+rc;
  46. dens 1 ;
  47. vu2 = 1000. 0. 0. ;
  48. vu3 = 0 100. 100. ;
  49. *
  50. long = 0.92 ;
  51. larg = 0.04 ;
  52. *
  53. * coordonnees des 6 sommets initiaux
  54. *
  55. p1 = 0. 0. 0. ;
  56. p2 = 0. .24 0. ;
  57. p3 = 0. .68 0. ;
  58. p4 = 0. .92 0. ;
  59. *
  60. p5 = 0. 0. larg ;
  61. p6 = 0. .24 larg ;
  62. p7 = 0. .68 larg ;
  63. p8 = 0. .92 larg ;
  64. *
  65. *
  66. * creation du point de traction du bec
  67. *
  68. af = p4 plus p8 ; af = af /2 ;
  69. *
  70. *
  71. * creation lignes s'appuyant sur les sommets
  72. *
  73. *
  74. l1 = droi 2 p1 p5 ;
  75. l2 = droi 2 p2 p6 ;
  76. l3 = droi 2 p3 p7 ;
  77. l4 = droi 2 p4 p8 ;
  78. *
  79. * creation des milieux des lignes l2 et l3
  80. *
  81. m2 = 0. .24 ( larg/2 ) ;
  82. m3 = 0. .35 ( larg/2 ) ;
  83. m4 = 0. .46 ( larg/2 ) ;
  84. m5 = 0. .57 ( larg/2 ) ;
  85. m6 = 0. .68 ( larg/2 ) ;
  86. *
  87. *
  88. *
  89. si complet ;
  90. na = 9;nb = 12; nc = 6;
  91. sinon;
  92. aa = 7; nb = 8; nc = 4;
  93. finsi;
  94. c1b = droi 9 p1 p2 ;
  95. c2b = droi 12 p2 p3 ;
  96. c3b = droi 6 p3 p4 ;
  97. *
  98. c1h = droi 9 p5 p6 ;
  99. c2h = droi 12 p6 p7 ;
  100. c3h = droi 6 p7 p8 ;
  101. *
  102. * creation maillage
  103. *
  104. s1 = dall l1 c1b l2 c1h ;
  105. s2 = dall l2 c2b l3 c2h ;
  106. s3 = dall l3 c3b l4 c3h ;
  107. *
  108. sbec = s1 et s2 et s3 ;
  109. *trac sbec;
  110. sbec = sbec coul jaun ;
  111. pf = sbec poin proc af ;
  112. *
  113. *********************************************************************
  114. * *
  115. * impression des donnees *
  116. * *
  117. *********************************************************************
  118. *
  119. mess ' angle inclinaison inclin ' inclin 'degres' ;
  120. mess ' angle rotation teta ' teta 'degres' ;
  121. mess ' angle = 90 plaque a plat ' ;
  122. mess ' temperature temp ' temp ;
  123. mess ' force de traction en N eft ' eft ;
  124. mess ' epaisseur U epu ' epu ;
  125. mess ' longueur plaque long ' long ;
  126. mess ' largeur plaque larg ' larg ;
  127. mess ' epaisseur bec epc ' epc ;
  128. mess ' ro equivalent rcf ' rcf ;
  129. *
  130. *
  131. *********************************************************************
  132. *********************************************************************
  133. *
  134. ********** modele du materiau et caracteristiques
  135. *
  136. *
  137. mo = mode sbec MECANIQUE ELASTIQUE coq4 ;
  138. matx = mate mo 'YOUN' 35E+09 'RHO' rcf 'ALPHA' 2.5E-06 'NU' 0.3 ;
  139. ca = cara mo EPAI epc ;
  140. mat1 = matx et ca ;
  141. *
  142. ********** matrice de masse ; acceleration pesanteur ; poids
  143. *
  144. m = mass mo mat1 ;
  145. *
  146. stet = sin teta ; ctet = cos teta ;
  147. cic = cos inclin ; sic = sin inclin ;
  148. gx = +9.81 * stet * cic ;
  149. gy = +9.81 * sic ;
  150. gz = -9.81 * ctet * cic ;
  151. accel = manu chpo sbec 6 ux gx uy gy uz gz mx 0. my 0. mz 0. ;
  152. *
  153. poids = m * accel ;
  154. *
  155. ********** condition aux limites
  156. *
  157. *
  158. * encastrement = bloq rota et depl
  159. *
  160. *
  161. * appui = bloq depl ; variante = bloq ux et-ou uy et-ou uz
  162. *
  163. ********** blocage rotation et deplacement l1
  164. *
  165. bl1rd = bloq rota depl l1 ;
  166. *
  167. ********** blocage deplacement [ ux uz ] l4
  168. *
  169. bl4d = bloq l4 ux uz ;
  170. *
  171. blqtot = bl1rd et bl4d ;
  172. *
  173. ********** creation force de traction
  174. *
  175. fo1 = force (0. eft 0.) pf ;
  176. *
  177. ********** RESOLUTION PAR RESO
  178. *
  179. sc = poids + fo1 ;
  180. ******************************** RESOLUTION PAR NONLIN
  181. *
  182. *
  183. ********** creation chargement
  184. *
  185. scf = fo1 ; scp = poids ;
  186. ev1 = evol manu xx1 ( prog 0. 1.1 ) ff1 ( prog 0. 1.1 ) ;
  187. ch1 = char meca scf ev1 ;
  188. ch2 = char meca scp ev1 ;
  189. cht = ch1 et ch2 ;
  190. sctot = scf + scp;
  191. *
  192. **********
  193. *
  194. tgd = table ;
  195. tgd.grands_deplacements = vrai ;
  196. tgd.k_sigma = vrai ;
  197. tgd.'REAC_GRANDS'=100.;
  198. * liste des temps a calculer
  199. si complet ;
  200. tfin = 0.37;
  201. sinon;
  202. tfin = 0.07;
  203. finsi;
  204. tcal = prog 0. 0.0001 0.0003 0.001 0.003 0.01 0.02
  205. 0.04 pas 0.03 tfin;
  206.  
  207. tgd.blocages_mecaniques = blqtot;
  208. tgd.caracteristiques = mat1;
  209. tgd.modele = mo;
  210. tgd.temps_calcules = tcal;
  211. tgd.chargement = cht;
  212. tgd.hypothese_deformations = 'LINEAIRE' ;
  213.  
  214.  
  215. pasapas tgd ;
  216.  
  217. ********** resultats
  218. *
  219. dep1 = peche tgd deplacements tfin;
  220. aa = maxi ( exco dep1 UX);
  221. si complet;
  222. valcomp = 2.0938e-03;
  223. sinon;
  224. valcomp = 1.24824E-03;
  225. finsi;
  226. err = ( aa - valcomp ) / valcomp abs;
  227. err100 = err * 100;
  228. mess ' ';
  229. mess ' ';
  230. mess ' ';
  231. mess ' ';
  232. mess ' la fleche maximum attendue est ' valcomp;
  233. mess ' ';
  234. mess ' ';
  235. mess ' valeur de la fleche trouvée : ' aa ;
  236. mess ' ';
  237. mess ' ';
  238. mess ' erreur commise en % ' err100;
  239. mess ' ';
  240. mess ' ';
  241. mess ' ';
  242. mess ' ';
  243. si ( err > 0.005) ; erreur 5; finsi;
  244. fin;
  245.  
  246.  
  247.  
  248.  
  249.  
  250.  
  251.  
  252.  
  253.  
  254.  
  255.  
  256.  
  257.  
  258.  
  259.  
  260.  

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