Télécharger la plaquette Cast3M


fr en

Télécharger nsclim.procedur

Retour à la liste

Numérotation des lignes : 

  1. * NSCLIM    PROCEDUR  MAGN      05/07/25    21:15:01     5133           
  2. 'DEBPROC' NSCLIM                                                       ;
  3.  
  4.  ARGU tit*'MOT' KIMPR/'LOGIQUE' NOMQ*'MOT'                             ;
  5.  
  6.  Si (Exist KIMPR)                                                      ;
  7.  TIMPR = KIMPR                                                         ;
  8.  Sinon                                                                 ;
  9.  TIMPR = VRAI                                                          ;
  10.  FinSi                                                                 ;
  11.  
  12.  Si TIMPR                                                              ;
  13.  mess '#### NSCLIM #######################################'
  14.       '#########################'                                      ;
  15.  FinSi                                                                 ;
  16.  
  17.  lnomq = mots 'DEBIT' 'VITESSE' 'TEMPERATURE' 'ADHERENCE'              ;
  18.  Si (NON (EXIST lnomq nomq))                                           ;
  19.  Mess '*********************************************************'      ;
  20.  Mess ' ERREUR   ERREUR   ERREUR   ERREUR   ERREUR   ERREUR     '      ;
  21.  Mess '                                                         '      ;
  22.  Mess 'L entree ' NOMQ   ' n existe pas '                              ;
  23.  Mess 'Choisir : VITESSE ou DEBIT ou TEMPERATURE ou ADHERENCE   '      ;
  24.  Mess '*********************************************************'      ;
  25.  erreur 21                                                             ;
  26.  quitter NSCLIM                                                        ;
  27.  Finsi                                                                 ;
  28.  
  29.  Si((EGA NOMQ 'VITESSE') ou (EGA NOMQ 'DEBIT')
  30.                          ou (EGA NOMQ 'TEMPERATURE') )     ;
  31.   ARGU Udebt/'EVOLUTION'                                   ;
  32.   Si(NON (Exist Udebt))                                    ;
  33.   ARGU Udeb/'FLOTTANT'                                     ;
  34.   FinSi                                                    ;
  35.  
  36.   ARGU    NCO*'MOT'             Tps*'FLOTTANT'             ;
  37.  
  38.   SWIRL=FAUX                                               ;
  39.   Si(EGA NCO 'SWIRL')                                      ;
  40.   SWIRL=VRAI                                               ;
  41.   asw = tps                                                ;
  42.   ARGU    NCO*'MOT'             Tps*'FLOTTANT'             ;
  43.   Finsi                                                    ;
  44.  
  46.   $mf*'MMODEL'   $sentr*'MMODEL'    clim/'CHPOINT'  NOMG*'MOT'
  47.   NOMM/'MOT'  sgm/'FLOTTANT' portee/'FLOTTANT'             ;
  48.  
  49.   Si (EXIST NOMM)                                          ;
  50.   NOMT = NOMM                                              ;
  51.   Sinon                                                    ;
  52.   NOMT='XXXX'                                              ;
  53.   FinSi                                                    ;
  54.  FinSi                                                     ;
  55.  
  56.  Si(EGA NOMQ 'ADHERENCE')                                 ;
  57.   ARGU    NCO*'MOT'
  58.   $mf*'MMODEL'   $sentr*'MMODEL'    clim/'CHPOINT'  NOMG*'MOT' ;
  59.  FinSi                                                     ;
  60.  
  61.  lnomg = mots 'LARGEMENT' 'STRICTEMENT'                     ;
  62.  Si (NON (EXIST lnomg nomg))                                ;
  63.  Mess '*********************************************************' ;
  64.  Mess ' ERREUR   ERREUR   ERREUR   ERREUR   ERREUR   ERREUR     ' ;
  65.  Mess '                                                         ' ;
  66.  Mess 'L entree ' NOMG   ' n existe pas '                         ;
  67.  Mess 'Choisir : LARGEMENT ou STRICTEMENT                       ' ;
  68.  Mess '*********************************************************' ;
  69.  erreur 21 ;
  70.  quitter NSCLIM ;
  71.  Finsi ;
  72.  
  73. *
  74. *  NOMQ = DEBIT        Débit imposé    (Débit entrant >0)
  75. *         VITESSE      Vitesse imposée (Vitesse entrante >0)
  76. *         <SWIRL> alfa  swirl
  77. *         TEMPERATURE  Température imposée (ou scalaire)
  78. *         ADHERENCE    u=v=w=0     imposées
  79. *         FPAROI       u.n=0  et Fparoi
  80. *         SYMETRIE     u.n=0 et rien si vitesse rien si Température
  81. *         SORTIE       -> p=0
  82. *         PRESSION     -> p=p0
  83. *  NCO    nom de la composante
  84. *
  85. *  NOMG = STRICTEMENT
  86. *         LARGEMENT
  87. *
  88. *  NOMT = SGE
  89. *         XXXXXXXXX
  90. *
  91. *
  92.  
  93. *!!!!!!! Petites vérifications !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  94.  
  95. * 1/ sentr doit être inclus dans mf
  96. * 2/ Faire le 2D
  97. * 3/ Traiter les cas des TRI7 -> ELNO QUAD ou MSOMMET
  98. * 5/ Objet Chargement
  99. * 6/ Completer les CL FPAROI SYMETRIE SORTIE PRESSION
  100. * 7/ Dans le cas Strictement le débit est faux
  101.  
  102. *!!!!!!! Petites vérifications !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  103.  
  104. *##############################################################
  105. * Préparation géométrie calcul normales et tangentes
  106.  AT = doma $sentr 'TABLE'                                   ;
  107.  Si (EXIST AT 'nx' )                                        ;
  108.  
  109.  
  110.    nx = AT.'nx'                                             ;
  111.    ny = AT.'ny'                                             ;
  112.  
  113.    tx = AT.'tx'                                             ;
  114.    ty = AT.'ty'                                             ;
  115.  
  116.   Si(EGA (VALE DIME) 3)                                     ;
  117.    nz = AT.'nz'                                             ;
  118.    tz = AT.'tz'                                             ;
  119.  
  120.    ax = AT.'ax'                                             ;
  121.    ay = AT.'ay'                                             ;
  122.    az = AT.'az'                                             ;
  123.   FinSi                                                     ;
  124.  
  125.  Sinon                                                      ;
  126.  
  127.  mess tit ': Préparation géométrie : calcul normales et tangentes';
  128.  
  129.    vol   = doma $mf maillage                                ;
  130.    Sentr = doma $sentr maillage                             ;
  131.    volr = vol elem 'APPUYE' 'LARGEMENT' Sentr               ;
  132.    sent = Sentr                                             ;
  133.    inefmd = (doma $mf 'TABLE').'INEFMD';                    ;
  134.    mess 'inefmd=' inefmd ;
  135.    Si (ega inefmd 3)                                        ;
  136.     $volr = model volr 'NAVIER_STOKES' 'MACRO'              ;
  137.     volr = doma $volr maillage                              ;
  138.     $sentr= model sentr 'NAVIER_STOKES' 'MACRO'             ;
  139.     sent = doma $sentr maillage                             ;
  140.     volr = volr elem 'APPUYE' 'LARGEMENT' Sent              ;
  141.    FinSi                                                    ;
  142.    mo   = mode volr mecanique elastique                     ;
  143.    chp = pres mass mo Sent  1.                              ;
  144.    Si (EGA (vale mode) 'AXIS' )                             ;
  145.    fx = 'FR'                                                ;
  146.    fy = 'FZ'                                                ;
  147.    Sinon                                                    ;
  148.    fx = 'FX'                                                ;
  149.    fy = 'FY'                                                ;
  150.    finSi                                                    ;
  151.    chnm = (((exco chp fx)**2)+((exco chp fy)**2))           ;
  152.   Si(EGA (VALE DIME) 3)                                     ;
  153.    chnm = chnm + ((exco chp fz)**2)                         ;
  154.   FinSi                                                     ;
  155.    chnm = chnm **0.5                                        ;
  156.    chn = chp / chnm                                         ;
  157.  
  158.    nx = (exco fx chn) ;
  159.    ny = (exco fy chn) ;
  160.  
  161.   Si(EGA (VALE DIME) 3)                                     ;
  162.    nz = (exco fz chn) ;
  163. *
  164.    ax = (maxi (coor 1 sent)) - (mini (coor 1 sent));
  165.    ay = (maxi (coor 2 sent)) - (mini (coor 2 sent));
  166.    az = (maxi (coor 3 sent)) - (mini (coor 3 sent));
  167.    av = ax ay az;
  168.    nav = norm av ;
  169.    av = av * (1./nav) ;
  170.    ax = manu chpo sent 'SCAL' (coor 1 av);
  171.    ay = manu chpo sent 'SCAL' (coor 2 av);
  172.    az = manu chpo sent 'SCAL' (coor 3 av);
  173. *
  174.    tx = (ny * az) - (nz * ay) nomc 'UX' ;
  175.    ty = (nz * ax) - (nx * az) nomc 'UY' ;
  176.    tz = (nx * ay) - (ny * ax) nomc 'UZ' ;
  177.    cht = tx + ty + tz ;
  178.    ncht = PSCA cht cht (MOTS 'UX' 'UY' 'UZ') (MOTS 'UX' 'UY' 'UZ');
  179.    nicht = inve ncht ;
  180.    tx = (tx nomc 'SCAL') * nicht ;
  181.    ty = (ty nomc 'SCAL')* nicht ;
  182.    tz = (tz nomc 'SCAL') * nicht ;
  183. *
  184.    ax = (ny * tz) - (nz * ty) ;
  185.    ay = (nz * tx) - (nx * tz) ;
  186.    az = (nx * ty) - (ny * tx) ;
  187. *
  188.    nx = nomc 'UX' nx ;
  189.    ny = nomc 'UY' ny ;
  190.    nz = nomc 'UZ' nz ;
  191. *
  192.    tx = nomc 'UX' tx ;
  193.    ty = nomc 'UY' ty ;
  194.    tz = nomc 'UZ' tz ;
  195. *
  196.    ax = nomc 'UX' ax ;
  197.    ay = nomc 'UY' ay ;
  198.    az = nomc 'UZ' az ;
  199. *
  200.    AT.'nx' = nx                                  ;
  201.    AT.'ny' = ny                                  ;
  202.    AT.'nz' = nz                                  ;
  203.    MNC=1                                         ;
  204.    Si( ((ABS nx MAXI) > (ABS ny MAXI)) et
  205.    ((ABS nx MAXI) > (ABS nz MAXI)) ); MNC=1;FinSi;
  206.    Si( ((ABS ny MAXI) > (ABS nx MAXI)) et
  207.    ((ABS ny MAXI) > (ABS nz MAXI)) ); MNC=2;FinSi;
  208.    Si( ((ABS nz MAXI) > (ABS ny MAXI)) et
  209.    ((ABS nz MAXI) > (ABS nx MAXI)) ); MNC=3;FinSi;
  210.    AT.'MNC' = MNC                                ;
  211.  
  212.    AT.'tx' = tx                                  ;
  213.    AT.'ty' = ty                                  ;
  214.    AT.'tz' = tz                                  ;
  215.  
  216.    AT.'ax' = ax                                  ;
  217.    AT.'ay' = ay                                  ;
  218.    AT.'az' = az                                  ;
  219.  
  220.    AT.'chn'  = nx + ny + nz                      ;
  221.    AT.'cht'  = tx + ty + tz                      ;
  222.    AT.'cha'  = ax + ay + az                      ;
  223.   FinSi                                          ;
  224.  
  225.   Si(EGA (VALE DIME) 2)                          ;
  226.    tx = nomc 'UX' ny                             ;
  227.    ty = nomc 'UY' nx ;                           ;
  228.    nx = nomc 'UX' nx                             ;
  229.    ny = nomc 'UY' ny                             ;
  230.    MNC=1                                         ;
  231.    Si( ((MAXI nx) > (MAXI ny)) ); MNC=2; FinSi   ;
  232.    AT.'MNC' = MNC                                ;
  233.  
  234.    AT.'tx' = tx                                  ;
  235.    AT.'ty' = ty                                  ;
  236.    AT.'nx' = nx                                  ;
  237.    AT.'ny' = ny                                  ;
  238.  
  239.    AT.'chn'  = nx + ny                           ;
  240.    AT.'cht'  = tx + ty                           ;
  241.   FinSi                                          ;
  242.  
  243.    AT.'Sentr'=Sentr                              ;
  244.    AT.'Sentr1'= chan Sentr  'POI1'               ;
  245.    Si (EGA (VALE DIME) 3)                        ;
  246.    CSentr = cont Sentr coul VERTE                           ;
  247.    Sinon                                                    ;
  248.    CSentr = MANU 'SEG2'
  249.               (Sentr POINT 'INITIAL')  (Sentr POINT 'FINAL');
  250.    FinSi                                                    ;
  251.    CSentr =  CSentr COUL ROUGE                              ;
  252.    CSentr1= chan CSentr 'POI1' COUL ROUGE                   ;
  253.    Si('EGA' (VALE MODE ) 'AXIS')                            ;
  254.    xstr1= coor 1 CSentr1                                    ;
  255.    dr =((maxi xstr1) - (mini xstr1))/(nbel CSentr1)/100.    ;
  256.    pax = (CSentr1 'POINT' 'DROIT' (0. 0.) (0. 1.e5) dr)     ;
  257.    Csentr1= diff pax Csentr1                                ;
  258.    FinSi                                                    ;
  259.  
  260.    Sentro = diff CSentr1 AT.'Sentr1' COUL VERT              ;
  261.    AT.'Sentro' = Sentro                                     ;
  262.    AT.'Csentr1'= Csentr1                                    ;
  263.    AT.'Cgrav'  = bary (Sentro et  Csentr1)                  ;
  264.  
  265.    Si TIMPR                                                 ;
  266.    trace  (Sentro et Csentr1) TITR
  267.    (chai tit '  Intérieur et contour')                      ;
  268.    trace  (AT.'Cgrav' et Csentr1) TITR
  269.    (chai tit '  Centre de gravité et contour')              ;
  270.    FinSi                                                    ;
  271.  
  272.  Si(EGA (nbel Sentro) 0 )                                 ;
  273.  Mess '*********************************************************' ;
  274.  Mess ' ERREUR   ERREUR   ERREUR   ERREUR   ERREUR   ERREUR     ' ;
  275.  Mess '                                                         ' ;
  276.  Mess tit '  L intérieur est vide                               ' ;
  277.  Mess '*********************************************************' ;
  278.  erreur 21 ;
  279.  quitter NSCLIM                                             ;
  280.  Finsi ;
  281.  
  282. * On en déduit que connaissant q on a u = q*q1              ;
  283.    Setr   = mesu Sentr                                      ;
  284.    AT.'Setr' = Setr                                         ;
  285. * Calcul du facteur géométrique EF
  286. *  list AT.'chn';
  287.    v1 = kcht $Sentr vect sommet (redu (AT.'chn') Sentro)    ;
  288. *  v1 = kcht $Sentr vect sommet (0. 0. 1.)   ;
  289. *  list v1 ;
  290. *  ung = vect v1 0.1 ux uy uz jaune ;
  291. *  trace ung (doma $Sentr maillage)  ;
  292. *
  293.    q1 = dbit $Sentr v1 'IMPR'                               ;
  294. *  trace (doma $Sentr maillage) TITR 'Sentr';
  295. *  list q1 ;
  296.    AT.'facgeo' = Setr/q1 'ABS'                              ;
  297.    mess tit ' Surface = ' Setr
  298.                      ' Facteur géométrique '   AT.'facgeo'  ;
  299.    AT.'Sentro' = Sentro                                     ;
  300.    AT.'CSentr' = CSentr                                     ;
  301.  
  302.  FinSi                                                      ;
  303. *##############################################################
  304.  
  305.    Sentr   = AT.'Sentr'                                     ;
  306.    Setr    = AT.'Setr'                                      ;
  307.  
  308.  Si(NON(EGA NOMQ 'ADHERENCE'))                              ;
  309.      IK1=0                                                  ;
  310.    Si (NON (EXIST Udebt))                                   ;
  311.      Ud = Udeb                                              ;
  312.      IK1 = IK1 + 1                                          ;
  313.    FinSi                                                    ;
  314.    Si (EXIST Udebt)                                         ;
  315.      t  = extr Udebt 'ABSC'                                 ;
  316.      u  = extr Udebt 'ORDO'                                 ;
  317.      Ud = 'IPOL' tps t  u                                   ;
  318.      IK1 = IK1 + 1                                          ;
  319.    FinSi                                                    ;
  320.  
  321.    Si (NON (EGA IK1 1 ))                                    ;
  322.      ERREUR 5                                               ;
  323.      QUITTER NSCLIM                                         ;
  324.    FinSi                                                    ;
  325.  FinSi                                                      ;
  326.  
  327.  
  328.  Si(EGA NOMQ 'DEBIT')                                       ;
  329.    Qd = Ud                                                  ;
  330.    Ud = Ud / Setr                                           ;
  331.    Si (EGA NOMG 'STRICTEMENT'); Ud = Ud * AT.'facgeo';FinSi ;
  332.  FinSi                                                      ;
  333.  Si(EGA NOMQ 'VITESSE')                                     ;
  334.    Qd = Ud * Setr                                           ;
  335.  FinSi                                                      ;
  336.  Si(EGA NOMQ 'TEMPERATURE')                                 ;
  337.    Qd = ' '                                                 ;
  338.  FinSi                                                      ;
  339.  
  340.  Si TIMPR                                                   ;
  341.  mess '**** Conditions limites ***** ' tit '  ' NOMQ ' ' NCO
  342.                               ' Ud = ' Ud  'Qd = ' Qd       ;
  343.  FinSi                                                      ;
  344.  
  345. ***************************************************************
  346. *  Imposition d'une fluctuation
  347.  
  348.  Si (EGA NOMT 'SGE')                                        ;
  349.  Si TIMPR                                                   ;
  350.  mess ' Modèle SGE : Imposition d une fluctuation :'
  351.  ' Ecart type ' sgm '  portée ' portee                      ;
  352.  FinSi                                                      ;
  353.  
  354. * Ne concerne que VITESSE et DEBIT
  355.  Si(NON ((EGA NOMQ 'VITESSE') ou (EGA NOMQ 'DEBIT')))       ;
  356.  Mess '*********************************************************' ;
  357.  Mess ' ERREUR   ERREUR   ERREUR   ERREUR   ERREUR   ERREUR     ' ;
  358.  Mess '                                                         ' ;
  359.  Mess 'L imposition d une fluctuation ne concerne que           ' ;
  360.  Mess ' VITESSE ou DEBIT                                        ' ;
  361.  Mess '*********************************************************' ;
  362.  erreur 21 ;
  363.  quitter NSCLIM                                             ;
  364.  Finsi ;
  365.  
  366. * La portee doit être > 0
  367.  Si (portee '&lt;EG' 0.)                                       ;
  368.  ERREUR 5                                                   ;
  369.  QUITTER NSCLIM                                             ;
  370.  FinSi                                                      ;
  371.  
  372.     sgp = sgm / 2.                                          ;
  373.  
  374.    Si (NON (EXIST AT 'unct1' ))                             ;
  375. * initialisation unct1,uncn1,unca1
  376.   unct2='ALEA' 'BANDES_TOURNANTES' Sentr 'EXPO'
  377.         'SIGMA' sgp 'MOYENNE' Ud 'LAMBDA' portee            ;
  378.   unct2=(unct2 - Ud)                                        ;
  379. *   trace unct2 Sentr TITR 'unct2'                          ;
  380.  
  381.   uncn2='ALEA' 'BANDES_TOURNANTES' Sentr 'EXPO'
  382.         'SIGMA' sgm 'MOYENNE' Ud 'LAMBDA' portee            ;
  383.  
  384. *   trace uncn2 Sentr TITR 'uncn2'                          ;
  385.  
  386.   unca2='ALEA' 'BANDES_TOURNANTES' Sentr 'EXPO'
  387.         'SIGMA' sgp  'MOYENNE' Ud 'LAMBDA' portee           ;
  388.   unca2=(unca2 - Ud)                                        ;
  389. *   trace unca2 Sentr TITR 'unca1'               ;
  390.  
  391.    AT.'unct2' = unct2                            ;
  392.    AT.'uncn2' = uncn2                            ;
  393.    AT.'unca2' = unca2                            ;
  394.    AT.'t1'    = 0.                               ;
  395.    AT.'t2'    = 0.1*Tps                          ;
  396.    FinSi                                           ;
  397.  
  398.     t1  = AT.'t1'                                ;
  399.     t2  = AT.'t2'                                ;
  400.     alf = (Tps - t1)/(t2 - t1)                   ;
  401.  
  402.    Si (ALF >EG 1.)                                 ;
  403.     alf = 0.                                     ;
  404.     AT.'unct1'=AT.'unct2'                        ;
  405.     AT.'uncn1'=AT.'uncn2'                        ;
  406.     AT.'unca1'=AT.'unca2'                        ;
  407.     t1 = Tps                                     ;
  408.     t2 = Tps + (portee / (abs Ud))               ;
  409.     AT.'t1'= t1                                  ;
  410.     AT.'t2'= t2                                  ;
  411.  
  412.   unct2='ALEA' 'BANDES_TOURNANTES' Sentr 'EXPO'
  413.         'SIGMA' sgp 'MOYENNE' Ud 'LAMBDA' portee ;
  414.   unct2=(unct2 - Ud);
  415. *   trace unct2 Sentr TITR 'unct2'               ;
  416.  
  417.   uncn2='ALEA' 'BANDES_TOURNANTES' Sentr 'EXPO'
  418.         'SIGMA' sgm 'MOYENNE' Ud 'LAMBDA' portee ;
  419.  
  420. *   trace uncn2 Sentr TITR 'uncn2'               ;
  421.  
  422.   unca2='ALEA' 'BANDES_TOURNANTES' Sentr 'EXPO'
  423.         'SIGMA' sgp 'MOYENNE' Ud 'LAMBDA' portee ;
  424.   unca2=(unca2 - Ud);
  425. *   trace unca2 Sentr TITR 'unca2'               ;
  426.  
  427.    AT.'unct2' = unct2                            ;
  428.    AT.'uncn2' = uncn2                            ;
  429.    AT.'unca2' = unca2                            ;
  430.  
  431.    FinSi ;
  432.  
  433.  unct1=AT.'unct1'                                ;
  434.  unct2=AT.'unct2'                                ;
  435.  uncn1=AT.'uncn1'                                ;
  436.  uncn2=AT.'uncn2'                                ;
  437.  unca1=AT.'unca1'                                ;
  438.  unca2=AT.'unca2'                                ;
  439.  
  440.  unct = (alf * unct1) + ((1. - alf) * unct2);
  441.  uncn = (alf * uncn1) + ((1. - alf) * uncn2);
  442.  unca = (alf * unca1) + ((1. - alf) * unca2);
  443.  
  444.    Si TIMPR                                            ;
  445.  tt=chai tit ' uncn' ' Tps=' Tps ' t1=' t1 ' t2=' t2   ;
  446.  trace uncn Sentr TITR tt                              ;
  447.  tt=chai tit ' unca'             ' Alf=' alf           ;
  448.  trace unca Sentr TITR tt                              ;
  449.  trace unct Sentr TITR (chai tit ' unct')              ;
  450.    FinSi                                               ;
  451.  Sinon                                                 ;
  452. ***************************************************************
  453.   uncn = Ud                                            ;
  454.   unct = 0.                                            ;
  455.   unca = 0.                                            ;
  456.  FinSi                                                 ;
  457. ***************************************************************
  458.  
  459.  
  460. Si((EGA NOMQ 'VITESSE') ou (EGA NOMQ 'DEBIT'))         ;
  461. NC1 = chai 1 NCO                                       ;
  462. NC2 = chai 2 NCO                                       ;
  463. NC3 = chai 3 NCO                                       ;
  464. LNC = MOTS NC1 NC2                                     ;
  465. UCxyzn=         (nomc NC1 (uncn*nx) 'NATU' 'DISCRET')
  466.              et (nomc NC2 (uncn*ny) 'NATU' 'DISCRET')  ;
  467.  
  468. UCxyzt=         (nomc NC1 (unct*tx) 'NATU' 'DISCRET')
  469.              et (nomc NC2 (unct*ty) 'NATU' 'DISCRET')  ;
  470.  
  471.  
  472.   Si(EGA (VALE DIME) 3)                                ;
  473. UCxyza=         (nomc NC1 (unca*ax) 'NATU' 'DISCRET')
  474.              et (nomc NC2 (unca*ay) 'NATU' 'DISCRET')  ;
  475.  
  476. UCxyzn=UCxyzn et (nomc NC3 (uncn*nz) 'NATU' 'DISCRET') ;
  477. UCxyzt=UCxyzt et (nomc NC3 (unct*tz) 'NATU' 'DISCRET') ;
  478. UCxyza=UCxyza et (nomc NC3 (unca*az) 'NATU' 'DISCRET') ;
  479. LNC = LNC et (MOTS NC3)                                ;
  480. Si SWIRL                                               ;
  481.  mess ' SWIRL : ' asw ' * Ud'                          ;
  482.   xyzw = extr uncn maillage                            ;
  483.   xw   = (coor 1 xyzw) - (coor 1 AT.'Cgrav') ;
  484.   yw   = (coor 2 xyzw) - (coor 2 AT.'Cgrav') ;
  485.   zw   = (coor 3 xyzw) - (coor 3 AT.'Cgrav') ;
  486.   GM   =        (nomc 'UX' xw 'NATU' 'DISCRET')
  487.              et (nomc 'UY' yw 'NATU' 'DISCRET')
  488.              et (nomc 'UZ' zw 'NATU' 'DISCRET')        ;
  489.   AM   =        (nomc 'UX' ax 'NATU' 'DISCRET')
  490.              et (nomc 'UY' ay 'NATU' 'DISCRET')
  491.              et (nomc 'UZ' az 'NATU' 'DISCRET')        ;
  492.   TM   =        (nomc 'UX' tx 'NATU' 'DISCRET')
  493.              et (nomc 'UY' ty 'NATU' 'DISCRET')
  494.              et (nomc 'UZ' tz 'NATU' 'DISCRET')        ;
  495.   xa = PSCA GM (MOTS 'UX' 'UY' 'UZ') AM (MOTS 'UX' 'UY' 'UZ');
  496.   ya = PSCA GM (MOTS 'UX' 'UY' 'UZ') TM (MOTS 'UX' 'UY' 'UZ');
  497.   R    = (xw*xw)+(yw*yw)+(zw*zw) **0.5                       ;
  498.   R0   = maxi R                                              ;
  499.   sw = (Ud/R0)*asw                                           ;
  500. UCxyza= UCxyza et (nomc NC1 (ya*ax*sw) 'NATU' 'DISCRET')
  501.                et (nomc NC2 (ya*ay*sw) 'NATU' 'DISCRET')
  502.                et (nomc NC3 (ya*az*sw) 'NATU' 'DISCRET')     ;
  503.  
  504. UCxyzt= UCxyzt et (nomc NC1 (xa*tx*(-1.)*sw) 'NATU' 'DISCRET')
  505.                et (nomc NC2 (xa*ty*(-1.)*sw) 'NATU' 'DISCRET')
  506.                et (nomc NC3 (xa*tz*(-1.)*sw) 'NATU' 'DISCRET') ;
  507.  
  508. FinSi                                                  ;
  509.   FinSi                                                ;
  510.  
  511.   UCxyz = UCxyzn et UCxyzt                             ;
  512.   Si(EGA (VALE DIME) 3)                                ;
  513.   UCxyz = UCxyz et UCxyza                              ;
  514.   FinSi                                                ;
  515. * UCxyz = (UCxyzn*0.) et UCxyzt et UCxyza;
  516.  
  517.  Si (EGA NOMG 'LARGEMENT')                             ;
  518. ***************************************************
  519. * NNC=chai AT.'MNC' NCO                                ;
  520. * UCc = EXCO (REDU UCxyzn (AT.'CSentr')) NNC NNC       ;
  521. * Sur la paroi seule la composante normale est retenue
  522. * Il semblerait plus simple et correct de faire comme suit
  523. ***************************************************
  524.   UCc =       REDU UCxyzn (AT.'CSentr')                ;
  525.   UCxyz=(REDU UCxyz (AT.'Sentro')) et UCc              ;
  526.  FinSi                                                 ;
  527.  
  528.   Si (EGA NOMG 'STRICTEMENT')                          ;
  529.   UCxyz= REDU UCxyz (AT.'Sentro')                      ;
  530.   FinSi                                                ;
  531.    Si TIMPR                                            ;
  532.    Si(EGA (VALE DIME) 3)                               ;
  533.    ung = vect UCxyz 0.05 NC1 NC2 NC3 jaune             ;
  534.    Sinon                                               ;
  535.    ung = vect UCxyz 0.05 NC1 NC2 jaune                 ;
  536.    FinSi                                               ;
  537.    trace ung (AT.'CSentr')      TITR tit               ;
  538.    FinSi                                               ;
  539.  
  540. FinSi                                                  ;
  541.  
  542. Si(EGA NOMQ 'TEMPERATURE')                             ;
  543. LNC = MOTS NCO                                         ;
  544. UCxyz= MANU 'CHPO' Sentr 'NATU' 'DISCRET' NCO uncn     ;
  545.   Si (EGA NOMG 'STRICTEMENT')                          ;
  546.   UCxyz = REDU UCxyz (AT.'Sentro')                     ;
  547.   FinSi                                                ;
  548.    Si TIMPR                                            ;
  549.    titt = chai tit ' composante : ' NCO                ;
  550.    trace UCxyz Sentr      TITR titt                    ;
  551.    FinSi                                               ;
  552. FinSi                                                  ;
  553.  
  554. Si(EGA NOMQ 'ADHERENCE')                               ;
  555. NC1 = chai 1 NCO                                       ;
  556. NC2 = chai 2 NCO                                       ;
  557. NC3 = chai 3 NCO                                       ;
  558. LNC = MOTS NC1 NC2                                     ;
  559. UCxyz= (MANU 'CHPO' Sentr 'NATU' 'DISCRET' NC1 0.)
  560.     et (MANU 'CHPO' Sentr 'NATU' 'DISCRET' NC2 0.)     ;
  561.   Si(EGA (VALE DIME) 3)                                ;
  562. UCxyz=UCxyz
  563. et (MANU 'CHPO' Sentr 'NATU' 'DISCRET' NC3 0.)         ;
  564. LNC = LNC et (MOTS NC3)                                ;
  565.   FinSi                                                ;
  566.   Si (EGA NOMG 'STRICTEMENT')                          ;
  567.   UCxyz = REDU UCxyz (AT.'Sentro')                     ;
  568.   FinSi                                                ;
  569. FinSi                                                  ;
  570.  
  571.  
  572.  Si(EXIST CLIM)                                        ;
  573.     unc=exco LNC (redu CLIM Sentr)      'NOID'         ;
  574.     unc = nomc LNC LNC unc ;
  575.     RESU = (CLIM - unc ) + UCxyz                       ;
  576.  Sinon                                                 ;
  577.     RESU = UCxyz                                       ;
  578.  FinSi                                                 ;
  579.  
  580.  Si TIMPR                                              ;
  581.  mess '#### FIN NSCLIM ###################################'
  582.       '#########################'                                      ;
  583.  FinSi                                                 ;
  584.  FINPROC RESU                                                          ;
  585.  
  586.  
  587.  
  588.  

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