Télécharger topoptim2.dgibi

Retour à la liste

Numérotation des lignes :

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

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