Vibroacoustic Simulation. Alexander Peiffer

Читать онлайн книгу.

Vibroacoustic Simulation - Alexander Peiffer


Скачать книгу
alt="StartLayout 1st Row 1st Column left-parenthesis Start 1 By 1 Matrix 1st Row normal upper Phi EndMatrix Superscript upper T Baseline Start 1 By 1 Matrix 1st Row upper K EndMatrix Start 1 By 1 Matrix 1st Row normal upper Phi EndMatrix minus omega squared Start 1 By 1 Matrix 1st Row normal upper Phi EndMatrix Superscript upper T Baseline Start 1 By 1 Matrix 1st Row upper M EndMatrix Start 1 By 1 Matrix 1st Row normal upper Phi EndMatrix right-parenthesis Start 1 By 1 Matrix 1st Row bold-italic q EndMatrix prime 2nd Column Blank 2nd Row 1st Column equals left-parenthesis Start 1 By 1 Matrix 1st Row upper K EndMatrix prime minus omega squared Start 1 By 1 Matrix 1st Row upper M EndMatrix prime right-parenthesis Start 1 By 1 Matrix 1st Row bold-italic q EndMatrix Superscript prime Baseline 2nd Column equals Start 1 By 1 Matrix 1st Row normal upper Phi EndMatrix Superscript upper T Baseline Start 1 By 1 Matrix 1st Row bold-italic upper F EndMatrix equals Start 1 By 1 Matrix 1st Row bold-italic upper F EndMatrix prime EndLayout"/> (1.124)

      Fn′ are the components of the modal force vector F′. Using the above defined orthogonality and normalisation we get:

       StartLayout 1st Row Start 4 By 4 Matrix 1st Row 1st Column omega 1 squared minus omega squared 2nd Column Blank 3rd Column Blank 4th Column Blank 2nd Row 1st Column Blank 2nd Column down-right-diagonal-ellipsis 3rd Column Blank 4th Column Blank 3rd Row 1st Column Blank 2nd Column Blank 3rd Column omega Subscript n Superscript 2 Baseline minus omega squared 4th Column Blank 4th Row 1st Column Blank 2nd Column Blank 3rd Column Blank 4th Column omega Subscript upper N Superscript 2 Baseline minus omega squared EndMatrix Start 4 By 1 Matrix 1st Row bold-italic q prime 1 2nd Row bold-italic q Subscript n Superscript prime Baseline 3rd Row vertical-ellipsis 4th Row bold-italic q Subscript upper N Superscript prime Baseline EndMatrix equals Start 4 By 1 Matrix 1st Row bold-italic upper F prime 1 2nd Row bold-italic upper F prime Subscript n 3rd Row vertical-ellipsis 4th Row bold-italic upper F prime Subscript upper N EndMatrix prime EndLayout (1.125)

      or in another form:

       bold-italic q prime Subscript n Baseline equals StartFraction bold-italic upper F prime Subscript n Baseline Over omega Subscript n Superscript 2 Baseline minus omega squared EndFraction (1.126)

      The response shape is reconstructed using (1.122). Thus, with modal coordinates expressed by independent degrees of freedom – the modes – the system decouples. This can by seen by the diagonal form of the matrix. For a system with N degrees of freedom, the modal solution is exact. The similarity to Equation (1.28) underlines that every mode can be interpreted as a single resonator. Thus, every dynamical system can be interpreted as a group of independent resonators of modal mass mn and stiffness ωn2mn In accordance with Equation (1.63) we may introduce a modal damping ηn for every mode. In that case the response in modal coordinates is:

       bold-italic q prime Subscript n Baseline equals StartFraction bold-italic upper F prime Subscript n Baseline Over omega Subscript n Superscript 2 Baseline left-parenthesis 1 plus j eta Subscript n Baseline right-parenthesis minus omega squared EndFraction (1.127)

      Finally, we conclude that the change in the coordinate base to modal coordinates keeps the general format of the equation of motion.

       left-parenthesis Start 1 By 1 Matrix 1st Row upper K EndMatrix prime minus omega squared Start 1 By 1 Matrix 1st Row upper M EndMatrix prime right-parenthesis Start 1 By 1 Matrix 1st Row bold-italic q EndMatrix Superscript prime Baseline equals Start 1 By 1 Matrix 1st Row bold-italic upper F EndMatrix prime (1.128)

      1.5 Random Process

      As clean harmonic signals are rare in technical systems we need a mathematical toolbox to deal with nonharmonic signals. Due to the Fourier theory (Nelson and Elliott, 1993) every deterministic time signal can be synthesized by a sum of harmonic components. This provides a theoretical link between periodic and harmonic but also transient signals. See appendix A.1 for a brief introduction.

      The same theory can also be applied to stochastic or random signals but in a slightly different manner. Fourier analysis and methods for investigating the random processes and the description of mechanical systems by impulse response or frequency response functions is an important toolset for the description of vibroacoustic systems. In addition and due to digitalisation most signals or spectra are given as a discrete, digital set of values. This discrete formulation creates some pitfalls that may also lead to misinterpretations. Even though the signal analysis is not a major subject in vibroacoustic simulation it is a very important especially when acoustic experiments are performed.

      Signals from technical processes are often not predictable and are generated randomly like pressure fluctuations in a turbulent flow, the impact of raindrops on a roof, or the stochastic combustion in a jet. For dealing with such signals the above described methods must be adapted. In addition we need formulations that allow for the definition of the statistics of random processes as there is no deterministic functionality between time or frequency and the physical quantity, for example force or displacement.

      1.5.1 Probability Function

      Imagine a random process creating signal sequences as shown in Figure 1.19. At each time the signal value f(t) may be different and has a certain continuous value. One option to characterize this signal is to define the probability that the signal value is less or equal to a specific value fk. Thus we define a probability for f(t) to be less than or equal to fk

       upper P left-parenthesis f right-parenthesis equals Prob left-bracket f less-than-or-equal-to f Subscript k Baseline right-bracket (1.129)

      Figure 1.19 Stochastic fluctuation with time. Source: Alexander Peiffer.

      Next, we are interested in the probability that the value of f(t) is in a range defined by Δf=f2−f1 meaning the probability Prob[f1<f≤f2]. Now we can define the probability density function as

       p left-parenthesis f right-parenthesis equals limit Underscript normal upper Delta f right-arrow 0 Endscripts StartFraction Prob left-bracket f 1 less-than f less-than-or-equal-to f 2 right-bracket right-bracket Over normal upper Delta f EndFraction (1.130)

       Скачать книгу