Access the full text.
Sign up today, get DeepDyve free for 14 days.
References for this paper are not available at this time. We will be adding them shortly, thank you for your patience.
This work addresses friction-induced modal interactions in jointed structures, and their eﬀects on the passive mitigation of vibrations by means of friction damping. Under the condition of (nearly) commensurable natural frequencies, the nonlinear character of friction can cause so-called nonlinear modal inter- actions. If harmonic forcing near the natural frequency of a speciﬁc mode is applied, for instance, another mode may be excited due to nonlinear energy transfer and thus contribute considerably to the vibration response. We inves- tigate how this phenomenon aﬀects the performance of friction damping. To this end, we study the steady-state, periodic forced vibrations of a system of two beams connected via a local mechanical friction joint. The system can be tuned to continuously adjust the ratio between the ﬁrst two natural frequencies in the range around the 1 : 3 internal resonance, in order to trigger or suppress the emergence of modal interactions. Due to the re-distribution of the vibration energy, the vibration level can in fact be reduced in certain situations. However, in other situations, the multi-harmonic character of the vibration has detrimen- tal eﬀects on the eﬀective damping provided by the friction joint. The resulting response level can be signiﬁcantly larger than in the absence of modal interac- Corresponding author Email addresses: krack@ila.uni-stuttgart.de (Malte Krack), lbergman@illinois.edu (Lawrence A. Bergman), avakakis@illinois.edu (Alexander F. Vakakis) Preprint submitted to Elsevier arXiv:2101.03232v1 [nlin.PS] 28 Dec 2020 tions. Moreover, it is shown that the vibration behavior is highly sensitive in the neighborhood of internal resonances. It is thus concluded that the condition of internal resonance should be avoided in the design of friction-damped systems. Keywords: friction damping, modal interactions, internal resonance, jointed structures 1. Introduction Friction damping is a well-known means of achieving vibration reduction. The damping eﬀect is due to the dissipative character of dry friction occurring in mechanical joints. For this purpose, friction interfaces may be either newly introduced to the structure or may already exist e. g. in the form of bolted or riveted joints. Friction damping is particularly suited for the passive vibration reduction of lightly damped ﬂexible structures. Various applications can be found in the ﬁeld of aerospace structures, combustion engines or turbomachin- ery blades [1, 2]. Dry friction is a nonlinear phenomenon: Depending on the vibration level of a jointed structure, the behavior in the local contact interfaces may range from sticking via micro-slip to macro-sliding. Hence, dry friction may induce the com- mon nonlinear phenomena. In the context of steady-state harmonically-forced vibrations, these phenomena include: (a) The dependence of the resonance fre- quency and the eﬀective damping on the vibration level, (b) the loss of stability of the fundamental periodic response giving rise to sub-harmonic, quasi-periodic or chaotic vibrations, and (c) the occurrence of so-called nonlinear modal inter- actions. Nonlinear modal interactions refer to energy exchanges between two or more of a structure’s modes of vibration caused by nonlinear eﬀects. Thus, a mode that is not directly forced by an external source, may be excited indirectly by the action of nonlinear forces, and therefore contribute signiﬁcantly to the overall vibration response. This occurs when the natural frequency of the indirectly forced mode is in a rational relation with the natural frequency of a directly 2 forced mode [3], i. e., when the two modes are in internal resonance. Under this condition, a suitable nonlinearity can initiate mode mixing caused by energy transfers between the modes, even if their frequencies are well separated. Note that the necessary condition of internal resonance does not have to hold for the natural frequencies of the linearized system, but this condition can be satisﬁed for larger vibration levels of the nonlinear system when the natural frequencies change accordingly. In the context of friction damping, numerical and experimental evidence of non- linear modal interactions was reported, e. g., in [4, 5, 6, 7, 8]. Speciﬁcally, an important eﬀect on the steady-state forced vibration response in the excitation frequency range around a particular natural frequency is commonly reported: The response of the indirectly forced mode assumes a local maximum. At the same time, the expected pronounced peak in the amplitude-frequency curve of the directly forced mode transforms to a double peak with a local minimum in between two local maxima. This phenomenon appears to be more profound if the joint’s normal load ﬂuctuates and the joint undergoes lift oﬀ/impacts as a consequence of vibrations [6, 7, 8]. Ferri et al. explored the energy transfers from low to high frequencies of friction damped systems under the condition of internal resonances, both in the steady-state forced response and the free vibrations [9, 10, 11]. Although nonlinear modal interactions are known to occur in friction-damped systems, their inﬂuence on the vibration reduction performance is unknown and yet has to be studied. It should be noted that nonlinear modal interactions are known to play the essential role in the vibration control concept of targeted energy transfer [12]. The key of this concept is that the vibration energy is passively and irreversibly transferred from a primary mode to either a strongly nonlinear attachment or a higher-frequency modes, where it is dissipated more eﬃciently. The purpose of the present study is thus to investigate how nonlinear modal interactions aﬀect the eﬃcacy of friction damping. To this end, we study a system composed of two cantilevered beams connected via a friction joint, as presented in Sec. 2. The system can be tuned so that the beams’ ﬁrst bending 3 Figure 1: Investigated model: (a) two beams connected via friction joint, (b) elastic Coulomb friction model; E = 6 · 10 = E , ρ = 2500 = ρ , l = 1.0 = l , b = 0.02 = b , h = 0.02 1 2 1 2 1 2 1 2 1 modes are in 1 : 3 internal resonance. It is demonstrated in Sec. 3 how and under what conditions the friction nonlinearity induces the intended nonlinear modal interactions. The consequences for the vibration reduction performance are quantiﬁed and explained in Sec. 4. This article ends with concluding re- marks in Sec. 5. 2. Model description 2.1. Problem formulation Consider the system of two beams connected via a friction joint, depicted in Fig. 1a. Bending vibrations with the deﬂection w(x, t) in y-direction are con- sidered. For convenience, the problem is divided into two sections, one for each beam, with local coordinates x and x , respectively. In accordance with the 1 2 classical Euler-Bernoulli beam theory , the equations of motion of the system small deﬂections, plane sections remain plane, constant length neutral axis, linear-elastic material 4 can be written as, ′′ ′′ (E I w (x , t)) + ρ A w¨ (x , t) = 0 , i = 1, 2 , 0 < x < ℓ (1) i i i i i i i i ′ ′ w (0, t) = 0 = w (0, t) , w (ℓ, t) = 0 = w (ℓ, t) (2) 1 1 2 2 ′ ′ ′′ ′′ ′′ ′′ (E I w (ℓ, t)) + f (t) = f = − (E I w (0, t)) , w (ℓ, t) = 0 = w (0, t)(3) 1 1 exc j 2 2 1 2 1 2 f = f [(w (ℓ) − w (0))] (4) j c 1 2 k dw |f + k dw | ≤ µN c c c c c f [w ] : df = , (5) c c c 0 otherwise where Eq. (3) denotes the force balance at the joint coupling the two beams, and Eqs. (4)-(5) provide the constitutive force-deformation law for the frictional con- nection. Note that local spatial coordinates are used for each of the two beams. In Eqs. (1)-(5), prime denotes diﬀerentiation with respect to local coordinate x or x , and overdot denotes diﬀerentiation with respect to time, and A = b h 2 i i i b h and I = are the area and the moment of inertia of the rectangular cross section, respectively. Note that the elastic Coulomb friction law is considered at the frictional connection (governed by Eq. (5)). The elastic Coulomb law relates the friction force f and the relative displacement w . This relationship depends c c on the joint stiﬀness k and the limiting friction force µN , which is the product of the friction coeﬃcient 0 ≤ µ ≤ 1 and the joint normal force N > 0. The hysteretic dependence of the friction force on the relative joint deformation w is expressed by the operator notation [w ]. An example for a simple hysteresis cycle is depicted in Fig. 1b. The joint normal force is well known to have an essential inﬂuence on the friction damping performance. Therefore, its value was varied during this work, as detailed later. 2.2. Spatial discretization System (1)-(5) is spatially discretized by a ﬁnite element procedure. Each beam section is equidistantly divided into N = N = 20 elements. Standard e,1 e,2 beam elements with consistent element matrices were utilized with two nodes per element and two degrees of freedom per node, one for the deﬂection and 5 one for the inclination. This results in a set of second-order ordinary diﬀerential equations governing the vector of generalized coordinates u(t), ¨ ˙ Mu(t) + C u(t) + Ku(t) + f [u(t)] = p e cos Ωt . (6) mod c exc w (ℓ) The total number of generalized coordinates is N = 2 (N + N ) = 40. A har- e,1 e,2 monic external forcing of magnitude p and frequency Ω is considered in this exc study. e is the unit vector with all components zero except the one associ- w (ℓ) ated with the displacement w (ℓ) at location x = ℓ. A linear viscous damping 1 1 term was introduced in Eq. (6) to take into account the material damping eﬀect. It is deﬁned as modal damping, Φ C Φ = diag (2D ω ) , (7) mod i i in the modal space associated to the linear system without joint. The quantities ω and D are the i-th modal frequency and damping ratio, respectively, and i i h i Φ = is the matrix of the N mass-normalized eigenvectors deﬁned φ , ·,φ 1 N in terms of the usual orthonormality conditions, T T 2 Φ MΦ = I , Φ KΦ = diag ω . (8) A modal damping ratio of D = 1%, i = 1, . . . , N was speciﬁed for all modes, if not explicitly stated otherwise. 2.3. Tuning The left beam and the right beam have identical material and geometrical properties except for their height h , as listed below Fig. 1. While the left beam’s height h remained constant throughout this study, the height h of the right 1 2 beam was varied in order to vary the natural frequencies of its modes, i. e., its component modes in the absence of the joint. Note that all natural frequencies of the beam are proportional to its height. Thus, any ratio between speciﬁc natural frequencies of the left and the right beam can be realized by continuously varying h . In this study, we focused on the 1 : 3 internal resonance between 6 the ﬁrst bending mode of the left beam and the ﬁrst bending mode of the right beam. To this end, h has to be speciﬁed to be close to 3h . For convenience, 2 1 a detuning parameter Δ is introduced, h (Δ) = 3h (1 + Δ) . (9) 2 1 For Δ = 0, the corresponding modes of the left and the right beam are in 1 : 3 internal resonance, while for Δ > 0 (Δ < 0) the natural frequency ratio between the right and the left beam is larger (smaller) than three. 2.4. Model reduction To reduce the computational eﬀort of the nonlinear dynamic analysis, a conventional modal truncation was utilized. Based on a preliminary convergence h i study, the six modes Φ = associated with the six lowest red φ . . . φ 1 6 natural frequencies of the system without joint served as reduced basis. In the neighborhood of the 1 : 3 internal resonance, these six modes are the leading four modes of the left beam, and the leading two modes of the right one. The associated natural frequencies are ω = 1.0, ω = 6.27, ω = 17.5 and ω = 34.4 1 3 4 6 for the left beam, and, assuming h = 3h , ω = 3.0 and ω = 18.8 for the right 2 1 2 5 beam. Eq. (6) was projected accordingly, resulting in the following reduced set of equations of motion, 2 T T q¨ (t)+2D ω q˙ (t)+ω q (t)+φ f [Φ q(t)] = p φ e cos Ωt , i = 1, . . . , 6 , i i i i i c red exc w (ℓ) i i i 1 (10) h i in terms of the modal coordinates q = . q . . . q 1 6 2.5. Inﬂuence of the linearized joint behavior N µ For suﬃciently small magnitudes |w | ≤ of the joint displacement coor- dinate w , the friction contact is only sticking and the joint behaves as a linear spring with stiﬀness k . The larger k is, the stronger the coupling between the c c left and the right beam becomes. For k = 0, the two beams are decoupled, and the system modes are isolated to either the left or the right beam, see Fig. 2a. For increasing k , the modes of the linearized system, transform from isolated 7 Figure 2: Vibrational deﬂection shapes of the ﬁrst and second linear modes: (a) decoupled 3EI with k = 0, ω = 1.00, ω = 3.00, (b) coupled with k = k := 1.25 , ω = 1.45, c 1 2 c weak 1 c 3 3EI ω = 3.07, (c) coupled with k = k strong := 125 · , ω = 2.56, ω = 4.70 2 c 1 2 c 3 modes to system modes with increasingly integral character; that is, the vibra- tion energy is more evenly distributed among the two beams, see Fig. 2b-c. In the nonlinear regime, a larger k also leads to a stronger inﬂuence of the nonlin- earity on the overall system response. In this study, two diﬀerent values of k 3EI were considered: The case of relatively a weak coupling with k weak = 1.25 , 3EI strong and the case of a relatively strong coupling with k = 125 · . It was found that the mechanism based on nonlinear modal interactions and most of the qualitative results obtained in this work are identical for both weak and strong coupling. Accordingly, most results are presented only for the case of 3EI weak coupling with k = 1.25 weak . Only in Sec. 4, where the resulting fric- tion damping performance is analyzed, results are also presented for the case of strong coupling, owing to the observed important diﬀerences between the weakly and strongly coupled cases. 3. Friction-induced nonlinear modal interactions The condition of internal resonance is necessary but not suﬃcient for the oc- currence of an interaction between the corresponding modes. As is demonstrated in this section, the friction nonlinearity can in fact cause this modal interaction and lead to pronounced eﬀects on the forced response. In Subsec. 3.1, it is inves- tigated in what dynamic regime this interaction takes place and where it is most prominent. In Subsec. 3.2, the underlying vibration mechanism is analyzed in terms of the dynamic power ﬂow through the system. Throughout this section, 8 1.6 1.6 1.6 1.6 ∗ ∗ Figure 3: Normalized steady-state response of the modal amplitudes qˆ , qˆ of the left 1 2 and the right beam’s ﬁrst bending mode: (a)-(b) results for left and right beam for p ∈{0.02,0.07,0.2,0.5,0.7,0.9,1.0,1.2,1.4,3.5 exc variation of the excitation level, μN = 1.0, , (c)-(d) results for left and right beam for variation of the limit friction force, μN ∈ {0.3, 0.75, 0.9, 1.0, 1.2, 1.5, 2.0, 4.5, 15, 45}, p = 1.0 exc the detuning parameter is set to Δ = 0; i. e., the ﬁrst bending modes of the left and the right beam are exactly in internal resonance condition. 3.1. Inﬂuence of excitation level and limit friction force on the intensity of the nonlinear modal interactions The nonlinear, steady-state forced response is considered in the excitation frequency range in the neighborhood of ω , i. e., the natural frequency of the left beam’s ﬁrst bending mode. The vibration response was computed numerically using the well-established harmonic balance method. In the computations, only 9 the ﬁrst ﬁve harmonics were retained in the Fourier expansion of the response of each beam. This number was deemed suﬃcient in accordance with a prelimi- nary convergence study, where the results were compared with direct numerical time integration. Note that the highest considered natural frequency ω is much larger than the ﬁfth harmonic of the driving frequency Ω, 5Ω ≈ 5.0 ≪ 34.4 ≈ ω , cf. Subsec. 2.4. Hence, the discarded modes are far from the driving frequency and its ﬁrst ﬁve harmonics. The excitation force p and limit friction force µN exc were varied over a wide range of values. The results for the modal amplitudes qˆ , qˆ of the left and the right beam’s ﬁrst bending mode are depicted in Fig. 3 1 2 in terms of their ﬁrst and third harmonic components. The third harmonic of q is so small that it lies outside the depicted amplitude range in Fig. 3a and c. The results are normalized with respect to the resonant case for p = 1.0 exc and sticking contact conditions. More speciﬁcally, the depicted quantities were divided by the magnitude of the ﬁrst harmonic of q for these conditions. Consider ﬁrst the results depicted in Fig. 3a-b for varying excitation level. As stated above, the nonlinear friction joint acts as a spring k for small vibrations, so in that limit the system exhibits linear behavior. Since friction damping is absent in this case, the resonance peak is relatively sharp and located at the corresponding natural frequency Ω = 1.45ω of the ﬁrst system mode. On the other hand, for large vibrations, the friction joint is mostly slipping, but the limited friction force becomes small compared the other forces in the system. Hence, for large vibrations the inﬂuence of the friction joint decreases, and the nonlinear system behavior approaches the linear behavior of the system with- out a joint. Thus, the resonance peak sharpens again and resonance frequency approaches the natural frequency Ω = ω . In between these (quasi-)linear lim- iting cases, the friction nonlinearity causes a considerable damping eﬀect. As a consequence, the resonance peak broadens signiﬁcantly, in particular for lower to moderate excitation levels. This behavior is well known for friction-damped systems. At a certain range of the excitation level, considerable eﬀects of nonlinear modal interactions can be ascertained in our results: The third harmonic of the right 10 beam’s response assumes a prominent maximum. In this regime, the third har- monic clearly dominates the response of the right beam. In the same range, the peak of the left beam’s response transforms into a double-peak. More specif- ically, the amplitude-frequency curve of the fundamental harmonic component assumes two local maxima with an intermediate local minimum. This phe- nomenon is a characteristic feature of nonlinear systems in internal resonance. Provided that the elastic Coulomb nonlinearity is the only nonlinearity present in the system, the following scaling property relates the response u(t) of the nonlinear system, the excitation level p and the limit friction force µN [13], exc u| (t) = γ u| (t) , (11) γN µ ,γp N µ ,p exc exc for any positive real constant γ. If the response u(t) is known for a set of val- ues for µN and p , Eq. (11) allows to determine the response for accordingly exc scaled parameters γµN and γp without the need for re-computation. As a exc consequence, the variation of the excitation level p has a similar eﬀect as exc the variation of the limit friction force µN . Hence, the amplitude-frequency curves depicted in Fig. 3c-d are strictly correlated with the ones in Fig. 3a-b, and they are merely presented here to more clearly demonstrate the inﬂuences of the individual parameters on the steady-state forced responses. If the ex- pected excitation level is approximately constant and known a priori, the limit friction force µN can be designed to minimize the vibration level. In this case, a reduction by a factor ten can be achieved, compared to the linear case with only sticking joint. The vibration level is less sensitive for values µN larger than this nominal optimum, which is usually preferred in order to improve the reliability with respect to uncertainties [13]. It should be noted that the qualitative dependence of the maximum vibration level on p and µN is not changed due to the nonlinear modal interactions exc considered here. Still, the regime of pronounced inﬂuence of the modal inter- actions is close to the optimum design point, and therefore appears relevant to the design process. 11 external loading material material damping damping transfer transfer joint-right beam left beam-joint friction 1.6 damping (a) (b) 1.6 1.6 (d) (c) Figure 4: Energetic analysis of the steady-state vibration response: (a) Schematic illustration of the power ﬂow, (b) power balance in left beam, (c) power balance in joint, (d) power balance in right beam 3.2. Analysis of the dynamical mechanism governing the nonlinear modal inter- actions In this subsection, we examine the vibration mechanism that leads to emer- gence of the phenomena associated with nonlinear modal interactions. To this end, we consider the steady-state forced response for a speciﬁc excitation level p and a speciﬁc limit friction force µN such that µN/p = 1.0. For each exc exc system component, i. e., the left beam, the joint and the right beam, we analyze the following averaged mechanical powers: The power input from the external forcing, the power dissipated by material or friction damping and the power transferred to or received from components. These quantities are illustrated in Fig. 4. A positive value indicates a power input to, whereas a negative value 12 indicates a power output from the considered component. The orientation of the arrows in the schematic overview Fig. 4a takes into account the actual sign of the power. As can be seen in Fig. 4b, the left beam is driven by the power of the external forcing. Some of this power is directly dissipated by material damping as a consequence of the structural vibration, but the largest part of the input power is transmitted directly to the frictional joint. In the joint, most of this power is dissipated by dry sliding friction, see Fig. 4c. Owing to the nonlinear char- acter of friction, the periodic joint force consists not only of the fundamental frequency component but also of higher frequency components, among which the third harmonic is particularly large. In the neighborhood of the internal resonance, the right beam exhibits a comparatively high dynamic compliance with respect to the third harmonic of the excitation frequency. Hence, power is also transmitted from the joint to the right beam. Although this power is small compared to the power dissipated in the joint, it is still large enough to maintain the vibration of the right beam. The power transmitted to the right beam is, of course, dissipated by the inherent material damping. From a structural dynamics point of view, two eﬀects of the nonlinear modal interactions are of particular importance: (a) The spatial re-distribution of the vibration energy within the system: While in the absence of modal interactions, the vibration is localized in the left beam, the vibration energy is more evenly distributed between the two beams when nonlinear modal interactions occur. (b) The temporal re-distribution of the vibration energy: While in the absence of modal interactions, the fundamental frequency component clearly dominates the vibration response, considerable energy is transferred to higher frequency components in the case of nonlinear modal interactions. 13 4. Resulting friction damping performance and its dependence on frequency detuning In this section, we investigate how a frequency detuning of the system around the considered 1 : 3 internal resonance aﬀects the friction damping performance. To this end, the parameter Δ is varied, which controls the height and therefore the stiﬀness of the right beam. From a vibration control perspective, the ob- jective is to maximize the structural reliability and therefore to minimize the mechanical stress in the system. We therefore consider the normal stress σ in xx the axial direction in either of the coupled beams. At a given spatial location within the structure, the stress varies with time t. The steady-state, periodic forced response generally exhibits multi-harmonic character. We thus consider its root-mean-square (RMS) value of the stress, u 1 σˆ (x , y) = σ (x , y, t) dt i = 1, 2 . (12) i i xx where x is the local axial coordinate for each beam, y is the transverse coordi- 2π nate, and T = is the fundamental period of oscillation. ′′ In the absence of normal loading, σ (x , y, t) = Ew (x , t)y. The critical xx i i transversal location is thus always the top or the bottom of the rectangular cross-sectional area, y = ±h /2, of each beam. The critical stress within each beam is thus σˆ = max σˆ x , , i = 1, 2 . (13) crit,i i 0≤x ≤ℓ 2 The critical cross-sectional area is typically located at the clamped end, i. e., at x = 0 and x = ℓ for the right and the left beam, respectively. Hence, the 1 2 critical stress within the system is the largest critical component stress, σˆ = max (σˆ , σˆ ) . (14) crit crit,1 crit,2 Owing to the comparatively small vibration level of the right beam, σˆ is crit always attained in the left beam, that is σˆ = σˆ in all cases presented in crit crit,1 this study. In the case of steady-state vibrations in the presence of harmonic 14 external forcing, the critical stress σˆ is a function of the excitation frequency crit Ω, and the maximum σˆ within the considered frequency range is denoted by crit σˆ , crit,max σˆ = max σˆ (Ω) . (15) crit,max crit This quantity is used to assess the eﬀect of frictional damping on the fatigue of the coupled structure, and is thus deﬁned as the objective function in the subsequent numerical studies. For suﬃciently large |Δ|, the nonlinear modal interactions with the right beam are suppressed, and the response level approaches an (almost) constant value. This detuned case serves as a reference. Stress values are accordingly normalized by the maximum (resonant) stress σˆ in this case. As it turns out, crit,max,detuned some quantitative and qualitative diﬀerences between weak and strong coupling can be ascertained. The case of weak coupling is considered ﬁrst. 4.1. Weak coupling The detuning parameter Δ was varied in the neighborhood of Δ = 0. For several values of the detuning, the critical stress σˆ is shown as a function of crit the excitation frequency in Fig. 5. In all four cases, the critical stress response exhibits two local maxima and an intermediate minimum. This is similar to the results depicted in Fig. 3a and c for the case Δ = 0. The location of the local minimum and the two local maxima depends on Δ. As Δ is increased, the lo- cation of the minimum is shifted to larger frequencies. The domain of inﬂuence of this phenomenon is limited to a comparatively small frequency window. For larger |Δ|, as in Figs. 5a and d, the dominant peak occurs around Ω = 1.03ω , and approaches the level of the detuned case. The overall eﬀect of the nonlinear modal interactions is most pronounced when the local minimum is located at Ω ≈ 1.03ω , i. e., for Δ ≈ 0.03. Not only the location of the secondary peaks varies, but also their height de- pends on Δ. For smaller Δ, the right secondary peak is larger than the left one. For a certain value of Δ = 0.017ω both peaks are equal. For larger Δ, the left peak is larger than the right one. This can be seen well in Fig. 6, which 15 Figure 5: Steady-state response of the normalized critical stress for diﬀerent detuning values: (a) Δ = −0.09ω , (b) Δ = 0.017ω =: Δ , (c) Δ = 0.05ω =: Δ , (d) Δ = 0.16ω 1 1 best 1 worst 1 Figure 6: Dependence of the maximum critical stress on detuning: (a) Global maximum of the steady-state critical stress, composed of the two local maxima as indicated in Fig. 5, (b) inﬂuence of the type of the linear viscous damping 16 shows the relationship σˆ (Δ) between the maximum stress level σˆ crit,max crit,max within the considered frequency range and the detuning parameter Δ. For a certain Δ, the stress level is minimized. The associated Δ = Δ is referred to best as optimum tuning. This corresponds to the point where the secondary peaks have equal stress level, see Fig. 5b. Since the critical frequency point switches at this point from the right to the left critical frequency, the characteristic has a kink at Δ = Δ . The height of the right maximum shrinks monotonically best for increasing Δ. In contrast, the left maximum grows quickly with Δ up to a certain point Δ , and then shrinks slightly again. Remarkably the stress worst level at Δ is larger than in the detuned case. In other words, the presence worst of an internal resonance has detrimental eﬀects on the friction damping perfor- mance. An explanation for this possibly unexpected phenomenon is given in Subsec. 4.3. The nonlinear modal interactions rely on the vibration energy transfer from low to high frequencies. The energy dissipation provided by material damping, which is assumed to be of viscous linear type in this work, depends on the vi- bration frequency: The higher the vibration frequency, the larger is the amount of energy dissipated within a given time span. That is why viscous damping is more eﬃcient for higher vibration frequencies. While the fundamental fre- quency component is largely conﬁned in the left beam’s ﬁrst mode, the higher frequency component is largely conﬁned in the right beam’s ﬁrst mode. The damping ratios associated with these two modes can have a signiﬁcant inﬂuence on the overall dynamic behavior of the system. A modal damping ratio D = 1% for all modes was assumed, as deﬁned in Eq. (7). In order to demonstrate that the qualitative behavior shown in Fig. 6a is invariant under the variation of the assumed damping model, two diﬀerent cases are considered: Mass-proportional damping with C = α M, and stiﬀness-proportional damping with C = β K. d d Generally, mass-proportional damping results in lower damping of higher fre- quency components, while the opposite holds for stiﬀness-proportional damp- ing. The coeﬃcients α , β were speciﬁed to maintain equal modal damping d d ratio of the ﬁrst mode D = 1% regardless of the damping model. Hence, 17 the damping ratio D of the interacting higher-frequency mode depends on the damping model: For Δ = 0, one has D = 1/3 for mass-proportional damping, D = D = 1% for equal modal modal damping, and D = 3% for stiﬀness- 2 1 2 proportional damping. The results are depicted in Fig. 6b. It can be clearly seen that the qualitative dependence on Δ does not change. However, quanti- tative diﬀerences exist that are caused by two diﬀerent eﬀects: (a) The better the higher frequency component is damped, the lower is the worst stress level (that is the largest value of σˆ ). (b) A stronger damping of the higher crit,max frequency component also reduces the response of the higher frequency mode, and may even suppress the nonlinear modal interactions. As a consequence, the stress-detuning curve ﬂattens out for larger damping of the higher frequency components. 4.2. Strong coupling 3EI We now consider the case of a much larger joint stiﬀness k strong = 125· = c 3 100k weak, which is 100 times larger than the coupling stiﬀness considered in the previous sections. A larger joint stiﬀness results in a stronger coupling between the two beams, and, hence, a more pronounced inﬂuence of the nonlinearity. As in the case of k weak, two quasi-linear limit cases exist for small and for large exci- tation levels, respectively, as can be ascertained from Fig. 3. The shift between the natural frequencies for sticking and free contact conditions is much larger for strong k , and the dynamics in the intermediate regime between the quasi-linear limit cases becomes much more complicated. The ﬁrst 21 Fourier terms were required in the harmonic balance method, compared to the ﬁve terms in the previous case, in order to achieve satisfactory agreement with direct numerical time integration. We would like to emphasize that the ‘noisiness’ in the inter- mediate regime is not due to the harmonic truncation but was also predicted by the numerical time integration. Instead, this behavior is considered intrinsic to the system. Due to the strong coupling, the non-smooth character of the stick-slip transitions is no longer ﬁltered out by the structure, but propagates to the global vibration behavior. For moderate and high excitation levels, the 18 Figure 7: Steady-state response of the beams’ tip displacement amplitudes for variation of the excitation level and strong coupling, Δ = 0: (a) Left beam, (b) right beam system again exhibits the most pronounced response in the neighborhood of Ω ≈ ω . Here, the dynamic behavior is comparatively regular and dominated by the ﬁrst component mode of the left and the right beam, associated with the modal coordinates q and q , respectively. The vibration energy is mainly con- 1 2 ﬁned to the ﬁrst and third harmonic, as in the previous case, while the response of the higher harmonics (not illustrated) is comparatively small. In the regime of higher excitation levels, it is thus still reasonable to investigate the nonlinear interaction between the ﬁrst two modes. We now consider the third-largest excitation level depicted in Fig. 7. The detun- ing parameter Δ was varied, and the steady-state forced response was computed in the frequency range 0.5ω ≤ Ω ≤ 2.8ω . The resulting maximum stress level 1 1 σˆ is depicted in Fig. 8. The qualitative dependence of σˆ on Δ dif- crit,max crit,max fers considerably from the results obtained for the case of weak coupling shown in Fig. 6. There is a qualitative diﬀerence, since the nonlinear modal interac- tions are only detrimental regarding the maximum stress level attained in the system. Hence, in the case of strong coupling, the harmonically forced system cannot be tuned to perform better due to the 1 : 3 internal resonance. More- over, there are also quantitative changes, as the nonlinear modal interactions cause a variation of the stress level σˆ of ≈ ±15 − 20% for k weak, whereas crit,max 19 Figure 8: Dependence of the maximum critical stress on detuning for strong coupling the stress level is more than doubled for k strong. The quantitative discrepancy is obviously caused by the increased relevance of the nonlinearity on the overall dynamics, owing to the larger joint stiﬀness k . The qualitative diﬀerence is less obvious. This is addressed in the following subsection. 4.3. On the detrimental eﬀects of the nonlinear modal interactions As was ascertained both for the cases of weak and strong coupling, the non- linear modal interactions caused by the 1 : 3 internal resonance between the two coupled beams can have detrimental eﬀects on the friction damping per- formance. This, probably, is an unexpected outcome. Apparently, the presence of the higher frequency component (mainly in the right beam) in the vibration response reduces the eﬀective damping of the fundamental frequency compo- nent (mainly concerning the left beam). Note that damping of the fundamental frequency component is most important, since the critical response level was always reached in the left beam, whose steady-state response is dominated by the fundamental frequency component. In what follows, the eﬀective damping provided by the elastic Coulomb friction element is investigated in more detail, in order to conﬁrm the aforementioned assertions. 20 First, we deﬁne a measure for the eﬀective damping of the fundamental fre- quency component. To this end, consider a periodic, two-harmonic joint dis- placement w (t), iΩt iδ 3iΩt w (t) = ℜ{a e + a e e } . (16) c 1 3 Herein, a and a are the amplitudes of the ﬁrst and third harmonic, and δ 1 3 31 is the possible phase lag between the harmonic components. The following eﬀective damping coeﬃcient is deﬁned, 1 −iΩt (iΩa ) f [w (Ωt)] e dΩt 1 c c 2π (2π) d = . (17) eﬀ,1 Ω a Herein, denotes the complex conjugate, and f [w (Ωt)] is the nonlinear force c c characteristic of the elastic Coulomb element. It can be veriﬁed that if f was replaced by a linear viscous damper force with coeﬃcient d, Eq. (17) would yield the expected value d = d. eﬀ,1 The elastic Coulomb model has two parameters, namely µN and k . The model possesses scaling properties with respect to these parameters as follows, f [w ]| = γ f w , (18) c c c c γN µ N µ f [w ]| = f [γw ]| , (19) c c c c γk k c c for arbitrary γ ∈ R . In order to obtain a complete picture of the behavior of d , it is thus suﬃcient to keep k and µN ﬁxed and to determine the de- eﬀ,1 c pendence of this measure only on a , a and δ . From the resulting functional 1 3 31 relationship, d can then be reproduced for any given value for k and µN by eﬀ,1 c utilizing Eqs. (18)-(19). The dependence of the eﬀective damping measure d on the parameters a , eﬀ,1 1 and δ was studied in detail. The investigation was limited to those pa- rameter ranges that are relevant for the results reported in this work. To this end, the forced response depicted in Fig. 6 and Fig. 8 was considered and the corresponding values for a , and δ were determined. The results of this 1 31 analysis are not presented here for the sake of brevity, and the measure d eﬀ,1 was computed within the relevant parameter ranges using Eq. (17). The results 21 Figure 9: Eﬀective damping of the elastic Coulomb element provided to the ﬁrst harmonic, versus its amplitude a : (a) variation of the phase lag between the harmonic components 0 ≤ δ ≤ π, for ﬁxed = 0.075; (b) variation of the amplitude ratio between the harmonic components 0 ≤ ≤ 0.2, for ﬁxed δ = 0.85π; the left and right zooms correspond to the relevant amplitude ranges for weak and strong coupling, respectively are illustrated in Fig. 9a-b. The evolution versus a qualitatively resembles previously reported results ob- tained in the absence of nonlinear modal interactions, e. g., [14, 15, 16]. We therefore limit the discussion to the inﬂuence of the higher harmonic compo- nent. The left and right zooms in each sub-ﬁgure in Fig. 9 correspond to the relevant ranges of the amplitude a for weak and strong coupling, respectively. Apparently, diﬀerent properties of the vibration response cause the detrimental reduction of d in the cases of weak and strong coupling. eﬀ,1 • In the case of weak coupling, the phase lag δ of the third harmonic plays the central role, almost independent of its magnitude a . Under the condition of in-phase (anti-phase) joint motion, in a multi-harmonic sense, the elastic Coulomb element is most (least) eﬀective. • In contrast, the third harmonic amplitude of the input joint displacement, a , is crucial in the case of strong coupling, regardless of its phase. The larger a , the smaller the eﬀective damping of the elastic Coulomb element. It is well-known that the amplitude ratio commonly undergoes considerable 22 variation in the forced response in the neighborhood of an internal resonance. Recent investigations by the authors suggest that the same is true for the so- called motion complexity, i. e., the phase diﬀerence between the coordinates of the nonlinearly interacting modes [17, 18]. 5. Conclusions The concept of friction damping relies on the inherently nonlinear eﬀects of dry friction. The nonlinear forces can trigger nonlinear interactions among (nearly) internally resonant modes, in spite of the considerable damping that is present in the system. The inﬂuence on the steady-state forced vibrations is particularly pronounced in the neighborhood of the optimum design point, and therefore appears relevant for design of friction-damped systems. Depending on the system properties, it is indeed possible to take advantage of the energy re- distribution caused by the modal interactions, and to optimally tune the system to achieve a minimum response level. In the numerical study, a reduction of the maximum response level on the order of 10% was achieved, which is deemed rel- atively small. More importantly, the vibration behavior is much more sensitive, and thus susceptible to inherent uncertainties, as compared to the detuned case in the absence of modal interactions. Moreover, parameter regimes exist, for which the nonlinear modal interactions are detrimental. Apparently, the pres- ence of higher frequency components in the relative joint displacement can di- minish the eﬀective damping provided by the friction joint. In the last example, the response level was almost doubled. It is thus concluded that the condition of internal resonance should be avoided in the design of friction-damped systems. References [1] A. A. Ferri, Friction damping and isolation systems, Journal of Vibration and Acoustics 117 (B) (1995) 196–206. [2] K. Popp, L. Panning, W. Sextro, Vibration damping by friction forces: Theory and applications, Journal of Vibration and Control 9 (3-4) (2003) 419–448. 23 [3] A. H. Nayfeh, Nonlinear Interactions: Analytical, Computational and Ex- perimental Methods, John Wiley & Sons, New York, 2000. [4] W. E. Whiteman, A. A. Ferri, Multi-mode analysis of beam-like structures subjected to displacement-dependent dry friction damping, Journal of Sound and Vibration 207 (3) (1997) 403–418. [5] M. Berthillier, C. Dupont, R. Mondal, J. J. Barrau, Blades forced response analysis with friction dampers, Journal of Vi- bration and Acoustics 120 (2) (1998) 468–474. [6] J. J. Chen, B. D. Yang, C. H. Menq, Periodic forced response of structures having three-dimensional frictional constraints, Journal of Sound and Vibration 229 (4) (2000) 775–792. [7] J. J. Chen, C. H. Menq, Periodic response of blades having three-dimensional nonlinear shroud constraints, Journal of Engineering for Gas Turbines and Power 123 (4) (2001) 901–909. [8] M. Krack, L. Panning-von Scheidt, J. Wallaschek, A method for nonlin- ear modal analysis and synthesis: Application to harmonically forced and self-excited mechanical systems, Journal of Sound and Vibration 332 (25) (2013) 6798–6814. [9] A. A. Ferri, B. S. HECK, Analytical investigation of damping enhancement using active and passive structural joints: Journal of guidance, control, and dynamics, Journal of Guidance, Control, and Dynamics 15 (5) (1992) 1258– [10] N. Do, A. A. Ferri, Energy transfer and dissipation in a three-degree-of- freedom system with stribeck friction, in: ASME 2005 International Me- chanical Engineering Congress and Exposition, 2005, pp. 195–204. [11] N. Do, A. A. Ferri, Properties of energy transfer within a dry fric- tion damped structural system, in: 49th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference, 16th 24 AIAA/ASME/AHS Adaptive Structures Conference,10th AIAA Non- Deterministic Approaches Conference, 9th AIAA Gossamer Spacecraft Fo- rum, 4th AIAA Multidisciplinary Design Optimization Specialists Con- ference, Structures, Structural Dynamics, and Materials and Co-located Conferences, American Institute of Aeronautics and Astronautics, 2008. [12] A. F. Vakakis, O. V. Gendelman, G. Kerschen, L. A. Bergman, D. M. McFarland, Y. S. Lee, Nonlinear Targeted Energy Transfer in Mechanical and Structural Systems, Springer Berlin Heidelberg, 2008. [13] M. Krack, S. Tatzko, L. Panning-von Scheidt, J. Wallaschek, Reliability optimization of friction-damped systems using nonlinear modes, Journal of Sound and Vibration 333 (2014) 2699–2712. [14] W. E. Whiteman, A. A. Ferri, Displacement-dependent dry friction damping of a beam-like structure, Journal of Sound and Vibration 198 (3) (1996) 313–329. [15] D. Laxalde, F. Thouverez, Complex non-linear modal analysis for mechanical systems application to turbomachinery bladings with friction interfaces, Journal of Sound and Vibration 322 (4-5) (2009) 1009–1025. [16] M. Krack, Nonlinear modal analysis of nonconservative systems: Extension of the periodic motion concept, Computers and Structures 154 (2015) 59– [17] M. Krack, L. A. Bergman, A. F. Vakakis, Motion complexity in a non- classically damped system with closely spaced modes: From standing to traveling waves, in press, Proceedings of the Institution of Mechanical En- gineers, Part K: Journal of Multibody Dynamics (2015) 1–13. [18] M. Krack, L. A. Bergman, A. F. Vakakis, Global complexity eﬀects due to local damping in a nonlinear system in 1:3 internal resonance, submitted to the Archive of Applied Mechanics.
Nonlinear Sciences – arXiv (Cornell University)
Published: Dec 28, 2020
You can share this free article with as many people as you like with the url below! We hope you enjoy this feature!
Read and print from thousands of top scholarly journals.
Already have an account? Log in
Bookmark this article. You can see your Bookmarks on your DeepDyve Library.
To save an article, log in first, or sign up for a DeepDyve account if you don’t already have one.
Copy and paste the desired citation format or use the link below to download a file formatted for EndNote
Access the full text.
Sign up today, get DeepDyve free for 14 days.
All DeepDyve websites use cookies to improve your online experience. They were placed on your computer when you launched this website. You can change your cookie settings through your browser.