* fichier : reso_asy.dgibi ************************************************************************ * Section : Mathematiques Fonctions ************************************************************************ * 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