Télécharger la plaquette Cast3M


fr en

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