background image
Code_Aster
®
Version
4.0
Titrate:
AHLV101 Guides wave at anechoic exit
Date:
12/01/98
Author (S):
F. STIFKENS, G. ROUSSEAU
Key:
V8.22.101-A
Page:
1/12
Manual of Validation
V8.22 booklet: Harmonic accoustics
HP-51/96/094 - Ind A
Organization (S):
EDF/EP/AMV
Manual of Validation
V8.22 booklet: Harmonic accoustics
V8.22.101 document
AHLV101 - Guide wave at anechoic exit
Summary:
A rectilinear guide of wave at anechoic exit, with rigid walls, whose propagation medium is air
“normal”, is excited by a harmonic incidental wave, normal with the face of input. One calculates the field of
acoustic pressure of the harmonic response by using the élasto-accoustics formulation in
pressure-displacement-potential of displacements.
The tests relate to 3 different modelings (finite elements élasto-acoustics three-dimensional,
two-dimensional and axisymmetric), they make it possible to validate the matrices of rigidity, mass, impedance and
vector source for 3 modelings.
The result of reference comes from an analytical calculation.
background image
Code_Aster
®
Version
4.0
Titrate:
AHLV101 Guides wave at anechoic exit
Date:
12/01/98
Author (S):
F. STIFKENS, G. ROUSSEAU
Key:
V8.22.101-A
Page:
2/12
Manual of Validation
V8.22 booklet: Harmonic accoustics
HP-51/96/094 - Ind A
1
Problem of reference
1.1 Geometry
wave
plane
X
rigid side surfaces
Z
y
.
With
B
C
D
O
exit
anechoic
Tube with rectangular section:
length:
L = L
X
= 1.0 m
height:
H = L
y
= 0.1 m
width:
L = L
Z
= 0.2 m
Co-ordinates of the points (in m):
With
B
C
D
X
0.
0.
1.00
1.00
y
0.
0.05
0.
0.05
Z
0.20
0.10
0.20
0.10
1.2
Properties of materials
Air:
= 1.3 kg. m
­ 3
C
= 343. Mr. S
­ 1
1.3
Boundary conditions and loading
()
Pressure of normal incidental wave at the entry
with
Frequency
Impedance at the end CD
=
P
P
I
T
P
AP
F
Hz
Z
C
Kg m
S
I
=
=
=
=
- -
0
0
2
1
1 0
500
445 9
* exp
.
.
.
.
background image
Code_Aster
®
Version
4.0
Titrate:
AHLV101 Guides wave at anechoic exit
Date:
12/01/98
Author (S):
F. STIFKENS, G. ROUSSEAU
Key:
V8.22.101-A
Page:
3/12
Manual of Validation
V8.22 booklet: Harmonic accoustics
HP-51/96/094 - Ind A
2
Reference solution
2.1
Method of calculation used for the reference solution
The frequencies of the excitation are rather low and jointly the guide of wave is sufficiently
length compared to its side dimensions so that one limits oneself to the plane waves: the phenomenon is
then identical in all points of a plan of wave, i.e. does not depend on the co-ordinates describing
points of this plan, y and Z for example.
One gives on this assumption the well-known general solution of the equations of accoustics for
the two sizes pressure
p
and acoustic speed
v
:
v
F
G
=
-




+
+




T
X
C
T
X
C
éq 2.1-1
p
F
G
=
-




-
+








C
T
X
C
T
X
C
éq 2.1-2
The guide is supposed to be closed at the end of X-coordinate L on an impedance
Z
L
; there is one
reflection on the level of this impedance, which gives a wave of return
G
.
In each point of the guide, there is then superposition of the two functions
F
and
G
; by definition even
final impedance
Z
L
impose on the point of X-coordinate L, enters
p
and
v
the relation.
p
v
Z
L
L
L
=
In the harmonic case
F
and
G
are written:
F
G
T
X
C
I E
T
X
C
R E
I
T X
C
I
T X
C
-




=
+




=
-




+




where
I
and
R
are determined by the boundary conditions.
In the calculation of the impedance
Z
p
v
=
in any point
X
the variable time this time is eliminated,
in accordance with the calculation even of the impedances and is written:
()
Z X
Z I E
R E
I E
R E
I
X
C
I
X
C
I
X
C
I
X
C
=
-
+
-
-
0
background image
Code_Aster
®
Version
4.0
Titrate:
AHLV101 Guides wave at anechoic exit
Date:
12/01/98
Author (S):
F. STIFKENS, G. ROUSSEAU
Key:
V8.22.101-A
Page:
4/12
Manual of Validation
V8.22 booklet: Harmonic accoustics
HP-51/96/094 - Ind A
The final impedance becomes:
Z
Z I E
R E
I E
R E
L
I
L
C
I
L
C
I
L
C
I
L
C
=
-
+
-
-
0
One calls
Z
C
0
=
iterative impedance.
On the fluid border at the entry of the guide the condition limits of incidental wave type imposed on
P
P E
I
I T
=
0
, is obtained by writing at the border the following linear relation:
p
C
P
N
I
-
=
v
éq 2.1-3
where
v
N
=
v. N
is speed according to the unit normal
N
outgoing of the fluid.
One imposes moreover on the exit of the guide a value of final impedance
Z
Z
L
=
0
who does one of them
anechoic end.
The final impedance is equal to the iterative impedance
Z
0
when
R
=
0
, i.e. when there is not
no the wave of return; one then has a pure travelling wave in the direction of the incidental wave, that is to say:
v
p
=
=
-




-




I E
C IE
I
T X
C
I
T X
C
thus the relation of imposed incidental wave [éq 2.1-3] is written:
(
)
(
)
p
v
p
-
=
= +
= =
C
X
cv X
C IE
N
I T
0
0
2
from where one identifies
2
C IE
P
I T
I
=
; one deduces the expression from it from the travelling wave of pressure in
guide when one imposes
P
I
with the input of the guide:
p
=
=
-
-




P E
P E
I
I
X
C
I
T X
C
2
2
0
background image
Code_Aster
®
Version
4.0
Titrate:
AHLV101 Guides wave at anechoic exit
Date:
12/01/98
Author (S):
F. STIFKENS, G. ROUSSEAU
Key:
V8.22.101-A
Page:
5/12
Manual of Validation
V8.22 booklet: Harmonic accoustics
HP-51/96/094 - Ind A
2.2
Results of reference
Pressure at the points A, B, C, D (for modelings A, B, C).
2.3
Uncertainty on the solution
Analytical solution.
2.4 References
bibliographical
[1]
F. STIFKENS “Introduction into
Code_Aster
of condition limits of incidental wave type in
vibroacoustic - Ratio HP-61/95/026/
background image
Code_Aster
®
Version
4.0
Titrate:
AHLV101 Guides wave at anechoic exit
Date:
12/01/98
Author (S):
F. STIFKENS, G. ROUSSEAU
Key:
V8.22.101-A
Page:
6/12
Manual of Validation
V8.22 booklet: Harmonic accoustics
HP-51/96/094 - Ind A
3 Modeling
With
3.1
Characteristics of modeling
Pressure-potential formulation of elements displacements
“3d_FLUIDE”
(
MEFL_HEXA20
and
MEFL_FACE8
)
vis-a-vis impedance
imposed
vis-a-vis wave
incidental imposed
Z
X
C
With
B
D
y
Cutting =
15
meshs HEXA20 according to the x axis
2
meshs HEXA20 according to the y axis
2
meshs HEXA20 according to the axis of Z
Limiting conditions:
ONDE_FLUI:
(GROUP_MA: Input
NEAR: 1.0)
IMPE_FACE:
(GROUP_MA: Exit
IMPE: 445.9)
Name of the nodes
With = No1
B = No780
C = No751
D = No763
3.2
Characteristics of the mesh
A number of nodes:
471
A number of meshs and types:
60
HEXA20
8
QUAD8
3.3 Functionalities
tested
Controls
Keys
AFFE_MODELE
“MECHANICAL”
“3D”
GROUP_MA
[U4.22.01]
DEFI_MATERIAU
FLUID
RHO
[U4.23.01]
CELE_R
AFFE_CHAR_MECA
ONDE_FLUI
NEAR
GROUP_MA
[U4.25.01]
IMPE_FACE
IMPE
CALC_MATR_ELEM
“RIGI_MECA”
MODEL
[U4.41.01]
“MASS_MECA”
CHAM_MATER
“IMPE_MECA”
CHARGE
“ONDE_FLUI”
CALC_VECT_ELEM
“CHAR_MECA”
MODEL
[U4.41.02]
CHAM_MATER
CHARGE
DYNA_LINE_HARM
[U4.54.02]
background image
Code_Aster
®
Version
4.0
Titrate:
AHLV101 Guides wave at anechoic exit
Date:
12/01/98
Author (S):
F. STIFKENS, G. ROUSSEAU
Key:
V8.22.101-A
Page:
7/12
Manual of Validation
V8.22 booklet: Harmonic accoustics
HP-51/96/094 - Ind A
4
Results of modeling A
4.1 Values
tested
Localization
Sizes
Reference
Aster
% difference
With
p (real)
p (imag)
0.5
0.0
0.499997
1.2 10
­ 5
6 10
­ 4
-
B
p (real)
p (imag)
0.5
0.0
0.499997
1.2 10
­ 5
6 10
­ 4
-
C
p (real)
p (imag)
­ 0.482466
­ 0.131252
­ 0.482352
­ 0.131670
2.4 10
­ 2
3.2 10
­ 1
D
p (real)
p (imag)
­ 0.482466
­ 0.131252
­ 0.482352
­ 0.131670
2.4 10
­ 2
3.2 10
­ 1
4.2 Parameters
of execution
Version: 3.05.10
Machine: CRAY C90
System:
UNICOS 8.0
Overall dimension memory:
8 megawords
Time CPU To use:
64.08 seconds
background image
Code_Aster
®
Version
4.0
Titrate:
AHLV101 Guides wave at anechoic exit
Date:
12/01/98
Author (S):
F. STIFKENS, G. ROUSSEAU
Key:
V8.22.101-A
Page:
8/12
Manual of Validation
V8.22 booklet: Harmonic accoustics
HP-51/96/094 - Ind A
5 Modeling
B
5.1
Characteristics of modeling
Formulation pressure potential of elements displacements
“2d_FLUIDE”
(
MEFLSE3
and
MEFLQU8
)
vis-a-vis
impedance
imposed
vis-a-vis wave
incidental imposed
X
C
With
B
D
y
Cutting =
15
meshs QUAD8 according to the x axis
2
meshs QUAD8 according to the y axis
Limiting conditions:
ONDE_FLUI:
(GROUP_MA: Input
NEAR: 1.0)
IMPE_FACE:
(GROUP_MA: Exit
IMPE: 445.9)
Name of the nodes
With = No1
B = NO3
C = No751
D = No153
5.2
Characteristics of the mesh
A number of nodes:
125
A number of meshs and types:
30
QUAD8
4
SEG3
5.3 Functionalities
tested
Controls
Keys
AFFE_MODELE
“MECHANICAL”
“2d_FLUIDE”
GROUP_MA
[U4.22.01]
DEFI_MATERIAU
FLUID
RHO
[U4.23.01]
CELE_R
AFFE_CHAR_MECA
ONDE_FLUI
NEAR
GROUP_MA
[U4.25.01]
IMPE_FACE
IMPE
CALC_MATR_ELEM
“RIGI_MECA”
MODEL
[U4.41.01]
“MASS_MECA”
CHAM_MATER
“IMPE_MECA”
CHARGE
“ONDE_FLUI”
CALC_VECT_ELEM
“CHAR_MECA”
MODEL
[U4.41.02]
CHAM_MATER
CHARGE
DYNA_LINE_HARM
[U4.54.02]
background image
Code_Aster
®
Version
4.0
Titrate:
AHLV101 Guides wave at anechoic exit
Date:
12/01/98
Author (S):
F. STIFKENS, G. ROUSSEAU
Key:
V8.22.101-A
Page:
9/12
Manual of Validation
V8.22 booklet: Harmonic accoustics
HP-51/96/094 - Ind A
6
Results of modeling B
6.1 Values
tested
Localization
Sizes
Reference
Aster
% difference
With
p (real)
p (imag)
0.5
0.0
0.499997
1.2 10
­ 5
6 10
­ 4
-
B
p (real)
p (imag)
0.5
0.0
0.499997
1.2 10
­ 5
6 10
­ 4
-
C
p (real)
p (imag)
­ 0.482466
­ 0.131252
­ 0.482352
­ 0.131670
2.4 10
­ 2
3.2 10
­ 1
D
p (real)
p (imag)
­ 0.482466
­ 0.131252
­ 0.482352
­ 0.131670
2.4 10
­ 2
3.2 10
­ 1
6.2 Parameters
of execution
Version: 3.05.10
Machine: CRAY C90
System:
UNICOS 8.0
Overall dimension memory:
8 megawords
Time CPU To use:
62.57 seconds
background image
Code_Aster
®
Version
4.0
Titrate:
AHLV101 Guides wave at anechoic exit
Date:
12/01/98
Author (S):
F. STIFKENS, G. ROUSSEAU
Key:
V8.22.101-A
Page:
10/12
Manual of Validation
V8.22 booklet: Harmonic accoustics
HP-51/96/094 - Ind A
7 Modeling
C
7.1
Characteristics of modeling
Pressure-potential formulation of elements displacements
“AXIS_FLUIDE”
(
MEAXFLS3
and
MEAXFLQ8
)
vis-a-vis wave
incidental imposed
X
C
With
B
D
vis-a-vis impedance
imposed
y
Cutting =
15
meshs QUAD8 according to the y axis
2
meshs QUAD8 according to the x axis
Limiting conditions:
ONDE_FLUI:
(GROUP_MA: Input
NEAR: 1.0)
IMPE_FACE:
(GROUP_MA: Exit
IMPE: 445.9)
Name of the nodes
With = No1
B = NO3
C = No151
D = No153
7.2
Characteristics of the mesh
A number of nodes:
125
A number of meshs and types:
30
QUAD8
4
SEG3
7.3 Functionalities
tested
Controls
Keys
AFFE_MODELE
“MECHANICAL”
“AXIS_FLUIDE”
GROUP_MA
[U4.22.01]
DEFI_MATERIAU
FLUID
RHO
[U4.23.01]
CELE_R
AFFE_CHAR_MECA
ONDE_FLUI
NEAR
GROUP_MA
[U4.25.01]
IMPE_FACE
IMPE
CALC_MATR_ELEM
“RIGI_MECA”
MODEL
[U4.41.01]
“MASS_MECA”
CHAM_MATER
“IMPE_MECA”
CHARGE
“ONDE_FLUI”
CALC_VECT_ELEM
“CHAR_MECA”
MODEL
[U4.41.02]
CHAM_MATER
CHARGE
DYNA_LINE_HARM
[U4.54.02]
background image
Code_Aster
®
Version
4.0
Titrate:
AHLV101 Guides wave at anechoic exit
Date:
12/01/98
Author (S):
F. STIFKENS, G. ROUSSEAU
Key:
V8.22.101-A
Page:
11/12
Manual of Validation
V8.22 booklet: Harmonic accoustics
HP-51/96/094 - Ind A
8
Results of modeling C
8.1 Values
tested
Localization
Sizes
Reference
Aster
% difference
With
p (real)
p (imag)
0.5
0.0
0.499997
1.2 10
­ 5
6 10
­ 4
-
B
p (real)
p (imag)
0.5
0.0
0.499997
1.2 10
­ 5
6 10
­ 4
-
C
p (real)
p (imag)
­ 0.482466
­ 0.131252
­ 0.482352
­ 0.131670
2.4 10
­ 2
3.2 10
­ 1
D
p (real)
p (imag)
­ 0.482466
­ 0.131252
­ 0.482352
­ 0.131670
2.4 10
­ 2
3.2 10
­ 1
8.2 Parameters
of execution
Version: 3.05.12
Machine: CRAY C98
System:
UNICOS 8.0
Overall dimension memory:
8 megawords
Time CPU To use:
62.77 seconds
background image
Code_Aster
®
Version
4.0
Titrate:
AHLV101 Guides wave at anechoic exit
Date:
12/01/98
Author (S):
F. STIFKENS, G. ROUSSEAU
Key:
V8.22.101-A
Page:
12/12
Manual of Validation
V8.22 booklet: Harmonic accoustics
HP-51/96/094 - Ind A
9
Summary of the results
The discretization is strong since it is approximately 45 nodes by wavelength. This is why us
let us obtain results of a high precision: pressure calculated by
Code_Aster
at the point it
less favorable differs from the theoretical value from less than 1%.
It should be also noted that all modelings used give identical results.