Télécharger vitetfor.procedur

Retour à la liste

Numérotation des lignes :

  1. * VITETFOR PROCEDUR AM 96/09/09 21:15:09 2269
  2. 'DEBPROC' VITETFOR MTA*'TABLE' ETAB12*'TABLE'
  3. FCORU*'CHPOINT' ;
  4. *
  5. * MTA table de la rigidite associee au contact
  6. *
  7. FLAG = VRAI ;
  8. MAN = ETAB12.'FMAN' ;
  9. MAN = 'CHAN' 'ATTRIBUT' MAN 'NATURE' 'DISCRET' ;
  10. *
  11. MDI = 'DIME' MTA;
  12. 'SI' ( 'NEG' MDI 0) ;
  13. A = MTA.7 ;
  14. LISEA = MTA.6;
  15.  
  16.  
  17. MALI = 'MAXI' LISEA;
  18. 'SI' ( 'EGA' MALI 1 ) ;
  19. GEO1 = 'EXTR' A 'MAIL' 'MULT' ;
  20. *
  21. * on recupere quelques quantites necessaires
  22. *
  23. H = ETAB12.'DT' ;
  24. H2 = H * H ;
  25.  
  26. *------
  27. *
  28. * calcul de l'operateur
  29. *
  30. ATOT = ETAB12.'MASSE' 'ET' A ;
  31. *
  32. * calcul du 1-er second membre et de COE1
  33. *
  34. FU = ( -0.25 * H2 ) * MAN;
  35.  
  36. VDT = ( H * ETAB12.'VITI') + ETAB12.'DEPI' ;
  37. VDT = 'CHAN' 'ATTRIBUT' VDT 'NATURE' 'DIFFUS' ;
  38. FL = 'EXCO' ( A * VDT) 'FLX' 'FLX' ;
  39. FL = FL - ( 'REDU' FCORU GEO1 ) ;
  40. FT = FL 'ET' FU ;
  41. COE1 = 'RESO' ATOT FT 'NOUNIL' ;
  42. *
  43. * calcul du 2-eme second membre et de COE2
  44. *
  45. FL = 'EXCO' ( A * ETAB12.'VITI' ) 'FLX' 'FLX' ;
  46. FU = ( -0.5 * H ) * MAN;
  47. FT = FL 'ET' FU ;
  48. COE2 = 'RESO' ATOT FT 'NOUNIL' ;
  49. *
  50. * calcul des vecteurs VADD et FADD
  51. *
  52. DTL = (( 2. / H ) * COE1 ) - COE2 ;
  53. DTLL = 'REAC' A DTL ;
  54. VADD = 'RESO' ETAB12.'MASSE' DTLL ;
  55. FADD = (-4./ H2 ) * ( 'REAC' A COE1) ;
  56. 'SI' ( 'EXIS' ETAB12 'AMORTISSEMENT' ) ;
  57. FADD = FADD - (ETAB12.'AMORTISSEMENT' * VADD );
  58. 'FINSI' ;
  59. FLAG = FAUX ;
  60. 'FINSI';
  61. 'FINSI';
  62. *
  63. 'SI' FLAG ;
  64. VADD = 0. * ETAB12.'VITI' ;
  65. FADD = 0. * MAN ;
  66. 'FINSI' ;
  67. VADD = 'ENLE' VADD 'LX' ;
  68. *
  69.  
  70. 'FINPROC' VADD FADD ;
  71.  

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