* fichier :  fsi4.dgibi
*           Test Fsi4.dgibi: Jeux de données          *
*           --------------------------------          *
*                                                     *
* CAS TEST DU 91/10/04     PROVENANCE : PETI
OPTI ECHO 0;
SAUT PAGE;
*
*******************************************************
*                                                     *
*                        TEST FSI4                    *
*                                                     *
*           CONCENTRIC SHELLS COUPLED BY FLUID        *
*                                                     *
*  Calculation of the frequencies of two concentric   *
*  shells with springs coulped by water in between.   *
*  The inner radius is 0.5m and the outer radius is   *
*  1.0m                                               *
*                                                     *
*                                                     *
*                P5          P7                       *
*                 |         |                         *
*                 |         |                         *
*                 |         |                         *
*                 |  fluid  |                         *
*                 |         |                         *
*                 |         |                         *
*                 |         |                         *
*                 |         |                         *
*          ^^^^^^^ ^^^^^^^^^ ^^^^^^^                  *
*            K   P6    K   P8   K                     *
*                                                     *
*  The boundary conditions are                        *
*                                                     *
*      dp  |                                          *
*     ---- |                 = 0.                     *
*      dz  | h = 0. , 1.0                             *
*                                                     *
*      dp  |                                          *
*     ---- |                 = 0.                     *
*      dr  | r = Ri , Ro                              *
*******************************************************
*
TEMPS;
*
OPTI DIME 2;
OPTI MODE FOUR 1  ;
OPTI ELEM QUA4;
*
GRAPH = 'N';
*
*=====================================================*
*               geometry                              *
*=====================================================*
*
* Dimensions en metres
*
* Points
*
N1 = 1        ; N2 = 1           ;
P5 = 500.0 1000.0 ; P6 = 500.0 0.0 ;
P7 = 1000. 1000.0 ; P8 = 1000. 0.0 ;
*
*
*
S5 =P5 D N1 P6 ; S6 = P8 D N1 P7 ;
*
STEEL = S5 ET S6                 ;
WATER = (S5 PLUS (2.5 0.))  TRAN N2 (495.0 0. )    ;
RAC1  = RACC 3.0   WATER STEEL   ;
*
* OPTIO FOR TRACE
*
SI (NEG GRAPH 'N');
  TITR ' FSI4 : MAILLAGE';
  TRAC QUAL (WATER ET (0 0));
FINSI;
*
*=====================================================*
*    modele - affecte - materiau - rigidite - masse   *
*=====================================================*
*
MODCOQ=MODL STEEL MECANIQUE COQ2 ;
MODLIQ=MODL WATER LIQUIDE  LQU4;
MODRAC=MODL RAC1 MECANIQUE LIQUIDE RACO;

*
MATCOQ=MATR MODCOQ RHO 8.E-9  YOUN 2.E5  NU 0.3 ;
MATLIQ = MATR (MODLIQ ET MODRAC)
         RHO 1.E-9  RORF 1.E-9  CSON 1435.0E3 CREF 
         1435.0E3 LCAR 1000. G 0.;
CARCOQ=CARB MODCOQ EPAI 5.0  ;
CARRAC=CARB MODRAC LIQU WATER;
MATCOQ=MATCOQ ET CARCOQ;
MATLIQ=MATLIQ ET CARRAC;
*
K    =3.7458E4/2.0;
PP   =PROG K (-1*K) K  ;
PPNEW = PROG K ( -1*K) ( -1*K) K;
RIG1 = RIGI  (MODCOQ ET MODLIQ ET MODRAC)
             (MATCOQ ET MATLIQ)   ;
MAS4 = MASS  (MODRAC ET MODLIQ  )(MATLIQ)  ;
*
*
*=====================================================*
*                boundary conditions                  *
*=====================================================*
*
RES1 = BLOQ UZ RT  (S5 ET S6)        ;
RIG2 = (APPU UR UT K P8 ) ET (APPU UR UT K P6 ) ;
RIG3 = ( MANU 'RIGI' type 'RIGIDITE' (P6 D 1 P8) 
       (MOTS UR)  PPNEW )
    ET ( MANU 'RIGI' type 'RIGIDITE' (P6 D 1 P8) 
       (MOTS UT)  PPNEW );
RIG4 = RIG1 ET RIG2         ET RIG3 ;               ;
*
*=====================================================*
* theritical value  of frequency
*
FRE1  = TABLE ;
FRE1.1 = 18.88 ;
FRE1.2 = 55.278;
LIST1 = PROG FRE1.1 FRE1.2 ;             ;
*
*=====================================================*
*              calculation of the frequencies         *
*                            and                      *
*                extraction of some results           *
*=====================================================*
*
*
* Use of the operator VIBR. (option PROC)
*
*
RESUL = VIBR PROC LIST1        (RES1 ET RIG4 )  MAS4  ;
*
*=====================================================*
*                   results                           *
*=====================================================*
*
MESS ' RESULTATS ';
MESS ' --------- ';
SAUT 1 LIGN;
*
*
FRE2 = TABLE; MOD = TABLE; DEF = TABLE;ERG = TABLE;
*
I = 0;
REPETER BLOC1 2;
  I = I + 1;
  FRE2.I = TIRE RESUL FREQ RANG I;
  ERG.I = 100 * (ABS ((FRE1.I - FRE2.I) / FRE1.I));
  MESS ' MODE '  I     ;
  MESS ' ----------';
  MESS ' Frequence theorique :' FRE1.I 'Hz';
  MESS ' Frequence calculee  :' FRE2.I 'Hz';
  MESS '    Soit un ecart de : ' ERG.I '%';
  SAUT 1 LIGN;
FIN BLOC1;
*
*=====================================================*
*                   code validation                   *
*=====================================================*
*
ERGMAX = MAXI (PROG ERG.1 ERG.2 );
*
SI (ERGMAX <EG 5.);
   ERRE 0;
SINON;
   ERRE 5;
FINSI;
*
SAUT 1 LIGN;
TEMPS;
SAUT 1 LIGN;
*
FIN;
*