********************************************** * Plaque de plâtre BA20 * * * ********************************************** * Longueur L = 500 mm * épaisseur plâtre e= 20 mm * Flexion 3 pts * courbe déplacement, contrainte ... *----------Donnée----------------- jeuini = -0.0005; Lo= 0.500; ep = 0.005; la = 0.01; nl = 100; nep = 4; nla = 6; TITR 'FLEXION 3 pts BA20'; OPTI DIME 3 MODE TRID ELEM CUB8; OPTI TRAC OPEN; *----------Maillage----------------- *surface 1 A1= 0. 0. 0.; A2= (Lo/2) 0. 0.; A3= Lo 0. 0.; A4= Lo ep 0.; A5 = (Lo/2) ep 0.; A6 = 0. ep 0.; L1 = DROI nl A1 A3; L2 = DROI nep A3 A4; L3 = DROI nl A6 A4; L4 = DROI nep A1 A6; L5 = DROI nep A2 A5; 2L = INVE L2; 3L = INVE L3; 4L = INVE L4; 5L = INVE L5; S1 = DALL L1 L2 L3 L4; S1 = S1 coul bleu; A7 = poin S1 proc ((Lo/4) 0. 0.); *surface 2 B1= 0. 0. la; B2= (Lo/2) 0. la; B3= Lo 0. la; B4= Lo ep la; B5 = (Lo/2) ep la; B6 = 0. ep la; L7 = DROI nl B1 B3; L8 = DROI nep B3 B4; L9 = DROI nl B6 B4; L10 = DROI nep B1 B6; L11 = DROI nep B2 B5; 8L = INVE L8; 9L = INVE L9; 10L = INVE L10; 11L = INVE L11; S2 = DALL L7 L8 L9 L10; S2 = S2 coul bleu; B7 = poin S2 proc ((Lo/4) 0. la); * volume de la poutre V1 = VOLU nla S1 S2; * surface où sera appliqué le blocage pc1 = poin v1 proc (0.245 0. 0.); pc2 = poin v1 proc (0.255 0. 0.); pc3 = poin v1 proc (0.255 0. la); pc4 = poin v1 proc (0.245 0. la); l1c = droi 1 pc1 pc2; l2c = droi nla pc2 pc3; l3c = droi 1 pc3 pc4; l4c = droi nla pc4 pc1; sc1= dall l1c l2c l3c l4c; elim sc1 v1 0.0001; * surfaces d'application du chargement thermique * FE = Face exposée p1 = poin S1 proc (0. 0. 0.); p2 = poin S1 proc (lo 0. 0.); p3 = poin S2 proc (lo 0. la); p4 = poin S2 proc (0. 0. la); Lt1 = DROI nl p1 p2; Lt2 = DROI nla p2 p3; Lt3 = DROI nl p3 p4; Lt4 = DROI nla p4 p1; FE = dall Lt1 Lt2 Lt3 Lt4; FE = FE coul roug; elim FE V1 0.0001; * FNE = Face non-exposée p1 = poin S1 proc (0. ep 0.); p2 = poin S1 proc (lo ep 0.); p3 = poin S2 proc (lo ep la); p4 = poin S2 proc (0. ep la); Lt1 = DROI nl p1 p2; Lt2 = DROI nla p2 p3; Lt3 = DROI nl p3 p4; Lt4 = DROI nla p4 p1; FNE = dall Lt1 Lt2 Lt3 Lt4; FNE = FNE coul vert; elim FNE V1 0.0001; *-----------modèle de comportement du matériaux platre------ M1 = mode V1 MECANIQUE ELASTIQUE ISOTROPE; M2 = mode V1 THERMIQUE ISOTROPE; *alpha = f(T) pt1 = PROG 20. 50. 100. 150. 200. 250. 300. 350. 400. 450. 500. 550. 600. 650. 700. 750. 800. 860. 900. 950. 1000. 1050. 1100. 1150. 1200.; pk1=PROG 0 (3.684e-4 / 50.) (9.984e-4 / 100.) (16.484e-4 / 150.) (23.184e-4 / 200.) (30.084e-4 / 250.) (37.184e-4 / 300.) (44.484e-4 / 350.) (51.984e-4 / 400.) (59.684e-4 / 450.) (67.584e-4 / 500.) (75.684e-4 / 550.) (83.984e-4 / 600.) (92.484e-4 / 650.) (101.184e-4 / 700.) (110e-4 / 750.) (110e-4 / 800.) (0.011 / 860.) (118e-4 / 900.) (128e-4 / 950.) (138e-4 / 1000.) (148e-4 / 1050.) (0.0158 / 1100.) (0.0168 / 1150.) (178e-4 / 1200.); eva1 = EVOL MANU 'T' pt1 'ALPH' pk1; * Young T2 = prog 20. 100. 200. 300. 400. 500. 600. 700. 800. 900. 1000. 1100.; pk2=PROG 2.06E11 2.06E11 1.85E11 1.65E11 1.44E11 1.23E11 6.38E10 2.68E10 1.85E10 1.39E10 9.26E9 4.63E9 ; eve1 = EVOL MANU 'T' T2 'YOUN' pk2 ; mat1 = mate M1 YOUN 2.82E9 nu 0.3 ALPH eva1 RHO 845. ; * conductivite thermique pt3=PROG 0. 20. 50. 100. 150. 200. 250. 300. 350. 400. 450. 500. 550. 600. 650. 700. 750. 800. 1200.; pk2=PROG 53.3 53.3 52.3 50.7 49 47.3 45.7 44 42.3 40.7 39 37.3 35.7 34 32.4 30.7 29 27.3 27.3; evk2 = EVOL MANU T pt3 K pk2 ; * chaleur specifique pt4 = PROG 0. 20. 50. 100. 150. 200. 250. 300. 350. 400. 450. 500. 550. 600. 650. 700. 725. 735. 750. 785. 800. 850. 900. 1200. ; pc2 = PROG 439.8 439.8 459.7 487.6 510.4 529.8 547.3 564.7 583.7 605.9 632.9 666.5 708.3 760.2 813.8 1008.2 1666.2 5000 1482.9 875 803.3 694.7 650. 650. ; evc2 = EVOL MANU T pt4 'C' pc2 ; mat2 = mate M2 'K' evk2 'C' evc2 'RHO' 845. ; * mode de convection M3= MODE FE CONVECTION CONS CON1; M4= MODE FNE CONVECTION CONS CON3; modt= M1 ET M2 ET M3 ET M4; M5 = MODE FE RAYONNEMENT CONS CON2; M6 = MODE FNE RAYONNEMENT CONS CON4; modtr= M5 ET M6; **coefficients d'échange** **selon Eurocode 1** *convection matc1 = MATE M3 'H' 25.; matc2 = MATE M4 'H' 4.; *rayonnement matr1=MATE M5 'EMIS' 0.9; matr2=MATE M6 'EMIS' 0.9; matot=matc1 et matc2 et matr1 et matr2 et mat1 et mat2; *-----------------conditions limites mécanique--------------- lc1 = DROI nla A1 B1; lc2 = DROI nla A6 B6; lc3 = DROI nla A3 B3; lc4 = DROI nla A4 B4; 2lc = INVE lc2; 4lc = INVE lc4; Scl1 = dall lc1 L10 2lc 4L; Scl2 = dall lc3 L8 4lc 2L; elim Scl1 V1 0.0001; elim Scl2 V1 0.0001; CL1 = BLOQ DEPL lc1; CL2 = BLOQ Uy Uz lc3; CL = CL1 et CL2; *----------------Chargement thermique------------- * Face exposée LI1 = prog 0. 350. 600.; LI2 = prog 0. 1500. 300.; ev2 = evol manu 'temps' LI1 'T' LI2; chan1=manu chpo FE 1 T 1. 'NATURE' 'DISCRET'; chart1=CHAR TECO ev2 chan1; chart2=CHAR TERA ev2 chan1; *face non exposée tnexp=prog 0. 450. 900.; tempnexp=prog 20. 20. 20.; nexp=evol manu 't' tnexp 'T' tempnexp; chan2=manu chpo FNE 1 T 1. 'NATURE' 'DISCRET'; chart3=CHAR TECO nexp chan2; chart4=CHAR TERA nexp chan2; chart = chart1 et chart2 et chart3 et chart4; * dess ev2 titre 'flot= P.F(t)'; *------------------------- calcucl pas à pas----------------- TAB1 = table; TAB1.'MODELE'=modt; TAB1.'CARACTERISTIQUES'=matot; TAB1.'BLOCAGES_MECANIQUES'= CL; TAB1.'CHARGEMENT'=chart; TAB1.procedure_perso1 = vrai; TAB1.RAYONNEMENT = table ; TAB1.RAYONNEMENT . 1 = table ; TAB1.RAYONNEMENT . 1 . TYPE = 'INFINI' ; TAB1.RAYONNEMENT . 1 . MODELE = modtr ; *TAB1.'PROCEDURE_THERMIQUE'=NONLINEAIRE; TAB1.'PROCEDURE_THERMIQUE'=DUPONT; TAB1.'CTE_STEFAN_BOLTZMANN'=5.67E-8; TAB1.'TEMPERATURES'=table; TAB1.'TEMPERATURES' . 0 =manu chpo V1 1 T 20.; TAB1.'CELCIUS'= vrai; TAB1.'TEMPS_CALCULES'=prog 0 pas 0.5 600.; * *------------ procédure utilisateur perso1 -------------------- * debproc perso1 TAB1*table; WTAB = TAB1.'WTABLE'; dep = TAB1.'ESTIMATION'.'DEPLACEMENTS'; i = 1; j = 0.; k = 0.; delta = 0.0025; nb3 = 100; lbloq = prog; repe BO3 nb3; xj = (Lo/4) + j; pcontrol = Poin S1 proc (xj 0. (la / 2.)); ypc = extr dep UY pcontrol; SI (ypc '<' jeuini); mess 'Doit être bloqué'; * lbloq = inse lbloq i xj; mess ' *************************'; mess ' mise a jour du au contact entre le seuil et le pt pc:'; mess ' *************************'; * création du nouveau blocage * CL3 = BLOQ 'MINI' UY pcontrol; CL3 = BLOQ 'MINI' UY sc1; * création du nouveau chargement tim = DEPI CL3 (jeuini); EVF1 = EVOL MANU 'Temps' LI1 'x' (prog 1. 1. 1.); CHA2 = CHAR 'DIMP' tim EVF1; * actualisation des blocages mécaniques WTAB.'BLOCAGES_MECANIQUES' = (WTAB. 'BLOCAGES_MECANIQUES' ) ET CL3; * actualisation du chargement mécanique WTAB.'CHARGEMENT' = (WTAB.'CHARGEMENT' ) ET CHA2; * Tab1.ok = vrai; sinon mess 'OK'; finsi; i = i + 1; j = j + delta; fin BO3; finproc ; PASAPAS TAB1; * Post traitement JC **pour avoir une idée des déplacements en E1 et E** B = poin V1 proc (0.25 0. (la / 2.)); nbi=dime TAB1.TEMPS; time=prog; depe=prog; repe bouci nbi; indic = &bouci - 1 ; timi=TAB1.TEMPS.indic; time=time inse &bouci timi; toto=TAB1.DEPLACEMENTS.indic; depee=extr toto UY B; depe=depe inse &bouci depee; fin bouci; evdepe=evol vert manu 'temps' time 'depE' depe; TAB2 = table ; TAB2 . 1 = 'MOT' 'TIRR MARQ M PLEIN CARR'; TAB2 . 'TITRE' = table ; TAB2 . 'TITRE' . 1 = 'MOT' 'depl au pt B'; dess (evdepe) 'LEGE' TAB2 'TITX' 'temps (s)' 'TITY' 'deplacement (m)' 'DATE' 'GRIL'; toto0=defo (tab1.deplacements.20) 1. v1 blanc; toto=defo (tab1.deplacements.90) 1. v1 vert; toto1=defo (tab1.deplacements.120) 1. V1 roug; toto2=defo (tab1.deplacements.180) 1. v1 jaun; trac (toto0 et toto et toto1 et toto2); * déplacemet point B NN = DIME (TAB1.TEMPS) ; TEM = PROG ; TT = PROG ;I = 0 ; REPE BO1 (NN) ; T1 = TAB1.TEMPS.(&BO1 - 1) ; TEM = INSE TEM (&BO1) T1; T2 = TAB1.DEPLACEMENTS.(&BO1 - 1) ; T2 = EXTR T2 UY B; TT = INSE TT (&BO1) T2; FIN BO1 ; EV1 = EVOL MANU 'abs' TEM 'cor' TT ; @EXCEL1 EV1 'depp_uy_B.csv' ;