Télécharger tc3bired.dgibi

Retour à la liste

Numérotation des lignes :

  1. * test snap back dans autopilot
  2. * on verifie principalement que le pilotage fonctionne correctement
  3. * on peut accessoirement verifier que la liste des temps calcules n'est pas monotone
  4.  
  5.  
  6.  
  7. opti echo 1 trac ps;;
  8.  
  9. opti echo 1 ;
  10.  
  11.  
  12.  
  13. opti dime 2 elem QUA8 mode plan defo ;
  14. h1 = 1. ; h2 = 1. ; h3 = 8. ;
  15. e1 = 5e-2 ; e2 = 1.20 ; e3 = 5e-2 ;
  16.  
  17. zz = 0. ;
  18. p1 = 0 zz ; p2 = e1 zz ; p3 = (e1 + e2) zz ; p4 = (e1 + e2 + e3) zz ;
  19. x_moy = (e1 + e2 + e3) / 2. ;
  20. zz = h1 ;
  21. p5 = 0 zz ; p6 = e1 zz ; p7 = (e1 + e2) zz ; p8 = (e1 + e2 + e3) zz ;
  22. zz = h1 + h2 ;
  23. p9 = 0 zz ; p10 = e1 zz ; p11 = (e1 + e2) zz ; p12 = (e1 + e2 + e3) zz;
  24. zz = h1 + h2 + h3 ;
  25. p13 = 0 zz ; p14 = e1 zz ; p15 = (e1 + e2) zz; p16 = (e1 + e2 + e3) zz;
  26.  
  27. ds1 = 75e-3 ; ds2 = 15.e-2 ; ds3 = 225e-3 ;
  28. ds1 = 600e-3 ; ds2 = 120.e-2 ; ds3 = 1350e-3 ;
  29. ds1=ds1/2; ds2=ds2/2; ds3=ds3/2;
  30. ds1=ds1/2; ds2=ds2/2; ds3=ds3/4;
  31. l12 = p1 d p2 dini ds1 dfin ds1 ;
  32. l26 = p2 d p6 dini ds3 dfin ds3 ;
  33. l65 = p6 d p5 dini ds1 dfin ds1 ;
  34. l51 = p5 d p1 dini ds3 dfin ds3 ;
  35. s12 = dall l12 l26 l65 l51 ;
  36. l23 = p2 d p3 dini ds2 dfin ds2 ;
  37. l37 = p3 d p7 dini ds3 dfin ds3 ;
  38. l76 = p7 d p6 dini ds2 dfin ds2 ;
  39. l62 = inve l26 ;
  40. s23 = dall l23 l37 l76 l62 ;
  41. l34 = p3 d p4 dini ds1 dfin ds1 ;
  42. l48 = p4 d p8 dini ds3 dfin ds3 ;
  43. l87 = p8 d p7 dini ds1 dfin ds1 ;
  44. l73 = inve l37 ;
  45. s34 = dall l34 l48 l87 l73 ;
  46. l56 = inve l65 ;
  47. l610 = p6 d p10 dini ds3 dfin ds3 ;
  48. l109 = p10 d p9 dini ds1 dfin ds1 ;
  49. l95 = p9 d p5 dini ds3 dfin ds3 ;
  50. s56 = dall l56 l610 l109 l95 ;
  51. l67 = inve l76 ;
  52. l711 = p7 d p11 dini ds3 dfin ds3 ;
  53. l1110 = p11 d p10 dini ds2 dfin ds2 ;
  54. l106 = inve l610 ;
  55. s67 = dall l67 l711 l1110 l106 ;
  56. l78 = inve l87 ;
  57. l812 = p8 d p12 dini ds3 dfin ds3 ;
  58. l1211 = p12 d p11 dini ds1 dfin ds1 ;
  59. l117 = inve l711 ;
  60. s78 = dall l78 l812 l1211 l117 ;
  61. l910 = inve l109 ;
  62. l1014 = p10 d p14 dini ds3 dfin ds3 ;
  63. l1413 = p14 d p13 dini ds1 dfin ds1 ;
  64. l139 = p13 d p9 dini ds3 dfin ds3 ;
  65. s910 = dall l910 l1014 l1413 l139 ;
  66. l1011 = inve l1110;
  67. l1115 = p11 d p15 dini ds3 dfin ds3 ;
  68. l1514 = p15 d p14 dini ds2 dfin ds2 ;
  69. l1410 = inve l1014;
  70. s1011 = dall l1011 l1115 l1514 l1410 ;
  71. l1112 = inve l1211;
  72. l1216 = p12 d p16 dini ds3 dfin ds3 ;
  73. l1615 = p16 d p15 dini ds1 dfin ds1 ;
  74. l1511 = inve l1115;
  75. s1112 = dall l1112 l1216 l1615 l1511 ;
  76.  
  77. couch1 = (s12 et s23 et s34) coul bleu ;
  78. couch2 = (s56 et s67 et s78) coul blan ;
  79. couch3 = (s910 et s1011 et s1112) coul roug ;
  80. su = couch1 et couch2 et couch3 ;
  81. xsubto = mesu su surf ;
  82. mess 'xsubto=' xsubto 'm2' ;
  83. titr 'Maillage' (nbno su) 'noeuds' ;
  84. trac su ;
  85. *trac noeud su ; trac elem su ; trac (cont su) ;
  86.  
  87. lgau = (l139 et l95 et l51) coul vert ;
  88. lbas = (l12 et l23 et l34) coul jaun ;
  89. ldro = (l48 et l812 et l1216) coul turq ;
  90. lhau = (l1615 et l1514 et l1413) coul bleu ;
  91. titr 'Conditions aux limites' ;
  92. *trac (lgau et lbas et ldro et lhau) ;
  93.  
  94. fer1 = (l26 et l610 et l1014) coul blan ;
  95. fer2 = (l37 et l711 et l1115) coul rose ;
  96. fers = fer1 et fer2 ;
  97. titr 'Ferraillage' ;
  98. *trac fers ;
  99. trac ((cont su) et fers) ;
  100.  
  101. voile = su et fers ; tass voile ; mena ;
  102. titr 'Maillage complet' (nbel voile) 'elements' ;
  103. *trac ecla voile ;
  104.  
  105. mo123 = mode su mecanique elastique isotrope plastique
  106. drucker_prager ; "
  107. mof1 = mode(changer (l26 et l610 et l37 et l711) seg2)
  108. mecanique elastique isotrope plastique cinematique coq2 ;
  109. mof2 = mode ( changer (l1014 et l1115) seg2)
  110. mecanique elastique isotrope plastique cinematique coq2 ;
  111. mof = mof1 et mof2 ;
  112. mo123f = mo123 et mof ;
  113.  
  114.  
  115. maf1 = mate mof1 youn 200000e6 nu 0.3 rho 7800. sigy 400.e6 h 100000e6 ;
  116. xep1 = 5. * pi * ((39.e-3 / 2.) ** 2) ;
  117. mess 'Epaisseur nappe 1' xep1 'm2' ;
  118. caf1 = cara mof1 epai xep1 ;
  119. rif1 = rigi mof1 (maf1 et caf1);
  120. mmf1 = mass mof1 (maf1 et caf1);
  121. maf2 = mate mof2 youn 200000e6 nu 0.3 rho 7800. sigy 400.e6 h 100000e6 ;
  122. xep2 = 5. * pi * ((25.e-3 / 2.) ** 2) ;
  123. mess 'Epaisseur nappe 2' xep2 'm2' ;
  124. caf2 = cara mof2 epai xep2 ;
  125. rif2 = rigi mof2 (maf2 et caf2);
  126. mmf2 = mass mof2 (maf2 et caf2);
  127.  
  128. maf = maf1 et maf2 ;
  129. caf = caf1 et caf2 ;
  130. rif = rif1 et rif2 ;
  131. mmf = mmf1 et mmf2 ;
  132.  
  133. cl = (bloq depl lbas) et (bloq rota lbas) ;
  134.  
  135. chprof = (manu chpo lgau 1 scal (h1 + h2 + h3)) - (coor 2 lgau) ;
  136. chpres = chprof * (1000. * 9.81) ;
  137. fopres = pres mass mo123 chpres ;
  138. vepres = vect fopres 1e-4 fx fy bleu ;
  139. xr = extr (resu fopres) fx ((extr (resu fopres) mail) poin 1) ;
  140. xrt= (h1 + h2 + h3) * 1000. * 9.81 * ((h1 + h2 + h3) / 2.) ;
  141. mess 'press ef=' xr 'press th=' xrt 'ec=' ((xr - xrt) / xrt) ;
  142. titr 'Pression, resultante fx=' xr 'N' ;
  143. trac vepres (cont su) ;
  144.  
  145. xyo = 35000e6 ; xcp = 50e6 ; xtr = 0.6e6 + (0.06 * xcp) ;
  146. xft = xtr ; xfc = -1 * xcp ;
  147. xal = (xfc + xft) / (xfc - xft) ;
  148. xbe = 1. ;
  149. xk = (2. * xfc * xft) / (xfc - xft) ;
  150. xga = xal ;
  151. xde = xbe ;
  152. xet = xal ;
  153. xmu = xbe ;
  154. xkl = xk ;
  155. xyo2 = xyo / -2. ;
  156. xh = ((xal + 1.) ** 2) * ((xyo * xyo2) / (xyo - xyo2)) ;
  157. opti echo 0 ;
  158. mess ' ' ;
  159. mess 'Beton : trac=' xtr 'comp=' xcp 'Pa' ;
  160. mess ' Youn=' xyo 'Pa' ;
  161. mess '------------------------------------------------------';
  162. mess 'Param : xft=' xft 'xfc=' xfc 'Pa' ;
  163. mess '------------------------------------------------------';
  164. mess 'Criteres :' ;
  165. mess ' xal=' xal 'xbe=' xbe 'xk=' xk 'Pa' ;
  166. mess ' xga=' xga 'xde=' xde ;
  167. mess ' xet=' xet 'xmu=' xmu 'xkl=' xkl 'Pa' ;
  168. mess ' xh=' xh 'Pa' ;
  169. opti echo 1 ;
  170.  
  171. ma123 = mate mo123 youn xyo nu 0.2 rho 2500.
  172. alfa xal beta xbe k xk gamm xga delt xde
  173. eta xet mu xmu kl xkl h xh ;
  174.  
  175. ri123 = rigi mo123 ma123 ;
  176. mm123 = mass mo123 ma123 ;
  177.  
  178. mm = mm123 et mmf ;
  179. gg = manu chpo voile 1 uy -9.81 ;
  180. fopoid = mm * gg ;
  181. vepoid = vect fopoid 2.e-4 fx fy jaun ;
  182. xr = extr (resu fopoid) fy ((extr (resu fopoid) mail) poin 1) ;
  183. titr 'Poids , resultante fy=' xr 'N' ;
  184. trac vepoid (cont su) ;
  185. titr 'Evolutions du chargement en temps' ;
  186. ev011 = evol bleu manu 'temps' (prog 0 1 2 10) 'F' (prog 0 1 1 1) ;
  187. chr1 = char meca fopoid ev011 ;
  188. ev001 = evol roug manu 'temps' (prog 0 1 2 4 ) 'F' (prog 0 0 1 2) ;
  189. dess (ev011 et ev001) mima ;
  190. chr2 = char meca fopres ev001 ;
  191. chr = chr1 et chr2 ;
  192.  
  193. ma123f = ma123 et maf et caf ;
  194. ri = (ri123 et rif) ;
  195.  
  196. ma = mm123 ;
  197. cp = manu chpo (ma extr 'MAIL' 'NOMU') 1 ux 1;
  198. po = ma * cp;
  199. po2 = po chan cham mo123 gravite; po = po2 chan chpo mo123;
  200.  
  201. cpv = (po exco fx) ** 0.5;
  202. **optio trac x isov surf;
  203. trac su cpv;
  204. *list cpv;
  205.  
  206. tab1 = table ;
  207. tab1.plastique = vrai ;
  208. tab1.precision = 1e-4;
  209. tab1.precdecharge = 1e-4;
  210. tab1.autoresu = 1;
  211. tab1.grands_deplacements = vrai ;
  212. tab1.'K_SIGMA' = vrai ;
  213. tab1.masse = ma;
  214. tab1.modele = mo123f;
  215. tab1.caracteristiques=ma123f;
  216. tab1.chargement=chr;
  217. tab1.BLOCAGES_MECANIQUES=cl;
  218. ltps = prog 0 1.5 1.6 1.65 pas 0.3 1.80 ;
  219. *ltps = prog 0 1.5 1.6 1.65 pas 0.03 1.85;
  220. *ltps = prog 0 1.5 1.6 pas 0.01 1.65 pas 0.006 1.85;
  221. *ltps = prog 0 1.5 1.6 1.65 pas 0.1 1.85;
  222. tab1.temps_calcules=ltps;
  223. tab1.maxdefor=1e-3;
  224. tab1.'PRECISION'=1E-5;
  225. tab1.deltaiter=29;
  226. ripv = ri et cl;
  227. tab1.automatique = vrai ;
  228. tab1.autopas = 5000 ;
  229. tab1.'MAXSOUSPAS'=5000;
  230. pasapas tab1 ;
  231.  
  232. * on se contente de verifier que le calcul est bien passe.
  233.  
  234. uu = tab1.deplacements.38;
  235. xx = extr uu ux p16;
  236. valref = 3.22079E-03;
  237. list xx ;
  238. si (abs (xx - valref) > 1e-3); erreur 5; finsi;
  239. list tab1.temps;
  240. fin;
  241.  
  242.  
  243.  
  244.  
  245.  
  246.  
  247.  
  248.  
  249.  
  250.  
  251.  
  252.  
  253.  

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