Télécharger la plaquette Cast3M


fr en

Télécharger lapn.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  lapn.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ******************************************************************
  5. ******************************************************************
  6. * CALCUL DU LAPLACIAN EN VF                                      *
  7. * LAPL(T)=0                                                      *
  8. *                                                                *
  9. * On controlle que l'algorithme est lineaire exact               *
  10. * A. BECCANTINI LTMF    JUILLET    2001                          *
  11. ******************************************************************
  12. ******************************************************************
  13.  
  14.  'OPTION'  'DIME' 2 
  15.            'ISOV' 'SULI' 'TRAC' 'X'
  16.            'ELEM' 'QUA4' 'ECHO'  0 ;
  17.  
  18.  GRAPH = FAUX ;
  19. * GRAPH = VRAI ;
  20.  
  21. *
  22. **** Domaine
  23. *
  24.  
  25.  A1 = 0.0 0.0 ;
  26.  A2 = 1.0 0.0 ;
  27.  A3 = 1.0 1.0 ;
  28.  A4 = 0.0 1.0 ;
  29.  
  30.  A1A2 = A1 'DROIT' 4 A2 ;
  31.  A2A3 = A2 'DROIT' 5 A3 ;
  32.  A3A4 = A3 'DROIT' 6 A4 ;
  33.  A4A1 = A4 'DROIT' 5 A1 ;
  34.  
  35.  DOMCON = A1A2 'ET' A2A3 'ET' A3A4 'ET'
  36.           A4A1 ;
  37.  DOMTOT = 'SURF' DOMCON 'PLAN' ;
  38.  
  39. *
  40. *** L'objet modele Euler
  41. *
  42.  
  43.  
  44.  $DOMTOT = 'MODELISER' DOMTOT 'EULER';
  45.  
  46.  $DOMCON = 'MODELISER'  DOMCON 'EULER';
  47.  $A1A2 = 'MODELISER'  A1A2 'EULER';
  48.  $A2A3 = 'MODELISER'  A2A3 'EULER';
  49.  $A3A4 = 'MODELISER'  A3A4 'EULER';
  50.  $A4A1 = 'MODELISER'  A4A1 'EULER';
  51.  
  52.  TDOMTOT = 'DOMA'  $DOMTOT 'VF';
  53.  
  54.  TDOMCON = 'DOMA'  $DOMCON 'VF';
  55.  TA1A2 = 'DOMA'  $A1A2 'VF';
  56.  TA2A3 = 'DOMA'  $A2A3 'VF';
  57.  TA3A4 = 'DOMA'  $A3A4 'VF';
  58.  TA4A1 = 'DOMA'  $A4A1 'VF';
  59.  
  60.  MDOMCON = TDOMCON . 'QUAF' ;
  61.  MA1A2 = TA1A2 . 'QUAF' ;
  62.  MA2A3 = TA2A3 . 'QUAF' ;
  63.  MA3A4 = TA3A4 . 'QUAF' ;
  64.  MA4A1 = TA4A1 . 'QUAF' ;
  65.  
  66. **** old stuff  $DOMTOT = 'DOMA' DOMTOT ;
  67.  
  68.  MDOMTOT = TDOMTOT . 'QUAF' ;
  69.  
  70.  'ELIMINATION'  (MDOMTOT ET MDOMCON) 0.0001 ;
  71.  'ELIMINATION'  (MDOMTOT ET MA1A2) 0.0001 ;
  72.  'ELIMINATION'  (MDOMTOT ET MA2A3) 0.0001 ;
  73.  'ELIMINATION'  (MDOMTOT ET MA3A4) 0.0001 ;
  74.  'ELIMINATION'  (MDOMTOT ET MA4A1) 0.0001 ;
  75.  
  76.  
  77. ************************************
  78. ************************************
  79. *** Conditions initiales et aux ****
  80. *** limites                     ****
  81. ************************************
  82. ************************************
  83. *
  84. *
  85. *  A4     -      A3
  86. *
  87. *  
  88. *  |             
  89. *
  90. *  A1     -      A2
  91. *
  92.  
  93.  
  94. *
  95. * TEST1 
  96. * Champ lineaire
  97. * Conditions aux limites de type Dirichlet
  98. *
  99.  
  100.  XXLIM YYLIM = 'COORDONNEE' DOMCON ;
  101.  ACOEF = 1.0 ;
  102.  BCOEF = 7.0 ;
  103.  TLIM = (ACOEF '*' XXLIM) '+' (BCOEF '*' YYLIM) ;
  104.  
  105.  XXCEN YYCEN = 'COORDONNEE' ('DOMA' $DOMTOT 'CENTRE') ;
  106.  TN = (ACOEF '*' XXCEN) '+' (BCOEF '*' YYCEN) ;
  107.  
  108. *
  109. *** Graphique des c.i.
  110. *
  111.  
  112.  MOD1     =  'MODELISER'  ('DOMA' $DOMTOT 'MAILLAGE') 'THERMIQUE' ;
  113.  
  114.  'SI' GRAPH ;
  115.     CHM_TN  =  'KCHA' $DOMTOT 'CHAM' TN ;
  116.     'TRAC' CHM_TN  MOD1 'TITR' ('CHAINE' 'TN at t= ' 0.0);
  117.  'FINSI' ;
  118.  
  119.  
  120. *
  121. **** Solutions
  122. *
  123.  
  124.  GRADTN MCHAM = 'PENT'  $DOMTOT  'FACE' 'DIAMANT' 
  125.           TN 'CLIM' TLIM ;
  126.  
  127. * Le reste ne serve pas, mais il faut l'initialiser
  128.  
  129.  
  130.  MU = 0.0 ;
  131.  CV = 1.0 ;
  132.  KAPPA = 1.0 ;
  133.  
  134.  RN = 'MANUEL' 'CHPO'  ('DOMA' $DOMTOT 'CENTRE') 1 'SCAL' 1.0 ;
  135.  
  136.  VN = 'MANUEL' 'CHPO'  ('DOMA' $DOMTOT 'CENTRE') 2 'UX' 0.0 'UY' 0.0 ;
  137.  GRADVN = 'MANUEL' 'CHPO'  ('DOMA' $DOMTOT 'FACE') 4
  138.       'P1DX' 0.0 'P1DY' 0.0
  139.       'P2DX' 0.0 'P2DY' 0.0 ;
  140.  LINC = 'MOTS' 'RN' 'UX' 'UY' 'TN' ;
  141.  
  142. *
  143. **** Le calcul
  144. *
  145.  
  146.  IJACO IRESI1 DT = 'LAPN' 'VF' 'PROPCOST' 'RESI' 'EXPL'
  147.     $DOMTOT  MU KAPPA CV RN VN TN GRADVN  GRADTN
  148.       LINC  ;
  149.  
  150.  'SI' (('MAXIMUM' IRESI1 'ABS') '>' 1.0D-5) ;
  151.     'MESSAGE' 'Probleme 1' ;
  152.     'ERREUR' 5 ;
  153.  'FINSI' ;
  154.  
  155. *
  156. **** TEST 2
  157. *    Champ lineaire
  158. *    Condition limite de type Neumann
  159. *
  160.  
  161.  QLIM = (-1.0 '*' KAPPA) '*'
  162.   ('MANUEL' 'CHPO' ('DOMA' $DOMCON 'CENTRE') 2 'UX' ACOEF 'UY' BCOEF);
  163.  
  164.  GRADTN MCHAM = 'PENT'  $DOMTOT  'FACE' 'DIAMANT' 
  165.           TN  ;
  166.  
  167.  IJACO IRESI1 DT = 'LAPN' 'VF' 'PROPCOST' 'RESI' 'EXPL'
  168.     $DOMTOT  MU KAPPA CV RN VN TN GRADVN  GRADTN
  169.       LINC  'QIMP' QLIM ;
  170.  
  171.  'SI' (('MAXIMUM' IRESI1 'ABS') '>' 1.0D-5) ;
  172.     'MESSAGE' 'Probleme 1' ;
  173.     'ERREUR' 5 ;
  174.  'FINSI' ;
  175.  
  176.  'FIN' ;
  177.  
  178.  
  179.  
  180.  
  181.  
  182.  
  183.  
  184.  
  185.  
  186.  
  187.  
  188.  
  189.  
  190.  
  191.  
  192.  

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