Code_Aster
®
Version
7.4
Titrate:
Operator
CALC_INTE_SPEC
Date:
04/02/05
Author (S):
S. CAMBIER, G. JACQUART
Key
:
U4.36.03-G
Page
:
1/6
Instruction manual
U4.3- booklet: Function
HT-66/05/004/A
Organization (S):
EDF-R & D/AMA, EDF-Pole Industry/CNPE of Tricastin
Instruction manual
U4.3- booklet: Function
Document: U4.36.03
Operator
CALC_INTE_SPEC
1 Goal
To calculate a matrix interspectrale starting from a function of time. The matrix is the average
arithmetic of a certain number of matrices interspectrales calculated on various temporal blocks
of a function of the time which one wants to know the interspectre.
Product a concept of the tabl_intsp type.
Code_Aster
®
Version
7.4
Titrate:
Operator
CALC_INTE_SPEC
Date:
04/02/05
Author (S):
S. CAMBIER, G. JACQUART
Key
:
U4.36.03-G
Page
:
2/6
Instruction manual
U4.3- booklet: Function
HT-66/05/004/A
2 Syntax
int [tabl_intsp] = CALC_INTE_SPEC
(
INST_INIT
=
/
II
[R]
/
0
[DEFECT]
INST_FIN
= yew
[R]
DUREE_ANALYZE
=
da
[R]
DUREE_DECALAGE
=
dd
[R]
NB_POIN =
Np
[I]
FUNCTION
=
fo
[function,
formulate]
TITRATE
=
titrate
[l_Kn]
INFORMATION =
/
1
[DEFECT]
/2
);
Code_Aster
®
Version
7.4
Titrate:
Operator
CALC_INTE_SPEC
Date:
04/02/05
Author (S):
S. CAMBIER, G. JACQUART
Key
:
U4.36.03-G
Page
:
3/6
Instruction manual
U4.3- booklet: Function
HT-66/05/004/A
3 Operands
3.1 Operands
INST_INIT
/
INST_FIN
INST_INIT = II
First value of the parameter for which the signals will be used for the calculation of the matrix
interspectrale (urgent initial).
INST_FIN = yew
Last value of the parameter for which the signals will be used for the calculation of the matrix
interspectrale (urgent final).
Note:
The functions will be calculated with the mode of interpolation which was associated to them. It is
advised not to have a problem of discretization which the functions have one
authorized linear interpolation.
F (T)
INST_INIT
INST_FIN
T
DUREE_DECALAGE
fenestrate 2
fenestrate 3
T
1
T
2
T
3
fenestrate 1
DUREE_ANALYZE
Appear 3.1-a: Analyzes and calculation of interspectre on 3 windows with recovery
3.2 Operands
DUREE_ANALYZE
/
DUREE_DECALAGE
DUREE_ANALYZE = da
The functions will be cut out in several windows of duration of analysis
da
. For each one of
these windows a matrix interspectrale is calculated. The matrix interspectrale result of
the operator will be the arithmetic mean of the calculated matrices.
DUREE_DECALAGE = dd
Allows during the function division according to the duration of analysis in windows, to shift
each window one compared to the other one duration
dd
. If
T
K
is the initial moment of
K
ième
fenestrate, the initial moment of
K
1
+
ième
window will be
T
K
+
dd
Code_Aster
®
Version
7.4
Titrate:
Operator
CALC_INTE_SPEC
Date:
04/02/05
Author (S):
S. CAMBIER, G. JACQUART
Key
:
U4.36.03-G
Page
:
4/6
Instruction manual
U4.3- booklet: Function
HT-66/05/004/A
If
[]
[]
X K
y K
and
are two discrete temporal signals.
Are
[]
[]
X K
Y K
and
their transforms of FOURIER discrete, then [bib1] the matrix
interspectrale
[]
[]
[]
[]
[]
S K is worth
S
K
S
K
S
K
S
K
xx
xy
xy
yy
*
where
[]
[]
[]
[]
[]
[]
S
K
p N T
X K
X K
S
K
p N T
X K
Y K
xx
I
p
xy
I
p
.
.
.
.
*
*
=
=
=
=
1
1
1
1
where
N
is the number of points per block,
p
is the number of blocks.
Caution:
This average adapted perfectly to the “real” signals, results of a measurement is not appropriate
not without precaution for functions close to a sine (the frequency of the average must
to be much higher than the frequency of the signal.
Note:
If the treated signals come from operator GENE_FONC_ALEA via possibly it
calculation of a dynamic response (operator DYNA_TRAN_MODAL for example), then it is
advised to treat each hard copy of GENE_FONC_ALEA independently. In this case,
it is necessary to choose durations of analysis and shift equal to the duration of each hard copy
GENE_FONC_ALEA (cf GENE_FONC_ALEA [U4.36.05]).
3.3 Operand
NB_POIN
NB_POIN = Np
A number of points of the parameter for one duration of analysis. For each point the functions will be
calculated according to the type of definite interpolation and prolongation. The number of points must be
a power of 2 (calculation of the fast transform of Fourier).
Note:
If the signals consist of a number (power of two) sufficient of points with one
not constant, it is preferable to choose this number to avoid interpolations which can
to generate artefacts. In particular, if the treated signals come from the operator
GENE_FONC_ALEA via possibly the calculation of a dynamic response (operator
DYNA_TRAN_MODAL for example), this number will correspond to the double of the number of points
informed in GENE_FONC_ALEA key word NB_POIN or obtained by INFO=2 in
GENE_FONC_ALEA (cf GENE_FONC_ALEA [U4.36.05].
3.4 Operand
FUNCTION
FUNCTION
=
List names of the functions of concept of the type
function
, which one wishes to calculate the matrix
interspectrale.
Code_Aster
®
Version
7.4
Titrate:
Operator
CALC_INTE_SPEC
Date:
04/02/05
Author (S):
S. CAMBIER, G. JACQUART
Key
:
U4.36.03-G
Page
:
5/6
Instruction manual
U4.3- booklet: Function
HT-66/05/004/A
3.5 Operand
TITRATE
TITRATE
=
titrate
is the title of the concept
tabl_intsp
to print at the head results [U4.03.01].
3.6 Operand
INFORMATION
INFORMATION
=
Specify the options of impression on the file
MESSAGE
.
1
print the initial frequency, the final frequency and the pitch in frequency.
2
like
1
more for each autospectre and interspectre, a criterion of convergence in
function of the number of random hard copies.(a random hard copy corresponds to a window
of analysis).
4
Phase of checking
It is checked if the number of points
Np
is a power of 2.
5 Remarks
of use
In addition to the concept of the type
tabl_intsp
product, the operator creates the concepts of the type
fonction_C
who
constitute the matrix interspectrale. This square matrix being, definite positive, functions
complexes defining the higher triangular part of the matrix are enough.
These functions can be printed using operator IMPR_TABLE [U4.91.03].
6 Example
FONC1=RECU_FONCTION (RESU_GENE=DYNAMODE, NOM_CHAM=' DEPL',
NOEUD=' N51', NOM_CMP=' DY',
INTERPOL=' LIN')
FONC2=RECU_FONCTION (RESU_GENE=DYNAMODE, NOM_CHAM=' DEPL',
NOEUD=' N52', NOM_CMP=' DY',
INTERPOL=' LIN')
INTERS=CALC_INTE_SPEC (INST_INIT=0.,
INST_FIN=10.24,
DUREE_ANALYZE=1.024,
DUREE_DECALAGE=1.024,
NB_POIN=1024,
FONCTION= (FONCT1, FONCT2,)
)
7 Bibliography
[1]
Note DER HP-61/93-067 - Generation of random signals of spectral concentration given -
G. JACQUART
Code_Aster
®
Version
7.4
Titrate:
Operator
CALC_INTE_SPEC
Date:
04/02/05
Author (S):
S. CAMBIER, G. JACQUART
Key
:
U4.36.03-G
Page
:
6/6
Instruction manual
U4.3- booklet: Function
HT-66/05/004/A
Intentionally white left page.