Professional Documents
Culture Documents
Virtual SEA Mid-Frequency Structure-Borne Noise Modeling Based On Finite Element Analysis
Virtual SEA Mid-Frequency Structure-Borne Noise Modeling Based On Finite Element Analysis
net/publication/273903885
CITATIONS READS
15 497
1 author:
Laurent Gagliardini
Groupe PSA
80 PUBLICATIONS 748 CITATIONS
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
All content following this page was uploaded by Laurent Gagliardini on 20 February 2017.
L. Petrinelli
GECI France-2, Esplanade Grand Siècle-F78000- Versailles
G. Borello
INTERAC-10, Impasse Borde-Basse, ZA la Violette-F31240- L’UNION, [email protected], www.interac.fr
• Increasing frequency lets the structures be hyper- Figure2: Front central part of a minivan floor structure defining 4
sensitive to small uncertainties in material subsystems.
properties, geometrical details. Such a sensitivity,
inherent to mass production objects such as cars,
cannot be ignored in the medium frequency range.
In this case only space and frequency averaged
responses (i.e. energetic responses) can be Figure 3: Matrix of Energetic Transfer Function between the 4 sub-
systems as defined on figure 2, computed (red- 100 Hz-2 kHz) and
predicted with a convenient precision.
measured (blue- 20Hz-10 kHz).
{ }
∆f , Ω i
H ij = +
(1) related to a reinforcing plate, unfortunately weld to the
Re v F j j prototype. This problematic is shown on purpose, in
∆f
order to emphasis the difficulties inherent to a design
process: the body design is possibly changing locally,
where indicates a frequency average, and a resulting in the fact that any modeling, even the most
{ }
∆f Ω accurate, remains uncertain regarding the design
process. Work is in progress in order to account for such
space average, and Pi = Re v +j F j is the input
modeling uncertainties.
∆f
power at the excitation point (subscript + indicates the
conjugate transposed vector or matrix). LIMITATIONS OF CLASSICAL SEA
(a)
(b)
Figure 5: 1kHz third octave averaged velocities map measured over
(a) the floor; bare (a) and trimmed (b).
FLOOR STRUCTURE
Figure 4. Top (a) and bottom(b) views of a minivan floor structure Few iterations of the Virtual SEA process are necessary
mesh to compute a robust sub-structuring. Figures 6 and 8
show these results on the “NVH meshes”. A different
color marks each SEA subsystem. Figures 7 and 9
present the respective AutoSEA1 model deduced from
the two sub-structuring analysis and associated coupling
MEASUREMENTS
schemes.
Experiments have been carried out with an artificial
At 500 Hz the numerical tool predicts a structure divided
excitation (electro-dynamic shaker) applied to the back
into 13 SEA subsystems, whereas 14 subsystems are
hook supporting the exhaust line. Only the panels
detected to model the minivan floor at 2000 Hz: the
vibration were measured. A LASER vibrometer scanning
differences of sub-structuring clearly appear. Moreover it
confirms the number of subsystems increases with the
frequency.
The same DLFs have been used for both analytical and
Virtual SEA modeling. Only two distinct DLF were
considered: one for the frame components and one for
the panels. These DLFs both decrease as a constant
power of the frequency.
TRIMMED FLOOR
Figure 11. Comparisons of velocity for subsystem S47 Figures 15, 16, 17 and 18 present the comparisons
between the experimental data and the Virtual SEA
predictions. Again, the effect of the trim on the vibratory
behavior is well predicted by the Virtual SEA technique.
REFERENCES
Contrary to experimental SEA, where CLF can be This system of equations can produce a few negative
computed independently from DLF, Virtual SEA imposes CLF values in some frequency bands. This indicates that
the equivalent masses to be computed before the CLF. the proposed model –including coupling scheme,
These masses are used to compute the Energy Transfer estimated subsystems energy and SEA assumptions- is
matrix using expression (A3). not physical. This difficulty is overpass first by taking
advantage of the multiple excitation of each subsystem,
The CLF calculation is performed from an analysis of and formulating a least square problem, and second by
exchanges between subsystems. The balance equation using a Monte-Carlo inverse process as mentioned for
for subsystem i, subjected to excitation, is expressed as: the equivalent masses determination.