Télécharger topoptim_03.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier topoptim_03.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. ************************************************************************
  6. ** Topology optimization of a simple 2D structure subjected to a
  7. ** thermal loading.
  8. **
  9. ** Author:
  10. ** Guenhael Le Quilliec (LaMe - Polytech Tours)
  11. **
  12. ** Version:
  13. ** V2.0 2017/04/18 Updated to make it compatible with TOPOPTIM V2.0
  14. ** V1.0 2014/04/28 Original version compatible with TOPOPTIM V1.0
  15. ************************************************************************
  16.  
  17. graph0 = FAUX ;
  18.  
  19. * General options
  20. OPTI 'DIME' 2 'ELEM' QUA4 ;
  21.  
  22. * Mesh
  23. p0 = 0.0 0.0 ;
  24. p1 = 0.0 40.0 ;
  25. p2 = 40.0 0.0 ;
  26. lgn0 = DROI 40 p0 p1 ;
  27. msh0 = TRAN lgn0 40 p2 ;
  28. p2 = msh0 POIN 'PROC' p2 ;
  29. lgn1 = lgn0 PLUS p2 ;
  30. ELIM (CONT msh0) lgn1 0.01 ;
  31. lgn2 = DROI 9 (0.0 18.0) (0.0 22.0) ;
  32. ELIM lgn0 lgn2 0.01 ;
  33.  
  34. * Model and material
  35. mod0 = MODE msh0 'THERMIQUE' 'CONDUCTION' ;
  36. mat0 = MATE mod0 'K' 50.0 ;
  37.  
  38. * Boundary conditions and loading
  39. bc0 = BLOQ 'T' lgn2 ;
  40. load0 = SOUR mod0 0.01 msh0 ;
  41.  
  42. * Finite element model table
  43. mdl0 = TABL ;
  44. mdl0.'MODELE' = mod0 ;
  45. mdl0.'CARACTERISTIQUES' = mat0 ;
  46. mdl0.'BLOCAGES_THERMIQUES' = bc0 ;
  47. mdl0.'CHARGEMENT' = load0 ;
  48.  
  49. * Optimization table
  50. tab0 = TABL ;
  51. tab0.'RESOLUTION_LINEAIRE' = mdl0 ;
  52. tab0.'RAPPORT_RAIDEURS_THERMIQUES' = 26.0e-3 / 50.0 ;
  53. tab0.'TRAC' = graph0 ;
  54.  
  55. * Optimization
  56. TOPOPTIM tab0 ;
  57.  
  58. * Plot to screen
  59. topo0 = tab0.'TOPOLOGIE'.(tab0.'CYCLE') ;
  60. topomsh0 = tab0.'MAILLAGE'.(tab0.'CYCLE') ;
  61. SI graph0 ;
  62. TRAC (REDU topo0 topomsh0) (REDU mod0 topomsh0)
  63. (PROG 0.0 'PAS' (1.0 / 56.0) 1.0) ;
  64. FINS ;
  65.  
  66. FIN ;
  67.  
  68.  
  69.  
  70.  

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