Code_Aster ®
Version
2.3
Titrate:
Elements vibroacoustic
Date:
02/11/92
Author (S):
Fe WAECKEL
Key:
R4.02.02-A
Page:
1/12
Organization (S): EDF/EP/AMV
Handbook of Référence
R4.02 booklet: Accoustics
Document: R4.02.02
Vibroacoustic elements
Summary:
The elements described here make it possible to carry out calculations of the frequencies and clean modes of a structure
coupled to a fluid. They also allow the acoustic calculation of answer.
After the formulation of the problem of coupling fluid-structure, this document describes the step followed for
to implement in Aster the new finite elements.
Handbook of Référence
R4.02 booklet: Accoustics
HI-75/7952/A
Code_Aster ®
Version
2.3
Titrate:
Elements vibroacoustic
Date:
02/11/92
Author (S):
Fe WAECKEL
Key:
R4.02.02-A
Page:
2/12
Contents
1 Introduction ............................................................................................................................................ 3
2 vibroacoustic coupling ................................................................................................................ 4
2.1 Presentation .................................................................................................................................... 4
2.2 Formulation of the vibroacoustic problem ..................................................................................... 4
2.2.1 Description of the structure ..................................................................................................... 4
2.2.1.1 Conservation equation of the momentum ........................................... 5
2.2.1.2 Relation of compatibility ............................................................................................ 5
2.2.1.3 Law of behavior in isotropic linear elasticity .................................................. 5
2.2.2 Description of the fluid ............................................................................................................... 5
2.2.2.1 Conservation equation of the momentum ........................................... 5
2.2.2.2 Conservation equation of the mass ...................................................................... 5
2.2.2.3 Law of behavior ................................................................................................. 6
2.2.2.4 Equation of propagation waves or equation of Helmholtz ................................ 6
2.2.3 Description of the interaction fluid-structure ............................................................................ 6
2.2.4 Formulation of the problem coupled ........................................................................................... 6
2.3 Variational equations associated the problem .......................................................................... 7
2.3.1 Variational equations associated the structure ................................................................ 7
2.3.2 Variational equation associated the equation of the fluid .......................................................... 7
2.4 Discretization by finite elements ..................................................................................................... 8
2.5 Choice of an additional variable for the description of the fluid ................................................. 8
2.5.1 Formulation of the new problem ........................................................................................ 8
2.5.2 Variational formulation associated the problem (P') ........................................................... 9
2.6 Discretization by finite elements ..................................................................................................... 9
2.7 Calculations of acoustic answer ..................................................................................................... 10
2.7.1 Speed imposed on the fluid .................................................................................................... 10
2.7.2 Impedance imposed on a wall of the fluid ........................................................................ 10
2.7.3 Discretization by finite elements .......................................................................................... 11
3 Integration in Aster ......................................................................................................................... 11
4 Bibliography ........................................................................................................................................ 12
Handbook of Référence
R4.02 booklet: Accoustics
HI-75/7952/A
Code_Aster ®
Version
2.3
Titrate:
Elements vibroacoustic
Date:
02/11/92
Author (S):
Fe WAECKEL
Key:
R4.02.02-A
Page:
3/12
1 Introduction
The vibratory behavior of a structure is often modified if this one is in the presence of a fluid:
it is what is called the vibroacoustic coupling. One distinguishes the cases from coupling into two
categories: either the fluid is infinite (it is the case of the immersed structures), or the fluid is contained
in a limited medium (it is the case of the tanks more or less filled with fluid).
The finite elements described here make it possible to solve the problems of coupling with a fluid of
finished dimension.
General notations:
P:
instantaneous total pressure in a point of the fluid,
p0:
pressure at rest,
p:
acoustic pressure,
:
instantaneous total density in a point of the fluid,
T
0:
density at rest,
:
acoustic density,
:
density of the structure,
S
X:
acoustic displacement,
U:
displacement of the structure,
:
gradient of acoustic displacements,
, F: pulsation, frequency,
C:
speed of sound in the fluid,
, K: wavelength, numbers wave,
:
tensor of the constraints of the structure,
:
tensor of the structural deformations,
C:
tensor of elasticity of the structure,
T:
tensor of the constraints of the fluid.
Handbook of Référence
R4.02 booklet: Accoustics
HI-75/7952/A
Code_Aster ®
Version
2.3
Titrate:
Elements vibroacoustic
Date:
02/11/92
Author (S):
Fe WAECKEL
Key:
R4.02.02-A
Page:
4/12
2
Vibroacoustic coupling
2.1 Presentation
That is to say an elastic structure defined in a field S which vibrates in the presence of a true fluid, not
weighing, compressible, in isentropic evolution defined in a field F. One indicates by
=
F
S, their common surface.
One notes N, the normal external with the fluid field F.
NS
Sf
N
F
Known
S
At a given moment, the state of the fluid is defined by its field of pressure P and that of the structure by
its field of displacement U.
It is considered that the coupled system is subjected to small disturbances around its state
of balance where the fluid and the structure are at rest.
As follows:
P = p0 + p
U = U (u0 = 0)
The problem of interaction fluid-structure then consists in solving two problems simultaneously:
· one in the structure subjected, on, to a field of pressure p imposed by the fluid
· the other in the fluid subjected to a field of displacement U of the wall.
2.2
Formulation of the vibroacoustic problem
2.2.1 Description of the structure
Assumption:
The structure is homogeneous and obeys the laws of linear elasticity.
Taking into account this assumption, one can write the various following equations controlling the state of
the structure [bib2].
Handbook of Référence
R4.02 booklet: Accoustics
HI-75/7952/A
Code_Aster ®
Version
2.3
Titrate:
Elements vibroacoustic
Date:
02/11/92
Author (S):
Fe WAECKEL
Key:
R4.02.02-A
Page:
5/12
2.2.1.1 Conservation equation of the momentum
The conservation equation of the momentum is written, in the absence of voluminal forces
others that inertias:
d2 U
-
I = 0
ij, J
S
dt2
[éq 2.2.1.1-1]
where:
S
is the density of the structure,
U
is displacement,
ij
is the tensor of the constraints.
2.2.1.2 Relation of compatibility
= 1 U + U
kl
K, L
L, K
2
[éq 2.2.1.2-1]
where kl is the tensor of the deformations
2.2.1.3 Law of behavior in isotropic linear elasticity
= C
ij
ijkl
kl
[éq 2.2.1.3-1]
with the moduli of elasticity Cijkl checking the identities: Cijkl = Cklij = Cjikl = Cijlk
C being the tensor of elasticity.
2.2.2 Description of the fluid
Assumption: the fluid obeys the laws of linear accoustics.
Taking into account this assumption, the equations controlling the state of the fluid are:
2.2.2.1 Conservation equation of the momentum
The conservation equation of the momentum is written, in the absence of sources:
2
D X
T
-
I = 0
ij, J
0
2
dt
[éq 2.2.2.1-1]
where:
T
is the tensor of the constraints in the fluid,
kl
is the density of the fluid in a natural state,
0
X
is the field of displacement of a particle of fluid.
2.2.2.2 Conservation equation of the mass
With the first command and in the absence of acoustic sources, it is expressed by the relation:
p
+
X
0 div
0
T
T =
[éq 2.2.2.2-1]
Handbook of Référence
R4.02 booklet: Accoustics
HI-75/7952/A
Code_Aster ®
Version
2.3
Titrate:
Elements vibroacoustic
Date:
02/11/92
Author (S):
Fe WAECKEL
Key:
R4.02.02-A
Page:
6/12
2.2.2.3 Law of behavior
T = - p
ij
ij
[éq 2.2.2.3-1]
The fluid is supposed in evolution barotrope (the pressure p is, for the fluid given, a function
known only density) p = C
20;
where c0 is the speed of sound in the fluid at rest.
2.2.2.4 Equation of propagation waves or equation of Helmholtz
One deduces it by combination from the conservation equations from the mass [éq 2.2.2.2-1] and from
momentum [éq 2.2.2.1-1] written in harmonic mode, with the pulsation:
2
p + K p = 0
[éq 2.2.2.4-1]
where K =/C is the number of wave.
2.2.3 Description of the interaction fluid-structure
With the interface fluid-structure (), the fluid being nonviscous, it does not adhere to the wall. One thus writes:
· the continuity of the normal constraints:
.n = T .n = - p. N
ij
I
ij
I
ij
I
[éq 2.2.3-1]
· the continuity normal speeds:
dx
i.n =
i.n
I
I
dt
dt
[éq 2.2.3-2]
2.2.4 Formulation of the coupled problem
Ultimately, the formulation of the problem of vibroacoustic in terms of displacements for
structure and of pressure in the fluid led to the equations of the harmonic problem (P):
C
.u
+ 2 U = 0
ijkl
K, lj
S I
in S
[éq 2.2.4-1]
p + k2p = 0 in F
C
.u
.n = - p. N
ijkl
K, L
I
ij
I
on
1
p
U .n
I
I =
2
N on
0
Handbook of Référence
R4.02 booklet: Accoustics
HI-75/7952/A
Code_Aster ®
Version
2.3
Titrate:
Elements vibroacoustic
Date:
02/11/92
Author (S):
Fe WAECKEL
Key:
R4.02.02-A
Page:
7/12
2.3
Variational equations associated the problem
One solves the problem coupled by using the finite element method starting from the weak formulation
problem.
2.3.1 Variational equations associated the structure
Either U, kinematically acceptable in S, the equation [éq 2.2.4-1] can be written in the form
integral:
C
.u
U + 2 U U FD = 0
ijkl
K, Li
I
S
I
I
S
After integration by part, one obtains the weak formulation:
C
.u
U
- 2 U U FD -
C
U U
.ns dS = 0
ijkl
K, L
I, J
S
I
I
ijkl
I
K, L
I
S
Maybe, if one takes into account the boundary condition [eq 2.2.3 - 1]:
C
.u
U
- 2 U U FD -
p U .n dS = 0
ijkl
K, L
I, J
S
I
I
I
I
S
[éq 2.3.1-1]
2.3.2 Variational equation associated the equation of the fluid
That is to say p, kinematically acceptable in F. One writes in variational form the equation
[éq 2.2.2.4-1]:
2
p p + K p p FD = 0
F
div (p
.p) FD -
p. p
FD + K 2
p p
FD = 0
F
F
F
dp
p
.
ds -
p. p
FD + K 2
p p
FD =
0
F
N
F
F
Maybe, if one takes into account the boundary condition [éq 2.2.3-1]
1
p. p - k2 p. p FD - U .p dS = 0
2
N
0
F
[éq 2.3.2-1]
Handbook of Référence
R4.02 booklet: Accoustics
HI-75/7952/A
Code_Aster ®
Version
2.3
Titrate:
Elements vibroacoustic
Date:
02/11/92
Author (S):
Fe WAECKEL
Key:
R4.02.02-A
Page:
8/12
2.4
Discretization by finite elements
The approximation by finished parts of the complete problem leads then to the system:
K 2M
- C
U
= 0
Q
- CT
H
-
p
.2
.c2
0
0
maybe,
K
- C
M
0
U
U
- 2
= 0
p
Q
p
0
H
.CT
0
c2
where:
K and M are the matrices of stiffness and mass of the structure,
H and Q are the fluid matrices obtained respectively starting from the bilinear forms:
p. p FD
p. p FD
F
and F
C is the matrix of coupling obtained starting from the bilinear form
p. U dS
N
F
The choice of the formulation led to a nonsymmetrical system matric what does not make it possible to use
traditional algorithms of resolution.
2.5
Choice of an additional variable for the description of the fluid
2.5.1 Formulation of the new problem
To obtain a symmetrical problem, one associates the variable pressure, a variable
additional. This new variable is, that is to say the potential of displacement of the fluid such as X =
grad [bib1], [bib4], [bib5], that is to say the variable [bib3] such as = -. The variable makes it possible to take
in account directly fluids with variable density. However, it does not represent anything
physically. This is why, the potential of displacements is preferred to him.
One thus replaces displacement X of the fluid by grad in the equations of the problem (P) [§ 2.2.4].
One thus obtains the new problem to be solved (P'):
C
.u
+ 2 U = 0
ijkl
K, lj
S
I
in S
2 + k2p = 0
0
in F
[éq 2.5.1-1]
p = 2
0
in F
[éq 2.5.1-2]
C
.u
.n = - 2 .n
ijkl
K, L
I
0
ij
I on
[éq 2.5.1-3]
U .n
I
I = on
[éq 2.5.1-4]
N
Handbook of Référence
R4.02 booklet: Accoustics
HI-75/7952/A
Code_Aster ®
Version
2.3
Titrate:
Elements vibroacoustic
Date:
02/11/92
Author (S):
Fe WAECKEL
Key:
R4.02.02-A
Page:
9/12
2.5.2 Variational formulation associated the problem (P')
One applies to the equation [éq 2.5.1-1] the formula of GREEN:
2
2
+ K p FD = 0
a.c.
0
F
grad grad
[K 2 p -
2
2
0
] FD +
ds
0
=
0
a.c.
F
F
N
Maybe, if one takes into account the boundary condition [éq 2.5.1-3]
p
FD -
grad grad FD +
U dS = 0
a.c.
[éq 2.5.2 - 1]
N
c2
F
0
F
Moreover, one writes in weak form the equation [éq 2.5.1 - 2] for Q a.c.,
pq
Q
p - 2 Q FD = 0 Q a.c.
FD - 2
FD = 0
0
c2
c2
F
F
0
F
[éq 2.5.2 - 2]
By summoning the equations [éq 2.5.2 - 1] and [éq 2.5.2 - 2], one obtains the associated variational equation
with the fluid:
pq
Q + p
FD - 2
FD -
grad grad FD +
.u dS = 0
0
N
c2
c2
F
0
F
0
F
(Q,) a.c.
[éq 2.5.2 - 3]
2.6
Discretization by finite elements
While proceeding with the same step as that used in [§ 2.3], one is led to the system
matric according to:
U
U
K
O
0
M
0
.M
0
M
M
0
F
0
p
- 2
0
0
fl
p
= 0
c2
c2
0
T
T
M
0
0
0
.M
fl
.H
0
0
c2
Handbook of Référence
R4.02 booklet: Accoustics
HI-75/7952/A
Code_Aster ®
Version
2.3
Titrate:
Elements vibroacoustic
Date:
02/11/92
Author (S):
Fe WAECKEL
Key:
R4.02.02-A
Page:
10/12
K and M being respectively matrices of stiffness and mass of the structure,
M being the matrix of coupling fluid-structure obtained starting from the bilinear form:
U dS and
M,
,
F Mfl H being fluid matrices, respectively obtained starting from the bilinear forms:
p2 FD,
p FD and of
grad 2 FD
F
F
F
2.7
Calculations of acoustic answer
2.7.1 Speed imposed on the fluid
On a part v of the fluid border F, one can impose a limiting condition of normal speed type
v0.
The term of edge of the fluid is written then:
- 2
.u dS = - 2
.u dS + I
v dS
0
N
0
N
0
0
-
F
F
v
v
2.7.2 Impedance imposed on a wall of the fluid
On a part Z of the fluid border F, one can impose a limiting condition of type impedance Z.
p = Z vn
where vn is the outgoing normal speed of the fluid.
By deferring this condition in the equation translating the conservation of the momentum
[éq 2.2.2.1 - 1] and by taking account of the law of behavior of the fluid [éq 2.2.2.3 - 1], one a:
0
grad p + p = 0
Z
[éq 2.7.2 - 1]
To preserve the symmetry of the system, one expresses the equation [eq 2.7.2 - 1] according to the potential of
displacement of the fluid, one a:
3
0
grad! +
= 0
Z 3
T
Maybe, in harmonic:
0
grad + I
= 0
Z
The term of edge of the fluid is written then:
2
2
2
3
0
0
.
ds =
0
.
ds + I
ds
-
F
N
F Z
N
v Z
Ultimately, to impose an impedance on a wall of the fluid amounts introducing into the system one
term of damping.
Handbook of Référence
R4.02 booklet: Accoustics
HI-75/7952/A
Code_Aster ®
Version
2.3
Titrate:
Elements vibroacoustic
Date:
02/11/92
Author (S):
Fe WAECKEL
Key:
R4.02.02-A
Page:
11/12
2.7.3 Discretization by finite elements
If one imposes limiting conditions of imposed speed type or impedance of wall imposed on the fluid,
one is led to solve the following matric system:
K
0
0
M
0
.M
0
S
0
0
0
0
M
- 2
M
- I 3
0
F
0
0
0
fl
0
0
0
= I V
.c2
c2
0
2
0
T
0
0
0
0
.
M
0
0
Q
MT
fl
.H
0
S
0
Z
c2
Q being the matrix obtained starting from the bilinear form:
2 dS and V, the vector obtained from
Z
v dS
0
0
v
.
3 Integration
in
Aster
The elements described previously belong, for the fluid part, with modeling
“3d_FLUIDE” of phenomenon MECANIQUE and, for the interface fluid-structure, with modeling
“FLUI_STRU” of the same phenomenon.
They lead to voluminal or surface elements, for the fluid part, in pressure-potential of
displacement and with surface elements for the interface fluid-structure in potential of
displacement of the fluid-displacement of the structure.
Handbook of Référence
R4.02 booklet: Accoustics
HI-75/7952/A
Code_Aster ®
Version
2.3
Titrate:
Elements vibroacoustic
Date:
02/11/92
Author (S):
Fe WAECKEL
Key:
R4.02.02-A
Page:
12/12
4 Bibliography
[1]
T. BALANNEC, S. COURTIER-ARNOUX, E. LUZZATO, P. THOMAS: Confrontation of tools
numerical of department AMV on a test of coupling fluid-structure. Internal report
E.D.F. D.E.R. HP-61/91.012
[2]
P. GERMAIN: Introduction to the mechanics of the continuous mediums. Masson
[3]
R.J. GIBERT: Vibrations of the structures. Interactions with the fluids. Sources of excitation
random. Collection of Direction of Etudes and Recherche of E.D.F. n°69.
[4]
R. OHAYON, R. VALID: True symmetric formulations off free vibrations off fluid-structure
interaction. Applications and extensions. International conference one numerical method for
coupled systems
[5]
P. THOMAS: Formulation of coupling fluid-structure for the modal analysis in SIVA.
Application to the three-dimensional elements. Internal report E.D.F.D.E.R. HP-35.82/259
Handbook of Référence
R4.02 booklet: Accoustics
HI-75/7952/A
Outline document