* fichier :  fsi3.dgibi
* CAS TEST DU 91/10/04     PROVENANCE : PETI                                    
OPTI ECHO 0;                                                                    
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;
*
GRAPH = 'N';
*
*================================================================*
*               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;
*
*================================================================*
*      modl - materiau - rigidite - masse                        *
*================================================================*
*
MODCOQ=MODL STEEL MECANIQUE COQ2  ;
MODLIQ=MODL WATER LIQUIDE LQU4   ;
MODRAC=MODL RAC1 MECANIQUE LIQUIDE RACO;

*
MATCOQ=MATR MODCOQ RHO 7.8E3 YOUN 2.E11 NU 0.3 ;
MATLIQ = MATR (MODLIQ ET MODRAC)
         RHO 1.E3 RORF 1.E3 CSON 1435. CREF 1435. LCAR 1. G 0.;
CARCOQ=CARB MODCOQ EPAI 0.002;
CARRAC=CARB 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   = TIRE RESUL FREQ RANG 1;
  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  = TIRE RESUL DEPL RANG 1;
     DEF0 = DEFO  STEEL  MOD   0.  ;
     DEF  = DEFO STEEL  MOD   VERT;
     TRAC (DEF   ET DEF0);
  FINSI;
*
*==============================================================*
*                   code validation                            *
*==============================================================*
*
ERGMAX = ERG      ;
*
SI (ERGMAX <EG 5.);
   ERRE 0;
SINON;
   ERRE 5;
FINSI;
*
SAUT 1 LIGN;
TEMPS;
SAUT 1 LIGN;
*
FIN;
*