Télécharger la plaquette Cast3M


fr en

Télécharger mfil.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier mfil.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. ************************************************************************
  6. ** Filtering of a field using a rigidity matrix genrated by MFIL.
  7. **
  8. ** Author:
  9. ** Guenhael Le Quilliec (LaMe - Polytech Tours)
  10. **
  11. ** Version:
  12. ** 1.0 2021/01/04 Original version compatible with MFIL V1.0
  13. ** 1.1 2021/01/07 Version compatible with MFIL V1.1
  14. ************************************************************************
  15.  
  16. * General options
  17. OPTI 'DIME' 2 'MODE' 'PLAN' 'CONT' 'ELEM' QUA4 ;
  18. graph0 = FAUX ;
  19.  
  20. * Filtering radius
  21. r0 = 5.0 ;
  22.  
  23. * Mesh
  24. p0 = 0.0 0.0 ;
  25. p1 = 0.0 100.0 ;
  26. p2 = 100.0 0.0 ;
  27. lgn0 = DROI 100 p1 p0 ;
  28. msh0 = TRAN lgn0 100 p2 ;
  29.  
  30. * Model
  31. mod0 = MODE msh0 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' ;
  32. * Material
  33. mat0 = MATE mod0 'YOUN' 200.0e9 'NU' 0.3 'RHO' 1.0 ;
  34.  
  35. * Volume node field (as rho = 1)
  36. vol0 = EXCO 'FX' ((MASS mod0 mat0) * (MANU 'CHPO' msh0 1 'UX' 1.0)) ;
  37. * Rigidity matrix for filtering
  38. filter0 = MFIL vol0 r0 ;
  39.  
  40. * Create a node field to be filtered
  41. pnt0 = ((msh0 POIN 'DROI' (100.0 0.0) (60.0 100.0) 10.0) 
  42.               POIN 'DROI' (0.0  60.0) (100.0 60.0) 30.0)
  43.      ET (msh0 POIN 'SPHE' (40.0 30.0) (40.0 50.0) 5.0)
  44.      ET (msh0 POIN 'SPHE' (-50.0 150.0) (50.0 150.0) 15.0) ;
  45. field0 = (MANU 'CHPO' msh0 1 'SCAL' 0.0 'NATURE' 'DIFFUS')
  46.        + (MANU 'CHPO' pnt0 1 'SCAL' 1.0 'NATURE' 'DIFFUS') ;
  47.  
  48. * Filter the node field
  49. field1 = field0 * filter0 ;
  50.  
  51. * Plot the field to screen before filtering
  52. SI graph0 ;
  53.     TRAC field0 (msh0 CHAN 'POI1') ;
  54.     TRAC field0 msh0 ;
  55. FINS ;
  56.  
  57. * Plot the field to screen after filtering
  58. SI graph0 ;
  59.     TRAC field1 (msh0 CHAN 'POI1') ;
  60.     TRAC field1 msh0 ;
  61. FINS ;
  62.  
  63. FIN ;
  64.  
  65.  
  66.  
  67.  

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