Download reso_asy.dgibi

Back to the list

* fichier :  reso_asy.dgibi
* CAS TEST DU 91/06/18     PROVENANCE : PLAF
* si GRAPH = N, les graphiques ne sont pas affich{s
* si GRAPH diff{rent de N, tous les graphiques sont affich{s
*
GRAPH = 'N' ;
*
SAUT PAGE ;
*--------------------------------------------------------------------*
*                            RESO_ASY                                *
*--------------------------------------------------------------------*
*
OPTION ECHO 0 ;
TEMPS ;
*
OPTI DIME 2 ELEM SEG2 ;
*
P1 = 0 0 ;
P2 = 1 0 ;
M1 = D 1 P1 P2 ;
*
L_MOT = MOTS 'UX' 'UY' ;
*
*   Matrice antisymétrique
*
M_A = MANU RIGI type 'MASSE' M1 L_MOT 'ANTI'  ( PROG  0    -2    -1   0
                                                2    0    -2   -1
                                                1    2    0   -2
                                                0    1    2    0  ) ;
*
*   Matrice symétrique
*
M_S = MANU RIGI type 'MASSE' M1 L_MOT ( PROG 100   0    0    0
                                         0  200    0     0
                                         0    0  300    0
                                         0    0    0  400  ) ;
*
*   Matrice totale
*
M_M = M_S ET M_A ;
*
*   Vecteur force donn{
*
CHP_F = ( MANU CHPO P1 2 'FX' 1. 'FY' 2. 'NATURE' 'DISCRET') ET
        ( MANU CHPO P2 2 'FX' 3. 'FY' 4. 'NATURE' 'DISCRET') ;
*
*   Appel au solveur
*
DEP1 = RESO  M_M  CHP_F ;
*
TEMPS ;
*
*     Code de bon fonctionnement
*
P1_UX = EXTRAIRE DEP1 P1 'UX' ;
P1_UY = EXTRAIRE DEP1 P1 'UY' ;
P2_UX = EXTRAIRE DEP1 P2 'UX' ;
P2_UY = EXTRAIRE DEP1 P2 'UY' ;
R_P1_UX = 1.03006E-02 ;
R_P1_UY = 1.00463E-02 ;
R_P2_UX = 9.96486E-03 ;
R_P2_UY = 9.92506E-03 ;
RES1 = ABS ( ( P1_UX - R_P1_UX ) / R_P1_UX ) ;
RES2 = ABS ( ( P1_UY - R_P1_UY ) / R_P1_UY ) ;
RES3 = ABS ( ( P2_UX - R_P2_UX ) / R_P2_UX ) ;
RES4 = ABS ( ( P2_UY - R_P2_UY ) / R_P2_UY ) ;
SI ( RES1 <EG 1E-2 ) ;
   ERRE 0 ;
SINON ;
   ERRE 5 ;
FINSI ;
SI ( RES2 <EG 1E-2 ) ;
   ERRE 0 ;
SINON ;
   ERRE 5 ;
FINSI ;
SI ( RES3 <EG 1E-2 ) ;
   ERRE 0 ;
SINON ;
   ERRE 5 ;
FINSI ;
SI ( RES4 <EG 1E-2 ) ;
   ERRE 0 ;
SINON ;
   ERRE 5 ;
FINSI ;
*
FIN ;



 

© Cast3M 2003 - All rights reserved.
Disclaimer