Simulation of Transient Nonlinear Friction-Induced Vibrations Using Complex Interface Modes: Application to the Prediction of Squeal Events

Sinou, J.-J.; Besset, S.

doi:https://doi.org/10.1155/2017/6107251

Shock and Vibration

On this page

Abstract Introduction Conclusion Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2017 | Article ID 6107251 | https://doi.org/10.1155/2017/6107251

Simulation of Transient Nonlinear Friction-Induced Vibrations Using Complex Interface Modes: Application to the Prediction of Squeal Events

J.-J. Sinou^1,2,3and S. Besset^1,2

Academic Editor: Tony Murmu

Received25 Sept 2016

Accepted30 Nov 2016

Published15 Jan 2017

Abstract

During the past decades, the problem of friction-induced vibration and noise has been the subject of a huge amount of works. Various numerical simulations with finite elements models have been largely investigated to predict squeal events. Although a nonlinear analysis is more predictive than Complex Eigenvalues Analysis, one of the main drawbacks of the time analysis is the need of large computational efforts. In view of the complexity of the subject, this approach appears still computationally too expensive to be used in industry for finite element models. In this study, the potential of a new reduced model based on a double modal synthesis (i.e., a classical modal reduction via Craig and Bampton plus a condensation at the frictional interface based on complex modes) for the prediction of self-excited vibrations of brake squeal is discussed. The effectiveness of the proposed modal reduction is tested on a finite element model of a simplified brake system. It will be shown that numerical results of times analysis by applying the proposed reduction correlate well with those of the nonlinear analysis based on a reference model, hence demonstrating the potential of using adapted modal reductions to predict the squeal propensity and to estimate self-excited vibrations and noise.

1. Introduction

Brake squeal is the result of friction-generated vibration and noise. It corresponds to a very common problem in industrial application. Thereby friction-induced vibration has become one of the most costly problems and important concerns in automotive engineering. Indeed this instability which results in a harsh noise, especially in the frequency range between 0.1 kHz and 15 kHz, is annoying to passengers’ hearing which leads to consistent customer complaints. Comprehensive overviews on friction-induced vibration and much classic understanding of the mechanism of the disc brake squeal phenomenon can be found in [1–7].

Finite element models and numerical simulations have been extensively investigated to detect brake squeal instabilities and to predict amplitudes during squeal events. Nowadays three kinds of analysis can be performed in the design and optimization process of brake systems in order to avoid or at least to reduce the squeal occurrence: eigenvalue analysis to detect frequency instabilities, time analysis, and nonlinear methods [8] to determine self-excited vibrations during the squeal event and acoustic analysis to estimate acoustic radiation. While the Complex Eigenvalue Analysis (CEA) represents a fundamental tool that is widely used in industry to predict unstable vibration modes, this approach may lead to under- or overestimation of unstable modes [9–14], and more realistic results can be achieved via the time nonlinear analysis [15]. While a brake squeal is characterized primarily as an acoustic phenomenon, the prediction of the acoustic response during squeal events has been usually ignored. There have been only few attempts to calculate and analyze the acoustic radiation of brake squeal [16–19]. Indeed the calculation of acoustic events necessarily involves an estimate of the nonlinear amplitudes and velocities of brake squeal via a time analysis. One of the major drawbacks of the nonlinear time simulation, nevertheless, is its high computational cost which explains that this approach based on time analysis has mostly been neglected in the past. As a consequence, it is difficult to calculate and track the acoustic radiation during squeal events. So the development of numerical methods to facilitate the prediction of amplitude vibrations during squeal events is one of the crucial issues and challenges.

As previously explained working in the time domain and performing simulations on a full finite element model require prohibitive computational costs. In order to overcome these issues, one solution is to perform model reduction. The use of appropriate reduction bases for the prediction of squeal frequencies and self-excited vibrations during the squeal event and its associated numerical implementation are complex tasks. Recently some efforts have been spent to overcome this problem by considering the potential of using reduction techniques of finite element models. Vermot des Roches [20] proposed a reduction method adapted to large models with local nonlinearities. The reduction method allows keeping the exact system real modes. Even if the computations remain costly, the numerical results were very encouraging. These results contributed greatly to new developments on the use of different modal reductions for the community of brake squeal. Then Loyer et al. [21] proposed using spatial model reductions for friction destabilized systems in order to study the nonlinear behaviors of an elastic layer. They concluded that the inclusion of static modes in the reduction basis is mandatory to approximate the self-excited vibrations of the mechanical system. They also indicated that the model reduction may lead to an alteration of the transient evolution, often leading to an advance phenomenon.

Recently, Besset and Sinou [22] proposed a new efficient reduction method based on a generalized double modal synthesis method. This numerical technique can significantly reduce the size of the original finite element model of the brake system by combining a classical modal reduction and a condensation at the frictional interface via a reduced complex mode basis. One of the main limitations of this previous study is the restriction of the problem on the prediction of frequency instabilities. One objective of this paper is, therefore, to address the use of the generalized double modal synthesis method in conducting rigorous numerical simulations of brake squeal for both the detection of squeal frequencies and the calculation of self-excited vibrations during squeal events. This paper consists of four major sections: introduction, finite element brake model and reduction method adapted to numerical brake system, numerical simulations and results based on the Complex Eigenvalue Analysis and the transient analysis, and conclusions.

2. Numerical Brake Model

In this section, the finite element model of the brake system is first presented. Then the double modal reduction based on the Craig and Bampton reduction and the interface reduction via complex eigenmodes is briefly introduced.

2.1. Finite Element Model

A simplified model of a brake system (i.e., “pad-on-disc”: a steel pad in sliding contact with an annular disc) is depicted in Figure 1. The finite element model used in this paper was explained previously in detail in [22]. At the upper surface, the pad is constrained in the tangential and radial directions (with no rotational degrees of freedom). The inner surface of the disc is clamped. The finite frictional surface-to-surface contact is composed of 220 contact nodes. It can be noted that a convergence study versus the interface refinement (i.e., the number of contact points at the frictional interface) has been carried on in order to ensure the independence of the simulation results versus the interface meshing. Details on the geometric and material properties of the finite element model are given in Table 1.

(a)

(b)

The nonlinear equation of the numerical model can be written in the following form:where , , and are the classical mass, stiffness, and damping matrices, respectively. A modal damping (with the damping coefficients fixed at 0.1% for each mode) is considered. is the generalized displacement vector while the dot denotes derivative with respect to time. represents the hydraulic pressure that is directly applied to the back-plate of the pad. The pressure force applied on the upper surface of the pad is equal to Pa in our case. A velocity field corresponding to the disc rotation is imposed on the disc surface. contains contributions from both the contact nonlinear forces and frictional forces at the pad/disc interface with a classical Coulombs law. Considering previous experimental investigations [23], we assume that the nonlinearities at the friction interface between the pad and the disc are due to the cubic nonlinear terms and the possible loss of contact. The components of the nonlinear contact force vectors denoted by and ( and , resp.) along the normal direction (along the tangential direction, resp.) are given at contact node bywhere corresponds to the relative displacement between the pad and the disc for the contact node. and define the linear and nonlinear stiffness at the friction interface. Values of these linear and nonlinear stiffness are fixed at N·m⁻¹ and N·m⁻³, respectively. is the relative velocity between the pad and the disc and defines the friction coefficient that is assumed to be constant.

2.2. Reduction Method Adapted to Numerical Brake Model

In order to reduce the finite element model and more specifically to condense the pad and disc, the Craig and Bampton technique is first used. This reduction technique is one of the most classical processes to reduce the finite element by retaining all the contact nodes at the pad/disc interface and generalized degrees of freedom. The Craig and Bampton technique [24] consists in building a reduced model using fixed interface modes and constraint modes. This general process with the application on the brake model under study was explained previously in detail in [22]. Convergence of Craig and Bampton method for a similar finite element brake system has also been discussed in [25].

Secondly, a reduction method using complex interface modes is applied in order to condense the pad/disc section (i.e., the remaining physical interface degrees of freedom that are kept unreduced via the Craig and Bampton method). This reduction method is now briefly explained.

In the process, the mass and stiffness matrices (denoted by and , resp.) are firstly implicitly reduced. They are obtained by using the transfer matrix :where and define the interface degrees of freedom and the generalized degrees of freedom, respectively (see [22] for more details).

By considering the reduced mass and stiffness matrices via the Craig and Bampton method (denoted by and , resp.) and the transfer matrix , we obtain

For the reader comprehension, the reduced stiffness matrix via the Craig and Bampton method contains not only the structural stiffness matrix but also the Jacobian matrix of the linearized expression of the nonlinear contact and friction forces at the disc/pad interfaces around the nonlinear equilibrium (see [22] for more details).

Then, right and left eigenvectors (denoted by and ) are calculated by resolving the following equations:

By introducing the matrices and of the first modes and , respectively, and the two matrices and ,

Finally we obtain the reduced brake model:where defines the generalized degrees of freedom of the reduced brake system. We have

Algorithm in universal formal pseudocode is provided in [22].

3. Numerical Simulations

The numerical simulations and results are decomposed into two main parts: the stability analysis via CEA and the transient nonlinear analysis.

The stability analysis that is the first step in order to predict the squeal propensity of a numerical brake model is only briefly presented in this study. This is explained by the fact that a comprehensive study was conducted previously by the authors on this subject in [22].

In the second part of this section, numerical simulations of the transient nonlinear self-excited vibrations will be presented. More specifically, efficiency of the reduction technique based on complex modes will be illustrated. All the results will be analyzed and compared to a reference calculation based on the finite element model without any interface condensation.

3.1. Stability Analysis

The stability analysis of the numerical brake model can be performed by using the CEA on the reduced finite element model (see (7)). As a reminder, the nonlinear steady-state operating point and the Jacobian matrix of the linearized expression of the nonlinear contact and friction forces at the disc/pad interfaces are taken into account in this calculation.

The stability of a brake system is determined by the sign of real part of eigenvalues: an equilibrium is asymptotically stable if all eigenvalues have negative real parts; it is unstable if at least one eigenvalue has positive real part. In this last case instability is generated from the equilibrium point and self-excited vibration occurs.

Figures 2(a) and 2(b) display the evolutions of real parts and frequencies as a function of the friction coefficient, respectively. As previously demonstrated in [22], the coalescence pattern is well predicted by the reduction process based on complex interfaces modes (in comparison with the reference numerical model). The curves highlight the phenomenon referred to as mode coupling: as the friction coefficient increases, the real parts increase or decrease slowly and the Hopf bifurcation point occurs for a friction value of . The associated frequencies of the two modes involved in the coalescence (i.e., the stable and unstable modes) tend to get closer. It is observed that the two modes reach not exactly the Hopf bifurcation point (for ) due to the non-equally damped configuration between the two modes (see [26, 27] for more details on the effect of damping for the coalescence patterns). After the Hopf bifurcation point, the frequencies of the two modes remain close and their associated real parts become opposite: one mode is stable (i.e., the associated real part is negative and decreases versus the increase of the friction coefficient) whereas the other one is unstable (i.e., the associated real part is positive and increases versus the increase of the friction coefficient).

(a)

(b)

3.2. Transient Simulation and Self-Excited Vibrations

Based on the previous CEA results (see Figure 2 and Section 3.1), the following friction coefficients were chosen for the transient nonlinear simulation and analysis:(i) with no instability.(ii) and with one instability.

Firstly the time series (transient dynamics and self-excited vibrations), then the spectral content by using short time Fourier transform (to visualize faster time-varying signals), and finally the limit cycles and the Fast Fourier Transform (to represent the nonlinear stationary signal in the frequency domain) are studied in detail. The displacements selected to illustrate the results correspond to(i)one degree-of-freedom of the disc in the direction normal to the disc surface (i.e., the direction normal to the frictional interface),(ii)one degree-of-freedom of the disc in the direction tangential to the disc surface.

For each case, results from the reduction process based on complex interfaces modes are compared with those from the reference model. Convergence versus the number of interface modes has been conducted (results not presented in the following in order not to overload the paper content). In the following, it has been chosen to show results with 15 interface modes.

Figures 3(a), 3(b), 4(a), and 4(b) depict the time series and the spectrograms for for one degree-of-freedom of the disc in the direction normal and tangential to the contact surface, respectively. It appears that the initial vibration amplitude decreases until no more movement occurs. The same behavior is obtained for the reduction process based on complex interfaces modes and the reference numerical model. It is also noted that the time simulation agrees with the CEA results that predict no instability for .

(a)

(b)

(c)

(d)

(e)

(f)

(a)

(b)

(c)

(d)

(e)

(f)

As the friction coefficient is increased to , the vibration amplitudes grow (see Figures 3(c), 3(d), 4(c), and 4(d)). It can be seen that a fast increase of the transient oscillations (i.e., the initial growing transient oscillations with the highest amplitudes) is followed by a slight decrease until the nonlinear stationary periodic motions. Evolution of the transient vibrations and the associated spectrograms during the transient oscillations and stationary self-excited vibrations for the reference model and the reduced system are very similar (compare Figures 3(c) and 3(d) and Figures 4(c) and 4(d), resp.). Moreover, Figures 5(a), 5(b), 6(a), and 6(b) give the limit cycles and Fast Fourier Transform on the stationary self-excited vibrations for both the normal and tangential directions. It is clearly shown that the results based on the reduction technique are in perfect agreement with the reference performed via the finite element model without any interface condensation. As indicated in Figures 5(b) and 6(b), the FFTs exhibit two significant peaks: the dominant peak is related to the contribution of the fundamental frequency = 394 Hz. It is noted that this frequency for the stationary self-excited oscillations is very close (but not equal) to the frequency of the unstable mode predicted by the CEA at , as indicated in Table 2. Also the frequency of the second predominant peak (at 788 Hz) is equal to a multiple (i.e., a superharmonic) of the fundamental frequency and corresponds to the second harmonic component (i.e., ). The contribution of 2x harmonic component is also visible on the limit cycles with the appearance of a “small-inside loop.” Moreover it is also observed that the FFTs exhibit the superharmonics (with = 3, 4, and ) with smaller contributions to the global nonlinear vibrations. Thus it appears that new frequency peaks that are not predicted by the stability analysis are visible in the transient and stationary signals during squeal events. As previously explained in [9, 10], CEA may lead to underestimation of resonances and more complex and realistic results can be achieved via the nonlinear analysis of the numerical brake model in order to predict the contribution of harmonic components. It can be noted that the fundamental frequency and the resonances of 2x and 3x harmonic components are also highlighted during the transient simulations (see Figures 3(c), 3(d), 4(c), and 4(d)).

(a)

(b)

(c)

(d)

(a)

(b)

(c)

(d)

Increasing the friction coefficient further to does not show major differences, as illustrated in Figures 3(e), 3(f), 4(e), and 4(f). The squeal event corresponds to nonlinear harmonic self-excited friction-induced vibrations of the mode-coupling mechanism with the participation of one single instability: we have only one fundamental frequency (at = 367 Hz) and its superharmonics participate in the nonlinear oscillations. However, it appears that the nonlinear behavior is more complex with the major contributions of not only the 1x and 2x harmonic components (at about 367 Hz and 734 Hz) but also a strong participation of the 3x harmonic component (at about 1101 Hz) during the transient time simulation of squeal events (see Figures 5(d) and 6(d)). Considering more specifically the nonlinear steady-state friction-induced motions, the FFTs exhibit other superharmonics . However, all the combinations with are generally less important than the fundamental frequencies as well as the second and third harmonic components (i.e., and ). Considering more specifically the FFT for in the tangential direction (see Figure 6(d)), two additional comments can be made: firstly, the major contribution comes from the 3x harmonic component. Secondly, the contributions of the 5x and 2x harmonic components are greater than that of the fundamental frequency. It can also be noted that the fundamental frequency of the stationary self-excited oscillations is slightly different from the frequency of the unstable mode predicted by the CEA at , as indicated in Table 2. These results highlight the need to consider time analysis not only to estimate the contribution of each harmonic components of the self-excited vibrations but also to predict the value of the fundamental frequency and the main frequency contribution during squeal events. It can be noted that the values of the fundamental frequency calculated from the reference model and the reduced system for the transient time analysis are very close (as shown in Table 2 for and ).

Finally, all the results via the reduction technique based on complex modes (for the time series, the spectral content, and the limit cycles with the Fast Fourier Transform) are in perfect agreement with the reference performed via the finite element model without any interface condensation. This illustrates the efficiency and robustness of the proposed reduction method. The computation time and the number of degrees of freedom of the numerical model are extensively reduced by using the reduction process based on complex interfaces modes. This demonstrates the potential of the proposed method (see Table 3 that highlights the less computational cost of the proposed method).

4. Conclusion

Extension of the double modal synthesis approach [22] has been proposed and used to study the self-excited vibrations and time analysis during squeal events. By conducting a direct comparison with a reference model, the proposed modal reduction strategy has been found to be highly effective not only to predict the transient and stationary amplitudes but also to estimate the nonlinear vibration signature during squeal events.

The proposed methodology has the potential to drastically improve the computational costs. It is not dependent on the particular choices of modeling (geometry of the brake system, friction model, or local nonlinearities in the system) and it appears easily implementable on a finite element model.

In addition, numerical results highlight that a nonlinear analysis is required to predict accuracy of the squeal events. The applicability of the proposed modal reduction strategy to a full industrial brake system will be the next step in a near future. Then appropriate comparisons and correlations with experiments could be investigated.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

J.-J. Sinou acknowledges the support of the Institut Universitaire de France.

References

J. T. Oden and J. A. C. Martins, “Models and computational methods for dynamic friction phenomena,” Computer Methods in Applied Mechanics and Engineering, vol. 52, no. 1–3, pp. 527–634, 1985.
View at: Publisher Site | Google Scholar | MathSciNet
N. M. Kinkaid, O. M. O'Reilly, and P. Papadopoulos, “Automotive disc brake squeal,” Journal of Sound and Vibration, vol. 267, no. 1, pp. 105–166, 2003.
View at: Publisher Site | Google Scholar
R. A. Ibrahim, “Friction-induced vibration, chatter, squeal and chaos. Part II: dynamics and modelling,” Applied Mechanics Review, vol. 47, pp. 227–253, 1994.
View at: Publisher Site | Google Scholar
R. A. Ibrahim, “Friction-induced vibration, chatter, squeal, and chaos—part I: mechanics of contact and friction,” Applied Mechanics Review, vol. 47, no. 7, pp. 209–226, 1994.
View at: Publisher Site | Google Scholar
D. A. Crolla and A. M. Lang, “Brake noise and vibration—state of art,” Vehicule Tribology, Tribology Series, vol. 18, pp. 165–174, 1991.
View at: Google Scholar
H. Ouyang, W. Nack, Y. Yuan, and F. Chen, “Numerical analysis of automotive disc brake squeal: a review,” International Journal of Vehicle Noise and Vibration, vol. 1, no. 3-4, pp. 207–231, 2005.
View at: Publisher Site | Google Scholar
A. Papinniemi, J. C. S. Lai, J. Zhao, and L. Loader, “Brake squeal: a literature review,” Applied Acoustics, vol. 63, no. 4, pp. 391–400, 2002.
View at: Publisher Site | Google Scholar
N. Coudeyras, S. Nacivet, and J.-J. Sinou, “Periodic and quasi-periodic solutions for multi-instabilities involved in brake squeal,” Journal of Sound and Vibration, vol. 328, no. 4-5, pp. 520–540, 2009.
View at: Publisher Site | Google Scholar
J.-J. Sinou, “Transient non-linear dynamic analysis of automotive disc brake squeal—on the need to consider both stability and non-linear analysis,” Mechanics Research Communications, vol. 37, no. 1, pp. 96–105, 2010.
View at: Publisher Site | Google Scholar
J.-J. Sinou, A. Loyer, O. Chiello et al., “A global strategy based on experiments and simulations for squeal prediction on industrial railway brakes,” Journal of Sound and Vibration, vol. 332, no. 20, pp. 5068–5085, 2013.
View at: Publisher Site | Google Scholar
F. Massi, L. Baillet, O. Giannini, and A. Sestieri, “Brake squeal: Linear and nonlinear numerical approaches,” Mechanical Systems and Signal Processing, vol. 21, no. 6, pp. 2374–2393, 2007.
View at: Publisher Site | Google Scholar
S. Oberst and J. C. S. Lai, “Statistical analysis of brake squeal noise,” Journal of Sound and Vibration, vol. 330, no. 12, pp. 2978–2994, 2011.
View at: Publisher Site | Google Scholar
S. Oberst and J. C. S. Lai, “Chaos in brake squeal noise,” Journal of Sound and Vibration, vol. 330, no. 5, pp. 955–975, 2011.
View at: Publisher Site | Google Scholar
S. Oberst and J. C. S. Lai, “Nonlinear transient and chaotic interactions in disc brake squeal,” Journal of Sound and Vibration, vol. 342, pp. 272–289, 2015.
View at: Publisher Site | Google Scholar
A. R. AbuBakar and H. Ouyang, “Complex eigenvalue analysis and dynamic transient analysis in predicting disc brake squeal,” International Journal of Vehicle Noise and Vibration, vol. 2, no. 2, pp. 143–155, 2006.
View at: Publisher Site | Google Scholar
S. Oberst, J. C. S. Lai, and S. Marburg, “Guidelines for numerical vibration and acoustic analysis of disc brake squeal using simple models of brake systems,” Journal of Sound and Vibration, vol. 332, no. 9, pp. 2284–2299, 2013.
View at: Publisher Site | Google Scholar
K. Soobbarayen, S. Besset, and J.-J. Sinou, “Noise and vibration for a self-excited mechanical system with friction,” Applied Acoustics, vol. 74, no. 10, pp. 1191–1204, 2013.
View at: Publisher Site | Google Scholar
K. Soobbarayen, J.-J. Sinou, and S. Besset, “Numerical study of friction-induced instability and acoustic radiation—effect of ramp loading on the squeal propensity for a simplified brake model,” Journal of Sound and Vibration, vol. 333, no. 21, pp. 5475–5493, 2014.
View at: Publisher Site | Google Scholar
K. Soobbarayen, S. Besset, and J.-J. Sinou, “A simplified approach for the calculation of acoustic emission in the case of friction-induced noise and vibration,” Mechanical Systems and Signal Processing, vol. 50, pp. 732–756, 2015.
View at: Google Scholar
G. Vermot des Roches, Frequency and time simulation of squeal instabilities. Application to the design of industrial automotive brakes [Ph.D. thesis], École Centrale Paris, Châtenay-Malabry, France, 2011.
A. Loyer, J.-J. Sinou, O. Chiello, and X. Lorang, “Study of nonlinear behaviors and modal reductions for friction destabilized systems. Application to an elastic layer,” Journal of Sound and Vibration, vol. 331, no. 5, pp. 1011–1041, 2012.
View at: Publisher Site | Google Scholar
S. Besset and J.-J. Sinou, “Modal reduction of brake squeal systems using complex interface modes,” Mechanical Systems and Signal Processing, vol. 85, pp. 896–911, 2017.
View at: Publisher Site | Google Scholar
J.-J. Sinou, N. Coudeyras, and S. Nacivet, “Study of the nonlinear stationary dynamic of single and multi-instabilities for disk brake squeal,” International Journal of Vehicle Design, vol. 51, no. 1-2, pp. 207–222, 2009.
View at: Publisher Site | Google Scholar
M. C. C. Bampton and R. R. Craig Jr., “Coupling of substructures for dynamic analyses,” AIAA Journal, vol. 6, no. 7, pp. 1313–1319, 1968.
View at: Publisher Site | Google Scholar
M. Monteil, S. Besset, and J.-J. Sinou, “A double modal synthesis approach for brake squeal prediction,” Mechanical Systems and Signal Processing, vol. 70-71, pp. 1073–1084, 2016.
View at: Publisher Site | Google Scholar
O. N. Kirillov, “Destabilization paradox due to breaking the Hamiltonian and reversible symmetry,” International Journal of Non-Linear Mechanics, vol. 42, no. 1, pp. 71–87, 2007.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
G. Fritz, J.-J. Sinou, J.-M. Duffal, and L. Jézéquel, “Effects of damping on brake squeal coalescence patterns—application on a finite element model,” Mechanics Research Communications, vol. 34, no. 2, pp. 181–190, 2007.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2017 J.-J. Sinou and S. Besset. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1036

Downloads

891

Citations