* 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