Procedures | Cast3M

   1 : * IN_MINI   PROCEDUR  BAZE      97/01/30    21:15:12     2468           
   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 :   HH= 'MANU' 'CHPO' SEGCOU 1 'Q'  ('PROG'  DVI DV ) NATURE DISCRET ;
  30 :   'SI' ('EGA' I 1);                                                             
  31 :   HTH  = HH ;                                                           
  32 :   'SINON' ;                                                                     
  33 :   HTH1= HTH 'ET' HH ;                                                           
  34 :  'DETR' HTH ;HTH= HTH1 ;  'DETR' HH;                                    
  35 :   'FINSI' ;                                                             
  36 :    'SI' ('EGA' IMENA 100 ) ; 'MESS' ' menage' ;                         
  37 :  'MENAGE'  ;imena = 0 ; 'FINSI' ;                                       
  38 :    'FIN' BLOCALC ;                                                      
  39 : ********                                                                
  40 :     OB1 = FCAN MODE THERMIQUE ISOTROPE BARR ;                           
  41 :     MAT1 = MATE OB1 'K' 1. ;                                            
  42 :     CAR1 = CARA  OB1  'SECT'  1. ;                                      
  43 :     RIGT = CONDUC OB1 (MAT1 ET CAR1 ) ;                                 
  44 : ********                                                                
  45 :   III = 'INDEX' B_ANTI ;INI = 'DIME' III ;                              
  46 :      IZ = 0 ;                                                           
  47 :             'REPETER' BENS  INI ;                                       
  48 :    IZ = IZ + 1 ;STN = B_ANTI.IZ; 
  49 :    
  50 :    'SI' ('NON' ('EXISTE'  STN 'MTYP' )) ;
  51 :    
  52 :        'SI' ('EXISTE' STN  'IMPOSE' )  ;
  53 :  'MESS' ' APPLICATION POTENTIEL IMPOSE  '  ;
  54 :    TBB = 'BLOQUER' (EXTR (STN . 'IMPOSE') 'MAILLAGE' )  'T';                    
  55 :    FFF = DEPIMP   TBB  (STN.'IMPOSE') ;                                   
  56 :    HTH = HTH ET FFF ;
  57 :    RIGT = RIGT ET TBB ;                                                   
  58 :        'SINON'  ;    
  59 : *  ON A PASSE  DIRECTEMENT UN OBJET DE TYPE RAIDEUR                             
  60 :    RIGT = RIGT 'ET' ( STN.'COND'  ) ;                                   
  61 :        'FINSI' ;     
  62 :    
  63 :    'SINON' ;
  64 :  MMMM = STN.'MTYP' ;                                                    
  65 :        'SI' ( 'EGA' MMMM  'TBLOQ' ) ;                                   
  66 :    GEOG =STN.'LGEO' ;                                                    
  67 :  'MESS' ' APPLICATION COND. LIM. CALCUL DU SAUT ' MMMM ;     
  68 :    TBB =    'BLOQUER' 'T' GEOG   ; ;                                    
  69 :    RIGT = RIGT 'ET' TBB                    ;                            
  70 :                                                      
  71 :        'FINSI' ;                                                        
  72 :        'SI' ( 'EGA' MMMM  'RENSE' ) ;                                   
  73 :    GEOG =STN.'LGEO' ;                                                    
  74 : 'MESS' ' APPLICATION COND. LIM. CALCUL DU SAUT ' MMMM ;    
  75 :    RIGT = RIGT 'ET' ( 'RELA' 'ENSE' 'T' GEOG ) ;                        
  76 :        'FINSI' ;
  77 :  RIGT = RIGT 'ET' ('BLOQUER' 'T' ORIG) ;                
  78 :    'FINSI' ;
  79 :              'FIN' BENS ;                                              
  80 :                              
  81 :   HHHH= 'RESOU'  RIGT   HTH ;                                           
  82 :   HHHH= ('ENLEVER' HHHH LX ) 'NOMC' 'SCAL';                             
  83 :   'FINPROC' HHHH ;                                                      
  84 :                                                       
  85 :