Test name
dyna7
Calculation type
MECHANICS DYNAMICS 2D FOURIER
Finite element type
LQUA4, LSU2
Topic
Vibrating mode shapes of a thin fluid web.
The cavity is supposed undeformable. The gravity effect is taken into consideration.
The shaking mode shapes are computed in an annular space with Fourier modes number 1 and 3.
The three options of the VIBR operator are tested in this example.
Goal
Find the frequencies of the first and third Fourier shape.
Reference
Analytical solution.
Version
97' customer version
Model description
Test dyna7 Results
Results
CASTEM figures
* Test Dyna7.dgibi: Jeux de données *
* --------------------------------- *
* *
OPTI ECHO 0;
GRAPH='N';
SAUT PAGE ;
*
*******************************************************
* DYNA7
* CAS TEST OPERATEUR VIBRATION
*
*
* Cet exemple permet de tester les trois options de
* l'opérateur VIBR
*
* On calcule les modes de ballotement d'une lame mince
* de fluide dans un espace annulaire en mode de
* Fourier 1 et 3.
*
* On suppose la cavité indéformable. L'effet de la
* gravité est pris en compte.
*
*
*******************************************************
TITRE 'MODES DE BALLOTEMENT D UNE LAME MINCE DE FLUIDE';
OPTI DIME 2 ELEM QUA4;
GRAPH='N';
TEMPS ;
*--------------------------- MAILLAGE -----------------
*----------------- ON NE MAILLE QUE LE FLUIDE ---------
*
A1=0.995 0.;A2=1.005 0.;A3=1.005 2.;A4=0.995 2.;
A5=1.005 1.2;A6=0.995 1.2;
L1=D 1 A1 A2;
L2=A2 D 24 A5 D 16 A3;
L3=D 1 A3 A4;
L4=A4 D 16 A6 D 24 A1;
FLUID=DALL L1 L2 L3 L4;
*
SI (NEG GRAPH 'N');
TRAC FLUID ;
FINSI;
*
OPTI ELEM SEG2;
SURFA=D 1 A3 A4;
*
SI (NEG GRAPH 'N');
TRAC SURFA ;
FINSI;
*
MAIL=FLUID ET SURFA;
*
SI (NEG GRAPH 'N');
TRAC MAIL;
FINSI;
*
*---------- DEFINITION DE QUELQUES CONSTANTES ---------
*
ROREF=1.;CEREF=1.;LCARA=1.;
RO0=1000.;CS=1200.;GRA=9.81;
*
*-- MODELE ET OBJET AFFECTE FLUIDE ET SURFACE LIBRE ---
*-------- POUR L'HARMONIQUE DE FOURIER 1 --------------
*
OPTI MODE FOUR 1;
MODLIQ1=MODL FLUID LIQUIDE LQU4;
MODLIQ2=MODL SURFA LIQUIDE LSU2;
MODLIQ=MODLIQ1 ET MODLIQ2;
*
*--------------- MATERIAU LIQUIDE ---------------------
*
MAT=MATR MODLIQ RHO RO0 RORF ROREF CSON CS CREF CEREF
LCAR LCARA G GRA;
*
*--------------- RIGIDITE ET MASSE --------------------
*
RAI=RIGI MODLIQ MAT;
MAS=MASS MODLIQ MAT;
*
*______________________________________________________
* ON CALCULE LE PREMIER MODE AVEC L'OPTION INTERVALLE
* DE VIBR. ON RENORMALISE LE CHAMP DE DEPLACEMENT
* (ET DE PRESSION) DE MANIERE A CE QUE LA PRESSION SOIT
* EGALE A 1 A L'ALTITUDE 0 ET ON EXTRAIT LES PRESSIONS
* AUX POINTS A2, A5 ET A3
*______________________________________________________
*
MOD1=VIBR INTE 0. 1. RAI MAS;
F1=TIRE MOD1 FREQ RANG 1;
PRE=TIRE MOD1 DEPL RANG 1;
PBAS=EXTR PRE P A2;
PP=PRE/PBAS;
X1=EXTR PP P A2;Y1=EXTR PP P A5;Z1=EXTR PP P A3;
*
*______________________________________________________
* ON CALCULE LE PREMIER MODE AVEC L'OPTION
* PROCHE DE VIBR
*______________________________________________________
*
MOD2=VIBR PROC (PROG 0.5) RAI MAS;
F2=TIRE MOD2 FREQ RANG 1;
PRE=TIRE MOD2 DEPL RANG 1;
PBAS=EXTR PRE P A2;
PP=PRE/PBAS;
X2=EXTR PP P A2;Y2=EXTR PP P A5;Z2=EXTR PP P A3;
*
*______________________________________________________
* ON CALCULE LE PREMIER MODE AVEC L'OPTION SIMULTANE
* DE VIBR
*______________________________________________________
*
MOD3=VIBR SIMU 0.5 1 RAI MAS;
F3=TIRE MOD3 FREQ RANG 1;
PRE=TIRE MOD2 DEPL RANG 1;
PBAS=EXTR PRE P A2;
PP=PRE/PBAS;
X3=EXTR PP P A2;Y3=EXTR PP P A5;Z3=EXTR PP P A3;
*
*______________________________________________________
* ON RECOMMENCE LES OPERATIONS CETTE FOIS
* AVEC L'HARMONIQUE 3
*______________________________________________________
*
OPTI MODE FOUR 3;
MODLIQ1=MODL FLUID LIQUIDE LQU4;
MODLIQ2=MODL SURFA LIQUIDE LSU2;
*
MODLIQ=MODLIQ1 Et MODLIQ2;
*
MAT=MATR MODLIQ RHO RO0 RORF ROREF CSON CS CREF CEREF
LCAR LCARA G GRA;
*
RAI=RIGI MODLIQ MAT;
MAS=MASS MODLIQ MAT;
*
MOD1=VIBR INTE 0. 1. RAI MAS;
*------------------------------------------------------
G1=TIRE MOD1 FREQ RANG 1;
PRE=TIRE MOD1 DEPL RANG 1;
PBAS=EXTR PRE P A2;
PP=PRE/PBAS;
T1=EXTR PP P A2;U1=EXTR PP P A5;V1=EXTR PP P A3;
*------------------------------------------------------
MOD2=VIBR PROC (PROG 0.9) RAI MAS;
*
G2=TIRE MOD2 FREQ RANG 1;
PRE=TIRE MOD2 DEPL RANG 1;
PBAS=EXTR PRE P A2;
PP=PRE/PBAS;
T2=EXTR PP P A2;U2=EXTR PP P A5;V2=EXTR PP P A3;
*------------------------------------------------------
MOD3=VIBR SIMU 0.9 1 RAI MAS;
G3=TIRE MOD3 FREQ RANG 1;
PRE=TIRE MOD3 DEPL RANG 1;
PBAS=EXTR PRE P A2;
PP=PRE/PBAS;
T3=EXTR PP P A2;U3=EXTR PP P A5;V3=EXTR PP P A3;
*------------------------------------------------------
SAUT PAGE; SAUT 2 LIGNE ;
MESS '*';
MESS '******** ON VERIFIE LES FREQUENCES';
MESS '*';
MESS ' ';
MESS
' ANALYTIQUE * INTERVALLE * PROCHE
* SIMULTANE';
MESS
' -------------*-----------------*-----------
-----*-------------';
MESS ' HARM 1 0.489 * ' F1 '*' F2 '* ' F3;
MESS ' HARM 3 0.863 * ' G1 '*' G2 '* ' G3;
MESS ' ';
MESS '*';
MESS
'******** ON VERIFIE LE CHAMP DE PRESSION AUX POINTS
A1,A5 ET A3';
MESS '*';
MESS ' ';
MESS
' ANALYTIQUE * INTERVALLE * PROCHE
* SIMULTANE';
MESS
' -------------*-----------------*-----------
-----*-------------';
MESS 'HARM 1 1.00 * ' X1 '*' X2 '* ' X3 ;
MESS 'HARM 1 1.81 * ' Y1 '*' Y2 '* ' Y3 ;
MESS 'HARM 1 3.76 * ' Z1 '*' Z2 '* ' Z3 ;
MESS 'HARM 3 1.00 * ' T1 '*' T2 '* ' T3 ;
MESS 'HARM 3 18.31 * ' U1 '*' U2 '* ' U3 ;
MESS 'HARM 3 201.72 * ' V1 '*' V2 '* ' V3 ;
TEMPS ;
*------------- code de bon fonctionnement ------------*
SAUT 1 PAGE;
ERR=TABLE;
ERR.1=100*(ABS(0.489 -F1)/0.489);
ERR.2=100*(ABS(0.489-F2)/0.489);
ERR.3=100*(ABS(0.489 -F3)/0.489);
ERR.4=100*(ABS(0.863 -G1)/0.863);
ERR.5=100*(ABS(0.863 -G2)/0.863);
ERR.6=100*(ABS(0.863 -G3)/0.863);
ERR.7=100*(ABS(1-X1)/1);
ERR.8=100*(ABS(1-X2)/1);
ERR.9=100*(ABS(1-X3)/1);
ERR.10=100*(ABS(1.81 -Y1)/1.81);
ERR.11=100*(ABS(1.81 -Y2)/1.81);
ERR.12=100*(ABS(1.81 -Y3)/1.81);
ERR.13=100*(ABS(3.76 -Z1)/3.76);
ERR.14=100*(ABS(3.76 -Z2)/3.76);
ERR.15=100*(ABS(3.76 -Z3)/3.76);
ERR.16=100*(ABS(1-T1));
ERR.17=100*(ABS(1-T2));
ERR.18=100*(ABS(1-T3));
ERR.19=100*(ABS(18.31 -U1)/18.31);
ERR.20=100*(ABS(18.31 -U2)/18.31);
ERR.21=100*(ABS(18.31 -U3)/18.31);
ERR.22=100*(ABS(201.72-V1)/201.72);
ERR.23=100*(ABS(201.72-V2)/201.72);
ERR.24=100*(ABS(201.72-V3)/201.72);
I=0;
REPETER BOUC1 24;
I=I+1;
SI (ERR.I <EG 5);
SAUT 1 LIGN;
ERRE 0;
SINON;
ERRE 5;
MESS 'erreur valeur ' I ' du tableau';
FINSI;
FIN BOUC1;
*-----------------------------------------------------*
FIN;
Test dyna7 Comments
MODLIQ1 = MODE FLUID LIQUIDE LQU4;
MODLIQ2 = MODE SURFA LIQUIDE LSU2;
The operator MODEL with keyword LIQUIDE defines a liquid formulation. The type of finite elements depends of the type of formulation. LQU4 is a four node massive quadrangle and LSU2 is a two node element of free surface.
MAT=MATE MODLIQ RHO RO0 RORF ROREF CSON CS CREF CEREF LCAR LCARA G GRA;
The material characteristics are defined with MATER operator and use the following parameters :
PRE = TIRE MOD1 DEPL RANG 1;
PBAS = EXTR PRE P A2;
The object PRE is a ponctual field and one of its components is the pressure P. The data PBAS is the pressure at point A2.