JOURNAL OF THEORETICAL AND APPLIED MECHANICS 43, 3, pp. 501-509, Warsaw 2005 ADVANCED MODELLING OF VIBRATORY MACHINES Mariusz Giergiel AGH University of Science and Technology. Cracow e-mail: giergiel@agh.edu.pl Józef Giergiel Rzeszow University of Technology Along with increasing dimensions, speed, and output of modern machi- nes,manynewlyappearingproblemsarebeingencountered indesignand manufacturingof suchmachines.Vibrationsoccuring in certain technolo- gical processes create undesirable phenomenawhich limit the durability, cause excessive dynamic reactions with surroundings or increase level of emitted sound. For a certain class of machines, however, vibration con- stitutes the primary factor ensuring that a desired technological process is performed correctly. This class is hereafter referred to as vibratory machines. Based on the considerationsmade in this paper, a conclusion can be drawn that from the point of view of mechatronical design of vi- bratorymachines it is necessary to determine precisemodels ofmachines including the driving system, especially the induction motor. Key words: vibratorymachines, machine dynamics, modelling 1. Introduction Machines which fulfil a manufacturing or a transport process based on transmission of vibrations of the body of amachine to themachinedmedium, i.e. so-called vibrating machines, are the only, and perhaps the most advan- tageous form of the realization of a manufacturing process in many branches of industry. The mode of working of such a machine, which is based on cau- sing intensive polytonal vibrations, often results in transmission of vibration to the ground, which, in turn, has a disturbing effect to the surroundings and is detrimental to health of people and buildings they live in. Amplitudes of 502 M.Giergiel, J.Giergiel dynamic forces transmitted to the ground achieving up to 20kkN, or evenmo- re, together with frequencies up to 100Hz are the major source of vibration in industrial plants. CAEsystemsbecame themost necessary and essential instrument of today designer’s work. Computer and applicable software offers a lot of helpful func- tions, especially in analysis, modeling and simulation. Today design engineer with CAE system can create a model of a machine and observe it by perfor- ming necessary analysis and simulations as well. What is more, it is possible to make durability calculations, analysis and even dynamical simulations of the working machine. To realize all these purposes, in the design of vibratory machines it is necessary to use precise models of dynamic events in the electromechanical system: vibratory machine - inductive drive. 2. The mathematical model of a vibratory machine 2.1. Model of the electrical part Thepropulsionof vibratorymachines is usually builtwith single run three- phased asynchronous engines. It is indispensably to use suitablemodel of pro- pulsion, when performing simulation. The asynchronous machine is the most popular electromechanical converter used in systems with various dynamical structures. Dynamical advances that exist during the starting phase, braking, voltage-hesitation are the basic gauge during mechanical, electrical and ther- mal selection of engines for designed propulsion. Nonlinear differential systems determine the mathematical model of an engine. The nonlinearity results from relationships between generalized coor- dinates describingmechanical and electric quantities. The structure of engine electromagnetic configuration classifies the type of a machine to holonomi- cal systems, where equations of constraints are integrable. The dependence between electromagnetic induction in a slit and voltage is linear. It is possible to use double or single-frame models of an asynchronous machine. What is more, we are able to explore the Kloss model or a sta- tic characteristic as well. However, from exchanged models, only single and double-framemodels suit dynamics of the asynchronous machine. The use of a double-frame model permits one to describe a wide group of asynchronous machines. This model gives extraordinary precision of calcula- tions including static characteristics. The single-framed model has quite big Advanced modelling of vibratory machines 503 precision for annular and single-framedmachines with cages ring-shaped bars about round rods (Puchała, 1977). Approaching the problem of mathematical modelling of an asynchronous machine, it is necessary to take care of basic purposes, which should be execu- ted.Thefirst thing is thatwe expectmathematicalmodels in formof relations between voltages, currents, torques and rotation speed. We are especially in- terested in dynamical conditions like the starting phase, braking, charge and discharge phases, changes during the switching of voltages and resistance or finally short-circuit.Themost important conditions required to obtain a useful mathematical model are the verifiability of its adequacy and the possibility to fix parameters throughmeasurements. The double-frame model of an asynchronous machine described in terms of stator streams and two frames divided by the carriage voltage, denoted by Φs, Φ1, Φ2, can be expressed by following equations (Rams, 1977) Φ̇s =−as ( Φs− σ2 σ Φ1− σ1 σ Φ2 ) +jωxΦs+U Φ̇1 =−a1 ( − σ2 σ Φs+ σ∆s+σs2 σ Φ1− σ∆s σ Φ2 ) − j(ωx−ωe)Φ1 Φ̇2 =−a2 ( − σ1 σ Φs+ σ∆s+σs1 σ Φ2− σ∆s σ Φ1 ) − j(ωx−ωe)Φ2 (2.1) Mel =−C · Im [ Φs (σ2 σ Φ1+ σ1 σ Φ2 )] C = 2Muω 2 op Iw where Mel – electromagnetic torque Mu – breakdown torque σw,σ1,σ2,σs1,σs2,σ∆s – parameters of the motor model C – electromechanical constant of the motor p – number of pairs of poles of the asynchronous engine U – normalized excitation voltage ωo – synchronous speed ωe – electrical angular velocity ωx – angular velocity of the stator coordinate system Iw – mass moment of inertia of the rotor Because of difficulty in the identification of parameters of the system and restricted precision, especially when interested in electromechanical worth, we should use the simplest model as it is possible. The parameters describing the double-framemodelmust be consistent with engine characteristics and its 504 M.Giergiel, J.Giergiel structural parameters. Usually, even if all data of the engine are known, it is necessary to use optimization methods to obtain right parameters. Moreover, a double-frame machine has parameters which are not defined in literature, as it goes with single-frame models. Summing up, there is no possibility to get accurate parameters even if catalogue data of the engine are well known. The simplicity and possibility of using catalogue data of engines are themain reasons why the single-framemodel has become themost popular. The stator of amachine has p pole pairs of clutches and the rotor is single-frame. Such a model is described by equations (Puchała, 1977) Φ̇s =−as ( Φs− 1 kw ∓w ) +U Φ̇w =−aw ( Φw ks kw − 1 kw Φs ) − jωeΦw (2.2) Mel =−C · Im(ΦsΦw) C = 2Muω 2 op Iw The model of an asynchronous machine is expressed by equations, where the linked streams of the rotor and stator divded by the worth of carriage voltage, marked asΦw and Φs, are condition co-ordinates. The following denote kw,ks – dissipation of the rotor and the stator, respectively aw,as – proportional parameters of the engine model 2.2. The model of the mechanical part For investigation purposes, the mathematical models describing relations between themachine body and its propulsionwere used for digital simulation of dynamical effects (ultra-resonancemachine). Thesemodelsmake it possible to examine amajor part of difficult dynamical problemsunable to explore ana- lytically. These dynamical effects have considerable influence on the rotation speed of the system. Somepopularvibratorymachines are equippedwith thependularvibrator, which makes them incline by more degrees of latitude than it is observed in machineswith strong axis vibrator. Ifwe consider amachinewith thependular vibrator, we should remember that the model of the mechanical part include flat motion of the body, pendulum and pulp of the vibrator. Moreover, we should take care of relations between the string suspension and the body of the machine. Our calculation takes also into account the susceptibility of the Advanced modelling of vibratory machines 505 Fig. 1. A physical model of a vibratorymachine with the pendulum inertial vibrator clutch that keeps the propulsion-vibrator system tighter (M+mh+mn)x+(mh+mn)Hα+(mhr+mnR)β−mnesinγ= =−n ∑ [kxx(x−γiα)+kxy(y+µiα)+kxαα]− bxx+mneγ 2 cosγ (M+mh+mn)y+mneγ cosγ = =−n ∑ [kyx(x−γiα)+kyy(y−µiα)+kyαα]− byy+mneγ 2 sinγ (mh+mn)Hx+[Ic+(mh+mn)H 2]α+(mhr+mnR)Hβ−mneHγ sinγ = =−n ∑ {kαx(x−γiα)+kαy(y−µiα)+kααα+ +µi[kyx(y−γiα)+kyy(y−µiα)+kyαα]+ −γi[kxx(x−γiα)+kxy(y−µiα)+kxαα]}+ −bα+mneHγ 2 cosγ+ b1h 2(β−α)+k1h 2(β−α) (2.3) 506 M.Giergiel, J.Giergiel (mhr+mnR)x+(mhr+mnR)Hα+(mhr 2+mnR 2+ Ih)β−mneRγ sinγ= =mneRγ 2 cosγ− b1h 2(β−α)−k1h 2(β−α)+ −mnesinγ+mnecosγ−mneHαsinγ−mneRβ sinγ+(In+mne 2)γ= =−ksm ( γ−φ d1 d2 ) − bsm ( γ−φ d1 d2 ) − boγ 2 sgn(γ−β) [ Iw+ I1+ I2 d21 d22 ] ϕ=−ksm ( φ d1 d2 −γ )d1 d2 − bsm ( φ d1 d2 −γ )d1 d2 + +Mel−Mo sgnϕ ReΦ̇s =−as ( ReΦs− 1 kw ReΦw ) +ReU ImΦ̇s =−as ( ImΦs− 1 kw ImΦw ) +ImU ReΦ̇w =−aw ( ReΦw ks kw − 1 kw ReΦs ) − ϕ̇ipImΦw (2.4) ImΦ̇w =−aw ( ImΦw ks kw − 1 kw ImΦs ) − ϕ̇ipReΦw Mel =−C(ReΦsImΦw− ImΦsReΦw)Iw Equations (2.3)1,2,3 describe flat motion of the machine body. Equation (2.3)4 describes motion of the pendulum,while equation (2.3)5 expresses mo- tion of the vibrator. Equation (2.3)6 describes rotarymotion of the engine and equations (2.4) define dynamical effects in the propulsion system. The following denote: x,y,α – horizontal coordinates of themass center of mass and angle of deviation from the machine body level β – angle of inclination of the pendulummeasured with respect to machine level ϕ – angle of rotation of the engine γ – angle of rotation of the vibrator Φs,Φw – fluxes of the stator and rotor divided by reference voltage M – mass of the machine body mh – mass of the pendulum mn – unbalanced mass of the vibrator I – mass moment of inertia of the driving system with respect to the rotor axis Advanced modelling of vibratory machines 507 Ic,Ih,In – central moment of inertia of the machine body, the pendu- lum and the unbalanced mass, respectively Iw – moment of inertia of the rotor I1,I2 – moments of belt pulleys H – distance between the axle of the pendulum and the center of mass of the machine body h – distance between the point of fastening of the pendulum stabilization system and its rotation axis e – eccentricity of the vibrator R – distance between the pendulum axis and vibrator axis r – distance between the center of mass of the pendulum and its rotation axis d1,d2 – diameters of belt pulleys µi,νi – coordinates of the fixing points of elastic elements in the central coordinate system Cµν k1,b1 – coefficients of elasticity and viscous damping of the pendu- lum stabilization system ksm,bsm – coefficients of elasticity and damping of the coupling αs,αw – proportional parameters of the engine model n – number of elastic elements in the set Mo – anti-torque of the driving system i – transmission ratio of the belt drive between the motor and vibrator bo – resistance coefficient of the bearings. During simulation carried out on the presented model of the ultra- resonance vibrationmachine, we take care of precision in the determination of parameters of the string suspension system. Temporary resonance simulations demand high standards of the string suspension system and are dependent on flexibility and damping features of the system aswell as relations between thr co-ordinates. The model of the ultra-resonance vibration machine, which is still in the building phase, enables one to describe the string-damping reaction of the suspension. If we are able to qualify the elasticity matrix K, we can also generally describe the string element in a matrix form as follows    kxx kxy kxα kyx kyy kyα kαx kαy kαα       xi yi αi    =    pxi pyi pαi    (2.5) where 508 M.Giergiel, J.Giergiel xi,yi,αi – displacements of the upper end of the elastic element pxi,pyi,pαi – components of forces andmoments applied to it. Because of constant character of relations between the co-ordinates, we may use an equivalent viscosity damping for motion along co-ordinates mar- ked as: bx, by, bα. In most of popular types of the string suspension (steel strings with mounted endings or a multi-layer metal –rubber block, the ma- trix of damping and elasticity coefficients should be automatically generated in a special model, which concerns Clayperon’s (Michalczyk, 1995) influence of string forces on the flexibility in all angular and linear directions of the deformation. 3. Conclusions To realize all the above mentioned purposes in the design of vibratory machines, it is necessary to use a precise model of dynamic effects in the electromechanical system: vibratory machine – inductive drive. The model was verified by M.Giergiel (2002) and to the best author’s knowledge it can be used for analysis and simulation of dynamics of vibratory machines in the design process, virtual prototyping andmechatronical design. References 1. Banaszewski T., 1991,Przesiewacze, Ossolineum,Kraków 2. Giergiel M., 2002, Komputerowe wspomaganie w projektowaniu maszyn wi- bracyjnych, Wydawnictwo IGSMiE PAN, Studia, rozprawy, monografie, Kra- ków 3. Giergiel M., Sapińska-Wcisło A., 2004, Zastosowanie elementów inteli- gentnych w zawieszeniach nadrezonansowych maszyn wibracyjnych, Przegląd Mechaniczny, 6 4. Michalczyk J., 1995, Maszyny wibracyjne: obliczenia dynamiczne, drgania i hałas, WNT,Warszawa 5. PuchałaA., 1977,Dynamika maszyn i układów elektromechanciznych, PWN, Warszawa 6. Rams W., 1977, Synteza dynamicznych własności maszyn asynchronicznych, Praca doktorska, AGH Advanced modelling of vibratory machines 509 Zaawansowane modelowanie maszyn wibracyjnych Streszczenie Wraz ze zwiększeniem wymiarów, prędkości oraz wydajności współczesnych ma- szyn i urządzeń występują coraz większe problemy w projektowaniu i konstrukcji, wynikające z występującymi w nich drganiami. W procesie projektowania maszyn wibracyjnychniezwykle interesujący jest aspektmożliwościwydatnego skrócenia fazy projektowejmaszynpoprzez szczegółoweokreśleniew fazie koncepcji projektowejwła- sności ruchowych maszyny. W pracy przedstawiono zagadnienia związane z budową modelu fizycznegomaszynywibracyjnej zgodnego z obiektem rzeczywistym i użytecz- nego w procesie komputerowego wspomagania prac inżynierskich. Szczególną uwagę zwrócononamodelowanie silnikanapędzającegowibrator, co z jednej stronyzapewnia prawidłowe odzwierciedlenie istotnych zjawisk fizycznych, z drugiej możliwość iden- tyfikacji parametrówmodelu. Manuscript received April 18, 2005; accepted for print June 2, 2005