THE RHEOLOGICAL-DYNAMICAL HARMONIC OSCILLATOR UDC (045)=111. Dragan D. Milašinović

Transcription

1 FACTA UNIVERSITATIS Series: Mechanics, Automatic Control and Robotics Vol. 6, N o, 7, pp THE RHEOLOGICAL-YNAMICAL HARMONIC OSCILLATOR UC 53.5(45)= ragan. Milašinović Faculty of Civil Engineering, Subotica, Serbia Abstract. In this paper the rheological-dynamical theory of the analytical dynamics of discrete visco-elasto-plastic system is presented. The rheological-dynamical analogy (RA) has been developed on the basis of mathematical-physical analogy between the rheological model and the discrete dynamical model with viscous damping and is aimed to be used for the analysis of inelastic deforming of materials and structures. In this presentation, the coupled initial conditions of the stress-strain state of the rheological visco-elasto-plastic model are applied for the study of discrete dynamics system. Mechanical systems that have their masses and elastic forces distributed, such as cables, rods, beams, plates, etc., rather than lumped together in concentrated masses by springs belong to this class of vibrations of continuous media. These systems consist of an infinitely larger number of particles, and hence require an infinitely large number of coordinates to specify their configurations. These notes give an example illustrating how discrete model can be derived from special limits of the continuum model using the principle of analogy. This technique is useful because discrete model is often much easier to deal with than continuum model, both conceptually and computationally. ey words: RA analogy, RA harmonic oscillator, viscous damping ratio. INTROUCTION Any system, discrete or continuous, capable of vibrating must have mass and stiffness. Mass implies that once the system starts moving it will possess inertia which will make the system continue moving. On the other hand, stiffness implies that change in the configuration of the system due to some internal or external disturbance will be accompanied by a change in potential energy. However, such a system cannot show vibration characteristics until it is disturbed from the position of equilibrium. A vibrating system would be conservative if its potential energy is interchangeable with the kinetic energy so that over any cycle of operation (at the end of which the configuration of the system is the same as it was at the beginning) the work done is zero. A non-conservative system implies gradual loss of total energy content with time through forces, which do not depend on a potential. Traditionally the energy dissipation in a system, caused by losses in the surrounding medium, feedback in the foundation soil, internal friction of the construction materials, hysteretic behavior of the members of the Received februar 7

2 4.. MILAŠINOVIĆ system and the friction at the connections, is represented by an idealized viscous damping force, i.e. a force directly proportional to the velocity of the system. The external forces are generally induced by several factors such as the earthquake, wind, rotating machines and etc. When a linear analysis is used to predict the dynamic response of a structure, viscous damping is the only factor that takes account of the energy dissipation of the structure and it is well understood that the viscous damping coefficients are of major importance in the dynamic analysis. The structure mass and stiffness matrices remain constant during the analysis and they satisfy the well known orthogonality conditions, Clough and Penzien []. If the damping matrix also satisfies orthogonality, the incremental equations of motion for a discretised multi-degree-of-freedom structure can be decoupled into n independent equations, one for each normal mode of the structure. It is equivalent to assuming that the normal modes of the damped system are identical to those of the undamped system; this is basically true for low values of damping. It is obvious that if the damping matrix is a linear combination of orthogonal matrices, such as mass and stiffness matrices, it will satisfy orthogonality. The strain induced in a purely elastic linear material is proportional to the stress that produces the deformation. When a linear visco-elastic material (hereditary medium) is subjected to time-dependent variations of stress and strain, the fundamental shear or volume deformation of the material is no longer related to stress by a simple constant of proportionality-the shear modulus or bulk modulus. Rather, the relation between stress and strain is most generally represented by a linear partial diffrential equation of arbitrary order. Mechanics of hereditary medium is presented in scientific literature by the array of monograpf: Reiner [] and Goroshko and Hedrih (Stevanović) []. Inelastic deforming of materials is based on the mathematical physical analogy between rheological model and dynamical model with viscous damping, the so-called RA that was suggested in the monograph from Milašinović [4]. RA modeling technique for material behavior of axially cyclically loaded rods has already been explained by Milašinović [5], where RA was used to predict the buckling behavior of slender columns. In the second paper, Milašinović [6], the author demonstrates that RA is also capable of modeling the fatigue behavior. The efficient numerical implementation of RA and its practical application was also studied by Milašinović [7] in the viscoelasto-plastic behavior of metalic rods where the RA loading function for Hencky s theory is derived. The objective of the present paper is to explain the physical mechanism of viscous damping, which is based on analogy. In such an analysis, the structure is usually represented by a simple vibration model. Taking into account RA, rheologic behavior of bar can only be characterized by a single parameter, i.e. the dynamic time of retardation T = / ω, or a characteristic time for which a wave at the velocity c takes to propagate a length of the bar l. In the stage of visco-elasticity, the RA model has the same phase angle as a simple single-degree-of-freedom spring mass system with damping in the steady state vibration and from that the viscous damping ratio is obtained.. THE ISCRETE VISCO-ELASTO-PLASTIC RA SYSTEM Many vibrating objects vibrate normally in a way that establishes standing waves in the object. A standing wave is formed, for instance, in a stretched flexible string when

3 The Rheological-ynamical Harmonic Oscillator 5 one end is fixed and vibratory motion is imparted to other end. The points of a standing wave at which the amplitude is equal to zero are called its nodes; the points at which the amplitude reaches its maximum value are called antinodes. The distance between two adjacent nodes and between two adjacent antinodes are the same and equal to half the wave-length λ of the traveling waves. This distance is called the length of the standing wave: λ st = λ /. We shall now consider the bar of Fig. a) of length l loaded by the axial force P(t), which retains its magnitude and direction as the bar deflects. Suppose that length of the bar under compressive critical load at the point of elasticity E is equal to the length of the standing wave (l = λ st ). In this way we obtained a physical significance for the characteristic length l. Let A be the constant cross section area of the uniform bar where the lateral dimensions remains small in comparison with the wavelength. In this case the radial motion of the bar (Poisson s ratio effect) may be neglected and the bar may be regarded as long bar. The vibration of the bar of significant lateral dimensions is determined from a corrected theory given by Love [3]. Fig.. iscrete visco-elasto-plastic RA system Let y(x,t) be the longitudinal motion of the bar yxt (, ) =ε ( tx ) () The total strain ε(t) is modeled according to the principle of RA by Milašinović [5] d ε dε d σ dσ + + ε = + + σ σ Y m c k m c k k () dt dt dt dt

4 6.. MILAŠINOVIĆ where: λ λ N m= = kt γ Eλ N + H λ, c = = kt γ kt H m =, E c = k( +ϕ ) T E H k = k( +ϕ ), E H k = k( ϕ ϕ ) (3) E H EH A k = l ϕ =ϕ ( + i), i =,,,... T = l EH ρ i is the level for visco-plastic yielding at any stage of dynamic equilibrium. φ * represents the visco-elastic creep coefficient of the material of the bar, see Milašinović [5]. T represents a characteristic time (the dynamic time of retardation) for which a wave at the velocity c = E H / ρ takes to propagate the distance l. The total visco-elasto-plastic strain was solved theoretically on the harmonic excitation σ(t) = σ +σ A sin(ω σ t) by Milašinović [9]. The strain ε h, under constant stress σ < σ Y, with initial conditions (4) is only visco-elastic strain σ ε =ε el =, E H dε ϕ iϕ =σ + ( σ σ ) (4) dt E T E T Y H H The cyclic visco-elasto-plastic steady strain is given by σ ε = ( +ϕ ) (5) h E H σ ( +ϕ ) +δ δϕ A ε p =ε () t = sin ω t arctan σ EH +δ +δ +ϕ (6) The constant and variable components of cycle are σ and σ A respectively, while σ Y is the yield stress. δ represents the relative frequency of the RA model, δ = ω σ /ω, where ω σ = ω F. The RA system (Fig. b) consists from a main mass M, which is periodically excited and its resonance amplitudes should be minimized by the influence of the material of the bar-like visco-elasto-plastic substance under cyclic strain. When gusts occur, the

5 The Rheological-ynamical Harmonic Oscillator 7 cyclic stress causes cyclic visco-elasto-plastic strain to the supstance, which then absorbs energy. The motion y(l,t) of the damper P is described by the equation similar to Eq. where Ad y Ady A d P dp m + c + k y = m + c + kp k PY (7) l dt l dt l dt dt y = y( l, t) =ε ( t) l P = P() t =σ () t A (8) According to the Alambert s principle, the following equation could be written for the dynamic equilibrium of forces acting on the RA system Thus d y M + Pt () = Ft () (9) dt d y d y Pt () = Ft () M, P dt Y = FY M dt Y 3 dp df d y = M () 3 dt dt dt 4 d P d F d y = M dt dt dt 4 The motion of the RA system is described by the equation with the highest derivative of the fourth order with respect to time t or 4 3 A d y A dy A d F d y df d y m + c + k y = m M c M l dt l dt l dt dt dt dt d y d y Y + k F M k F Y M dt dt d y d y d y d y A d y A dy A mm + cm + km k M + m + c + k y = dt dt dt dt dt 4 3 Y 4 3 l l dt l d F df = m + c + kf k F dt dt Y () ()

6 8.. MILAŠINOVIĆ where: m T k = + ϕ, m EHT k = + ϕ c +ϕ = k +ϕ T c EHT, k = + ϕ (3) k EH k = + ϕ, k ϕ ϕ = k + ϕ A mechanical longitudinal wave propagates in an elastic solid metals at the finite velocity, ~5m/s. Accordingly, the characteristic time T to be small number, so that the terms, which have multiplication by (T ) become negligible. For these oscillatory RA systems the third order differential equation with respect to time t may be composed 3 +ϕ d y d y ϕ ϕ d yy T dy 3 +ϕ dt dt +ϕ dt +ϕ dt +ϕ T M + M M + k + ky = +ϕ df ϕ ϕ = T + F F +ϕ dt +ϕ Y (4) In the stage of visco-elasticity φ = φ *, the Eq. 4 takes the form of 3 +ϕ d y d y T dy +ϕ df T M + M + k + ky = T 3 + F +ϕ dt dt +ϕ dt +ϕ +ϕ dt (5) 3. THE VISCO-ELASTIC RA OSCILLATOR amping is basically a dissipation of energy, which occurs in the damper P. Metallic dampers mainly consisting of low carbon steel. Two equivalent expressions for the critical (minimal) damping of the bar defined as mass and/or stiffness proportional can be obtained using by the basic RA equations see Milašinović [8] = = E H + = cmin ccr m m λ λn T (6) λ λ N cmin = ccr = k + = kt E H (7)

7 The Rheological-ynamical Harmonic Oscillator 9 Thus the minimal viscous damping ratio of the RA model is ξ min = (8) Owing to the visco-elastic nature of the material of the damper, the measure of phase angle between the applied stress and subsequent strain is more importantly δ ϕ tan α= (9) + δ +ϕ The representative computations of the previous equation are shown in Fig.. The curves on this figure have been computed for the values of creep coefficient: φ * =.5,,, 3, 4, 5, and 6., tanα,8,6,4, creep coef.=.5 creep coef.= creep coef.= creep coef.=3 creep coef.=4 creep coef.=5 creep coef.=6, logδ Fig.. Frequency dependence of the phase angle of the RA model The RA model has the same natural frequency as a simple single-degree-of-freedom spring mass system and because of that we can form a new single RA system with one lumped element of mass (Fig. 3), where natural and relative frequencies are: keq k eq ω = = =ω m+ M M m ( + ϕ )( +η) + m ωf ωf δ = = ( +ϕ )( +η ) =δ ( +ϕ )( +η) ω ω () ()

8 3.. MILAŠINOVIĆ Fig. 3. A new single RA visco-elastic system with one lumped element of mass k eq is the permanent stiffness k eq k = () + ϕ η is the mass ratio M η = (3) m m is mass of the damper. δ is the relative frequency of the RA model where ω σ δ = = ω ω σ T F F (4) ω =ω (5) Also, the RA model has the same phase angle as a simple single-degree-of-freedom spring mass system with damping in the linear steady state vibration and from that we can form equality, from which yields the viscous damping ratio of single RA system with one lumped element of mass δϕ ξδ ϕ δ( δ ) =, ξ= +δ +ϕ δ ( +δ +ϕ ) δ The viscous damping ratio, which is formulated using the principle of analogy, is a function of various relative parameters: φ *, η and δ. It is known that for the base isolation, the fundamental frequency of the whole structure is dominated by the natural frequency of the base isolator. Consequently, it is very important to evaluate the viscous damping ratio for the isolator, which is the function of mentioned parameters. (6)

9 The Rheological-ynamical Harmonic Oscillator ξ logδ M/m= M/m= M/m= M/m=4 Fig. 4. Frequency (δ) dependence of the viscous damping ratio of single RA system with one lumped element of mass for creep coefficient φ * = The viscous damping ratio is shown in Fig. 4 for the creep coefficient φ * =, where all values are positive. Owing to the frequency dependence of the viscous damping ratio ξ it is useful to consider separately the situations arising when the ξ is positive (system is stable) and when it is negative. Negative damping ratio means that the complementary solution of the response would not die away (the system is unstable because of the factor e ξ ω t ). However, it may be important to analyze the dependence of the viscous damping ratio on the relative frequency δ * (Fig. 5), because this frequency is the function of mass ratio η ξ logδ* M/m= M/m= M/m= M/m=4 Fig. 5. Frequency (δ * ) dependence of the viscous damping ratio of single RA system with one lumped element of mass for creep coefficient φ * =

10 3.. MILAŠINOVIĆ The purpose of the energy dissipation devices is to increase the natural period of the structure so that the acceleration response of the structure is decreased during the gust. The effective period T eqd for the RA system can be written as or T eqd * π π π ( +ϕ )( +η) * * d = = = ω ω ξ ω ξ T eq π π = = +ϕ +µ * ω ω where T eq is the natural period of the structure. The effective period ratio can be represented as * ( )( ) (7) (8) T eqd T = * ( + ϕ )( +η) ξ (9) or T eq T * = ( +ϕ )( +η ) (3) where T is the the period of the energy dissipation devices π T = = πt (3) ω The effective period ratios of the RA system, shown in Fig. 6 are constants for all mass ratios while viscous damping ratios have positive values. However, the effective period ratio increases rapidly as the applied mass ratio increases. Teqd/T. logδ* M/m= M/m= M/m= M/m=4 Fig. 6. Frequency (δ * ) dependence of the effective period ratio of single RA system with one lumped element of mass for creep coefficient φ * =

11 The Rheological-ynamical Harmonic Oscillator 33 This time considering also a new single RA system with one lumped element of mass with effective stiffness k eq =k / (+φ * ) and effective damping c eq =ξ k eq (m+m), under the harmonic excitation F(t)=F A sin(ω F t), we obtain the well known solution y y ξ ω t y y d ω d yωd + ξω + ξω yt () = y + e cos ω t arctan + F A ξδ + sin ωf t arctan k eq ( ) ( ) δ + ξδ δ (3) where: y = F l, FA ϕ y = l A E T A AEH ϕδ( δ ) ξ= δ +δ +ϕ H ( ) ω =ω ( + ϕ )( +η), ω d =ω ξ (33) δ =δ ( +ϕ )( + η ) k eq k =, + ϕ ϕδ( δ ) k ceq = ( m+ M) δ ( +δ +ϕ ) +ϕ The amplitude of the steady state response of the single RA system with one lumped element of mass can be written in the form FA A = k eq ( δ ) + ( ξδ ) (34) The static deflection F A /k eq is multiplied by the dynamic magnification factor = ( δ ) + ( ξδ ) (35) 4. THE VERIFICATION OF RESULTS OF VISCO-ELASTIC RA OSCILLATOR ON LOW CARBON STEEL BAR Extensive work in the several years has been done on analogy; see Milašinović [6] and [7]. The test material was low carbon steel. The low carbon steels are an economic yet effective solution for the seismic retrofit of highway bridges. These readily available materials have a yielding stress of 3 ksi (.64Mpa).

12 34.. MILAŠINOVIĆ The proportional stress (σ P = 4MPa) or reaction-stress of the clamped bar under axially fatigue process is obtained using physical characteristics of low carbon steel only, like: specific heat, coefficient of linear thermal expansion and mass density. Elasticity stress (σ E = 87MPa) of the bar under compression is obtained as Euler s critical stress from which we have that elasticity stress also becomes dependent upon the dimensions of the bar (its length and diameter) and thus is no more a physical characteristics of the material only. In ductile materials like metals, some difference between the point of proportionality and the point of elasticity produce the deviation from perfect elasticity with dissipation of mechanical energy through quasi-viscous flow or visco-elastic creep. This is the reason of the drop in the stress-strain curve in the average σ-ε diagram with upper and lower-yield points, see Milašinović [7]. In order to clarify the influence the relative ratios: η = M/m, φ * and δ = T ω F (ω F = πf F ) on the viscous damping ratio and the dynamic magnification factor, the steel bar was analyzed under four values of η (,,, 3683), φ * = and excitation frequency, f F = 5Hz (δ=.94). Table gives values of viscous damping ratios and dynamic magnification factors, calculated from Eqs. (6) and (35) for relative frequencies varying from δ =. to δ = and for the maximal mass ratio, η = c/ccr,6,5,4,3,, -, -, -,3 -,4 -,5 -,6,, logδ* M/m= M/m= M/m= M/m=3683 Fig. 7. The dependence of the viscous damping ratio on relative frequencies δ * Figure 7 gives variations of the viscous damping ratio on the relative frequency δ * for four values of η (,,, 3683). It can be seen from Fig. 7 that if the mass ratio is smaller, the viscous damping ratio will be larger and also that if the relative frequency is larger, the viscous damping ratio will be smaller. It is interesting to note that for η =, positive viscous damping ratio decreasing from.6 to.4. This confirms similar experimental observations made on steel, which show.>c / c cr >.3.

13 The Rheological-ynamical Harmonic Oscillator 35 Table. Values of viscous damping ratios and dynamic magnification factors δ=t ω F Prototype T =.967s f F =5Hz ω F =πf F =94.5rad/s δ=.94,,,3,4,5,6,7,8,9,94,,,3,4,5,6,7,8,9,,,3,4,5,6,7,8, δ * = δ (+φ * )(+η) Prototype σ P =4Mpa F A =(9 π/4)4=46n M=F A /g=44kg η=m/m=44/.4=3683 δ * =.9584,58,57,35385,454,5564,6377,735899,848,94656,9584,585,57 3, ,54 5,5645 6,377 7, ,479 9,46564,585,57 3, ,54 5, ,77 73, ,79 94,6564 5,85,57 35,3855 4,54 55,645 63,77 735, ,79 946,564 5,85,57 353,855 45,4 556,45 637,7 7358,994 84,79 946,564 5,85 ξ,336,33,855,6,95,99,454,98,33,6 -,33 -,84 -,836 -,587 -,8437 -,84 -,686 -,63 -,799 -,346 -,384-3,589-5,3995-8,8388 -,69-4,7574-8,485 -,3477-6,797-6,7-78, ,588-93, ,44-99,63 -,4 -,374 -,66-4,346-4,779-4,93-5, -5,4-5,64-5,79-5,9-5,97,75,465,45,485,3878,66763,83 3, ,5496, 9,5549,933,749,5999,3753,576,8793,43,77,9,48,988,549,346,37,7,9, 8,9E-5,96E-5 8,99E-6 5,E-6 3,4E-6,4E-6,78E-6,38E-6,9E-6 8,88E-7,5E-7 E-7 5,65E-8 3,6E-8,5E-8,85E-8,4E-8,E-8 9,5E-9

14 36.. MILAŠINOVIĆ The dynamic magnification factors for the above four mass ratios are shown in Fig. 8. It is seen that all dynamic magnification factors increase with respect to relative frequencies.<δ *. Note, however, that all factors decreasing for negative values of ξ (δ * ). The negative values of ξ are examined only as a hypothetical case because the complementary solution of the response associated with negative viscous damping ratio would not die away. The system is unstable because of the factor e ξ ω t. 3,,, logδ* M/m= M/m= M/m= M/m=3683 Fig. 8. The dependence of the dynamic magnification factor on relative frequencies δ * The corresponding value for maximal mass ratio η = 3683 was =., while for η = it is found that = 37.. Note however that the dynamic magnification factor has maximal value of = 3.76 for η =. In this case the relative frequency δ * = is very close to one, which is the case of resonance. Resonance occurs when the frequency of the excitation ω F is equal to the natural frequency of the RA system ω = ω ( +ϕ )( + η) Finally it is useful to analyze the dependences of the viscous damping ratio (Fig. 9) and the dynamic magnification factor (Fig. ) on the relative frequency δ of the RA model, becouse this frequency is independent of mass ratio η.

15 The Rheological-ynamical Harmonic Oscillator 37, logc/ccr,,,,, logδ M/m= M/m= M/m= M/m=3683 Fig. 9. The dependence of the viscous damping ratio on relative frequencies δ * log..... logδ M/m= M/m= M/m= M/m=3683 Fig.. The dependence of the dynamic magnification factor on relative frequencies δ The exemplary calculations for one degree of freedom visco-elastic RA system are performed based on the formulas derived in the former chapter. Calculations are done for the harmonic excitation F() t = F + FAsin( ω Ft) where F = σ A, (σ =E H ε ) Solution under constant permanent force F has the form of y + yξω ξω t y + yξω F d ω k d y ωd eq yt () = y + e cos ω t arctan + The first part of the solution, the complementary, usually called the transient, will eventually die out. Also, the rate of decay and natural frequency of the system depend on the system parameters only, while the amplitude of vibration and phase angle are determined by the initial conditions.

16 38.. MILAŠINOVIĆ Some values of parameters of a single RA system with one lumped element of mass of steel bar appear in Table. the results of other parameters that must be used for computational results of longitudinal vibration are presented in Table. The bar is loaded with cyclic sinusoidal load in symmetrical cycle and frequency f = 5Hz (ω F = πf = 94.5rad/s). Table. The RA model and single RA system parameters Φ=9mm l =5mm σ =4MPa F =(9 π/4)x4 =46N k=e H A/l =9.8 MN/m m=ρal =.4kg M=F /g=44kg ω F =πf F =π5= =94.5rad/s RA model Single RA system Characteristic time T =l / E H /ρ[s].967 Creep coefficient φ *. Initial strain ε =σ /E H.676 y =ε l =F /k [m].338 (displacement) Initial strain rate dε /dt=σ /λ = 3.98 dy /dt=(dε /dt)l 6.99 (velocity) =ε φ * /T [/s] [m/s] Mass ratio η=m/m 3683 Angular frequency ω=/t [/s] 34 ω * =ω/ (+φ * )(+η) [rad/s] Relative frequency δ=ω F T.94 δ * =ω F /ω *.9584 Viscous damping ξ. ξ.6 ratio amped frequency ω * d = ω * - ξ The solution with parameters from Table takes the form of y y + y ξω + = ωd =.7m * ξω t t.5575 t e = e = e y + yξω arctan = arctan.566 y = rad ωd F 46 ( +ϕ ) = ( + ) = =.m 6 k 9.8 yt e t.5575 t ( ) =.7 cos( ) +.

17 The Rheological-ynamical Harmonic Oscillator y(t) [m] ω dxt [rad] Fig.. The transient vibration amplitudes with permanent force F (<t<.4s) Figure (<t<.4s) and Fig. (8<t<s) display the transient vibration amplitudes. From the beginning (Fig. ), the RA oscillator moves periodically and then alight to a constant regime (Fig. ). The complementary solution die out after 3s ( = 95rad) y(t) [m] ω dxt [rad] Fig.. The transient vibration amplitudes with permanent force F (8<t<s) Let us consider the periodically force F(t) (F A =(9 π / 4)x4=46N), which is applied slowly F() t = FAsin( ω Ft) where f=5hz (ω F =πf = 94.5rad/s). The solution of the RA oscillator has the form of

18 4.. MILAŠINOVIĆ where F A ξδ yt () = sin ωf t arctan k eq ( ) ( ) δ + ξδ δ F A ( +ϕ ) k ( δ ) + ( ξδ ) (.9584 ) + ( ) = =.4 3 =.37m ξδ arctan arctan.68rad = = δ.9584 Finally, the solution is yt ( ) =.37 sin( t.68) Figure 3 display the steady vibration amplitudes, which are found for three types of vibrations: the RA visco-elastic, the RA oscillator (elastic solution) and the RA oscillator (visco-elastic solution). The RA visco-elastic results are obtained by solving Eq. 6. y(t)[m],5,3,,9,7,5,3, -, -,3 -,5 -,7 -,9 -, -,3 -,5,,4,6,8 Time[s] RA oscillator: elastic solution RA visco-elastic solution RA oscillator: visco-elastic solution Fig. 3. Steady state response of forced longitudinal vibration of prototype The elastic solution according to the RA oscillator has the maximal amplitude of.4m, which is greater for. (.4/.338=.) from maximal amplitude under constant force F. It is in accordance with dynamic magnification factor (see Fig., M/m = 3683). The RA oscillator give the visco-elastic maximal amplitude of.37m, which is greater for. (.37/.4=.) from maximal amplitude obtained using the RA visco-elastic model (see Eq. 6).

19 The Rheological-ynamical Harmonic Oscillator 4 5. VISCO-ELASTO-PLASTIC RA OSCILLATOR The only difference to be noted for a bar that is elasically deformed and one that is inelastically deformed is the stress-strain relation for the material of the bar. Inelasticity theory would be greatly simplified if one function could be found to approximate the stress-strain diagram over both the elastic and inelastic range. This function may be obtained as isochronous octahedral shearing stress-octahedral shearing strain RA diagram (see [7], Section 8). Since about 997 (Milašinović [4]), mathematical-phisycal analogy between the visco-elasto-plastic rheological model and the dynamical model was proposed to an explicit form to predicting many different inelastic and time dependent problems related to prismatic bars, such as buckling, fatigue etc. According to the principle of analogy, the visco-elasto-plastic oscilator has the same solution as the visco-elastic one (see Chapter 3), if we chose ϕ =ϕ ( + i), i =,,,... (36) where i is the level for visco-plastic yielding at any stage of dynamic equilibrium. Also, the forces: F max, F A and F Y are not independent. The following expressions may be derived, for the known: F P, F E, F and φ * (see [7], Section 5.) F F ( + iϕ ) + F ( +ϕ )[ i( + iϕ ) ], F iϕ ( + iϕ ) () i P E max = r () i () i F Fmax () i Fmax =, F () i () i () i A Fmax Y FP = + F E, ϕ r =. (37) In this way the RA oscillator has the advantage that all the calculations of the strainto-stress type (which is the way usually needed in computer calculation) may be carried out explicitly, i.e., without the need for any step-by-step or iterative integration procedure in the loading steps. 6. CONCLUSION In this paper the rheological-dynamical theory of the analytical dynamics of discrete visco-elasto-plastic system is presented. The coupled initial conditions of the stress-strain state of the rheological visco-elasto-plastic model are applied for the study of single RA system. The motion of the RA system is described by the equation with the highest derivative of the fourth order with respect to time t. In the case of linear visco-elasticity, RA model has the same natural frequency as a simple single-degree-of-freedom spring mass system and because of that a new single RA system with one lumped element of mass is formed. Also, the RA model has the same phase angle as a simple single-degree-of-freedom spring mass system with damping in the linear steady state vibration and from that the viscous damping ratio of the RA system is obtained. The viscous damping ratio is a function of various relative parameters of the RA model like: creep coefficient φ *, mass ratio η and relative frequency δ. It is known that for the base isolation, the fundamental frequency of the whole structure is dominated by the natural frequency of the base isolator. Consequently, it is very important to evaluate the viscous damping ratio for the isolator, which is the function of

20 4.. MILAŠINOVIĆ mentioned relative parameters. The purpose of the energy dissipation devices is to increase the natural period of the structure so that the acceleration response of the structure is decreased during the gust. The effective period ratio increases rapidly as the applied mass ratio increases. The visco-elastic RA oscillator is confirmed and validated by some of the published experimental data on the monotonic and the cyclic loading of mild-steel. While studying the discrete visco-elasto-plastic RA system in Chapter, the fourth order differential equation with the small parameter before the highest derivative is given as a basic result. For such equations, as known according to A. I. Tykhonov theorems, the solution of singular equation does not always converge to initial. So we cannot approve that the third order differential equation describes the initial system correctly. Nevertheless, visco-elasto-plastic RA oscillator has been generated from mathematical-phisycal analogy, through an elegant framework that bypasses most of the questions arising from the vibration of visco-elasto-plastic bodies. Acknowledgement: This work was supported by the Ministry of Science and Environmental Protection - Republic of Serbia, Project: REFERENCES. Clough, R. W. and Penzien, J., ynamics of Structures, McGraw-Hill, New York, 975. Goroshko, O. A. and Hedrih (Stevanivić),., Analytical ynamics of iscrete Hereditary Systems, University of Niš, iew-niš (in Serbian), Love, A. E. H., A Treatise on the Mathematical Theory of Elasticity, 4 th Ed., over Publications, New York, Milašinović,.., The Finite Strip Method in Computational Mechanics, Faculties of Civil Engineering: Subotica, Budampest, Belgrade, Milašinović,.., Rheological-dynamical analogy: prediction of buckling curves of columns, International Journal of Solids and Structures 37(9) pp ,. 6. Milašinović,.., Rheological-dynamical analogy: modeling of fatigue behavior, International Journal of Solids and Structures 4() pp. 8-7, Milašinović,.., Rheological-dynamical analogy: visco-elasto-plastic behavior of metallic bars, International Journal of Solids and Structures 4(6-7) pp , Milašinović,.., Rheological-dynamical Analogy: Frequency ependence of the System Parameters of Internally amped Bars, Multidiscipline Modeling in Materials and Structures () pp. 43-7, Milašinović,.., Rheological-ynamical Theory of Visco-Elasto-Plasticity and Fatigue: Part, Multidiscipline Modeling in Materials and Structures () pp. -9, 6.. Reiner, M., Rhe'ologie the'orique, unod, Paris, France, 955. REOLOŠO-INAMIČI HARMONIJSI OSCILATOR ragan. Milašinović U radu je predstavljena reološko-dinamička teorija analitičke dinamike diskretnih visko-elastoplastičnih sistema. Reološko-dinamička analogija (RA) je razvijena na bazi matematičko-fizičke analogije između reološkog modela i diskretnog dinamičkog modela sa viskoznim prigušenjem sa ciljem da bude korištena u analizi neelastičnog deformisanja materijala i konstrukcija. U ovoj presentaciji, povezani početni uslovi naponsko-deformacijskog stanja reološkog visko-elastoplastičnog modela su primjenjeni u studiji diskretnog dinamičkog sistema. Mehanički sistemi kao što su kablovi, štapovi, grede, ploče i drugi, čije su mase i elastične sile raspoređeni različito od

21 The Rheological-ynamical Harmonic Oscillator 43 oprugama spojenih koncentrisanih masa, spadaju u grupu vibracija kontinualnih sistema. Ovi sistemi imaju beskonačno veliki broj djelića i zahtjevaju beskonačno veliki broj koordinata za specificiranje njihove konfiguracije. Ovaj tekst daje primjer i ilustruje način kako diskretan model može biti izveden, korištenjem principa analogije, iz specifičnih granica kontinualnog modela. Ova tehnika je korisna jer diskretan model u konceptualnom i proračunskom smislu često je mnogo jednostavniji za analizu od kontinualnog modela. ljučne reči: RA analogija, RA harmonijski oscilator, viskozno relativno prigušenje