Télécharger la plaquette Cast3M


fr en

Télécharger visufour1.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : visufour1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ************************************************************************
  5. ************************************************************************
  6. *                                                                      *
  7. *             Exemple d'utilisation de four2tri                        *
  8. *         Visualisation 3D de résultats de calcul Fourier              *
  9. *               avec possibilité de recombinaison des harmoniques      *
  10. *                                                                      *
  11. *                      D. Combescure - Janvier 2007                    *
  12. *                      Laboratoire DYN - CEA Saclay                    *
  13. ************************************************************************
  14. ************************************************************************
  15. *
  16. * FLAGVISU = VRAI; OPTI TRAC PSC;
  17. FLAGVISU = FAUX;
  18. *
  19. opti dime 2 elem tri3 mode four 'NOHA';
  20. *
  21. ************************************************************************
  22. * MAILLAGE 
  23. ************************************************************************
  24.  
  25. p1 = 10 0.01;
  26. p2 = 11 0.01;
  27. p3 = 20 0.01;
  28. p4 = 30 0.01;
  29. p5 = 30 1;
  30. v1 = 0 1;
  31. *
  32. lig1 = d 3 p1 p2;
  33. mesh2d = tran 3 lig1 v1;
  34. opti elem seg2;
  35. lig2 = d 3 p2 p3;
  36.  
  37. ************************************************************************
  38. * MODELE
  39. ************************************************************************
  40.  
  41. mod1 =  modeli mesh2d mecanique elastique;
  42. mod2 =  modeli lig2 mecanique elastique coq2;
  43. mod3 =  modeli (p5 et p4) mecanique elastique cerc;
  44. modt = mod1 et mod2;
  45. mesht = mesh2d et lig2;
  46.  
  47. * >>>  GENERATION DU MAILLAGE 3D SEULEMENT  <<<
  48. tab1 = TABLE;
  49. mesh3d = FOUR2TRI modt tab1 (PROG 0. pas 10. 360.);
  50. SI FLAGVISU; TRAC mesh3d; FINSI;
  51.  
  52.  
  53. ************************************************************************
  54. * CHPO  2D Fourier -> 3D
  55. ************************************************************************
  56. * listreel pour les isovaleurs
  57. izo1 = prog -1. PAS 0.1 1.;
  58. izo2 = 0.001 * izo1;
  59.  
  60.  
  61. * creation des tables de chpoint
  62. U2D = tabl; U2DZ = tabl;
  63. F2D = tabl; F2DZ = tabl;
  64. PRE1 = tabl;
  65.  
  66. REPE bnf 3; nf = &bnf - 1;
  67.  
  68.   opti dime 2 elem tri3 mode four nf;
  69.   U2D  . nf = manu chpo mesht  3 'UR' 1. 'UT' -1. 'UZ' 0. ;
  70.   U2DZ . nf = manu chpo mesht  3 'UR' 0. 'UT'  0. 'UZ' 1. ;
  71.   F2D  . nf = manu chpo mesht  2 'FR' 0.001 'FT' 0.001 ;
  72.   F2DZ . nf = manu chpo mesht  1 'FZ' 1. ;
  73.   PRE1 . nf = manu chpo mesht  2 'P' 0.001 'EPAI' 0.001 ;
  74.  
  75. FIN bnf;
  76.  
  77. opti mode four 'NOHA';
  78.  
  79. * >>> FOUR2TRI <<<
  80.  
  81. REPE bnf 3; nf = &bnf - 1;
  82.  
  83.   tab1 = TABLE;
  84.   tab1.'DEPLACEMENTS' = TABLE;
  85.   tab1.'DEPLACEMENTS'. 1 = U2D . nf;
  86.   tab1.'DEPLACEMENTS'. 2 = U2DZ . nf;
  87.   mesh3d = FOUR2TRI (mod1 et mod2) tab1 nf ;
  88.   dep3D1 = (tab1.deplacements_3d. 1);
  89.   dep3D2 = (tab1.deplacements_3d. 2);
  90.   mo3d = mots 'UX' 'UY' 'UZ';
  91.   SI FLAGVISU;
  92.     TITR (chai 'n_{H}='nf ' UR=-UT=+1');
  93.     u1 = psca dep3D1 dep3D1 mo3d mo3d;
  94.     trac (0 1E4 -1E4) (defo mesh3D dep3D1 u1) cach izo1;
  95.     TITR (chai 'n_{H}='nf ' UZ=+1');
  96.     u2 = psca dep3D2 dep3D2 mo3d mo3d;
  97.     trac (0 1E4 -1E4) (defo mesh3D dep3D2 u2) cach izo1;
  98.   FINSI;
  99.  
  100.   OUBL tab1 'DEPLACEMENTS';
  101.   tab1.'CHAMPS_POINTS_SYMETRIQUES' = TABLE;
  102.   tab1.'CHAMPS_POINTS_SYMETRIQUES'. 1 = PRE1 . nf;
  103.   mesh3d = FOUR2TRI (mod1 et mod2) tab1 nf ;
  104.   SI FLAGVISU;
  105.     trac mesh3d (exco (tab1.'CHAMPS_POINTS_SYME_3D'. 1) 'P')
  106.     TITR (chai 'n_{H}='nf ' p=+0.001') izo2;
  107.   FINSI;
  108.  
  109. FIN bnf;
  110.  
  111.  
  112. ************************************************************************
  113. * Combinaison d'harmoniques
  114. ************************************************************************
  115. *
  116. * creation d'une force composee de plusieurs harmoniques
  117. opti mode four 1;
  118. FPONC = manu chpo p3 1 FR 1. nature discrete;
  119. opti mode four 2;
  120. FPONC = FPONC et (manu chpo p3 1 FR 1. nature discrete);
  121. opti mode four 3;
  122. FPONC = FPONC et (manu chpo p3 1 FR 1. nature discrete);
  123. opti mode four 4;
  124. FPONC = FPONC et (manu chpo p3 1 FR 1. nature discrete);
  125. opti mode four 5;
  126. FPONC = FPONC et (manu chpo p3 1 FR 1. nature discrete);
  127. opti mode four 6;
  128. FPONC = FPONC et (manu chpo p3 1 FR 1. nature discrete);
  129. opti mode four 'NOHA';
  130.  
  131. * >>> FOUR2TRI <<<
  132.  
  133. tab1 = table;
  134. tab1.'EFFORTS' = TABLE;
  135. tab1.'EFFORTS'. 1 = FPONC;
  136.  
  137. F3D = VIDE 'CHPOINT';
  138. REPE bnf 6; nf = &bnf ;
  139.   mesh3d = four2tri (mod1 et mod2) tab1 nf (prog 0. pas 10. 360.);
  140.   F3D = F3D + tab1. 'EFFORTS_3D'. 1;
  141.   FVECT = vect F3D 'FX' 'FY' 'FZ' ROSE;
  142.   SI FLAGVISU;
  143.     trac (0 -1E4 0) mesh3d FVECT 
  144.     TITR (chai 'F^{3D}=\S_{n=1..'nf'} F_{n}');
  145.   FINSI;
  146. FIN bnf;
  147.  
  148. *
  149. TEMP 'IMPR' 'MAXI' 'CPU';
  150.  
  151. FIN ;
  152. *
  153.  
  154.  
  155.  
  156.  
  157.  
  158.  
  159.  
  160.  
  161.  
  162.  

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