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