* fichier : ther62.dgibi ************************************************************************ ************************************************************************ SAUT PAGE ; *-------------------------------------------------------------------* * TEST THER62 * * * TEST DES CL DE TEMPERATURE POUR coq4 ET coq3 * * ------------------------------------------------------------ * * * * TEMPERATURE IMPOSEE + CONVECTION * * + FLUX * * + SOURCE * * * * Ce test permet de v{rifier le bon fonctionnement des divers * * op{rateurs thermiques de CASTEM2000 dans le cas des {l{ments * * coq4 et coq3 . * * * * UNE PLAQUE RECTANGULAIRE CONSTITUEE D'{L{ments COQ3 et coq4 est * * soumise @ une temp{rature impos{e @ une de ses extr{mit{s et suc- * * essivement @ une condition de convection forc{e, @ une condition * * de flux impos{ et @ une condition de source volumique impos{e. * * * * Les r{sultats sont compar{s a ceux obtenus * * par un calcul {quivalent utilisant les {l{ments massifs. Le mail- * * -lage n'utilise que des cub8 et est oriente differement pour sim- * * -plifier le calcul. * * On compare les temperatures sur la face superieure; le long de la * * ligne AS2 DS2 * * Test elabore par Mr Dowlatyari 03/91. * * Test modifie le 14/05/91 * * * * * * * *-------------------------------------------------------------------* ******************************************************************** **************** 1/ calcul en massif ******************************* ******************************************************************** *qua3 ; graph = 'N'; TEMPS ; * *------------- CREATION DE LA GEOMETRIE: POINTS SUPPORTS DES ELEMENTS -- * As1= 0.25 0. 0.; AS2 =0.25 1.5 0. ; AS3 = 0.25 3. 0.; BS1= 0.25 0. 2.; BS2 =0.25 1.5 2.; BS3 = 0.25 3. 2.; CS1= 0.25 0. 4.; CS2 =0.25 1.5 4; CS3 = 0.25 3. 4.; DS1= 0.25 0. 6.; DS2 =0.25 1.5 6.; DS3 = 0.25 3. 6.; *------------- CREATION DES SURFACES ET VOLUMES ------------------ - * DAD = AS3 d 3 DS3; DDD = DS3 D 2 DS1; DDA = DS1 D 3 AS1; DAA = AS1 D 2 AS3; V = -0.5 0. 0.; OEIL = 5. 1.5 3.; si ( NEG GRAPH 'N' ); trac oeil qtot; trac oeil scon; trac oeil sinf; finsi; *--------------- DONNEES DU PROBLEME DE THERMIQUE ----------------------- * *-------------- MODELISATION ------------------------------------------- * * *--------------- DONNEES DES CARACTERISTIQUES DU MATERIAU --------------- * * *--------------- CREATION DES MATRICES DE CONDUCTIVITE ------------------ * CND1 = CONDUCTIVITE MOD1 MAT1 ; CND2 = CONDUCTIVITE MOD2 MAT2 ; * *--------------- TEMPERATURES IMPOSEES: BLOQUE + DEPI ------------------- * BB1 = BLOQUE SINF 'TINF' ; BB2 = BLOQUE SINF 'TSUP' ; BB3 = BLOQUE SINF 'T ' ; BB = BB1 ET BB2 ET BB3 ; EE = EE1 ET EE2 ET EE3 ; * *--------------- FLUX EQUIVALENTS A LA CONDITION DE CONVECTION ---------- * CNV1 = CONVECTION MOD2 MAT2 'T' 200. ; * *--------------- FLUX IMPOSE -------------------------------------------- * * *--------------- SOURCE DE CHALEUR -------------------------------------- * S1 = SOURCE MOD1 MAT1 1000. SSOU ; * *---------------- ASSEMBLAGE DES TERMES DE CONDUCTIVITE ----------------- * CCC1 = CND1 ET CND2 ET BB ;CCC2 = CND1 ET BB ; * *---------------- ASSEMBLAGE DES TERMES DE FLUX EQUIVALENTS ------------- * FFF1 = EE ET CNV1 ;FFF2 = EE ET FLU1 ;FFF3 = EE ET S1 ; * *---------------- 1ERE RESOLUTION : CONVECTION -------------------------- * CHTER1 = RESOUDRE CCC1 FFF1 ; * *---------------- 2EME RESOLUTION : FLUX -------------------------------- * CHTER2 = RESOUDRE CCC2 FFF2 ; * *---------------- 3EME RESOLUTION: SOURCE ------------------------------- * CHTER3 = RESOUDRE CCC2 FFF3 ; * *---------------- POST-TRAITEMENT: - *---------------- EXTRACTION DES VALEURS CARACTERISTIQUES ----------- * *LISTE CHTER1 ; *LISTE CHTER2 ; *LISTE CHTER3 ; * *LIST LCON; LIST LFLU; LIST LSOU; ******************************************************************* ********************* 2/ CALCUL EN COQ3 ET COQ4 ******************* ******************************************************************* TITRE 'C.L. DE TEMPERATURE - ELEMENTS COQ3 COQ4' ; *TRI3 ; TEMPS ; * *------------- CREATION DE LA GEOMETRIE: POINTS SUPPORTS DES ELEMENTS -- * A1 = 0. 0. 0. ; A2 = 0. 1.5 0. ; B1 = 1.4142 0. 1.4142 ; B2 = 1.4142 1.5 1.4141 ; C1 = 2.8284 0. 2.8284 ; C2 = 2.8284 1.5 2.8284 ; D1 = 4.2426 0. 4.2426 ; D2 = 4.2426 1.5 4.2426 ; A3 = 0. 3. 0. ; B3 = 1.4142 3. 1.4142 ; C3 = 2.8284 3. 2.8284 ; D3 = 4.2426 3. 4.2426 ; *------------- CREATION DES QUADRANGLES ET TRIANGLES ------------------- * * QTOT = QQ01 ET QQ02 ET QQ03 ET QQ04 ET QQ05 ET QQ06 ET QQ07 ET QQ08 ; SI ( NEG GRAPH 'N'); OEIL = -1000. 0. 1000. ; TRACE OEIL QTOT ; FINSI; QSOUR = QQ04 ET QQ08 ; * DD3 = D 1 A3 A2 ; DD4 = D 1 A2 A1 ; L_INF = DD3 ET DD4 ; * S_COT = QQ04 ET QQ08 ; * * ATTENTION : il faut r{orienter les {l{ments (OPERATEUR ORIENT) * si on n'est pas sur que tous les elements sont orientes d'une * facon coh{rante * *--------------- DONNEES DU PROBLEME DE THERMIQUE ----------------------- * *-------------- MODELISATION ------------------------------------------- * * *--------------- DONNEES DES CARACTERISTIQUES DU MATERIAU --------------- * * *--------------- CREATION DES MATRICES DE CONDUCTIVITE ------------------ * CND1 = CONDUCTIVITE MOD1 MAT1 ; CND2 = CONDUCTIVITE MOD2 MAT2 ; * *--------------- TEMPERATURES IMPOSEES: BLOQUE + DEPI ------------------- * BB1 = BLOQUE L_INF 'TINF' ; BB2 = BLOQUE L_INF 'TSUP' ; BB3 = BLOQUE L_INF 'T ' ; BB = BB1 ET BB2 ET BB3 ; EE = EE1 ET EE2 ET EE3 ; * *--------------- FLUX EQUIVALENTS A LA CONDITION DE CONVECTION ---------- * CNV1 = CONVECTION MOD2 MAT2 'TSUP' 200.; * *--------------- FLUX IMPOSE -------------------------------------------- * * *--------------- SOURCE DE CHALEUR -------------------------------------- * S1 = SOURCE MOD1 MAT1 1000. QSOUR ; * *---------------- ASSEMBLAGE DES TERMES DE CONDUCTIVITE ----------------- * CCC1 = CND1 ET CND2 ET BB ;CCC2 = CND1 ET BB ; * *---------------- ASSEMBLAGE DES TERMES DE FLUX EQUIVALENTS ------------- * FFF1 = EE ET CNV1 ;FFF2 = EE ET FLU1 ;FFF3 = EE ET S1 ; * *---------------- 1ERE RESOLUTION : CONVECTION -------------------------- * CHTER1 = RESOUDRE CCC1 FFF1 ; * *---------------- 2EME RESOLUTION : FLUX -------------------------------- * CHTER2 = RESOUDRE CCC2 FFF2 ; * *---------------- 3EME RESOLUTION: SOURCE ------------------------------- * CHTER3 = RESOUDRE CCC2 FFF3 ; * *---------------- POST-TRAITEMENT: TRACE DES CHAMPS RESULTATS ----------- *---------------- ET EXTRACTION DES VALEURS CARACTERISTIQUES ----------- * *LISTE CHTER1 ; *LISTE CHTER2 ; *LISTE CHTER3 ; * OEIL= -1000. 0. 1000. ; * SI ( NEG GRAPH 'N'); TITRE 'ELEMENTS COQ3,COQ4 : CONVECTION' ; TRACER OEIL QTOT CHTR1I ; TITRE 'ELEMENTS COQ3,COQ4 : FLUX IMPOSE' ; TRACER OEIL QTOT CHTR2I ; TITRE 'ELEMENTS COQ3,COQ4 : SOURCE' ; TRACER OEIL QTOT CHTR3I ; FINSI; * * ************************************************************* ********************* 3/ comparaison entre les 2 calculs **** ************************************************************* si ( NEG graph 'N'); TAB = 'TABLE' ; TAB.2 = 'MARQ ETOI' ; TITRE 'COQ3 ET COQ4 ;TEMPERATURE + CONVECTION ;+ : MASSIF,* : COQUE' ; DESSIN TAB (EV1R ET EV1C ) ; DESSIN TAB (EV2R ET EV2C ) ; TITRE 'COQ3 ET COQ4 ;TEMPERATURE + SOURCE ;+ MASSIF,* : COQUE '; DESSIN TAB (EV3R ET EV3C ) ; TEMPS ; finsi; ************************************************************* *************** 4/ TEST DE FONCTIONNEMENT ******************* ************************************************************* ERCA= ABS (( T1A-T1Aq)/T1A); ERCB= ABS (( T1B-T1Bq)/T1B); ERCC= ABS (( T1C-T1Cq)/T1C); ERCD= ABS (( T1D-T1Dq)/T1D); ERFA= ABS (( T2A-T2Aq)/T2A); ERFB= ABS (( T2B-T2Bq)/T2B); ERFD= ABS (( T2D-T2Dq)/T2D); ERSA= ABS (( T3A-T3Aq)/T3A); ERSB= ABS (( T3B-T3Bq)/T3B); ERSC= ABS (( T3C-T3Cq)/T3C); ERSD= ABS (( T3D-T3Dq)/T3D); ERSC ERSD; temps; SI ( ERMAX <EG 0.05 ); SINON; FINSI; fin;
© Cast3M 2003 - Tous droits réservés.
Mentions légales