********************************************************************************************************************** * * TITR'Test_ajout_élément'; OPTI ECHO 0 DIME 3 ELEM CUB8; nb_mail = 1; * COORDONNNEES P1 = 0. 0. 0.; P2 = 1. 0. 0.; P3 = 1. 0. 3.; P4 = 0. 0. 3;; P5 = 0. 2. 0.; P6 = 1. 2. 0.; P7 = 1. 2. 3.; P8 = 0. 2. 3; * FACE AVANT L1L2 = P1 DROI nb_mail P2 ; L2L3 = P2 DROI nb_mail P3 ; L3L4 = P3 DROI nb_mail P4; L4L1 = P4 DROI nb_mail P1; Contavan = L1L2 et L2L3 et L3L4 et L4L1; Surfavan = SURF PLAN Contavan COUL ROUGE ; * FACE ARRIERE L6L5 = P6 DROI nb_mail P5 ; L5L8 = P5 DROI nb_mail P8 ; L8L7 = P8 DROI nb_mail P7; L7L6 = P7 DROI nb_mail P6; Contarri = L6L5 et L5L8 et L8L7 et L7L6; Surfarri = SURF PLAN Contarri COUL ROUGE ; *TRAC (Surfavan et Surfarri); * FACE GAUCHE L4L8 = P4 DROI nb_mail P8 ; L8L5 = P8 DROI nb_mail P5 ; L5L1 = P5 DROI nb_mail P1; L1L4 = P1 DROI nb_mail P4; Contgauc = L4L8 et L8L5 et L5L1 et L1L4; Surfgauc = SURF PLAN Contgauc COUL BLEU ; *TRAC Surfgauc; * FACE DROITE L7L3 = P7 DROI nb_mail P3 ; L3L2 = P3 DROI nb_mail P2 ; L2L6 = P2 DROI nb_mail P6; L6L7 = P6 DROI nb_mail P7; Contdroi = L7L3 et L3L2 et L2L6 et L6L7; Surfdroi = SURF PLAN Contdroi COUL BLEU ; * FACE INFERIEURE L2L1 = P2 DROI nb_mail P1 ; L1L5 = P1 DROI nb_mail P5 ; L5L6 = P5 DROI nb_mail P6; L6L2 = P6 DROI nb_mail P2; Continfe = L2L1 et L1L5 et L5L6 et L6L2; Surfinfe = SURF PLAN Continfe COUL VERT ; * FACE SUPERIEURE L4L3 = P4 DROI nb_mail P3 ; L3L7 = P3 DROI nb_mail P7 ; L7L8 = P7 DROI nb_mail P8; L8L4 = P8 DROI nb_mail P4; Contsupe = L4L3 et L3L7 et L7L8 et L8L4; Surfsupe = SURF PLAN Contsupe COUL VERT ; Surftot = (Surfavan et Surfarri et Surfgauc et Surfdroi et Surfinfe et Surfsupe); *TRAC Surftot; Surftot2 = ELIM Surftot 0.001; Votot = Surftot2 VOLU; Trac Votot; ******************************************************************* * DEFINITION DES MODELES ******************************************************************* **************************************************** * A ) Beton **************************************************** Motot= MODE Votot MECANIQUE ELASTIQUE ISOTROPE; ******************************************************************* ******************************************************************* ** DEFINITION DES MATERIAUX **** ******************************************************************* ******************************************************************* ********************************************* *A ) BETON ********************************************* Matot = MATE Motot 'YOUN' 31.E9 'NU' 0.2 ; ******************************************************************* * DEFINITION DES CONDITIONS AUX LIMITES ******************************************************************* * surfaces latérales : pression de confinement CLgauc = BLOQ UX Surfgauc ; CLdroi = BLOQ UX Surfdroi ; CLavan = BLOQ UY Surfavan ; CLarri = BLOQ UY Surfarri ; * surfaces supérieurs et inférieures : déplacement imposé CLinfe = BLOQ UZ Surfinfe ; CLsupe = BLOQ UZ Surfsupe ; CLtot = CLinfe et CLsupe; ******************************************************************* * DEFINITION DU CHARGEMENT (pression augmentée jusqu'à X kPa puis déplacement vertical imposé * Remarque : cette technique permet de déduire la contrainte max mais pas la déformation max ******************************************************************* LTPSCHA = PROG 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ; FORCCHA = PROG 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20; DEPLCHA = PROG 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39; *pilotage déplacement IncrDepl = -1.E-4; F1 = DEPI CLSupe IncrDepl; *evolution pour le chargement par déplacement imposé EVCHA = EVOL MANU 'TEMPS' LTPSCHA 'DEPLACEMENT' DEPLCHA ; CHSTAI1 = CHAR F1 EVCHA 'MECA' ; ********************************************************************* *pilotage effort IncrForc = 5.0E-3; F2 = FSUR 'MASS' Motot Surfgauc (0. 0. IncrForc); F4 = FSUR 'MASS' Motot Surfavan (0. 0. IncrForc); F6 = F2 et F4; *evolution pour le chargement par effort imposé EVCHA = EVOL MANU 'TEMPS' LTPSCHA 'FORCE' FORCCHA ; CHSTAI2 = CHAR F6 EVCHA 'MECA' ; *DESS EVCHA ; FORCCHA2 = FORCCHA * IncrForc; EVCHA2 = EVOL MANU 'TEMPS' LTPSCHA 'DEPLACEMENT' FORCCHA2 ; DESS EVCHA2 ; ************************************************************************* CHSTAI3 = CHSTAI1; LTPSCHA2 = LTPSCHA; compteur = 1; VAL_TOL = 1.E-4; NVERIF = 0; PTAB3 = TABLE ; PTAB3.'TEMPS_CALCULES' = LTPSCHA ; PTAB3.'TEMPS_SAUVES' = LTPSCHA2 ; PTAB3.'BLOCAGES_MECANIQUES' = CLtot ; PTAB3.'CARACTERISTIQUES' = MaTot ; PTAB3.'CHARGEMENT' = CHSTAI3; PTAB3.'MODELE' = MoTot ; PTAB3.'PRECISION' = VAL_TOL ; PTAB3.GRAND_DEPLACEMENTS=VRAI; PTAB3.FTOL = VAL_TOL; 5E-3; PTAB3.ECONOMIQUE = VRAI ; PTAB3.PROCESSEURS = 'AUTOMATIQUE'; PTAB2 = TABLE ; PTAB3.'MES_SAUVEGARDES'= PTAB2; PTAB2.DEFTO = VRAI ; PTAB2.DEFIN = VRAI; PASAPAS PTAB3 ; *********************************************** * Deformées *********************************************** NDIM = DIME PTAB3.DEPLACEMENTS ; SNAP_DEF = 60; GEOTOT = EXTR PTAB3.'MODELE' 'MAIL'; DEP0 = PTAB3.DEPLACEMENTS.0; DEP1 = PTAB3.DEPLACEMENTS.60; DEF0 = DEFO GEOTOT DEP0 0. BLAN; DEF1 = DEFO GEOTOT DEP1 30. ROUG; TRAC (DEF0 ET DEF1);