Numérotation des lignes :

* fichier :  prot.dgibi*********************************************************************************************************************************************************************************************************************                                                                  **    Test de la procedure PROT : Projection de Temperature d'un    **    massif sur une coque                                          **                                                                  **********************************************************************OPTI ECHO 0;GRAPH ='N';*OPTI DIME 3 ;OPTI ELEM SEG2;**==================================================================**             geometrie : maillage                                 **==================================================================*** Dimension en metres** Points*P1 = -3. -1. 5.; P1B = -3. -1. 5.5; P1C = -3. -1. 4.5;P2 =  0. -1. 5.; P2B =  0. -1. 5.5; P2C =  0. -1. 4.5;P3 =  3. -1. 5.; P3B =  3. -1. 5.5; P3C =  3. -1. 4.5; P4 = -3.  1. 5.; P4B = -3.  1. 5.5; P4C = -3.  1. 4.5; P5 =  0.  1. 5.; P5B =  0.  1. 5.5; P5C =  0.  1. 4.5; P6 =  3.  1. 5.; P6B =  3.  1. 5.5; P6C =  3.  1. 4.5;  L1 = P1 D 5 P2 D 3 P5 D 5 P4 D 3  ;L2 = P6 D 5 P5 D 3 P2 D 5 P3 ;L3 = P3 D 3 P6  ;ELIM 1.E-3 (L1 ET L2);L1B = P1B D 5 P3B D 3 P6B D 5 P4B D 3  ;L2B = P3B D 3 P6B  ;ELIM 1.E-3 (L1B ET L2B);*L1C = P1C D 5 P3C D 3 P6C D 5 P4C D 3 ;L2C = P1C D 3 P4C  ;*ELIM 1.E-3 (L1C ET L2C); *  Creation des objets OPTI ELEM TRI6 ;PLAN1 = SURF  L1 PLANE ;OPTI ELEM QUA8 ;PLAN2 = SURF (L2 ET L3) PLANE ; OPTI ELEM TRI3 ; PLAN1B = SURF L1B PLANE ;PLAN1C = PLAN1B 'PLUS' (0. 0. -1.) ;ELIM PLAN1C L2C 0.00001;*PLAN1C = SURF L1C PLANE ; GEO  = (PLAN1 ET PLAN2) COUL ROUG ;ELIM GEO 1.E-7; OPTI ELEM PRI6 ; GEOV = PLAN1B VOLU PLAN1C ; * Definition du modele coque*---------------------------- MO1 = MODELI PLAN1 THERMIQUE COQ6 ;MO2 = MODELI PLAN1 THERMIQUE COQ6 ;MO3 = MODELI PLAN2 THERMIQUE COQ8 ;MO4 = MODELI PLAN2 THERMIQUE COQ8 ; MO13 =  MO1 ET MO3 ;MO24 =  MO2 ET MO4 ;MO = MO13 ET MO24 ; * Definition du materiau et de l'epaisseur MA1 = MATE MO13 EPAI 0.5 EXCE 2.5E-1 ;MA2 = MATE MO24 EPAI 0.5 EXCE -2.5E-1 ;MA  = MA1 ET MA2 ; * Definition du modele massif*----------------------------- MOV = MODELI GEOV THERMIQUE  ; * Definition du materiau et de l'epaisseur MAV = MATE MOV K 15. 'C' 460. RHO 8000. ; * Definition des conditions aux limites*--------------------------------------- CL1 = BLOQ T L2B ;CL2 = BLOQ T L2C ; RIT =  CL1 ET CL2; * Temperature imposee*--------------------- TEM1 = DEPI CL1 500.;TEM2 = DEPI CL2 0.; * Matrice conductivite*---------------------- CON1 = CONDUCTIVITE MOV MAV ;CONT = CON1 ET RIT ; * Second membre*--------------- TEMT = TEM1 ET TEM2 ; * Resolution du systeme*----------------------- CHT1 = RESO CONT TEMT ; CHT2 = CHAN CHAM CHT1 GEOV ; *trace CHT1 (PLAN1B ET PLAN1C) ; * Projection sur la coque*------------------------- *CHT3 = PROI MO MA CHT2 ; CHT31 = CALP  MO13 MA1 CHT2 ;CHT32 = CALP  MO24 MA2 CHT2 ;CHP31 = CHANGE 'CHPO'  MO13 CHT31;                        CHP32 = CHANGE 'CHPO'  MO24 CHT32;                         CHT31T = EXCO CHT31 T ; CHT31TI = EXCO CHT31 TINF ; CHT31TS = EXCO CHT31 TSUP ; CHT32T = EXCO CHT32 T ; CHT32TI = EXCO CHT32 TINF ; CHT32TS = EXCO CHT32 TSUP ; *SI (NEG GRAPH 'N') ;  trace CHT32TI MO13 ;  trace CHT32T MO13 ;  trace CHT32TS MO13 ;  trace CHT31TI MO13 ;  trace CHT31T MO13 ;  trace CHT31TS MO13 ;FINSI ;*T31S3 = EXTR CHP31 TSUP P3;T32I1 = EXTR CHP32 TINF P1;** Les temperatures TINF sur la coque associee au modele MO13 doit* etre nulle au point P3 ( Condition  au limite), et la temperature* TSUP sur la coque associee au modele MO24 doit etre egale a 500 en P1*DIF = (ABS(T32I1)) + (ABS(T31S3 - 500.)) ;SI (DIF &lt;EG 1.E-2) ;  ERRE 0 ;SINO ;  ERRE 5 ;FINSI ; FIN ;

© Cast3M 2003 - Tous droits réservés.
Mentions légales