background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
1/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A
Organization (S):
EDF-R & D/MFTT, CS















Manual of Reference
R4.07 booklet: Coupling fluid-structure
Document: R4.07.02



Modeling of the turbulent excitations




Summary:

One describes the modeling of the turbulent excitations available in Code_Aster and the way in which these
last are taken into account in a calculation of dynamics. The turbulent excitations are characterized
by a spectral concentration of efforts, specified using the operator
DEFI_SPEC_TURB
[U4.44.31]. Their catch in
count in a calculation of dynamics is done by projection of the spectrum on the basis of modal structure of which
one wants to calculate the answer. The operations of projection are carried out using the operator
PROJ_SPEC_BASE
[U4. 63.14].
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
2/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A
Count
matters
1
Principle of calculation ................................................................................................................................... 3
1.1
Determination of a modal base of the system under flow and projection of the excitation ....... 3
1.2
Calculation of the response to the turbulent excitation: frequential resolution .......................................... 4
1.2.1
Introduction ............................................................................................................................. 4
1.2.2
Calculation of the interspectres of modal excitations ..................................................................... 4
1.2.3
Calculation of the interspectres of modal answer ......................................................................... 4
1.2.4
Recombination on physical basis ........................................................................................ 5
1.2.5
Statistical elements .............................................................................................................. 5
1.3
Calculation of the response to the turbulent excitation: temporal resolution ............................................. 6
1.3.1
Factorization of the density interspectrale ............................................................................... 6
1.3.2
Generation of the random modal excitations ..................................................................... 6
1.3.3
Amendment of a modal base and projection ........................................................................ 6
1.3.4
Definition of the obstacles ......................................................................................................... 6
1.3.5
Dynamic resolution ............................................................................................................ 6
1.3.6
Projection of Ritz ................................................................................................................... 6
2
Models of turbulent excitation applicable to the telegraphic structures ....................................................... 7
2.1
General principles .......................................................................................................................... 7
2.1.1
Assumptions ............................................................................................................................. 7
2.1.2
Calculation of the interspectres of modal excitations ..................................................................... 7
2.2
Spectra of the type “length of correlation” ..................................................................................... 9
2.2.1
Key words ................................................................................................................................ 9
2.2.2
Definition of the model ............................................................................................................... 9
2.2.2.1
Density interspectrale ................................................................................................ 9
2.2.2.2
Modeling of the spectrum of turbulence by an expression with variables séparées..10
2.3
Model of turbulent excitation distributed .......................................................................................... 17
2.3.1
Key words .............................................................................................................................. 17
2.3.2
Decomposition on a family of functions of form ......................................................... 17
2.3.3
Setting in equations ................................................................................................................ 18
2.3.3.1
Application of a turbulent excitation distributed ....................................................... 18
2.3.3.2
Turbulent excitation identified on model GRAPPE1 ..................................... 18
2.3.3.3
Projection of the excitation on modal basis .............................................................. 19
2.4
Model of localized turbulent excitation ........................................................................................ 21
2.4.1
Key words .............................................................................................................................. 21
2.4.2
Bases .......................................................................................................................... 21
2.4.3
Setting in equations ................................................................................................................ 22
2.4.3.1
Application of a localized turbulent excitation ..................................................... 22
2.4.3.2
Turbulent excitation identified on model GRAPPE2 ..................................... 24
3
Bibliography ........................................................................................................................................ 25
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
3/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A
1
Principle of calculation
1.1 Determination of a modal base of the system under flow and
projection of the excitation
The calculation of the dynamic response of a system to a turbulent excitation induced by a flow
fluid is carried out by respecting the following stages:
1) initially, one calculates the modal base of the system except flow using
the operator
MODE_ITER_SIMULT
[U4.52.03],
2) one defines then the characteristics of the studied configuration, for taking into account of
phenomenon of coupling fluid-structure, using the operator
DEFI_FLUI_STRU
[U4.25.01].
This operator allows for example to inform the profiles speed associated with the areas
of fluid excitation, for configurations of the type “beam of tubes under flow
transverse ". It produces a concept of the type
[type_flui_stru]
intended to be used by
operators implemented downstream in the command file,
the 3) modal characteristics of the system under flow are then calculated using
the operator
CALC_FLUI_STRU
[U4.66.02]. One has at exit a modal base for
each rate of flow,
4) the definition of the turbulent excitation is done then by a call to the operator
DEFI_SPEC_TURB
[U4.44.31]. Modelings available are as follows:
·
spectra of the type “length of correlation”, specific of the configurations of the type
“beam of tubes under transverse flow”, for the application to the vibrations of
tubes of Steam Generator. The key words corresponding factors are
SPEC_LONG_COR_1
,
SPEC_LONG_COR_2
,
SPEC_LONG_COR_3
and
SPEC_LONG_COR_4
. These spectra are
preset; however, the user can adjust the parameters of them. This part is
developed with the paragraph [§2.2],
·
model of turbulent excitation distributed. The key word factor corresponding is
SPEC_FONC_FORME
. The spectrum of excitation is defined by its decomposition on one
family of functions of form while providing, on the one hand a matrix interspectrale, and
in addition a list of functions of form associated with this matrix. Concepts
[interspectre]
and
[function]
associated must be generated upstream. In the case
component “control rod”, the user can also use a spectrum of
turbulence preset, identified on model GRAPPE1. This part is developed with
paragraph [§2.3],
·
model of localized turbulent excitation. The key word factor corresponding is
SPEC_EXCI_POINT
. It is used in the case of a spectrum of excitation associated with one or
several specific forces and moments. The definition of the excitation is done then in
providing:
-
a matrix interspectrale of excitations (the concept
[interspectre]
associated
must be generated upstream),
-
the list of the nodes of application of these excitations,
-
the nature of the excitation applied of each one of these nodes (force or moment),
-
directions of application of the excitations thus defined.
This part is developed with the paragraph [§2.4].
5) The projection of the spectrum of turbulent excitation previously definite, on the basis of modal
structure under flow, is then carried out using the operator
PROJ_SPEC_BASE
[U4.63.14].
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
4/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A
1.2 Calculation of the response to the turbulent excitation
: resolution
frequential
1.2.1 Introduction
The calculation of the frequential response of the structure or the system coupled fluid-structure is done in
three stages:
1) calculation of the interspectres of modal excitations,
2) calculation of the interspectres of modal answer,
3) recombination on the physical basis.
Initially, one introduces for each mode the transfer function of the mechanical system
(structure alone or system coupled fluid-structure). Each of the three stages above is then
detailed.
1.2.2 Calculation of the interspectres of modal excitations
Interspectres of modal excitations
(
)
S
F U
QiQj
,
are determined by projection of the spectrum
of turbulent excitation on the basis of modal system mechanical (structure alone or coupled system
fluid-structure). This stage of projection is detailed in paragraph [§2] for the various models
applicable to telegraphic structures.
1.2.3 Calculation of the interspectres of modal answer
Interspectres of modal displacements
(
)
S
F U
qiqj
,
result then from the interspectres
modal excitations
(
)
S
F U
QiQj
,
using the following relation:
(
)
(
)
(
) (
)
U
F
H
U
F
S
U
F
H
U
F
S
J
QiQj
I
Q
Q
J
I
,
,
,
,
*
=
éq
1.2.3-1
where
(
)
H
F U
I
*
,
indicate the combined complex of the transfer function
(
)
H F U
I
,
system
mechanics considered. Being given a frequency F and a rate of flow
U
, the function of
transfer
(
)
H F U
I
,
mechanical system for mode I is defined by:
H F U
M
F
F
J
F
F
I
I I
I
I
I
(,)
=
- +


+


1
2
1
2
2
éq
1.2.3-2
where
M
I
indicate the modal mass of the mode
I
F
I
I
,
and
indicate respectively, at the speed
U
,
pulsation and the Eigen frequency of the mode
I
I
,
indicate, at the speed
U
, the reduced damping of
mode
I
, and
J
indicate the complex number such as
J
2
1
= -
.
The calculation of the interspectres of modal displacements starting from the interspectres of modal excitations
and of the transfer functions is carried out using the operator
DYNA_SPEC_MODAL
[U4.53.23].
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
5/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A
One deduces in particular from [éq 1.2.3-2] the relation binding the autospectres modal displacements to
autospectres of modal excitations:
(
)
(
)
(
)
U
F
S
U
F
H
U
F
S
QiQi
I
qiqi
,
,
,
2
=
éq
1.2.3-3
where
(
)
H F U
I
,
2
indicate the square of the module of
(
)
H F U
I
,
1.2.4 Recombination on physical basis
Being given a rate of flow
U
, the interspectre of physical displacement
(
)
F
X
X
S
U
U
,
,
2
1
2
1
at the points of X-coordinates
X
X
1
2
and
, at the frequency
F
, is obtained by modal recombination. This
operation is written:
(
)
(
)
S
X X F
X
X S
F U
U U
I
J
NR
I
NR
J
Q Q
I J
1 2
1
2
1
1
1
2
,
,
() ()
,
=
=
=
éq
1.2.4-1
Where
NR
indicate the number of modes of the base;
()
I
K
X
is the component at the point of discretization
X
K
deformation of I
ème
mode following the direction of space considered.
The recombination on physical basis is carried out using the operator
REST_SPEC_PHYS
[U4.63.22]. The direction of space considered is specified at the time of the call to this operator.
1.2.5 Elements
statistics
The modal variance
()
I
U
2
, associated at the speed
U
, expresses itself as follows:
(
)
I
Q Q
U
S
F U df
I J
2
0
2
()
,
=
éq
1.2.5-1
At the rate of flow
U
, value RMS
()
RMS
X
of response in an item X of the structure is
data by:
)
(
)
(
)
(
2
1
2
U
X
X
I
NR
I
I
RMS
=
=
éq
1.2.5-2
Where
NR
indicate the number of modes of the base and
()
I
X
is the component at the point
X
deformation of I
ème
mode following the direction of space considered.
This operation is carried out by the operator
POST_DYNA_ALEA
[U4.84.04].
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
6/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A
1.3
Calculation of the response to the turbulent excitation: temporal resolution
The temporal resolution proceeds according to the sequence of the following operations:
1.3.1 Factorization of the density interspectrale
The operator
FACT_INTE_SPEC
[U4.36.04] the factorization of the density interspectrale realizes
modal excitations
(
)
S
F U
QiQj
,
, before application of the method of Monte Carlo.
1.3.2 Generation of the random modal excitations
The operator
GENE_FONC_ALEA
[U4.36.05] generates random modal excitations
()
Q T
I
in
carrying out hard copies by the method of Monte Carlo. The operator
RECU_FONCTION
[U4.32.03]
allows to recover each evolution
()
Q T
I
.
1.3.3 Amendment of a modal base and projection
The operator
MODI_BASE_MODALE
[U4.66.21] the modal base of the structure in substituent modifies with
initial characteristics those obtained for a rate of flow considered.
The operator
PROJ_MATR_BASE
[U4.63.12] the projection of the matrices of mass and stiffness allows
assembled on the new modal basis previously definite.
1.3.4 Definition of the obstacles
The definition of the geometry of the obstacles is carried out, if necessary, using the operator
DEFI_OBSTACLE
[U4.44.21].
1.3.5 Resolution
dynamics
Transitory dynamic calculation for the mode
(
)
I
I
NR
1
is carried out using a diagram
of numerical integration with the operator
DYNA_TRAN_MODAL
[U4.53.21].
M Q T
C Q T
K Q T
Q T
II I
II I
II I
I
& & ()
& ()
()
()
+
+
=
éq
1.3.5-1
Where
MR. C
K
II
II
II
,
and
indicate respectively the generalized mass, damping and stiffness
associated I
ème
mode;
()
()
Q T
Q T
I
I
and
displacement and the excitation indicate respectively
generalized associated I
ème
mode.
1.3.6 Projection of Ritz
The restitution on physical basis is carried out using a projection of Ritz:
()
U X
U X
,
() ()
T
Q T
I
I
I
NR
=
=
1
éq
1.3.6-1
()
U X, T
indicate the assembled vector of physical displacements;
()
U X
I
is the assembled vector
defining I
ème
modal form and
()
Q T
I
generalized displacement following I
ème
mode.
This last operation is carried out using the operator
REST_BASE_PHYS
[U4.63.21].
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
7/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A
2
Models of turbulent excitation applicable to the structures
telegraphic
2.1 Principles
Generals
2.1.1 Assumptions
One supposes that the linear excitation induced on the telegraphic structure by turbulence of the flow
can be modelized in the form of a stationary process random ergodic Gaussian of
null average. This turbulent excitation thus is entirely characterized by its density
interspectrale
(
)
S
X X
F
1
2
,
,
, where
X
X
1
2
and
are two unspecified points of the beam and
indicate
the pulsation. The turbulent excitation applied to the structure is thus characterized by its density
interspectrale
S
F
.
Moreover, one supposes that the turbulent forces are independent of the movement of the structure.
The turbulent excitation is identified in experiments on a model of reference. It is
then applicable to any real component in geometrical similarity with the model of reference.
2.1.2 Calculation of the interspectres of modal excitations
One indicates by
()
F X S
T
,
linear density of turbulent excitation exerted on the beam;
X
is
the current X-coordinate of a point of the beam and S the complex pulsation (variable of Laplace). They are done
additional assumptions following H1 and H2:
H1. The excited length
L
E
is lower than the overall length
L
beam.
H2. The expression of
()
F X S
T
,
does not depend on the origin of the excited area
X
E
; that is translated
by
()
(
)
F X S
F X X S
T
T
E
,
,
=
-
.
In this case, one can express the linear density F
T
in the following form:
()
F X S
U D C
D
D
D
L S
T
F
H
E
R
,
,
,
, Re
=




1
2
2
éq
2.1.2-1
with:
= -
=
=
X X
L
S
sD
U
UD
E
E
R
Re
Where
indicate the density of the fluid,
U
is the mean velocity of flow of the fluid,
D
and
D
H
are respectively the diameter of the structure and the hydraulic diameter,
C
F
represent it
adimensional coefficient of turbulent force,
X
is the current X-coordinate of a point of the beam,
X
E
indicate the X-coordinate of the origin of the excited area,
L
E
represent the excited length,
is the variable
of space reduced,
S
is the complex pulsation (variable of Laplace),
S
R
is the complex pulsation
reduced,
is the kinematic viscosity of the fluid, finally “
Re
“the Reynolds number indicates.
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
8/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A
By geometrical assumption of similarity of the real component with the model of reference, one
obtains:
()
(
)
F X S
U D C
S
T
F
R
,
, Re
=
1
2
2
éq
2.1.2-2
Thus, the modal turbulent excitation
()
Q S
I
can be written in the field of Laplace (assumption H2):
()
(
)
Q X
F X S
X dx
L
F
L S
L
X D
I
T
I
E
T
E
I
E
E
X
X
L
E
E
E
()
,
()
,
(
)
=
=
+
+
0
1
éq
2.1.2-3
where
()
I
X
is the component of I
ème
modal deformation according to the direction of space in which
acts the turbulent excitation.
By means of the expression [éq 2.1.2-2], one deduces:
(
)
Q S
U DLL C
S
L
X D
I
E
F
R
I
E
E
()
, Re (
)
=
+
1
2
2
0
1
éq
2.1.2-4
The densities interspectrales of modal turbulent excitations are expressed then in the form:
(
)
(
)
S
F U
U DLL
D
U
F
L
X
L
X D D
Q Q
E
T
R
I
E
E
J
E
E
I J
,
,
,
, Re (
) (
)
=


+
+
1
2
2
2
1
2
0
1
0
1
1
2
1
2
éq 2.1.2-4
with
1
I J
NR
,
, where
NR
is the number of modes selected to determine the answer of
structure;
T
: interspectre of
C
F
enter
1
2
and
;
F
fD
U
R
=
: reduced frequency.
Note:
In what follows, one preserves the assumptions H1 and H2 and one notes
(
)
I
F
ij
R
, Re
the integral:
(
)
(
)
2
1
2
1
1
0
1
0
2
1
)
(
)
(
Re
,
,
,
Re
,
D
D
X
L
X
L
F
F
I
E
E
J
E
E
I
R
T
R
ij
+
+
=
éq
2.1.2-5
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
9/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A
Using this notation, the interspectres of modal excitations are written:
(
)
(
)
S
F U
U DLL
D
U I
F
Q Q
E
ij
R
I J
,
, Re
=


1
2
2
2
éq
2.1.2-6
The expression of the autospectres of modal excitations is similar:
(
)
(
)
S
F U
U DLL
D
U I
F
Q Q
E
II
R
I I
,
, Re
=


1
2
2
2
éq
2.1.2-7
2.2
Spectra of the type “length of correlation”
2.2.1 Key words
The key words factors
SPEC_LONG_COR_i
(
I varying from 1 to 4) of the operator
DEFI_SPEC_TURB
[U4.44.31] give access spectra of the type “length of correlation”. These spectra,
specific of the configurations of the type “beam of tubes under transverse flow”, are
preset but the user can adjust the parameters of them.
2.2.2 Definition of the model

2.2.2.1 Density
interspectrale
In the case of spectra of the type “length of correlation”, the density interspectrale characterizing
the turbulent excitation is supposed to be able to be put in a form at separable variables such
that:
(
)
()
(
)
S X X
S
X X
I
1
2
0
0
1
2
,
,
,
=
éq
2.2.2.1-1
In this expression,
()
S
0
represent the autospectre turbulence and
(
)
0
1
2
X X
,
indicate one
function of space correlation defined by:
(
)
0
1
2
2
1
X X
X
X
C
,
exp
=
-
-






éq
2.2.2.1-2
where
X
X
1
2
and
the X-coordinates of two points of observation indicate and
C
represent the length of
correlation.
Four analytical expressions are available in the operator
DEFI_SPEC_TURB
[U4.44.31]. These
expressions correspond each one to a particular representation of
()
S
0
.
The user defines a spectrum of turbulence by choosing one of these analytical forms, of which it
can adjust the parameters.
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
10/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A

2.2.2.2 Modeling of the spectrum of turbulence by an expression with separate variables
·
General case
The function
tt
introduced into the relation is modelized by a form with variables
separated:
(
)
(
) (
)
T
R
N
N
R
N
NR
F
F
S
1
2
1
2
1
,
,
, Re
,
, Re
=
=
éq
2.2.2.2-1
Where
NR
S
indicate the degree of the base of the functions of form
N
N
and
is a function
independent of the variable of space. These two functions are stored in the base of
data and can be selected by the user.
The autospectres of modal excitations are given by [éq 2.1.2-7] while introducing:
(
)
(
)
I
F
L
F
II
R
nor
N
R
N
NR
S
, Re
, Re
=
=
2
1
éq
2.2.2.2-2
with:
(
) (
) (
)
L
L
X
L
X
D D
nor
N
I
E
E
I
E
E
2
0
1
1
2
1
2
1
2
0
1
=
+
+
,
éq
2.2.2.2-3
The principle of calculation is as follows: one calculates the values first of all of
L
nor
2
while realizing
the calculation of the double integrals; one calculates then
(
)
N
R
F, Re
for all the values of
N; one obtains finally the expression of
(
)
S
F U
Q Q
I I
,
using the equation [éq 2.1.2-4].
·
Particular case: model used for the tubes of steam generator
The particular case of the study of the tubes of Steam Generator corresponds to a particular case of the general case
presented previously while posing
NS
= 1
. The interspectre of turbulent excitation between two
points of reduced X-coordinates
1
2
and
is then given by:
(
)
(
)
T
R
C
E
R
F
L
F
1
2
1
2
,
,
, Re
exp
, Re
=
-
-




éq
2.2.2.2-4
where
C
represent the length of correlation of the turbulent forces and
L
E
is the length
excited. In general, one takes
C
about 3 to 4 times the diameter external of the tube.
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
11/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A
Spectra of autocorrelation of modal excitations, in the case of profiles speed and of
density constant, are given by:
(
)
(
)
Re
,
2
1
,
2
2
R
II
E
QiQi
F
I
U
D
DLL
U
U
F
S


=
éq 2.2.2.2-5
with:
(
)
(
)
(
) (
)
I
F
F
L
L
X
L
X
D D
II
R
R
C
E
I
E
E
I
E
E
, Re
, Re
.exp
.
.
=
-
-




+
+
2
1
1
0
1
0
1
2
1
2
éq 2.2.2.2-6
In the general case of profiles of density and unspecified rate of flow,
one a:
(
)
()
() ()
() () () ()
2
1
2
1
2
2
1
2
2
1
1
2
2
exp
2
1
,
dx
dx
X
X
X
U
X
U
X
X
X
X
F
S
U
D
D
U
F
S
I
I
E
E
L
X
X
E
E
C
L
X
X
R
QiQj
E
E
E
E
E
E




-
-


=
+
+
éq 2.2.2.2-7
Where
D
is the diameter of the structure,
L
E
is the length of the excited area,
E
X
is
the X-coordinate of the origin of the excited area,
U
is the mean velocity of the flow,
()
S F
R
is a spectral concentration of separate excitation the mean velocity of the flow
U X
X
,
1
2
and
are the curvilinear X-coordinates of two points of observation on the tube,
()
X
E
is the profile of density of the fluid along the tube,
()
U X
E
is the profile speed
transverse of the flow along the tube and
C
indicate the length of correlation.
Adimensional profiles of density and transverse speed of the flow
external are in the following way defined
:
()
E
X
indicating the evolution of the density of the external fluid along the area
immersed
L
imm
tube, one indicates by
density of the external fluid realized
on the immersed part of the tube
:
()
=
+
1
L
X dx
imm
E
X
X
L
imm
imm
imm
éq 2.2.2.2-8
One indicates by
()
R X
adimensional profile of density such as
E
X
R X
()
()
=
.
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
12/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A
()
U X
E
indicating the evolution rate of flow of the external fluid over the length
excited
L
E
tube, one indicates by
U
rate of flow of the fluid realized on
excited length of the tube:
()
U
L
U X dx
E
E
X
X
L
E
E
E
=
+
1
éq 2.2.2.2-9
One indicates by
()
U X
adimensional profile transverse speed of the external flow,
such as
U X
U U X
E
()
()
=
.
By introducing the average sizes and the adimensional profiles into the expression
[éq 2.2.2.2-7], one obtains:
(
)
()
() () () () () ()
S
F U
U D
D
U.S.F
X
X
X
X U
X U X
X
X dx dx
QiQj
R
C
X
X
L
X
X
L
E
E
E
E
I
J
E
E
E
E
E
E
,
.
exp
=


-
-




+
+
1
2
2
2
2
1
1
2
2
1
2
2
1
2
1
2
éq
2.2.2.2-10
After having noted
= -
X X
L
E
E
, it comes:
(
) (
)
(
) (
) (
) (
)
]
2
1
2
1
2
2
1
2
1
0
1
0
2
1
1
2
2
3
3
2
exp
)
(
4
1
)
,
(
D
D
X
L
X
L
X
L
U
X
L
U
E
X
E
L
R
E
X
E
L
R
X
X
F
S
L
D
U
U
F
S
E
E
I
E
E
I
E
E
E
E
C
R
E
QiQj
+
+
+
+


+
+




-
-
×
=
éq 2.2.2.2-11
Where
()
S F
R
represent the spectrum of turbulence, definite according to a reduced frequency
F
R
(a Strouhal number). For a tube in interaction with a transverse flow,
F
R
is written:
F
fD
U
R
=
where
F
is the dimensioned frequency,
D
is the diameter of the tube and
U
is speed
average of the flow.
The double integral of the expression [éq 2.2.2.2-11] is evaluated by the operator
PROJ_SPEC_BASE
[U4.63.14].
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
13/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A
·
Case of multiple areas of excitation
If there are several areas of excitations, the notations are introduced
additional following:
The area of excitation
K
being identified by its X-coordinate of beginning
X
K
and its length
L
K
, one
note
()
U X
K
profile speed transverse of the fluid flow on the level of this area.
average transverse speed on the area of excitation
K
is then given by:
U
L
U X dx
K
K
K
X
X
L
K
K
K
=
+
1
()
One deduces the adimensional profile from it transverse speed, standardized on the area
K
:
U X
U X
U
K
K
K
()
()
=
K
indicating the total number of areas of excitation, average transverse speed on
the whole of the areas of excitation is defined by:
U
K
U
K
K
K
=
=
1
1
If
V
gap
is the speed intertube at the entry of the Steam Generator (the range speeds retailers is defined
in
CALC_FLUI_STRU
[U4.66.02] using the key word
VITE_FLUI
), one proceeds to one
the second standardization; transverse speed in an item X located in the area of excitation
K
is given by:
V X
V
U X
U
V
U
U U X
K
gap
K
gap
K
K
()
()
()
=
=
Thanks to this standardization, the arithmetic mean transverse speed on all them
areas of excitation is equal at the speed intertube; one has indeed:
1
1
1
K
L
V X dx
V
K
K
X
X
L
K
K
gap
K
K
K
()
+
=


=
The calculation of the interspectres of modal excitations, realized by the operator
PROJ_SPEC_BASE
[U4.63.14], is done by adding the contributions with each area of excitation according to
the relation:
(
)
()
S
F V
D
D
V
L
S F
Q Q
gap
K
ij
K
rk
K
K
I J
,
=
×
×




=
1
2
2
1
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
14/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A
with:
V
V
U
U
K
gap
K
=
×
and
F
fD
V
rk
K
=
and
L
X
X
X
X V
X V
X
X
X dx dx
ij
K
C
E
X
X
L
X
X
L
E
K
K
I
I
K
K
K
K
K
K
=
-
-






+
+
exp
() ()
()
() () ()
2
1
1
2
2
1
2
2
1
2
1
2
that is to say:
L
V
X
X
X
X U X U X
X
X dx dx
ij
K
K
C
E
X
X
L
X
X
L
E
K
K
I
I
K
K
K
K
K
K
=
×
-
-




+
+
4
2
1
1
2
2
1
2
2
1
2
1
2
exp
() () () () () ()
One poses:
L
X
X
X
X U X U X
X
X dx dx
ijk
C
E
X
X
L
X
X
L
E
K
K
I
I
K
K
K
K
K
K
=
-
-






+
+
exp
() () () () () ()
2
1
1
2
2
1
2
2
1
2
1
2
The expression of the interspectres of modal excitations becomes then:
(
)
()
S
F V
D
D
V
V
L
S F
Q Q
gap
K
K
ijk
rk
K
K
I J
,
=
×
×
×




=
1
2
2
4
1
from where:
(
)
()
(
)
S
F V
D
V
L
S F
Q Q
gap
K
ijk
rk
K
K
I J
,
=
×
×
×
=
1
4
3
3
1
·
Analytical expressions of the spectra available for the user
Various analytical expressions of the spectra available in the operator
DEFI_SPEC_TURB
[U4.44.31] are as follows:
·
SPEC_LONG_COR_1
Each speed
U
I
defined by the user by discretizing the range speeds
[
]
U
U
min
max
-
explored is initially standardized in the form
U
ikn
while applying
the equation:
U
U UU
ikn
I
K
=
where
U
K
and
U
respectively indicate the speed realized on the area of excitation “
K
“, and
mean velocity on the unit of the areas of excitation.
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
15/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A
A Reynolds number “local”
R
eik
, associated the area “
K
“and at the speed
U
I
is then
calculated starting from the local characteristics of the flow:
Re
ik
ikn
U
D
=
The spectrum of turbulent excitation associated the area “
K
“and at the speed
U
I
is given under
the shape of a vector
S
ik
, having as many components as of points used for
to discretize the frequential interval
[
]
F
F
min
max
-
, support of the excitation. J-ième
component
S
jik
this vector is provided by the expression:
S
F
F
F
F
jik
rjik
rc
rjik
rc
=
-




+


0
2
2
2
2
1
4
éq
2.2.2.2-12
F
rjik
is provided by:
F
F D
U
rjik
J
ikn
=
where:
F
J
is the value of frequency associated with the j-ième component in the discretization with
the frequential interval
[
]
F
F
min
max
-
,
F
rc
is a cut-off frequency being worth 0.2;
O
,
,
depend on the Reynolds number according to equations' provided in the table
below:

R
eik
O
] -
; 1.5 10
4
]
2.83504 10
- 4
3
0.7
] 1.5 10
4
; 3.5 10
4
]
13 10
20 42 14 10
9 81 10
1197 10
35 95 10
34 69 10
4
4
8
2
12
3
17
4
22
5
.
.
.
.
.
.
-
-
+
-
+


-
-
-
-
-
-
R
R
R
R
R
eik
eik
eik
eik
eik
Idem Idem
] 3.5 10
4
; 5 10
4
] Idem
4 0.3
] 5 10
4
; 5.5 10
4
] 50.18975
10
­ 4
Idem
Idem
] 5.5 10
4
; +
]
Idem
4 0.6
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
16/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A
·
SPEC_LONG_COR_2
The spectrum of turbulent excitation is written:
()
S F
F
F
R
R
rc
=
+
0
1
éq
2.2.2.2-13
The default values of the parameters are as follows:
0
3
15 10
2 7
01
=
=
=
-
.
.
.
F
rc
·
SPEC_LONG_COR_3
The spectrum of turbulent excitation is written:
()
S F
F
R
R
=
0
éq
2.2.2.2-14
with:
0
0
=
=
(
)
(
)
F
F
rc
rc
The default values of the parameters are as follows:
F
rc
= 2
If
F
F
R
rc
, one a:
0
3
5 10
0 5
=
=
-
.
if not
0
5
4 10
35
=
=
-
.
·
SPEC_LONG_COR_4
The spectrum of turbulent excitation is written:
S F
F
R
R
v
()
=
0
éq 2.2.2.2-15
with:
0
2
1
6 8 10
10
=
-
.
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
17/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A
The other parameters are defined by:
=
-
-
-
=
=
With
B
C
D
v
v
v
v
0 5
1 5
2 5
3 5
2
4
.
.
.
.
v
indicate the rate of vacuum;
v
is the volume throughput defined by
v
m
U
=
;
m
is the flow
mass and
U
indicate the mean velocity of the flow. Values of the coefficients of
polynomial in
v
are as follows:
With
B
C
D
=
= -
=
=
24 042
50 421
63 483
33 284
.
.
.
.
2.3
Model of turbulent excitation distributed
2.3.1 Key words
The key word factor
SPEC_FONC_FORME
of the operator
DEFI_SPEC_TURB
[U4.44.31] allows to define
a spectrum of excitation by its decomposition on a family of functions of form. The user with
possibility of defining the spectrum by providing a matrix interspectrale and a list of functions of
form associated. Concepts
[interspectre]
and
[function]
must then be generated in
upstream. In the case of the component “control rod”, the user can also use one
preset spectrum of turbulence, identified on model GRAPPE1.
2.3.2 Decomposition on a family of functions of form
The model of turbulent excitation distributed supposes that the instantaneous linear density of the forces
turbulent
()
F X T
T
,
can be broken up on a family of functions of form
()
K
X
of
dimension
K
in the following way:
=
=
K
K
K
K
T
T
X
T
X
F
1
)
(
)
(
)
,
(
éq
2.3.2-1
Coefficients
()
K
T
at every moment define the decomposition of the turbulent excitation on
family of functions of form.
The density interspectrale of turbulent excitation between two points of the telegraphic structure of X-coordinates
X
1
and
X
2
is written then:
S
X X
X
X S
F
K
L
K
K
K
K
K
L
(,
,)
() ()
()
1
2
1
2
1
1
=
=
=
éq
2.3.2-2
This formulation makes it possible to take into account an excitation whose space distribution is
unspecified.
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
18/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A
2.3.3 Setting in equations

2.3.3.1 Application of a turbulent excitation distributed
The length of application
L
is characterized in an intrinsic way by the field of definition of
functions of form associated with the excitation. The area of application is determined by the data of
name of the node around of which it is centered.
X
N
indicating the X-coordinate identifying this node, the turbulent excitation is imposed on
field
[
]
X
L
X
L
N
N
-
+
2
2
,
.
The turbulent excitation being able, in addition, to be developed in a way correlated in both
directions
Y
and
Z
orthogonal with the axis of the telegraphic structure, the functions of form are a priori
vectors with two components.
One thus informs, by convention, these functions about the interval
[
]
0 2
, L
, the field
[]
0, L
being
associated the direction
Y
and the field
[
]
L L
, 2
being associated the direction
Z
.

2.3.3.2 Turbulent excitation identified on model GRAPPE1
Functions of form
K
are the first 12 modal deformations of bending of the structure
identified in experiments, distributed according to two orthogonal directions' with the principal axis of
beam. The general analytical expression of these deformations is as follows:
R

K
Yk
Zk
X
X
X
()
()
()
=


éq
2.3.3.2-1
with:
Yk
Yk
Yk
Yk
Yk
Yk
Yk
Yk
Yk
X
With
N
L X
B
N
L X
C
CH nL X D HS nL X
()
cos
sin
=




+




+


+


éq
2.3.3.2-2
Zk
Zk
Zk
Zk
Zk
Zk
Zk
Zk
Zk
X
With
N
L X
B
N
L X
C
CH nL X D HS nL X
()
cos
sin
=




+




+


+


éq 2.3.3.2-3
where
N
Yk
and
N
Zk
numbers of waves indicate, L is the length of application of the excitation and them
coefficients
With
Yk
,
B
Yk
,
C
Yk
,
D
Yk
,
With
Zk
,
B
Zk
,
C
Zk
,
D
Zk
are constant real coefficients
characteristics of the function of form considered.
The first 6 functions of form are associated the direction
Y
and
With
Zk
,
B
Zk
,
C
Zk
,
D
Zk
are
thus null, for 1
K 6.
The 6 last functions of form are associated the direction
Z
and
With
Yk
,
B
Yk
,
C
Yk
,
D
Yk
are
thus null, for 7
K 12.
This family of functions of form is thus characterized by 5x12 = 60 real coefficients.
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
19/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A
The turbulent excitation identified on model GRAPPE1 is homogeneous in the two directions
orthogonal with the axis of the telegraphic structure, turbulence being décorrélée between these two directions.
The matrix interspectrale
[
]
S
K
L
identified on model GRAPPE1 is thus a matrix of
dimension 12x12, consisted of two identical diagonal blocks of dimension 6:
[
]
()
[
]
[]
[]
()
[
]
S
S
S
K
O
O
L
=


0
0
By square property of symmetry, this matrix is entirely defined by the data of the part
triangular higher (or lower) of
()
[
]
S
O
, that is to say 21 interspectres. For each one among them, them
characteristic parameters are the level of plate, the cut-off frequency and the slope of the spectrum
beyond this frequency.
The matrix interspectrale of turbulent excitation identified on model GRAPPE1 is thus
characterized by 63 real coefficients (3x21).
Note:
Excitations GRAPPE1 are available to two flows of reference. The whole of
data characterizing these excitations thus represents 246 real coefficients ([60+63] x2).

2.3.3.3 Projection of the excitation on modal basis
One notes:
I
I
I
X
DY X
DZ X
()
()
()
=


the modal deformed i-éme of the structure.
Are
ik
co-ordinates of the modal i-éme deformed of the structure on the basis of the functions of
form
K
X
()
:
I
ik
K
K
K
X
X
()
()
=
=
1
éq
2.3.3.3-1
Interspectres of modal excitations
S
Q Q
I J
()
applied to the structure are written then:
S
S
Q Q
ik
K
K
jl
K
K
I J
K L
()
()
=
=
=
1
1
éq
2.3.3.3-2
For each mode I of the structure, coefficients
ik
are given by integrating the equation
[éq 2.3.3.3-1] prémultipliée by the functions
J
, on the applicability of the excitation. One
obtains as follows:
J
X
L
X
L
I
ik
K
K
J
X
L
X
L
K
X L
X dx
X L
X L
dx
(
)
()
(
)
(
)
+
=
+
+
-
+
=
-
+
2
2
2
0
0
0
0
2
2
1
2
2
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
20/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A
()
J
X
L
X
L
I
ik
K
K
J
L
K
X L
X dx
X
X dx
I J
(
)
()
()
()
,
+
=
-
+
=
2
0
0
2
2
1
0
éq
2.3.3.3-3
For each
I
, the equation [éq 2.3.3.3-3] is written in matric form:
[]
()
()
has
B
jk
ik
ij
=
éq
2.3.3.3-4
with:
has
X
X dx
jk
J
L
K
=
()
()
0
that is to say:
(
)
has
X
X
X
X
dx
jk
Yj
Yk
Zj
Zk
L
=
+
()
()
()
()
0
and
B
X L
X dx
ij
J
X
L
X
L
I
=
+
-
+
(
)
()
2
0
0
2
2
that is to say:
(
)
B
DY X
X L
DZ X
X L
dx
ij
I
Yj
I
Zj
X
L
X
L
=
+
+
+
-
+
()
(
)
()
(
)
2
2
0
0
2
2
The resolution of each linear system of equations leads to
ik
.
The calculation of the scalar products is carried out in the operator
PROJ_SPEC_BASE
[U4. 63.14].
Note:
1)
functions
K
X
()
represent, in practice, the modal deformations raised on
model. The system (), with dominating diagonal, is thus well conditioned. In
private individual, when the telegraphic structure model has a homogeneous linear density, them
functions
K
X
()
are orthogonal and the matrix
[]
has
jk
is diagonal.
2) Tests comparing the applicability of the excitation with the field of definition of
the structure are carried out.
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
21/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A
2.4
Model of localized turbulent excitation
2.4.1 Key words
The key word factor
SPEC_EXCI_POINT
of the operator
DEFI_SPEC_TURB
[U4.44.31] is used in
cases of a spectrum of excitation associated with one or more forces and moments specific. The user
can define the spectrum while providing:
·
a matrix interspectrale of excitations (the concept
[interspectre]
associated must be
generated upstream),
·
the list of the nodes of application of these excitations,
·
the nature of the excitation applied of each one of these nodes (force or moment),
·
directions of application of the excitations thus defined.
It can also use a preset spectrum of turbulence, identified on model GRAPPE2.
2.4.2 Bases
The model of localized turbulent excitation is a particular case of the model of turbulent excitation
distributed. Thus, one supposes just as in paragraph [§2.3.2] that instantaneous linear density
turbulent forces
F X T
T
(,)
can be broken up on a family of functions of form
K
X
()
in the following way:
F X T
X
T
T
K
K
K
K
(,)
()
()
=
=
1
éq
2.4.2-1
Coefficients
K
T
()
at every moment define the decomposition of the turbulent excitation on
family of functions of form.
The density interspectrale of turbulent excitation between two points of the telegraphic structure of X-coordinates X
1
and X
2
is written then:
S
X X
X
X
S
F
K
L
K
L
K
K
K L
(,
,)
()
()
()
1
2
1
1
2
1
=
=
=
éq
2.4.2-2
The characteristic of the model of localized turbulent excitation is due to the specificity of the functions of
form
K
X
()
:
(
)
K
K
X
X X
()
=
-
allows to represent a specific force applied to the point of X-coordinate
X
K
(
)
K
K
X
X X
()
= -
allows to represent one specific moment applied to the point of X-coordinate
X
K
(
)
X X
K
-
and
(
)
-
X X
K
the distribution of Dirac function and the derivative indicate respectively of
distribution of Dirac function at the point of X-coordinate
X
K
.
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
22/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A
Taking into account the specificity of the functions of form, the projection of a localized turbulent excitation
on modal basis is much simpler than in the general case (excitation distributed), since one
can analytically calculate the expression of the projected excitation.
2.4.3 Setting in equations

2.4.3.1 Application of a localized turbulent excitation
One considers a turbulent excitation applied to a structure telegraphic and made up of forces and of
specific moments. This excitation is entirely characterized by the following data:
·
list nodes of application of the forces and moments specific,
·
nature of the excitation applied in each node (force or moment),
·
direction of the excitation applied in each node.
(
)
(
)
Thus
F
N
N
T
K
K
K
K
K
m
m
M
m
m
X T
F S
X X
M
S
X X
(,)
()
()
=
-
-
-
=
=
1
1
éq 2.4.3.1-1
is the expression of a located turbulent excitation, characterized by
K
forces and
M
moments
specific, applied respectively to the nodes of X-coordinates
X
K
and
X
m
in the directions
N
K
and
RN
m
.
One a:
()
()
N
K
K
K
=




0
cos
sin
and
N
m
defined in a similar way.
represent the azimuth giving the direction of application of the force (or the moment) in the plan
P
orthogonal with neutral fiber with the node of application, such as defined in the figure [Figure 2.4.3.1-a]
below:
P
Z
F
X
neutral fiber
Nap
node of application
Appear 2.4.3.1-a: Definition of the direction of application
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
23/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A
The generalized excitation associated I
ème
mode of the structure,
Q S
I
()
, being defined by:
Q S
X
X T dx
I
I
L
T
()
()
(,)
=
0
F
éq
2.4.3.1-2
where
L
represent the length of the beam and
I
X
()
deformation of the mode
I
, one obtains, taking into account
expression [éq 2.4.3.1-1]:
Q S
F S
X
M
S
X
I
K
K
K
I
K
K
m
m
M
I
K
m
()
()
(
)
()
(
)
=
-
=
=
1
1
N
N
éq
2.4.3.1-3
The calculation of the interspectres of modal excitations leads then to:
(
)
(
)
(
)
S
S
F F
I X
J X
F M
I X
J X
MR. F
X
J X
Q Q
S
K
K
K
K
K
K
K
K
S
K
K
m
M
K
K
m
m
S
I
m
m
K
K
m
M
K
I J
K
K
K
m
m
K
S
S
S
()
(
).
.
(
).
(
).
.
(
).
(
).
.
(
()
()
()
'
=




+




+
=
=
=
=
=
=
1
2
2
1
1
2
2
1
1
2
2
1
1
1
1
1
2
2
1
1
1
1
2
2
1
1
1
1
2
N
N
N
N
N
(
)
).
(
).
.
(
).
()
'
N
N
N
K
S
I
m
m
m
M
m
M
m
m
S
MR. M
X
J X
m
m
2
1
1
1
1
2
2
1
2
2
1




+




=
=
éq
2.4.3.1-4
Note:
When the user defines the spectrum of turbulent excitation, it must inform the matrix
interspectrale of the specific excitations whose terms intervene above. This
stamp has as a dimension K+M (a number of forces and specific moments applied).
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
24/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A

2.4.3.2 Turbulent excitation identified on model GRAPPE2
The turbulent excitation identified on model GRAPPE2 is represented by a force and one
moment resulting, applied in the same node following the two orthogonal directions to the axis of
structure. The linear density of this excitation has as an expression:
(
)
() (
)
[
]
F
T
H
p
T
R
p
T
R
X S
U D L
F S
X X
L
M S
X X
(,)
()
=
-
-
-




1
2
0
0
1
2
0
2
0
éq 2.4.3.2-1
Where
is the density of the fluid,
U
is the mean velocity of the flow,
D
H
is it
hydraulic diameter,
L
p
is the thickness of the plate of housing (corresponding to the length
excited),
X
0
is the X-coordinate of the point of application of the excitation,
S
S D
U
R
=
is the complex frequency
reduced,
F S
T
R
()
and
()
M S
T
R
are the adimensional coefficients representing the force and the moment
resulting.
Sizes
,
U
,
D
H
and
L
p
allow to dimension the excitation.
In substituent the expression [éq 2.4.3.2-3] in the relation [éq 2.4.3.1-4] defining the modal excitation
Q S
I
()
, one obtains:
()
() ()
Q S
U D L
F S
X
L
M S
X
I
H
p
T
R
I
p
T
R
I
()
()
'
=




+








1
2
0
0
1
0
0
1
2
0
2
0
éq 2.4.3.2-2
The specific force and moment identified on model GRAPPE2 being décorrélés, the calculation of
interspectres of modal excitations leads finally to:
()
()
()
()
()
()
S
U D
D
U
L
X
X
S
S
L
X
X
S
S
Q Q
H
p
I
J
F F
R
p
I
J
MR. M
R
I J
T T
T
T
=






×






+




×






1
2
0
1
1
0
1
1
0
1
1
0
1
1
2
2
2
0
0
4
0
0
'
'
éq
2.4.3.2-3
In this expression,
D
is the diameter external of the structure,
()
S
S
F F
R
T T
and
()
S
S
MR. M
R
T
T
the adimensional autospectres of force and moment represent respectively identified on
model GRAPPE2. The operator
PROJ_SPEC_BASE
[U4.63.14] calculates the interspectres excitations
modal according to the relation [éq 2.4.3.2-3] above.
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
25/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A
Note:
1) Autospectres adimensional GRAPPE2 are usable to simulate it
behavior of any structure in similarity with the model; one then utilizes
geometrical parameters structural feature to dimension
the excitation. Model GRAPPE2 having been built in similarity with
configuration jet engine, the following reports/ratios are fixed and characteristic of this
geometry:
D
D
L
D
H
p
and
It is reminded the meeting that
D
H
and
D
the hydraulic diameter and it indicate respectively
diameter external of the structure;
L
p
is the thickness of the plate of housing,
corresponding to the excited length.
The data of
,
U
and
D
is thus sufficient to dimension in manner
univocal the turbulent excitation starting from the autospectres adimensional.
2) Adimensional autospectres
()
S
S
F F
R
T T
and
()
S
S
MR. M
R
T
T
one and the other being defined by
three real coefficients (level of plate, reduced frequency of cut and slope beyond
of this frequency), only six constants make it possible to characterize the excitation
turbulent adimensional identified on model GRAPPE2.
Four configurations having been studied (ascending flow or going down, stem of
order centered or offset), the whole of the data characterizing the excitations
GRAPPE2 thus represents 24 real coefficients.


3 Bibliography
[1]
NR. GAY, T. FRIOU: Resorption of software FLUSTRU in ASTER. HT32/93/002/B
[2]
S. GRANGER, NR. GAY: Software FLUSTRU Version 3. Note principle. HT32/93/013/B
[3]
L. PEROTIN, Mr. LAINET: Integration of a general model of turbulent excitation in
Code_Aster: specifications. HT32/96/003/A
background image
Code_Aster
®
Version
7.4
Titrate:
Modeling of the turbulent excitations
Date:
05/04/05
Author (S):
A. ADOBES, L. VIVAN
Key
:
R4.07.02-B
Page
:
26/26
Manual of Reference
R4.07 booklet: Coupling fluid-structure
HT-66/05/002/A

























Intentionally white left page.