Télécharger fornod.procedur

Retour à la liste

Numérotation des lignes :

  1. * FORNOD PROCEDUR WANG 96/09/17 21:16:05 2278
  2.  
  3. 'DEBPROC' FORNOD LIG1*'MAILLAGE' COM1*'MOT' FO*'FLOTTANT' ;
  4.  
  5. *---------------------------------------------------------------------
  6. *
  7. * POUR CALCULER LES FORCES NODALES EQUIVALENTES D'UNE FORCE
  8. * SURFACIQUE Fo
  9. * ELEMENT HORIZONTAL DE TYPE SEG3
  10. * MODE DE CALCUL : 'PLANDEFO' 'AXIS' ET 'FOUR 0, 1'
  11. *
  12. * LIG1 'MAIILAGE' ligne horizontale sur laquelle on
  13. * applique la force surfacique
  14. * COM 'MOT' nom de la composante de la force
  15. * PLAN DEFO : FX, FY
  16. * FOUR 1 : FR ,FZ ,FT
  17. *
  18. * FO 'FLOTTANT' valeur de la composante de la force
  19. * surfacique
  20. *
  21. *---------------------------------------------------------------------
  22.  
  23. TYP_MODE = 'MOT' ('VALEUR' 'MODE') ;
  24. 'SI' ('EGA' TYP_MODE 'PLANDEFO') ;
  25. ITYP = 1 ;
  26. 'SINON' ; 'SI' ('EGA' TYP_MODE 'FOUR') ;
  27. ITYP = 2 ;
  28. 'SINON' ;
  29. 'MESSAGE' ' ***** ERREUR *****' ;
  30. 'MESSAGE' ' On ne traite pas les modèles tridimensionnelles.'
  31. 'MESSAGE' ' Arrêt de la procédure.' ;
  32. 'QUITTER' BOUCERR ;
  33. 'FINSI' ; 'FINSI' ; 'FINSI' ;
  34. *
  35. * Boucle sur les éléments
  36. *
  37. NB_ELEM = 'NBEL' LIG1 ;
  38. I_EL = 0 ;
  39. 'REPETER' BOUCK NB_ELEM ;
  40. I_EL = I_EL + 1 ;
  41. ELEM_I = LIG1 'ELEMENT' I_EL ;
  42. ELEM_B = 'CHANGER' ELEM_I 'POI1' ;
  43. NOEU1 = ELEM_B 'POINT' 1 ;
  44. NOEU2 = ELEM_B 'POINT' 2 ;
  45. NOEU3 = ELEM_B 'POINT' 3 ;
  46. R1 = 'COOR' NOEU1 1 ;
  47. R3 = 'COOR' NOEU3 1 ;
  48. LONG = 'ABS' (R3 - R1) ;
  49.  
  50. *** MODE DEFO PLAN ***
  51.  
  52. 'SI' ('EGA' ITYP 1) ;
  53. F1 = LONG * FO / 6. ;
  54. F2 = F1 * 4. ;
  55. F3 = F1 ;
  56.  
  57. 'SI' ('EGA' I_EL 1) ;
  58. FOR1 = ( 'FORCE' COM1 F1 NOEU1 ) 'ET' ( 'FORCE' COM1 F2 NOEU2 )
  59. 'ET' ( 'FORCE' COM1 F3 NOEU3 ) ;
  60. 'SINON' ; FOR1 = FOR1 'ET' ( 'FORCE' COM1 F1 NOEU1 ) 'ET'
  61. ( 'FORCE' COM1 F2 NOEU2 ) ET ( 'FORCE' COM1 F3 NOEU3 ) ;
  62. 'FINSI' ;
  63. 'FINSI' ;
  64.  
  65. *** MODE AXIS ET FOURIER ***
  66.  
  67. 'SI' ( 'EGA' ITYP 2 ) ;
  68. F1 = PI * LONG * FO * R1 / 6. ;
  69. F2 = PI * LONG * FO * ( R1 + R3 ) / 3. ;
  70. F3 = PI * LONG * FO * R3 / 6. ;
  71. 'SI' ('EGA' I_EL 1) ;
  72. FOR1 = ( 'FORCE' COM1 F1 NOEU1 ) 'ET' ( 'FORCE' COM1 F2 NOEU2 )
  73. 'ET' ( 'FORCE' COM1 F3 NOEU3 ) ;
  74. 'SINON' ;
  75. FOR1 = FOR1 'ET' ( 'FORCE' COM1 F1 NOEU1 ) 'ET'
  76. ( 'FORCE' COM1 F2 NOEU2 ) 'ET' ( 'FORCE' COM1 F3 NOEU3 ) ;
  77. 'FINSI' ;
  78. 'FINSI' ;
  79. 'FIN' BOUCK ;
  80. 'FINPROC' FOR1 ;
  81. ***********************************************************************
  82.  
  83.  
  84.  

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