Télécharger la plaquette Cast3M


fr en

Télécharger visucoq.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : visucoq.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ************************************************************************
  5. ************************************************************************
  6. *                                                                      *
  7. *             Exemple d'utilisation de coq2mas                         *
  8. *         Visualisation 3D de résultats de calcul coque multicouche    *
  9. *                                                                      *
  10. *                      D. Combescure - Juillet 2007                    *
  11. *                      Laboratoire DYN - CEA Saclay                    *
  12. ************************************************************************
  13. ************************************************************************
  14. *
  15. *FLAGVISU = VRAI;
  16. *FLAGLONG = VRAI;
  17. *
  18. FLAGVISU = FAUX;
  19. *FLAGLONG = FAUX;
  20. *
  21. opti dime 3 elem tri3;
  22. R = 10.;
  23. H = -1.;
  24. ep1 = 1.0;
  25. ep2 = 0.20;
  26. p1 = R 0. 0.;
  27. p2 = 0. R 0.;
  28. p3 = ((-1.)*R) 0. 0.;
  29. p4 = 0. ((-1.)*R) 0.;
  30. p0 = 0. 0. H;
  31. p0b = 0. 0. (-2.*H);
  32. *
  33. n1 = 20;
  34. lig0 = (CER3 n1 p1 p2 p3) 
  35.     ET (CER3 n1 p3 p4 p1); 
  36. sur0 = SURF lig0 SPHE p0;
  37. sur0b = SURF lig0 SPHE p0b;
  38. *
  39. mod1 = MODE sur0 mecanique elastique dkt cons 'Couche 1';
  40. mat1  = MATE mod1 YOUN 30000.D6 NU 0.25 RHO 2400.
  41.                   EPAI ep1 EXCE (1.*ep1);
  42. mod2 = MODE sur0 mecanique elastique dkt cons 'Couche 2';
  43. mat2  = MATE mod2 YOUN 30000.D6 NU 0.25 RHO 2400.
  44.                   EPAI ep1 EXCE ((-1.)*ep1);                   
  45. mod3 = MODE sur0 mecanique elastique dkt cons 'Couche 3';
  46. mat3  = MATE mod3 YOUN 30000.D6 NU 0.25 RHO 2400.
  47.                   EPAI ep1 ; 
  48. mod4 = MODE sur0b mecanique elastique dkt cons 'Couche 4';
  49. mat4  = MATE mod4 YOUN 30000.D6 NU 0.25 RHO 2400.
  50.                   EPAI ep2 ;                   
  51. modtot = mod1 et mod2 et mod3 et mod4;
  52. mattot = mat1 et mat2 et mat3 et mat4;
  53. rigtot = rigi modtot mattot;
  54. mastot = mass modtot mattot;
  55. *
  56. tabtot = table;
  57. mesh3D = coq2mas modtot mattot tabtot;
  58. *   
  59. SI FLAGVISU; 
  60.  trac (sur0 et sur0b) cach;
  61.  trac mesh3D cach; 
  62. FINSI;
  63. * 
  64. load = mastot * (manu chpo (sur0 et sur0b) 1 UZ -9.81);
  65. bl0 = bloq depl lig0;
  66. *
  67. *
  68. DEPJ = RESO (RIGTOT et BL0) load;
  69. CONJ = SIGMA DEPJ MODTOT MATTOT; 
  70. *
  71. TABTOT. 'DEPLACEMENTS' = TABLE;
  72. TABTOT. 'DEPLACEMENTS' . 1 = DEPJ;
  73. TABTOT. 'CONTRAINTES' = TABLE;
  74. TABTOT. 'CONTRAINTES' . 1 = CONJ;
  75. *
  76. MESH3D = COQ2MAS MODTOT MATTOT TABTOT;
  77. *
  78. SI FLAGVISU; 
  79.  dep3D = (TABTOT.'DEPLACEMENTS_VOLUMIQUE'. 1 . TOTAL);
  80.  trac (defo mesh3D dep3D) cach;
  81.  mesh3Df = TABTOT.'MAILLAGE_FIBRE_MOYENNE' .  TOTAL;
  82.  con3Df =  EXCO (TABTOT.'CONTRAINTES_FIBRE_MOYENNE'. 1 . TOTAL)
  83.           SMXX;
  84.  trac mesh3Df con3Df cach;
  85.  mesh3D = TABTOT.'MAILLAGE_VOLUMIQUE' .  TOTAL;
  86.  con3D =  EXCO (TABTOT.'CONTRAINTES_VOLUMIQUE'. 1 . TOTAL)
  87.           SMXX;
  88.  trac mesh3D con3D cach;
  89. * Pour relecture avec PARAVIEW
  90.  opti sort 'visucoq.inp';
  91.  sort avs mesh3D (TABTOT.'CONTRAINTES_VOLUMIQUE'. 1 . TOTAL);
  92. FINSI;
  93. *
  94. FIN;
  95. *
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  

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