******************************************************* * * * TEST ELAS11 * * * * SOLIDE (DALLE) SOUMIS A UN GRADIENT * * DE TEMPERATURE (Calcul axisymetrique 2D). * * * * Dalle soumise a un gradient de temperature * * s'exercant sur sa hauteur. * * Le champs de temperature s'exprime comme suit : * * * * T(x,y,z) = z * * * * L'ensemble est maille avec des elements massifs * * CUB8, d'une épaisseur de 5mm. * * * * On se propose de calculer la rotation de la section * * de dalle hors de son plan, afin d'établir * * un gradient thermique linéaire équivalent. * * * ******************************************************* * OPTI DIME 3 ELEM CU20 ; * *=====================================================* * geometrie : maillage * *=====================================================* * nel1 = 18 ; larg1 = 0.18 ; larg2 = 1.0 ; nlong = 10 ; l = 5 ; * *r1 = larg1/larg2 ; r1 = 2 ; n1 = nel1 ; n2 = nel1/r1 ; * Points de la section * po = 0. 0. 0. ; po1 = 0. 0. ((-1)*larg1/2) ; po2 = (larg2/2) 0. ((-1)*larg1/2) ; po3 = (larg2/2) 0. (larg1/2) ; po4 = 0. 0. (larg1/2) ; po5 = po3 SYME DROI po1 po4 ; po6 = po2 SYME DROI po1 po4 ; papp1 = 0. 0. ((-1)*larg1/2) ; papp2 = 0. l ((-1)*larg1/2) ; dtest = 'DROI' (n2/2) papp1 papp2 COUL CYAN ; *Droites du contour de la section * d6a2 = 'DROI' (n2/2) po6 po2 ; d2a3 = 'DROI' nel1 po2 po3 ; d3a5 = 'DROI' (n2/2) po3 po5 ; d5a6 = 'DROI' nel1 po5 po6 ; * Contour de la section de dalle * contd = (d6a2 ET d2a3 ET d3a5 ET d5a6 ) ; * TRAC contd ; cext1 = ('DALL' d6a2 d2a3 d3a5 d5a6 ) ; 'ELIM' 1e-6 cext1 ; tout4 = cext1 ; *Surface de la section - Dallage en Mailles * SURF1 = cext1 ; *SURF1 = 'CHAN' tout4 QUADRATIQUE ; 'TRAC' QUAL SURF1 ; *Sections de milieu et de bout à partir de la section *de la dalle * SURFM = SURF1 PLUS (0. (l/2) 0.) COUL JAUNE ; SURF2 = SURF1 PLUS (0. l 0.) COUL VERT ; *Sections d'appui de la dalle - conditions limites * Dapp1 = d6a2 PLUS (0. 0. 0.) COUL CYAN ; Sapp1 = Dapp1 TRAN 2 (0. 0.1 0.) ; Dapp2 = d6a2 PLUS (0. l 0.) COUL CYAN ; Sapp2 = Dapp2 TRAN 2 (0. -0.1 0.) ; TRAC (Sapp1 'ET' Sapp2) ; *Volume de la dalle * VOLD = SURF1 nlong VOLU TRAN (0. l 0.) ; *Modele entier de la dalle * DALLE = (VOLD 'ET' SURFM 'ET' SURF2 'ET' Sapp1 'ET' Sapp2 'ET' Dapp1 'ET' Dapp2 ) ; elim DALLE 0.001 ; TRAC QUAL DALLE ; *REPERE Rx1 = VECT (FORCE (1. 0. 0.) po) 0.2 FX FY FZ 'VERT'; Ry1 = VECT (FORCE (0. 1. 0.) po) 0.2 FX FY FZ 'BLEU'; Rz1 = VECT (FORCE (0. 0. 1.) po) 0.2 FX FY FZ 'ROUG'; Repxyz1 = (Rx1 ET Ry1 ET Rz1) ; TRAC DALLE Repxyz1 ; *=====================================================* * Conditions aux limites * *=====================================================* * * appui simple sur les extremitées de la dalle * CL1 = BLOQ UZ Sapp1 ; CL2 = BLOQ UZ Sapp2 ; * CL3 = BLOQ UX Dapp1 ; CL4 = BLOQ UX Dapp2 ; * CLI1 = CL1 'ET' CL2 'ET' CL3 'ET' CL4 ; *CLI1 = CL1 'ET' CL2 'ET' CL3 ; * *=====================================================* * modele - materiau * * caracteristique - rigidite * *=====================================================* * MODL1 = MODL VOLD MECANIQUE ELASTIQUE ISOTROPE ; MODL2 = MODL d5a6 MECANIQUE ELASTIQUE ISOTROPE ; MATR1 = MATR MODL1 YOUN 31E6 NU 0.2 ALPH 1.E-5; RI1 = RIGI MODL1 MATR1; * RIG1 = RI1 ET CLI1; * *=====================================================* * Chargement thermique * *=====================================================* * * Le solide est soumis a un gradient de temperature * lineaire, s'appliquant sur z . * La temperature est la suivante : * * T(x,y,z) = z ; * * Creation du champs de temperature * CHPX = COOR 1 VOLD ; CHPY = COOR 2 VOLD ; chpz = COOR 3 VOLD ; * CHP4 = CHPZ ; * * Transformation du sous-type "scal" en sous-type "t" * pour chp4. Utilisation de l'operateur "vari". * *LIS10 = PROG 0.18 0.17 0.16 0.15 0.14 0.13 0.12 0.11 0.10 0.09 0.08 0.07 0.06 0.05 0.04 0.03 0.02 0.01 0.0 ; *LIS20 = PROG 34.6 33.7 32.9 32.1 31.4 30.8 30.2 29.7 29.2 28.8 28.4 28.1 28. 27.8 27.6 27.4 27.2 27.1 27. ; * LIS10 = PROG 0.0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 0.11 0.12 0.13 0.14 0.15 0.16 0.17 0.18 ; LIS20 = PROG 27. 27.1 27.2 27.4 27.6 27.8 28. 28.1 28.4 28.8 29.2 29.7 30.2 30.8 31.4 32.1 32.9 33.7 34.6 ; EVOL1 = EVOL MANU T LIS20 y LIS10; DESS EVOL1 ; * CHP5 = VARI CHP4 EVOL1 T; * * Creation du champ de contraintes equivalentes * CHAM10 = THET MODL1 MATR1 CHP5; * * Creation du champ de forces equivalentes * CHAR1 = BSIG MODL1 CHAM10; * *=====================================================* * Resolution : champs de deplacements * *=====================================================* * DEP1 = RESO RIG1 CHAR1; * *=====================================================* * tenseur de déformations * *=====================================================* * CHDEF = EPSI DEP1 MODL1 MATR1 ; * ************************************************************** *opti donn 5 ; ************************************************************** * * Calcul du champs de contraintes totales. * CHAM1 = SIGM MODL1 MATR1 DEP1; * * On retranche les contraintes d'origine thermique * CHAM1 = CHAM1 - CHAM10; * * Interpolation aux noeuds du maillage. * CHAM2 = CHAN NOEUD MODL1 CHAM1; CHPO10 = CHAN CHPO MODL1 CHAM2; * **************************************************************************************************************************** * * POST TRAITEMENT * * Trace facultatif de la deformee et des reactions * AMPVEC = 5.E-2; REA1 = REAC DEP1 RIG1; VECT1 = VECT REA1 AMPVEC FX FY FZ BLEU; TRAC VOLD VECT1 ; * DEF0 = DEFO 0. DEP1 DALLE ; DEF1 = DEFO 5000. DEP1 DALLE ROSE ; TITR 'ELAS11 ; DEFORMEE, REACTIONS'; TRAC (VOLD 'ET' Dapp1 'ET' Dapp2) (DEF0 ET DEF1); TRAC VOLD DEF1 ; TRAC SURF1 (DEF0 ET DEF1); TRAC SURFM (DEF0 ET DEF1); TRAC SURF2 (DEF0 ET DEF1); * sur une partie du maillage * UUA = REDU DEP1 SURF1 ; DEFO1 = DEFO 0. UUA DALLE ; DEFO2 = DEFO -1000. UUA DALLE TURQ ; TRAC SURF1 (DEFO1 'ET' DEFO2); * UUDA = REDU DEP1 d5a6 ; DEFO3 = DEFO 0. UUDA d5a6 ; DEFO4 = DEFO 1000. UUDA d5a6 ROSE ; TRAC d5a6 (DEFO4 'ET' DEFO3) ; * déplacements sous forme de fleches * *OVEC = VECT DEP1 1.e-3; *TRAC OVEC SECT1 ; * iso-composante en déplacement sur section * OEIL1 = 0.0 1.0 0.0 ; * UUX = EXCO DEP1 UX ; TRAC OEIL1 UUX SURF1 ; * UUY = EXCO DEP1 UY ; TRAC OEIL1 UUY SURF1 ; TRAC OEIL1 UUY d2a3 ; * *TRAC OEIL1 UUY SURFM ; *TRAC OEIL1 UUY SURF2 ; * *=====================================================* * extraction d'une composante * * de déplacement et sauvegarde vers des tableaux * *=====================================================* * 'OPTI' 'SORT' 'UXUYUZ.dat' ; 'SORT' 'AVS' d5a6 UUDA ; 'OPTI' 'SORT' 'zz1' ; lelem1 = 'NBEL' d5a6 ('MOTS' SEG2 TRI3 QUA4) ; 'OPTI' 'ACQU' 'UXUYUZ.dat' ; * *************************************************************** *=====================================================* * extraction d'une composante * * de déformation et sauvegarde vers des tableaux * *=====================================================* * EPSx = EXCO CHDEF EPXX ; EPSxm = MAXI EPSx ; LIST EPSxm ; * EPSy = EXCO CHDEF EPYY ; EPSym = MAXI EPSY ; LIST EPSym ; 'OPTI' 'SORT' 'epsi.dat' ; 'SORT' 'AVS' VOLD CHDEF ; 'OPTI' 'SORT' 'zz1' ; lelem1 = 'NBEL' vold ('MOTS' SEG2 TRI3 QUA4) ; 'OPTI' 'ACQU' 'epsilon.dat' ; * opti donn 5 ; * *=====================================================* * extraction d'une contrainte * * et comparaison avec solution analytique * *=====================================================* * * On cherche la contrainte axiale sigzz au noeud A * * * Extraction de la contrainte recherchee. * SMZZA = EXTR CHPO10 SMxy po3; SMZZA = SMZZA / 1.E6; * SMZZA1 = -105; * * Calcul de l'ecart avec la solution de reference. * ERGZZA = 100 * (ABS ((SMZZA - SMZZA1) / SMZZA1)); * opti donn 5 ;