Télécharger la plaquette Cast3M


fr en

Télécharger @deduire.procedur

Retour à la liste

Numérotation des lignes : 

  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.  
  69.  

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