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