Télécharger la plaquette Cast3M


fr en

Télécharger in_mini.procedur

Retour à la liste

Numérotation des lignes : 

  1. * IN_MINI   PROCEDUR  JC220346  14/12/09    21:15:08     8314           
  2.   'DEBPROC' IN_MINI  FCAN*'MAILLAGE'   HS*'CHPOINT' ORIG*'POINT'        
  3.     B_ANTI*'TABLE' ;                                                         
  4. ************************************************************************
  5. *   integration   de v  par minimisation  moindres carres               
  6. *   en entree   maillage frontiere cote phi et hs sur ce maiilage       
  7. *   en sortie   V(b)= phi(b)-psi(b)  chpoint sur fcan                   
  8. *    avec psi(p1) = phi(p1)                                             
  9. ************************************************************************
  10.   HX='EXCO' HS BX ;HY='EXCO' HS BY ; HZ='EXCO' HS BZ ;                  
  11.   nbi = 'NBNO' ( fcan 'ELEM' 1 ) ;                                      
  12.   I= 0 ;                                                                
  13.   'MESS' ' SEPARATION ' (FCAN NBNO)  'POINTS'  (FCAN NBEL) 'ELEM ';     
  14.   IMENA = 0 ;                                                           
  15.   MESS ' CALCUL DU RHS ' ;                                              
  16.   'REPE' BLOCALC (NBEL FCAN);                                           
  17.   imena = imena + 1 ;                                                   
  18.   I= I + 1 ;                                                            
  19.   SEGCOU = FCAN 'ELEM' I;                                               
  20.   P1 = SEGCOU 'POIN' 1;                                                 
  21.   P2 = SEGCOU 'POIN' 2;                                                 
  22.   HX1 = 'EXTR'  HX  'SCAL' P1 ; HX2 = 'EXTR'  HX 'SCAL'  P2 ;           
  23.   HY1 = 'EXTR'  HY  'SCAL' P1 ; HY2 = 'EXTR'  HY 'SCAL'  P2 ;           
  24.   HZ1 = 'EXTR'  HZ  'SCAL' P1 ; HZ2 = 'EXTR'  HZ 'SCAL'  P2 ;           
  25.   HMOY =(( HX1 +  HX2) / 2.) ( ( HY1 + HY2) / 2.) ((HZ1 + HZ2) / 2.);   
  26.  
  27.   VL= P2 'MOINS' P1 ;DL= 'NORM' VL ;VL=VL / DL ;                        
  28.   DV =  VL 'PSCA'  HMOY  ;DVI= DV * -1.;    
  29.   SEGCOU1 = CHAN 'POI1' SEGCOU ;
  30.   HH= 'MANU' 'CHPO' SEGCOU1 1 'Q'  ('PROG'  DVI DV ) NATURE DISCRET ;
  31.   'SI' ('EGA' I 1);                                                             
  32.   HTH  = HH ;                                                           
  33.   'SINON' ;                                                                     
  34.   HTH1= HTH 'ET' HH ;                                                           
  35.  'DETR' HTH ;HTH= HTH1 ;  'DETR' HH;                                    
  36.   'FINSI' ;                                                             
  37.    'SI' ('EGA' IMENA 100 ) ; 'MESS' ' menage' ;                         
  38.  'MENAGE'  ;imena = 0 ; 'FINSI' ;                                       
  39.    'FIN' BLOCALC ;                                                      
  40. ********                                                                
  41.     OB1 = FCAN MODE THERMIQUE ISOTROPE BARR ;                           
  42.     MAT1 = MATE OB1 'K' 1. ;                                            
  43.     CAR1 = CARA  OB1  'SECT'  1. ;                                      
  44.     RIGT = CONDUC OB1 (MAT1 ET CAR1 ) ;                                 
  45. ********                                                                
  46.   III = 'INDEX' B_ANTI ;INI = 'DIME' III ;                              
  47.      IZ = 0 ;                                                           
  48.             'REPETER' BENS  INI ;                                       
  49.    IZ = IZ + 1 ;STN = B_ANTI.IZ; 
  50.  
  51.    'SI' ('NON' ('EXISTE'  STN 'MTYP' )) ;
  52.  
  53.        'SI' ('EXISTE' STN  'IMPOSE' )  ;
  54.  'MESS' ' APPLICATION POTENTIEL IMPOSE  '  ;
  55.    TBB = 'BLOQUER' (EXTR (STN . 'IMPOSE') 'MAILLAGE' )  'T';                    
  56.    FFF = DEPIMP   TBB  (STN.'IMPOSE') ;                                   
  57.    HTH = HTH ET FFF ;
  58.    RIGT = RIGT ET TBB ;                                                   
  59.        'SINON'  ;    
  60. *  ON A PASSE  DIRECTEMENT UN OBJET DE TYPE RAIDEUR                             
  61.    RIGT = RIGT 'ET' ( STN.'COND'  ) ;                                   
  62.        'FINSI' ;     
  63.  
  64.    'SINON' ;
  65.  MMMM = STN.'MTYP' ;                                                    
  66.        'SI' ( 'EGA' MMMM  'TBLOQ' ) ;                                   
  67.    GEOG =STN.'LGEO' ;                                                    
  68.  'MESS' ' APPLICATION COND. LIM. CALCUL DU SAUT ' MMMM ;     
  69.    TBB =    'BLOQUER' 'T' GEOG   ; ;                                    
  70.    RIGT = RIGT 'ET' TBB                    ;                            
  71.  
  72.        'FINSI' ;                                                        
  73.        'SI' ( 'EGA' MMMM  'RENSE' ) ;                                   
  74.    GEOG =STN.'LGEO' ;                                                    
  75. 'MESS' ' APPLICATION COND. LIM. CALCUL DU SAUT ' MMMM ;    
  76.    RIGT = RIGT 'ET' ( 'RELA' 'ENSE' 'T' GEOG ) ;                        
  77.        'FINSI' ;
  78.  RIGT = RIGT 'ET' ('BLOQUER' 'T' ORIG) ;                
  79.    'FINSI' ;
  80.              'FIN' BENS ;                                              
  81.  
  82.   HHHH= 'RESOU'  RIGT   HTH ;                                           
  83.   HHHH= ('ENLEVER' HHHH LX ) 'NOMC' 'SCAL';                             
  84.   'FINPROC' HHHH ;                                                      
  85.  
  86.  
  87.  
  88.  

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