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