Procedures | Cast3M

   1 : * @DEDUIRE  PROCEDUR  GOUNAND   05/10/10    21:15:01     5199           
   2 : debproc @deduire geo1*'MAILLAGE'   maianc*'MAILLAGE' 
   3 : mainouv*'MAILLAGE';
   4 : iidim= 'VALE' dime; 
   5 : 'SI' ( 'EGA' iidim 3) ;                                                     
   6 : xs ys zs = 'COOR' mainouv;                                                    
   7 : xv yv zv = 'COOR' maianc; 
   8 : 'SINO';
   9 : xs ys  = 'COOR' mainouv;                                                    
  10 : xv yv  = 'COOR' maianc; 
  11 : 'FINS';                                                 
  12 : mais = 'EXTR' xs 'MAILLAGE';                                                    
  13 : maiv = 'EXTR' xv 'MAILLAGE';                                                
  14 : nbp = 'NBEL' mais;   
  15 : 'SI' ( 'EGA' iidim 3);                                                          
  16 :  'REPETER' bouc3 nbp;i = 'POIN' &bouc3 mais;j = 'POIN' &bouc3 maiv;             
  17 :  xi yi zi = 'COOR' i;xj yj zj = 'COOR' j;tx = xi - xj;                      
  18 :  ty = yi - yj;tz = zi -zj; 
  19 :  mx = 'MANU' 'CHPO' j 1 'T' tx 'NATURE' 'DIFFUS';     
  20 :  my = 'MANU' 'CHPO' j 1 'T' ty 'NATURE' 'DIFFUS';                              
  21 :  mz = 'MANU' 'CHPO' j 1 'T' tz  'NATURE' 'DIFFUS';                              
  22 :  'SI' ( &bouc3 'EGA' 1);                                                     
  23 :   mtx = mx; mty = my ;mtz= mz; 'SINO';                                     
  24 :   mtx = mtx 'ET' mx; mty = mty 'ET' my; mtz = mtz 'ET' mz; 
  25 :  'FINS';               
  26 :  'FIN' bouc3; 
  27 : 'SINO';
  28 :  'REPETER' bouc2 nbp;i = 'POIN' &bouc2 mais;j ='POIN' &bouc2 maiv;  
  29 :  xi yi = 'COOR' i;xj yj  = 'COOR' j;tx = xi - xj;     
  30 :  ty = yi - yj;mx = 'MANU' 'CHPO' j 1 'T' tx 'NATURE' 'DIFFUS'; 
  31 :  my = 'MANU' 'CHPO' j 1 'T' ty 'NATURE' 'DIFFUS'; 
  32 :  'SI' ( &bouc2 'EGA' 1);                                                     
  33 :   mtx = mx; mty = my ; 'SINO' ;
  34 :   mtx = mtx 'ET' mx; mty = mty 'ET' my;
  35 :  'FINS';               
  36 :  'FIN' bouc2; 
  37 : 'FINS';
  38 : 
  39 : * il n'y a plus qu'a faire la resolution thermique                              
  40 : mot = 'MODE' geo1 'THERMIQUE' ;                                                 
  41 : mat = mot 'MATE' 'K' 2;                                                   
  42 : mth = 'CONDU' mot mat;                                               
  43 : bl = 'BLOQ' 'T' maianc;                                                         
  44 : mthb = mth 'ET' bl; 
  45 : 
  46 : * resolution pour les deplacements en X
  47 : ft = 'DEPI' bl mtx;                                                       
  48 : dx = 'RESO' mthb ft;                                                     
  49 : * resolution en Y                                                               
  50 : fty = 'DEPI' bl mty;                                                      
  51 : dy = 'RESO' mthb fty; 
  52 : * on additionne les deux chpoints 
  53 : ddx = ( 'EXCO' dx 'T') 'NOMC' 'UX';                                         
  54 : ddy = ( 'EXCO' dy 'T') 'NOMC' 'UY'; 
  55 : dep = ddx 'ET' ddy  ;
  56 : 
  57 : 'SI' ( iidim 'EGA' 3);                                               
  58 : ftz = 'DEPI' bl mtz;                                                      
  59 : dz = 'RESO' mthb ftz;                                                           
  60 : ddz = ( 'EXCO' dz 'T') 'NOMC' 'UZ';                                         
  61 : dep = dep 'ET' ddz;  
  62 : 'FINS';                                              
  63 : res = geo1 'PLUS' dep;                                                    
  64 : 'FINP' res;                                                 
  65 : 
  66 : 
  67 :  
  68 :