* fichier : dy_dev12.dgibi ************************************************************************ * * Cas-Test de la liaison palier de l'operateur DYNE * * * Masse ponctuelle sur 1 squeeze-film (NL) * Chargement statique = poids propre = -m*g --> Wcharg * + Chargement tournant = balourd = m*u*w^2 --> Wbal * avec ou sans centering spring de raideur k * * * ref : JI Inayat-Hussain et al, On the bifurcations of a rigid rotor * response in squeeze-film dampers, Journal of Fluids and * Structures 17 (2003) 433-459 * * * Auteur : BP, 2015-05-05 * * * Mots-clés : Vibrations non-lineaires, DYNE, squeeze-film, * dynamique des machines tournantes * ************************************************************************ * OPTI 'DIME' 3 'ELEM' SEG2 ; COMPLET = FAUX; * COMPLET = VRAI; * GRAPH = FAUX; GRAPH = VRAI; * FLCOUR=VRAI pour utiliser le modele de palier_court (sinon rodeli) * FLCOUR = FAUX; FLCOUR = VRAI; * devogelaere ou diff centrees ? FLDEVO = vrai; FLDEVO = faux; * donnees du squeeze : Lpalier = 0.05; Rarbre = 0.02; * mu = 0.02; rhofl = 900.; xjeu = 0.50E-3; xjeum1 = 1./xjeu; * adimensionne B = 0.015; valB = mot'0015'; * B = 0.09; valB = mot'009'; * B = 0.40; valB = mot'040'; * donnees du rotor (rigide) : M1 = 1.; *amplitude force statique g = 10.; * poids = M1*g; * poids propre adimensionne W = 0.05; valW = mot '005'; * W = 0.1; valW = mot '010'; * W = 0.5; valW = mot '050'; * on deduit la vitesse de rotation si (EGA W 0.); poids = 0.; OMEGHZ = 100. ; OMEGA = OMEGHZ * (2*pi); sinon; poids = M1*g; OMEGA = (g / (xjeu * W))**0.5; OMEGHZ = OMEGA/(2*pi); mess '\W = ' OMEGA 'rad/s = ' OMEGHZ 'Hz' ; * on deduit la viscosite mu = B * M1 * OMEGA * (xjeu**3) / (Rarbre * (Lpalier**3)); mess '\m = ' mu ' unite SI ' ; * amplitude balourd * Fbalad = 0.05; valU = mot'005'; Fbalad = 0.10; valU = mot'010'; Fbalad = 0.20; valU = mot'020'; Fbalad = 0.30; valU = mot'030'; Fbalad = 0.355; valU = mot'0355'; * Fbalad = 0.363; valU = mot'0363'; * Fbalad = 0.3635; valU = mot'03635'; * * Fbalad = 0.364; valU = mot'0364'; * * Fbalad = 0.3641; valU = mot'03641'; * * Fbalad = 0.370; valU = mot'0370'; * Fbalad = 0.385; valU = mot'0385'; * * Fbalad = 0.40; valU = mot'040'; * Fbal = Fbalad * M1 * xjeu * (OMEGA**2); Fbal = Fbalad * M1 * xjeu * (OMEGA**2); mess 'Fbal='Fbal 'N'; * centering spring * S = 0.10; valS = mot '010'; S = 0.30; valS = mot '030'; * S = 0.50; valS = mot '050'; k1 = M1 * ((S * OMEGA)**2); mess 'k='k1 'N/m'; Y0 = W / (S**2); mess 'static eccentricity ratio Y0='Y0; * valeurs adimensionnées du squeeze : cha1 = chai 'Squeeze-film : B=' FORMAT '(F5.3)' B ' - W=' FORMAT '(F5.3)' W ' - S=' FORMAT '(F5.3)' S ' - U=' FORMAT '(F5.3)' Fbalad ; mess cha1; TITR cha1; * prefixe si FLCOUR; mopal = mot 'court'; sinon; mopal = mot 'rodeli'; finsi; si FLDEVO; moalgo = mot 'DV'; sinon; moalgo = mot 'DC'; finsi; prefix = chai 'dy_dev12-'mopal'-'moalgo '-B'valB'-W'valW'-S'valS'-U'valU; iCI = 0; prefix = chai prefix '-CI0'; * iCI = 1; prefix = chai prefix '-CI1'; * iCI = 2; prefix = chai prefix '-CI2'; * iCI = 3; prefix = chai prefix '-CI3'; mess prefix ; ficxdr = chai prefix'.xdr'; ficps = chai prefix'.ps'; OPTI 'TRAC' 'PSC' 'POTR' 'HELVETICA_16' 'EPTR' 6 'FTRA' ficps; ************************************************ * Nombre de pas de temps, pas de temps ************************************************ si (COMPLET); si (((Fbalad > 0.36 ) et (Fbalad 1.E-6) ou ((extr zecart 1) < 0.1) ou ((extr zecart 2) < 0.1) ou ((extr zecart 3) < 0.1) ou ((extr zecart 4) > 1.E-6) ); ERRE 5; SINON; ERRE 0; FINSI; fin ;