Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
1/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
Organization (S):
EDF-R & D/AMA
Instruction manual
U2.06 booklet: Dynamics
Document: U2.06.07
Interaction ground-structure (ISS) in seismic analysis
with the Code_Aster interface - PROMISS 3D
Summary:
This document is a note of description of use of the Code_Aster interface - PROMISS 3D to treat them
problems of interaction ground-structure (ISS) in dynamic analysis: vibrations, seism… One presents to it like
case of application the standard case of a ship on erasing common subjected to a seismic excitation treaty with
various representations of the foundation: rigid or flexible with in this last case the taking into account
total or reduced of the modes of deformation of the foundation raft.
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
2/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
1 Introduction
Processing in seismic analysis of the problem of the behavior of ships posed on flexible foundation
require the taking into account of the interaction between the ground and the structure. However, traditional tools for
to treat this interaction (PARASOL and CLASSI) authorize only rigid foundations, or even
take into account only exclusively homogeneous grounds and foundation rafts of form not
unspecified (case of PARASOL). This is why code PROMISS 3D, developed at the Central School of
Paris, was chained by a procedure established in Code_Aster to make it possible to modelize with
time of the flexible foundations, the heterogeneous grounds - with an extension particular to the laminated grounds -
and of the foundation rafts of an unspecified form, and thus to allow calculations of dynamic interaction with one or
several unspecified structures.
On the one hand, the modeling of the structure of the ship as well as the loadings which are to him
applied, is realized with Code_Aster, and in addition, it is necessary to carry out the analysis stresses
dynamic obtained starting from the characteristics of the elements of the structure modelized using
this code. It is thus then necessary to constitute an interface between PROMISS 3D and Code_Aster
to connect the two preceding operations with the calculation of the linear dynamic evolution of
the ground-ships unit carried out by PROMISS 3D.
This document has thus as a matter to describe this interface consisted modules of calculation
developed around PROMISS 3D and by new specific controls of Code_Aster. One y
fact as a preliminary a simplified description of software PROMISS 3D, complete and detailed description
of its principle being made in the user's manual of PROMISS 3D - MISS 2D [bib1]. One presents to it
the case of standard application of a nuclear small island subjected to a seismic excitation treaty with different
case of foundation: rigid or flexible with or without reduction of modes of deformation of the foundation raft.
2
Description and principle of software PROMISS 3D
Software PROMISS 3D makes it possible to deal with the problems of propagation of wave in fields
rubber bands or fluids.
It uses the geometrical assumption of linearity and behavior: that is to say the equation of Navier
(conservation of the momentum) with the law of Hooke in the springy media and
the equation of the waves in the fluid environments.
This assumption makes it possible to apply a transformation of Fourier compared to the temporal variable
for the whole of the fields to be calculated and thus, to operate the resolution in the field of
frequencies. The return in the temporal field is carried out in postprocessing by the transformation of
Opposite Fourier.
Lastly, software PROMISS 3D rests on a method of under-structuring: the field of study is
broken up into under-fields coupled between them by interfaces. One applies a method to it of
resolution multi-fields and only the interfaces between fields require to be with a grid by
finite elements of border.
The resolution is carried out then on the borders of the under-fields and is founded on knowledge
elementary solutions, functions of Green, fields generated in an infinite field by one
specific stress. One can thus treat the case of the not limited fields, by avoiding any reflection
parasite on fictitious borders truncating the field of study. Moreover, one original extension and
economic was brought to the method by the introduction processing of laminated fields
implicitly taking into account various homogeneous layers of a field without having recourse to
a mesh of their interfaces.
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
3/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
3
Principle of the Code_Aster interface - PROMISS 3D
3.1
Case of the external field in PROMISS 3D
In the case of one or more ships subjected to a seismic excitation which one wants to study
the interaction ground-structure (ISS), one separates by an interface
the field of the “structure”
including/understanding the ships (but also possibly parts of ground not laminated like
fill) of the field “ground” (either laminated, or homogeneous or comprising even fluid parts
(e.g.: reserve of stopping)) modelized directly by PROMISS 3D [Figure 3.1-a].
Structure
small island
Interface
Ground layer 1
Ground layer 2
Appear Model 3.1-a: of interface ground - structure
The structure modelized by Code_Aster is regarded as an external field for
PROMISS 3D. In this case, one breaks up a displacement in this field on modes which,
reduced to the interface, can be null i.e. the dynamic clean modes of the structure on
base fixed
, or not null, i.e. static modes
:
U
I
I
J
I
has
B
=
+
The coefficients has
I
and B
J
are respectively the factors of participation of the dynamic modes and
statics. M and K are respectively the assembled matrices of mass and rigidity. Then, the writing
balance of the field “structure” within the meaning of virtual work provides the following system:
dd
ds
ds
S
dd
ds
ds
S
D
S
K
K
K
K
M
M
M
M
F
F
F
has
B
-
=
+
2
The matrices K
dd
and K
S
are the assembled rigidities projected respectively on the modes
dynamic
and statics
:
T
K
and
T
K
.
The matrices M
dd
, M
S
are the assembled masses projected respectively on the modes
dynamic
and statics
:
T
M
and
T
M
.
K
ds
and M
ds
are the link-words or products cross:
T
K
and
T
M
.
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
4/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
F
D
and F
S
are respective projections on the dynamic modes
and statics
vector forces
applied F to the field “structure”:
T
F and
T
F.
F
represent the action with the sign close to the field “structure” on the interface and can thus express itself with
to leave the preceding system, by eliminating the factor has, by: F
= F
eq
+ K * B
with F
eq
the vector of equivalent forces and K * the matrix of impedance of the field given by:
F
eq
= - F
S
+ (K
ds
-
2
M
ds
)
T
(K
dd
-
2
M
dd
)
- 1
F
D
K * = (K
S
-
2
M
S
) - (K
ds
-
2
M
ds
)
T
(K
dd
-
2
M
dd
)
- 1
(K
ds
-
2
M
ds
)
It is shown simply that one can be reduced to the resolution of 2 local problems to obtain
F
eq
and K *. Indeed F
eq
, solution of the 1
er
local problem, is the effort applied to the locked interface when one
apply the forces to the field “structure”. And when, in the 2
Nd
local problem, one imposes
displacements of unit static modes B on the interface without forces applied to the field
“structure”, one obtains K * B like effort applied to the interface.
The resolution is carried out on (S) the interface (S) which must (or must) be with a grid (S) with elements
surface directed towards the interior of the field “ground” and on this (or these) interface (S) it must y have
between 6 and 10 nodes per wavelength. In the field of the seism one is interested in a range of
frequency going of 0.1 Hz to 30 Hz.
3.2
Contents of the interface
The Code_Aster chaining - PROMISS 3D requires the transfer of following information.
In the Code_Aster direction towards PROMISS 3D, one transfers successively:
·
information concerning the mesh of the interface ground-structure (nodes and elements
surface) as well as the static modes of interface and the dynamic modes of the structure
reduced to the nodes of this interface and ordered according to its local classification,
·
information on the basis of static and modal dynamic mode independent of
local classification: modal masses and rigidities, modal factors of participation and
other terms of coupling between the dynamic modes
and static modes
and, for
each load interns with the structure, projections of the assembled vector corresponding, on
static and dynamic modes as well as the multiplicative function of time associated.
In direction PROMISS 3D towards Code_Aster, one recovers:
·
Evolutions of results of displacements, speeds and accelerations generalized (standard
“
TRAN_GENE
“)
(T),
'(T),
'' (T) on the one hand, and
(T),
'(T),
'' (T) in addition, projected respectively
starting from the dynamic modes
and of the static modes
. One will be able to then project, then
to combine these results on the physical basis: U (T) =
.
(T) +
.
(T).
In the case of a harmonic calculation one recovers complex evolutions by frequency
(
)
and
(
) (of type “
HARM_GENE
“) before projecting them to produce an evolution of it
harmonic on the physical basis: U (
) =
.
(
) +
.
(
).
(
) and
(
)
are always
complexes because of the form complexes impedance of ground: K.
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
5/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
3.3
Procedure of sequence of the programs
An automatic procedure was made up in order to encapsulate the programs necessary to
Code_Aster chaining - PROMISS 3D.
The user must first of all constitute the data of Code_Aster for the preprocessing like
those to describe the stratifications of the ground in module DOS2M of PROMISS 3D in order to calculate
functions of Green.
A first program, gtaster, starting from the file result of Code_Aster, generates the data of
PROMISS 3D including/understanding: mesh of the interface (suffix .mail), design data (.in), them
impedances of the structure (.ext), loadings on the structure (.cext), data of
postprocessing (.post).
A second program, ptaster, recover the movements of the interface calculated by PROMISS 3D in
displacements, speeds and accelerations generalized and transmits them to Code_Aster for
post-to treat.
It is possible in the same tender successively to launch the programs gtaster,
PROMISS 3D, ptaster on the central machine of processing by successive calls to
EXEC_LOGICIEL
contents in a macro-control of Code_Aster called
MACRO_MISS_3D
[U7.03.11].
4
Use of the Code_Aster interface - PROMISS 3D
The interface Code_Aster-PROMISS 3D follows the following diagram [Figure 4-a]:
Appear 4-a: the main stages of calculations of ISS
Signals of seismic excitation
IMPR_MISS_3D
Nom_étude.raster
Code_Aster
-
Calculation of the modes
dynamic
in base
embedded.
- Calculation of the modes
statics
-
M, K, C projected on the basis
modal:
and
- Mesh of the interface.
Code PROMISS 3D
Data of ground
Response of the structure,
determination of the factors
of participation modal
With
and
B
Code_Aster
Recombination
modal
SRO and functions
of transfer
Parameters of calculations
Nom_étude.sol
Nom_étude.optmiss
IMPR_MACR_ELEM
LIRE_MISS_3D
Order of Code_Aster
Data file PROMISS 3D
Representation
small island + foundation
Accélérogrammes
compatible with
spectrum of ground
Evolution in displacement, speed
and/or acceleration of the answers
MACRO_MISS_3D
RECU_FONCTION
CALC_FONCTION
1
era
Stage
2
čme
Stage
3
čme
Stage
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
6/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
4.1 Controls
Code_Aster for the preprocessing
Before using the controls of Code_Aster specific to the coupling, it is necessary to pass
by the following stages of modeling of the structure:
·
calculation of the clean modes by the control
MODE_ITER_SIMULT
[U4.52.03],
·
definition of the dynamic interface, producing the type (
CRAIGB
(recommended) or
MC-NEAL
) and it
mesh of the interface ground-structure, by the control
DEFI_INTERF_DYNA
[U4.64.01],
·
definition of a modal base supplements by the control
DEFI_BASE_MODALE
[U4.64.02]:
it calculates the static modes of the dynamic interface previously definite and complete
the base of the clean modes if the base is of type
CONVENTIONAL
. It is also possible of
to calculate unspecified static modes directly without dynamic interface (in it
case one does not use the control
DEFI_INTERF_DYNA
) by the control
MODE_STATIQUE
[U4.52.14] and to define a modal base then supplements of type RITZ by the control
DEFI_BASE_MODALE
. The interest of the modal bases of RITZ is to be able to assemble
modes calculated with boundary conditions different. For example, modes
dynamic calculated in base embedded with all ddls of the interface locked and them
static modes of interface calculated in various ways:
-
maybe with a solid condition of connection which confers on the interface a movement of body
rigid,
-
maybe with modes of unspecified interface of pace calculated like clean modes
(via the control
MODE_ITER_SIMULT
) of the structure on carpet of springs of ground;
·
assembly of the dynamic macronutrient by the control
MACR_ELEM_DYNA
[U4.65.01], with
to leave the modal base previously definite.
The data to be transferred from Code_Aster towards PROMISS 3D are obtained by the use of
2 specific controls:
·
the first control
IMPR_MACR_ELEM
[U7.04.33] allows to produce the mesh of
the interface ground-structure and modes static and dynamic reduced to this interface. These
terms are used to establish the contribution of the structure on the impedance.
·
evolution in time or frequency whose user will provide the moments or the frequencies
of restitution.
If the evolution is frequential, the frequencies of calculation, defined in PROMISS 3D, must
to be coherent with that defined in
IMPR_MISS_3D
(generally one uses the range of
frequency of study between 0 and 30 Hz at the time D `studies seismic).
If the evolution is temporal, the range of restitution in time defined in
IMPR_MISS_3D
must be included/understood in the duration of the seismic excitation.
One can, for the calculation of this evolution, to define at the same time loadings in the structure
by the key word factor
EXCIT
, including/understanding the vector assembled agent with each one of
these loads, and of the loadings coming from the ground (signals of seismic excitation) by
key word factor
EXCIT_SOL
, including/understanding the type of excitation and the direction of each one of these
loads. In each one of these key words factors, one also gives the signal in frequency
or in time associated with the definite load. Loadings given by
EXCIT
or
EXCIT_SOL
in the same call to
IMPR_MISS_3D
are combined in only one evolution calculated by
PROMISS 3D.
IMPR_MISS_3D
is thus répétable to obtain several calculated evolutions.
One initially encloses the data file of Code_Aster by the word
END
then one
carries out. One creates thus concepts results (modes and dynamic macronutrient
in particular) on a basis of data stored on the central machine of processing.
One thus gives the hand to PROMISS 3D in order to calculate by a resolution in the field
frequencies each evolution previously definite. According to the strategy of
restitution (in time or frequency) of this evolution, one will apply or not one
opposite transformation of Fourier.
It is possible to treat in the same tender without writing on a basis of data
controls of preprocessing, the launching of PROMISS 3D by
MACRO_MISS_3D
[U7.03.11] and controls of postprocessing. But except for the problems of small
cut, it is strongly advised, for reasons of size memory and control of
parameter time in the classes of tender, to continue to split the study into 3
time and to work with a data base Aster.
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
7/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
4.2 Controls
Code_Aster for postprocessing
The control
LIRE_MISS_3D
[U7.02.31] allows to recover an evolution of its choice among
those calculated by PROMISS 3D, the choice being done by the data of a logical unit. It is
necessary to point out the type of the evolution, transient or harmonic. In this last case, one
recover at the same time the real part and the imaginary part data for each frequency of calculation of
displacements, speeds and accelerations generalized. By the data of the modal base supplements
projection, via the dynamic macronutrient, one then obtains the transitory or harmonic evolution on
the physical base of the structure modelized by Code_Aster.
One can then carry out the conventional postprocessing of a seismic study:
·
Extraction of the temporal evolutions of fields of acceleration or displacement to various
levels of ground or structure by the control
RECU_FONCTION
[U4.32.03].
·
Calculation of the spectra of response in these same levels of ground or structure by
order
CALC_FONCTION
[U4.32.04] and the operand
SPEC_OSCI
.
4.3
Controls specific to MISS 3D
4.3.1 Controls
It is possible, except standard use of the interface Code_AsterPROMISS 3D, to use tools of
calculations specific to MISS 3D and to recover the useable results or not by Code_Aster.
Thus, while acting on the files of preprocessing PROMISS 3D (cf [§3.3] and [Figure 4-a]) related to
optimization of calculations (extension .optmis) and/or with the design data (extension .in), one can
to profit from the following options (cf [§ 5.1.4] to have examples of use):
·
Calculation with variable pitch of frequency. The resolution of the equation of the waves (cf [§ 2]) is done
then in the field of the frequencies with a more or less coarse pitch according to tapes'
frequencies. That makes it possible to refine around the interesting frequencies and to be less
precis elsewhere [§ 5.1.4.1].
·
Definition of points of control. The points of control make it possible to recover
information, in particular on the incidental fields and the fields diffracted by (S)
interface (S) (cf [Figure 3-a]), anywhere in the ground.
For that one must write a file of instructions MISS 3D (his name and its extension is
completely free) which will make it possible to extract towards an output file from the evolutions from
fields incidental or diffracted for each point of control and in each direction of
space [§ 5.1.4.2].
·
In the case of buried foundation (cf [§ 5.1.3.2]) of fictitious resonances appear in
certain configurations: soft grounds, wide foundation raft. They are due to the resonance of
started from ground excavated and are located at a frequency close to
H
V
F
p
fictitious
4
=
where
p
V
is
the speed of the wave of compression and
H
depth of the excavation.
An option makes it possible to be freed some by using control RFIC in the file related to
the optimization of calculations (extension .optmis) [§ 5.1.4.3].
·
Simultaneous calculation of the impedances of ground and a transitory and/or harmonic answer of
structure.
For that one must use a file of instructions MISS 3D which will make it possible to write in one
output file values of the impedances of ground or the seismic forces according to
frequency for all ddls of the interface [§ 5.1.4.4].
4.3.2 The files MISS 3D
Two files of MISSES 3D, located in the index associated with the study on the server dedicated to
Code_Aster and with software PROMISS 3D, are interesting to check or control which types of
calculations are carried out at the time of the study and time associated with each controls MISS 3D.
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
8/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
4.3.2.1 The data file main
It has necessarily the suffix .in (nom_étude.in for example) and contains all the controls
used by MISS 3D at the time of the study. It can call upon auxiliary data files. It
break up schematically into three parts:
·
definition of the data,
·
stages of calculation,
·
postprocessing,
these various phases being able to be connected and repeated by respecting the logic of the program.
The execution can be carried out in several phases with resumptions of the various stages of calculation
4.3.2.2 Auxiliary data files
The whole of the data necessary to the definition of a complex problem led to a file of
order big size in which the hierarchy of information tends to disappear. Of
more, often of similar calculations data files have which different only from some
lines, common parts being able to be consigned in the same file. In order to allow of such
cuttings, it is possible in certain menus to disconnect the reading of the data on a file
auxiliary by means of key word FICP.
4.3.2.3 The output file
It has necessarily the suffix .out (nom_étude.out for example) and gives an indication of the unit
controls read by the program, as well as the times CPU spent in each phase
of calculation. The information printed during the various phases of the program is detailed
for each key word.
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
9/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
5
Case of application of taking into account of the ISS by the interface
Code_Aster - PROMISS 3D
To take into account the ISS that amounts representing the ground by a mechanical system are equivalent. Two
methods are currently used:
·
lawful method within the competences of ground for which stiffnesses of the system of
springs of ground independent of the frequency and are adjusted on the first mode of
swinging and first mode of pumping of the system coupled ground-ship (S),
·
frequential method of coupling where the impedance of the ground evolves/moves according to
frequency.
It is the frequential method of coupling which is implemented in PROMISS 3D [bib1] for
to modelize the ISS. However, this method makes it possible to determine the stiffnesses of the system of springs
of ground of the lawful method thanks to a specific option of the chaining
Code_Aster/PROMISS 3D (option
MISS_IMPE
control
MACRO_MISS_3D
). In this case, it
calculation follows the following diagram:
Mesh
Order of Code_Aster
Calculation of the dynamic modes
structure in embedded base
(Dx=Dy=Dz=Drx=Dry=Drz=0)
DEFI_BASE_MODALE
MACR_ELEM_DYNA
Impression of the data for
the calculation of ISS per MISS 3D
AFFE_CHAR_MECA
MACRO_MATR_ASSE
POST_ELEM
MACRO_MODE_MECA
LIRE_MAILLAGE
DEFI_MATERIAU
AFFE_MATERIAU
AFFE_MODELE
AFFE_CARA_ELEM
Calculation of the 6 static modes
structure in rigid foundation
(LIAISON_SOLIDE) and PO blocking
(Dx=Dy=Dz=Drx=Dry=Drz=0)
Definition bases modal
and
Projection of M, K, C on the basis
Launching of MISS 3D
IMPR_MACR_ELEM
IMPR_MISS_3D
AFFE_CHAR_MECA
MACRO_MATR_ASSE
MODE_STATIQUE
MACRO_MISS_3D
(OPTION: MISS_IMPE)
End
Definition of the model
Appear Synoptic 5-a: of the calculation of the stiffnesses within the competence of equivalent ground
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
10/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
Method of adjustment stiffnesses of the system of springs of ground on the first mode of
swinging and the first mode of pumping of the system coupled ground-ship (S) is described in
document [bib3].
As case of application of the Code_Aster chaining/PROMISS 3D, one takes as example one
complex structure [Figure 5-b], [Figure 5-c] resting on a cruciform foundation raft [bib4].
The interest of this case is that it makes it possible to consider the various modes of representation of the interface
ground structure. Thus, the foundation can be considered either rigid, or flexible with the totality of the modes
statics, is flexible with some modes of foundation chosen according to a method of reduction
modal [bib2].
Appear 5-b: Mesh of the structure
Appear 5-c: Mesh of the foundation of the structure
Section E-W
Section NS
Structures
Interns
Chambers
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
11/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
5.1
chaining
Code_Aster - PROMISS 3D: Stages and parameters
Calculations are carried out by implementing the Code_Aster chaining - PROMISS 3D. The different ones
parameters and stages necessary to the calculation of ISS are described below:
5.1.1 Data transmitted by Code_Aster (cf [Figure 4-a])
Code_Aster transmits to PROMISS 3D: the mesh of the interface enters the ground and the structure (face
lower of the foundation raft), the modal base combining the dynamic modes
structure and modes
statics
as well as projection on this basis of the matrices M, K and C. This operation is carried out
via the operator
IMPR_MACR_ELEM
of Code_Aster.
Code_Aster transmits also the data relating to the seismic excitation
: they are them
accélérogrammes acc1, acc2 and acc3. That is done thanks to the control
IMPR_MISS_3D
Code_Aster.
The whole of the data transmitted via the operators
IMPR_MACR_ELEM
and
IMPR_MISS_3D
are
printed in the file nom_étude.raster result of the type:
libr
on the unit
26
by defect.
5.1.2 Data specific to PROMISS 3D for its calculation in the field
frequential
PROMISS 3D is based on the assumption of linearity as well from the geometrical point of view as of the point of
sight of the behavior of the materials [bib1]. This assumption of linearity makes it possible to solve them
problems in the frequential field. It is in the file: nom_étude.optmiss, which is in
the index associated with the Aster study, which one defines the parameters necessary to calculations in
field of the frequencies. One finds there in particular the range of frequency [Fmin, Fmax] in which
will be carried out the calculation and the pitch of sampling dF. The file nom_étude.optmiss is given in
appendix 1 of the document.
The rules of adjustment of the parameters are recalled in the document [bib1]. For our part us
took the following parameters for the study of a nuclear small island:
Fmax
The range of frequency depends on the nature of the ground. The maximum frequency
reserve is 20 Hz (ground means-slackness).
dF = 0.1 Hz
The sampling rate retained is identical to that of
accélérogrammes acc1.c2, acc2.c2 and acc3.c2 used for the seismic excitation.
Fmin = 0.1 Hz
This parameter is a function of the two precedents. It is necessary that the report/ratio
dF
F
F
Max
min
-
that is to say an entirety.
Z0 =-11.60 m
One must re-enter the dimension of the base of the foundation. For our study the foundation is
surface. Axis OZ of the model must always be vertical and the normals with
plans of the mesh of the foundation are obligatorily directed towards the interior of
field of the ground.
5.1.3 Data relating to the ground
It is in the file nom_étude.sol, which is in the index associated with the Aster study, that one
described the data relating to the ground. The constitution of the ground laminated there is indicated like
the site of the chamber of the seismic excitation and parameters of sampling of the functions of
Green. The file nom_étude.sol is given in Appendix 1 of the document.
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
12/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
5.1.3.1 Data of ground
They are given layer by layer. One described the mechanical characteristics there (Young modulus,
Poisson's ratio, density, reduced damping) of materials constitutive of
layers and their thicknesses.
Recalls:
·
The propagation velocity of the waves of compression is given by
:
(
)
G
2
1
1
2
-
-
=
p
.
·
The propagation velocity of the waves of shearing is given by:
G
=
S
.
·
The modulus of rigidity
G
, the modulus of elasticity
E
(Young) and the coefficient of
Poisson
are connected by the relation:
(
)
+
=
1
2
E
G
.
5.1.3.2 Space discretization of the ground and Geometry of the stratification
The discretization in voluminal finite elements of the half space infinite ground is not accessible. It is
conventional for not limited fields of R
3
to have recourse to a formulation by equation
integral based on the knowledge of a fundamental solution which for (S) the medium (X) considered (S)
in ISS is called: functions of Green. This solution is then discretized by finite elements of
border what makes it possible to limit the space discretization of the field and thus to net only it (or
) the interface (S) (cf [§5.1.3.3]).
Recalls:
The functions of Green give the “impulse” response of the medium to a source
specific, in the absence of any reflective surface. The solution makes it possible to reproduce
the signal source with a shift corresponding to the time of way source-receiver, and
a decrease of the level proportional to the distance source-receiver.
Adapted, the functions of Green can take into account all or some
boundary conditions on obstacles. The use of these functions is particularly
useful for the integral formulations of the problems of radiation by the structures and
of diffraction by obstacles.
The finite elements of border built and used by PROMISS 3D are generated from
“interfacing” of the mesh of the foundation belonging to the external field (Code_Aster). It
mesh must be carried out to leave D `linear or quadratic surface elements to the normals
imperatively directed towards the interior of the ground. However, it is to be announced that the presence of elements
quadratic on the interface nothing brings. Indeed, for PROMISS 3D, fields being used for calculation
functions of Green are constant by side of element (the presence of node S intermediaries does not have
no interest for PROMISS 3D). Moreover, it should be noted that PROMISS 3D has a method of
calculation original which avoids netting the interface between the various layers of the ground.
Two cases of figure can arise:
·
The foundation is surface: In this case, it is enough to only one level source and receiver
located on the free face (at the higher level of the soil horizon in contact with the air).
Concretely, in the file nom_étude.sol one specifies that only one source is
necessary, thanks to the key word
SOURCE
, by indicating if one solves a problem in
geometry 2D or 3D. The key word
RECEP
, which announces the position of the receiver for calculation
functions of Green, must appear on the description of the 1čre layer.
stresses being null on the free face, it is useless to calculate them what leads to
to use the option
ALGO_DEPL
for the calculation algorithm of the functions of Green.
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
13/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
One gives below an example of file nom_étude.sol commented on:
TITR
GROUND ONGBABY
* Name which one gives to the characteristics ground used
MATERIAL 2
RO E
2150 4400.E06
2070 1421.E06
* A number of Materials associated with the stratification with the ground
NAKED BETA
ETA
0.45
0.08 0.
0.45
0.114 0.
* Description of the soil mechanics characteristics
2 SLEEP
* Numbers stratification taken into account during calculation
25.0 SUBDUE 1
RECEP
* Thickness and material associated with the soil horizon.
* receiving is placed at the top of the 1
era
sleep
25.0 2 SUBDUE
* Layer 2 is not in contact with the foundation. One does not place there
* of receiver
SUBS SUBDUE 2
* The substratum is located, in the example, under the 2
čme
sleep
SOURCE 1 3D
* Shallow foundation of structure => Only one source
FORCE HORIZ
POSI 1
* The source is applied to the node of the 1
era
sleep
ALGO DEPL
* Because shallow foundation
SPEC CAR
* Automatic management of the sampling of the functions of Green
OFFSET 110/440
* Parameter of horizontal sampling of the functions of Green
Layer 1
Layer 2
Interface:
Elements
surface
linear.
Normals
are directed
towards the interior
ground medium
Appear 5.1.3.2-a: Representation of the ground with shallow foundation
·
The foundation is buried: It previously was seen, the finite elements of border
built and used by PROMISS 3D are generated starting from the mesh of the foundation
belonging to the external field (Code_Aster). In the case of a mesh cutting it
volume of laminated space, it is appropriate to have several levels sources and receivers
to cover the whole of the mesh. Thus, on the side part of the buried foundation, one
place:
-
a source in the center of gravity of each element, as on the level of the base of
foundation (cf [Figure 5.1.3.2-b]),
- a receiver on points of Gauss of each element. The rule stated in
PROMISS 3D is to lay out with more the 6 receivers on each element but it is
recommended to use only 4 of them placed at the node of the element like all them
quarters length of the element (cf [Figure 5.1.3.2-b]).
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
14/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
The provision of the receivers and the sources in the file nom_étude.sol is generic it is appropriate
thus to use for the side part of the foundation in contact with the ground a regulated mesh where all
the elements have even height (cf [Figure 5.1.3.2-b]).
Concretely, in the file nom_étude.sol one specifies the number of sources necessary, thanks to
key word
SOURCE
, by indicating if one solves a problem in geometry 2D or 3D. The word is applied
key
RECEP
to announce the position of the receivers, necessary to the calculation of the functions of Green, on
ad hoc soil horizons. At the time of the taking into account of the burial of a foundation, PROMISS 3D
implements formulations regularized to improve the precision of calculations and to limit the effort
of integration. In the case of a buried foundation it is thus imperative to use the option
ALGO_REGU
for the calculation algorithm of the functions of Green.
One gives below an example of file nom_étude.sol commented on. One supposes, for the example
considered, that the foundation is buried of 20m and that one has two elements on the height of
foundation (cf [Figure 5.1.3.2-b]):
TITR
GROUND OTU *
Name which one gives to the characteristics ground used
MATERIAL 2 *
A number of materials associated with the stratification with the ground
RO
E
NAKED
BETA
ETA
2150
4400.E06
0.45 0.08 0.
*
Description of the characteristics
2070 1421.E06
0.45 0.114 0. *
soil mechanics
10 SLEEP *
Numbers stratification taken into account during calculation
2.5 1 RECEP SUBDUE *
The 1
era
lay down ground is divided into 9 underlayers.
2.5 1 RECEP SUBDUE *
A receiver is placed at the top of each underlayer.
2.5 1 RECEP SUBDUE *
the first 8 underlayers are in opposite with the foundation.
2.5 1 RECEP SUBDUE *
2.5 1 RECEP SUBDUE *
2.5 1 RECEP SUBDUE *
2.5 1 RECEP SUBDUE *
2.5 1 RECEP SUBDUE *
5.0 1 RECEP SUBDUE *
Remainder of 1
era
horizon soil located in lower part of the foundation
25.0 2 SUBDUE *
Layer 2 is not in contact with the foundation. One does not place there
*
of receiver
SUBS SUBDUE 2 *
The substratum is located, in the example, under the 2
čme
sleep
SOURCE 3 3D *
3 (2+1) Sources placed at the center of gravity of each element
*
(2) and on the level of the base of the foundation (1)
FORCE HORIZ POSI 3 *
The source is applied to the node of the 3
čme
underlayer
FORCE HORIZ POSI 7 *
The source is applied to the node of the 7
čme
underlayer
FORCE HORIZ POSI 9 *
The source is applied to the node of the 9
čme
underlayer
ALGO REGU *
Because buried foundation
SPEC CAR *
Automatic management of the sampling of the functions of Green
OFFSET 110/440 *
Parameter of horizontal sampling of the functions of Green
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
15/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
Layer 1
Layer 2
Receiver
Receiver and
Source
Interface:
Elements
surface
linear.
Normals
are directed
towards the interior
ground medium
Appear 5.1.3.2-b: Representation of the ground with buried foundation
Note:
If on part of the burial, there is no rigid connection between the wall of the ship
buried and the ground, then it is necessary to define the surface elements of this wall without connection
rigid ground-ship by the key word
GROUP_MA_SOL_SOL
of the operator
IMPR_MACR_ELEM
.
Moreover, it is necessary to direct these surface elements towards the outside of the ground medium.
5.1.3.3 Parameters of sampling of the functions of Green
The calculation of the functions of Green follows the following stages:
·
decomposition of the solution in plane or cylindrical waves elementary,
·
resolution of the problem of the elementary waves by the methods of the coefficients of reflection
transmission,
·
synthesis of the solution in Cartesian space (space field) by transform of Fourier
opposite.
The functions of Green are thus sampled
.
The parameter of horizontal sampling: OFFSET
OFFSET = Xmax/NR
With Xmax length wraps with the biggest length of the foundation. For the small island the length of
to erase is 110 m => X max = 110 Mr.
With NR, a number of points of sampling. It is given starting from the average size of
elements of the foundation. This length to erase it nuclear small island is approximately 3 Mr. One chooses
12 points of sampling per element.
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
16/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
The distance between 2 samples is then given by: D = 3.00/12 = 0.25m.
The number of samples is thus equal to: NR = 110/0.25 = 440.
The parameter of spectral sampling: SPEC
SPEC = CAR PROMISS 3D manages the parameters of sampling.
5.1.4 The parameter setting of postprocessing specific to MISS 3D
It is possible, except standard use of the interface Code_AsterPROMISS 3D, to use tools of
calculations specific to MISS 3D and to recover the useable results or not by Code_Aster
(cf [§4.3]).
5.1.4.1 Calculation MISS 3D with variable pitch of frequency
To increase the precision of a calculation of ISS, one can re-enter a list of frequencies to the hand.
It is in the file:nom_étude.optmis (cf [§ 5.1.2]) that one defines the number and the list of
frequencies on which will carry calculation.
The control MISS 3D for the list of frequency is: LFREQ to which one must associate a continuation of
values corresponding to the frequencies for which one will carry out calculation.
LFREQ = nf
with nf a number of frequencies retained for calculation
f1 f2 f3…. fi…. fn
with fi frequency for which one wishes to carry out calculation
One gives below an example of file nom_étude.optmis commented on. In the example
considered, the foundation surface, is excavated 11.60m and one looks at the answer on 6
precise frequencies:
*
* File nom_etude.21 (.optmiss)
*
*
LFREQ 6
* One carries out a calculation on 6 frequencies.
7.1.7.2.7.3 12.3 15.7 18.9
*
List frequencies of calculation
Z0
- 11.6
*
Foundation
excavated
with
11.60m
SURFING *
Foundation
surface
Note:
In the file nom_étude.optmis one must choose the option either Liste of Frequency LFREQ or
that is to say sampling of the Frequency: FREQ. Options LFREQ and FREQ are not
compatible.
5.1.4.2 Calculation MISS 3D on points of control.
The points of control make it possible to recover information on the incidental fields, diffracted
or radiated by (S) the interface (S), anywhere in the ground.
To carry out calculations on points of control, one indicates, in the file: nom_étude.optmis
(cf [§ 5.1.2]), their number (instruction CONT) like their geometrical co-ordinates (X, y, Z) in
ground.
The calculation of the fields starting from the points of control is post-treaty as well in time in
frequency by MISS 3D because it does not intervene in the resolution of the problem coupled between the different ones
under-fields. One must thus define, in the nom_étude.optmis (cf [§ 5.1.2]), a data file
auxiliaries (cf [§ 4.3.2.2]) (instruction FICP) which must preferably (that is strongly to advise)
to reside in the index associated with the study on the server dedicated to Code_Aster and the software
PROMISS 3D.
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
17/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
In this file of postprocessing one must, in the order:
1)
To stipulate if one wishes to carry out a calculation in the temporal field (instruction TIME)
rather than into frequential (default option).
Note:
It is necessary to specify the end of calculations in time (return in the frequential field) by
the FINTime instruction.
2)
To see (instruction LIRA) the file with the extension .sign which contains the transform of
Fourier of the signal used as excitation in far fields during coupled calculation
Code_Aster PROMISS 3D. This file is in the index associated with the study on
server dedicated to Code_Aster and software PROMISS 3D,
Note:
During postprocessing, for calculations in time or frequency, in order to
to recover a coherent temporal signal with the excitation, the FFT of the signal of excitation
in far fields must be filtered (instruction FILTERS). In the frequential field
MISS 3D multiplies the FFT of the signal of excitation by a window of ampitude 1 on
all the frequency band. At the time of the passage in the temporal field that
amounts carrying out the product of convolution according to
:
-
=
-
T
D
T
sign
signal
T
excitation
D
signal
FFT
)
(
*
)
.
(
)
(
1
who allows
to completely describe the signal of excitation in the temporal field.
3)
To give the name of a file temporal result (or frequential) which will contain them
displacements, speeds or accelerations calculated at the point of control (instruction
FICH).
4)
To define which type of field (incidental, radiated or diffracted) will be used for calculation.
table below points out the whole of the results which can be obtained in
postprocessing starting from point of control per MISS 3D:
Instruction MISS 3D
Results
Associated field *
CUI
Displacements, speeds, accelerations
Incidental field
UCTR
Displacements, speeds, accelerations
Radiated field
CTOT **
Displacements, speeds, accelerations
Diffracted field
CSOL **
Displacements, speeds, accelerations
Fields incident+diffracté
* Incidental field, radiated or diffracted by the interfaces of the field
** Attention: For fields diffracted CTOT and CSOL before defining thanks to instruction FICP
the file of post processing it is imperatively necessary to insert in the file nom_étude.optmis the word
key CHMI (Field with the Interface). One gives below an example of file nom_étude.optmis
commented on:
* File Nom_etude.21 (.optmis)
FREQ .01 35.01 .25
* Calculation between 0.1 Hz and 35.01Hz by pitch of 0.25 Hz
ZO 0.
* Not excavated foundation
SURFING
* Shallow foundation
CONT 2
* 2 points of control will be used for the post-
* processing
0. 100. 0.
* Co-ordinates (X, y, Z) of the 1
er
not control
0. 100. 5.4
* Co-ordinates (X, y, Z) of the 2
Nd
not control
CHMI
* Taking into account of the field diffracted by the interface
FICP/home/gubonva/uaster/Br/fichier.post
* File of postprocessing of the points of control
* described and commented on hereafter
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
18/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
One can refer to the example of postprocessing of points of control with accompanying notes page
following to obtain an example of file of post processing (fichier.post) by taking guard
however to replace control CUI (incidental field) by CTOT or CSOL.
Note:
By defect the instructions of the table above deliver results in
displacement. To obtain the results of speed or in acceleration it is necessary to use
the operand
CAPTION
follow-up of the key word
QUICKLY
or
ACCE
. For example: CUI CAPTIONS
ACCE
in this case MISS 3D calculates the response in acceleration to the point of control.
Caution:
It is imperative to be coherent on the type of calculation to realize compared to the file
.sign (cf .2) exploited for postprocessing. If this file comes, for example, of one
signal in acceleration it is appropriate to seek a result in acceleration.
5)
As MISS 3D works in the frequential field, one must specify, among all them
sampled frequencies, on which numbers of frequency one wishes to carry it out
postprocessing (instruction FREQ).
Note:
In the case of a postprocessing in time cf 1., it is advisable to select all
the frequencies (instructions FREQ ALL), if not a filtering will be
automatically applied.
6)
To specify starting from which type of fields of excitation one will carry out postprocessing
(instruction FIELD). One can thus use one, two or the three fields relating to the waves
of pressure and shearing in the ground of the incidental far field.
7)
To indicate on which degrees of freedom one wishes post-to treat (instruction DDL).
Note:
in 3D:
DDL 1 corresponds to direction X
DDL 2 corresponds to the direction y
DDL 3 corresponds to direction Z
8)
To define the check-point on which one carries out postprocessing (instruction NOT).
One gives below an example of file nom_étude.optmis commented on. For the example
considered, the foundation surface, is not excavated and one carries out postprocessings on two
points of control:
* File Nom_etude.21 (.optmis)
FREQ .01 35.01 .25
* Calculation between 0.1 Hz and 35.01Hz by pitch of 0.25 Hz
ZO 0.
* Not excavated foundation
SURFING
* Shallow foundation
CONT 2
* 2 points of control will be used for the post-
* processing
0. 100. 0.
* Co-ordinates (X, y, Z) of the 1
er
not decontrōle
0. 100. 5.4
* Co-ordinates (X, y, Z) of the 2
Nd
not decontrōle
FICP/home/gubonva/uaster/Br/fichier.post
* File of postprocessing of the points of control
* described and commented on hereafter
One gives below an example of file of postprocessing of points of control with accompanying notes.
For the example considered, one seeks, accelerations in the temporal field caused by
incidental fields in two points of control located at a distance of 100 m of the foundation and at
respective depths of 0.m and 5.40m (cf [Figure 5.1.4.2-a]).
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
19/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
Layer 1
Layer 2
y
X
Z
100m
Not control 1
Not control 2
5.40m
Appear 5.1.4.2-a: Example of points of control
* File/home/gubonva/uaster/Br/fichier.post
* The file must be placed on the Aster server.
TIME NT=1024 TMAX=10.24
* One carries out a return in time to carry it out
* calculation. NT=nombre of pitch of time.
* TMAX=longuor of the window of time.
FILTER LIRA RÉPONSE.01.SIGN
* The incidental FFT of the fields is recovered
* calculated by MISS 3D starting from accéléro.
FICH réponse.01.p1.t
* results obtained from the 1
er
not
* control will be recorded in the file
* réponse.01.p1.t
CUI CAPTIONS ACCE
* Calculations of accelerations associated with the field
* incidental.
FREQ ALL
* One calculates displacements from all
* sampled frequencies.
FIELD ALL
* horizontal components and vertical of
* incidental field will intervene during calculation of
* displacements.
DDL 1
* One calculates here only the component in X of
* displacement.
POINT 1
* Calculations previously described will be done on
* the 1
er
not control.
FICH réponse.01.p2.t
* results obtained from the 2
Nd
not
* control will be recorded in the file
* réponse.01.p2;T
CUI CAPTIONS ACCE
* Calculations of displacements associated with the field
* incidental.
FREQ ALL
* One calculates displacements from all
* sampled frequencies.
FIELD ALL
* horizontal components and vertical of
* incidental field will intervene during calculation of
* displacements.
DDL 1
* One calculates here only the component in X of
* displacement.
POINT 2
* Calculations previously described will be done
* to the 2
Nd
not control.
FINT
* End of calculations in time
EOF
* End of file
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
20/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
5.1.4.3 Calculation MISS 3D in foundation buried with suppression of fictitious resonances
The formulation of the elements of borders applied to the dynamic problem of the interaction ground
structure can reveal frequencies of fictitious resonances which correspond to
frequencies of excitation of the finished field hidden (foundation) in a semi-infinite field (ground). These
fictitious resonances appear in certain configurations: soft grounds, large foundation raft (cf [§ 4.3.1]).
To free itself some one uses control RFIC in the file: nom_étude.optmis (cf [§ 5.1.2]).
Calculation with the elimination of fictitious resonances can be is carried out directly at the time from the 1
er
calculation
MISS 3D is post-treaty by MISS 3D after having realized of the phenomenon in comparison with
results.
Note:
Control RFIC is very greedy in calculating times. To use can go until doubling
calculating times for MISS 3D. It is thus advised to use control RFIC rather in
postprocessing.
The syntax of control RFIC of MISS 3D is as follows:
RFIC telemfd with telemfd size in meter of the smallest element of the mesh of the foundation
One gives below an example of file nom_étude.optmis (cf [§ 5.1.2]) commented on. In
the example considered, the foundation is buried a depth of 6.05m. the size of smallest
element of the foundation is 3 Mr. One looks at the answer on 6 precise frequencies in
being been free from the problem of fictitious resonances:
*
* File nom_etude.21 (.optmiss)
*
*
LFREQ 6
* One carries out a calculation on 6 frequencies.
7.1.7.2.7.3 12.3 15.7 18.9
*
List frequencies of calculation
Z0
6.05
*
Foundation
buried
with
6.05m
RFIC 3
* Elimination of fictitious resonances
5.1.4.4 Calculation MISS 3D with simultaneous search of the impedances of ground
The impedances of ground express the dynamic rigidity of the field. They are expressed in the form of
square matrix depending on the frequency. Each line and each column of this matrix
corresponds to a particular mode, a term of the matrix being the virtual work exerted by one of these
modes on another mode.
Note:
During a calculation in time (instruction TIME) no impedance of ground can be calculated.
MISS 3D allows to calculate the impedances of ground (instruction IMPDC) at the same time as
response of the foundation to the seism. For this calculation of the impedances of ground, MISS 3D carries them out
same calculations as a postprocessing with the operand
MISS_IMPE
control
MACRO_MISS 3D
of Code_Aster.
The result of the calculation of the impedances (matrices) is stored in a file result (instruction IMPE)
who must imperatively reside in the index associated with the study on the server dedicated to
Code_Aster and with software PROMISS 3D.
One gives below an example of file nom_étude.optmis (cf [§ 5.1.2]) commented on. In
the example considered, the foundation surface is excavated a depth of 11.60m.
* File Nom_etude.21 (.optmis)
FREQ .01 20. 0.1
* Calculation between 0.1 Hz and 20Hz by pitch of 0.1 Hz
ZO 11.60
* Excavated foundation with 11.60m
SURFING
* Shallow foundation
IMPE
/home/gubonva/uaster/Br/réponse.01.impe
* the calculated impedances will be recorded in
* file réponse.01.impe
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
21/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
5.1.5 Output data of PROMISS 3D
Following calculation by PROMISS 3D, Code_Aster recovers the modal factors of participation for
to obtain displacement, speed and the acceleration in any point of the structure by recombination
modal. This recombination can S `write in the form:
()
() ()
() ()
+
=
With
B
B
B
With
With
X
T
X
T
T
X
U
,
with
()
T
X
U,
: field of displacement of the structure
B
With
,
: dynamic and static modes
B
With
,
: modal factors of participation
For calculations of the spectra of answer, the file result which contains the factors of participations
modal names nom_étude.nn.t (where N corresponds to the sequence number of the coming loading
of Code_Aster is the occurrence of the call to the control
IMPR_MISS_3D
). It is created
automatically under the central machine of processing in the index indicated in the control
MACRO_MISS_3D
.
For the calculation of the transfer functions one applies like loading a harmonic excitation of
modulate 1. That is carried out in Code_Aster thanks to the control
IMPR_MISS_3D
. The file in
left PROMISS 3D comprising the dynamic response structure complexes is named then
nom_étude.nn.h (where N corresponds to the sequence number of the loading coming from Code_Aster is
occurrence of the call to the control
IMPR_MISS_3D
). It is created automatically under the machine
of processing in the index indicated in the macro-control
MACRO_MISS_3D
.
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
22/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
5.2
Method of calculation on rigid foundation
The rigid condition of foundation is obtained by imposing on the foundation raft of the ship a movement of
solid body (the number of static modes associated is then tiny room to 6). For that one uses them
linear relations of type
LIAISON_SOLIDE
in the control
AFFE_CHAR_MECA
[U4.44.01] enters
the whole of the node S of the foundation (gathered in the group of node: SRADIER).
The whole of the node S of the foundation is then connected to the central node:
PO
. The unfolding of
calculation is done in the following way [Figure 5.2-a]:
Mesh small island
RECU_FONCTION
CALC_FONCTION
Order of Code_Aster
Calculation of the dynamic modes
structure in embedded base
(Dx=Dy=Dz=Drx=Dry=Drz=0)
Definition of the model
DEFI_BASE_MODALE
MACR_ELEM_DYNA
Impression of the data for
the calculation of ISS per MISS 3D
LIRE_MISS_3D
AFFE_CHAR_MECA
MACRO_MATR_ASSE
POST_ELEM
MACRO_MODE_MECA
LIRE_MAILLAGE
DEFI_MATERIAU
AFFE_MATERIAU
AFFE_MODELE
AFFE_CARA_ELEM
Calculation of the 6 static modes
structure in rigid foundation
(LIAISON_SOLIDE) and PO blocking
(Dx=Dy=Dz=Drx=Dry=Drz=0)
Definition bases modal
and
Projection of M, K, C on the basis
Launching of MISS 3D
IMPR_MACR_ELEM
IMPR_MISS_3D
AFFE_CHAR_MECA
MACRO_MATR_ASSE
MODE_STATIQUE
MACRO_MISS_3D
Restitution on the basis of physical
transitory answer resulting MISS 3D
Calculation of the spectra
and of the transfer functions
End
Be reproduced Synoptic 5.2-a: of calculation on rigid foundation
The associated command files are given in Appendix 2.
1
era
Stage
2
čme
Stage
3
čme
Stage
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
23/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
5.3 Method of calculation on flexible foundation with totality of
static modes
The interface between the foundation and the ground is in the case of the ship already studied represented by
group nodes
SRADIER
on which only the translations are locked (Dx=Dy=Dz=0) the number
static modes to calculate amounts to 1731 which corresponds to the number of node S on the foundation
(577 nodes) multiplied by the number of associated degrees of freedom.
The method implemented is as follows [Figure 5.3-a]:
Mesh small island
RECU_FONCTION
CALC_FONCTION
Order of Code_Aster
Calculation of the dynamic modes
structure in locked base
(Dx=Dy=Dz=0)
Definition of the model
Boundary condition
DEFI_BASE_MODALE
MACR_ELEM_DYNA
Impression of the data for
the calculation of ISS per MISS 3D
LIRE_MISS_3D
MACRO_MATR_ASSE
POST_ELEM
MACRO_MODE_MECA
LIRE_MAILLAGE
DEFI_MATERIAU
AFFE_MATERIAU
AFFE_MODELE
AFFE_CARA_ELEM
AFFE_CHAR_MECA
Calculation des1731modes static
structure in flexible foundation
and bases locked
(Dx=Dy=Dz=0)
Definition bases modal
and
Projection of M, K, C on the basis
Launching of MISS 3D
IMPR_MACR_ELEM
IMPR_MISS_3D
MODE_STATIQUE
MACRO_MISS_3D
Restitution on the basis of physical
transitory answer resulting MISS 3D
Calculation of the spectra
and of the transfer functions
End
Be reproduced Synoptic 5.3-a: of calculation on flexible foundation with the totality of the static modes
The command files are given in Appendix 2.
1
era
Stage
2
čme
Stage
3
čme
Stage
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
24/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
5.4
Method of calculation on flexible foundation with some modes of
foundation selected
5.4.1 Main stages of the method
For the calculation of ISS, PROMISS 3D needs a base made up at the same time of null clean modes on
the interface ground-structure and other nonnull modes on this interface. This base corresponds to one
interface of type CRAIG_BAMPTON. For the first modes, one generally takes the modes
clean of the structure obtained by locking displacements on the interface (bases locked) and for
the seconds, one takes static the modes known as “constrained”, successively obtained while imposing
a unit displacement of each degree of freedom of each node of the interface (
foundation).
The principle of the method implemented here consists in replacing the constrained static modes
plethoric by clean modes of foundation in small number calculated on carpet of springs of ground
and selected according to an established criterion.
Several stages are necessary to conclude calculation:
5.4.2 Determination of the carpet of springs to be placed under the foundation.
The values of the stiffnesses of the springs equivalent on the laminated ground of the ship are given with
through a calculation of the transfer functions under harmonic stress of module 1 at the time of the study
with rigid foundation (cf [§5.2]).
One obtains the 6 values of total stiffness within the competence of the laminated ground: Kx (NR/m), Ky (NR/m), Kz (NR/m),
K
X (N.m), K
y (N.m), K
Z (N.m). These stiffnesses, independent of the frequency, are distributed with
proportion of surfaces of the elements around the nodes of the foundation thanks to the operand
RIGI_PARASOL
control
AFFE_CARA_ELEM
[U4.42.01] of Code_Aster.
5.4.3 The calculation of the dynamic modes of the structure
This calculation is carried out on basis embedded with the control
MODE_ITER_SIMULT
(one applies to
all nodes of the foundation the following boundary condition: Dx=Dy=Dz=Drx=Dry=Drz=0).
5.4.4 The calculation of the clean modes of foundation on carpet of spring
During calculation, one dissociates the modes with nonnull displacements of the infrastructure (to erase)
modes of the superstructure (ships…) by considering that only the foundation raft is heavy. This is
realized while applying, with the elements not modelizing the foundation, a material of which mass
voluminal is null. One avoids thus, during the construction of the modal base gathering the modes of
foundation and dynamic of the structure, to consider the clean modes twice of
superstructure.
One enriches then the modal base established with [§5.4.3], via the control
DEFI_BASE_MODALE
, by
first calculated modes which all are of the modes of foundation since are rejected towards the high ones
frequencies all modes of the superstructure.
5.4.5 Selection of the modes
While reducing considerably the number of constrained modes of foundation one can manage to find
the solution in answer and frequency of resonance obtained with the preceding method putting in
work the totality of the static modes (cf [§5.3]) and allowing a saving of time of substantial calculation.
One judges that the method of reduction is interesting, in term of saving of time, when the number
clean modes of foundation on carpet of spring is with most equal to the third of the number of modes
statics on flexible foundation (for this study, the method is interesting if the number of modes
of foundation on carpet of spring is lower than 1731/3
500 modes cf [§5.3]).
To refine the selection of the modes, one can use the method recommended by E. Balmes [bib2] which
consist in retaining only the modes of foundation whose Eigen frequency remains lower than twice
the cut-off frequency used during the calculation of the dynamic modes [Figure 5.4.5-a].
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
25/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
Cut-off frequency used
during dynamic calculation
For the ship fc=23 Hz
Classification of the ddl of the structure
total starting from the classification of
the modal base of Ritz Q
Construction of the modal base of Ritz Q
containing the dynamic modes and them
modes of foundation with:
N
= a number of dynamic modes
N
= a number of modes of foundation
(N
must be with most equal to 500 for the ship)
N= N
+n
= an optimal number of modes
DEFI_BASE_MODALE
NUME_DDL_GENE
PROJ_MATR_BASE
MODE_ITER_SIMULT
REST_BASE_PHYS
Determination of the optimal number
modes of foundation
N
= N
such as NR counts all the modes
whose frequency is at the maximum
equalize with 2xfc (for the ship 46 Hz)
Projection of the matrices of mass
and of rigidity associated with calculation with
dynamic modes on the modal basis Q
Calculation of the modes generalized with the matrices
of mass and rigidity projected and restitution
modes orthogonalized on the physical basis
Appear 5.4.5-a: Optimization of the number of modes of foundation
The course of complete calculation with reduction of the modes of foundation is made way
following [Figure 5.4.5-b]:
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
26/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
Warnings:
* For the modes of foundation, contrary to the calculation of the dynamic modes, the nodes of
to erase do not have any stress on their displacements. To calculate the macronutrient
dynamics of the structure it is thus necessary to adapt two classifications related to limiting conditions
different (blocking with the interface for the dynamic modes or carpet from specific springs for
modes of foundation).
** One avoids, during the construction of the modal base (gathering the modes of foundation and
dynamic of the structure), to twice take into account the clean modes of main
under structures: internal chambers, structures,… by considering that only the foundation raft is heavy and
while applying to the elements modelizing the superstructure of materials with density
quasi null (
=10
- 3
kg/m
3
).
Mesh small island
RECU_FONCTION
CALC_FONCTION
Order of Code_Aster
Calculation of the dynamic modes
structure in embedded base
(Dx=Dy=Dz=Drx=Dry=Drz=0)
Definition of the model
Boundary condition
DEFI_BASE_MODALE
(NUME_REF:num_dyn)
MACR_ELEM_DYNA
Impression of the data for
the calculation of ISS per MISS 3D
LIRE_MISS_3D
Calculation of the spectra
transfer functions
LIRE_MAILLAGE
DEFI_MATERIAU
AFFE_MATERIAU
AFFE_MODELE
AFFE_CARA_ELEM
AFFE_CHAR_MECA
Calculation of the modes of foundations
on carpet of spring while considering
that only the foundation raft is heavy.
No Boundary conditions on
the foundation
Definition bases modal
and
Projection of M, K, C on the basis
Launching of MISS 3D
IMPR_MACR_ELEM
IMPR_MISS_3D
MACRO_MISS_3D
Restitution on the basis of physical
transitory answer resulting MISS 3D
End
MACRO_MATR_ASSE
(NUME_DDl:num_dyn
CHAMP_MATER:mat_dyn)
POST_ELEM
MACRO_MODE_MECA
AFFE_MATERIAU
AFFE_CARA_ELEM
(RIGI_PARASOL)
MACRO_MATR_ASSE
(NUME_DDl:num_fon
CHAMP_MATER:mat_fon)
POST_ELEM
MACRO_MODE_MECA
Be reproduced Synoptic 5.4.5-b: of calculation with modes of foundation on carpet of spring
*
**
*
**
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
27/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
One can carry out a parametric study on the number of modes of foundation. Method put in
work was evaluated [bib4] and consisted in making a continuation on a basis containing the 171 modes
dynamic and 500 modes of foundation and to retain at the time of the definition of the modal base only 80,
then 30 and finally 20 modes of foundation. This operation was carried out by modifying in continuation it
a number of modes of foundation by the operand
NMAX_MODE
control
DEFI_BASE_MODALE
,
then by starting again the channel of call to the successive operators starting from this operator until
postprocessing.
It should be noted that, on this study, when one uses the method of reduction of the modes, one finds that only
forty modes of foundation are necessary to reproduce the effect induced by
1731 static modes.
6 Bibliography
[1]
D. CLOUTEAU: User's manual of PROMISS 3D - MISS 2D, overhaul 6.3, by (LMSSM
Central school of Paris)
[2]
E. BALMES: Use generalized interfaces off dismantle off freedom in component synthesis mode
IMAC 1996
[3]
V. GUYONVARH - G. DEVESA: Methods of calculation of the seismic excitations to the works
CP N4. HP-52/99/006/A
[4]
V. GUYONVARH - G. DEVESA: Methods to consider the interaction ground-structure on the small island
nuclear power EPR with Code_Aster and MISS 3D. HP-62/00/007/A
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
28/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
Appendix 1 files of setting in data of PROMISS 3D
*
* File miss01.21 (.optmiss)
*
*
FREQ 0.1 20. 0.1
Z0 - 11.6
SURFING
*
* File miss01a.22 (.sol)
*
TITR
GROUND PENLY
MATERIAL 4
RO E NAKED ETA BETA
2150 4480.E06 0.40 0.08 0.
2070 1421.E06 0.45 0.114 0.
2150 1305.E06 0.45 0.16 0.
2400 6000.E06 0.45 0.06 0.
3 SLEEP
43.9 1 RECEP SUBDUE
31 2 SUBDUE
38.5 3 SUBDUE
SUBS SUBDUE 4
SOURCE 1 3D
FORCE HORIZ POSI 1
ALGO DEPL
* SPEC 0.12/16384
SPEC CAR
OFFSET 110/400
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
29/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
Appendix 2 Command files Aster
Initial common Left A2.1
# **********************************************************************
# Command file = miss01.comm
# **********************************************************************
BEGINNING ();
#
mall = LIRE_MAILLAGE ();
#
#
# -------------------------------------------------------------
# definitions of the model and the group of meshs
# -------------------------------------------------------------
#
#
INCLUDE (UNIT = 11);
#
# ------------------------------------------------------------
# definition of materials
# ------------------------------------------------------------
#
INCLUDE (UNIT = 12);
#
# -----------------------------------------------------------
# definition of the characteristics of the elements
# -----------------------------------------------------------
#
INCLUDE (UNIT = 13);
#
# -----------------------------------------------------------
# definition of the foundation
# -----------------------------------------------------------
&mail = DEFI_GROUP (
MESH = mall,
CREA_GROUP_NO =_F (
GROUP_MA = ' SRADIER',
),
);
A2.2 rigid Case Foundation
#
# ************************************************************
# CONDITION OF RIGIDITY OF THE FOUNDATION RAFT
# BOUNDARY CONDITION BASES BLOQUEE IN DYNAMICS
#
ch_cldyn = AFFE_CHAR_MECA (
MODEL = model,
DDL_IMPO =_F (GROUP_NO = ' SRADIER',
DX = 0.,
DY = 0.,
DZ = 0.,
DRX = 0.,
DRY = 0.,
DRZ = 0.,
),
);
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
30/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
# ************************************************************
# RESOLUTION OF THE DYNAMIC MODAL PROBLEM
#
# CALCULATIONS ELEMENTARY of the matrices of mass and stiffness
# CLASSIFICATION of the unknown factors of the problem
# ASSEMBLY of the matrices of mass and stiffness
#
#
MACRO_MATR_ASSE (model MODELE=,
CHAM_MATER= MATER,
CARA_ELEM= elem,
CHARGE= ch_cldyn,
NUME_DDL= num_dyn,
MATR_ASSE=_F (MATRICE= matrigi,
OPTION= “RIGI_MECA”),
MATR_ASSE=_F (MATRICE= matmass,
OPTION= “MASS_MECA”),
MATR_ASSE=_F (MATRICE= matamor,
OPTION= “AMOR_MECA”),
);
#
#-----------------------------------------
# calculation of the masses
#-----------------------------------------
#
#
#
masses = POST_ELEM (
MODEL = model,
INFORMATION = 1
CHAM_MATER = to subdue,
CARA_ELEM = elem,
MASS_INER =_F (
ALL = “YES”
),
);
#
#
# ----------------------------------------------------------
# calculation of the clean modes by successive tapes
# ----------------------------------------------------------
#
#
mod_dyn = MACRO_MODE_MECA (MATR_A= matrigi, MATR_B= matmass,
CALC_FREQ =_F (
FREQ = (0.1, 7., 10., 12., 14. , 16., 17., 19.,
21., 23.,),
),
NORM_MODE=_F (MASS_INER=masses),
FILTRE_MODE=_F (SEUIL= 1.D-3),
IMPRESSION=_F (),
);
#
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
31/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
# *****************************************************************
# CONDITION OF RIGIDITY OF THE FOUNDATION RAFT
# BOUNDARY CONDITION BASES BLOQUEE IN STATICS
# *****************************************************************
ch_clsta = AFFE_CHAR_MECA (
MODEL = model,
DDL_IMPO =_F (GROUP_NO = ' PO',
DX
=
0.,
DY = 0.,
DZ = 0.,
DRX = 0.,
DRY = 0.,
DRZ = 0.,
),
LIAISON_SOLIDE =_F (GROUP_NO = ' SRADIER'),
);
# *****************************************************************
# RESOLUTION OF THE STATIC MODAL PROBLEM
#
# CALCULATIONS ELEMENTARY of the matrices of mass and stiffness
# CLASSIFICATION of the unknown factors of the problem
# ASSEMBLY of the matrices of mass and stiffness
#
MACRO_MATR_ASSE (model MODELE=,
SOLVEUR=_F (METHODE= “MULT_FRONT”),
CARA_ELEM= elem,
CHARGE= ch_clsta,
CHAM_MATER= MATER,
NUME_DDL= num_sta,
MATR_ASSE=_F (MATRICE= rigistat, OPTION= “RIGI_MECA”),
(MATRICE= massetat, OPTION= “MASS_MECA”),
);
#
#
# CALCULATION OF THE DYNAMIC MACRONUTRIENT =
#---------------------------------------------
#
mod_sta = MODE_STATIQUE (MATR_RIGI= rigistat,
MATR_MASS= massetat,
DDL_IMPO=_F (GROUP_NO= “PO”,
TOUT_CMP= “YES”),
);
basmo = DEFI_BASE_MODALE (
RITZ=_F (MODE_MECA= mod_dyn),
RITZ=_F (MODE_STAT= mod_sta,
NMAX_MODE= 6),
NUME_REF= num_dyn);
#
mael = MACR_ELEM_DYNA (BASE_MODALE= basmo,
MATR_RIGI= matrigi,
MATR_MASS= matmass,
OPTION= “RITZ”);
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
32/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
A2.3 flexible Case Foundation with all the static modes
#
# ************************************************************
# CONDITION OF RIGIDITY OF THE FOUNDATION RAFT
# BOUNDARY CONDITION BASES BLOQUEE IN DYNAMICS
#
ch_cldyn = AFFE_CHAR_MECA (
MODEL = model,
DDL_IMPO =_F (GROUP_NO = ' SRADIER',
DX = 0.,
DY = 0.,
DZ = 0.,
),
);
# ************************************************************
# RESOLUTION OF THE DYNAMIC MODAL PROBLEM
#
# CALCULATIONS ELEMENTARY of the matrices of mass and stiffness
# CLASSIFICATION of the unknown factors of the problem
# ASSEMBLY of the matrices of mass and stiffness
#
#
MACRO_MATR_ASSE (model MODELE=,
CHAM_MATER= MATER,
CARA_ELEM= elem,
CHARGE= ch_cldyn,
NUME_DDL= num_dyn,
MATR_ASSE=_F (MATRICE= matrigi,
OPTION= “RIGI_MECA”),
MATR_ASSE=_F (MATRICE= matmass,
OPTION= “MASS_MECA”),
MATR_ASSE=_F (MATRICE= matamor,
OPTION= “AMOR_MECA”),
);
#
#-----------------------------------------
# calculation of the masses
#-----------------------------------------
#
#
#
masses = POST_ELEM (
MODEL = model,
INFORMATION = 1,
CHAM_MATER = to subdue,
CARA_ELEM = elem,
MASS_INER =_F (
ALL = “YES”
),
);
#
#
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
33/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
# ----------------------------------------------------------
# calculation of the clean modes by successive tapes
# ----------------------------------------------------------
#
#
mod_dyn = MACRO_MODE_MECA (MATR_A= matrigi, MATR_B= matmass,
CALC_FREQ =_F (
FREQ = (0.1, 7., 10., 12., 14., 16., 17., 19.,
21., 23.),
),
NORM_MODE=_F (MASS_INER=masses),
FILTRE_MODE=_F (SEUIL= 1.D-3),
IMPRESSION=_F (),
);
#
# CALCULATION OF THE DYNAMIC MACRONUTRIENT =
#---------------------------------------------
#
#
mod_sta = MODE_STATIQUE (MATR_RIGI= matrigi,
MATR_MASS= matmass,
DDL_IMPO=_F (GROUP_NO= “SRADIER”,
AVEC_CMP= (“DX”, “DY”, “DZ”),
),
);
basmo = DEFI_BASE_MODALE (
RITZ=_F (MODE_MECA= mod_dyn),
# 577 nodes has the interface
# => 577 * 3=1731 static modes
RITZ =_F (NMAX_MODE= 1800,
MODE_STAT=mod_sta),
NUME_REF= num_dyn);
#
mael = MACR_ELEM_DYNA (BASE_MODALE= basmo,
MATR_RIGI= matrigi,
MATR_MASS= matmass,
OPTION= “RITZ”);
A2.4 flexible Case Foundation with reduction of modes
#
# ************************************************************
# CONDITION OF RIGIDITY OF THE FOUNDATION RAFT
# BOUNDARY CONDITION BASES BLOQUEE IN DYNAMICS
#
ch_cldyn = AFFE_CHAR_MECA (
MODEL = model,
DDL_IMPO =_F (GROUP_NO = ' SRADIER',
DX = 0.,
DY = 0.,
DZ = 0.,
DRX = 0.,
DRY = 0.,
DRZ = 0.,
),
);
#
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
34/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
# ************************************************************
# RESOLUTION OF THE DYNAMIC MODAL PROBLEM
#
# CALCULATIONS ELEMENTARY of the matrices of mass and stiffness
# CLASSIFICATION of the unknown factors of the problem
# ASSEMBLY of the matrices of mass and stiffness
#
#
MACRO_MATR_ASSE (model MODELE=,
CHAM_MATER= MATER,
CARA_ELEM= elem,
CHARGE= ch_cldyn,
NUME_DDL= num_dyn,
MATR_ASSE=_F (MATRICE= matrigi,
OPTION= “RIGI_MECA”),
MATR_ASSE=_F (MATRICE= matmass,
OPTION= “MASS_MECA”),
MATR_ASSE=_F (MATRICE= matamor,
OPTION= “AMOR_MECA”),
);
#
#-----------------------------------------
# calculation of the masses
#-----------------------------------------
#
#
#
masses = POST_ELEM (
MODEL = model,
INFORMATION = 1,
CHAM_MATER = to subdue,
CARA_ELEM = elem,
MASS_INER =_F (
ALL = “YES”
),
);
#
#
# ----------------------------------------------------------
# calculation of the clean modes by successive tapes
# ----------------------------------------------------------
#
#
mod_dyn = MACRO_MODE_MECA (MATR_A= matrigi, MATR_B= matmass,
CALC_FREQ =_F (
FREQ = (0.1, 7., 10., 12., 14., 16., 17., 19.,
21., 23.),
),
NORM_MODE=_F (MASS_INER=masses),
FILTRE_MODE=_F (SEUIL= 1.D-3),
IMPRESSION=_F (),
);
#
#
#
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
35/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #
# DEFINITION OF CHARACTERISTICS OF GROUND TAKEN INTO ACCOUNT
# BY A CARPET OF SPRINGS
# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #
# GROUND HOMOGENEOUS PENLY
# ==================
#
elemb = AFFE_CARA_ELEM (….(idem other cases),…
RIGI_PARASOL=_F (
GROUP_MA= “SRADIER”,
COEF_GROUP= 1.,
GROUP_NO_CENTER= “PO”,
CARA= “K_TR_D_N',
VALE= (5.4E11, 5.4E11, 6.0E11,
6.5E14, 6.5E14, 1.09E15),
),
);
#
#
# *****************************************************************
# RESOLUTION OF THE MODAL PROBLEM OF FOUNDATION
# =========================================
#
# CALCULATIONS ELEMENTARY of the matrices of mass and stiffness
# CLASSIFICATION of the unknown factors of the problem
# ASSEMBLY of the matrices of mass and stiffness
#
#
MACRO_MATR_ASSE (model MODELE=,
SOLVEUR=_F (METHODE= “MULT_FRONT”),
CARA_ELEM= elemb,
# CHARGE= ch_clsta,
CHAM_MATER= materb,
NUME_DDL= num_stab,
MATR_ASSE=_F (MATRICE= rigistat, OPTION= “RIGI_MECA”),
(MATRICE= massetat, OPTION= “MASS_MECA”),
);
#
#
#-------------------------------------------------
# calculation of the masses reduced to the foundation
#-------------------------------------------------
#
#
#
masseb = POST_ELEM (
MODEL = model,
INFORMATION = 1,
CHAM_MATER = materb,
CARA_ELEM = elemb,
MASS_INER =_F (
ALL = “YES”
),
);
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
36/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
#
#
# ----------------------------------------------------------
# calculation of the clean modes by successive tapes
# ----------------------------------------------------------
#
# The calculation of the modes of “foundation” is carried out from
# The structure on which one keeps only the mass of the foundation raft
# because one wants that the modes with displacement not no one of the foundation
# and one do not want to recover second once the local modes of
# structure (mode of swinging of the chambers, of IF,….).
#
# The method of reduction of the modes is interesting
# when the number of modes of “foundation” is with most equal to the third
# of the number of static modes of reference. For the small island, one has
# 577 nodes on the interface is 577 * 3=1731 static modes.
# Using IMPR_STURMs one established the frequency band which
# enables us to keep 500 modes.
#
mod_sta = MACRO_MODE_MECA (MATR_A= rigistat, MATR_B= massetat,
CALC_FREQ =_F (
FREQ = (0.1, 60., 100., 130., 160., 200., 300.,
2000., 4000., 6000.),
),
NORM_MODE=_F (MASS_INER=masseb),
# FILTRE_MODE=_F (SEUIL= 1.D-3),
IMPRESSION=_F (),
);
#
# One defines our modal base while combining during our 1st test
# 171 dynamic modes of the structure in base encastree and 80
# “static” modes (foundation) with foundation on carpet of spring.
#
basmo = DEFI_BASE_MODALE (
RITZ=_F (MODE_MECA= mod_dyn),
RITZ=_F (MODE_MECA= mod_sta,
NMAX_MODE= 80),
NUME_REF= num_dyn);
mael = MACR_ELEM_DYNA (BASE_MODALE= basmo,
MATR_RIGI= matrigi,
MATR_MASS= matmass,
OPTION= “RITZ”);
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
37/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
Common Left A2.5: writing of the file .raster and launching of
PROMISS 3D
#--------------------------------------------------------
# READING OF the ACCELEROGRAMS (Penly),
#--------------------------------------------------------
#
#-----------------------------
# DEFI_FONCTION = UL 31
# Accelerogram = acc1.c2
# (resulting from the spectrum of ground Penly),
#-----------------------------
# DEFI_FONCTION = UL 32
# Accelerogram = acc2.c2
# (resulting from the spectrum of ground Penly),
#-----------------------------
# DEFI_FONCTION = UL 33
# Accelerogram = acc3.c2
# (resulting from the spectrum of ground Penly),
#
INCLUDE (UNIT = 31);
#
acce_x = CALC_FONCTION (
COMB =_F (
FUNCTION = acce1,
COEFF = 9.81,
),
);
#
INCLUDE (UNIT = 32);
#
acce_y = CALC_FONCTION (
COMB =_F (
FUNCTION = acce2,
COEFF = 9.81,
),
);
#
INCLUDE (UNIT = 33);
#
acce_z = CALC_FONCTION (
COMB =_F (
FUNCTION = acce3,
COEFF = 6.54,
),
);
#
# PREPROCESSING OF PROMISS 3D=
#--------------------------
# MESH AND IMPEDANCES Of INTERFACE
#
IMPR_MACR_ELEM (
MACR_ELEM_DYNA = mael,
FORMAT = “MISS_3D”,
SOUS_TITER = “miss01”,
IMPR_MODE_STAT= “YES”,
IMPR_MODE_MECA= “NOT”,
AMOR_REDUIT =_F (
#
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
5.00000E-02, 5.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
38/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 5.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 5.00000E-02, 5.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 5.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
5.00000E-02, 5.00000E-02, 7.00000E-02, 7.00000E-02, 7.00000E-02,
7.00000E-02, 7.00000E-02, 7.00000E-02, 5.00000E-02, 7.00000E-02,
7.00000E-02,
),
GROUP_MA_INTERF = “SRADIER”,
);
#
# EXCITATION HARMONIC OF MODULE 1
#
FO1 = DEFI_FONCTION (NOM_PARA= “FREQ”
VALE= (0., 1., 100., 1.) );
#
#
#
#
IMPR_MISS_3D (MACR_ELEM_DYNA= mael,
FREQ_INIT= 0., FREQ_FIN= 20. PAS= 0.1,
EXCIT_SOL=_F (DIRECTION= (1., 0., 0.), NOM_CHAM= “ACCE”,
FONC_SIGNAL= fo1),
EXCIT_SOL=_F (DIRECTION= (1., 0., 0.), NOM_CHAM= “ACCE”,
FONC_SIGNAL= fo1),
EXCIT_SOL=_F (DIRECTION= (1., 0., 0.), NOM_CHAM= “ACCE”,
FONC_SIGNAL= fo1),
);
#
#
#
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
39/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
# TRANSMISSION OF GIVE ASTER A PROMISS 3D
# CREATION OF the FILE = nom_etude.RASTER ON UNIT 26
IMPR_MISS_3D (
MACR_ELEM_DYNA = mael,
INST_INIT = 0.,
INST_FIN = 20.,
NOT = 0.01,
EXCIT_SOL =_F (
DIRECTION = (1., 0., 0.),
NOM_CHAM = “ACCE”,
FONC_SIGNAL = acce_x,
),
EXCIT_SOL =_F (
DIRECTION = (1., 0., 0.),
NOM_CHAM = “ACCE”,
FONC_SIGNAL = acce_y,
),
EXCIT_SOL =_F (
DIRECTION = (1., 0., 0.),
NOM_CHAM = “ACCE”,
FONC_SIGNAL = acce_z,
),
);
END ();
#
# PROCESSING OF PROMISS 3D BY EXEC_LOGICIEL
#---------------------------------------
#
# Response transitory of the structure
# subjected to the seismic loading
#
CONTINUATION ();
MACRO_MISS_3D (
OPTION =_F (TOUT= “YES”),
PROJET= “miss01”,
REPERTOIRE= “./uaster/miss01/”,
UNITE_IMPR_ASTER= 26,
UNITE_OPTI_MISS= 21,
UNITE_MODELE_SOL= 22,
);
END ();
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
40/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
Common Left A2.6: postprocessing of PROMISS 3D and examination
#
# ==================================================================
# Programs = miss01a_FDT.comm
# CALCULATION OF THE TRANSFER FUNCTIONS
# ==================================================================
#
CONTINUATION ();
#
# ==================================================================
# CALCULATION OF THE TRANSFER FUNCTIONS
# ON the NUCLEAR SMALL ISLAND
# One gives here as example that postprocessing on the Structures
interns with 1.50m
# ==================================================================
#
#
dyna = LIRE_MISS_3D (MACR_ELEM_DYNA= mael,
TYPE_RESU= “HARMO”,
TITER= “HARM_ACCE_EPR”,
UNITE=28);
#
#
#-----------------------------------
# STRUCTURE INTERNAL
#-----------------------------------
#
# Z=1.50 m
# =========
#
Hsi1x = RECU_FONCTION (
RESULTAT= dyna,
GROUP_NO= “NSIEZ3”,
NOM_CHAM= “ACCE”,
NOM_CMP= “DX”);
#
#
IMPR_COURBE (FORMAT = “AGRAF”,
FILE = “AGRAF”,
EXIT = “COLOR”,
TITER_GRAPHIQUE= “Functions IF 1.50m' has,
ECHELLE_X = “FLAX”,
ECHELLE_Y = “FLAX”,
LABEL_X = “frequency (Hz)”,
PRESENTATION = ' PAYSAGE',
DATE = ' OUI',
CURVE =_F (
FUNCTION = Hsi1x),
CURVE =_F (
PARTIE= “IMAG”,
FUNCTION = Hsi1x),
);
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
41/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
# ==================================================================
# Programs = miss01a_trans.comm
# CALCULATION OF THE SPECTRA OF ANSWER
# ==================================================================
#
CONTINUATION ();
#
#-------------------------------------------------------
#
CONTINUATION AFTER CALCULATION MISS
# ------------------------------------------------------
#
# =================================================================
# CALCULATION OF THE TRANSITORY ANSWERS
# ON the NUCLEAR SMALL ISLAND
# One gives as example calculations on the internal structures
# ==================================================================
#
resugene = LIRE_MISS_3D (MACR_ELEM_DYNA = mael,
TYPE_RESU = “TRANS”,
TITRATE = “TRANSIT”,
UNITE=28);
#
l_freq = (
0.200, 0.350, 0.500, 0.650, 0.950,
1.100, 1.250, 1.400 1.550, 1.700,
1.850, 2.000, 2.150, 2.300, 2.450, 2.600,
2.750, 2.900, 3.075, 3.300, 3.525, 3.800,
4.100, 4.400, 4.700, 5.000, 5.375,
5.750, 6.125, 6.500, 6.875,
7.250, 7.625, 8.000, 8.750, 9.500, 10.250,
11.000, 11.750, 12.500, 13.250, 14.000, 14.750,
16.000, 17.500, 20.000, 23.500, 28.000, 32.500,
37.000, 41.500, 46.000, 50.500, 56.000, 62.000,
74.000, 80.000, 86.000, 92.000, 98.000);
l_amor_s = (0.04);
#---------- Internal structures A 1.50 m ----------
#
SIAZdXr = RECU_FONCTION (
RESULT = resugene,
NOM_CHAM = “ACCE”,
TITER= “ABSOLUTE ACCELERATION IF Z=1.50 m O DEGR EXT. IN X”,
GROUP_NO = “NSIAZ3”
NOM_CMP = “DX”, INTERPOL = “FLAX”,
);
# ================================================================
# CALCULATION OF THE SPECTRA
# ================================================================
SIAZdXs = CALC_FONCTION (
SPEC_OSCI =_F (
FUNCTION = SIAZdXr,
FREQ = l_freq,
AMOR_REDUIT= l_amor_s));
Code_Aster
®
Version
6.4
Titrate:
Interaction ground-structure with the interface Code_Aster-PROMISS 3D
Date
:
09/06/04
Author (S):
G. DEVESA, V. GUYONVARH
Key
:
U2.06.07-B
Page
:
42/42
Instruction manual
U2.06 booklet: Dynamics HT-66/04/004/A
Intentionally white left page.