************************************************************************ * Modelisation du cuvelage inox d'une galerie * ************************************************************************ * ******************************* * * SOMMAIRE DU PROGRAMME * * * **************** * * * * * * 1/Construction tole et sinu * * * 2/ Galerie * * * 3/ Caractér et charges * * * 4/ Traitements * * * 5/ Resulats * * ******************************* ************************************************************************ *option opti echo 1; opti dime 3 elem tri3 ; *données en m larg = 2.5e-2; long = 1.2; ray = 0.95; ray2 = 0.99; ep = 1.5e-3 ; rayg = 1. ; *petite longueur en m pl = 0.06 ; *pas des sinusoides en m ps = 0.10; *hauteur sinusoide en m hs = 0.04; *projection long en m pro = 0.04; *densité densi larg ; *pression en kg/m²* p = 0.04e7; *young en MPA you1 = 2.e11; *Coeff poisson nu1 = 0.3; *Contrainte limite elastique en MPa sig1 = 320. ; *Poids volumique en kg/m3 rho1 = 7700 ; *section tole sect1 = pi*(((ray + ep)**2)-(ray**2)) ; *section sinusoide sect2 = ep * larg ; ************************************************************************ * 1/ Construction des connecteurs * ************************************************************************ o1 = 0. 0. 0. ; o2 = 0. long 0. ; *connecteur1 p0 = 0. 0. ray ; p1 = 0. pl ray ; p2 = 0. (pl + (2* pro) + ps ) ray ; p3 = 0. (pl + (2* pro) + (2* ps)) ray ; p4 = 0. (pl + (4* pro) + (3* ps)) ray ; p5 = 0. (pl + (4* pro) + (4* ps)) ray ; p6 = 0. (pl + (6* pro) + (5* ps)) ray ; p7 = 0. (pl + (6* pro) + (6* ps)) ray ; p8 = 0. (pl + (8* pro) + (7* ps)) ray ; p9 = 0. ((2* pl) + (8* pro) + (7* ps)) ray ; p11 = 0. (pl + pro) ray2 ; p12 = 0. (pl + pro + ps) ray2 ; p31 = 0. (pl + (3* pro) + (2* ps)) ray2 ; p32 = 0. (pl + (3* pro) + (3* ps)) ray2 ; p51 = 0. (pl + (5* pro) + (4* ps)) ray2 ; p52 = 0. (pl + (5* pro) + (5* ps)) ray2 ; p71 = 0. (pl + (7* pro) + (6* ps)) ray2 ; p72 = 0. (pl + (7* pro) + (7* ps)) ray2 ; p00 p01 p02 p03 p04 p05 p06 p07 p08 p09 p011 p012 p031 p032 p051 p052 p071 p072 = p0 p1 p2 p3 p4 p5 p6 p7 p8 p9 p11 p12 p31 p32 p51 p52 p71 p72 tour 1.5 o1 o2; l1 = droi p0 p00 ; l2 = droi p1 p01 ; l3 = droi p2 p02 ; l4 = droi p3 p03 ; l5 = droi p4 p04 ; l6 = droi p5 p05 ; l7 = droi p6 p06 ; l8 = droi p7 p07 ; l9 = droi p8 p08 ; l10 = droi p9 p09 ; l11 = droi p11 p011 ; l12 = droi p12 p012 ; l13 = droi p31 p031 ; l14 = droi p32 p032 ; l15 = droi p51 p051 ; l16 = droi p52 p052 ; l17 = droi p71 p071 ; l18 = droi p72 p072 ; l01 = droi p0 p1 ; l02 = droi p1 p11 ; l03 = droi p11 p12 ; l04 = droi p12 p2 ; l05 = droi p2 p3 ; l06 = droi p3 p31 ; l07 = droi p31 p32 ; l08 = droi p32 p4 ; l09 = droi p4 p5 ; l010 = droi p5 p51 ; l011 = droi p51 p52 ; l012 = droi p52 p6 ; l013 = droi p6 p7 ; l014 = droi p7 p71 ; l015 = droi p71 p72 ; l016 = droi p72 p8 ; l017 = droi p8 p9 ; l018 = droi p00 p01 ; l019 = droi p01 p011 ; l020 = droi p011 p012 ; l021 = droi p012 p02 ; l022 = droi p02 p03 ; l023 = droi p03 p031 ; l024 = droi p031 p032 ; l025 = droi p032 p04 ; l026 = droi p04 p05 ; l027 = droi p05 p051 ; l028 = droi p051 p052 ; l029 = droi p052 p06 ; l030 = droi p06 p07 ; l031 = droi p07 p071 ; l032 = droi p071 p072 ; l033 = droi p072 p08 ; l034 = droi p08 p09 ; s1 = dall l1 l01 l2 l018 cyli o1 o2 ; s2 = dall l2 l02 l11 l019 plan ; s3 = dall l11 l03 l12 l020 cyli o1 o2 ; s4 = dall l12 l04 l3 l021 plan ; s5 = dall l3 l05 l4 l022 cyli o1 o2 ; s6 = dall l4 l06 l13 l023 plan ; s7 = dall l13 l07 l14 l024 cyli o1 o2 ; s8 = dall l14 l08 l5 l025 plan ; s9 = dall l5 l09 l6 l026 cyli o1 o2 ; s10 = dall l6 l010 l15 l027 plan ; s11 = dall l15 l011 l16 l028 cyli o1 o2 ; s12 = dall l16 l012 l7 l029 plan ; s13 = dall l7 l013 l8 l030 cyli o1 o2 ; s14 = dall l8 l014 l17 l031 plan ; s15 = dall l17 l015 l18 l032 cyli o1 o2 ; s16 = dall l18 l016 l9 l033 plan ; s17 = dall l9 l017 l10 l034 cyli o1 o2 ; *assemblage sinusoide sinu = s1 et s2 et s3 et s4 et s5 et s6 et s7 et s8 et s9 et s10 et s11 et s12 et s13 et s14 et s15 et s16 et s17; trac sinu titr'Sinusoide'; *connecteur 2 à 15 planp0 l201 = l1 tour 24 o1 o2 ; l301 = l1 tour 48 o1 o2 ; l401 = l1 tour 72 o1 o2 ; l501 = l1 tour 96 o1 o2 ; l601 = l1 tour 120 o1 o2 ; l701 = l1 tour 144 o1 o2 ; l801 = l1 tour 168 o1 o2 ; l901 = l1 tour 192 o1 o2 ; l1001 = l1 tour 216 o1 o2 ; l1101 = l1 tour 240 o1 o2 ; l1201 = l1 tour 264 o1 o2 ; l1301 = l1 tour 288 o1 o2 ; l1401 = l1 tour 312 o1 o2 ; l1501 = l1 tour 336 o1 o2 ; c1 = c l1 o1 l201 ; c12 = c l201 o1 l301 ; c13 = c l301 o1 l401 ; c14 = c l401 o1 l501 ; c15 = c l501 o1 l601 ; c16 = c l601 o1 l701 ; c17 = c l701 o1 l801 ; c18 = c l801 o1 l901 ; c19 = c l901 o1 l1001 ; c110 = c l1001 o1 l1101 ; c111 = c l1101 o1 l1201 ; c112 = c l1201 o1 l1301 ; c113 = c l1301 o1 l1401 ; c114 = c l1401 o1 l1501 ; c115 = c l1501 o1 l1 ; cer1 = (c1 et c12 et c13 et c14 et c15 et c16 et c17 et c18 et c19 et c110 et c111 et c112 et c113 et c114 et c115); *connecteur 2 à 15 planp1 l202 = l2 tour 24 o1 o2 ; l302 = l2 tour 48 o1 o2 ; l402 = l2 tour 72 o1 o2 ; l502 = l2 tour 96 o1 o2 ; l602 = l2 tour 120 o1 o2 ; l702 = l2 tour 144 o1 o2 ; l802 = l2 tour 168 o1 o2 ; l902 = l2 tour 192 o1 o2 ; l1002 = l2 tour 216 o1 o2 ; l1102 = l2 tour 240 o1 o2 ; l1202 = l2 tour 264 o1 o2 ; l1302 = l2 tour 288 o1 o2 ; l1402 = l2 tour 312 o1 o2 ; l1502 = l2 tour 336 o1 o2 ; o3 = 0. pl 0.; c2 = c l2 o3 l202 ; c22 = c l202 o3 l302 ; c23 = c l302 o3 l402 ; c24 = c l402 o3 l502 ; c25 = c l502 o3 l602 ; c26 = c l602 o3 l702 ; c27 = c l702 o3 l802 ; c28 = c l802 o3 l902 ; c29 = c l902 o3 l1002 ; c210 = c l1002 o3 l1102 ; c211 = c l1102 o3 l1202 ; c212 = c l1202 o3 l1302 ; c213 = c l1302 o3 l1402 ; c214 = c l1402 o3 l1502 ; c215 = c l1502 o3 l2 ; cer2 = (c2 et c22 et c23 et c24 et c25 et c26 et c27 et c28 et c29 et c210 et c211 et c212 et c213 et c214 et c215); *connecteur 2 à 15 planp2 l203 = l3 tour 24 o1 o2 ; l303 = l3 tour 48 o1 o2 ; l403 = l3 tour 72 o1 o2 ; l503 = l3 tour 96 o1 o2 ; l603 = l3 tour 120 o1 o2 ; l703 = l3 tour 144 o1 o2 ; l803 = l3 tour 168 o1 o2 ; l903 = l3 tour 192 o1 o2 ; l1003 = l3 tour 216 o1 o2 ; l1103 = l3 tour 240 o1 o2 ; l1203 = l3 tour 264 o1 o2 ; l1303 = l3 tour 288 o1 o2 ; l1403 = l3 tour 312 o1 o2 ; l1503 = l3 tour 336 o1 o2 ; pl4 = pl + 2 * pro + ps; o4 = 0. pl4 0.; c3 = c l3 o4 l203 ; c32 = c l203 o4 l303 ; c33 = c l303 o4 l403 ; c34 = c l403 o4 l503 ; c35 = c l503 o4 l603 ; c36 = c l603 o4 l703 ; c37 = c l703 o4 l803 ; c38 = c l803 o4 l903 ; c39 = c l903 o4 l1003 ; c310 = c l1003 o4 l1103 ; c311 = c l1103 o4 l1203 ; c312 = c l1203 o4 l1303 ; c313 = c l1303 o4 l1403 ; c314 = c l1403 o4 l1503 ; c315 = c l1503 o4 l3 ; cerr3 = (c3 et c32 et c33 et c34 et c35 et c36 et c37 et c38 et c39 et c310 et c311 et c312 et c313 et c314 et c315); *connecteur 2 à 15 planp3 l204 = l4 tour 24 o1 o2 ; l304 = l4 tour 48 o1 o2 ; l404 = l4 tour 72 o1 o2 ; l504 = l4 tour 96 o1 o2 ; l604 = l4 tour 120 o1 o2 ; l704 = l4 tour 144 o1 o2 ; l804 = l4 tour 168 o1 o2 ; l904 = l4 tour 192 o1 o2 ; l1004 = l4 tour 216 o1 o2 ; l1104 = l4 tour 240 o1 o2 ; l1204 = l4 tour 264 o1 o2 ; l1304 = l4 tour 288 o1 o2 ; l1404 = l4 tour 312 o1 o2 ; l1504 = l4 tour 336 o1 o2 ; pl5 = pl + 2*pro + 2* ps ; o5 = 0. pl5 0.; c4 = c l4 o5 l204 ; c42 = c l204 o5 l304 ; c43 = c l304 o5 l404 ; c44 = c l404 o5 l504 ; c45 = c l504 o5 l604 ; c46 = c l604 o5 l704 ; c47 = c l704 o5 l804 ; c48 = c l804 o5 l904 ; c49 = c l904 o5 l1004 ; c410 = c l1004 o5 l1104 ; c411 = c l1104 o5 l1204 ; c412 = c l1204 o5 l1304 ; c413 = c l1304 o5 l1404 ; c414 = c l1404 o5 l1504 ; c415 = c l1504 o5 l4 ; cer4 = (c4 et c42 et c43 et c44 et c45 et c46 et c47 et c48 et c49 et c410 et c411 et c412 et c413 et c414 et c415); *connecteur 2 à 15 planp4 l205 = l5 tour 24 o1 o2 ; l305 = l5 tour 48 o1 o2 ; l405 = l5 tour 72 o1 o2 ; l505 = l5 tour 96 o1 o2 ; l605 = l5 tour 120 o1 o2 ; l705 = l5 tour 144 o1 o2 ; l805 = l5 tour 168 o1 o2 ; l905 = l5 tour 192 o1 o2 ; l1005 = l5 tour 216 o1 o2 ; l1105 = l5 tour 240 o1 o2 ; l1205 = l5 tour 264 o1 o2 ; l1305 = l5 tour 288 o1 o2 ; l1405 = l5 tour 312 o1 o2 ; l1505 = l5 tour 336 o1 o2 ; pl6 = pl + 4 * pro + 3* ps ; o6 = 0. pl6 0.; c5 = c l5 o6 l205 ; c52 = c l205 o6 l305 ; c53 = c l305 o6 l405 ; c54 = c l405 o6 l505 ; c55 = c l505 o6 l605 ; c56 = c l605 o6 l705 ; c57 = c l705 o6 l805 ; c58 = c l805 o6 l905 ; c59 = c l905 o6 l1005 ; c510 = c l1005 o6 l1105 ; c511 = c l1105 o6 l1205 ; c512 = c l1205 o6 l1305 ; c513 = c l1305 o6 l1405 ; c514 = c l1405 o6 l1505 ; c515 = c l1505 o6 l5 ; cer5 = (c5 et c52 et c53 et c54 et c55 et c56 et c57 et c58 et c59 et c510 et c511 et c512 et c513 et c514 et c515); *connecteur 2 à 15 planp5 l206 = l6 tour 24 o1 o2 ; l306 = l6 tour 48 o1 o2 ; l406 = l6 tour 72 o1 o2 ; l506 = l6 tour 96 o1 o2 ; l606 = l6 tour 120 o1 o2 ; l706 = l6 tour 144 o1 o2 ; l806 = l6 tour 168 o1 o2 ; l906 = l6 tour 192 o1 o2 ; l1006 = l6 tour 216 o1 o2 ; l1106 = l6 tour 240 o1 o2 ; l1206 = l6 tour 264 o1 o2 ; l1306 = l6 tour 288 o1 o2 ; l1406 = l6 tour 312 o1 o2 ; l1506 = l6 tour 336 o1 o2 ; pl7 = pl + 4* pro + 4* ps ; o7 = 0. pl7 0.; c6 = c l6 o7 l206 ; c62 = c l206 o7 l306 ; c63 = c l306 o7 l406 ; c64 = c l406 o7 l506 ; c65 = c l506 o7 l606 ; c66 = c l606 o7 l706 ; c67 = c l706 o7 l806 ; c68 = c l806 o7 l906 ; c69 = c l906 o7 l1006 ; c610 = c l1006 o7 l1106 ; c611 = c l1106 o7 l1206 ; c612 = c l1206 o7 l1306 ; c613 = c l1306 o7 l1406 ; c614 = c l1406 o7 l1506 ; c615 = c l1506 o7 l6 ; cer6 = (c6 et c62 et c63 et c64 et c65 et c66 et c67 et c68 et c69 et c610 et c611 et c612 et c613 et c614 et c615); *connecteur 2 à 15 planp6 l207 = l7 tour 24 o1 o2 ; l307 = l7 tour 48 o1 o2 ; l407 = l7 tour 72 o1 o2 ; l507 = l7 tour 96 o1 o2 ; l607 = l7 tour 120 o1 o2 ; l707 = l7 tour 144 o1 o2 ; l807 = l7 tour 168 o1 o2 ; l907 = l7 tour 192 o1 o2 ; l1007 = l7 tour 216 o1 o2 ; l1107 = l7 tour 240 o1 o2 ; l1207 = l7 tour 264 o1 o2 ; l1307 = l7 tour 288 o1 o2 ; l1407 = l7 tour 312 o1 o2 ; l1507 = l7 tour 336 o1 o2 ; pl8 = pl + (6* pro) + (5* ps) ; o8 = 0. pl8 0.; c7 = c l7 o8 l207 ; c72 = c l207 o8 l307 ; c73 = c l307 o8 l407 ; c74 = c l407 o8 l507 ; c75 = c l507 o8 l607 ; c76 = c l607 o8 l707 ; c77 = c l707 o8 l807 ; c78 = c l807 o8 l907 ; c79 = c l907 o8 l1007 ; c710 = c l1007 o8 l1107 ; c711 = c l1107 o8 l1207 ; c712 = c l1207 o8 l1307 ; c713 = c l1307 o8 l1407 ; c714 = c l1407 o8 l1507 ; c715 = c l1507 o8 l7 ; cer7 = (c7 et c72 et c73 et c74 et c75 et c76 et c77 et c78 et c79 et c710 et c711 et c712 et c713 et c714 et c715); *connecteur 2 à 15 planp7 l208 = l8 tour 24 o1 o2 ; l308 = l8 tour 48 o1 o2 ; l408 = l8 tour 72 o1 o2 ; l508 = l8 tour 96 o1 o2 ; l608 = l8 tour 120 o1 o2 ; l708 = l8 tour 144 o1 o2 ; l808 = l8 tour 168 o1 o2 ; l908 = l8 tour 192 o1 o2 ; l1008 = l8 tour 216 o1 o2 ; l1108 = l8 tour 240 o1 o2 ; l1208 = l8 tour 264 o1 o2 ; l1308 = l8 tour 288 o1 o2 ; l1408 = l8 tour 312 o1 o2 ; l1508 = l8 tour 336 o1 o2 ; pl9 = pl + (6* pro) + (6 * ps) ; o9 = 0. pl9 0.; c8 = c l8 o9 l208 ; c82 = c l208 o9 l308 ; c83 = c l308 o9 l408 ; c84 = c l408 o9 l508 ; c85 = c l508 o9 l608 ; c86 = c l608 o9 l708 ; c87 = c l708 o9 l808 ; c88 = c l808 o9 l908 ; c89 = c l908 o9 l1008 ; c810 = c l1008 o9 l1108 ; c811 = c l1108 o9 l1208 ; c812 = c l1208 o9 l1308 ; c813 = c l1308 o9 l1408 ; c814 = c l1408 o9 l1508 ; c815 = c l1508 o9 l8 ; cer8 = (c8 et c82 et c83 et c84 et c85 et c86 et c87 et c88 et c89 et c810 et c811 et c812 et c813 et c814 et c815); *connecteur 2 à 15 planp8 l209 = l9 tour 24 o1 o2 ; l309 = l9 tour 48 o1 o2 ; l409 = l9 tour 72 o1 o2 ; l509 = l9 tour 96 o1 o2 ; l609 = l9 tour 120 o1 o2 ; l709 = l9 tour 144 o1 o2 ; l809 = l9 tour 168 o1 o2 ; l909 = l9 tour 192 o1 o2 ; l1009 = l9 tour 216 o1 o2 ; l1109 = l9 tour 240 o1 o2 ; l1209 = l9 tour 264 o1 o2 ; l1309 = l9 tour 288 o1 o2 ; l1409 = l9 tour 312 o1 o2 ; l1509 = l9 tour 336 o1 o2 ; pl10 = pl + (8* pro) + (7* ps) ; o10 = 0. pl10 0.; c9 = c l9 o10 l209 ; c92 = c l209 o10 l309 ; c93 = c l309 o10 l409 ; c94 = c l409 o10 l509 ; c95 = c l509 o10 l609 ; c96 = c l609 o10 l709 ; c97 = c l709 o10 l809 ; c98 = c l809 o10 l909 ; c99 = c l909 o10 l1009 ; c910 = c l1009 o10 l1109 ; c911 = c l1109 o10 l1209 ; c912 = c l1209 o10 l1309 ; c913 = c l1309 o10 l1409 ; c914 = c l1409 o10 l1509 ; c915 = c l1509 o10 l9 ; cer9 = (c9 et c92 et c93 et c94 et c95 et c96 et c97 et c98 et c99 et c910 et c911 et c912 et c913 et c914 et c915); *connecteur 2 à 15 planp9 l210 = l10 tour 24 o1 o2 ; l310 = l10 tour 48 o1 o2 ; l410 = l10 tour 72 o1 o2 ; l510 = l10 tour 96 o1 o2 ; l610 = l10 tour 120 o1 o2 ; l710 = l10 tour 144 o1 o2 ; l810 = l10 tour 168 o1 o2 ; l910 = l10 tour 192 o1 o2 ; l1010 = l10 tour 216 o1 o2 ; l1110 = l10 tour 240 o1 o2 ; l1210 = l10 tour 264 o1 o2 ; l1310 = l10 tour 288 o1 o2 ; l1410 = l10 tour 312 o1 o2 ; l1510 = l10 tour 336 o1 o2 ; *rotation sinu sinu2 = sinu tour 24 o1 o2 ; sinu3 = sinu tour 48 o1 o2 ; sinu4 = sinu tour 72 o1 o2 ; sinu5 = sinu tour 96 o1 o2 ; sinu6 = sinu tour 120 o1 o2 ; sinu7 = sinu tour 144 o1 o2 ; sinu8 = sinu tour 168 o1 o2 ; sinu9 = sinu tour 192 o1 o2 ; sinu10 = sinu tour 216 o1 o2 ; sinu11 = sinu tour 240 o1 o2 ; sinu12 = sinu tour 264 o1 o2 ; sinu13 = sinu tour 288 o1 o2 ; sinu14 = sinu tour 312 o1 o2 ; sinu15 = sinu tour 336 o1 o2 ; sinutot = (sinu et sinu2 et sinu3 et sinu4 et sinu5 et sinu6 et sinu7 et sinu8 et sinu9 et sinu10 et sinu11 et sinu12 et sinu13 et sinu14 et sinu15) coul rouge; *assemblage tole lij = ((2* pro) + ps); ce1 = cer1 trans (0. pl 0.) ; ce2 = cer2 trans (0. lij 0.) ; ce3 = cerr3 trans (0. ps 0.) ; ce4 = cer4 trans (0. lij 0.) ; ce5 = cer5 trans (0. ps 0.) ; ce6 = cer6 trans (0. lij 0.) ; ce7 = cer7 trans (0. ps 0.) ; ce8 = cer8 trans (0. lij 0.) ; ce9 = cer9 trans (0. pl 0.) ; tole = (ce1 et ce2 et ce3 et ce4 et ce5 et ce6 et ce7 et ce8 et ce9) coul bleu; elim (tole et sinutot); titr 'Tole Inox avec Sinusoides'; RP1 = @REPERE O1 (prog 0.2 0.2 0.2) ; trac (tole et sinutot et rp1); *Conditions aux limites* ************************ *sinusoides bas* **************** pat1 = (pl + pro) ; o11 = 0. pat1 0. ; climtr1 = l11 tour 168 o1 o2 ; climtr2 = l11 tour 192 o1 o2 ; clim = c climtr1 o11 climtr2 ; clim1 = clim trans (0. (-1. * pat1) 0.) ; clim2 = clim trans (0. (long) 0.) ; climtot = clim1 et clim2 ; elim (climtot et sinu8 et sinu9) ; sinucl = (sinu8 et sinu9) elem appuy strict climtot ; pat2 = (pl + (3* pro) + (2* ps )) ; o31 = 0. pat2 0. ; climtr12 = l13 tour 168 o1 o2 ; climtr22 = l13 tour 192 o1 o2 ; clim2 = c climtr12 o31 climtr22 ; clim12 = clim2 trans (0. (-1. * pat2) 0.) ; clim22 = clim2 trans (0. (long) 0.) ; climtot2 = clim12 et clim22 ; elim (climtot2 et sinu8 et sinu9) ; sinucl2 = (sinu8 et sinu9) elem appuy strict climtot2 ; pat3 = (pl + (5* pro) + (4* ps )) ; o51 = 0. pat3 0. ; climtr13 = l15 tour 168 o1 o2 ; climtr23 = l15 tour 192 o1 o2 ; clim3 = c climtr13 o31 climtr23 ; clim13 = clim3 trans (0. (-1. * pat3) 0.) ; clim23 = clim3 trans (0. (long) 0.) ; climtot3 = clim13 et clim23 ; elim (climtot3 et sinu8 et sinu9) ; sinucl3 = (sinu8 et sinu9) elem appuy strict climtot3 ; pat4 = (pl + (7* pro) + (6* ps )) ; o71 = 0. pat4 0. ; climtr14 = l17 tour 168 o1 o2 ; climtr24 = l17 tour 192 o1 o2 ; clim4 = c climtr14 o71 climtr24 ; clim14 = clim4 trans (0. (-1. * pat4) 0.) ; clim24 = clim4 trans (0. (long) 0.) ; climtot4 = clim14 et clim24 ; elim (climtot4 et sinu8 et sinu9) ; sinucl4 = (sinu8 et sinu9) elem appuy strict climtot4 ; stot = (sinucl et sinucl2 et sinucl3 et sinucl4 ); trac ((stot coul jaune) et tole et sinutot) titr 'Blocages inferieurs'; cls = bloq stot ux uy uz rx ry rz; *galerie* *points galerie* **************** pg0 = 0. 0. rayg ; pg1 = 0. pl rayg ; pg2 = 0. (pl + (2* pro) + ps ) rayg ; pg3 = 0. (pl + (2* pro) + (2* ps)) rayg ; pg4 = 0. (pl + (4* pro) + (3* ps)) rayg ; pg5 = 0. (pl + (4* pro) + (4* ps)) rayg ; pg6 = 0. (pl + (6* pro) + (5* ps)) rayg ; pg7 = 0. (pl + (6* pro) + (6* ps)) rayg ; pg8 = 0. (pl + (8* pro) + (7* ps)) rayg ; pg9 = 0. ((2* pl) + (8* pro) + (7* ps)) rayg ; pg11 = 0. (pl + pro) rayg ; pg12 = 0. (pl + pro + ps) rayg ; pg31 = 0. (pl + (3* pro) + (2* ps)) rayg ; pg32 = 0. (pl + (3* pro) + (3* ps)) rayg ; pg51 = 0. (pl + (5* pro) + (4* ps)) rayg ; pg52 = 0. (pl + (5* pro) + (5* ps)) rayg ; pg71 = 0. (pl + (7* pro) + (6* ps)) rayg ; pg72 = 0. (pl + (7* pro) + (7* ps)) rayg ; pg00 pg01 pg02 pg03 pg04 pg05 pg06 pg07 pg08 pg09 pg011 pg012 pg031 pg032 pg051 pg052 pg071 pg072 = pg0 pg1 pg2 pg3 pg4 pg5 pg6 pg7 pg8 pg9 pg11 pg12 pg31 pg32 pg51 pg52 pg71 pg72 tour 1.5 o1 o2; lg1 = droi pg0 pg00 ; lg2 = droi pg1 pg01 ; lg3 = droi pg2 pg02 ; lg4 = droi pg3 pg03 ; lg5 = droi pg4 pg04 ; lg6 = droi pg5 pg05 ; lg7 = droi pg6 pg06 ; lg8 = droi pg7 pg07 ; lg9 = droi pg8 pg08 ; lg10 = droi pg9 pg09 ; lg11 = droi pg11 pg011 ; lg12 = droi pg12 pg012 ; lg13 = droi pg31 pg031 ; lg14 = droi pg32 pg032 ; lg15 = droi pg51 pg051 ; lg16 = droi pg52 pg052 ; lg17 = droi pg71 pg071 ; lg18 = droi pg72 pg072 ; *connecteur 2 à 15 planp0 lg201 = lg1 tour 24 o1 o2 ; lg301 = lg1 tour 48 o1 o2 ; lg401 = lg1 tour 72 o1 o2 ; lg501 = lg1 tour 96 o1 o2 ; lg601 = lg1 tour 120 o1 o2 ; lg701 = lg1 tour 144 o1 o2 ; lg801 = lg1 tour 168 o1 o2 ; lg901 = lg1 tour 192 o1 o2 ; lg1001 = lg1 tour 216 o1 o2 ; lg1101 = lg1 tour 240 o1 o2 ; lg1201 = lg1 tour 264 o1 o2 ; lg1301 = lg1 tour 288 o1 o2 ; lg1401 = lg1 tour 312 o1 o2 ; lg1501 = lg1 tour 336 o1 o2 ; cg1 = c lg1 o1 lg201 ; cg12 = c lg201 o1 lg301 ; cg13 = c lg301 o1 lg401 ; cg14 = c lg401 o1 lg501 ; cg15 = c lg501 o1 lg601 ; cg16 = c lg601 o1 lg701 ; cg17 = c lg701 o1 lg801 ; cg18 = c lg801 o1 lg901 ; cg19 = c lg901 o1 lg1001 ; cg110 = c lg1001 o1 lg1101 ; cg111 = c lg1101 o1 lg1201 ; cg112 = c lg1201 o1 lg1301 ; cg113 = c lg1301 o1 lg1401 ; cg114 = c lg1401 o1 lg1501 ; cg115 = c lg1501 o1 lg1 ; cegg1 = (cg1 et cg12 et cg13 et cg14 et cg15 et cg16 et cg17 et cg18 et cg19 et cg110 et cg111 et cg112 et cg113 et cg114 et cg115); *connecteur 2 à 15 planp1 lg202 = lg2 tour 24 o1 o2 ; lg302 = lg2 tour 48 o1 o2 ; lg402 = lg2 tour 72 o1 o2 ; lg502 = lg2 tour 96 o1 o2 ; lg602 = lg2 tour 120 o1 o2 ; lg702 = lg2 tour 144 o1 o2 ; lg802 = lg2 tour 168 o1 o2 ; lg902 = lg2 tour 192 o1 o2 ; lg1002 = lg2 tour 216 o1 o2 ; lg1102 = lg2 tour 240 o1 o2 ; lg1202 = lg2 tour 264 o1 o2 ; lg1302 = lg2 tour 288 o1 o2 ; lg1402 = lg2 tour 312 o1 o2 ; lg1502 = lg2 tour 336 o1 o2 ; o3 = 0. pl 0.; cg2 = c lg2 o3 lg202 ; cg22 = c lg202 o3 lg302 ; cg23 = c lg302 o3 lg402 ; cg24 = c lg402 o3 lg502 ; cg25 = c lg502 o3 lg602 ; cg26 = c lg602 o3 lg702 ; cg27 = c lg702 o3 lg802 ; cg28 = c lg802 o3 lg902 ; cg29 = c lg902 o3 lg1002 ; cg210 = c lg1002 o3 lg1102 ; cg211 = c lg1102 o3 lg1202 ; cg212 = c lg1202 o3 lg1302 ; cg213 = c lg1302 o3 lg1402 ; cg214 = c lg1402 o3 lg1502 ; cg215 = c lg1502 o3 lg2 ; cegg2 = (cg2 et cg22 et cg23 et cg24 et cg25 et cg26 et cg27 et cg28 et cg29 et cg210 et cg211 et cg212 et cg213 et cg214 et cg215); *connecteur 2 à 15 planp2 lg203 = lg3 tour 24 o1 o2 ; lg303 = lg3 tour 48 o1 o2 ; lg403 = lg3 tour 72 o1 o2 ; lg503 = lg3 tour 96 o1 o2 ; lg603 = lg3 tour 120 o1 o2 ; lg703 = lg3 tour 144 o1 o2 ; lg803 = lg3 tour 168 o1 o2 ; lg903 = lg3 tour 192 o1 o2 ; lg1003 = lg3 tour 216 o1 o2 ; lg1103 = lg3 tour 240 o1 o2 ; lg1203 = lg3 tour 264 o1 o2 ; lg1303 = lg3 tour 288 o1 o2 ; lg1403 = lg3 tour 312 o1 o2 ; lg1503 = lg3 tour 336 o1 o2 ; pl4 = pl + 2 * pro + ps; o4 = 0. pl4 0.; cg3 = c lg3 o4 lg203 ; cg32 = c lg203 o4 lg303 ; cg33 = c lg303 o4 lg403 ; cg34 = c lg403 o4 lg503 ; cg35 = c lg503 o4 lg603 ; cg36 = c lg603 o4 lg703 ; cg37 = c lg703 o4 lg803 ; cg38 = c lg803 o4 lg903 ; cg39 = c lg903 o4 lg1003 ; cg310 = c lg1003 o4 lg1103 ; cg311 = c lg1103 o4 lg1203 ; cg312 = c lg1203 o4 lg1303 ; cg313 = c lg1303 o4 lg1403 ; cg314 = c lg1403 o4 lg1503 ; cg315 = c lg1503 o4 lg3 ; cegg3 = (cg3 et cg32 et cg33 et cg34 et cg35 et cg36 et cg37 et cg38 et cg39 et cg310 et cg311 et cg312 et cg313 et cg314 et cg315); *connecteur 2 à 15 planp3 lg204 = lg4 tour 24 o1 o2 ; lg304 = lg4 tour 48 o1 o2 ; lg404 = lg4 tour 72 o1 o2 ; lg504 = lg4 tour 96 o1 o2 ; lg604 = lg4 tour 120 o1 o2 ; lg704 = lg4 tour 144 o1 o2 ; lg804 = lg4 tour 168 o1 o2 ; lg904 = lg4 tour 192 o1 o2 ; lg1004 = lg4 tour 216 o1 o2 ; lg1104 = lg4 tour 240 o1 o2 ; lg1204 = lg4 tour 264 o1 o2 ; lg1304 = lg4 tour 288 o1 o2 ; lg1404 = lg4 tour 312 o1 o2 ; lg1504 = lg4 tour 336 o1 o2 ; pl5 = pl + 2*pro + 2* ps ; o5 = 0. pl5 0.; cg4 = c lg4 o5 lg204 ; cg42 = c lg204 o5 lg304 ; cg43 = c lg304 o5 lg404 ; cg44 = c lg404 o5 lg504 ; cg45 = c lg504 o5 lg604 ; cg46 = c lg604 o5 lg704 ; cg47 = c lg704 o5 lg804 ; cg48 = c lg804 o5 lg904 ; cg49 = c lg904 o5 lg1004 ; cg410 = c lg1004 o5 lg1104 ; cg411 = c lg1104 o5 lg1204 ; cg412 = c lg1204 o5 lg1304 ; cg413 = c lg1304 o5 lg1404 ; cg414 = c lg1404 o5 lg1504 ; cg415 = c lg1504 o5 lg4 ; cegg4 = (cg4 et cg42 et cg43 et cg44 et cg45 et cg46 et cg47 et cg48 et cg49 et cg410 et cg411 et cg412 et cg413 et cg414 et cg415); *connecteur 2 à 15 planp4 lg205 = lg5 tour 24 o1 o2 ; lg305 = lg5 tour 48 o1 o2 ; lg405 = lg5 tour 72 o1 o2 ; lg505 = lg5 tour 96 o1 o2 ; lg605 = lg5 tour 120 o1 o2 ; lg705 = lg5 tour 144 o1 o2 ; lg805 = lg5 tour 168 o1 o2 ; lg905 = lg5 tour 192 o1 o2 ; lg1005 = lg5 tour 216 o1 o2 ; lg1105 = lg5 tour 240 o1 o2 ; lg1205 = lg5 tour 264 o1 o2 ; lg1305 = lg5 tour 288 o1 o2 ; lg1405 = lg5 tour 312 o1 o2 ; lg1505 = lg5 tour 336 o1 o2 ; pl6 = pl + 4 * pro + 3* ps ; o6 = 0. pl6 0.; cg5 = c lg5 o6 lg205 ; cg52 = c lg205 o6 lg305 ; cg53 = c lg305 o6 lg405 ; cg54 = c lg405 o6 lg505 ; cg55 = c lg505 o6 lg605 ; cg56 = c lg605 o6 lg705 ; cg57 = c lg705 o6 lg805 ; cg58 = c lg805 o6 lg905 ; cg59 = c lg905 o6 lg1005 ; cg510 = c lg1005 o6 lg1105 ; cg511 = c lg1105 o6 lg1205 ; cg512 = c lg1205 o6 lg1305 ; cg513 = c lg1305 o6 lg1405 ; cg514 = c lg1405 o6 lg1505 ; cg515 = c lg1505 o6 lg5 ; cegg5 = (cg5 et cg52 et cg53 et cg54 et cg55 et cg56 et cg57 et cg58 et cg59 et cg510 et cg511 et cg512 et cg513 et cg514 et cg515); *connecteur 2 à 15 planp5 lg206 = lg6 tour 24 o1 o2 ; lg306 = lg6 tour 48 o1 o2 ; lg406 = lg6 tour 72 o1 o2 ; lg506 = lg6 tour 96 o1 o2 ; lg606 = lg6 tour 120 o1 o2 ; lg706 = lg6 tour 144 o1 o2 ; lg806 = lg6 tour 168 o1 o2 ; lg906 = lg6 tour 192 o1 o2 ; lg1006 = lg6 tour 216 o1 o2 ; lg1106 = lg6 tour 240 o1 o2 ; lg1206 = lg6 tour 264 o1 o2 ; lg1306 = lg6 tour 288 o1 o2 ; lg1406 = lg6 tour 312 o1 o2 ; lg1506 = lg6 tour 336 o1 o2 ; pl7 = pl + 4* pro + 4* ps ; o7 = 0. pl7 0.; cg6 = c lg6 o7 lg206 ; cg62 = c lg206 o7 lg306 ; cg63 = c lg306 o7 lg406 ; cg64 = c lg406 o7 lg506 ; cg65 = c lg506 o7 lg606 ; cg66 = c lg606 o7 lg706 ; cg67 = c lg706 o7 lg806 ; cg68 = c lg806 o7 lg906 ; cg69 = c lg906 o7 lg1006 ; cg610 = c lg1006 o7 lg1106 ; cg611 = c lg1106 o7 lg1206 ; cg612 = c lg1206 o7 lg1306 ; cg613 = c lg1306 o7 lg1406 ; cg614 = c lg1406 o7 lg1506 ; cg615 = c lg1506 o7 lg6 ; cegg6 = (cg6 et cg62 et cg63 et cg64 et cg65 et cg66 et cg67 et cg68 et cg69 et cg610 et cg611 et cg612 et cg613 et cg614 et cg615); *connecteur 2 à 15 planp6 lg207 = lg7 tour 24 o1 o2 ; lg307 = lg7 tour 48 o1 o2 ; lg407 = lg7 tour 72 o1 o2 ; lg507 = lg7 tour 96 o1 o2 ; lg607 = lg7 tour 120 o1 o2 ; lg707 = lg7 tour 144 o1 o2 ; lg807 = lg7 tour 168 o1 o2 ; lg907 = lg7 tour 192 o1 o2 ; lg1007 = lg7 tour 216 o1 o2 ; lg1107 = lg7 tour 240 o1 o2 ; lg1207 = lg7 tour 264 o1 o2 ; lg1307 = lg7 tour 288 o1 o2 ; lg1407 = lg7 tour 312 o1 o2 ; lg1507 = lg7 tour 336 o1 o2 ; pl8 = pl + (6* pro) + (5* ps) ; o8 = 0. pl8 0.; cg7 = c lg7 o8 lg207 ; cg72 = c lg207 o8 lg307 ; cg73 = c lg307 o8 lg407 ; cg74 = c lg407 o8 lg507 ; cg75 = c lg507 o8 lg607 ; cg76 = c lg607 o8 lg707 ; cg77 = c lg707 o8 lg807 ; cg78 = c lg807 o8 lg907 ; cg79 = c lg907 o8 lg1007 ; cg710 = c lg1007 o8 lg1107 ; cg711 = c lg1107 o8 lg1207 ; cg712 = c lg1207 o8 lg1307 ; cg713 = c lg1307 o8 lg1407 ; cg714 = c lg1407 o8 lg1507 ; cg715 = c lg1507 o8 lg7 ; cegg7 = (cg7 et cg72 et cg73 et cg74 et cg75 et cg76 et cg77 et cg78 et cg79 et cg710 et cg711 et cg712 et cg713 et cg714 et cg715); *connecteur 2 à 15 planp7 lg208 = lg8 tour 24 o1 o2 ; lg308 = lg8 tour 48 o1 o2 ; lg408 = lg8 tour 72 o1 o2 ; lg508 = lg8 tour 96 o1 o2 ; lg608 = lg8 tour 120 o1 o2 ; lg708 = lg8 tour 144 o1 o2 ; lg808 = lg8 tour 168 o1 o2 ; lg908 = lg8 tour 192 o1 o2 ; lg1008 = lg8 tour 216 o1 o2 ; lg1108 = lg8 tour 240 o1 o2 ; lg1208 = lg8 tour 264 o1 o2 ; lg1308 = lg8 tour 288 o1 o2 ; lg1408 = lg8 tour 312 o1 o2 ; lg1508 = lg8 tour 336 o1 o2 ; pl9 = pl + (6* pro) + (6 * ps) ; o9 = 0. pl9 0.; cg8 = c lg8 o9 lg208 ; cg82 = c lg208 o9 lg308 ; cg83 = c lg308 o9 lg408 ; cg84 = c lg408 o9 lg508 ; cg85 = c lg508 o9 lg608 ; cg86 = c lg608 o9 lg708 ; cg87 = c lg708 o9 lg808 ; cg88 = c lg808 o9 lg908 ; cg89 = c lg908 o9 lg1008 ; cg810 = c lg1008 o9 lg1108 ; cg811 = c lg1108 o9 lg1208 ; cg812 = c lg1208 o9 lg1308 ; cg813 = c lg1308 o9 lg1408 ; cg814 = c lg1408 o9 lg1508 ; cg815 = c lg1508 o9 lg8 ; cegg8 = (cg8 et cg82 et cg83 et cg84 et cg85 et cg86 et cg87 et cg88 et cg89 et cg810 et cg811 et cg812 et cg813 et cg814 et cg815); *connecteur 2 à 15 planp8 lg209 = lg9 tour 24 o1 o2 ; lg309 = lg9 tour 48 o1 o2 ; lg409 = lg9 tour 72 o1 o2 ; lg509 = lg9 tour 96 o1 o2 ; lg609 = lg9 tour 120 o1 o2 ; lg709 = lg9 tour 144 o1 o2 ; lg809 = lg9 tour 168 o1 o2 ; lg909 = lg9 tour 192 o1 o2 ; lg1009 = lg9 tour 216 o1 o2 ; lg1109 = lg9 tour 240 o1 o2 ; lg1209 = lg9 tour 264 o1 o2 ; lg1309 = lg9 tour 288 o1 o2 ; lg1409 = lg9 tour 312 o1 o2 ; lg1509 = lg9 tour 336 o1 o2 ; pl10 = pl + (8* pro) + (7* ps) ; o10 = 0. pl10 0.; cg9 = c lg9 o10 lg209 ; cg92 = c lg209 o10 lg309 ; cg93 = c lg309 o10 lg409 ; cg94 = c lg409 o10 lg509 ; cg95 = c lg509 o10 lg609 ; cg96 = c lg609 o10 lg709 ; cg97 = c lg709 o10 lg809 ; cg98 = c lg809 o10 lg909 ; cg99 = c lg909 o10 lg1009 ; cg910 = c lg1009 o10 lg1109 ; cg911 = c lg1109 o10 lg1209 ; cg912 = c lg1209 o10 lg1309 ; cg913 = c lg1309 o10 lg1409 ; cg914 = c lg1409 o10 lg1509 ; cg915 = c lg1509 o10 lg9 ; cegg9 = (cg9 et cg92 et cg93 et cg94 et cg95 et cg96 et cg97 et cg98 et cg99 et cg910 et cg911 et cg912 et cg913 et cg914 et cg915); *assemblage galerie lij = ((2* pro) + ps); ceg1 = cegg1 trans (0. pl 0.) ; ceg2 = cegg2 trans (0. lij 0.) ; ceg3 = cegg3 trans (0. ps 0.) ; ceg4 = cegg4 trans (0. lij 0.) ; ceg5 = cegg5 trans (0. ps 0.) ; ceg6 = cegg6 trans (0. lij 0.) ; ceg7 = cegg7 trans (0. ps 0.) ; ceg8 = cegg8 trans (0. lij 0.) ; ceg9 = cegg9 trans (0. pl 0.) ; gal = (ceg1 et ceg2 et ceg3 et ceg4 et ceg5 et ceg6 et ceg7 et ceg8 et ceg9) ; elim gal ; titr'galerie' ; trac gal ; titr'galerie, tole et sinusoides' ; trac (gal et tole et sinutot) ; *CL galerie* ************ kg = bloq gal ux uy uz rx ry rz ; *CL, relations de deplacements galerie sinusoide* ************************************************* sinup = (sinu et sinu2 et sinu3 et sinu4 et sinu13 et sinu14 et sinu15) ; sinupg = (sinu13 et sinu14 et sinu15) ; sinupd= (sinu et sinu2 et sinu3 et sinu4) ; sinum = (sinu5 et sinu6 et sinu7 et sinu8 et sinu9 et sinu10 et sinu11 et sinu12) ; sinumg = (sinu9 et sinu10 et sinu11 et sinu12) ; sinumd = (sinu5 et sinu6 et sinu7 et sinu8) ; kx1d = rela 'MAXI' 'DEPL' 'DIRECTION' (0.010 0. 0.) sinupd ; kx1g = rela 'MAXI' 'DEPL' 'DIRECTION' (-0.010 0. 0.) sinupg ; kz1 = rela 'MAXI' 'DEPL' 'DIRECTION' (0. 0. 0.010) sinup ; kx2d = rela 'MAXI' 'DEPL' 'DIRECTION' (0.010 0. 0.) sinumd ; kx2g = rela 'MAXI' 'DEPL' 'DIRECTION' (-0.010 0. 0.) sinumg ; kz2 = rela 'MAXI' 'DEPL' 'DIRECTION' (0. 0. -0.010) sinum ; *Relation de déplacements tole* ******************************* *surface de tole positive droite tot1 = ( c1 et c12 et c13 et c14 ) ; totr1 = tot1 trans (0. pl 0.) ; tot2 = ( c2 et c22 et c23 et c24 ) ; totr2 = tot2 trans (0. lij 0.) ; tot3 = ( c3 et c32 et c33 et c34 ) ; totr3 = tot3 trans (0. ps 0.) ; tot4 = ( c4 et c42 et c43 et c44 ) ; totr4 = tot4 trans (0. lij 0.) ; tot5 = (c5 et c52 et c53 et c54 ) ; totr5 = tot5 trans (0. ps 0.) ; tot6 = ( c6 et c62 et c63 et c64 ) ; totr6 = tot6 trans (0. lij 0.) ; tot7 = ( c7 et c72 et c73 et c74 ) ; totr7 = tot7 trans (0. ps 0.) ; tot8 = ( c8 et c82 et c83 et c84 ) ; totr8 = tot8 trans (0. lij 0.) ; tot9 = ( c9 et c92 et c93 et c94 ) ; totr9 = tot9 trans (0. pl 0.) ; tott = (totr1 et totr2 et totr3 et totr4 et totr5 et totr6 et totr7 et totr8 et totr9); elim (tott et tole) ; tolpd = tole elem appuy strict tott ; *surface de tole negative droite totmd1 = ( c15 et c16 et c17 et c18 ) ; totrmd1 = totmd1 trans (0. pl 0.) ; totmd2 = ( c25 et c26 et c27 et c28 ) ; totrmd2 = totmd2 trans (0. lij 0.) ; totmd3 = ( c35 et c36 et c37 et c38) ; totrmd3 = totmd3 trans (0. ps 0.) ; totmd4 = ( c45 et c46 et c47 et c48 ) ; totrmd4 = totmd4 trans (0. lij 0.) ; totmd5 = (c55 et c56 et c57 et c58 ) ; totrmd5 = totmd5 trans (0. ps 0.) ; totmd6 = ( c65 et c66 et c67 et c68 ) ; totrmd6 = totmd6 trans (0. lij 0.) ; totmd7 = ( c75 et c76 et c77 et c78 ) ; totrmd7 = totmd7 trans (0. ps 0.) ; totmd8 = ( c85 et c86 et c87 et c88 ) ; totrmd8 = totmd8 trans (0. lij 0.) ; totmd9 = ( c95 et c96 et c97 et c98 ) ; totrmd9 = totmd9 trans (0. pl 0.) ; tottmd = (totrmd1 et totrmd2 et totrmd3 et totrmd4 et totrmd5 et totrmd6 et totrmd7 et totrmd8 et totrmd9); elim (tottmd et tole) ; tolmd = tole elem appuy strict tottmd ; *surface de tole negative gauche totmg1 = ( c19 et c110 et c111 et c112 ) ; totrmg1 = totmg1 trans (0. pl 0.) ; totmg2 = ( c29 et c210 et c211 et c212 ) ; totrmg2 = totmg2 trans (0. lij 0.) ; totmg3 = ( c39 et c310 et c311 et c312) ; totrmg3 = totmg3 trans (0. ps 0.) ; totmg4 = ( c49 et c410 et c411 et c412 ) ; totrmg4 = totmg4 trans (0. lij 0.) ; totmg5 = (c59 et c510 et c511 et c512 ) ; totrmg5 = totmg5 trans (0. ps 0.) ; totmg6 = ( c69 et c610 et c611 et c612 ) ; totrmg6 = totmg6 trans (0. lij 0.) ; totmg7 = ( c79 et c710 et c711 et c712 ) ; totrmg7 = totmg7 trans (0. ps 0.) ; totmg8 = ( c89 et c810 et c811 et c812 ) ; totrmg8 = totmg8 trans (0. lij 0.) ; totmg9 = ( c99 et c910 et c911 et c912 ) ; totrmg9 = totmg9 trans (0. pl 0.) ; tottmg = (totrmg1 et totrmg2 et totrmg3 et totrmg4 et totrmg5 et totrmg6 et totrmg7 et totrmg8 et totrmg9); elim (tottmg et tole) ; tolmg = tole elem appuy strict tottmg ; *surface de tole positive gauche totpg1 = ( c113 et c114 et c115 ) ; totrpg1 = totpg1 trans (0. pl 0.) ; totpg2 = ( c213 et c214 et c215 ) ; totrpg2 = totpg2 trans (0. lij 0.) ; totpg3 = ( c313 et c314 et c315 ) ; totrpg3 = totpg3 trans (0. ps 0.) ; totpg4 = ( c413 et c414 et c415) ; totrpg4 = totmg4 trans (0. lij 0.) ; totpg5 = (c513 et c514 et c515) ; totrpg5 = totpg5 trans (0. ps 0.) ; totpg6 = ( c613 et c614 et c615 ) ; totrpg6 = totpg6 trans (0. lij 0.) ; totpg7 = ( c713 et c714 et c715 ) ; totrpg7 = totpg7 trans (0. ps 0.) ; totpg8 = ( c813 et c814 et c815 ) ; totrpg8 = totpg8 trans (0. lij 0.) ; totpg9 = ( c913 et c914 et c915 ) ; totrpg9 = totpg9 trans (0. pl 0.) ; tottpg = (totrpg1 et totrpg2 et totrpg3 et totrpg4 et totrpg5 et totrpg6 et totrpg7 et totrpg8 et totrpg9); elim (tottpg et tole) ; tolpg = tole elem appuy strict tottpg ; *surface de tole dessous ad12 = c18 ; atd12 = ad12 trans (0. pl 0.) ; ad22 = c28 ; atd22 = ad22 trans (0. lij 0.) ; ad32 = c38 ; atd32 = ad32 trans (0. ps 0.) ; ad42 = c48 ; atd42 = ad42 trans (0. lij 0.) ; ad52 = c58 ; atd52 = ad52 trans (0. ps 0.) ; ad62 = c68 ; atd62 = ad62 trans (0. lij 0.) ; ad72 = c78 ; atd72 = ad72 trans (0. ps 0.) ; ad82 = c88 ; atd82 = ad82 trans (0. lij 0.) ; ad92 = c98 ; atd92 = ad92 trans (0. pl 0.) ; attt =(atd12 et atd22 et atd32 et atd42 et atd52 et atd62 et atd72 et atd82 et atd92); riu = (atd12 et atd22 et atd32 et atd42 et atd52 et atd62 et atd72 et atd82 et atd92 ); elim 0.01 (riu et tole) ; riu2 = tole elem appuy largement riu ; titr 'Surface de contact dessous'; trac riu2; riu2 = riu2 coul rouge ; kt1d = rela 'MAXI' 'DEPL' 'DIRECTION' (0.050 0. 0.) tolpd ; kt1g = rela 'MAXI' 'DEPL' 'DIRECTION' (-0.050 0. 0.) tolpg ; kt2d = rela 'MAXI' 'DEPL' 'DIRECTION' (0.050 0. 0.) tolmd ; kt2g = rela 'MAXI' 'DEPL' 'DIRECTION' (-0.050 0. 0.) tolmg ; kt1z = rela 'MAXI' 'DEPL' 'DIRECTION' (0.050 0. 0.) tolpd ; kt2z = rela 'MAXI' 'DEPL' 'DIRECTION' (0.050 0. 0.) tolpg ; kt3z = rela 'MAXI' 'DEPL' 'DIRECTION' (-0.050 0. 0.) tolmd ; kt4z = rela 'MAXI' 'DEPL' 'DIRECTION' (-0.05 0. 0.) tolmg ; ktdz = rela 'MAXI' 'DEPL' 'DIRECTION' (-0.05 0. 0.) riu2 ; ************************************************************************ * 3/ Caracteristiques matériaux et charges * ************************************************************************ *caractéristiques* ****************** mod1 = mode (tole et sinutot) mecanique elastique dkt ; mate1 = mate mod1 young you1 nu nu1 rho rho1 epai ep; *poids propre* ************** mass1 = mass mod1 mate1 ; pop1 = manu chpo (tole et sinutot) 1 uz -9.81 nature diffus ; fo1 = mass1 * pop1 ; vectfo1 = vect fo1 fx fy fz jaune ; *liquide beton* *************** *cote droit* ************ att1 = ( c1 et c12 et c13 et c14 et c15 et c16 et c17 ) ; atr1 = att1 trans (0. pl 0.) ; att2 = ( c2 et c22 et c23 et c24 et c25 et c26 et c27 ) ; atr2 = att2 trans (0. lij 0.) ; att3 = ( c3 et c32 et c33 et c34 et c35 et c36 et c37 ) ; atr3 = att3 trans (0. ps 0.) ; att4 = ( c4 et c42 et c43 et c44 et c45 et c46 et c47 ) ; atr4 = att4 trans (0. lij 0.) ; att5 = (c5 et c52 et c53 et c54 et c55 et c56 et c57 ) ; atr5 = att5 trans (0. ps 0.) ; att6 = ( c6 et c62 et c63 et c64 et c65 et c66 et c67 ) ; atr6 = att6 trans (0. lij 0.) ; att7 = ( c7 et c72 et c73 et c74 et c75 et c76 et c77 ) ; atr7 = att7 trans (0. ps 0.) ; att8 = ( c8 et c82 et c83 et c84 et c85 et c86 et c87 ) ; atr8 = att8 trans (0. lij 0.) ; att9 = ( c9 et c92 et c93 et c94 et c95 et c96 et c97 ) ; atr9 = att9 trans (0. pl 0.) ; atrt = (atr1 et atr2 et atr3 et atr4 et atr5 et atr6 et atr7 et atr8 et atr9); elim 0.01 (atrt et tole) ; resu1 = tole elem appuy largement atrt ; *zto1 = (coor 3 resu1) + (ray *(manu chpo tole 1 scal 1.)); *cote gauche* ************* att12 = ( c19 et c110 et c111 et c112 et c113 et c114 et c115 ) ; atr12 = att12 trans (0. pl 0.) ; att22 = ( c29 et c210 et c211 et c212 et c213 et c214 et c215) ; atr22 = att22 trans (0. lij 0.) ; att32 = ( c39 et c310 et c311 et c312 et c313 et c314 et c315) ; atr32 = att32 trans (0. ps 0.) ; att42 = ( c49 et c410 et c411 et c412 et c413 et c414 et c415) ; atr42 = att42 trans (0. lij 0.) ; att52 = ( c59 et c510 et c511 et c512 et c513 et c514 et c515) ; atr52 = att52 trans (0. ps 0.) ; att62 = (c69 et c610 et c611 et c612 et c613 et c614 et c615) ; atr62 = att62 trans (0. lij 0.) ; att72 = ( c79 et c710 et c711 et c712 et c713 et c714 et c715) ; atr72 = att72 trans (0. ps 0.) ; att82 = (c89 et c810 et c811 et c812 et c813 et c814 et c815) ; atr82 = att82 trans (0. lij 0.) ; att92 = ( c99 et c910 et c911 et c912 et c913 et c914 et c915) ; atr92 = att92 trans (0. pl 0.) ; atrt2 = (atr12 et atr22 et atr32 et atr42 et atr52 et atr62 et atr72 et atr82 et atr92); elim 0.01 (atrt2 et tole) ; resu2 = tole elem appuy largement atrt2 ; *zto2 = (coor 3 resu2) + (ray *(manu chpo tole 1 scal 1.)); *dessous* ********* atrd2 = (atd12 et atd22 et atd32 et atd42 et atd52 et atd62 et atd72 et atd82 et atd92 et atr12 et atr22 et atr32 et atr42 et atr52 et atr62 et atr72 et atr82 et atr92 et atr1 et atr2 et atr3 et atr4 et atr5 et atr6 et atr7 et atr8 et atr9); elim 0.01 (atrd2 et tole) ; resu3 = tole elem appuy largement atrd2 ; *calcul surface de contact droite* ********************************** mess 'Hauteur de betonnage en m'; mess ' Choix ?'; zv = lireflot -1. 1.; ppress = p * (zv + 1.) ; poin1 = sinu8 poin initial ; poin2 = sinu8 poin final ; zmind = coor 3 poin1 ; ensemp = sinu8 poin droi poin1 poin2 ; dmind = droi poin1 poin2 ; elim (dmind et tole) ; pmaxdd = sinu2 poin initial; pmaxfd = sinu2 poin final ; zmaxd = coor 3 pmaxdd ; ensemp2 = sinu2 poin droi pmaxdd pmaxfd ; dmaxd = droi pmaxdd pmaxfd ; elim (dmaxd et tole) ; zu1 = (zv < zmind) ; zu2 = (zv > zmaxd) ; si zu1 ; pk1 = l801 poin final ; hu0 = c pk1 o1 poin1 ; hu1 = trans hu0 (0. pl 0.) ; hu2 = trans l802 (0. lij 0.) ; hu3 = trans l803 (0. ps 0.) ; hu4 = trans l804 (0. lij 0.) ; hu5 = trans l805 (0. ps 0.) ; hu6 = trans l806 (0. lij 0.) ; hu7 = trans l807 (0. ps 0.) ; hu8 = trans l808 (0. lij 0.) ; hu9 = trans l809 (0. pl 0.) ; surfcald = (hu1 et hu2 et hu3 et hu4 et hu5 et hu6 et hu7 et hu8 et hu9) ; elim 0.01(surfcald et tole) ; res = surfcald elem appuy strict tole ; elim 0.01 (res et tole) ; titr 'Surface d etude cote droit' ; trac res ; res1 = orie res (1. 0. 0.) ; mod2 = mode (res1) mecanique elastique dkt ; mate2 = mate mod2 young you1 nu nu1 rho rho1 epai ep ; chpress = press coqu mod2 ppress (-1. 0. 0.) ; sinon ; si zu2 ; surfcald = resu1 ; res1 = surfcald elem appuy largement tole ; elim (res1 et tole); titr 'Surface d etude cote droit' ; trac res1 ; mod2 = mode (res1) mecanique elastique dkt ; mate2 = mate mod2 young you1 nu nu1 rho rho1 epai ep; res1 = orie res1 'DIRECTION' (1. 0. 0.) ; chpress = press coqu mod2 ppress (-1. 0. 0.) ; sinon ; ab = ((ray ** 2 ) - ((zv)** 2))** (1/ 2) ; pinter = ab 0. zv ; youyou = pl10 + pl ; pin = ab youyou zv ; pmax1 = resu1 poin proc pinter ; x1 = coor 1 pmax1 ; y1 = coor 2 pmax1 ; z1 = coor 3 pmax1 ; x2 = coor 1 poin1 ; y2 = coor 2 poin1 ; z2 = coor 3 poin1 ; t1 = x1 pl z1 ; t3 = x1 pl4 z1 ; t5 = x1 pl5 z1 ; t7 = x1 pl6 z1 ; t9 = x1 pl7 z1 ; t11 = x1 pl8 z1 ; t13 = x1 pl9 z1 ; t15 = x1 pl10 z1 ; t17 = x1 youyou z1 ; t2 = x2 pl z2 ; t4 = x2 pl4 z2 ; t6 = x2 pl5 z2 ; t8 = x2 pl6 z2 ; t10 = x2 pl7 z2 ; t12 = x2 pl8 z2 ; t14 = x2 pl9 z2 ; t16 = x2 pl10 z2 ; t18 = x2 youyou z2 ; tl1 = droi pmax1 t1 ; tl2 = droi t1 t3 ; tl3 = droi t3 t5 ; tl4 = droi t5 t7 ; tl5 = droi t7 t9 ; tl6 = droi t9 t11 ; tl7 = droi t11 t13 ; tl8 = droi t13 t15 ; tl9 = droi t15 t17 ; tl10 = droi poin1 t2 ; tl11 = droi t2 t4 ; tl12 = droi t4 t6 ; tl13 = droi t6 t8 ; tl14 = droi t8 t10 ; tl15 = droi t10 t12 ; tl16 = droi t12 t14 ; tl17 = droi t14 t16 ; tl18 = droi t16 t18 ; tc1 = c pmax1 o1 poin1 ; tc2 = c t1 o3 t2 ; tc3 = c t3 o4 t4 ; tc4 = c t5 o5 t6 ; tc5 = c t7 o6 t8 ; tc6 = c t9 o7 t10 ; tc7 = c t11 o8 t12 ; tc8 = c t13 o9 t14 ; tc9 = c t15 o10 t16 ; oo = 0. youyou 0. ; tc10 = c t17 oo t18 ; si1 = dall tc1 tl10 tc2 tl1 cyli o1 o2 ; si2 = dall tc2 tl11 tc3 tl2 cyli o1 o2 ; si3 = dall tc3 tl12 tc4 tl3 cyli o1 o2 ; si4 = dall tc4 tl13 tc5 tl4 cyli o1 o2 ; si5 = dall tc5 tl14 tc6 tl5 cyli o1 o2 ; si6 = dall tc6 tl15 tc7 tl6 cyli o1 o2 ; si7 = dall tc7 tl16 tc8 tl7 cyli o1 o2 ; si8 = dall tc8 tl17 tc9 tl8 cyli o1 o2 ; si9 = dall tc9 tl18 tc10 tl9 cyli o1 o2 ; marre = (si1 et si2 et si3 et si4 et si5 et si6 et si7 et si8 et si9) ; elim 0.01 (marre et tole) ; res = marre elem appuy largement tole ; res1 = orie res 'DIRECTION' (1. 0. 0.) ; mod2 = mode (res1) mecanique elastique dkt ; mate2 = mate mod2 young you1 nu nu1 rho rho1 epai ep; chpress = press coqu mod2 ppress (-1. 0. 0.) ; vcp = vect chpress fx fy fz jaune ; finsi ; finsi ; *calcul surface de contact gauche* ********************************** poin3 = sinu9 poin initial ; poin4 = sinu9 poin final ; zming = coor 3 poin3 ; ensemp = sinu9 poin droi poin3 poin4 ; dming = droi poin3 poin4 ; elim (dming et tole) ; pmaxdg = sinu15 poin initial; pmaxfg = sinu15 poin final ; zmaxg = coor 3 pmaxdg ; ensemp3 = sinu15 poin droi pmaxdg pmaxfg ; dmaxg = droi pmaxdg pmaxfg ; elim (dmaxg et tole) ; zu1 = (zv < zming) ; zu2 = (zv > zmaxg) ; si zu1 ; pkg2 = l901 poin final ; hug0 = c pkg2 o1 poin3 ; hug1 = trans hug0 (0. pl 0.) ; hug2 = trans l902 (0. lij 0.) ; hug3 = trans l903 (0. ps 0.) ; hug4 = trans l904 (0. lij 0.) ; hug5 = trans l905 (0. ps 0.) ; hug6 = trans l906 (0. lij 0.) ; hug7 = trans l907 (0. ps 0.) ; hug8 = trans l908 (0. lij 0.) ; hug9 = trans l909 (0. pl 0.) ; sog = (hug1 et hug2 et hug3 et hug4 et hug5 et hug6 et hug7 et hug8 et hug9) ; elim 0.01 (sog et tole) ; resg = sog elem appuy strict tole ; elim 0.01 (resg et tole); titr 'Surface d etude cote gauche' ; trac resg ; resg1 = orie resg (1. 0. 0.) ; mod3 = mode (resg1) mecanique elastique dkt ; mate3 = mate mod3 young you1 nu nu1 rho rho1 epai ep; chpressg = press coqu mod3 ppress (1. 0. 0.) ; sinon ; si zu2 ; surfcalg = resu1 ; res = surfcalg elem appuy strict tole ; elim 0.01 (res et tole); titr 'Surface d etude cote gauche' ; trac resg ; resg1 = orie resg (1. 0. 0.) ; mod3 = mode (resg1) mecanique elastique dkt ; mate3 = mate mod3 young you1 nu nu1 rho rho1 epai ep; chpressg = press coqu mod3 ppress (1. 0. 0.) ; sinon ; o3o = 0. 0. 15. ; pin = syme pmax1 'PLAN' o1 o2 o3o; pmaxg = resu2 poin proc pin ; x1 = coor 1 pmaxg ; y1 = coor 2 pmaxg ; z1 = coor 3 pmaxg ; x2 = coor 1 poin3 ; y2 = coor 2 poin3 ; z2 = coor 3 poin3 ; tg1 = x1 pl z1 ; tg3 = x1 pl4 z1 ; tg5 = x1 pl5 z1 ; tg7 = x1 pl6 z1 ; tg9 = x1 pl7 z1 ; tg11 = x1 pl8 z1 ; tg13 = x1 pl9 z1 ; tg15 = x1 pl10 z1 ; tg17 = x1 youyou z1 ; tg2 = x2 pl z2 ; tg4 = x2 pl4 z2 ; tg6 = x2 pl5 z2 ; tg8 = x2 pl6 z2 ; tg10 = x2 pl7 z2 ; tg12 = x2 pl8 z2 ; tg14 = x2 pl9 z2 ; tg16 = x2 pl10 z2 ; tg18 = x2 youyou z2 ; tgl1 = droi pmaxg tg1 ; tgl2 = droi tg1 tg3 ; tgl3 = droi tg3 tg5 ; tgl4 = droi tg5 tg7 ; tgl5 = droi tg7 tg9 ; tgl6 = droi tg9 tg11 ; tgl7 = droi tg11 tg13 ; tgl8 = droi tg13 tg15 ; tgl9 = droi tg15 tg17 ; tgl10 = droi poin3 tg2 ; tgl11 = droi tg2 tg4 ; tgl12 = droi tg4 tg6 ; tgl13 = droi tg6 tg8 ; tgl14 = droi tg8 tg10 ; tgl15 = droi tg10 tg12 ; tgl16 = droi tg12 tg14 ; tgl17 = droi tg14 tg16 ; tgl18 = droi tg16 tg18 ; tgc1 = c pmaxg o1 poin3 ; tgc2 = c tg1 o3 tg2 ; tgc3 = c tg3 o4 tg4 ; tgc4 = c tg5 o5 tg6 ; tgc5 = c tg7 o6 tg8 ; tgc6 = c tg9 o7 tg10 ; tgc7 = c tg11 o8 tg12 ; tgc8 = c tg13 o9 tg14 ; tgc9 = c tg15 o10 tg16 ; oo = 0. youyou 0. ; tgc10 = c tg17 oo tg18 ; sig1 = dall tgc1 tgl10 tgc2 tgl1 cyli o1 o2 ; sig2 = dall tgc2 tgl11 tgc3 tgl2 cyli o1 o2 ; sig3 = dall tgc3 tgl12 tgc4 tgl3 cyli o1 o2 ; sig4 = dall tgc4 tgl13 tgc5 tgl4 cyli o1 o2 ; sig5 = dall tgc5 tgl14 tgc6 tgl5 cyli o1 o2 ; sig6 = dall tgc6 tgl15 tgc7 tgl6 cyli o1 o2 ; sig7 = dall tgc7 tgl16 tgc8 tgl7 cyli o1 o2 ; sig8 = dall tgc8 tgl17 tgc9 tgl8 cyli o1 o2 ; sig9 = dall tgc9 tgl18 tgc10 tgl9 cyli o1 o2 ; marreg = (sig1 et sig2 et sig3 et sig4 et sig5 et sig6 et sig7 et sig8 et sig9) ; elim 0.01 (marreg et tole) ; resg = marreg elem appuy strict tole ; resg1 = orie resg (1. 0. 0.) ; mod3 = mode (resg1) mecanique elastique dkt ; mate3 = mate mod3 young you1 nu nu1 rho rho1 epai ep; chpressg = press coqu mod3 ppress (1. 0. 0.) ; vcpg = vect chpressg fx fy fz jaune ; finsi ; finsi ; *calcul surface de contact dessous* *********************************** *on se contente d appliquer la force sur riu2 *titr *'Surf réelle d application de la pres fct hauteur de betonnage cote' ; *trac (res1 et resg1) ; *resii = resg1 et res1 et attt ; *titr *'Surf réelle d application de la pres fct hauteur de betonnage dessous'; *trac (resii) ; resu41 = riu2 elem appuy strict resu3 ; trac (resu41 et resg1 et res1); resu42 = resg1 moins resu41 ; resu4 = resu42 moins res1 ; trac resu4; trac (resg1 et resu4 et res1); mod4 = mode (resu4) mecanique elastique dkt ; mate4 = mate mod4 young you1 nu nu1 rho rho1 epai ep; zu3 = (zv ega 0.) ; zu4 = (zv < 0.) ; zu5 = (zv > 0.) ; si zu3 ; pp = long * 3.14159 * p * ray * ray *0.5 * 0.25 ; sinon; si zu4; pm = cos (abs (zv)) ; pom = 90. - pm ; pp = p * long * 0.5 * (abs (zv)) * (pom / ray) * ((ray ** 2) - (zv ** 2))** (1/2 ) * 2 ; sinon; pm = cos (abs (zv)) ; pom = 90. - pm ; pp0 = 3.14159 * 0.25 * (ray ** 2) * 0.25 ; pp1 = p * 0.5 * zv * pm * ((ray ** 2) - (zv ** 2))** (1/2 ) * 2 ; pp = (pp0 + pp1) * long ; chprd = press coqu mod4 pp (0. 0. 1.) ; vcpd = vect chprd fx fy fz jaune ; finsi ; finsi ; ************************************************************************ * 4/ Traitements * ************************************************************************ *deplacements* ************** oeil1 = 0. 1000. 0. ; *poids propre* ************** k1 = rigi mod1 mate1 ; k = k1 et cls ; dep = resou (kx1d et kx1g et k et kz1 et kx2d et kx2g et kz2 et kt1d et kt1g et kt2d et kt2g et kt1z et kt2z et kt3z et kt4z et ktdz ) (fo1 ); def0 = defo (tole et sinutot ) dep 0. vectfo1 bleu ; def1 = defo (tole et sinutot ) dep rouge ; *liquide beton* *************** *droite dep2 = resou (kx1d et kx1g et k et kz1 et kx2d et kx2g et kz2 et kt1d et kt1g et kt2d et kt2g et kt1z et kt2z et kt3z et kt4z et ktdz ) (chpress) ; defo0 = defo (tole et sinutot ) dep2 0. vcp bleu ; defo1 = defo (tole et sinutot ) dep2 rouge ; *gauche dep3 = resou (k et kltg et klsg) (chpressg) ; defog0 = defo (tole et sinutot ) dep3 0. vcpg bleu ; defog1 = defo (tole et sinutot ) dep3 rouge ; *dessous dep4 = resou (kx1d et kx1g et k et kz1 et kx2d et kx2g et kz2 et kt1d et kt1g et kt2d et kt2g et kt1z et kt2z et kt3z et kt4z et ktdz ) (chprd) ; defd0 = defo (tole et sinutot ) dep4 0. vcpd bleu ; defd1 = defo (tole et sinutot ) dep4 rouge ; *somme des chargements* *********************** dept = resou (kx1d et kx1g et k et kz1 et kx2d et kx2g et kz2 et kt1d et kt1g et kt2d et kt2g et kt1z et kt2z et kt3z et kt4z et ktdz ) (fo1 et chpress et chpressg ) ; defot0 = defo (tole et sinutot) dept 0. vectfo1 bleu ; defot1 = defo (tole et sinutot) dept rouge ; *Contraintes* ************* chp1 = sigma mod1 mate1 dep ; chp2 = sigma mod1 mate1 dept ; *animation* *********** ani = anime 30 (tole et sinutot) (dep) 30. ; ani2 = anime 30 (tole et sinutot ) (dept) 30. ; ************************************************************************ * 5/ Résultats * ************************************************************************ *deformations* ************** titr 'deformee au cours du chargement poids propre ' ; trac (def0 et def1) ; trac oeil1 (def0 et def1) ; trac ani oscil ; titr 'deformee au cours du chargement béton droite ' ; trac (defo0 et defo1) ; trac oeil1 (defo0 et defo1) ; titr 'deformee au cours du chargement béton gauche ' ; trac (defog0 et defog1) ; trac oeil1 (defog0 et defog1) ; titr'deformee au cours du chargement beton desous' ; trac (defd0 et defd1) ; trac oeil1 (defd0 et defd1) ; titr ' deformee avec poids propre et chargement total béton' ; trac (defot0 et defot1) ; trac oeil1 (defot0 et defot1) ; trac ani2 oscil ; *igcast 'DEFO' 'RESU' vol dept ; *contraintes* ************* titr 'contrainte sous poids propre' ; trac mod1 mate1 (chp1) ; titr 'contrainte sous chargement béton et poids propre' ; trac mod1 mate1 (chp2) ; fin ;