Télécharger reso_asy.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : reso_asy.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * CAS TEST DU 91/06/18 PROVENANCE : PLAF
  5. * si GRAPH = N, les graphiques ne sont pas affich{s
  6. * si GRAPH diff{rent de N, tous les graphiques sont affich{s
  7. *
  8. GRAPH = 'N' ;
  9. *
  10. SAUT PAGE ;
  11. *--------------------------------------------------------------------*
  12. * RESO_ASY *
  13. *--------------------------------------------------------------------*
  14. *
  15. OPTION ECHO 0 ;
  16. TEMPS ;
  17. *
  18. OPTI DIME 2 ELEM SEG2 ;
  19. *
  20. P1 = 0 0 ;
  21. P2 = 1 0 ;
  22. M1 = D 1 P1 P2 ;
  23. *
  24. L_MOT = MOTS 'UX' 'UY' ;
  25. *
  26. * Matrice antisymétrique
  27. *
  28. M_A = MANU RIGI type 'MASSE' M1 L_MOT 'ANTI' ( PROG 0 -2 -1 0
  29. 2 0 -2 -1
  30. 1 2 0 -2
  31. 0 1 2 0 ) ;
  32. *
  33. * Matrice symétrique
  34. *
  35. M_S = MANU RIGI type 'MASSE' M1 L_MOT ( PROG 100 0 0 0
  36. 0 200 0 0
  37. 0 0 300 0
  38. 0 0 0 400 ) ;
  39. *
  40. * Matrice totale
  41. *
  42. M_M = M_S ET M_A ;
  43. *
  44. * Vecteur force donn{
  45. *
  46. CHP_F = ( MANU CHPO P1 2 'FX' 1. 'FY' 2. 'NATURE' 'DISCRET') ET
  47. ( MANU CHPO P2 2 'FX' 3. 'FY' 4. 'NATURE' 'DISCRET') ;
  48. *
  49. * Appel au solveur
  50. *
  51. DEP1 = RESO M_M CHP_F ;
  52. *
  53. TEMPS ;
  54. *
  55. * Code de bon fonctionnement
  56. *
  57. P1_UX = EXTRAIRE DEP1 P1 'UX' ;
  58. P1_UY = EXTRAIRE DEP1 P1 'UY' ;
  59. P2_UX = EXTRAIRE DEP1 P2 'UX' ;
  60. P2_UY = EXTRAIRE DEP1 P2 'UY' ;
  61. R_P1_UX = 1.03006E-02 ;
  62. R_P1_UY = 1.00463E-02 ;
  63. R_P2_UX = 9.96486E-03 ;
  64. R_P2_UY = 9.92506E-03 ;
  65. RES1 = ABS ( ( P1_UX - R_P1_UX ) / R_P1_UX ) ;
  66. RES2 = ABS ( ( P1_UY - R_P1_UY ) / R_P1_UY ) ;
  67. RES3 = ABS ( ( P2_UX - R_P2_UX ) / R_P2_UX ) ;
  68. RES4 = ABS ( ( P2_UY - R_P2_UY ) / R_P2_UY ) ;
  69. SI ( RES1 <EG 1E-2 ) ;
  70. ERRE 0 ;
  71. SINON ;
  72. ERRE 5 ;
  73. FINSI ;
  74. SI ( RES2 <EG 1E-2 ) ;
  75. ERRE 0 ;
  76. SINON ;
  77. ERRE 5 ;
  78. FINSI ;
  79. SI ( RES3 <EG 1E-2 ) ;
  80. ERRE 0 ;
  81. SINON ;
  82. ERRE 5 ;
  83. FINSI ;
  84. SI ( RES4 <EG 1E-2 ) ;
  85. ERRE 0 ;
  86. SINON ;
  87. ERRE 5 ;
  88. FINSI ;
  89. *
  90. FIN ;
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  

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