Télécharger topoptim2.dgibi

Retour à la liste

Numérotation des lignes :

  1. ************************************************************************
  2. ************************************************************************
  3. ************************************************************************
  4. ** Topology optimization of a simple 2D structure subjected to a
  5. ** mechanical loading, with penalty factor and GSF.
  6. **
  7. ** Guenhael Le Quilliec (LMR - Polytech Tours)
  8. ** 2014/04/28
  9. ************************************************************************
  10.  
  11. * General options
  12. OPTI 'DIME' 2 'MODE' 'PLAN' 'CONT' 'ELEM' QUA4 ;
  13. OPTI 'TRAC' 'PSC';
  14.  
  15. * Mesh
  16. nelx0 = 90 ;
  17. nely0 = 30 ;
  18. p0 = 0.0 0.0 ;
  19. p1 = 0.0 (FLOT nely0) ;
  20. p2 = (FLOT nelx0) 0.0 ;
  21. lgn0 = DROI nely0 p0 p1 ;
  22. msh0 = TRAN lgn0 nelx0 p2 ;
  23. p2 = msh0 POIN 'PROC' p2 ;
  24.  
  25. * Boundary conditions and loading
  26. bc0 = (BLOQ 'UX' lgn0) ET (BLOQ 'UY' p2) ;
  27. load0 = FORC (0.0 -1.0) p1 ;
  28.  
  29. * Optimization table
  30. tab0 = TABL ;
  31. tab0.'MAILLAGE' = msh0 ;
  32. tab0.'BLOCAGES_MECANIQUES' = bc0 ;
  33. tab0.'CHARGEMENT' = load0 ;
  34. tab0.'TAUX_FILTRAGE' = 10 ;
  35.  
  36. * Start value of P factor
  37. tab0.'P_FACTEUR' = 1.0 ;
  38. * Num of iters to hold P constant from start
  39. tab0.'P_STABLE' = 5 ;
  40. * Increment by this amount
  41. tab0.'P_INCREMENT' = 0.2 ;
  42. * Increment every 'P_NITER' iters
  43. tab0.'P_NITER' = 1 ;
  44. * Max value of 'P_NITER'
  45. tab0.'P_MAX' = 3.0 ;
  46.  
  47. * Start value of GSF
  48. tab0.'GSF_FACTEUR' = 1.0 ;
  49. * Num of iters to hold the GSF constant from start
  50. tab0.'GSF_STABLE' = 25 ;
  51. * Increment by this amount
  52. tab0.'GSF_INCREMENT' = 0.05 ;
  53. * Increment every 'GSF_NITER' iters
  54. tab0.'GSF_NITER' = 1 ;
  55. * Max value of 'GSF_NITER'
  56. tab0.'GSF_MAX' = 5.0 ;
  57.  
  58. * Optimization
  59. TOPOPTIM tab0 ;
  60.  
  61. * Plot to screen
  62. TRAC (REDU tab0.'TOPOLOGIE_CH' tab0.'TOPOLOGIE_MAIL')
  63. (REDU tab0.'MODELE' tab0.'TOPOLOGIE_MAIL')
  64. (PROG 0.0 'PAS' (1.0 / 56.0) 1.0) ;
  65.  
  66. FIN;
  67.  
  68.  
  69.  
  70.  
  71.  

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