Télécharger phasage.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : phasage.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. complet=vrai;
  5. opti dime 3 elem cu20 mode trid;
  6. am=-1; bm = -0.5; cm=-0.275;cp=0.275;bp=0.5;ap=1.;
  7.  
  8. *** nca2=10;nca=nca2*2;nd1=2; nd2=4;
  9. si complet ;
  10. nca2=5;nca=nca2*2;nd1=1; nd2=2;
  11. sinon;
  12. nca2=4;nca=nca2*2;nd1=1; nd2=2;
  13. finsi;
  14. pa = am am 20;pb=bm am 20; pc= cm am 20; pd = cp am 20;
  15. pe= bp am 20; pf = ap am 20;
  16. liba = pa droi nd1 pb droi nd1 pc droi nd2 pd droi nd1 pe droi nd1 pf;
  17. su = liba tran nd1 ( 0 bp 0) trans nd1 ( 0 0.225 0.) trans nd2 (0 0.55 0) tran nd1 ( 0 0.225 0) tran nd1 ( 0 BP 0);
  18.  
  19. *trac su;
  20.  
  21. vol1= su volu trans nca2 ( 0 0 -10);
  22. su2 = vol1 face 2;
  23. pceny = point su2 droite ( 0 0 10) ( 1 0 10) 0.51;
  24. pcenxy = point pceny droite ( 0 0 10) ( 0 1 10) 0.51;
  25.  
  26. su21= elem su2 appu strict pcenxy coul rouge;
  27.  
  28. *trac ( su2 et su21);
  29.  
  30. vol2 = su21 volu trans nca2 (0 0 -10.);
  31. su3 = vol2 face 2;
  32.  
  33. voltot = vol1 et vol2;
  34.  
  35. *trac voltot;
  36.  
  37. pc1d = cm cp 0.; pc1f = cm cp 20;
  38. cab1 = pc1d droi nca pc1f chan seg2;
  39. pc2d = cp cm 0.; pc2f = cp cm 20;
  40. cab2 = pc2d droi nca pc2f chan seg2;
  41. pc3d = cp cp 0.; pc3f = cp cp 20;
  42. cab3 = pc3d droi nca pc3f chan seg2;
  43. pc4d = cm cm 0.; pc4f = cm cm 20;
  44. cab4 = pc4d droi nca pc4f chan seg2;
  45. pc5d = 0 0 0.; pc5f = 0 0 20;
  46. cab5 = pc5d droi nca pc5f chan seg2;
  47.  
  48. cable= cab1 et cab2 et cab3 et cab4 et cab5 coul bleu;
  49. fer= cable plus (0 0 0 );
  50. list ( nbno cable);
  51.  
  52. *elim voltot cable 0.01;
  53.  
  54.  
  55. *trac ( vol2 et cable);
  56.  
  57.  
  58. *
  59. * description des modeles
  60. *
  61. mobet1 = model vol1 mecanique elastique;
  62. mobet2 = model vol2 mecanique elastique;
  63. mobeton = mobet1 et mobet2;
  64.  
  65. mocab1= model cab1 mecanique elastique barr;
  66. mocab2= model cab2 mecanique elastique barr;
  67. mocab3= model cab3 mecanique elastique barr;
  68. mocab4= model cab4 mecanique elastique barr;
  69. mocab5= model cab5 mecanique elastique barr;
  70.  
  71. mofer = model fer mecanique elastique barr;
  72. *
  73. * description des materiaux et sections
  74. *
  75.  
  76. mabet1= mate mobet1 young 40000e6 nu 0.2 rho 2500.;
  77. mabet2= mate mobet2 young 40000e6 nu 0.2 rho 2500.;
  78. mabeton = mabet1 et mabet2;
  79.  
  80. macab1= mate mocab1 young 193000e6 nu 0.3 rho 7859 sect 0.0025;
  81. macab2= mate mocab2 young 193000e6 nu 0.3 rho 7859 sect 0.0025;
  82. macab3= mate mocab3 young 193000e6 nu 0.3 rho 7859 sect 0.0025;
  83. macab4= mate mocab4 young 193000e6 nu 0.3 rho 7859 sect 0.0025;
  84. macab5= mate mocab5 young 193000e6 nu 0.3 rho 7859 sect 0.0025;
  85.  
  86. mafer = mate mofer young 100000e6 nu 0.3 rho 7859 sect 0.00001;
  87. *
  88. * condition aux limites de deplacement
  89. *
  90.  
  91. pp= point voltot droite (0 0 0 )( 0 0 20) 0.01;
  92. blxy = bloqu 'UX' 'UY' pp;
  93. px= point su3 proche (0.5 0 0 );
  94. bly = bloqu 'UY' px;
  95. blz = bloqu 'UZ' su3;
  96. claccr = rela accro fer voltot 0.01;
  97. cltot = blxy et bly et blz et claccr;
  98.  
  99. *
  100. * chargement du au poids
  101. *
  102.  
  103. gamma = -9.81;
  104. mmas = masse mobeton mabeton;
  105. vz = manu chpo voltot 1 'UZ' gamma;
  106. fpoi= mmas * vz;
  107. evpoi = evol manu 'TEMPS' (prog 0. 10000.) 'COEF' (PROG 1. 1.);
  108. charpoid= chargement fpoi evpoi;
  109.  
  110. *
  111. * description des caracteristioques generales pour les pertes quasi instantanee
  112. *
  113.  
  114. coefprec=table;
  115.  
  116. coefprec . 'FF' = 0.16 ;
  117. coefprec . 'PHIF' = 0.0015;
  118. coefprec . 'GANC' = 0.012;
  119. coefprec . 'RMU0' = 0.45 ;
  120. coefprec . 'FPRG' = 1940.e6 ;
  121. coefprec . 'RH10' = 2.3;
  122.  
  123. *
  124. * creation de la table des etapes pour le cas 2 il n'y a q'une etape
  125. *
  126. cas3etap=table;
  127. cas2=table;
  128. cas3etap . 1 = cas2 ;
  129. cas2 . 'TPS' = 300. ;
  130.  
  131. ********************************** premier groupe
  132.  
  133. group1_1 = table;
  134. cas2 . 1 = group1_1 ;
  135.  
  136. group1_1 . 'GEOMETRIE1' = manu poi1 pc2d;
  137. group1_1 . 'MODELE' = mocab2;
  138. group1_1 . 'MATERIAU' = macab2;
  139. group1_1 . 'FORCE' = 3.75e6;
  140. group1_1 . 'COEF_PREC' = coefprec;
  141. group1_1 . 'TYPE_CAB' = '1EXT';
  142.  
  143. group1_2= table;
  144. cas2 . 2 = group1_2 ;
  145.  
  146. group1_2 . 'GEOMETRIE1' = manu poi1 pc1d;
  147. group1_2 . 'MODELE' = mocab1;
  148. group1_2 . 'MATERIAU' = macab1;
  149. group1_2 . 'FORCE' = 3.75e6;
  150. group1_2 . 'COEF_PREC' = coefprec;
  151. group1_2 . 'TYPE_CAB' = '1EXT';
  152.  
  153.  
  154. ***************************** deuxieme groupe
  155. si complet;
  156.  
  157. cas3=table;
  158. cas3etap . 2 = cas3;
  159.  
  160. cas3.'TPS'= 450.;
  161. gr31=table;
  162. cas3 . 1 = gr31;
  163. gr31. 'GEOMETRIE1' = manu poi1 pc3d ;
  164. gr31. 'MODELE' = mocab3;
  165. gr31. 'MATERIAU' = macab3;
  166. gr31. 'FORCE' = 3.75e6;
  167. gr31. 'COEF_PREC' = coefprec;
  168. gr31. 'TYPE_CAB' = '1EXT';
  169.  
  170. gr32=table;
  171. cas3 . 2 = gr32;
  172. gr32. 'GEOMETRIE1' = manu poi1 pc4d ;
  173. gr32. 'MODELE' = mocab4;
  174. gr32. 'MATERIAU' = macab4;
  175. gr32. 'FORCE' = 3.75e6;
  176. gr32. 'COEF_PREC' = coefprec;
  177. gr32. 'TYPE_CAB' = '1EXT';
  178.  
  179. **************************** troisieme groupe
  180.  
  181. cas4=table;
  182. cas3etap . 3 = cas4;
  183. cas4.'TPS'= 600.;
  184. gr41=table;
  185. cas4 . 1 = gr41;
  186. gr41. 'GEOMETRIE1' = manu poi1 pc5d ;
  187. gr41. 'GEOMETRIE2' = manu poi1 pc5f ;
  188. gr41. 'MODELE' = mocab5;
  189. gr41. 'MATERIAU' = macab5;
  190. gr41. 'FORCE' = 3.75e6;
  191. gr41. 'COEF_PREC' = coefprec;
  192. gr41. 'TYPE_CAB' = '2EXT';
  193.  
  194.  
  195. finsi;
  196.  
  197. *
  198. * appel a la procedur tension
  199. *
  200.  
  201. cas3etap = tension cas3etap;
  202. *
  203. * donnees relatives au levee
  204. *
  205. tlev = table;tlev1=table;tlev2=table;
  206. tlev . 1 = tlev1;tlev . 2 = tlev2;
  207. tlev1.'INSTANT' = 0.; tlev2.'INSTANT' = 150.;
  208. tlev1.'MODELE' = mobet1;tlev2.'MODELE' = mobet2;
  209. tlev1.'MATERIAU' = mabet1;tlev2.'MATERIAU' = mabet2;
  210. tlev1.'COEF1'=70.;tlev2.'COEF1'=70.;
  211. tlev1.'COEF2'= MANU CHML mobet1 'TAUX' 9.e-3 'STRESSES';
  212. tlev2.'COEF2'= MANU CHML mobet2 'TAUX' 9.e-3 'STRESSES';
  213. raysec1 = manu chml mobet1 'EPAI' 50. 'STRESSES';
  214. raysec2 = manu chml mobet2 'EPAI' 100. 'STRESSES';
  215. raysec = raysec1 et raysec2;
  216. tlev1.'SECHAGE' = raysec1;tlev2.'SECHAGE' = raysec2;
  217.  
  218.  
  219. *
  220. * appel a phasage
  221. *
  222.  
  223. tphas= table;
  224.  
  225. tphas.'FLUAGE' = 'BPEL99';
  226. tphas.'RETRAIT' = 'BPEL99';
  227. tphas.'PRECONTRAINTE'= cas3etap;
  228. tphas.'LEVEES'= tlev;
  229. tphas. 'TEMPS_FINAL' = 3000.;
  230. tphas.'BLOCAGES' = cltot;
  231. tphas.'DEPOU'= vrai;
  232. tphas.'MOD_RESTE' = mofer;
  233. tphas.'MAT_RESTE' = mafer;
  234. *option debug ;
  235. *option echo 2;
  236.  
  237.  
  238.  
  239.  
  240. tphas = phasage TPHAS;
  241.  
  242.  
  243. *opti sauv 'complet.res';
  244.  
  245. *sauv tphas;
  246.  
  247. *fin;
  248.  
  249. tsui = tphas.table_suite;
  250. ttem= tsui.temps;
  251. tdep= tsui.deplacements;
  252. tcon= tsui.contraintes;
  253. na = (dime ttem) - 1;
  254. maxdepz= maxi( abs (exco tdep . na UZ));
  255. maxcoz=maxi ( abs ( exco tcon . na SMZZ));
  256. si complet;
  257. message ' calcul complet = vrai';
  258. verdep=1.1017; vericont=1.86069e+07;
  259. sinon ;message ' calcul complet = faux';
  260. verdep= 1.0949 ; vericont= 8.37572e+06;
  261. finsi;
  262. mess ' deplacemment en Z maxi ' maxdepz ' comparé à ' verdep;
  263. mess ' contrainte en Z maxi ' maxcoz ' comparée à ' vericont;
  264.  
  265.  
  266. erdep= abs((maxdepz - verdep) / verdep );
  267. ercont = abs ((maxcoz - vericont) / vericont);
  268.  
  269. si ( erdep > 0.001) ;
  270.  
  271. mess ' erreur deplacement ' erdep ' comparé à 0.001';
  272. erreur 5;
  273. finsi;
  274. si ( ercont > 0.001 );
  275.  
  276. mess ' errreur contrainte ' ercont ' comparé à 0.001';
  277. erreur 5;
  278. finsi;
  279.  
  280. fin;
  281.  
  282.  
  283.  
  284.  
  285.  
  286.  
  287.  
  288.  
  289.  
  290.  
  291.  
  292.  
  293.  
  294.  
  295.  
  296.  
  297.  
  298.  
  299.  
  300.  
  301.  
  302.  
  303.  
  304.  
  305.  
  306.  
  307.  
  308.  
  309.  
  310.  

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