Numérotation des lignes :

* fichier :  fsi3.dgibi************************************************************************************************************************************************* CAS TEST DU 91/10/04     PROVENANCE : PETI ********************************************    Test fsi3.dgibi: jeux de données     **    ---------------------------------    ******************************************** * SI GRAPH = N PAS DE GRAPHIQUE AFFICHE* SINON SI GRAPH DIFFERENT DE N TOUS* LES GRAPHIQUES SONT AFFICHES GRAPH = 'N' ; SAUT PAGE;SI (NEG GRAPH 'N') ;  OPTI ECHO 1   ;  OPTI TRAC X   ;SINO ;  OPTI ECHO 0   ;FINSI ; SAUT PAGE;********************************************                                         **               TEST FSI3                 **                                         **    CYLINDRICAL TANK CONTAINING WATER    **                                         **  Calculation of the frequencies of      **  cylindrical tank containing water, for **  m = 30.  The radius is 1.43m and       **  height is 1.039m.                      **                                         **                                         **                            P2           **                           +             **                           |             **                           |             **                fluid      |             **                           |             **                           |             **                           |             **             --------------+             **             P4              P3          **                                         ** The boundary condition :    p3 is fixed **                                         **      dp  |                              **     ---- |                 = 0.         **      dz  | h = 0. , 1.43                **                                         **      dp  |                              **     ---- |                 = 0.         **      dr  | r = R                        ******************************************** TEMPS; OPTI DIME 2;OPTI MODE FOUR 10 ;OPTI ELEM QUA4; *=========================================**               geometry                  **=========================================** Dimensions en metres * Points N1 = 20       ; N2 = 30          ;P2 = 1.43 1.039  ;P3 = 1.43 0.0    ;P4= 0. 0.        ;P5 = 1.43 1.039  ;P6 = 1.43 0.0    ;P7= 0. 0.        ;P8 = 0. 1.039    ; S2 = P3 D N1 P2 ; S3 =P3 D N2 P4 ;S4 = P6 D N1 P5 ; S5 =P7 D N2 P6 ;S6 = P7 D N1 P8 ; S7 =P8 D N2 P5 ; STEEL = S2 ET S3                 ;WATER = DALL S4 S5 S6 S7  QUELC   ;RAC1  = RACC 0.002 WATER STEEL   ; * OPTIO FOR TRACE SI (NEG GRAPH 'N');  TITR ' FSI3 : MAILLAGE'; TRAC QUAL (WATER ET (0 0));FINSI; *=========================================**      MODE - materiau - rigidite - masse **=========================================**MODCOQ = MODE STEEL MECANIQUE COQ2  ;MODLIQ = MODE WATER LIQUIDE   LQU4   ;MODRAC = MODE RAC1  MECANIQUE LIQUIDE RACO; MATCOQ = MATE MODCOQ RHO 7.8E3 YOUN 2.E11                      NU 0.3 ;MATLIQ = MATE (MODLIQ ET MODRAC)         RHO 1.E3 RORF 1.E3 CSON 1435.          CREF 1435. LCAR 1. G 0.; CARCOQ = CARA MODCOQ EPAI 0.002;CARRAC = CARA MODRAC LIQU WATER;MATCOQ = MATCOQ ET CARCOQ;MATLIQ = MATLIQ ET CARRAC; RIG1 = RIGI  (MODCOQ ET MODLIQ ET MODRAC)             (MATCOQ ET MATLIQ)  ;MAS1 = MASS  (MODCOQ ET MODLIQ ET MODRAC)             (MATCOQ ET MATLIQ)   ; *=========================================**                boundary conditions      **=========================================* RES1 = BLOQ DEPL ROTA P3 ;**=========================================** theritical value  of frequency FRE1  = 106.64 ;LIST1 = PROG FRE1 ; *=========================================**     calculation of the frequencies      **                   and                   **       extraction of some results        **=========================================* * Use of the operator VIBR. (option PROC) RESUL = VIBR PROC LIST1  (RES1 ET RIG1 )  MAS1 ; *=========================================**             results                     **=========================================* MESS ' RESULTATS ';MESS ' --------- ';SAUT 1 LIGN;   FRE2  = RESUL . MODES . 1 . FREQUENCE;  ERG   = 100 *         (ABS ((FRE1   - FRE2  ) / FRE1  ));  MESS ' MODE '  1  ;  MESS ' ----------';  MESS 'Frequence theorique :' FRE1   'Hz';  MESS 'Frequence calculee  :' FRE2   'Hz';  MESS '    Soit un ecart de : ' ERG   '%';  SAUT 1 LIGN; * mode shape extraction and optional trace   SI (NEG GRAPH 'N');     MOD  = RESUL . MODES . 1 . DEFORMEE_MODALE;     DEF0 = DEFO STEEL  MOD   0.  ;     DEF  = DEFO STEEL  MOD   BLEU;     TRAC (DEF   ET DEF0);  FINSI; *=========================================**           code validation               **=========================================* ERGMAX = ERG      ; SI (ERGMAX &lt;EG 5.);   ERRE 0;SINON;   ERRE 5;FINSI; SAUT 1 LIGN;TEMPS;SAUT 1 LIGN; FIN;

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