SYSTUS RefManAnalysisVol3 en PDF
SYSTUS RefManAnalysisVol3 en PDF
SYSTUS RefManAnalysisVol3 en PDF
2019
(VOL.3)
www.esi-group.com
SYSTUS 2019
REFERENCE ANALYSIS MANUAL
(VOL.3)
© 2019 by ESI Group
Published: January 2019
GL/SYST/18/08/03/A
This documentation is the confidential and proprietary product of ESI Group.
Any unauthorized use, reproduction, or transfer of this documentation is strictly prohibited.
This documentation is subject to limited distribution and restricted disclosure.
All PAM- and SYS- product names as well as other products belonging to ESI's portfolio are
tradenames or trademarks of ESI Group, unless specifically mentioned. All other trademarks are the
property of their respective owners.
Specifications are subject to change without prior notice.
Contents
INDEX i
elementary matrix files with a DATA file, for subsequent utilization with the DYNAMIC or
TRANSIENT routine.
General syntax:
SOLVE (labels)
l Computation check
o MATRICES
o RESULTS
o PARTIAL
o STIFFNESS
o MASS
o DAMPING or BUCKLING
the contribution of each element or node to the left-hand side (stiffnesses Ki) and right-hand
side (forces Fi).
b) Execution of routines which:
l assemble the basic stiffness matrices and load vectors Ki and Fi into the global stiffness
containing vectors, the components of which are infinitely small for the free
degrees of freedom (residues) and reactions for the constrained degrees of freedom.
The largest residue should be negligibly small compared with the resultant of
applied loads, the order of magnitude of residues compared with that of applied
loads being a measure of the precision of simultaneous equation solution.
c) Execution of routines which, according to the type of structure, produce force, stress and
strain files.
The data on results files is stored by load case.
SOLVE
The algorithm uses an internal numbering system. This minimizes the number of terms of the
matrices. The working files are stored automatically in core memory or on disk, and this
algorithm is compatible with all types of computer.
The triangulated matrix can be stored for use with other load cases.
If the triangulated matrix is not stored, the optional label WITHOUT is used to optimize the
space on the hard disk.
For more information about this command and the various solvers available see Chapter 21.
SEARCH DATA 1
SOLVE
SAVE MATRIX RESULTS 1
SEARCH DATA 2
SEARCH MATRIX 1
SOLVE STIFFNESS RESULTS
SAVE RESULTS 2
Forces in the elements are calculated from a displacement file. No other operation is required
apart from the calculation of stresses and strains at the Gaussian points, where these exist. If the
EXTRACT label is present, displacements and forces are printed.
The MODAL and TRANSIENT labels (see 11.2.1.3.1 and 12.3.5) are used to integrate
dynamic loads in beams.
Example:
SEARCH DATA 1
SOLVE
SAVE STRESSES 1
AVERAGING STRESSES STRAINS
ELEMENTS DIMENSION 3
RETURN
SAVE RESULTS 1
SEARCH DATA i
SOLVE (STIFFNESS) (MASS) (DAMPING or BUCKLING)
SAVE (STIFFNESS) (MASS) (DAMPING or BUCKLING) j
Buckling matrices are computed from displacements for 3D continuous media and shell
elements, and from forces for beam elements. In axisymmetric shell option, buckling matrices
are computed from forces contained in the STRESS file.
Example 1: normal utilization
SOLVE
SOLVE BUCKLING
Notes:
l If the structure contains at the same time beams and shells elements, it is necessary to do SOLVE
STRESSES before SOLVE STIFFNESS BUCKLING.
l If the imposed load used for non linear computation contains temperatures or imposed strains, the
user must verify the coherence between displacements and the load used to compute buckling
matrices (time function under SOLVE are not taken into account).
Note
l Elimination methods are more efficient with sparse mode matrix storage. These methods are not
available for non-symmetrical matrices. The EXTRACT 100 label can be used to obtain more fully
detailed printout.
SOLVE HARMONIC
In this case, it is not possible to execute individually the calculation steps of the solution
process.
l The file STIFFNESS holds for each element as many stiffness matrices and load vectors as
modified for each harmonic, and a further contribution is made to the results file.
l During the calculation of forces, each harmonic is analyzed in turn, by selection of the
This command transforms a displacement file, assumed to hold temperature values, into a file of
temperature cards which may be used in thermoelastic analysis and eventually creates the
corresponding thermoelastic loads. The number of temperature cards is equal to the number of
loads of the initial displacement file.
GENERATE This label is used to chain a thermal transient analysis with a
thermoelastic analysis by generating thermoelastic loads from thermal
transient file. This label can also be used with a result file issued from a
thermal static analysis.
CARD Label CARD defines a list of the thermal computation cards to be taken
into account. Each card generates a load in the DATA file.
The command TEMPERATURE must precede the associated SOLVE command.
It is supposed that a displacement file (temperatures previously calculated by a thermal analysis)
which is compatible with the thermoelastic analysis (i.e. is for the same mesh) is currently in
existence.
Note In automatic processing of a temperature file,
l the file NODAL STRESS must not have been previously accessed or created, as it would be treated
as a file of temperatures,
l the temperature data held on the displacements file must be compatible with the thermoelastic
analysis.
c) TEMPERATURE TRANSIENT SHELL command
The results of a non-linear thermal transient analysis on a structure meshed with thermal shell
elements (code 24XX) with the THERMAL ENTHALPY option (see Heat Transfer
Reference manual - Chapter 6), can be used for a thermo- mechanical analysis with the
SPATIAL SHELL option.
The interface procedure linearizes the temperature field in the thickness of the element
(computation of mean temperature and temperature gradient in the thickness) and computes the
mean temperature at the nodes of the elements. Input data is as follows:
This procedure requires the DATA input file of the mechanical analysis and the TRANSIENT
input file of the thermal analysis. Thermomechanical loads must be defined by the usual syntax:
THERMAL i.
Label LINEAR obtains linearization of the temperature field in the thickness of the element.
Label CARD defines a list of the thermal computation cards to be taken into account. Each
card generates two loads, the first corresponding to mean temperature, and the second to the
temperature gradient of each node.
9.5.3 Examples
Thermal Static Analysis followed by a Thermoelastic Analysis
The two following input data illustrate the use of the command TEMPERATURE. The
difference between the two input data is emphasized by the use of bold characters.
DEFINITION DEFINITION
THERMAL TEST THERMAL TEST
OPTION THERMAL PLANE OPTION THERMAL PLANE
RECTANGLE RECTANGLE
X Y 10 I J 10 X Y 10 I J 10
MATERIAL PROPERTIES MATERIAL PROPERTIES
M / KX KY 10 M / KX KY 10
CONSTRAINTS CONSTRAINTS
LINE / TT LINE / TT
LINE YY 10 / TT LINE YY 10 / TT
LOAD LOAD
1 TEMPERATURE / NOTHING 1 TEMPERATURE / NOTHING
LINE YY 10 / TT = 200 LINE YY 10 / TT = 200
RETURN RETURN
SOLVE SOLVE
SAVE DATA RESU 1 SAVE DATA RESU 1
DEFINITION DEFINITION
THERMOMECHANICAL TEST THERMOMECHANICAL TEST
OPTION TWO-DIME TRANSLATION OPTION TWO-DIME TRANSLATION
THERMO THERMO
RESTART GEOMETRY RESTART GEOMETRY
MATERIAL PROPERTIES MATERIAL PROPERTIES
E / E 210000 NU 0.3 LX LY LZ E / E 210000 NU 0.3 LX LY LZ
0.001 0.001
CONSTRAINTS CONSTRAINTS
LINE/ UX UY LINE/ UX UY
LOAD LOAD
1 FORCE / NOTHING 1 FORCE / NOTHING
LINE YY 10 / FY -10 LINE YY 10 / FY -10
2 THERMOELASTIQUE / THERMAL 1 RETURN
RETURN
The thermal input data is followed by a static thermal analysis. In the thermo-mechanical input
data you can notice:
l Two loads in the first case: a mechanical load plus a thermal load on which it is specified that
the thermal distribution will be retrieved from a thermal result file. If this result file has several
load cases it would be necessary to define a line for each load.
l Only the mechanical load is defined in the second case.
The command TEMPERATURE precedes the command SOLVE, and does not have the same
action:
l in the first case, the temperature distribution is retrieved from the result file (RESU 1) and is
as a thermal load. The temperature distribution is also retrieved from the result file.
DEFINITION
THERMAL TEST
OPTION THERMAL PLANE
RECTANGLE
X Y 10 I J 10
MATERIAL PROPERTIES
E / KX KY 10
CONSTRAINTS
LINE / TT
LINE YY 10 / TT
LOAD
1 T VARIABLE / NOTHING
LINE YY 10 / TT = 1 FT 1
TABLE
1 / 1 0 0 1000 1000
RETURN
TRANSIENT NONLINEAR
ALGORITHM OPTIMIZE -1000
PRECISION ABSOLUTE
TIME
200 STEP 1 / STORE 1
RETURN
DEFINITION
THERMOMECHANICAL TEST
OPTION TWO-DIM TRANSLATION THERMOELASTIC
RESTART GEOMETRY
MATERIAL PROPERTIES
E / E 210000 NU 0.3 LX LY LZ 0.001
CONSTRAINTS
LINE / UX UY
LOAD
1 FORCE / NOTHING
LINE YY 10 / FY -10
RETURN
SEARCH TRAN 1
TEMPERATURE GENERATE CARD 1 TO 200
SOLVE
SAVE DATA RESU 3
l one load, which will be a unit force at each node. Its directions is that of the contact forces,
9.6.2 Files
Input DATA,
l initial problem loads ("real" loads),
l contact unit loads.
RESULTS,
l file obtained with the previous data (or coming from a "MERGE" procedure).
l output RESULTS, if storage has been required.
CONTACT Command
This command initializes the procedure.
EXTRACT controls extraction level:
<no> =0 or 1 extraction of values introduced by the user under CONTACT procedure.2
(default) extraction of RESULTS used by CONTACT procedure.
LOAD Command
<list of numbers> list gives numbers of contact loads.
Note
l When the number of the contact loads is upper than 200, the SIZEn label on the CONTACT
command, with n= number of contact loads, must be used.
GAP
Optional command used to introduce gaps.
<no> is the number of the contact load.
<value> is the value of this contact gap.
CALCULATE Command
Defines combinations of the "REAL" loads whose effect must be calculated.
PRINT prints results during computation
STORE specifies that a local result file will be created at the end of the
CONTACT procedure.
<combination of "real" loads> = <title> / <value> <no> ...
<no> number of the "real" load
<value> load coefficient in the current combination.
9.6.4 Example
CHAPTER 10 SUBSTRUCTURING
10.1 INTRODUCTION
10.1.1 Terminology
The finite element method leads to the subdivision of the structure to be analyzed into
subdomains of simple geometric form (triangle, quadrilateral, tetrahedron etc ...) known as finite
ELEMENTS, which are interconnected by NODES. The purpose of defining a mesh in this
way is to enable the construction of a discrete system of equations describing structural
behavior, by means of integration over regions of simple geometric form.
An alternative representation which corresponds more closely to reality is that in which the
structure is considered as an assembly of components of arbitrary geometric form known as
SUBSTRUCTURES, each of which may itself be an assembly of lower level substructures.
l In the following discussion, by LEVEL ZERO SUBSTRUCTURES, or
SUPERELEMENT, is meant a part of the structure which is modeled in the usual way by
elements, with which are associated nodes, material properties etc ...
l A LEVEL 1 SUBSTRUCTURE is a part of the structure formed by the assembly of several
PRIMARY NODES or INTERFACE of the substructure, the remaining nodes being known
as SECONDARY NODES.
The representation of a structure so defined may be symbolized by a tree in which there are
links between substructures of different levels.
Note that the highest order substructure, which represents the complete structure, has no
primary nodes. In the following text, SE denotes substructure and SE(ij) substructure i, of
level j.
10.1.2 Example
In Figure 10.1 are shown two substructures of level zero (types 1 and 2).
+ secondary node
0,+ primary node
The assembly of two substructures of type 1 and of two others of type 2 gives two substructures
of level 1 as shown in Figure 10-2.
The assembly of the substructures SE100 and SE200 so obtained, together with a third
substructure of level zero, forms the level 2 substructure which represents the structure shown in
Figure 10-3.
Figure 10-1
Figure 10-2
Figure 10-3
K.U = F (1)
where :
K : stiffness matrix of the substructure
U : nodal displacement vectors
F : load vectors (equivalent nodal forces)
Separation of primary nodes (p) and secondary nodes (s) gives the following system:
Where
Kss . Us + Ksp . Up = Fs (2)
Elimination of degrees of freedom at secondary nodes gives the following system of equations:
(4)
(5)
This procedure finally results in the calculation of displacements at all nodes of zero level
substructures. Forces and stresses are then computed from these displacements in the same way
as in analysis without substructuring.
Figure 10-4
Figure 10-5
Figure 10-6
l the basic load vectors Fi and initial strains SSi which depend upon applied loads.
Matrix
This file is created after the elimination of unknowns at secondary nodes ; it consists of two
parts :
l the matrix Kss-1.Ksp (see equation 4, paragraph 10.2) required in the calculation of
Load vectors
The structure of this file is the same as that of the matrix file. On it are held the load vectors of
equations 4 and 5, Chapter 10.2.
Displacements
Solution of the system of equations gives displacements of nodes of each substructure ; these
are saved on the displacements file.
INPUT FILES.
l the procedure consists of the creation, using data read from input files, of one or more files of
Figure 10-7
The purposes of this mode of operation are to enable modular programming and also to permit
the execution of the analysis in several stages, with the possibility of stopping and later
restarting the calculation at the beginning of any of the following calculation steps.
Calculation of zero level substructure stiffnesses
Input file : DATA
Output file : STIFFNESS
Calculation of zero level substructure matrix
Input files : DATA and STIFFNESS
Output files : MATRIX and LOAD VECTOR
Calculation of matrix of substructure of level greater than or equal to 1
Input files : DATA for the substructure whose matrix is to be calculated.
DATA
MATRIX
for each associated lower level substructure
LOAD
VECTOR
Output files: MATRIX and LOAD VECTOR
For the highest level substructure, the output files are MATRIX and DISPLACEMENTS.
Calculation of displacements of a substructure
Input files : DATA, MATRIX, LOAD VECTOR.
DATA and DISPLACEMENTS of the associated higher level substructure.
Output files : DISPLACEMENTS
Figure 10-8
In the example shown in 10.4.1 Re-Use of Substructures, substructure 4 is modified after the
initial analysis. The files STIFFNESS, MATRIX and LOAD VECTOR of all other
substructures are re-used.
It may be seen from this example that by judicious subdivision of the structure, a significant
amount of computer time may be saved. The way in which substructures are assembled is also
of importance, as is shown by the example above in which only substructures 4 and 100 are
affected by the modification.
The alternative means of assembly represented by the tree shown in 10.4.1 Re- Use of
Substructures would have resulted in all substructures on the path linking substructure 4 to the
highest level substructure 100 being affected :
Figure 10-9
Figure 10-13
Example:
The structure shown in Figure 10- 13 consists of a plate (SE10) loaded in the direction
perpendicular to its plane, connected to another plate (SE20). Under these conditions SE10 is
analyzed as a plate subject to bending (option 2), SE20 as a plate subject to plane stress
(option 7), and the shared nodes as if they were nodes of structure type shell (option 3).
The degrees of freedom of the plate SE10 correspond to the unknowns UZ, RX and RY. The
degrees of freedom of the plate SE20 subject to plane stress are UX and UY, or UY and UZ
after transformation of coordinates.
Equations are therefore obtained at shared nodes in terms of the unknowns UY, UZ, RX and
RY. It is therefore necessary to introduce a third substructure the nodes of which are those
shared by SE10 and SE20, and are subject to fictitious constraints in the directions UX and RZ.
This modeling strategy divides the size of the problem by at least 2 by calculating at most nodes
either 2 or 3 unknowns instead of 6, and results in a significant reduction of computing time. It
should also be noted that the analysis of an assembly of plane membranes in space as a three
dimensional shell requires the stabilization by fictitious constraints of all nodes not located at the
intersection of several membranes.
The possibility of locally representing each membrane as a plate subject to plane stress therefore
results in both a reduction in computing time and simplified data preparation.
input data file either previously created by the routine DEFINITION or copied from a
permanent file by the command SEARCH DATA.
In the command DEFINITION, the presence of the label SUPER ELEMENT distinguishes
input data of a substructure from that for a problem without substructuring.
n number of the superelement. (1 < N < 9999)
(FILE) i suppresses the creation of the specified output files and all associated
computation, where existing files are to be re-used. Permissible labels are
STIFFNESS, MATRIX and LOAD VECTOR.
i is the number of the file
(TRANS) (R) geometric transformation to be applied to data held on the files MATRIX
and LOAD VECTOR of the superelement, during assembly of a higher
level substructure.
Example:
RETURN
SUPER ELEMENT 11 STIFFNESS MATRIX LOAD VECTOR 10
SUPER ELEMENT 12 STIFFNESS 10 SYMMETRY LINE THETA =45
SUPER ELEMENT 13 STIFFNESS 10 ROTATION PSI = 90
Figure 10-14
Figure 10-15
Note
l for θ rotation round the z axis,
l for ϕ rotation round the y axis,
l This is contrary to common practice.
COMPONENTS ei ej ek ...
CONNECTION
n/(ni ei) (nj ej) (nk ek)
...
LOAD
1 TITLE OF LOAD CASE
ci ei cj ej ....
2 ...
RETURN
On the command SUBSTRUCTURE is given the same data as for the definition of a zero level
substructure, that is:
n number of the substructure
(FILE) i suppresses the creation of the specified output files and all associated
computation, where existing files are to be re-used. Permissible labels are
MATRIX and LOAD VECTOR.i is the number of the file.
(TRANS) (R) geometric transformation as described in Chapter 10.5.1.
Load cases are each specified as the sum of several load cases of associated substructures.
The command LOAD is followed by a series of pairs of lines of input data.
The first line of a pair begins at the first character with the load case number, which must be
sequential. The load case number is followed by the title of the load case.
On the second line of a pair, which must begin with a blank in column 1, are specified pairs of
numbers (ci ei) where ci is the number of the load case and ei that of the substructure, for each
contribution to the total load that is to be taken into account. A substructure number ei may
appear one time only on this line.
Where the input data defining connections and/or load cases of several substructures are
identical apart from the numbers of associated substructures, commands need not be repeated if
preceded by the command EQUIVALENCE.
The command EQUIVALENCE defines a correspondence between the number of an
associated substructure and the number referred to in the data defining connections or load.
In this command therefore are specified pairs of values (ei fi) where ei is the number of the
substructure for which connections or loads are to be defined, and fi is the number by which the
substructure is identified in subsequent data.
The command RETURN terminates the input data defining the substructure.
Notes
l If no error is found by the routine, the input data file of the substructure is automatically saved.
l The permanent input data files of all associated substructures are updated by the addition of data
essential for the execution of later stages of the calculation.
l Each substructure is assigned a status code which has the following signification:
l 0: isolated substructure (not connected to any other)
l 1: connected substructure
l 2: STIFFNESS file created
l 3: MATRIX file created
l 4: LOAD VECTOR file created
l 5: DISPLACEMENTS file created
l 6: RESULTS file created (end of analysis)
l In command CONNECTION a couple (ni, ei) can appear only once.
l There "waiting" connections can exist with only one couple (ni, ei), which will be used to define a
higher level substructure.
l For the higher level substructure loads, the edition labels (Chap.4.8.2) can be used. They will be
taken into account by editing the results of every superelements of the substructure.
Example:
order of level,
l assembly and reduction of load vectors of higher level substructures, in ascending order of
level. For the highest level substructure, this operation yields nodal displacements,
l calculation of nodal displacements of each substructure, in descending order of level,
Caution:
l Before using CALCULATE command data must be load in memory.
The syntax of the command CALCULATE which enables these operations to be carried out is
as follows:
CALCULATE (STEP) n
(STEP) is a label specifying the stage of analysis it is desired to reach for substructure n as
follows:
STIFFNESS calculation of stiffnesses (zero level substructure)
MATRIX calculation of triangularization matrix and reduced load vectors
FORCE or calculation of forces or stresses (zero level substructure)
RESULTS
Each calculation step may be requested individually by including in the command
CALCULATE the label ONLY. In this case, the user must check that all intermediate files
(input files) required during the step are available. Syntax is as follows:
Notes
l Permanent input data files of substructures are updated upon completion of calculation steps, so
that a record of the status of each substructure is maintained. These files must not be modified by
the user during the solution process.
l If the requested sequence of calculation is terminated prematurely because for example the job
time limit is exceeded or the system goes down, only the calculation step in progress at the time
the job aborted is lost. If the same job is later resubmitted, the calculation is automatically
restarted at the beginning of the step that was aborted.
l Calculation steps that create files that are to be re-used by several substructures must be executed
first.
l If thermal loading is applied to zero level substructures by using a file of temperatures, then the
calculation of stiffnesses and forces or stresses must be requested individually.
l For a zero level substructure, the COMPUTE DISPLA n (ONLY) generates a DISPLA n file
necessary to forces computation. Physic displacements are then contains in the RESU n file.
Examples:
SEARCH TEMPERATURE 1
CALCULATE STIFFNESS 10
SEARCH TEMPERATURE 2
CALCULATE FORCES 20
File
As with standard substructure commands, input and output files are automatically managed.
The user must never modify these files.
Input file
l DATA: data files of every superelement of the substructure
Output file
l DATA: modified data files: connections are now defined.
l LINK: eventually alphanumerical files where connections are described in the SYSTUS
input data format (command SUBSTRUCTURE, subcommand CONNECTION)
Input Data
CONNECTION <number>
LINK EPSILON <epsilon> FILE
RELEASE
n/p
CHECK (ONLY) <number>
RETURN
CONNECTION command
This is the first level command for the automatic link generation
<number > is the number of the highest level substructure on which "connection"
will be applied. If the whole substructure has just been defined, this
number can be omitted.
LINK Command (k)
EPSILON Coupling tolerance. This is the largest distance at which two nodes are
<epsilon>: considered to be connected. Default value is equal to 10 percent of the
minimum of the leading size of the elements, for all the elements. This
leading size is the length of the longest size of an element.
The default value is intended to ensure that no nodes associated with
the same element will be connected.
FILE : specifies that connections will be described on formatted data files as
in the standard "CONNECTION" subcommand of the
"SUBSTRUCTURE" command:
CONNECTION
<NUMBER > / Ni Ei Nj Ej
with:
<number> : external number of connection
(Ni Ei): number of the Ni node of superelement Ei
- One file, whose name is BINDk.TIT, is created for each superelement and can be used in
SUBSTRUCTURE command by ASSIGN and READ command.
- These files, whose name is BINDk.TIT, can of course, be used under SUBSTRUCTURE
command, with $ASSIGN and $READ utilities.
CHECK command (k)
Extracts substructure data for the tree-node k (level, connections, etc., see 10.5.2). If no tree-
node is specified, the complete tree-structure is examined.
(ONLY) : Only one tree-node is examined.
RELEASE command
This command permits the release of a connection generated by the A.L.G.
n : tree-node number (substructure number)
p : absolute number of the connection (given by CHECK)
The named connection is released throughout the tree-structure and renumbering is carried out
by the A.L.G.
Example
The following sequence presents a simple use of the A.L.G.
DEFINITION
SIMPLE EXAMPLE OF A.L.G
OPTION TWO DIMENTIONAL
GEOMETRY
NODES
1/00
2/10
3/11
4/01
ELEMENTS
1/1234
RETURN
SUPERELEMENT 1
SUPERELEMENT 2 XX = 1
SUPERELEMENT 3 XX = 1 YY = 1
SUPERELEMENT 4 YY = 1
SUBSTRUCTURE 12
COMPONENTS 1 2
RETURN
SUBSTRUCTURE 123
COMPONENTS 12 3
RETURN
SUBSTRUCTURE 1234
COMPONENTS 123 4
RETURN
CONNECTIONS
LINK
CHECK
RETURN
10.7 CONCLUSIONS
Substructuring gives SYSTUS increased problem solving power and greater flexibility. The
most significant advantages of the method are summarized below.
l There is no absolute limit to the size of a problem that may be analyzed (except memory
limitation, A.L.G.).
l Discretization of the structure is simplified by its subdivision into independent components,
on entry and postprocessing. Freedom to choose the system of axes most appropriate to each
substructure is an important advantage.
l Subdivision of the solution process into calculation steps leads to more efficient use of the
computer, with shorter runs and quicker recovery from incidents such as system failure.
l Economies of computing time can be very significant if advantage is taken of various
CHAPTER 11 DYNAMIC
11.1 INTRODUCTION
The DYNAMIC module determines eigenfrequencies and eigenmodes for a physical problem,
and computes dynamic response or sensitivity using a modal superimposition process. The
basic matrices required for the problem (STIFFNESS, MASS, DAMPING and BUCKLING)
must be computed beforehand using the SOLVE command (Chapter 9), as also the static results
which are needed in certain cases.
The computation modules can be classified in four categories:
l Eigenfrequencies and eigenmodes
l Dynamic responses
l Sensitivity analysis
l Utility modules
The stiffness matrix can take account of geometrical stiffness due to a state of static stress, the
mass matrix, and added mass resulting from interaction between the structure and a non-
viscous, incompressible fluid (limited to certain computation options).
Direct analysis methods are suitable for medium size problems, and modal synthesis for large
size problems.
l Thermal domain: conductivity and thermal capacity matrices.
l Buckling domain: stiffness and geometrical stiffness matrices. Critical load factors
A special module can be used to study harmonic response with any damping matrix.
Sensitivity analysis
This method can be used to analyze the influence of isolated stiffness, mass or damping
changes, on eigenfrequencies and the form of real eigenmodes.
Utility modules
The 4 utility modules supply relevant data to the computation modules (modal damping and
response spectra), or are used for special processing (effective mass and dynamic transient
values).
l complex eigenmodes.
l Furthermore, only this algorithm can be used for computation of complex eigenfrequencies
and eigenmodes.
algorithms apart from INVERSE POWER, can be used to indicate the value of α
(recommended value: 10 to 20% of where ω1 is the first elastic body circular
eigenfrequency).
l Identification of eigenfrequencies in a given frequency band. FREQUENCY (CIRCULAR
FREQUENCY) labels are used in the ITERATIVE or BLOCK SUBSPACE algorithm, and
the LIMIT label where the INVERSE POWER algorithm, to indicate the position of the
center frequency sought. Matrix K' is not necessarily defined positive (α negative),
explaining its limitation to certain algorithms. Due to their antagonistic function, the SHIFT
and FREQUENCY (CIRCULAR FREQUENCY) labels cannot coexist. In case of
buckling, the shift sign must be the sign of the coefficient, which characterizes the critical
load searched.
Ahead of the DYNAMIC module, the SUM procedure can be used to include geometrical
stiffnesses (BUCKLING file) resulting from a state of static stress, in the STIFFNESS file.
- Modification of the mass matrix:
The mass matrix is normally consistent. Experience shows that eigenvalues are obtained with
sufficient precision for large vibratory problems, if the matrix is concentrated on the diagonal.
Computation time and memory space gains can be considerable. The DIAGONAL label is
available with all algorithms (except GIVENS), and allows the use of a diagonal mass matrix.
For computation options with transnational and rotational degrees of freedom, and with certain
algorithms, the inertia values assigned to rotational degrees of freedom are then ignored (care is
required with certain local modes, such as straight beam torsion).
Utilization of the DIAGONAL label is not relevant for buckling problems.
- Simultaneous modification of stiffness and mass matrices:
GUYAN's elimination method, activated by the SELECT or SUPPRESS label, can be used to
assign stiffness and mass matrices to a reduced geometry. The size of the problem is thus
reduced for search of the eigenvalues. The reduced stiffness is obtained by static condensation,
and reduced mass by limited development of the dynamic equilibrium equation. The degrees of
freedom which should be retained on a priority basis are those which, when blocked, would
increase the eigenfrequencies of the system to the greatest extent. The GUYAN method is
useful with the GIVENS algorithm, which computes all eigenfrequencies, and is also used in
association with the LANCZOS BAND SUBSPACE algorithm, with which the reduced mass
and stiffness matrices can be retained for assignment to special elements (code X9XX).
Note
l In the case of harmonic options, the stiffness and mass matrices depend on the harmonic. There
are consequently as many eigenvalue problems to be solved as there are harmonic taken into
account.
also referred to as the K-orthogonality matrix φTKφ. These matrices must be diagonal, and
the diagonal term values depend on the vector norm.
l as for static analysis, computation of residues R = (λM - K) φ which are stored in the reaction
file.
l Computation of Rayleigh quotients. The difference with respect to the eigenfrequency gives
an idea of the accuracy of eigenmode computation.
These operations are always carried out with a consistent mass matrix. The residues depend on
the vector norm. If the structure has rigid constraints, their order of magnitude should be
compared with that of the reactions.
For free-free structures, the number of zero eigenfrequencies must be equal to the number of
eigenmodes of rigid bodies of the option used. If this is not the case with SUBSPACE
methods, it is advisable to increase the size of the computation subspace. In all cases, STURM
sequences can be used to check the number of eigenfrequencies situated in a given frequency
band.
Results and files
The eigenmode search module prints the eigenfrequencies, and in certain cases eigenmodes,
and creates three files which can be saved.
DATA number of loads, equal to the number of computed modes. Circular
eigenfrequencies and generalized masses associated with RESULT file
eigenmodes are contained in this file.
RESULT This standard file contains the types of result effectively computed
(displacement, reaction and force) for each eigenmode. The vector
norm in this file is that selected by the user:
elastic (ELASTIC label) φtKφ = 1
mass (MASS label) φtM φ= 1
maximum (MAX label): the largest component of each eigenvector is
1. This norm is taken by default.
MODE This file is required for modal method dynamic response computation
(XI.4). After computation of eigenmodes, this file contains the same types
of result as the RESULT file. The eigenvector norm is the elastic norm.
The procedures described in paragraph 11.2.1.3 are used to include other
types of result (forces, accelerations, stresses and strains) in the MODE
file, and where appropriate, to align rigid body modes with the MASS
norm.
Note
l Rigid body modes can not be normed with elastic norm. Their number is automatically calculated
from the shift value α and the eigenfrequencies.
If the ELASCTIC norm is selected, rigid body modes are normed (RESULT and MODE files) with
respect to the shifted stiffness matrix (φtMφ =1/α).
If the MASS or MAX norm is selected, rigid body modes are normed (MODE file) with respect to
the mass matrix (φtMφ =1).
Input Data
Commands common to all methods
Input files
- DATA
- STIFFNESSES (K)
- MASSES or BUCKLING (M)
MODE command:
This command is mandatory. It is followed by a label (GIVENS, INVERSE POWER or
SUBSPACE), specifying the selected eigenmode search method (GIVENS by default if the
model size is below 500 d.o.f, if not, use SUBSPACE LANCZOS). Selection of one of the 4
available SUBSPACE methods must be specified.
SHIFT required for free-free structures. An α suitable value is of the order of
0.1 to 0.2 , where ω1 is the first non-zero eigenfrequency of the
structure. When no value of α is specified, an automatic determination
method, based on stiffness and mass matrices trace, is used. In this case
the SHIFT label must be the last one of the line ; the order of size of
α with respect to must be checked.
BUCKLING required for a buckling eigenmodes and critical loads search. If the critical
loads are not all of the same sign, they are sorted by increasing absolute
value order, except with SUBSPACE BAND where they are sorted by
algebraic values. INVERSE POWER, SUBSPACE BLOCK and
ITERATIVE SUBSPACE may be used to search critical loads around a
given value.
DIAGONAL the mass matrix is diagonalized before computation. This label must
not be used for buckling problems, and is not taken into account with
the GIVENS method.
SIZE t this specifies the subspace size for the SUBSPACE method, but is not
taken into account with the SUBSPACE BLOCK method
By default the subspace SIZE is computed automatically.
For large systems:
. 2nv + 1 for LANCZOS method
. nv (1 + 0.16) for ITERATIVE SUBSPACE (mini nv +8, maxi nv
+32)
VECTORS command
This command is not essential. If it is absent, nv (number of computed eigenfrequencies and
eigenmodes) is taken as at value of 1. All eigenfrequencies are computed with the GIVENS
method.
NORM this specifies the norm for the eigenmodes printed, and stored in the
RESULT file. The MAX norm is used by default. There is a choice
of the following norms: MAX, MASS (φtMφ=1) or ELASTIC
(φtKφ=1).
PRINT i controls printout of eigenmodes.
i: number of eigenmodes printed
by default: no printout.
FORCES this label is used to obtain the forces associated with the eigenmodes
(modal forces). The accelerations are added to the results.
Note
l The eight generalized forces are obtained for the thick shell elements (2006, 2008, 2403, 2404
elements).
DYNAMIC
MODE GIVENS (SHIFT = α)
SUPPRESS ni, nj...
SELECT ni, nj...
(PRINT) VALUES n
VECTORS nv (NORM <LABEL>) (PRINT i) (FORCE)--
(ITERATION niter) (PRECISON ε)
RETURN
Eigenfrequency search for a free- free structure, using GUYAN's method of elimination.
Printout of the first 20 eigenfrequencies, computation of the first 5 eigenmodes with mass norm,
and storage of displacements and reactions in the RESULT file.
INVERSE POWER method
MODE command:
STRUCTURAL The standard sparse method is used. The mass matrix is consistent.
In the lack of STRUCTURAL label, the resolution method can be specified
by the METHOD command.
In the presence of DIAGONAL label mass matrix is diagonalized (the
method can be specified by the SCALING label).
PRECISION ε precision requested for computation of eigenmodes (default value: 10-3).
ITERATION n maximum number of iterations for computation of an eigenfrequency
(default value: 20).
If fo and f 1 are different, the SHIFT optional label is used to define the number of central
frequencies to be taken into account for the search of nv eigenfrequencies.
For buckling, the command LIMIT LOWER co (UPPER c1) is used to define central critical
loads.
METHOD Command: optional (useless in the presence of the STRUCTURAL label).
It permits to choose a type of solver of the linear equation systems. By default a direct method
of resolution of symmetrical systems is chosen. To search eigenfrequencies in a frequency band
(label LIMIT), GAUSS factorization is necessary. The NON SYMMETRICAL label must be
specified for the iterative method.
For more information about this command and the various solvers available see Chapter 21.
FILE Command: optional
This command is used to impose an initial deformation contained in a RESULT file (nload:
load number). An initial deformation suitable to centre frequency leads to the reduction of
computation time.
Note
l The CONDITIONING label on the MODE command is used to obtain a valuation of the initial
numerical conditioning of the linear system associated with the stiffness matrix (static
conditioning).
Example 1: Inverse power method
Search for 2 eigenfrequencies and eigenmodes round 2500, with special checks for the iterative
algorithm. Eigendisplacements (MAX norm) are stored in the RESULT file.
Example 2: Inverse power method : First critical load of a vertical post imposing the
initial deformed shape.
DEFINITION
:
:
LOAD
1 AXIAL FORCE/ NOTHING
N / FZ – 1000
2 TRANSVERSAL FORCE / GX GY 1. NOTHING
RETURN
SOLVE
SOLVE BUCKLING
DYNAMIC
MODE INVERSE POWER BUCKLING
VECTOR 1
FILE 2
METHOD …
RETURN
SEARCH DATA 1
SOLVE STIFFNESS MASS
DYNAMIC
MODE INVERSE POWER CONDI ITER 5
METHOD…
RETURN
Note:
l If the mass is null, a unity mass matrix is taken into account.
ITERATIVE SUBSPACE method
SEARCH DATA 1
RENUMBER
1 17 23
RETURN
SOLVE STIFFNESS MASS
DYNAMIC
MODE SUBSPACE ITERATIVE SHIFT
VECTOR 10 NORM MAX
RETURN
SAVE DATA DISP 2
The algorithm uses the band storage mode, and the data are reclassified before computation of
the stiffness and mass matrices. The first 10 eigenfrequencies and eigenmodes of a free-free
structure are computed, with automatic shift determination. The reclassified data and results are
stored.
Example 2: Iterative subspace method
The 5 frequencies closest to 50 Hz are computed with a scaling factor of 1000, using a diagonal
mass matrix and a special iteration check.
Example 3: Iterative subspace method (specific solver)
DYNAMIC
MODE SUBSPACE (SHIFT = α) (DIAGONAL)(SIZE t)
VECTORS nv (NORM <LABEL> ) (PRINT i) (FORCE) --
(PRECISON ε) (STURM (FREQUENCYf)or (CIRCULARFREQUENCYp))
RETURN
Search for 5 eigenfrequencies and 5 eigenmodes, normalized with respect to stiffness (elastic
norm).
b-An algorithm using the BAND storage mode is also available.
DYNAMIC
MODES SUBSPACE BAND (SHIFT = α) (DIAGONAL(SCALING) (SIZE t) --
VECTORS nv (NORM <LABEL> ) (PRINT i) (FORCE) --
(PRECISON ε) (STURM (FREQUENCY f )or (CIRCULAR FREQUENCY p))
RETURN
The data are the same as below. The DIAGONAL label can be precised by SCALING.
Example: Lanczos subspace method: BAND
DYNAMIC
MODES SUBSPACE BAND SIZE 15
VECTORS STURM FREQUENCY 50
RETURN
Search for frequency lower than 50 with a maximum subspace size 15 (min. of 101 and 15).
Lanczos SUBSPACE method with GUYAN option
This algorithm is identical to that described in (Chapter. 11.2.1. Input Data Part). The function
of DIAGONAL and SIZE labels is the same. The Lanczos algorithm, for the searching of
eigenvalues and eigenvectors is preceded by a GUYAN reduction of the problem. The
GUYAN's reduced matrices can be stored (REDUCTION label). Eigenfrequencies and
eigenmodes associated with reduced matrices are computed. The METHOD command is used
to choose the employed method for GUYAN's reduction.
DYNAMIC
MODE SUBSPACE (GUYAN)
SUPPRESS
10 TO 15 / UX RZ
10 TO 30 STEP 10 INT / UY
VECTORS 20
RETURN
DYNAMIC
MODE SUBSPACE DIAGONAL (GUYAN)
SELECT FREQUENCY 100
/ UX UY UZ
VECTOR 10 STURM FREQUENCY 15 INTERRUPTION
RETURN
Guyan's elimination - diagonal mass matrix. Elimination of all rotations, and selection of
translational degrees of freedom for which the ratio of the diagonal terms of K to those of M
exceeds (200.π)2 . Search for frequencies lower than 15 with a maximum subspace size 21. If
the minimum frequency associated with secondary degrees of freedom is lower than 15, the
computation will stop.
DYNAMIC
MODE SUBSPACE REDUCTION
SELECT
10 12 15 30 / UX UY UZ
RETURN
Guyan's elimination and creation of STIGEN and MASGEN files, containing reduced
stiffnesses and reduced masses associated with 4 nodes, with selection of translation degrees of
freedom. The first eigenfrequency and the correspondant eigenmodes, associated with the
reduced matrices are computed.
DYNAMIC
MODE SUBSPACE BLOCK nvpb (MEMORY) --
(FREQUENCY (or CIRCULAR FREQUENCY) f )--
(SHIFT α) (DIAGONAL (SCALING)) (EXTRACT j)
METHOD …
RETURN
BLOCK nvpb This label is mandatory where the Lanczos subspace method is used with
the block option.
nvpb is the number of vectors per block.
Recommended value: maximum order of multiplicity in the case of
multiple frequencies.
The default value nvpb (1), can be left out if no number appears on the
command after the BLOC label.
MEMORY For little models, the computations executed in the central memory are
quicker.
FREQUENCY f Computation of frequencies round a given frequency.
or CIRCULAR These two labels, which cannot coexist, are used to specify the centre
FREQUENCY frequency round which the frequencies are computed.
By default the lower frequencies are computed.
SHIFT this label cannot be used in association with FREQUENCY or
CIRCULAR FREQUENCY. For buckling α defines a central critical
load.
DIAGONAL this optional label SCALING modifies the mass condensation method
on the diagonal (Chapter 12.3.3).
EXTRACT j j = 10 ; detailed control of eigenvalues convergence.
VECTORS command: specific label
ITERATION niter maximum number of LANCZOS iterations. Subspace size is niter x nvpb
(by default niter=2nv+2).
maximum value of niter: min (niter, ndr /nvpb)
ndr : number of unconstrained degrees of freedom of the structure
nvpb : number of vectors per block
PRECISION rtol Convergence criterion used for the eigenvalue convergence test.
By default, rtol = 10–4 (rtol = λi - λi - 1/λi -1)
STURM is used to adjust the computed eigenmodes number.
nv begins the number of modes between 0 and f1.
If FREQUENCY label is present on MODE card nv is the number of
modes of [f-f1, f+f1].
initial nv and niter define the maximum subspace size.
METHOD Command: optional
This command is only available when the DIAGONAL label is used.
It permits to choose a type of solver of linear equation systems. By default a direct method of
resolution of symmetrical systems is chosen. To search eigenfrequencies in a frequency band
(label FREQUENCY or CIRCULAR FREQUENCY), the label NON SYMMETRICAL
must be specified for the iterative method.
For more information about this command and the various solvers available see Chapter 21.
Example 1: LANCZOS SUBSPACE method with BLOCK option
5 eigenfrequencies are computed round 200 Hz. The number of vectors per block is 2, the
maximum number of iterations is 30 (subspace size 60), and the elastic norm is selected.
Example 2: LANCZOS SUBSPACE method with BLOCK option
The number of eigenfrequencies between 0 and 50 Hz are searched. The eigenmodes obtained
are computed with a subspace maximum size 61. The factorization associated to STURM
sequences increases the computation time.
Optional labels:
VECTORS i TO Selection of vectors processed. By default, all vectors are processed.
j STEP k or The selection can be used to reduce the size of the modal base before
(CARDS) dynamic response calculations (Chapter 11.4). For this reason, the
vector number is this of their order in the selection list.
SEARCH nf If the MODES file has been made permanent by a SAVE command, it
can be accessed without the need to copy this file, simply by giving the
file number accompanied by the SEARCH label.
EXTRACT Printout of computed modal forces.
THICK In the presence of thick shells, this label leads to take into account the
shear forces.
The label THICK preserves the 8 components of force in the file
FORCE (in RESU*.TIT) and in the file STRESS (file STRE*.TIT),
instead of 6 usually stored components for the shells. These 8
components of force are, in order of storage:
NX NY NXY MX MY MXY TX TY
The shell elements to which this label is applied are:
- 2403, 2404
- 2406, 2408
- 2006, 2008
It is not taken into account in the following post-processing procedures:
- SHELL (ref Chapter 13.13)
- CONVERT RESULTS STRESS
- COMBINE FORCE/STRESS
Note
l If it exists, the local RESULT file is destroyed by this procedure.
Output files:
MODES containing displacements, reactions, forces and accelerations.
Example:
SEARCH DATA 1
SOLVE STIFFNESS MASSES
DYNAMIC
MODES SUBSPACE
VECTORS 8
RETURN
SAVE DATA RESULTS 2 MODES 3
SOLVE FORCES MODAL
SAVE MODES 2
SEARCH DATA 2
SOLVE FORCES MODAL SEARCH 3 CARDS 1 TO 7 STEP 2
SAVE MODES 3
SEQUENCE command
This command converts a local RESULTS file into a MODES file. The MODES file contains
all the quantities of the RESULTS file. The local data file must be the file resulting from
eigenmode computation. If the local displacement file exists, the MODES file can also contain
accelerations (ACCELERATION label).
For all options, this command can be used to align rigid body modes on the MASS norm,
which generally means that modal quantities of a homogeneous order of magnitude can be
obtained for all eigenmodes. (RIGID nvs label)
For the continuous media options (TWO-DIMENSIONAL and THREE-DIMENSIONAL),
elementary forces in the elements, and elementary stresses and strains at the nodes, obtained by
SOLVE FORCES and AVERAGE (cf. 12.2 ) can be replaced by Tresca and von Mises
equivalent stresses (EQUIVALENT STRESS label). This provides an alternative way of
obtaining equivalent stresses from a SPECTRAL ANALYSIS or STOCHASTIC DYNAMIC
study.
Input files:
DATA resulting from eigenmode computation.
RESULTS can contain the 5 normal local files (displacements, reactions, forces, node stresses
and node strains).
Input data:
Optional labels:
ACCELERATION If the local displacement file exists, the MODES file will contain
accelerations
RIGID nvs If nvs rigid body modes exist, this label is used to align these modes on
the MASS norm .
EQUIVALENT For the continuous media options (TWO-DIMENSIONAL and THREE-
STRESSES DIMENSIONAL), the elementary quantities (forces, node stresses and
node strains) are replaced by the corresponding Tresca and von Mises
criteria in the MODES file.
SHELLS In the SPATIAL SHELL, the RESULT file from the eigenmodes
computation, being treated by SOLVE FORCES and STRESS
(cf.Chapter 13.10), the elementary stresses on the three skins (superior,
average, lower) are replaced by the TRESCA and Von Mises criteria
corresponding, in the MODE file (six components).
BEAMS In the SPATIAL BEAM option and the RESULT file from the eigenmodes
computation being treated by SOLVE FORCES and STRESSES (cf.
Chapter 13.10), the normal stresses at four points of the computation of
each section are obtained at the two extremities of the beam (inertial forces
not taking into account).
VECTORS i TO j Selection of vectors processed. By default, all vectors are processed. The
STEP k selection can be used to reduce the size of the modal base before dynamic
response calculations (Chapter 11.4). For this reason, the vector number is
this of their order in the selection list.
Output files:
VECTORS can contain displacements, reactions, forces, node stresses, node strains
and accelerations.
Special Problems
Enumeration of eigenvalues (Sturm sequences)
Basically, a number of eigenfrequencies between zero and a given frequency can be obtained.
This method involves triangulation or Gaussian reduction of matrix K-λM where λ is the pivot
associated with the given value. The cost of this reduction is of the order of one static solution.
The number of eigenvalues between two pivots, λi and λj, can be obtained by two successive
computations.
The predetermined number of eigenvalues can differ slightly (in particular if the structure
contains multiple eigenvalues) from the number obtained by eigenvalue computation. This
corresponds to rounding errors.
Input data:
The DIAGONAL label can be used to take account of the effect of a DIAGONAL mass, if an
eigenvalue search method is to be used with this option.
The optional EXTRACT label permits to obtain the localization and the value of negative
pivots.
Vibration modes:
Pivot value is as follows, according to the selected label:
FREQUENCY λi = (2.π.f)2
CIRCULAR FREQUENCY λi = ω2
COEFFICIENT λi = c (c>0)
Buckling modes:
Only the label COEFFICIENT may be used to define the pivot (where c>0 or c<0).
The METHOD label is used to choose the utilized solver (see chapter 21). It must lead to a
total triangulation of GAUSS. Iterative solvers (incomplete triangulation) here are less accurate.
Example
SUMMATION
Eigenvalues equation is ( ):
The specific physical quantity is thus the problem eigenvalue of the problem defined by the
calorific conductivity and calorific capacity matrices K and C. For stable systems λ is a real
positive number. Classical output data gives the square root of the absolute value of λ. It is
possible to obtain the real value of λ by processing the problem like a buckling problem.
Example:
SEARCH DATA 1
SOLVE STIFFNESS MASS
SAVE STIFFNESS MASS 1
DYNAMIC ; square root
MODES SUBSPACE
VECTORS 10
RETURN
SEARCH DATA STIFFNESS 1
ASSIGN 33 MASS1.TIT BINARY
DYNAMIC ;λ
MODES SUBSPACE BUCKLING
VECTORS 10
RETURN
...
DYNAMIC
MODEFLUID STRUCTURE(CONDENSATION)(NONE)(VERIFY)(FACTORs)--
(SIZEt)(RESTART)(PRESSURE)(METHOD i)
FRONT FLUID
n1, n2 ... nj
FRONT STRUCTURE
i1, 12 ... ik
VECTORS nv (NORM <LABEL>) (PRINT i) (PRECISION ε) (FORCES)
RETURN
MODE command:
FLUID calls the LANCZOS algorithm adapted for FLUID-STRUCTURE
STRUCTURE problems
CONDENSATION used to stop computation after condensation of the fluid on the structure.
This gives a new mass file, which can be stored under the name MATRIX,
and which contains (M+Ma) . This file is required for SPECTRAL
ANALYSIS (see Chapter 11.4).
NONE is used to minimize the size of necessary files. This label can not
coexist with the CONDENSATION label.
At the end of the eigenmode calculation, the MATRIX file can not be
created.
VERIFY in the presence of this label, there is no computation.
Only
- the number of solid, wet and fluid nodes,
- the size of the fluid matrix,
- the global size for the fluid condensation on the structure (generally the
maximum size of the computation)
are computed.
This command can be used to try different fluid or structure fronts
(optimisation of the bandwidth).
FACTOR s definition of scaling factor s. By default, this factor is computed
automatically. It may be used in the presence of a large number of
constraint elements very stiff (for example: wet wall both quoted).
SIZE t is used to define the size of the subspace (by default this size is computed
automatically).
RESTART a first computation done and a MATRIX file saved, this label is used to
remake an eigenmodes computation without reform the added mass matrix.
The parameters which can be modify are the followings:
- the number of eigenmodes,
- the numerical parameters (size, factor, precision),
- the element material properties which doesn't have a influence on the
solid mass (method 1), or all element material properties of solid (method
2).
The geometry and the FRONTS (FLUID, STRUCTURE) can't be
modified.
PRESSURE in the presence of this label, the fluid potential UP and the fluid
reactions FP are printed for each fluid node joined.
METHOD i i=1
the added mass matrix is associated with all the solid nodes (method by
default).
i=2
the added mass matrix is associated with all the "wet" solid nodes.
i = 3 (SHELL FLUID option).
Added mass matrix is not explicitly formed. This method minimizes file
size and computation time when a usual number of eigenmodes is
searched. In the presence of the PARTICIPATION label, participation
factors corresponding with a uniform movement of supports are
computed. A mono-spectral analysis is possible (no representative
factor computation, for the reactions: the linked mass effects are not
taken into account).
FRONT command:
This command is used to renumber fluid and structure nodes to optimize the performance of the
band storage method employed for condensation and eigenmode search. The bandwidths
obtained are printed. If independent fluid volumes exist, at least one node of each volume must
appear in the list associated with FRONT FLUID, and likewise for the rigid structure (FRONT
RIGID). This command is similar to RENUMBER (see Chapter 6.3.2).
VECTOR command:
nv number of eigenvalue and eigenvector pairs to be computed.
NORM selection of the norm for the eigenvectors to be printed or stored in the
RESULT file with one of the following labels.
ELASTIC :normalization with respect to the stiffness matrix(φT K φ = 1)
MASS : normalization with respect to mass matrix (φT (M+Ma) φ = 1)
MAX : maximum vector component has the value 1 (default norm).
PRINT i eigenvector print check.
PRECISION ε the label PRECISION permits a larger precision in the eigenmodes
searching (by default ε = 10-2).
FORCE Modal forces computation
Results
Files created by the DYNAMIC routine are as follows:
l DATA
l RESULTS
l MODES
M + Ma : method 1
Ma : method 2
The eigenvectors form a result file with nv loads. Problem data are modified for compatibility
with the result file (number of loads).
The norm in the MODES file is elastic. For spectral analysis purposes, forces and modal
accelerations can be added to this file, using the FORCE label or the SOLVE FORCES
MODAL command (Chapter 11.2.1).
The DATA, RESULTS and MODES files are limited to the shell structure. They consequently
only contain the nodes and elements of the shell, with the nodes given in the renumbering order
obtained with the FRONT STRUCTURE command. Potential fluid displacements are
retained. The number of degrees of freedom for these files is not the same as for conventional
SHELL options (5 for SHELL HARMONIC, 7 for SHELL FLUID).
Note
l The file size much depends on the number of structure nodes.
In SHELL FLUID option
method 1: 5 000 structure nodes MATRIX file 3.6 Go
method 2: 5 000 "wet" structure nodes
The time of computation much depends on the number of fluid nodes and "wet" structure
nodes.
Example: 1- Initial computation
Calorific capacity
Notes
l There's no boundary conditions. Matrix HF is generally a singular matrix and the SHIFT command
has to be used. Only harmonic options with a harmonic number not equal to zero lead to a non-
singular matrix HF.
l In incompressible case (c = ∞), the number of eigenmodes is limited to the number of free surface
nodes.
DEFINITION
SLOSHING - THERMAL ANALOGY : COMPRESSIBLE FLUID
OPTION THERMAL AXISYM HARMONIC SYMETRIC 1
TWO-DIMENSIONAL
NODES
1 / 0.
11 / 18.288
2 / 0. 12.192
12 / 18.288 12.192
ELEMENT
101 / 1 11 12 2
1001 / 2 12
EDGE
2001 / 1 11 36
2101 / 2 12 36
3001 / 1 2 24
3101 / 11 12 24
MATERIAL PROPERTIES
E 101 TO 964 / KX KY KZ 1.RHO 1. C 1*-4
E 1001 TO 1036 / H 1. RHO 1. C 0.1019
VERIFY
RETURN
SOLVE STIFFNESS MASS
DYNAMIC
MODES SUBSPACE
VECTORS 10
RETURN
Analogy:
3D elements : K=1
2D elements : K=0
we obtain a generalized eigenvalue equation, with the same form as that given in Chapter
11.2.1:
(K* + λ M*) Ψ = 0
where K* is not defined positive
λ is a complex number (λ = α + j ω)
Problems can only be processed where the K, M and C physical matrices are symmetrical and
real.
Only one eigenfrequency and eigenmode search algorithm is available, this being the
generalized LANCZOS algorithm for SUBSPACE method with the BLOCK option
(11.2.1.2.6). The K matrice is assembled in memory, according to the choices done using the
METHOD command. By default, it is assembled in band form (renumbering is required). The
matrices M and C are on file. In the presence of the MEMORY label, they are assembled in
memory in band form (renumbering is required).
The M matrix can be condensed on its diagonal.
The matrices (K, M and C) are real, and the eigenfrequencies are real, or conjugated two by
two. Only eigenfrequencies with a positive or zero imaginary part are retained, together with
the corresponding eigenvectors.
Note
l The fact that the K* and M* matrices are double the size of the K, M and C matrices results in
high computation times for complex eigenfrequencies and eigenmodes. An approximate method,
described in Chapter11.6.3, can be used with matrices which can be non-symmetrical. This
method is interesting for large models, for which the real eigenfrequencies have already been
computed.
Input Data
Input files:
DATA
STIFFNESS SEARCH
MASS or SOLVE STIFFNESS MASS DAMPING command
DAMPING
DYNAMIC
MODE SUBSPACE BLOCK nvpb DAMPING (MEMORY)REAL --
((DIAGONAL (SCALING)) (EXTRACT j) (SHIFT (α))
VECTOR nv (ITERATION niter ) (PRECISION rtol)--
(NORM <LABEL>) (PRINT i) (FORCE)
METHOD
RETURN
The input data is very similar to that of a real eigenmode search using the SUBSPACE
algorithm with the BLOCK option (Chapter 11.2.1.input data part). Only the particularities are
described here. The use of STURM sequences is not possible.
MODE command:
DAMPING initiates the complex eigenmode search. In the case of high damping
system, some eigenvalue can be real (ω = 0). The optional label REAL is
used to obtain only the complex eigenmodes. If the number of the real
eigenvalues is important the size of the subspace must be controlled by
using the ITERATION label of the VECTOR command.
The SUBSPACE and BLOCK nvpb labels are required, as only SUBSPACE BLOCK
permits a complex eigenmode search. The frequency search round a given frequency is not
possible. For free-free structures α is approximately 10 to 20% of the first non-zero circular
eigenfrequency. When no value of α is specified, an automatic determination method, based on
stiffness and mass matrices trace, is used. In this case the SHIFT label must be the last one of
the line ; the order of size of α with respect to the first non-zero circular eigenfrequency must
be checked.
EXTRACT j j = 10 ; detailed control of eigenvalues convergence.
VECTOR command:
ITERATION niter maximum number of LANCZOS iterations. Subspace size is niter x
nvpb.
by default, niter = 4 nv + 2
Maximum value of niter: min (niter, 2 ndr/nvpb).
ndr : number of non-constrained degrees of freedom of the structure
nvpb : number of vectors per block
FORCE This label is used for computation of forces in beams, with the
following formulation:
F = -KU - λCU- λ2 MU
METHOD Command: optional
This command is only available when the DIAGONAL label is used.
It permits to choose a type of solver of linear equation systems. By default a direct method of
resolution of symmetrical systems is chosen. To search eigenfrequencies in a frequency band
(label FREQUENCY or CIRCULAR FREQUENCY), the label NON SYMMETRICAL
must be specified for the iterative method.
For more information about this command and the various solvers available see Chapter 21.
Results
Output files
l DATA: contains the names of the loads corresponding to the eigenmodes
l RESULTS: contains eigenmodes for the selected norm
l MODES: contains elastic norm eigenmodes (ΨT K* Ψ = 1)
Notes
l Eigenmodes Ψ are normalized with respect to matrices K*, M*.
l Only eigenmodes with a positive or zero imaginary part are retained, together with the
corresponding eigenvectors.
l Two cards in the MODES file (or loads in the RESULT file) correspond to each stored complex
eigenmode, such that:
l The real part corresponds to an odd number
l The imaginary part corresponds to the next even number.
Complex eigenfrequencies: λ = α + jω = 2π (α ' + jf)
If the system has damping greater than critical damping on an eigenmode, the corresponding
eigenfrequency is real (ω = 0). All real eigenfrequencies are retained.
If the system has damping less than critical damping, λ and its conjugate are
eigenvalues, and only λ is retained. Eigenfrequencies are stored in ascending value of ω.
A positive value for α corresponds to an unstable mode.
Rayleigh quotients:
Effective verification of eigenfrequency and eigenmode computation (λ' = λ).
In free-free structures, Rayleigh coefficients have no meaning for rigid body modes.
Example 1 : Complex eigenmodes
Six eigenfrequencies and eigenmodes are computed for a damped structure. The number of
vectors per block is 2, as we know that double frequencies exist. The result and eigenvector
files obtained contain 12 loads or cards, with odd numbers corresponding to real parts, and even
numbers to imaginary parts.
Example 2: Complex eigenmodes with integration of geometric stiffness
SEARCH DATA 1
SOLVE
SOLVE BUCKLING
SUMMATION
SAVE STIFFNESS 1
SOLVE STIFFNESS MASS DAMPING
SEARCH STIFFNESS 1
DYNAMIC
MODE SUBSPACE BLOCK 1 DAMPING
VECTOR 5 NORM MASS
RETURN
Caution must be taken as BUCKLING and DAMPING are on unit 33. Only the SOLVE
STIFFNESS MASS DAMPING command is reliable for damping problems defined in
elements (Chapter 11.2.3).
With
The eigenmode equation is solved by factorisation QR. All the eigenvalues are then obtained.
The eigenmode vectors are eventually obtained by inverse iteration algorithm.
Notes
l Unsymmetries
They can be included in the STIFFNESS, MASS, DAMPING files (element having X8XX, X9XX
codes), or can be described directly in the modulus for all types of element.
l Models size
This modulus is particularly suitable to study the model instabilities (positive real eigenvalue) with
weak size (some thousand degrees of freedom) and having large non- symmetries. For tall models,
the modulus PSEUDO MODAL (Chapter 11.6.3)for which the input data is very similar can be
used.
Input Data
Input files:
DATA
STIFFNESS SEARCH
MASS or SOLVE STIFFNESS MASS DAMPING command
DAMPING
DYNAMIC
MODESCOMPLEX(RATIO r)(ALL)(GYROSCOPIC ω0)
ELEMENTS
nd TO ne (STEP nf) (INTERNAL) / TABLE t3 t2 t1 (ELEMENT n)
ng / TABLE …
TABLE
t1 / (mnij, Kmnij)…
t2 / (mnij, Mmnij)…
t3 / (mnij, Cmnij)…
VECTOR nv (PRINT i) (RAYLEIGH) (THRESHOLD ε)(MODULUS ρ)
RETURN
SAVE DATA RESULTS (MODES) j
by default:
no VECTOR command r=10
VECTOR command without PRINT N=2nv+1
ALL optional label useful to save the zero imaginary part eigenvalues.
GYROSCOPIC In option SPATIAL BEAM, the beam being carried by Oy axis, this
ω0 command is used to take into account the anti-symmetrical damping
matrices due to the gyroscopic effects (element of1002 code, ponctual
mass). ω0 is the uniform rotation velocity (rd/s). In output, gyroscopic
damping matrices are stored on DAMPING file which can be used on
HARMONIC RESPONSE modulus (see Chapter.11.5).
Notes
l the elements concerned do not necessarily have X9XX codes,
l no reference frame can be defined with the ELEMENT command,
l matrices which can be non-symmetrical must be given in full, and are added to those contained in
the STIFFNESS, MASS and DAMPING files.
l non-symmetrical stiffness, mass and damping matrices are given with the TABLE procedure if
t3t2t1 >0 or in files t3t2t1 <0 (Chapter 11.6.3. Stiffness, Mass and Damping Matrices Stored in
Files Part).
TABLE command (optional)
This command defines the stiffness, mass and damping matrices in the way described in the
ELEMENT command.
VECTOR command (optional)
This command is used to release the computation of modal deformed shape.
nv number of complex modal deflections requested (by default nv =1).
PRINT i This label controls result printout.
In the absence of this label, no results are printed.
In the presence of this label:
- i>0, results are printed for the first i modes.
- i=0, results are printed for all modes.
RAYLEIGH In the special case where stiffness, mass and damping matrices are
symmetrical, the Rayleigh quotients give an indication of solution quality.
THRESHOLD is used to erase the stable eigenvalues for which the absolute value of the
imaginary part is lower than ε (choice of the eigenmodes).
MODULUS is used to erase the eigenvalues for which the modulus value is greather
than ρ (choice of the eigenmodes). By default ρ=106.
Results
In the lack of the VECTOR command, only the complex eigenvalues are printed.
If the VECTOR command is present, 3 files are created.
- DATA : contains the titles of the 2nv loads corresponding to the eigenmodes.
- RESULT (MODES): eigenmodes with the norm MAX.
Notes on complex eigenvalues (
l All complex eigenvalues with a positive part are retained because they correspond to a physically
unstable system.
l Only complex eigenvalues with a positive or zero imaginary part are retained. In the lack of the
ALL label, the eigenvalues with a very small imaginary part are not retained.
l The retained eigenvalues are stored in ascending value of size. Their critical damping fraction
is associated with.
SEARCH DATA 1
DYNAMIC
MODE COMPLEX RATIO 5 ALL
RETURN
SEARCH DATA 1
DYNAMIC
MODES COMPLEX
ELEMENTS
1 2 / TABLE 30201
TABLE
1 / 1233 k 2133 k
2/
3/ 1233 c
VECTOR 20 MODULUS 1 * 7
RETURN
SAVE DATA RESULTS (MODES) 2
The generalized coordinates (α,β) of each substructure become the unknown factors of the
problem, and are related by reduced mass (M r ) and reduced stiffness (K r ) matrices. These
matrices can be obtained from vibration tests, allowing the integration of experimental super-
elements. The substructures are assembled by writing the equality of the displacements and
physical forces on the constrained degrees of freedom, using Kr and Mr .
Selection of a method means choosing a function base, namely eigenmodes ( ), and static
modes (U). Static modes do not have an orthogonality property with respect to eigenmodes,
thus partially justifying the multiplicity of methods. The seven methods currently available are
described in Chapter 11.3.2.
There is only one substructuring level. Substructures are referred to as level 0 physical or
dynamic super- elements, and the complete structure is referred to as a level 1 dynamic
substructure.
Variants
Each free or mixed method has one or more variants, involving a number of supplementary
hypotheses:
l free and mixed ATT methods : 2 variants
The variants are less precise than the basic methods, which are used by default. Three levels of
precision are defined:
Precision 0
Attachment or residual eigenmodes are taken into account for computation of reduced stiffness
and reduced mass matrices. This level of precision concerns the free and mixed ATT and FRD
methods, but does not exist for the free and mixed FRF methods.
Precision 1
Attachment or residual eigenmodes are only taken into account for computation of the reduced
stiffness matrix, and are ignored for computation for the reduced mass matrix.
This level of precision concerns:
l free and mixed FRF methods
l variant no 1 for free and mixed ATT methods
l variant for free and mixed FRD methods
Precision 2
Attachment or residual eigenmodes are ignored for calculation of the reduced stiffness matrix
and reduced mass matrix.
This level of precision concerns:
l variant no 2 of the free and mixed ATT methods
This level does not exist for the free and mixed FRD methods.
In this case, static results corresponding to free constrained degrees of freedom are not
necessary. The free and mixed FRF methods allow the use of super-elements with rigid body
modes. The special input data for CALCULATE DYNAMIC MATRIX for these super-
elements, is described in Chapter 11.3.6.
The following table summarizes possible combinations between precision levels and free and
mixed methods.
Method Free and mixed Free and mixed Free and mixed
Precision Level ATT FRF FRD
by default (no
PRECISION = 0 by default
constrained stiffness)
PRECISION = 1 variant 1 by default (**) variant (**)
PRECISION = 2 variant 2(*) variant (*)
Notes
l An experimental type substructure (see Chapter 11.3.8) can use methods and variants marked (*)
or (**). A residual flexibility matrix restricted to free constrained degrees of freedom is not
necessary in cases marked (*), but must be supplied in cases marked (**)
l In all cases, the substructure can be without mass, or with mass and without eigenmodes.
l Precision level is selected using the RESIDUAL label of the CALCULATE DYNAMIC command
(Chapter 11.3.6).
Restrictions
Three methods which use frontier node displacements only as unknown quantities, have no
restrictions:
l fixed method
l free-FRD method
l mixed-FRD method
On the points of the structure to which more than two substructures are connected
For these points, a restriction is imposed on the type (fixed or free) of the constrained degrees of
freedom of the associated substructures:
l if free or mixed ATT or FRF methods are used, only one free constrained degree of freedom
Boundary Conditions
Connections with the exterior are described conventionally.
For constrained degrees of freedom between super-elements, there is a choice between fixed
interfaces, which stabilize the super-element and allow static modes, and free interfaces.
The degrees of freedom of fixed interfaces are constrained:
l stiffly
l elastically if constrained stiffnesses exist between super-elements.
The degrees of freedom of free interfaces are left free in all cases, and constrained stiffnesses,
where these exist, are defined in the physical properties of the substructure.
Static Loads
Normally only exterior loads for which dynamic response computation is to be executed for the
complete structure are defined. Loads included in the STIFFNESS file for the super-element
are integrated by the ASSEMBLE DYNAMIC STIFFNESS routine (Chapter 11.3.6).
Static loads corresponding to static modes are generated automatically by the program
(CALCULATE DYNAMIC RESULT command), and the RESULT file has the number of
the STIFFNESS file for the super-element.
Nevertheless, computation of static modes can be initiated manually by including loads in the
data file, and applying the following rules:
l unit displacement loads imposed on fixed interface degrees of freedom are defined before the
degrees of freedom is the same as the order of appearance of these degrees of freedom in the
definition of the substructure links.
In this case, no external load must be defined.
l MASS
If rigid body modes exist, they must be aligned with the MASS norm (see SEQUENCE
command: Chapter 11.2.1).
If computation is restarted, three other pre-existent files can be associated with the super-
element (these files are not recreated):
l RESULT: containing static modes
l MATRIX: 2 separate files containing reduced mass and stiffness matrices.
Step 2
Creation of data file for level 1 substructure. This step is essentially used to describe
connections between super-elements (position and type of link), but also:
l Type of static modes associated with free links
l Number of eigenmodes effectively integrated for each level 0 super-element
l Connection stiffness which can be associated with the free links.
Step 3
Automatic computation of static modes for each level 0 super-element.
Step 4
Computation, for each level 0 super-element, of reduced mass and reduced stiffness matrices
(Kr, Mr).
Step 5
Assembly of Kr and Mr matrices to obtain assembled reduced stiffness matrix KR and
assembled reduced mass matrix MR.
Step 6
The eigenvalue equation is solved: (KR - ω2 MR) VR = 0
This gives circular eigenfrequencies ω and eigenmodes VR.
The eigenmodes VR contain the generalized coordinates (α, β).
Circular eigenfrequencies ω are the circular eigenfrequencies of the assembled structure.
Step 7
The physical eigenmodes for each substructure are reconstituted from VR, eigenmodes and
static modes.
Step 8
The physical eigenmodes of the substructures are assembled, to obtain the eigenmodes of the
complete structure. It is also possible to assemble the stiffness and mass matrices, for dynamic
response computation on the complete structure.
Steps 3, 4, 5 and 7 must be executed for each super-element.
The different steps are initiated by four specific commands, described in detail in Chapter
11.3.6.
Command Step
SUPERELEMENT DYNAMIC 1
SUBSTRUCTURE DYNAMIC 2
CALCULATE DYNAMIC 3 to 7
ASSEMBLE DYNAMIC 8
INPUT Files
Each level 0 super- element is assumed to possess the following files before step 1 is
commenced:
DATA
MODE eigenmodes
MASS elementary masses
STIFFNESS elementary stiffnesses and loads
The number of the MODE file must be different from the number of the super-element, in order
to avoid erasing this file on completion of step 7.
File RIGD
This file contains the reduced stiffness matrix for a level 0 super-element (Kr). For a level 1
substructure, the file contains the assembled reduced stiffness matrix (KR).
File MASD
This file contains the reduced mass matrix for a level 0 super-element (Mr). For the level 1
substructure, the file contains the assembled reduced mass matrix (MR).
FileVECD
This file is created for the level 1 substructure only, and contains circular frequencies ω and
eigenmodes VR in generalized coordinates.
OUTPUT Files
Two OUTPUT files are created for each level 0 super-element on completion of step 7. These
are:
MODE
RESULT
These files have the same number as the level 0 super- element concerned, and contain
restriction of level 1 substructure modes to the level 0 super-element. These are standard
SYSTUS files, which can be handled by post-processing procedures.
Furthermore, the DATA file for the super-element is compatible with the MODE and RESULT
files.
These files are assembled, together with the STIFFNESS and MASS files where appropriated,
in step 8. The eigenmodes are normalized on the elastic norm for the complete structure, in the
MODE file.
8
Step 1 2 3 4 5 6 7
DATAP DATA1
DATA1
STIF0 DATA1
VECD1 DATA1
MASS0 DATA0
DATA1 DATA1 DATA0
MODE0 DATA0 RIGD0
Input DATA0 MODE0,S
files DATA0 STIF0 MASD0
DATA0 RIGD1 STIF0
(RESU0)* MASS0
STIF0 MASD1 MASS0 STIF0
(RIGD0)* RESU0
RESU0 MASS0
(MASD0)* MODE0
MODE0
DATA1 DATAS
DATA0 DATA0
DATA0 DATA0 RIGD1 DATA1 RESUS
Output DATA0 MASD1 MODES
files RIGD0 RESU0,S
DATA1 RESU0 VECD1
MASD0 MODE0,S
STIFS
MASS
l 1: level 1 substructure
l S: complete structure
Note
l DATA0 files are modified throughout computation.
l SUBSTRUCTURE DYNAMIC
l CALCULATE DYNAMIC
l ASSEMBLE DYNAMIC
RETURN
or
SEARCH DATA ig
.
SUPERELEMENT DYNAM n MODE i MASS k STIFFNESS l --
(RESULT j) (MATRIX m) (SUPPRESS nvsupp) --
(EXPERIMENTAL iexp)(XX x 0 YY y 0 ZZ z0)
SUPERELEMENT DYNAMIC
If the program detects no error, the DATA file for the level 0 super- element is created
automatically.
XX, YY, ZZ define translations taking into account at the ASSEMBLE DYNAMIC
level.
CONNECTIONS Command
This command describes the connections between the level 0 super-elements.
Each connection defines a link between two or more level 0 super-elements.
A level 0 super-element is involved in a connection by one degree of freedom only.
A degree of freedom for a level 0 super-element node must not be involved in more than one
connection.
The fixed method, and the free- FRD and mixed- FRD methods have no restrictions.
Restrictions for the other methods are described in Chapter 11.3.2.
M : connection number introduced
ei : level 0 super-element number introduced
ni : level 0 super-element node number introduced
dofi : constrained degree of freedom label.
The labels used are UX, UY, UZ, RX, RY and RZ
Type : type of link. The label is FIXED or FREE.
STIFFNESS Command
This optional command is used to introduce constrained stiffnesses on the free interface degrees
of freedom of level 0 super-elements, where the connection is not infinitely stiff.
ei : level 0 super-element number introduced
ni : level 0 super-element node number introduced
dofi : connected degree of freedom label
The labels used are UX, UY, UZ, RX, RY and RZ
ki : constrained stiffness value.
Note
l Constrained stiffnesses are introduced on definition of the substructure, and are taken into account
in step 4 (computation of RIGD for the level 0 super-element).
MODIFYCommand
Where the physical properties for certain super-elements are modified, and in order to avoid
repeat computation for unchanged super-elements, this command can be used to recompute
only those super-elements concerned.
ei, ej number introduced for a level 0 super-element for which the physical
properties have been modified, and for which computation must be
repeated.
Notes
l If the program detects no error, the DATA file for the level 1 substructure is created
automatically.
l The permanent data files for all level 0 super-elements are modified by the inclusion of information
essential for subsequent computation phases.
l Each level 0 super-element, and the level 1 substructure, has a status code which is interpreted as
follows:
l 0 : isolated super-element or substructure
l 1 : super-element constrained
l 3 : super-element or substructure RIGD or MASD files created
l 5 : level 0 super-element MOD and RESULT files created.
l Level 1 substructure VECD file created.
l (the code value can be obtained by the VERIFY DYNAMIC n command – Chap 7.2.4)
<LABEL >n This label indicates the computation step required for super-element
number n
The labels used are:
l RESULT
l MATRIX
l MODE
RESULT n indicates that only static results for level 0 super-element
number n (step 3), or all elements of level 1 substructure number n, are
to be computed.
MATRIX n indicates that files RIGD and MASD for level 0 super-
element number n (step 4), or level 1 super-element number n (step 5)
are to be created.
MODE n indicates that MODE and RESULT files for level 0 super-
element number n (step 7), or file VECD for level 1 substructure
number n (step 6) are to be created. In this last case and in the presence
of the WHOLE label, the step 7 will be done for all levels 0 super
element of the substructure number n.
VECTOR n1 n1 is the number of modes required.
For the level 1 substructure, by default (n1)1 = 1
For the level 0 super-elements: : (n1)0 ≤ (n1)1 by default : with
VECTOR label (n1)0 = 1
without VECTOR label (n1)0 = (n1)1
(with MODE label)
EXTRACT in the presence of this label, print RIGD and MASD matrices computed
during the step.
(with MATRIX label)
PRINT i This label is used to control printout of the following:
- static results for level 0 super-elements
- modes for level 0 super-elements or the level 1 substructure.
If the label does not exist, no modes are printed.
If i > 0, the first i modes are printed
If i is null or absent, all modes are printed.
It is advisable to use this label, in particular for computation of static
results, in order to avoid overloading the listing.
NOCHECK n3 Having computed the modes for the level 0 super-element (step 7) or
the level 1 substructure (step 6), this label authorizes recomputation of
the modes, changing n1, n4, n5, ε or i.
n3 = 1 authorizes recomputation of the modes for the level 0 super-element
(step 7), without repeating the preceding steps.
n3 = 10 authorizes recomputation of the modes for the level 1 substructure (step 6),
without repeating the preceding steps
(with MODE label).
OTHER If the free-ATT (or free-FRF) method is selected for computation of all
super-elements composing the substructure, the OTHER label
authorizes computation of level 0 super-element n using the free-FRF
(or free-ATT) method (with MATRIX label - step 4).
RESIDUAL j authorizes utilization of variants for the 3 free methods (XI.3.2.4):
j = 0 integration of residual stiffness and mass
j = 1 integration of residual stiffness only
j = 2 residual stiffness and mass ignored.
(With MATRIX label - step 4)
the following table indicates possible values according to selected free
method:
Free method ATT FRF FRD
j=0 j=0
by default by default
j=1 j=1 j=1
RESIDUAL j
by default
j=2 j=2
Nvs number of rigid body modes for the super-element (nvs ≤ 6).
m1,... m solid body mode generalized masses as supplied by the
nvs
PARTICIPATION module (cf. 11.4.2), if not set at 1 by using
SEQUENCE RIGID nvs (cf. 11.2.1.).
l Physical properties and stiffness and mass matrices for the complete structure (Chapter.
11.3.6.).
The fixed constraints of nodes forming part of a fixed link of the level 1 substructure are deleted
from the physical properties, and the associated constrained stiffnesses are cancelled, as these
have no physical reality in the assembled structure.
Following modal synthesis, the files obtained are equivalent to those created following
computation of modes for the complete structure.
However, constraints between super- elements corresponding to the links of the dynamic
substructure are not included in the STIFFNESS file, making conventional static computation
of the complete structure impossible, and restricting certain possibilities of the dynamic response
module (see Chapter 11.4). Unless LINK command can be used (Chapter 11.3.6.)
Assembly of data and modes
Input files:
DATA for level 1 substructure
DATA for level 0 super-elements
MODE for level 0 super- elements, with each MODE file containing the modes of the
structure, restricted to the degrees of freedom of the level 0 super-element nodes.
These are the result files for modal synthesis step 7.
Output files:
DATA for the complete structure, reduced to the following entities:
l Geometry
l Physical properties
NORM <label> defines the norm for the eigenmodes in the RESULT file. The
associated label can be MAX, MASS or ELASTIC. The MAX norm
is used by default.
PRINT j edition of the j first eigenmodes. In the lack of PRINT, no printing. In
the presence of PRINT without j value , edition of all eigenmodes.
Example: Static results generation for a spectral analysis with static correction.
VERIFY
RETURN
SOLVE
SAVE DATA RESULTS 2001.
STIFFNESS stiffness matrix for the complete structure, containing the specified second
members.
MASS mass matrix for the complete structure.
The MASS file for the complete structure contains the elementary mass matrices. The
STIFFNESS file for the complete structure contains the elementary stiffness matrices, and the
applied load vectors specified by the LOAD command.
Input data:
.
.
RETURN
LOAD command
This command is used to create the applied load vectors contained in the STIFFNESS file for
the complete structure:
j1 Title/
(ek nk ck) (ei ni ci)...
l Spectral analysis
l Stochastic analysis
The only difference occurs with computation of modal forces (PARTICIPATION module –
Chapter 11.4.2), and stems from the fact that it is not possible to create the static result file.
Multi-excitation computation is not possible. For computation with forces on the free nodes, the
PARTICIPATION command must be followed by the LOAD command (Chapter. 11.4.2).
If the command LINK can be used, (Chapter 11.3.6), the restrictions of the previous paragraph
doesn't exist any more.
Sensitivity computation (Chapter 11.6) is also possible.
DEFINITION
MODAL SYNTHESIS * FIXED-FIXED METHOD * SUPER1 *
OPTION PLANE BEAM
ONE-DIMENSIONAL
NODE
1/ 0. 0.
6/ 0. 1.
ELEMENTS
1/ 1 6 5
MATERIAL PROPERTIES
M 1 TO 5 /E 2.*11 NU 0.30 RHO 7800 AX 0.0005 IZ 2.6*-8
CONSTRAINTS
N 1 6 / UX UY RZ
RETURN
SAVE DATA 1
SOLVE STIFFNESS MASS
SAVE STIFFNESS MASS 1
DYNAMIC
MODE
VECTOR 6
RETURN
SAVE MODES 11
DEFINITION
MODAL SYNTHESIS * FIXED-FIXED METHOD * SUPER2 *
OPTION PLANE BEAM
ONE-DIMENSIONAL
NODE
6/ 0. 1.
11/ 0. 2.
ELEMENTS
6/ 6 11 5
MATERIAL PROPERTIES
M 6 TO 10 /E 2.*11 NU 0.30 RHO 7800 AX .0005 IZ 2.6*-8
CONSTRAINTS
N 6/ UX UY RZ
RETURN
SAVE DATA 2 ...
Example 1: continued
SEARCH DATA 1
SUPERELEMENT DYNAMIC 100 STIFFNESS MASS 1 MODES 11
SEARCH DATA 2
SUPERELEMENT DYNAMIC 200 STIFFNESS MASS 2 MODES 21
SUBSTRUCTURE DYNAMIC 1000
COMPONENTS 100 5 200 5
CONNECTIONS
1 / 100 6 UX FIXED 200 6 UX FIXED
2 / 100 6 UY FIXED 200 6 UY FIXED
3 / 100 6 RZ FIXED 200 6 RZ FIXED
RETURN
CALCULATE DYNAMIC RESULT 100
CALCULATE DYNAMIC MATRIX 100
CALCULATE DYNAMIC RESULT 200
CALCULATE DYNAMIC MATRIX 200
CALCULATE DYNAMIC MATRIX 1000
CALCULATE DYNAMIC MODE 1000 VECTORS 5 TOUT
ASSEMBLE DYNAMIC 1000 FILE 1100
ASSEMBLE DYNAMIC MASS 1000 FILE 1100
RETURN
SEARCH DATA 1
DEFINITION
MODAL SYNTHESIS * METHOD MIXED F.R.F. * SUPER1 *
OPTION PLANE BEAM
REPEAT
CONSTRAINTS
N 1/ UX UY RZ
RETURN
SAVE DATA 3
SOLVE STIFFNESS MASS
SAVE STIFFNESS MASS 3
DYNAMIC
MODE
VECTOR 6
RETURN
SAVE MODES 31
SEARCH DATA 2
DEFINITION
MODAL SYNTHESIS * METHOD MIXED F.R.F.* SUPER2 *
OPTION PLANE BEAM
REPEAT
CONSTRAINTS
N 6/ UX UY RZ
LOAD
1 FORCE AT END /
N 11 / FX -2
RETURN
SOLVE STIFFNESS MASS
SAVE DATA STIFFNESS MASS 4
SEARCH DATA 3
SUPERELEMENT DYNAMIC 300 STIFFNESS MASS 3 MODE 31
SEARCH DATA 4
SUPERELEMENT DYNAM 400 STIFF MASS 4 MODE 21 RESULT 2 MATR 200
SUBSTRUCTURE DYNAMIC 2000 ; (METHOD F.R.F.)
COMPONENTS 300 5 400 5
CONNECTIONS
1 / 300 6 UX FREE 400 6 UX FIXED
2 / 300 6 UY FREE 400 6 UY FIXED
3 / 300 6 RZ FREE 400 6 RZ FIXED
RETURN ...
DYNAMIC
COMBINATION 5 RELATIVE
PRINT MAXI MODULUS
DISPLACEMENT 11
REACTION 1
FORCE 1
CARD 0 TO 1000 STEP 100
CARD 1 TO 1000 STEP 1 STORE
RETURN
SAVE DATA TRANSIENT 2201
END
Added mass effects are taken into account by the two rigid body translation modes and rigid
body rotation mode. These modes are aligned with respect to the MASS mode.
DEFINITION
Beam with load at end - EL.1: manual computation of static modes
OPTION TWO-DIMENSIONAL
TWO-DIMENSIONAL
NODES
1/
2 / 10.
3 / 12.
11 / 0. 1.5
12 / 10. 1.5
13 / 12. 1.5
ELEMENTS
101 / 1 2 12 11
501 / 2 3 13 12
EDGES
1001 / 1 11 3
1011 / 1 2 20
1031 / 2 3 4
MATERIAL PROPERTIES
E 101 TO 160 501 TO 512 / E 2.*11 NU 0.3 RHO 7800. H 1.
CONSTRAINTS
N 1 11 1001 1002 / UX UY
LOADS
1 N 12 FX / NOTHING
N 12 / FX 1.
2 N 12 FY / NOTHING
N 12 / FY 1.
3 N 13 FX / NOTHING
N 13 / FX 1.
4 N 13 FY / NOTHING
N 13 / FY 1.
VERIFY
RETURN
SAVE DATA 1
SOLVE
SAVE RESULT 1
DYNAMIC
MODE SUBSPACE
VECTORS 10
RETURN
SAVE STIFFNESS MASS MODES 1
DEFINITION
BEAM WITH LOAD AT ONE END - EL.2 : EXPERIMENTAL ELEMENT -
OPTION TWO-DIMENSIONAL
GEOMETRY
NODES
12 / 10. 1.5
13 / 12. 1.5
VERIFY
RETURN
SAVE DATA 20
SEARCH DATA 1
SUPERELEMENT DYNAMIC 100 STIFF MASS 1 MODE 1 RESULT 1
SEARCH DATA 20
SUPERELEMENT DYNAMIC 220 EXPERIMENTAL 222
Note
l With some fixed links, practically necessary with an experimental super-element, the eigenmodes
order of magnitude, (normed with the MASS norm), and static modes can be different. This leads to
an eigenvalue problem not correctly solved and necessitates the use of a low offset which does not
change the circular and eigenmodes values.
CONNECTION command
DYNAMIC activate the D.A.L.G.
FREE the generated links are supposed to be free (by default FIXED).
TRANSLATION
the links are generated only for the translation degrees of freedom of
the considered option (the fluid structure options are not allowed).
num number of the static substructure
LINK command
FILE compulsory label to generate the ASCI LINKnum.TIT file including
the dynamic links.
Limitations:
l Maximum number of substructures connected into one point : 6,
SEARCH DATA 50
SUPER ELEMENT 1
SUPER ELEMENT 2 ZZ 5.
SUPERELEMENT 3 ZZ 10.
SUBSTRUCTURE 100
COMPONENTS 1 2 3
RETURN
CONNECTION DYNAMIC TRANSLATION 100
LINK FILE ; file generation LINK 100.TIT
RETURN
SEARCH DATA 50
DYNAMIC
MODES …
VECTORS nv
RETURN
SAVE STIFFNESS MASS MODES 50
SEARCH DATA 50
SUPER ELEMENT DYNAMIC 101 STIFNESS MASS MODES 50
SUPER ELEMENT DYNAM I C 102 STIFNESS MASS MODES 50
SUPERELEMENT DYNAMIC 103 STIFNESS MASS MODES 50
SUBSTRUCTURE DYNAMIC 1000
COMPONENTS 101 nv 102 nv 103 nv
LINK
$READ LINK 100.TIT
RETURN
END
The unknown quantities of the problem are U 1 , absolute displacements of free nodes, and R 2 ,
reactions on the constrained nodes.
Driving movement Ue corresponds to the static deflection imposed by the supports
(Ue = - K11-1 K12 U2). Relative movement Ur is defined with respect to this static deflection
(U1 = Ue + Ur). This verifies the equation:
The relative movement is projected into eigenmode base ϕ , computed with the constrained
degrees of freedom 2 (U 2 = 0), U r = ϕ q . Generalized coordinates q checks the dynamic
equilibrium equation:
In this case:
l generalized coordinates q verify independent equations,
l relative movement, reflecting dynamic behavior, depends on the imposed forces and
Generally, the behavior of the structure is linear, thus allowing linear accumulation of the
different excitations. A number of non-linearities (stops, friction) can be taken into account,
where the excitation is defined in temporal terms.
Solver limitations:
The number of eigenvectors of the modal base (before a possible selection) cannot exceed 640.
Definition of Excitations
Four different computation modules are available, corresponding to four different types of
excitation definition.
MODAL ANALYSIS module (Chapter 11.4.3)
Excitation is defined in terms of time. The modal transient response, followed by the physical
transient response are obtained in the relative or absolute reference frame. A number of non-
linearities (stops, friction) can be taken into account, by introducing pseudo-forces in the second
member of the equations.
HARMONIC RESPONSE module (Chapter 11.4.4)
Excitation amplitude is defined in terms of frequency. The stationary response to a harmonic
excitation (amplitude, phase) is obtained. Variation of the excitation frequency gives the
stationary or transient transfer function.
STOCHASTIC module (random vibrations - Chapter 11.4.5)
Associating its power spectrum density (PSD) with each ergodic zero mean value stationary
type excitation, the following are computed using transfer functions:
l response PSD.
If cross spectra exist between a number excitations, only the real part of these cross spectra can
be taken into account.
If damping is introduced for the complete structure, in BASILE form (C = αK + βM), the α
and β coefficients are transformed, for each circular eigenfrequency ω, into a critical damping
fraction:
DAMPING command
This command is included in the MODAL ANALYSIS module ( Chapter 11.4.3 ) and
HARMONIC RESPONSE module ( Chapter 11.4.4 ). Critical damping fractions γ ,
corresponding to a given physical damping factor, can be obtained for elastic eigenmodes using
the MODAL DAMPING routine (Chapter 11.7.1 ), or by complex frequency computation
(Chapter 11.2.3).
Notes
l The number of rigid body modes is defined by the command RIGID of the computation command.
For HARMONIC RESPONSE MODAL this label is not always necessary. Without this label, only
BETA type damping is taken into account.
l If the BETA and MODAL labels coexist, the damping of elastic body modes is defined with BETA
and the damping of rigid body modes is defined with MODAL.
l In case a variable damping is used (label MODAL) conjointly with a selection of modes, the selection command
must precede the command DAMPING. Modes are given by their serial number in the selected list.
Examples:
1 - Elastic modes with damping, rigid body modes without damping:
modes 1 to 6 : β = 0
modes 10 to 18 :
modes 19 and 20 :
2 - Elastic and rigid body modes with damping:
modes 1 to 3 : β = 1
modes 4 to 6 : β = 3
modes 7 to 20 :
l Imposed acceleration
o for a mono-excitation study
US : static deflection associated with an uniform acceleration type load (GX, GY or GZ).
In a multi-excitation context, the approximation of an uniform acceleration is most of the
time sufficient.
φk : eigenmode
There is a pseudo-mode for each seismic direction, and the existence of these pseudo-modes
explains why the PARTICIPATION module associated with the SPECTRAL ANALYSIS
module (Chapter 11.4.6 ) is different from that associated with the MODAL ANALYSIS
module (Chapter 11.4.3) and the HARMONIC RESPONSE module (Chapter 11.4.4).
Representativeness Criteria
Associated with analysis of a static deflection, U S , these criteria require the existence of a
RESULT file containing this deflection, except in a MONO-EXCITATION context where this
deflection is generated by the program.
"Support movement" type load
US : static deflection associated to the U2 movement of the supports.
Taking account of the orthogonality of the eigenmodes and residual modes:
we can write:
Input Data
* Participation factors for MODAL ANALYSIS or HARMONIC RESPONSE module
Input files:
DATA from static computation
STIFFNESS from static computation (if FORCE type loads)
RESULTS from static computation, required if MULTI is present or FORCE present without
RIGID label or LOAD command
MASSES
MODES
Input data:
PARTICIPATION command:
The printout list also contains the modal base representativeness criteria: effective masses,
displaced masses and strain energies.
Example 1: Mono-excitation of the supports (PLANE BEAM option)
The participation factors corresponding to uniform movements of the supports on the x and y
axes are computed. DATA 100 contains 2 modal loads.
Example 2: Mono-excitation of the supports and force (SHELL option)
The participation factors corresponding to the first form imposed displacement type loads and
the first two force type loads are computed.
DATA 100 contains 6 modal loads.
PARTICIPATION command
nv computation of MONO-SPECTRAL participation factors for the first nv
eigenmodes. These correspond to uniform unit accelerations of all
supports, following the translational degree of freedom for the option
concerned. In the presence of the DOUBLE label, computation are done in
double precision.
nv can be replaced by FREQUENCY fo or CIRCULAR FREQUENCY
ω0.
The eigenmodes of which the frequency is lower than fo (or the circular
frequency is lower than ω0) will be taken into account.
In the lack of nv, FREQUENCY or CIRCULAR FREQUENCY, the
number of eigenmodes taken into account is that of the MODE file.
DIRECTION optional label used to define, for a mono-spectral study, a seismic
direction. Only this direction is taken into account subsequently (GN
label).
CX CY CZ c1, c2, c3 are the director cosines for this direction.
RESIDUAL- TERMS command : if residual modes or static corrections are taken into
account.
Computation of residual terms, and inclusion of these modes in the MODE file. This MODE
file can only be used by a consistent SPECTRAL ANALYSIS, about level seismic directions
taken into account, with the RESIDUAL-TERMS command.
GN or GX, GY, Seismic directions
GZ
if the GN label is present, redefinition of director cosines c1, c2, c3 is
mandatory.
n or nx, ny, nz numbers of imposed acceleration loads defining static deflections.
MODE mc label and mandatory value
number of dynamic eigenmodes taken into account.
FREQUENCY fc optional label. The value of mc is recomputed, integrating all eigenmodes
for which the eigenfrequency is not greater than frequency of cutting fc.
<LABEL> at least one of the following result LABELS must be specified:
DISPLACEMENTS, REACTIONS, FORCES, STRESSES,
STRAINS or ACCELERATIONS. These results are printed if the
PRINT label is present.
In the presence of ACCELERATION label, the keyword
QUADRATIC starts up the calculation of residual modes in
acceleration with a quadratic summation of the contribution of
eigenmodes taken into account.
EQUIVALENT if the MODE file has been obtained by SEQUENCE (STRESS) (SHELL)
(BEAM) EQUIVALENT command (see section 11.2.1.3.2), these same
labels must be present to process FORCES, STRESSES, STRAINS
(contains transformation of the file static RESULT).
PARTICIPATION MULTI command: required in the case of a MULTI-SPECTRAL study.
nv computation of multispectral participation factors corresponding to the
first nv eigenmodes, and residual modes where appropriate.
m number of static data imposed displacement loads (movements of the
support) taken into account.
Notes
l Computation can be fractionized.
l If fluid-structure interaction is taken into account and if modes have been computed with the
method 1, the MASS file must be replaced by the MATRIX file. If modes have been computed with
the method 2, the MASS file must be recomputed. As a static RESULT file is not possible, the
RESIDUAL-TERMS and MULTI labels are not authorized.
l The RESIDUAL-TERMS and MULTI labels are not authorized after computation of eigenmodes by
the modal synthesis method (no static RESULT file).
l In the case of a MULTI-SPECTRAL study, the data file imposed displacement values are taken into
account for computation of participation factors, and pseudo-acceleration spectra must
consequently be defined.
l In a fluid-structure interaction computation, the effective mass of the rigidly constrained degrees
of freedom is not accounted for. Nevertheless, this is compatible with a spectral dynamic analysis.
Results
Files:
DATA with generalized masses and participation factors.
MODES including residual modes (if RESIDUAL-TERMS label present).
The printout list also contains the representativeness criteria associated with imposed
displacement type loads: effective masses.
Example 1: Mono-spectral study with static correction (PLANE BEAM or BIDIM option)
DEFINITION
:
LOAD
1DIRECTION X / GX 1. NONE
2 DIRECTION Y / GY 1. NONE
RETURN
:
SEARCH DATA MASS RESULT 1 MODES 2
DYNAMIC
PARTICIPATION FREQUENCY 33
RESIDUAL-TERMS GX 1 GY 2 FREQ 33 DISPLA FORCE
RETURN
SAVE DATA 100 MODES 2
DEFINITION
:
CONSTRAINT
N 1 101 / UX (UY UZ RX RY RZ)
LOAD
1 NODE 1 X / NONE
N 1 / UX 1.
2 NODE 101 X / NONE
N 101 / UX 1.
RETURN
:
SEARCH DATA MASS RESULT 1 MODES 2
DYNAMIC
PARTICIPATION 10 MULTI 2
RETURN
SAVE DATA 100
DEFINITION
:
CONSTRAINT
N 1 101 / UX UY
LOAD
1NODE1 Y / NONE
N 1 / UY 1.
2NODE1 Z / NONE
N 1 / UZ 1.
3NODE101 Y / NONE
N 101 / UY 1.
4NODE101 Y / NONE
N 101 / UZ 1.
5DIRECTION Y / GY 1. NONE
6DIRECTION Z / GZ 1. NONE
RETURN
Analysis
Time-related modal analysis assumes execution of a number of modules.
In the most comprehensive case (absolute response for a multi-excitation context), five main
steps can be distinguished:
A - Static analysis (SOLVE).
B - Modal base (DYNAMIC MODES).
C - Participation factors (PARTICIPATION).
D - Evolution of generalized coordinates (MODAL ANALYSIS).
E - Evolution of physical quantities (COMBINATION).
Step A: Although this step is only absolutely necessary for the supports in a MULTI-
EXCITATION context, it is recommended in the case of force excitations, in order to be able to
compute the representativeness criteria of the selected modal base. Supports with imposed
movement form part of the constrained nodes.
Step B: This step computes the eigenmodes of the structure. Constraint conditions are
identical to those for step A.
Step C: This step computes modal forces, and various criteria characterizing the
representativeness of the modal base. The PARTICIPATION module is described in
paragraph 11.4.2.
Step D: This step computes the time-domain evolution of the generalized coordinates, taking
account of non-linearities. These modal acceleration velocity and displacement responses, as
also forces induced by non- linearities, are stored in the MODAL file. The MODAL
ANALYSIS module is described in paragraph 11.4.3.
Step E: This step computes the physical components for the desired results on the selected
elements and nodes. Where appropriate, these are stored in a TRANSIENT file, reduced to the
selected nodes and elements. A DATA file (also reduced) is associated with the TRANSIENT
file. The COMBINATION module is described in section 11.4.3.4.
nv number of the highest order eigenmode taken into account for computation.
DISPLA nd the first nd mono-excitation type support participation factors are used:
nd ≤ 3 (cannot coexist with MULTI).
MULTI nd the first nd multi-excitation type support participation factors are used
(cannot coexist with DISPLA).
FORCE nf the first nf imposed force type participation factors are used.
RIGID nvs there are nvs rigid body modes (eigenfrequency zero) for which generalized
mass has been recomputed (PARTICIPATION procedure).
ik number of the table defining the acceleration function associated with the
kth support movement type load (1 ≤ k ≤ nd).
(jl , j2,...,jnf) numbers of tables associated with the first nf force type loads defined in
the data.
In the more general case, a triplet is required for each of the nd support movement type loads.
l For a multi-excitation study, the complete syntax form is recommended, only response in the
absolute reference frame having a physical significance.
TABLE command
All tables must be defined under this command, or under TABLE in the DEFINITION
procedure. In the latter case, the TABLE command must not appear at this point.
Sequential file type tables (FILE) are not authorized for modal analysis, and direct access files
must be used in their place (see Chapter 4.15.4 - tabulated file functions).
DAMPING command (optional)
A precise physical description of the parameters for this command is given in paragraph 11.4.1.
No damping is taken into account if this command is not present.
STOP command (PRINT n) (RELEASE) (optional)
If this command is present, one or more of the five types of non-linearity described below
(unidirectional or cylindrical stop, unidirectional or plane frictional stop and material stop) is
taken into account.
PRINT n non-linearity status and temporal non-linearities evolution are printed every
n computation steps
RELEASE contact terminates in the stop when the contact force is zero (F=(ku+c )).
By default, the contact terminates if differential displacement u is zero.
MODAL The prediction of the solution is made on the generalized coordinates.
METHOD m if the ALGORITHM VARIABLE STEP is used, m> 1 changes the
prediction method of the impact time.
Notes
l In the case of force loads only, or in a mono-excitation context, the gap criteria are expressed in
the relative reference frame.
l For multi-excitation, the gap criteria are in the absolute reference frame. In the case of contact
with the exterior, the "MULTI fd, fv, fγ" label specifies the numbers of the tables describing the
displacement, velocity and acceleration of the external medium in each direction where relevant.
l When a local reference frame is defined, the angles PSI, THETA and PHI are expressed in
degrees.
a - Unidirectional stop
min min. and max. gap values, counted algebraically on the axis defined by (U)
max
k stop stiffness
c stop damping
Note
l Node n1 is stopped by the external medium, or node n2 if
b - Cylindrical stop
Friction-free contact be two concentric cylinders with local axis z, the gap being measured
radially .
UNIDIRECTIONAL
The frictional stop is considered to be unidirectional
n1, n2 same definition as before
min, max, k , c same definition as for a unidirectional stop (acting in z or y direction with
respect to the option)
FORCE f if this label is present, the frictional stop is said to be decoupled. The
two gap parameters of the stop must be defined. Tangential force is
calculated from f, NUST, NUDY, according to the value of EPSI.
Only tangential force is taken into account in the computation.
Is this label is not present, the frictional stop is said to be coupled. The
four parameters of the stop must be defined. Tangential force is
calculated from stop normal force and static (NUST) and dynamic
(NUDY) friction coefficients according to the value of EPSI.
Tangential and normal forces are taken into account in the computation.
PLANE The frictional stop is considered to be plane. Data are identical to those for
a unidirectional frictional stop (MULTI label followed by 3 triplets of
values).
e - Material stop
The material stop is used to take into account a force-displacement relation of plastic material
type.
This stop acts either in traction or in compression depending on the sign of the gap. If the gap is
zero, it acts in traction and compression, it is equivalent to an unidirectional stop.
Due to the history effect, restart calculations are not possible when the material behavior is
becoming plastic.
HISTORY ik ik ≥ 0 the stop takes into account the history effect (restart calculation is
impossible when the stop turns non-linear).
ik < 0 the stop is bilinear.
COEF xc xc > 1 the threshold variable MIN or MAX is divided by this parameter
when the stop turns non-linear.
ALGORITHM command (optional)
This command is followed by one or two labels (which can coexist), and is used to modify
certain parameters of the integration algorithm:
ZERO initial modal accelerations have imposed zero values, except if the initial
conditions are picked up from the MODAL file.
VARIABLE stop impact time between two consecutive computation times is
STEP nR recomputed. The minimum step is defined by
nR ( , by default nR = 5)
if this command is not present:
l Initial accelerations are computed by balancing modal forces.
l Impact times are not computed. Numerical integration using constant step is quicker but can
Modal displacements and velocities are picked up from card i of the MODAL file, and modal
accelerations are computed from the equilibrium equation. This then becomes card (i+1) of the
MODAL file.
In the lack of the CARD label, the last card of the MODAL file is restarted. By default (lack of
the INITIAL label on the TIME command), the initial time of the computation is the one
corresponding to the chosen card.
Initial physical conditions (TRANSIENT or RESULT file)
Relative movement displacements and physical velocities can be defined from the following:
l TRANSIENT file (i : card number).
l RESULT file:no: file number. Velocity values are zero. Load 1 is taken into account.
l list syntax, taking precedence over the other definitions for the nodes listed.
Physical velocities and displacements are used to compute MODAL velocities and
displacements. MODAL accelerations are deduced from the equilibrium equation, or taken to
have zero values (ALGORITHM command).
Note
l If a TRANSIENT file generated by COMBINATION (Chapter 11.4.3.) is picked up by the
RESTART CARD command, displacement and velocity values for all nodes must have already
been computed and stored.
TIME command
Definition of computation times
l in the absolute reference frame, introducing displacement and velocity functions, in addition
Notes:
l Where any force type loads are considered, the relative and absolute reference frames are the
same.
l Where participation factors has been computed in a given direction (DIRECTION label), only
results in the relative reference frame can be obtained.
Input files :
Output files:
DATA
TRANSIENT
These files correspond to a geometry reduced to the selected nodes and elements (the DATA
file only contains geometry and material properties required for analysis of the TRANSIENT
file).
Listing:
Optional printout of requested results for nodes and elements, and selected cards. Printout of
maximum values and occurrence times, where appropriate.
Input data
DYNAMIC
COMBINATIONS nv (RELATIVE) (FLUID)
TITLE
title
SELECT i1 TO i2 STEP i3
(FUNCTION)((i11, i12, i13) ...(ind1, ind2, ind3))
(TABLE)
iT/ ...
PRINT (MAXI) (MODULUS) (COMPONENT)
(RESULT) i1 TO i2 STEP i3 (INTERNAL)
CARD n1 TO n2 STEP n3 (STORE )
...
RETURN
RELATIVE used to obtain recombination in the relative reference frame (by default,
recombination in the absolute reference frame).
FLUID this label is optional if the fluid-structure coupling (MODE FLUID
STRUCTURE) has been taken into account when calculating eigenmodes.
Reduced pressures of the fluid are not recombined.
TITLE command
This optional command is used to associate a title (at most 40 characters) with the sum of loads
defined by FUNCTION command. By default the title of load 1 is kept.
SELECT command
Optional mode selection command.
Only modes specified in the list are used for recombination. If only one mode is selected, the
response on this mode is obtained.
The list must be a subset of the list defined in the MODAL ANALYSIS module.
FUNCTION command (optional)
If a MODAL ANALYSIS has been executed defining only the support acceleration functions,
displacement and velocity functions can be defined with this command, to obtain results in the
absolute reference frame. All TABLES must be redefined.
TABLE command (optional)
This command is present if the FUNCTION command is also present (see Chap. 11.4.3.3).
PRINT command
No printout or storage occurs if this command is absent.
If the MAXI label is present, maximum values for the results requested for any of the selected
cards are printed, with their corresponding occurrence times.
The COMPONENT and MODULUS labels can coexist, and one of these labels must at least
exist.
If the COMPONENT label is present, maximum values are computed component by
component.
If the MODULUS label is present, maximum modulus values are computed.
RESULTS command
These commands specify list of nodes or element used to define a reduced geometry. The
requested results are computed for a reduced geometry. The type of results requested must exist
in the RESULTS and MODES files. A number of RESULT commands can exist
simultaneously. No computation is done in the absence of a RESULTS command.
The types of result are selected by means of the following labels:
(RESULTS) = DISPLACEMENT list of nodes
REACTION list of nodes connected to the external medium,or belonging to a
STOP
The REACTION label is used to obtain, in addition, the force or forces exerted by the nodes on
the non-linearities (action of n1 on the exterior, or n1 on n2).
RESULTS can also be replaced by NODES, synonymous with DISPLACEMENTS and
REACTIONS (if DISPLACEMENTS is used, VELOCITIES and ACCELERATIONS are
automatically computed).
(RESULTS) = FORCE list of elements
(if FORCES the reduced geometry contains the nodes of the selected elements).
Forces in beams are computed in the element local axes. The local damping matrices are
neglected. If local reference frames are defined in the material properties, they are ignored for
this computation. .
For elements other than beams, forces are computed by combination of modal forces. In these
cases, it is supposed that the relation between displacement and acceleration is still available
mode by mode . Modal forces must exist in the MODES file (SOLVE FORCES
MODAL - 11.2.1.3.3). For all elements, using SOLVE FORCES TRANSIENT (12.3.5) after
DYNAMIC COMBINATION is used to obtain forces precisely, by neglecting however for
beams the loads of span, but taking into account concentrated dampings possibly defined in the
characteristics .
CARD command
Card selection command for result printout and storage.
n1 TO n2 STEP list of card numbers for printout or storage.
n3
STORAGE this label indicates that the command selects cards to be stored.
There can be one or two CARD commands.
No storage or printout is executed in the absence of a CARD command.
If the CARD command does not have the STORE label, printout and storage are done
simultaneously.
If the STORE label exists for one of the two CARD commands:
l printout is executed in accordance with the command not containing the STORE label.
l storage is executed in accordance with the command containing the STORE label.
Note
l In the case of variable step computations (see 11.4.3.3 ALGORITHM command), if n3 is too high,
the computation will stop at the last card for which the generalized coordinates were computed.
DEFINITION
MODAL ANALYSIS NON LINEAR * STOP IN MULTI-EXCITATION CASE
OPTION PLANE BEAMS
ONE-DIMENSION
NODES
1/ 1. 0.
21/ 1. 0.5081
ELEMENTS
1/ 1 21 20
MATERIALS
/E=2.07*11 NU=0.3 IZ=8.32*-7 S=3.87*-3 RO=7483
CONSTRAINT
1/UX UY RZ
LOAD
1 EXCITATION HORIZONTAL /
NODE 1/UX=1
RETURN
SOLVE STIFFNESS MASS
SAVE STIFFNESS MASS DATA 2
SOLVE
SAVE RESULT 2
DYNAMIC
MODES GIVENS
VECTORS 10
RETURN
SAVE MODES 2
SEARCH DATA STIFFNESS MASS MODES RESULT 2
DYNAMIC
PARTICIPATION 10 MULTI 1
MODAL ANALYSIS 10 MULTI 1
FUNCTION MULTI 1 29
TABLE
1/2, (0.001,0.,0.) (0.011,0.423333*-5,0.127*-2) (0.016, 1.164167*-5, 0.127*-2)
2/2, (0.001,0.,0.) (0.011,0.127*-2,0.254) (0.016,0.127*-2,-0.254)
9 / 1, (.001,0), (.011 0.254) (0.016,-0.254)
STOP
21 / UX MIN=-1.27*-7 MAX=1.27*-7 STIFFNESS 2.50*5 MULTI 1 2 9
INITIAL TIME 0.001
0.015 STEP 0.000125/DISPLACEMENT 200 STORE 1
RETURN
SAVE DATA MODAL 22
...
Example: (continued)
This problem is considered as a physical variable problem in Chapter 11.5 (absolute movement,
any damping matrix). With modal methods, generalized coordinates q k , associated with each
eigenmode, and describing relative movement with respect to the static deflection imposed by
the supports, verify the analytical equation:
Transfer Functions
Transfer functions H(ω) are basic tools for the study of random vibrations. They describe the
dynamic amplification of the response for each excitation frequency. These functions are used
for stochastic dynamic studies in particular (cf. Chapter 11.4.5 ). There are two equivalent
mathematical definitions:
l H(ω) is the ratio of the Fourier transform of the response to the Fourier transform of the
excitation:
o Computation of harmonic response combined with a frequency scan, provides the transfer
function. The modal method is consequently particularly suitable for their determination.
l H(ω) is also the Fourier transform of unit impulse response h(t) to a Dirac delta function at
t=0.
which can be computed analytically for each eigenmode, and is obtained in the physical space
by combining the different modal contributions.
Input Data
Input files:
DATA obtained from PARTICIPATION, containing participation factors pk
MODES obtained from computation of modes, containing modal forces if this
type of result is requested (SOLVE FORCES MODAL )
RESULT label:
At least one or more labels, with one per data line, selected from the following:
DISPLACEMENT VELOCITY ACCELERATION REACTION FORCE
A list of nodes (ni) or elements (mi) is associated with each type of result.
A new reduced geometry is defined from the various lists of nodes and elements, comprising:
l Nodes for which at least one result is explicitly requested,
l Nodes belonging to an element for which forces are requested,
Computation is then executed for each type of result explicitly requested, for all nodes or
elements of the reduced geometry.
Results
Files:
DATA reduced geometry at nodes and elements, selected by
DISPLACEMENT, VELOCITY, ACCELERATION, REACTIONS
or FORCE commands.
TRANSIENT result cards numbered sequentially according to ascending order of
frequencies, in modulus and phase form:
Modulus (positive) odd card
Phase (between -π and π) even card
Notes
l In contrast to HARMONIC RESPONSE computation (11.5), results are in the relative reference
frame, except in a case where the ABSOLUTE label can be used with the TRANSFER command.
l In the DATA and TRANSIENT files, the order of loads is redefined by the computation modules.
The loads corresponding to movements of the support are followed by force type loads.
l TRANSIENT DYNAMIC computation (11.7.4) can be used to combine phase-offset excitations.
Results are given in time-related form.
SOLVE FORCES TRANSIENT (12.3.5) after TRANSIENT DYNAMIC is used to obtain coherent
forces with displacement, velocity and acceleration fields. This forces can locally be different from
those directly obtained by combination of modal forces, especially in the proximity of concentrated
loads.
l Transfer function units depend on the type of result and the nature of load (response/excitation
ratio).
l The conversion of a transient file, outcoming from multi-excitation in modal dynamic response, by
using the command CONVERT is not possible directly. It is necessary to copy the cards of a local
TRANSIENT file in a RESULT file with the command MERGE CARD.
DEFINITION
COMPARE RHD RHM ABSOLUTE REFERENCE FRAME
OPTION PLANE BEAMS
ONE-DIMENSION
NODES
1 / 1. 0.
21 / 1. 0.5081
ELEMENTS
1 / 1 21 20
MATERIAL
/ E = 2.07 * 11 NU = 0.3 IZ = 8 . 32 * - 7 S = 3 . 87 * - 3 RHO = 7483
CONSTRAINT
1 / UX UY RZ
LOAD
1 /
NODE 1 / UX = 1
RETURN
SOLVE STIFFNESS mass
SAVE STIFFNESS mass DATA 2
DYNAMIC
HARMONIC RESPONSE
FREQUENCY
30 STEP 10 / DISPLA 1
RETURN
SAVE DATA TRANSIENT 10
Note
l The direct method is described in Chapter 11.5.
the circular frequency: . The zero order moment is used to evaluate the
mean value of the response.
l in the case of Gaussian signals, analysis of statistical moments and computation of material
property parameters: bandwidth, mean number of times a given level is exceeded, statistical
data on absolute and local maximum values, and fatigue analysis.
Before this module is used, the transfer functions must be computed by HARMONIC
RESPONSE MODAL (Chapter 11.4.4). Each of the three steps is covered by a separate
submodule.
Input data
DYNAMIC
PSD (CROSS) (LOG) (EDIT)
FUNCTION n1,...,nk
CALCULATE (list of result labels)
TABLE (FREQUENCY)
n /1,(........)
CARD na TO nb ... STEP nc
RETURN
SAVE DATA TRANSIENT n
PSD command
This command initiates computation of the response PSD for the structure.
The CROSS PSD label indicates integration of the real part of the cross power spectral densities
between the different loads, in the computation.
The LOG label selects bi- logarithmic interpolation in the load spectrum tables (linear
interpolation by default). In the presence of this label, tables are necessary type 1. (Chapter
4.15.3).
The EDIT label is used to edit the response power spectral density (PSD).
FUNCTION n1, ... ,nk command (mandatory)
Omitting cross PSD values, k is equal to j, the number of loads for which transfer functions
have been computed.
If cross PSD values are taken into account, k is equal to j+j.(j-1)/2
List n1,..,nk is used to give a list of the load spectrum tables (classed in succession by imposed
displacement and force on the free DOF, and then in ascending order of load numbers
introduced), followed by a list of cross PSD tables, classified coherently by load pair with the
load classification.
Example with 4 loads
FUNCTION sp1 sp2 sp3 sp4isp12 isp13 isp14 isp23 isp24 isp34
with
spi = PSD for load i
ispij = cross PSD between loads i and j
Input data
DYNAMIC
INTEGRATION
CARD n1 TO n2 STEP n3
CALCULATE (list of result labels)
MOMENT( m1,..,mj) (PRINT)
RETURN
SAVE DATA DISPLA VELOCITY ACCELERATION REACTION FORCE n
l The trapeze integration method is used. Frequency scan must be sufficiently fine in the vicinity of
a relatively undamped eigenfrequency.
l Integration is done between 0 and +∞, which explains why the odd order moments are non-zero.
DYNAMIC
STOCHASTIC n1 n2 n3
CALCULATE (list of result labels)
BANDWIDTH
LEVEL CROSSING j
u1,...,uj
MAXIMA k
r1,...,rk
FATIGUE COEF c EXPONENT a
RETURN
SAVE DATA DISPLA VELOCITY ACCEL REACTION FORCE n
STOCHASTIC command
This command initiates computation of statistical parameters.
n1, n2, n3 are load numbers corresponding to order 0, 2 and 4 moments for integration
computation. By default of every list, ni = i
CALCULATE command (mandatory)
This command specifies the type of result computed (list of possible labels:
DISPLACEMENT, VELOCITY, ACCELERATION, REACTION, FORCE).
DATA
RESULT
These two files contain as many loads as there are computed parameters, in the same order as
for printout.
The listing contains all requested results.
DEFINITION
STOCHASTIC - POUTRE CRANDALL ET MARK
OPTION PLANE BEAM
GEOMETRY
NODE
1/
2 / 0. 4.
ELEMENT
1/12
MATERIAL CARACTERISTICS
ELEM/E10.3*6AX 0.0625IZ3.23*-4VZVPZ0.125
CONSTRAINTS
N 1 / UX UY RZ
N 2 / UY
MASS
N 2 / AX 7.28*-4
RETURN
; eigen modes
DYNAMIC
MODES
VECTOR 1 FORCE NORM ELASTIC
RETURN
STRESSES
SEQUENCY BEAM
DYNAMIC
PARTI
HARMONIC RESPONSE MODAL DISPLA
DAMPING GAMMA 0.01
FREQUENCY
2000. STEP 10.
MODE 20 RATIO 0.1
TRANSFER STATIONNARY
DISPLA 2
FORCE 1
RETURN
DYNAMIC
PSD
FUNCTION 1 2 2
CALCULATE DISPLA FORCES
TABLE FREQUENCY
1/1(0.0.)(19.90.)(20.5930.)(2000.5930.)(2000.10.)
2 / 1 (0. 0.) (2000.1 0.)
CARDS 1 TO 241 STEP 1
RETURN
SAVE DATA TRANSIENT 2
; statistical moments
DYNAMIC
INTEGRATION
CARDS 1 TO 241 STEP 1
CALCULATE DISPLA FORCES
MOMENT
RETURN
SAVE DATA RESULTS 3
; displacement analysis
DYNAMIC
STOCHASTIC
CALCULATE DISPLA
BANDWITH
LEVEL 2
0. 0.001
MAXIMA 3
10. 100. 1000.
FATIGUECOEFFICIENT0.31EXPONENT6.09
RETURN
SAVE DATA RESULTS 4
END
of seismic excitation, essentially according to the RG 1.92 rev. 1 rules of the US Nuclear
Regulatory Commission. The ZPA method defined in RG 1.92 rev. 3 is not available. For
SYSTUS, amplified modes are combined with the usual combinations available (SRSS,
CQC etc.) up to the frequency cutting of the analysis, i.e. up to F point defined on RG 1.92
curves. For these modes, no distinction is performed between periodic and rigid components.
All modes with a frequency higher than the frequency cutting are taken into account by using
a pseudo-mode. The contribution of these rigid modes is evaluated by using the structure
static deflection when it is subjected to a uniform acceleration load. The overall response is
then derived by taking the square root of the sum of the square of amplified response and
rigid response.
Absolute movement U1 of the free nodes is broken down into two parts:
U1 = Ue + Ur
Ue : static driving movement resulting from displacement U2 of the supports.
Ur : relative movement with respect to the static deflection imposed by the supports.
U r depends on acceleration U 2 of the supports. This movement is studied by projection in the
modal base corresponding to U 2 = 0. This modal base can also include the residual modes
(static corrections).
Two types of analysis are possible:
l If all the supports have the same movement, the analysis is said to be MONO-SPECTRAL.
All nodes constrained by the external medium are assumed to possess this movement. This
concerns the non-modeled end of an elastic constraint. Only the relative movement is
studied, according to the pseudo-acceleration spectra.
l If the supports have different movements, the analysis is said to be MULTI-SPECTRAL.
Imposed displacement type loads are used to identify the different groups of supports. These
loads are the first of the static computation. Relative movement and driving movement are
studied according to the pseudo-acceleration spectra, and maximum support displacement
data.
Analysis Procedure
Spectral analysis assumes execution of a number of modules. Chaining of the modules differs
according to whether the excitation is mono- or multi-spectral.
Distinction is made between four main steps:
A - Static analysis (SOLVE)
printing a certain number of criteria (effective modal masses, the sum of which is compared
with the effective total mass with unconstrained degrees of freedom, for each excitation
analyzed).
l To compute static corrections which allow for modes omitted from the analysis, to ensure
correct representation of the static component of the seismic response (not amplified by
inertial effects).
D - Analysis by spectral combination
This final step provides the seismic response physical results. This response is not obtained in a
fully deterministic manner.
It is the result of the following operations:
l Computation of the response of each mode for a given direction or excitation (value read in
the spectrum, multiplied by the modal participation factor).
l Combination of modal responses for a given direction or excitation. These combinations are
Special Options
There are two special options for this type of analysis. The SHELL AXISYMMETRICAL
HARMONIC FLUID option embodies the particularities of both.
FLUID-STRUCTURE options
The mass required for computation of modal participation factors must take account of the
added mass due to the fluid. For execution of the PARTICIPATION module (Chapter 11.4.2),
the MASS file must be replaced by the MATRIX file, created after computation of eigenmodes,
and taking account of fluid-structure interaction.
As a static RESULT file is impossible, only mono-spectral analysis is possible without static
corrections.
Only the combination of fluid reduced pressures has a physical meaning. For the fluid the
reduced pressures associated with the walls movement are always obtained (processing of
displacements or accelerations).
AXISYMMETRICAL HARMONIC options (SHELL or THREE DIMENSIONAL)
Multi-spectral analysis is not possible.
Support excitations are axisymmetrical. Axis of symmetry OY must be vertical, but the position
of the horizontal seismic directions is not imposed.
A vertical translational movement, parallel to OY, is axisymmetrical, and corresponds to
HARMONIC 0. For a unit translation, the components at each point are as follows:
UX = 0 UY = 1 UZ = 0 (RZ = 0).
Only the axisymmetrical modes are loaded, or more precisely those which are symmetrical with
respect to any meridian plane, such that UZ = 0 at any point. Antisymmetrical torsion modes,
where only UZ is non-zero, can only be loaded by rotation round the axis of symmetry.
A unit translation in direction X corresponds to a radial displacement equal to cos θ, and a
circumferential displacement equal to -sin θ, namely HARMONIC 1 SYMMETRICAL, with
the following support excitation components:
UX = 1 UY = 0 UZ = - 1 (RZ = 0)
Only the modes of this harmonic are loaded.
A unit translation in direction Z, deduced from the previous translation by rotating through
+π/2, obviously corresponds to HARMONIC 1 ANTISYMMETRICAL, and only loads the
modes of this harmonic, deduced from those of the symmetrical harmonic (rotation through
π/2). Radial displacement is equal to sin θ, and circumferential displacement to cos θ, giving
the components:
MONO-SPECTRAL Analysis
All supports have identical movements. Loads corresponding to support movements (unit
accelerations) are generated by the PARTICIPATION module.
Preliminary steps
l Step 0: static DATA file
constraints with the external medium
unit acceleration loads (GX, GY, GZ) if static corrections.
l Step 1: static results (SOLVE) if static corrections
l Step 2: modal analysis (Chapter 11.2)
DYNAMIC MODES
(SOLVE FORCES MODAL: if forces and accelerations).
l Step 3: participation factors (cf 11.4.2.3)
Note
l If DIRECTION and RESIDUAL-TERMS GN are used, the static load for step 0 must be in this
direction, by combination of the elementary accelerations defined by GX, GY, GZ
(for example GX 0.707 GY - 0.707).
Input data
Input files:
DATA from PARTICIPATION
MODES from MODES computation without static corrections
from PARTICIPATION with static corrections.
Input data
l FILE nf: dampings are searched in DAMPnf.DAT file created by DAMPING (11.7.1).
The labels of the DIRECTION command describe the methods for response combination in
each excitation direction. At least one of the following labels must be given; otherwise the label
by default is taken into account.
SRSS (by default)
ABSOLUTE
SEMI-SRSS(in this case, there must be 3 directions for the label to be meaningful).
TABLE command
Tables ix, iy and iz defining the pseudo-acceleration spectra:
l are type 1 (Chapter 4.15.3). If the first number of the table is different from 1, the pseudo-
acceleration values are multiplicated by this number mandatory positive (case of the table iy).
l comprise pairs (frequency in Hz, acceleration in units of acceleration),
It is necessary to use this command in the absence of FILE label on SPECTRUM command.
Results
Files:
DATA
the number of loads depends on selected combination methods
RESULTS
In the presence of the EDIT label, the listing contains all results, if not it contains only the load
titles generated by the analysis.
The quantities obtained are relative movement amplitudes.
Accelerations represent a special case. With the inclusion of static corrections, a good
approximation of absolute accelerations is obtained. Otherwise these are relative accelerations.
The acceleration storage must be explicitly requested (SAVE DATA RESULTS
ACCELERATION i commands).
The reactions at nodes linked to a local frame are also a particular case (CONSTRAINT
command), they are defined in the local frame of the node.
DEFINITION
MONO-SPECTRAL ANALYSIS WITH STATIC CORRECTIONS
OPTION SPATIAL BEAMS
ONE-DIMENSION
NODE
1/
11/ 10. 0. 0.
ELEMENT
1/ 1 11 10
MATERIALS
/E=2*11NU=0.3 RHO 7800 S 0.025 IX IY IZ 0.0002
CONSTRAINT
N 1 11 / UX UY UZ RX RY
LOAD
1 ACCELERATION GY / GY 1.
2 ACCELERATION GZ / GZ 1.
RETURN
1 / 1 0. 1 1. 3. 1. 5. 10. 25. 10. 30. 1. 100. 1.
2 / 1 0. 1 2. 3. 2. 6. 15. 30. 15. 33. 1.5 100. 1.5
RETURN
SAVE DATA RESULTS 1100
END
MULTI-SPECTRAL Analysis
The different sets of supports have different movements, and are identified by means of
imposed displacement type loads.
The main steps of the analysis are summarized in Table 11. 5.
Preliminary steps
l Step 0: static DATA file
o Constraints with the external medium. If local reference frames are defined for the nodes,
support movements are expressed in these reference frames. Generally, only identical
reference frames for all nodes are significant in a seismic context, and static corrections are
only strictly valid in this case.
o Imposed displacement loads defining the amplitude and direction of the excitation (in the
DYNAMIC
PARTICIPATION nv
RESIDUAL-TERMSGX nx GY ny GZnz MODES mc--
(FREQUENCY fc) <LABEL> (EQUIVALENT) (PRINT)
PARTICIPATION nv MULTI m
RETURN
Note
l If static corrections are not included, only the PARTICIPATION nv MULTI m command is
necessary.
Input data
Input files:
DATA from PARTICIPATION MULTI
RESULTS from SOLVE (static computation)
MODES from MODES no static corrections
from PARTICIPATION with static corrections
Input data
DYNAMIC
SPECTRAL ANALYSIS MULTI GX n1x ... nix GY n1y ... niy GZ n1z ... niz--
(RESIDUAL-TERMS) (EDIT)
COMBINE (LABELS)
STATIC (LABELS)
DISPLACEMENTS U1, U2, ..., Um
MODE n (FREQUENCY fc) (LABELS)
DURATION t
DAMPING (GAMMA γ ) (MODAL) (FILE nf)
j1 TO j2 STEP j3 / (GAMMA) γj
SPECTRUM i1 i2, ..., im(FILE ACCEL g1 g2 ... gm)
DIRECTION (LABELS)
TABLES
i1/1 (fn, San)
:
im/1 "
RETURN
The COMBINE, MODE, DURATION, DAMPING and TABLE commands are identical to
those for MONO-SPECTRAL analysis.
SPECTRAL ANALYSIS MULTIPLE command
GX n1x ... nix, GY n1y ... niy GZ n1z ... niz
nix, niy, niz are static load numbers (imposed displacements associated
with each direction X, Y, Z).
RESIDUAL- inclusion of static corrections
TERMS
The pseudo-acceleration associated with each residual mode is the last maximal table value
describing the excitation spectrum in its direction. This value is also used to correct reactions at
supports.
EDIT Complete printing of the results.
STATIC command:
This command is used to describe methods for combination of static results with the following
labels:
l SRSS (by default)
l ABSOLUTE
l MIXED
At least one of these labels must be given; otherwise the label by default is taken into account.
In the lack of this command, only the relative movement with respect to the static deflection
imposed by the supports is determinated.
DISPLACEMENT command:
U1, U2, ..., Um are multiplier factors, used to obtain maximum values for support physical
displacements from the displacement values specified in the m imposed displacement type static
loads (determination of driving movement due to the supports).
If the STATIC command is existed, it must be immediately followed by the DISPLACEMENT
command.
SPECTRUM command:
This command associates the table numbers with each imposed displacement type static load in
load order.
In the absence of FILE label, these numbers are table numbers defined by the TABLE
command. In the presence of RESIDUAL-TERMS label, the acceleration linked with static
correction is the last value of the table.
In the presence of FILE label, these numbers define SPECix.DAT files.
These files are of the same type of those created by GENERATION SPECTRA (Chapter
11.7.2).
Ground maximal accelerations have to be defined with the help of the ACCELERATION
label. They are used to correct reactions at support levels and in the presence of RESIDUAL-
TERMS label are associated with static correction. Dampings have to be defined for each
mode. PRINT label brings about the printing of spectra actually taken into account.
DIRECTION command
This command enables the user to manage the combination of the excitations due to different
supports in one direction and the combination of the basic directions (X, Y, Z)
Groups of supports in one direction. Two methods exist:
l Root-sum squares: QUADRATIC and SEMI-QUADRATIC labels,
l Simple summation of the response of each support: ABSOLUTE, MIX labels. At the end of
this stage, we have the dynamic responses for each direction (X, Y, Z). In the presence of the
CORRECTION label, the unit static correction (calculated in the CORRECTION module)
for the considered direction, is summed in a quadratic manner to the dynamic response. The
acceleration taken into account when calculating the static correction is the maximum of the
last values of the pseudo-acceleration spectra in each of the basic directions (X, Y, Z). If
pseudo-acceleration spectra are given on file, the value of the acceleration taken into account
in each direction must be given.
Summation of the basic directions
The end of the procedure is carried out by taking into account the response in each elementary
direction and taking into account one of the four possible labels: QUADRATIC, SEMI-
QUADRATIC, ABSOLUTE, MIX
At this stage of the calculation, we have relative responses compared to the static deformation
imposed by the absolute displacements of the various groups of supports. The quadratic sum
(mandatory) of the static deformation due to the supports displacements (Commands STATICS,
DISPLACEMENT) and of the relative dynamic responses due to the supports accelerations
enable to the obtain the final result.
Notes:
l If a stress analysis (or forces in continuous media with the bending stresses method, especially
used in the nuclear domain, has been carried out :
l The driving movement creates secondary stresses (Q),
l The relative movement with respect to the static deflection creates primary stresses (P).
l This remark justifies the distinction between different kinds of results.
l It is obvious that total stresses must be considered as primary stresses.
Results
Files:
DATA
the number of loads depends on the selected combination method
RESULTS
In the presence of the EDIT label, the listing contains all results, if not it contains only the load
titles generated by the analysis.
Quantity amplitudes are obtained for the following:
l driving movement (static),
l relative movement with respect to the static deflection imposed by the supports,
l absolute movement.
Accelerations represent a special case (no driving acceleration). With the inclusion of static
corrections, a good approximation of absolute accelerations is obtained, otherwise these are
relative accelerations. The accelerations storage must be explicitly requested (command SAVE
DATA RESULTS ACCELERATION i).
The reactions at nodes linked to a local frame are also a particular case (CONSTRAINT
command), they are defined in the local frame of the node.
Example: Multi-spectral analysis with static corrections
DEFINITION
MULTI-SPECTRAL ANALYSIS WITH STATIC CORRECTIONS
OPTION SPATIAL BEAMS
ONE-DIMENSION
NODE
1/
11/10.0.0.
ELEMENTS
1/1 1110
MATERIALS
/E=2*11NU=0.3 R0 7800 S 0.025 IX IY IZ 0.0002
CONSTRAINT
N 111 / UX UYUZRXRY
LOAD
1 NODE 1 - UY=1/ NONE
N 1 / UY 1.
2 NODE 1 - UZ=1/ NONE
N 1 / UZ 1.
3 NODE 11 - UY=1/ NONE
N 11 / UY 1.
4 NODE 11 - UZ=1/ NONE
N 11 / UZ 1.
5 ACCELERATION GY / GY 1. NONE
6 ACCELERATION GZ / GZ 1. NONE
RETURN
SOLVE
SAVE DATA RESULTS 9000
SOLVE MASS
SAVE STIFFNESS MASS 9000
DYNAMIC
MODES GIVENS
VECTORS 6 PRINT -1
RETURN
Example: (continued)
Transient Response
This modulus solves in time the modal equation of an elastic structure:
* The TABLE command must appear only once and defines all tables useful for the modulus.
PSEUDO MODAL TRANSIENT
initiates execution of the modulus.
ELEMENTS, optional commands of structure modification (§ Chapter 11.6.3).
TABLES
ALGORITHM Newmark’s parameter definition.
or α (default 0.05)
Harmonic Response
This modulus is used to solve, for a circular frequency or frequency given, the modal equation
of an elastic structure:
Input files:
OMEGA label
PSEUDO MODAL HARMONIC
initiates the module execution.
ELEMENTS, optional commands of the structure modification (Chapter 11.6.3).
TABLES
FREQUENCY (INITIAL f0)
defines the frequency scan.
This command follows necessary the commands of the structure
modification if they are existing.
MODE Command
Definition of the scan step in the vicinity of the eigenmodes contained in the minimum scan
window defined by the FREQUENCY command.
If fm is the eigenfrequency computed, the (2nm+1) frequencies obtained by dividing the interval
[ fm(1-r), fm(1+r) ] into 2 nm equal intervals, are added to the minimum scan window. If r = 0,
only eigenfrequencies processed are added. PRECISION command is used to suppress
additional frequencies which are too close (by default ε = 10-5).
A solution in the form is obtained by solving a linear system for each value
of circular frequency ω.
No damping (K - Mω2) X = F Y=0
The solution is infinite for the circular eigenfrequencies of the system. A real component, with
its sign, is obtained.
With damping
Notes
l Only one damping file exists as input file of this procedure. C0 is of viscous type (GENERALIZED
label) or structural (STRUCTURAL label).
l With usual notations for one degree of freedom systems:
DATA file created by DEFINITION, which must contain the load or loads for
which harmonic response is computed.
STIFFNESS files created from the DATA file by the command. SOLVE
MASS STIFFNESS MASS (DAMPING).
DAMPING This file is only required if GENERALIZED or STRUCTURAL type
damping is applied.
MODE This file is only required in the presence of the MODE command.
Input data
DYNAMIC
HARMONIC RESPONSE (VELOCITY) (ACCELERATION)(REACTION)
DAMPING (ALPHA α) (BETA β) (GAMMA γ) ((GENERALIZED) (STRUCTURAL))
(METHOD <label> ( COMPLEX))
CIRCULAR FREQUENCY INITIALw0
FREQUENCY
ω1 STEP dω1 / DISPLA n
ω2 STEP dω2/ DISPLA n
...
MODES nm (RATIO r) (PRECISION ε) (COMPLEX)
RETURN
RESPONSE command
If the VELOCITY and ACCELERATION labels are present, velocities and accelerations are
computed in amplitude and phase form as displacements.
In the presence of the REACTION label, reactions and residual in forces are computed. The
maximal residual amplitude force is printed.
DAMPING command
If the GENERALIZED or STRUCTURAL labels are present, damping Co is extracted from
the DAMPING file. The various damping forms are added. If STRUCTURAL label is present
Co is of structural type.
METHOD Command: optional
This command is used to choose solver type. In the presence of damping, a complex solver is
required.
<label>
without damping SKYLINE, OPTIMISED
with damping SKYLINE, SPARSE, ITERATIV (COMPLEX)
For more information about this command and the properties of the various solvers available see
Chapter 21.
By default, standard solvers are used.
CIRCULAR FREQUENCY or FREQUENCY command
Computation circular frequencies or frequencies are introduced after the command of the same
name, followed by the initial bound of the interval (default value 0). Each interval (ωi-1 , ωi ) is
then subdivided into an integer number of subintervals, each with length dωi.
The DISPLA n parameter forces printout of moduli and phases every n intervals. In the
absence of this label there is no printout. Cards are stored on TRANSIENT file at each
subinterval.
MODE Command: optional
In the presence of a MODE file, definition of the step of balayage near the eigenfrequencies
included in the minimum balayage window defined by the FREQUENCY command. In the
lack of damping, this command is ineffective.
If fm is the treated eigenfrequencies, the (2nm+1) frequencies obtained by dividing the interval
[fm(1-r), fm(1+r)] into 2 nm equal intervals are added to the minimum balayage0 If r is null, only
the treated eigenfrequencies are added. The PRECISION label is used to suppress the
additional frequencies too close (by default ε = 10-5).
If the MODE file contains complex eigenmodes (modulus SUBSPACE BLOC DAMP or
PSEUDO MODAL DAMP), the COMPLEX label is required. If the system is enough
damped, the insertion of complex eigenfrequencies is better.
11.5.3 Results
Results (displacement, velocity, acceleration) are obtained in the absolute reference frame.
Computation of reactions and forces is not available under complex form.
Output files:
DATA unmodified
TRANSIENT only complex results are stored in this file, in the form:
. with damping: modulus and phase,
. no damping: real component.
MODULUS (positive) or real component odd card
PHASE (between - π and π) even card
Notes
l TRANSIENT DYNAMIC (Chapter 11.7.4) can be used to obtain displacements, velocities and
accelerations under time-related form, with eventually combination of out of phase loads .
l SOLVE FORCE TRANSIENT can be used to associate dynamic forces to time-related response.
l This command is also used to associate with the dynamic forces (real or complex) the harmonic
response (real or complex).
Example 1:
Viscous damping
Example 2:
Viscous and structural damping
Example 3:
Generalized damping using iterative method and force calculation in a magnetodynamic
problem.
The TRANSIENT file number 2 contains the results stored in amplitude and phase while
TRANSIENT file number 3 contains results stored in real or imaginary part.
Example 4:
Harmonic response in a thermal problem.
SEARCH DATA 1
SOLVE STIFFNESS MASS
DYNAMIC
HARMONIC RESPONSE
DAMPING GENERALIZED
(METHOD COMPLEX <LABEL>)
CIRCULAR FREQUENCY
:
RETURN
<LABEL> SKYLINE
DIRECT (SPARSE)
ITERATIVE (SPARSE)
(for more precision, see chapter 21).
Input Data
Input files:
DATA
for the unmodified model
MODES
Input data
DYNAMIC
SENSITIVITY nv1 VECTOR nv2
MODIFICATION (list of modifications)
NODE (list of nodes)
ELEMENT (list of elements)
COUPLE (list of pairs of nodes)
RESPONSE
DISPLACEMENTS (list of nodes)
AXES ref (RESPONSE)(MODIFICATION)
RETURN
SENSITIVITY command
nv1 number of modes for which sensitivities are computed,
nv2 number of modes of the modal base in which computation is performed
(nv1 ≤ nv2).
MODIFICATION command
(list of list of labels indicating the type of modification applied to the structure
modifications)
AX, AY, AZ addition of mass in directions x, y, z
IX, IY, IZ addition of rotational inertia about x, y, z axes
This command must be followed by one or more of the NODE, ELEMENT or COUPLE
subcommands (COUPLE subcommands must appear after NODE or ELEMENT
subcommands).
NODE subcommand
List of nodes on which modifications are applied
Example:
NODE 1, 5, 10 TO 15 INTERNAL
ELEMENT subcommand
List of elements on which modifications are applied.
These are two-node elements.
Addition of a mass value between two nodes has no meaning.
Example:
ELEMENT 8, 4, 12 TO 18
COUPLE subcommand
List of pairs of nodes (not connected by existing elements), between which changes are applied.
Addition of a mass value between two nodes has no meaning.
Example:
COUPLE (4, 5), (10, 20)
RESPONSE command
This command initiates computation of sensitivities on mode shape. The list of nodes for the
DISPLACEMENT subcommand specifies where the eigendisplacement modifications are
computed.
AXES command
This command is used to define the reference frame, by means of its Bryan angles (Chapter
4.4.4), in which:
l modifications are defined (MODIFICATION label),
Results
All results appear in the listing. No file is created.
Sensitivity of circular eigenfrequencies (or eigenfrequencies)
Input Data
Input data files required:
DATA
MODES (real eigendisplacements)
DAMPING
STIFFNESS required for computation of modal matrices, reactions, forces in beams and
MASS Rayleigh quotients.
Input data
DYNAMIC
PSEUDO MODALnv(DAMPING ONLY)(ORTHOGONALITY)--
(NORM <LABEL>) (EXTRACT)
SELECT na TO nb STEP nc
ELEMENT
ndTO neSTEPnf(INTERNAL)/TABLEt3t2t1(ELEMENTn)
ng / TABLE ...
TABLE
t1 / (mnij, Kmnij), etc.
t2 / (mnij, Mmnij), etc.
t3 / (mnij, Cmnij), etc.
VECTORS ns (LABEL) (PRINT i) (RAYLEIGH)
RETURN
ORTHOGONALITY
The stiffness and mass modal matrices are stored, and used by default.
This label (optional) provides for recomputation from the STIFFNESS,
MASS and MODE files.
NORM This label (optional) indicates the MODE file real eigenmode norm to
be used: MAX, MASS or ELASTIC. The ELASTIC norm is used by
default.
EXTRACT This label (optional) provides for printout of the following:
-stiffness mass and damping modal matrices,
-global matrix of the system to be analyzed,
-eigenmodes obtained in modal space.
These data are printed in addition to circular eigenfrequencies, critical
damping fractions, logarithmic decrements and physical modal
deflections called for by the VECTOR command.
SELECT command (optional)
na TO nb STEP nc
selects the modes effectively taken into account. If omitted, the first nv modes are taken into
account.
ELEMENT command (optional)
This command indicates the elements to which modifications are applied. The modifications
are characterized by stiffness, mass and/or damping matrices which can be symmetrical or not.
nd, ne, nf, ng create a list of the numbers of these elements (previously defined with
DEFINITION).
Notes
l the elements concerned do not necessarily have X9XX codes,
l no reference frame can be defined with the ELEMENT command,
l matrices which can be non-symmetrical must be given in full, and are added to those contained in
the STIFFNESS, MASS and DAMPING files.
l non-symmetrical stiffness, mass and damping matrices are given with the TABLE procedure if
t3t2t1 >0 (cf. 7.8.1.2 a), or in files t3t2t1 <0 (Chapter 11.6.3).
TABLE command (optional)
This command defines the stiffness, mass and damping matrices in the way described in the
ELEMENT command.
Results
Output files:
DATA contains the titles of loads corresponding to the eigenmodes
RESULTS contains MAX norm eigenmodes
MODE contains MAX norm eigenmodes
Notes
l Only eigenfrequencies with a positive or zero imaginary part are retained, together with the
corresponding eigenvectors.
l For each complex eigenmode stored there are two cards in the MODE file (or two loads in the
RESULT file), such that:
l the real part corresponds to an odd number
l the imaginary part corresponds to the next even number.
l In the special case where the structure is not damped, and there are no elements which are non-
symmetrical in terms of stiffness or mass, the modes are real and stored by card.
l1 = 7 + 2*ls
o ne: super-element number
n0=0
o ls: dimension of matrix s
l the stiffness matrix associated with an element with nb nodes comprises *nb blocks, each
with kp*kp terms, where kp is the number of dof. Storage is by column, both for blocks and
the terms of each block.
Storage example for a matrix with (2 nodes and 3 dof):
ls=36
1 4 7 19 22 25
2 5 8 20 23 26
3 6 9 21 24 27
10 13 16 28 31 34
11 14 17 29 32 35
12 15 18 30 33 36
An identity matrix (2 nodes, 3 dof) can be written to the IDENTITY.DAT file in the following
form:
.
Examples
Example 1:
Computation of complex circular eigenfrequencies and eigenmodes for a plane beam. Damping
is symmetrical, and can be contained in the DAMPING file or entered directly using the
PSEUDO MODAL procedure.
DEFINITION
COMPLEX FREQUENCIES OF A DAMPED BEAM
OPTION PLANE BEAM
ONE-DIMENSIONAL
NODES
1 / 1. 0.
21 / 1. 0.5081
ELEMENTS
1 / 1 21 20
MATERIAL PROPERTIES
/ E 2.07*11 NU 0.3 AX 3.87*-3 IZ 8.32*-7 RHO 7483.
CONSTRAINTS
N 1 / UX UY RZ
LOAD
1
DAMPING
1 21 / AX AY IZ 1.
2 TO 20 / AX AY IZ 2.
VERIFY
RETURN
SAVE DATA 2
; direct computation
SEARCH DATA 2
DYNAMIC
MODE SUBSPACE BLOCK DAMPING
VECTORS 4
RETURN
Example 1 (contd.)
; pseudo-modal computation
SEARCH DATA 2
DYNAMIC
MODE SUBSPACE
VECTORS 100 STURM FREQUENCY 20000
RETURN
SAVE MODE 2
DEASSIGN 33 ;necessarytoavoidthesumofthetwotypesofdamping
Example 2
Computation of complex frequencies for a structure with non-symmetrical stiffness matrix.
DEFINITION
Complexfrequencies ofa structure withnon-symmetricalstiffness matrix
OPTION PLANE-BEAM
GEOMETRY
NODES
1/
2 / 1.
ELEMENTS
1 / 1 2 TYPE 1902
MATERIAL PROPERTIES
E 1 / TABLE 1
MASS
N 1 2 / AX 1.
TABLE
1 / 1111 2. 2211 2. 1211 1.
CONSTRAINTS
N 1 2 / UY RZ
LOAD
1
VERIFY
RETURN
SAVE DATA 5
characterized by its critical damping ratio ξi (ratio of damping to critical damping). Two
equivalent dampings, associated with E or T, are then computed for each mode k:
o If some dampings are defined in element material properties by ALPHA, BETA, EPSX
Input Data
Output files:
STIFFNESS, used for computation of eigenmodes
MASS
DATA, MODES obtained from computation of eigenmodes
Input data:
DYNAMIC
DAMPING nv (FILE nf (LOWER inf)(UPPER sup))
(SELECT i1 TO i2 STEP i3)
ELEMENTS j1 TO j2 STEP j3 (INTERNAL) / XSIE xsie
ELEMENTS k1 TO k2 STEP k3 (INTERNAL) / XSIE xsie XSIH xsih--
XSIV xsiv CIRCULAR p
ELEMENTS ...
...
RETURN
DAMPING Command nv
The critical damping ratios corresponding to the first nv eigenmodes are computed (nv ≤ 640).
If nv is absent, it is taken equal to the number of eigenmodes of the MODE file.
nv can although be defined by FREQUENCY fo (or CIRCULAR FREQUENCY ωo)
In the presence of FILE label, weighted dampings are saved in DAMPnf.DAT file (XSIE and
XSIT quantities). LOWER-UPPER optional labels can impose lower and upper values of
weighted equivalent dampings.
SELECT i1 TO i2 STEP i3 command optional
Mode selection. Computation is performed only for the modes selected from the first nv modes.
The eigenmodes of rigid bodies (null frequency) must in priori be suppressed.
ELEMENTS i1 TO i2 STEP i3 / XSIE xsie XSIH xsih XSIV xsiv CIRCULAR p
command
This command is used to define a component and its damping parameters. A component can be
defined with a number of ELEMENT commands, if necessary:
XSIE xsie critical damping ratio for component i (method 1)
XSIH xsih coefficient ξH - equation 3 (method 2)
Notes:
l Damping data are substitutive for the same element.
l The stiffness of elastic constraints is not taken into account for computation of strain energies (use
equivalent stiffness elements). Local masses are not taken into account for computation of kinetic
energies (use equivalent mass elements).
Results
All results appear in the listing. The critical damping ratio value is obtained for each
eigenmode.
Damping values defined by method 1 give XSIE1 (eq.1) and XSIT (eq.2). The user must
select weighting of strain energies or kinetic energies. Damping values defined by method 2
(XSIH, XSIV) give XSIVH (eq.4). XSIE column contains the summation of dampings
weighted by strain energies (XSIE1+XSIVH).
Example:
Reference:
1 - B.V. WINKEL and L.J. JULIK
Experience with non-uniform damping in the seismic analysis of nuclear plant components.
ASME Pressure Vessel and Piping Conference - June 19-23, 1983
Input Data
Accelerogramm defining by a table
Input file:
DATA A data file must exist and can contain the TABLES defining the
accelerograms (Chapter 4.15).
Input data
GENERATION command
the PRINT -1 label cancels printout.
ACCELEROGRAM command
i1 / TABLE nt TIME (INITIAL ti) FINAL tf STEP step (FILE nf)
i1 accelerogram number
nt1 number of TABLE containing the load 1.
ti initial time (0 by default).
tf final time.
step time step defined over [ti,tf]
The time step must be small enough for linear description of the
accelerogram by segments (mandatory parameter).
nf results are stored in SPECnf.DAT (20 by default).
i2 /.. other accelerogram (nf different)
DAMPING γ1, γ2,... γn command
Critical damping ratio values are introduced in sequence, with n < 10.
METHOD command
Selection of the integration method (DUHAMEL by default).
FREQUENCY INITIAL f0 command
This command defines frequency or circular frequency scan ranges
f0 : scan start frequency (f0 > 0)
One or more fp STEP dfp commands are used to define computation steps. The fp STEP dfp
card breaks down the interval into (fp - (fp - 1))/dfp intervals.
The DISPLA np label forces printout every np subintervals. By default, printout occurs at the
end of each interval. Card storage occurs on each subinterval.
Results
Output files
SPECnf.DAT (one file per accelerogram)
This sequential binary file contains the following for each critical damping ratio value,
according to frequency:
l pseudo-acceleration spectrum,
l pseudo-velocity spectrum,
l acceleration spectrum,
Results can be plotted using the CURVE FUNCTION SPECTRUM routine (see 14.8). The
abscissas of the curves are frequencies.
SPECnf.DAT File
The SPECnf.DAT file is a function tabulated on a sequential binary type file (see 4.15.4.1);
It can be directly used by SPECTRAL ANALYSIS (cf 11.4.6 SPECTRUM FILE command).
Only pseudo-accelerations are used. Critical damping fraction values γ 1 , γ 2 , γ n must at this
moment form a rising monotonous series.
It is organized as follows:
l 1 record containing NGAMMA + 3 words:
For each critical damping ratio value, four cards are stored ( SD, SV, SA and Tmax ), with
each card containing a pair of values:
Card1 to NGAMMA: frequency and associated pseudo-acceleration spectrum
Card NGAMMA+1 to 2 NGAMMA : frequency and velocity spectrum
Card 2 NGAMMA + 1 to 3 NGAMMA : frequency and acceleration spectrum
Card 3 NGAMMA +1 to 4 NGAMMA : frequency and time corresponding to maximum
displacement value.
The SPECnf.DAT file can be reread or created by a program of the following type:
;
; RESPONSE SPECTRUM AND GRAPHIC VERIFICATION
;
DEFINITION
HORIZONTAL AND VERTICAL RESPONSE SPECTRUM
OPTION TWO-DIMENSIONNAL
GEOMETRY
NODE
1/
TABLE
; TEMPORAL ABSOLUTE ACCELERATION
1 / 6 (2. 11) (5. 12) (3. 13) ; HORIZONTAL ACCELERATION
2 / 6 (1. 11) (2.2 12) (1.8 13); VERTICAL ACCELERATION
11/ 3 31.416 -1.5708 0. 20. ; SINUS - FREQUENCY 5 HERTZ
12 / 3 94.248 -1.5708 0. 20. ; SINUS - FREQUENCY 15 HERTZ
13 / 3 188.496 -1.5708 0. 20. ; SINUS - FREQUENCY 30 HERTZ
; GRAPHIC VERIFICATION OF RESPONSE SPECTRUM
31 / FILE 81 CARD 1 ; HORIZ.- GAMMA=0.02
32 / FILE 81 CARD 2 ; HORIZ.- GAMMA=0.04
33 / FILE 81 CARD 3 ; HORIZ.- GAMMA=0.06
41 / FILE 82 CARD 1 ; VERTI.- GAMMA=0.02
42 / FILE 82 CARD 2 ; VERTI.- GAMMA=0.04
43 / FILE 82 CARD 3 ; VERTI.- GAMMA=0.06
VERIFY
RETURN
SAVE DATA 1
DYNAMIC
GENERATION SPECTRUM
ACCELEROGRAM
1 / TABLE 1 TIME INIT 0. FINAL 20. STEP 0.0005 FILE 81
2 / TABLE 2 TIME INIT 0. FINAL 20. STEP 0.0005 FILE 82
DAMPING 0.02 0.04 0.06
FREQUENCY INITIAL 1.
50. STEP 1.
RETURN
The DATA and TRANSIENT files are taken from MODAL ANALYSIS. They contain 3
nodes in PLANE BEAM option. 6 spectra files will be generated by the computation (2
translation ddl by node). The files will have the number 61 to 66 (furthermore, a file
TABL99.DAT will be created).
l Unsupported mass,
Participation factors, and effective masses and inertia are computed for all degrees of freedom of
the option.
This procedure is a complement to the PARTICIPATION module, computing effective inertia
with respect to a given center of rotation. Modal response computation is not possible
following EFFECTIVE MASS.
Input Data
Input files:
DATA
MASS
MODES
There is no output file. Results appear in the listing.
DYNAMIC
EFFECTIVE MASS nv (NODE I1) (POINT x y z)
RETURN
In the lack of centre of rotation definition, nv is by default, equal to the eigenmodes number of
the MODE file.
Example
We obtain displacements (deformed mesh), velocities and accelerations. Forces in elements can
be obtained by SOLVE FORCE TRANSIENT procedure.
Mathematical formulation
. n elementary loads Fi (Ui) (DEFINITION)
. m circular frequencies (HARMONIC RESPONSE)
We obtain the response to a load of form
Input Data
Input files:
DATA a file contains:
-geometry (reduced)
-number of loads of the harmonic response computation
-(tables describing the evolution of loads amplitude and phase with circular
frequency or frequency)
TRANSIENT this file contains complex displacements stored in modulus and phase
form. It is created from:
HARMONIC RESPONSE with (see Chapter 11.5)
HARMONIC RESPONSE MODAL (see Chapter 11.4.4)
Input data:
TRANSIENT DYNAMIC
TITLE
…
PHASE(list)(DEGREES)(VARIABLE(FREQUENCY))
AMPLITUDE (list) (VARIABLE (FREQUENCY))
TABLE
t1/
:
tn/
CARD (list) (COMPLEX)
TIME INITIAL t 0
t 1 STEP dt 1
-----
TIME INITIAL t0
t1 STEP dt1
----
tn STEP dtn
PERIOD p NSTEP n
Results
Output files:
DATA modified file containing one load only,
TRANSIENT file containing displacements, velocities and accelerations for all times
requested.
Using the DATA and TRANSIENT files supplied by TRANSIENT DYNAMIC, the
procedure can be used to compute forces at each time (see Chapter 12.3.5).
Examples
Embedded-free beam under two harmonic forces.
l Transfer functions for a reduced geometry comprising nodes 1, 2, 6 and 11 and element 1.
DEFINITION
MULTIPLES AND PHASE-SHIFTED EXCITATIONS
OPTION PLANE BEAM
ONE-DIMENSIONAL
NODES
1/
11 / 2.
ELEMENTS
1 / 1 11 10
MATERIAL PROPERTIES
E 1 TO 10 / E 2.*11 NU 0.3 RHO 7800. AX 0.0005 IZ 2.6*-8
CONSTRAINTS
N 1 / UX UY RZ
LOADS
1 FORCE AT NODE 6 /
N 6 / FY 1000.
2 FORCE AT NODE 11 /
N 11 / FY 1000.
VERIFY
RETURN
SAVE DATA 1
SOLVE
SAVE RESULTS 1
DYNAMIC
MODES GIVENS
VECTORS 10 FORCES
RETURN
SAVE MODES 1
SEARCH DATA RESULTS 1
DYNAMIC
PARTICIPATION 10 FORCE 2
HARMONIC RESPONSE MODAL 10 FORCE 2
DAMPING GAMMA 0.04
FREQUENCY
5. STEP 5. / DISPL 1
30. STEP 25. / DISPL 1
TRANSFER STATIONARY
DISPL 6 11
FORCE 1
RETURN
SAVE DATA TRANSIENT 10
F6 = 1000 cos 10 πt
Note
l When TABLES have been defined in HARMONIC RESPONSE MODAL, TABLES associated to
VARIABLE label must be compulsorily defined on the reduced geometry.
CHAPTER 12 TRANSIENT
12.1 OVERVIEW
12.1.1 Introduction
In this chapter the linear and non-linear analyses of the response of structures to time-dependent
loads is described. The general equations considered are those for heat conduction,
K. θ + C . = Q (t)
where θ and are temperature and the first derivative of temperature with respect to time.
and for the dynamic analysis of structures,
KU+C. +M. = F (t))
where U, and are displacement, velocity and acceleration
Modal method (cf. Chapter 11.4.3) has some limitations concerning damping forms and non-
linearity types which are taken into account.
The direct method is stepwise integration of the physical ordinary differential equations. This
method is used for all viscous damping forms and for all non-linearities (material, geometry,
contact..). The time-related discretization of general equations uses a time-related scheme.
Authorized values for parameter ν vary from 0 to 1. If ν=0, we obtain the Euler algorithm. If
ν=1/2, we obtain the trapezoid algorithm with a velocity linear approximation between t and
t+dt. For ν ≥ 1/2 the time related scheme is unconditionally stable.
Dynamic analysis of structures
The time-related discretization is obtained from Newmark time-related scheme:
DEFINITION. For linear structural behavior however, several load cases may be specified
which will then be combined during the calculation into a single load case (FUNCTION
command).
l The method of solution (by default) is suited to band equations.In this case, it is necessary to
ensure that nodes are numbered in such a way as to minimize the bandwidth. This may be
done by means of the routine RENUMBER.
The METHOD command (see Chapter 21) is used to choose others equation solvers.
Computed results are essentially displacement or temperature cards, velocities, accelerations and
possibly reactions and forces.
These cards are written to the file TRANSIENT if requested.
12.1.4 Commands
The routines in which the direct method of calculation of dynamic response is implemented are
called by the command
TRANSIENT (parameters)
The parameters indicate the type of analysis to be carried out. Distinction is made between:
l Analysis of linear thermal or elastic behavior,
l Analysis of non-linear dynamic behavior,
The command TRANSIENT is followed by commands and data relevant to the chosen type of
analysis.
All of this data is described in detail in the rest of this chapter.
The command
RETURN
concludes the set of commands and data read under TRANSIENT, and initiates the
corresponding calculation of response.
or by the command :
LINEARIZE this label is used to transform quadratic elements (type 1003, 2006,
2008, 3010, 3015 and 3020) into linear elements without creation of
nodes. This transformation is internally performed. The stiffness
(conductivity) matrices calculation is carried out by an internal
assembly of subdivision of quadratic elements into linear elements.
Forces and stresses are computed on the original elements in order to
perform standard post-processing.
The routine for transient analysis is called by the command
or
The label SCALING specifies the use of an alternative method of diagonalization of the thermal
capacity matrix C. If it is omitted, the leading diagonal term is taken to be the sum of all the
terms of a line ; this procedure sometimes results in negative leading diagonal terms. If the label
SCALING is present, the diagonal matrix is obtained by multiplying the initial leading diagonal
by a scale factor, chosen so that there is no change in global thermal capacity.
Note
l If the STIFFNESS and MASS files does not exist during the input in the TRANSIENT module, they
are automatically computed. In the presence of quadratic elements (type 1003, 2006, 2008, 3010,
3015, 3020), the LINEARIZE label is automatically taken into account.
ALGORITHME IMPLICIT ν
where :
j is the number of load cases specified in the input data,
ik is the number of a table, defined by the routine TABLE, in which is
given the function of time fik by which the load cases k is to be
multiplied.
If Qk is the load defined by load case k of the input data, the total load Q(t) at any time is given
by :
If in the command FUNCTION the value of ik is given as zero, then fik (t) = 1.
The absence of the command FUNCTION is equivalent to the unity function.
The functions fik(t) must be defined by the routine TABLE, called in the usual way by a
command followed by groups of data :
TABLE
i / type (parameters)
i / FORTRAN
i / PROG
i / FILE n CARD j
The command
METHOD ...
is used to chose the solver for the linear equations system. The chapter 21 gives a detailed
description of this command.
The recursive algorithm requires the definition of initial temperatures.
By default, these temperatures are taken to be zero at all nodes.
Otherwise, initial temperatures may be specified by means of a routine called by the command
RESTART Indicates that initial temperatures are to be taken from the file
TEMPERATURE (or DISPLACEMENT), which must previously
have been defined.
CARD It is also possible to take initial temperatures from a card of a previously
created TRANSIENT file. In this case, the label RESTART must be
followed by CARD n, n being the number of the card non equal to 0 to
be selected. The results are then stored on the same TRANSIENT file,
after the n- card.
LAST Indicates that the last stored card is taken without specifying its number.
In this case, the label LAST must follow the label CARD in place of
the number n.
This procedure enables an analysis to be subdivided into several stages, with each stage
continuing from the point reached by the last. It also allows the possibility of repeating a stage
of the calculation, with a shorter time step for example.
Finally, initial temperatures having been set as described above (and in particular having been
set to zero if the label RESTART is omitted), it is possible to selectively reset nodal
temperatures to any specified values by following the command INITIAL CONDITIONS by
groups of data of the form
(nodes) / TT=θ
where :
(nodes) is a list of nodes
θ is the initial temperature of those nodes.
Examples:
The initial temperatures are read from card 8 except for nodes 10 and 20 where prescribed
temperature is 150 degrees.
INITIAL CONDITIONS
TT = 50
This sets the initial temperatures of all the nodes of the structure to 50 degrees.
The recursive solution algorithm also requires a definition of the time interval. This is done by
means of the command
where
INITIAL to to is the initial value of time. By default, to = 0.
COMPACT This label activates the compact numbering of the transient file. This
label can be used when some of the cards are not stored.
RESTART The initial time is taken from the specified card of the transient file, in
the case of a restart calculation.
SKIP This label avoids re-defining time steps in the case of a restart
calculation.
Time steps are specified by groups of data of the following form :
where :
l ti is the upper limit of the time interval [t(i-1), ti]
TIME
4 STEP 2 / TEMPERATURE 1
10 STEP 3 / TEMPERATURE
15 /TEMPERATURE
RETURN
This initiates the transient analysis, resulting in the calculation, printing and storage on file of the
requested temperature cards.
DEFINITION
THERMAL SHOCK
OPTION THERMAL SPATIAL
PARALLELEPIPED
X Y Z 0.1 I J K 5
ADD 1
MATERIAL PROPERTIES
E / KX = KY = KZ = 0.75 RHO = 1602 C = 840
CONSTRAINTS
PLANE ZZ 0 / TT
PLANE ZZ 0.1 / TT
LOAD
1 TEMPERATURE 1 / NOTHING
PLANE ZZ 0.1 / TT 1
RETURN
TRANSIENT LINEAR
FONCTION 1
TABLE
1 / 1 0 0 0.1 100
TIME
1000 STEP 100 / STORE 1 ; Transient analysis
1000 / STORE 1 ; Steady-state calculation
RETURN
The results shown in figures 12.1 and 12.2 below represent the temperature distribution at time
step 1000 for respectively transient analysis and steady-state.
12.3.2 General
Equation (1) is solved by a time-stepping procedure.
Given the values of U, et at time t satisfying equation (1), values of U, and at
time t + ∆t, are sought such that the equation is still satisfied.
To this end a recursive algorithm is used (see Chapter 12.1.2), in which at each time step a
system of equations whose matrix depends upon ∆t is solved.
Provided the time step ∆t does not change, it is not necessary to reconstruct and triangularize the
matrix at every step ; given the high cost of this operation, considerable reductions in computing
time may be achieved by adopting a constant time step.
In addition, to optimize the performance of the computer routine, the damping and mass
matrices C and M are diagonalized. Particularly in the case of damping matrices with the form
C = αK, the user must verify that the physical properties of the model are not too much
modified by the diagonalization operation; The diagonalisation does not take into account the α
K damping introduced on the material properties of elements or on the DAMPING command of
TRANSIENT ELASTIC.
or by the command
or
The label SCALING specifies the use of an alternative method of diagonalisation of the mass
and damping matrices. If it is omitted, the leading diagonal term is taken to be the sum of all the
terms of a line separately for translations and rotations. Rotational diagonal dampings are set to
zero. If the label SCALING is present, the diagonal matrix is obtained by multiplying the initial
leading diagonal by a scale factor, chosen so that there is no change in global mass or damping.
The label SCALING is required with quadratic elements.
Other commands, possibly followed by data, provide details of the analysis.
In the absence of any indication to the contrary, the structure is taken to be undamped (C = 0).
The optional ALGORITHM command is used to specify the values of the two parameters (β,
γ) of NEWMARK's time-related scheme (cf. Chapter 12.1.2.). Giving α parameter is a way of
automatically defining the βand γ parameters. If no label is present, a low numerical damping is
set (α = 0,05).
If α is given (with α ≥0, then and γ = 1/2 + α. If α ≤ 0 : no numerical damping
or
If β is given (β≥ 1/4), then γ = 2β if γ is not given
where:
j is the number of load cases specified in the input data.
Ik is the number of a table, defined by the routine TABLE, in which is given
the function of time fik by which the load case k is to be multiplied.
If Fk is the load defined by load case k of the input data, the total load F (t) at any time is given
by:
TABLE
i / type (parameters)
i / FORTRAN
i / PROG
i / FILE n CARD j
METHOD ...
is used to choose the solver for linear equations solver. The Chapter 21 gives a detailed
description of this command.
The recursive algorithm requires the definition of initial displacements and velocities.
By default, these displacements and velocities are taken to be zero at all nodes.
Otherwise, initial displacements and velocities may be specified by means of a routine called by
the command
RESTART The label RESTART indicates that initial conditions are issued from a
previous analysis. If the CARD label is not present, only the
displacements are to be taken from the local DISPLACEMENTS file.
In this case, they are supposed coming from a static computation; initial
velocities and accelerations are equal to zero.
CARD Initial displacements and velocities are to be taken from a card number
n of a TRANSIENT file previously computed. In this case, the label
RESTART must be followed by CARD n, n being the normal number
of the card to be selected. The results are then stored on the same
TRANSIENT file, after the n-th card.
If one of the PLASTIC, VISCOELASTIC or MODIFY labels is
present after the BEHAVIOR command, FORCES corresponding to
card number n are also retrieved.
If the number of recomputed cards is inferior to the initial number of
cards of the file, the last cards of the initial file are saved. To erase
these last cards, the SELECT TRANSIENT command must be used.
LAST Indicates that the last stored card is taken without specifying its number. In
this case, the label LAST must follow the label CARD in place of the
number n.
This procedure enables an analysis to be subdivided into several stages, with each stage
continuing from the point reached by the last. It also allows the possibility of repeating a stage
of the calculation, with a shorter time step for example.
Finally, initial displacements and velocities having been set as described above (and in
particular having been set to zero if the label RESTART is omitted), it is possible selectively to
reset nodal displacements and velocities by following the command INITIAL CONDITIONS
by groups of data of the form:
(nodes) / (U = u), (V = v)
where :
(nodes) is a list of nodes.
(U = u) are displacement labels chosen according to the type of computation option,
each followed by an initial value.
(V = v) are velocity labels chosen according to the type of computation option, each
followed by an initial value.
Example:
The recursive solution algorithm also requires a definition of the time interval. This is done by
means of the command
where
INITIAL to to is the initial value of time. By default, to = 0.
COMPACT This label activates the compact numbering of the transient file. This label
can be used when some of the cards are not stored.
RESTART The initial time is taken from the specified card of the transient file, in the
case of a restart calculation.
SKIP This label avoids re-defining time steps in the case of a restart calculation.
All the macro time steps for which the final time (ti) is lower to the initial
restart time step (to) are automatically skipped.
Time steps are specified by groups of data of the following form:
Where:
ti is the upper limit of the time interval ((ti-1), ti)
∆ti is the time step to be taken in the interval ((ti-1), ti).If necessary the routine
adjusts ∆ti to be an exact sub multiple of (ti - t(i-1)).
TIME
2 STEP 0.1 / DISPLA 4 STORE 1 REACT 2
10 STEP 0.5 / DISPLA 2 VELO REAC FORCE 3
RETURN
This initiates the transient analysis, resulting in the calculation, printing and storage on file of the
requested results.
12.3.4 Relations
This option, available under the command TRANSIENT ELASTIC is used for modeling
contacts either internal or external to the structure.
It is not allowed for static analysis.
During the transient response:
l Initially unconnected nodes may make contact, or vice versa;
l Initially unconstrained nodes may make contact with a support, or vice versa.
The linkage between two nodes of the structure, or with a support is represented by a linear
spring of finite stiffness normal to the surfaces which make contact. This structural component
comes into play when a limit of travel is exceeded, and ceases to be active if subsequently the
relative displacement falls below that limit.
Interpolation within a time step is automatically carried out, in order to determine exactly the
instant at which contact is made or broken.
Consequently, each change of state of a contact introduces shorter time steps, chosen so that the
change coincides with the end of a step.
Each change produces new cards of results. As a safeguard, if the number of time steps
considered exceeds twice the number requested, the analysis is terminated.
Relations must previously be defined under the command DEFINITION, by means of the
command RELATION (cf. 4.10).
This command is followed by groups of data as shown below:
RELATION
ni nj / (U) MIN = α MAX = β ABS = γ
ni / (U) MIN = α MAX = β ABS = γ
where :
ni, nj are the numbers of nodes which may make contact ;
(U) is a displacement label, chosen according to the type of computation option.
The first group of data specifies a contact interval within the structure, for which
α< U (ni) - U (nj) < β
The second group of data specifies a contact between the structure and a support, for which
α < U (ni) < β
The value γ of the stiffness of the spring which comes into play when contact is made must
always be given.
This stiffness acts in a direction defined by the local axes of the nodes as specified under the
command CONSTRAINT.
The analysis remains numerically stable, provided the stiffnesses of springs are of the same
order of magnitude as structural stiffnesses associated with the degrees of freedom involved.
can be used. This command, starting from the TRANSIENT file created in TRANSIENT
analysis, computes FORCES associated with the selected cards.
Files
Input
DATA
TRANSIENT including DISPLACEMENTS, VELOCITIES, ACCELERATIONS and
REACTIONS
Output
DATA
TRANSIENT including DISPLACEMENTS, VELOCITIES, ACCELERATIONS,
REACTIONS and FORCES
Optional labels
SEARCH f By default the input TRANSIENT file is copied on unit 30. The output
TRANSIENT file is created on unit 19 as usual (copy on unit 30 may
consume CPU time for large files). If the file TRANSIENT has previously
been created on a permanent medium by a SAVE command, the copy
operation may be omitted using the SEARCH label associated with the
TRANSIENT file identification number used on SAVE command.
COMPLEX and REAL
These labels are used to associate the dynamic forces (real or complex)
with the harmonic response (real or complex).
CARDS i TO j STEP k
The selected cards are copied to the output TRANSIENT file.(the number
of cards that may be selected is limited to 65 536, they conserve their initial
number).By default all cards are selected.
EXTRACT Control printout of results for each card.
FILTER The accelerations are recomputed by a natural physical formula
(derivation) from velocities. The forces in beam elements are computed
from filtered accelerations. This lead to acceleration and force levels
more physical than those obtained directly with accelerations of
Newmark's method.
Note
l After a calculation with the TRANSIENT NON LINEAR command, the size of the FORCE file
being reduced by SELECT TRANSIENT, the command SOLVE FORCES TRANSIENT (RESIDUE)
can be used.
Key :
optional
mandatory
Table 12.2 - Input and output files
12.4.1 Introduction
For non-linear dynamic problems, finite element analysis involves repeated solution of the
following type of equation:
ψ = F - φ - fa - fi = 0 (1)
Where
ψ = residual nodal force vector
F = vector of external forces applied to the structure
φ= vector of internal forces due to displacements of the structure
fa = vector of damping forces due to velocities of the structure
fi = vector of inertial forces due to accelerations of the structure.
All these vectors are functions of the state of the structure at time t, the state of the structure
being defined by:
U : nodal displacement vector
fa = C.
fi = M.
where:
K = stiffness matrix
C = damping matrix
M = mass matrix.
For linear behavior, matrices K, C and M are constant.
For non-linear analysis, equation (1) must be solved repeatedly using an iterative method.
where:
and a new residue Ψ i+1 , corresponding to ui+1, and so on until the residue is less than the
precision factor required by the user (PRECISION command):
Matrix , representing the left-hand side of the equations to be solved, is only used
to supply a convergence operator. It is not always necessary to recompute and retriangulate on
each iteration (OPTIMISE label of the ALGORITHM command). Where the BFGS method
(BFGS label of ALGORITHM command) is selected, the left-hand side matrix corresponds to
the cutting matrix between two successive iterations. In certain cases where convergence is
particularly difficult to obtain, in particular in the presence of large strains, it can be necessary to
compute tangent matrix with precision. A perturbation method can be used for
this purpose (PERTURBATION label of ALGORITHM command):
with
where n and m refer respectively to the nth and mth components of vectors ψ and U.
In all cases, the essential requirement is that residual forces decrease on each iteration, and tend
to cancel each other out.
Furthermore, to facilitate convergence, a prediction of the solution based on known quantities at
previous instants is made systematically for dynamic problems, for the first iteration of the
iterative process for solution of equation (1). For static problems, it is also possible to activate a
solution prediction algorithm (PREDICTION label of ALGORITHM command). This
optional algorithm should only be used for quasi-proportional and monotonic load cases, and
involves estimating the solution at the current time step, on the basis of solutions for the two
preceding time steps:
In contrast to linear computation, mass matrix M and damping matrix C are processed integrally
in this case (consistent, non-diagonal matrices).
In theory only the initial displacement and initial velocity are necessary to define the
problem. With Newmark algorithms γ≠2β, initial acceleration is necessary. In all case, an initial
acceleration which respects the equilibrium equation at initial state, is computed:
Code specifying whether matrices and load vectors are to be printed for
EXTRACT i
checking purposes.
i = 0 : no printout
i = 1 : printout of the left-hand side matrix of each element KD (cf.
12.4.2.2)
i = 2 : printout of right-hand side matrix for each element
i = 3 : printout of both matrices for each element
i = 4 : printout of force data for each element, the material of which is non-
linear (plasticity, creep, etc.), before and after correction of stresses.
Level 1 command TRANSIENT NON-LINEAR is followed by a set of level 2 commands
defining the following:
l Solution algorithm,
l Initial conditions,
l Time step.
RETURN
which starts computation of transient conditions, and printout of the required results.
Solution Algorithm
The command:
Notes:
l If the BFGS or GRAD label is present, it is sufficient to construct and triangularize the right-hand
side matrix once only, by specifying OPTIMIZE - n (where n is very large).
l However in certain cases, it can be appropriate to recompute (at start of increment) the matrix
during the analysis, if convergence cannot be obtained.
l In certain cases, it may be necessary to reconstruct and triangularize the left-hand side matrix
during the iterative process (total Lagrange formulation in large displacements, for example).
However it may be sufficient to execute this operation for the firs t few iterations only.
Example: ALGORITHM OPTIMISE 1 0 2
Triangularization is done on iteration 0 (start of increment), and then on iterations 1 and 2.
l If the PERTURBATION label is present, certain behavior models and modified updated Lagrange
formulation ( EULER label associated with BEHAVIOR command: see below) produce non-
symmetrical tangent matrices. It is therefore advisable to associate a non-symmetrical solver
(NON-SYMMETRICAL label associated with METHOD command: see below) with the
computation of tangent matrices by the PERTURBATION method.
l The PERTURBATION method leads to cumbersome computation, and should consequently be
reserved for cases where convergence cannot be obtained with conventional methods.
l The solution prediction method (PREDICTION label) is particularly efficient in the case of quasi-
proportional loads, namely where the direction of strain at any point varies little in time, and is
monotonic. However this method should be avoided in other cases.
l The RIKS method can be combined with the BFGS or OPTIMISE label on the one hand, and the
BERGAN or MEITERATION label on the other.
l In the case of the BERGAN or CRISFIELD method, time step ∆t or arc increment ∆s (if one of
these methods is combined with the RIKS method) are updated from the second increment.
Examples:
1/ RIKS method
The second method of constant arc-length incrementation is used between times 1 and 10.
Computation stops if the time exceeds 10, or if the number of increments performed since time
1 exceeds 30.
Note that time step DT = 0.5 in the second group of time step data determines the value of the
reference load, noted P in Chapter 19 (initial amplitude, noted , consequently has a value of
1).
2/ BERGAN method
Between times 1 and 10, time step DT (which had an initial value of 0.5), is automatically
updated at each increment according to variation of the BERGAN parameter.
For the CRISFIELD method, interaction between the data of the ALGORITHM command and
the time step data follows the same principle.
The command:
By default:
e = 10-6 for absolute precision
e = 10-5 for relative precision.
TOLERANCE t Required precision for relational equations, which must be expressed in
absolute values (even if the ABSOLUTE label is absent).
By default, t = 10-3.
NORM nor nor = 0 (default value) indicates the maximum norm:
Notes
l The convergence is dependent on the quantities specified by the user. If no criterion is defined,
the convergence is obtained with respect to the default force, displacement and energy precision
values.
l Convergence is obtained when the first criterion is met.
l The internal energy criterion is only available for non-linear mechanical problems.
l If the convergence is not obtained after the n iterations specified in the ALGORITHM command,
computation is stopped and an appropriate message is printed. However the partial results
obtained are printed and stored, for possible restart of the computation.
l When a convergence error occurs, the last card where the convergence is obtained is
automatically saved.
The command
is used to select the type of solver for linear equation systems. By default, a direct method of
solving symmetrical systems is selected. In some cases (contact media with friction, for
example), a non-symmetrical solver must be activated.
Iterative solution methods are generally suitable for processing systems with a large number of
equations.
Refer to Chapter 21 for a detailed description of the METHOD command and the different
solvers.
The rotational degrees of freedom are not used for convergence criteria for the shell spatial and
shell spatial multi options.
Structural Behavior
Material behavior and geometrical non-linearities
The command
BEHAVIOR...
is used to specify the material and geometrical non-linearities taken into account in a mechanical
analysis. If this command is absent, the materials are assumed to behave in accordance with
Hooke's law (linear or non-linear) with small displacements and small strains.
Behavior models, different geometrical non-linearity formulations and the contents of the
TRANSIENT file, are described in detail in Chapter 17.
Material non-linearities
PLASTIC Used to calculate a structure in which certain materials have:
- elastoplastic behavior
- elasto-visco-plastic behavior.
The label is used to activate the selected PLASTIC and/or CREEP
model type for each element.
Geometrical non-linearities
EULER Specifies integration of large displacement and large rotation
phenomena using eulerian variables (Cauchy stress tensor). This
formulation is compatible with material nonlinearities.
TOTLAGRANGE Specifies integration of large displacement and large strain phenomena
using the TOTAL LAGRANGE formulation.
UPDLAGRANGE Specifies integration of large displacement and large rotation
phenomena using an updated lagrangian formulation.
c/ TWO-DIMENSIONAL and THREE-DIMENSIONAL options
Material non-linearities
PLASTIC Used to compute a structure in which certain materials have:
- elastoplastic behavior
- elasto-visco-plastic behavior
- or viscoplastic behavior.
The label is used to activate the selected PLASTIC and/or CREEP model
type for each element.
VISCOELASTIC Used to assign a linear viscoelastic behavior to the structure.
HYPERELASTIC Indicates that certain materials have hyperelastic behavior.
Geometrical non-linearities
TOTLAGRANGE Specifies integration of large displacement and large strain phenomena
using the TOTAL LAGRANGE formulation.
This technique consists in dividing the imposed load at Gauss points in the elements, in several
load stages in order to compute the stress correction. The number of sub- increments is
calculated from the ratio: threshold/yield stress. It is limited to the value n entered by the user
(by default, 50).
Generalized plane strain
The command:
This command is used to impose a global behavior for a generalized plane strain formulation.
This option can be used to analyze a structure with translation geometry, subjected for example
to a heterogeneous temperature field in the cross- section of the structure, using a two-
dimensional cross-section model. The structure is then assumed to be situated between two
infinitely rigid planes, with mutual displacement of the planes in the form of a translation and
two rotations. One, two or three degrees of freedom can thus be added to the problem:
l translation on the z axis (DISPL label)
If none of the labels, (DISPL, ROTX or ROTY) are specified, the DISPL label is taken by
default. Furthermore, the ROTX and ROTY labels are only taken into account if the
translational degree of freedom is activated (DISPL label).
The user can also introduce a "bottom effect" in direction z, specifying total force f after the
BOTTOM label. A space function (defined by table g) can be associated with the bottom effect.
FUNCTION command
In contrast to linear analysis, only one load case may be specified in the input data for the
problem. The FUNCTION command is therefore generally omitted.
This command is still required however if a type 9 element is used (see 8.8 SPECIAL
ELEMENTS), the stiffness matrix and load vector of which are the result of static compression
using the sub-structuring method (see 7.2.5 EXTRACT BINARY).
FUNCTION i
i is used to assign the time function described in table i to the load of this element.
The various corresponding time functions are assigned directly to the loads concerned by the
DEFINITION procedure.
These time functions, as also the tables required for integration of the various non-linearities,
must be defined by a TABLE routine, called by the usual command, followed by its groups of
data:
TABLE
i / type, (parameters)
i / FORTRAN
i / PROG
i / FILEn CARDj
Initial Conditions
Recursive solution algorithms require definition of initial displacements and velocities.
By default, these displacements and velocities are taken to be zero at all nodes.
Otherwise, it is possible to specify initial displacements and velocities, using a routine called by
the command:
RESTART The label RESTART indicates that initial conditions are issued from a
previous analysis. If the CARD label is not present, only the
displacements are to be taken from the local DISPLACEMENTS file.
In this case, they are supposed coming from a static computation; initial
velocities and accelerations are equal to zero.
CARD Initial displacements and velocities are to be taken from a card number
n of a TRANSIENT file previously computed. In this case, the label
RESTART must be followed by CARD n, n being the number of the
card non equal to 0 to be selected. The results are then stored on the
same TRANSIENT file, after the n-th card.
If one of the PLASTIC, VISCOELASTIC or MODIFY labels is
present after the BEHAVIOR command, FORCES corresponding to
card number n are also retrieved.
If the number of recomputed cards is inferior to the initial number of
cards of the file, the last cards of the initial file are saved. To erase
these last cards, the SELECT TRANSIENT command must be used
LAST Indicates that the last stored card is taken without specifying its number.
In this case, the label LAST must follow the label CARD and in place
of the number n.
RIKS During a preliminary calculation performed with the RIKS method (cf.
12.4.3.2), the reference load (cf. 19.1) has been computed and saved in
the RIKS.TIT file. It is advisable to restart the calculation with this file
by adding the label RIKS on the command INITIAL CONDITIONS.
This procedure enables an analysis to be subdivided into several stages, with each stage
continuing from the point reached by the last. It also allows the possibility of repeating a stage
of the calculation, with a shorter time step for example.
Finally, initial displacements and velocities having been set as described above (and in
particular having been set to zero if the label RESTART is omitted), it is possible selectively to
reset nodal displacements and velocities by following the command INITIAL CONDITIONS
by groups of data of the form
where:
(nodes): list of nodes
(U = u) displacement labels selected according to the computation option, followed
by initial value at the nodes.
(V = v) velocity labels selected according to the computation option (elasticity
formulation only), followed by initial value at the nodes.
Example:
CONSTRAINT Activates the incremental release of rigid constraints. A blank character ‘ ‘
REMOVE must be preceded by the word CONSTRAINT.
Up to 10 lines can be defined.
GROUP list Group(s) of nodes to release.
The list can be a list of group names (between $ characters) and/or group
numbers.
Example: GROUP $grp1$ num2 num3 $grp4$ …
DURATION dt DURATION sets the duration of the time-dependent function to which is
multiplied the concentrated loads.
The linear table starts from value 1 at time tnc (time of the restart card) and
reaches the value 0 at time tnc + dt.
By default dt is equal to the difference between ending time of the transient
computation and restart time.
EXTRACT Lists information on the released nodes.
Note:
l The keyword CONSTRAINT REMOVE is available only for mechanical computations.
l Not all rigid constraints should be released to avoid rigid body motion.
l A node listed several times in the lists of nodes to release will be treated (constraint removal and
load assignment) only the first time encountered by the process.
Time Step
General
The recursive solution algorithm also requires definition of computation times. This definition
is made using a procedure called by the command:
where
INITIAL to is the initial value of time. By default, to = 0.
COMPACT This label activates the compact numbering of the transient file. This
label can be used when some of the cards are not stored.
RESTART The initial time is taken from the specified card of the transient file, in
the case of a restart calculation.
SKIP This label avoids re-defining time steps in the case of a restart calculation.
All the macro time steps for which the final time (ti) is lower to the initial
restart time step (to) are automatically skipped.
DTLAST In case of restart of a computation, this label overwrites the initial time
step with the previous one to reach first macro time.
The previous time step is computed from the difference between the
time of the restart card (defined in INITIAL CONDITION section) and
the one of the card that precedes it.
The overwriting of the time step is available only with automatic time
step, provided that the computed step is greater than the initial one.
Time steps are defined by the following groups of data:
where
ti = the upper limit of the time interval [t(i-1), ti]
∆ti = time step taken for computation in the interval [t(i-1), ti]. The value of
ti is adjusted by the routine, to a sub-multiple of [ti - t(i-1)].
By default, ∆ti = ti - t(i-1)
By convention, a time step equal to zero (ti = t(i-1)) leads to the computation of a static steady
state.
Printout and storage of results are performed according to information given in the input data
field.
TIME
2 STEP 0.1 / DISPLACEMENT 4 STORE 1
10 STEP 0.5 / DISPL 2 VELO 2 STORE 2 DISPLA 2
Starting from a time step ∆t which is already known (time step previously predicted or test time
step suggested by the user in the TIME data of the procedure TRANSIENT NON LINEAR),
the program tries to compute a solution at t+ ∆ t. If this solution cannot be obtained (non
convergence of computation), then a new time step is suggested:
∆tp = α∆t
where α is a user-defined parameter (by default, α=0.5). This operation is repeated until the
computation converges. ∆ t min being the minimum time step authorized by the user,
computation stops as soon as ∆tp < ∆tmin.
Having a solution at ( being the time step where convergence has been reached), the
general rules apply to the automatic time step adjustment are the following:
Checking the solution in time interval [t, t+∆t].
The solution obtained at t+∆t assumes predefined evolution of the displacement field between t
and t+∆t :
l linear variation for a static mechanical analysis,
and there is nothing to indicate that this evolution is compatible with load evolution.
When the residue at is acceptable, the program keeps the solution obtained at
and thus tries to take the next time step.
The prediction of the new time step is based on two estimations. The minimum value of these
two estimations is kept. The first estimation is based on the previous time step (between t and
) and is used to estimate the time step which could have been chosen (in place of )
while keeping an acceptable solution at the end of the time step. The second estimation is based
given for thermal, thermo-metallurgical and electromagnetism options to define the precision in
term of "displacement tolerable gap".
l By adding the label FORCE in the input data, this algorithm can be made predictive only. It
suppresses the control of the solution in the interval .
Input data
This algorithm is incompatible with the RIKS, BERGAN and CRISFIELD automatic
incrementation method.
with:
dtmin minimum time step authorized by the user (∆t < dtmin --> computation
stopped) - The DTMINI label must be given.
dtmax maximum time step authorized by the user (optional).
LABEL ε
acceptable precision for the control of the solution at and for the
prediction of the solution.
LABEL:
EPSD, if the acceptable precision is expressed in "DISPLACEMENT"
for the thermal and thermo-metallurgical options , ε = 5 degrees by
default.
EPSF, if the acceptable is expressed in "FORCE". By default,
where f is the precision in force, indicated on the command
PRECISION (ABSOLUTE or RELATIVE precision).
fc time step multiplication factor if convergence not obtained.
By default, fc=0.5.
fp maximum authorized time step multiplication factor. By default, fp=2.
nres number of successive time that the residue can increase during a card.
resm limiting value of the maximal residue.
The FORCE label avoids to control the solution at , thus the model becomes predictive
only.
The label WATCH allows to deactivate the automatic time stepping algorithm until a
convergence issue occurs. The time step is given by the user.
Note that the label FORCE is automatically activated in that case because it switches off the
control at .
DAMPING Command
DAMPING
This command is required to specify concentrated damping at the nodes (cf. 4.12), and damping
for each element (labels ALPHA and BETA for most elements - cf. 4.5) described in the
physical properties.
SIL Command
SIL processing_sil P1 P2 P3
This command executes the « sil_processing.cmd » process before and after each card with
three parameters:
P1: 0: call before the card, 1: call after the card.
//
// Example
//
ibeg = p1
time = p2
ncard = p3
if (ibeg = = 0) {
printf(« Execution before the card %d\n », ncard) ;
else {
printf(« Execution after the card %d\n », ncard) ;
}
RUNEND Command
The end of run can be controlled either by a final time (Cf. Time Step Section), or by stop
criteria. These stop criteria are managed by the command RUNEND.
"Card number" stopping criterion (thermal computation)
The "card number" criterion lets the user stop a transient thermal computation after a given
number of cards. The functionality is activated with keywords RUNEND CARD:
TRANSIENT NON-LINEAR
...
RUNEND CARD (nc) (INITALIZE) (EXTRACT)
...
RETURN
TRANSIENT NON-LINEAR
...
RUNEND THERMAL (GOTO n) (EXTRACT)
...
RETURN
l Material non-linearities
Other non-linearities can be combined, such as those associated with the different MEDIA
(Chapter 8.11), and contact media in particular, and those resulting from the introduction of
special RELATIONS (Chapter 12.4.5).
The ALGORITHM command defines:
l BFGS,
l NEWTON (OPTIMISE 1 label),
l Modified NEWTON method (no BFGS or OPTIMISE labels, with the tangent matrix
l BERGAN
l CRISFIELD
Table 12.3 indicates elements for selection of the most efficient algorithm according to the non-
linearity (or non-linearities) concerned.
The algorithms are assigned a priority index for each non-linearity.
Automatic incrementation methods can be associated with all algorithms defined above (BFGS,
Newton, etc).
Table 12.4 describes possible combinations of non-linearities and incrementation methods.
12.4.5 Relations
It is possible to integrate additional relations or inequalities between degrees of freedom
previously defined with RELATION commands (see 4.10) under the DEFINITION procedure.
Equations
Degrees of freedom at a number of nodes can be related by an equation with the form:
bij Uj = gi
An unknown Lagrange multiplier, interpreted as a force, is associated with this equation.
This equation must be described using the RELATION command under DEFINITION.
Two types of relation can be defined:
a) General relation
RELATION
(nodes) / TABLE (t)
where n is the node number, and j the number of the degree of freedom. The table comprises
pairs of values:
where k is an integer associated with node n and degree of freedom j by the relation k =
10.no+j, no being the position in which node n appears in the list of nodes in the RELATION
command.
Where appropriate, value g gives the right-hand side of the equation.
Example
RELATION
1 8 5 7 / TABLE 5
TABLE
5 / (12,1) (22,.5) (31, -7) (43,2.3) 10.
Note
l If all the nodes indicated are assigned the same system of axes by a CONSTRAINT command (see
4.6), the relation is expressed in the corresponding axes.
b) Displacement equalities
RELATION
(nodes) / (U)
Inequalities
a) Contact type
The TRANSIENT NON- LINEAR procedure can integrate contact phenomena, whether
internal or external to the structure.
During the dynamic response:
l Initially separate nodes can come into contact with each other, or vice versa,
l Initially connected nodes can come into contact with a support, or vice versa.
Nevertheless the iterative procedure for a given time value ensures contact status stability.
l No particular user action is required.
RELATIONS
ni nj / (U) MIN = α MAX = β (1)
ni / (U) MIN = α MAX = β (2)
where:
ni, nj numbers of nodes which can come into contact.
(U) displacement label according to the computation option.
The first group of data defines a contact internal to the structure:
α <U(ni) - U(nj) < β
The second group of data describes a contact between the structure and the support:
α <U(ni)<β
The inequality is expressed in the local system of axes of the nodes (defined by the
CONSTRAINT procedure).
b) Friction
Decoupled friction
There is no gap to be closed, and contact is assumed to be already established.
A conventional relation of equality between one or two degrees of freedom is introduced.
However the associated force exceeds a specified value, a quality is not longer maintained, the
force remains constant, and sliding occurs.
RELATION
ni n j / (U) ABS = - val
ni / (U) ABS = - val
In this case, the relation is only active if a contact relation acting on the same node or nodes and
another degree of freedom, is itself active.
Friction without sliding occurs as long as the associated force does not exceed a limit value,
dependent on contact force Fn.
RELATION
ni nj / (U) ABS = α
ni / (U) ABS = α
where :
ni nj numbers of two nodes between which friction can occur
(U) displacement label
α coefficient of friction (α > 0)
Sliding occurs when Ft ≥ α . Fn (Fn = normal contact force)
Example
RELATION
45 26 / UY MIN -.3 MAX 3.
45 26 / UX ABS = .5
RELATION
10 15 / UX MAX 4.
10 15 / UY ABS = .3
10 15 / UZ ABS = .3
c) Axes
Contact and friction relations operate in axes associated with the nodes by the CONSTRAINT
command in the DEFINITION procedure (see 4.7).
where
represents the number of iterations required to determine normal contacts (Lagrange
(i)
multipliers).
represents the number of iterations required to determine the Lagrange multipliers associated with
(j)
tangential displacement during sliding.
(N) number of relations.
(k) equilibrium iteration number for a given time card.
represents sliding displacement of a node between two equilibrium iterations i and i-1. The
∆bi
value (residue) must tend towards zero to obtain equilibrium for an iteration.
Relation status
Relation Stress Limit Code Force Node-component
1 0.000E+0 0.000+E0 1 0.321 52 2 63 2
2 -100. 0.000E+0 0 0.00E+0 52 2 63 2
3 -0.279E-01 0.600 3 0.192 52 1 63 1
Internal node
number where
Relation Imposed gap Relation Force at
ABS value component to
number or computed status nodes
which relation
applied
Interpretation of relation status codes:
0 Inactive
1 Active for coupled contact-friction relation
2 Intermittent friction stop without sliding (FT < Fn)
3 Intermittent sliding contact ( FT = Fn)
-1 Active for decoupled contact-friction or displacement equality relation
-2 Permanent friction stop without sliding
-3 Permanent sliding contact
Rather than maintaining a TRANSIENT file, SYSTUS enables the user to store minimum
results for subsequent reconstruction of any card of the TRANSIENT file, in a HISTORY file
(mass storage file).
The HISTORY file contains all displacements, in double precision, for all computation steps.
Furthermore, in the presence of "contact" type macro elements (CONTACT, GAP,
INTERFACE and PLANE), the HISTORY file contains the PRESSURE file, for all computed
time steps, in addition to displacements.
To store results in the HISTORY file, it is merely necessary to add the ARCHIVATE label to
the TRANSIENT NON-LINEAR command.
The HISTORY file can be manipulated with the usual SAVE, DELETE, SEARCH and
ASSIGN commands (unit 25).
SAVE HISTORY i
DELETE HISTORY i
SEARCH HISTORY i
ASSIGN 25 HISTi.TIT BINARY
The computation cards stored in the HISTORY file can be listed by the following command
sequence:
The user can then request storage of any precomputed card by searching the HISTORY file,
and specifying the REBUILD label on the TRANSIENT NON-LINEAR command. The
routine then reintegrates the material behavior law from the displacement history. This
recreation phase delivers a conventional SYSTUS TRANSIENT file.
The BEHAVIOUR and FUNCTION commands must contain the same information as for the
initial analysis.
The INITIAL CONDITIONS command defines either initial conditions (defined in initial
analysis) for integration of the behavior laws, or the card from which integration of the behavior
laws is made (recreation procedure restart). In this case, the TRANSIENT file containing the
card mentioned in the command is required at input. By default, recreation starts at card 0.
The CARD label defines the list of computation cards for which results must be stored.
This function is available for all computation options recognized by the TRANSIENT NON-
LINEAR procedure (Except for generalized plane strain option).
l Iteration for each time value, in which residual forces are minimized.
Time steps are used only to define the quasi-static application of loads to the structure. Velocity
and acceleration are not involved.
The procedure for minimization of residual forces is the same as described in Chapter 12.4.2.
Where the quasi-Newton is selected, convergence is substantially accelerated by modification of
the iterative procedure.
The essential factor is that residual forces decrease on each iteration, tending to zero.
12.5.3 Utilization
As mentioned above (Chapter 12.4.3), the STATIC label must appear in the TRANSIENT
command.
All other parameters and commands have the same meaning as before.
The recursive solution algorithm only requires definition of initial displacements in this case.
For mechanical problems, the RIKS, Bergan or Crisfield automatic incrementation methods can
be used. Furthermore, the RIKS method can also be used for structural instability problems.
Only the labels (U = u) appear in the groups of data following the INITIAL CONDITIONS
command.
In the same way, the only results printed and stored are DISPLACEMENTS, REACTIONS
and FORCES (no velocities or accelerations).
12.5.4 Substructuring
Introduction
The sub- structuring method has numerous advantages for dealing with large- scale linear
problems, including the sub-division of computation, possibility of elementary modifications,
optimum integration of structure symmetries, etc.
In a large number of non-linear mechanical problems, the structure has non-linear behavior
(plasticity, large strains, etc.) in a restricted area, which is generally clearly demarcated. It is
then of interest to isolate this area, and compute, once and for all, the elastic matrices and right-
hand side of equations for the other areas, which have a linear behavior (elasticity, small strains
and small displacements).
To take another example in the field of electromagnetism, there is a special medium with linear
behavior. It can be relevant to consider air as a special sub-structure, in order to concentrate
non-linear analysis on the fellow-magnetic areas.
Other examples can also be given, but in general terms, the use of sub-structuring in a non-
linear analysis leads both to a reduction of computation time, and an increase in the size of the
problems which can be handled.
It is obvious that the advantage which can be obtained from this method, depends directly on
the type of problem and the quality of sub-division into sub-structures (super elements).
On the other hand, all areas showing non-linear behavior must be grouped in a single, level 0
sub-structure.
DEFINITION
---
---
RETURN
SUPER ELEMENT nl
Computation Procedure
The conventional TRANSIENT NON-LINEAR procedure is used. Only the SUPER
ELEMENT nl label must be added.
No DATA file is needed at the start of the procedure, as the latter works directly on the data of
sub-structure n1.
The DATA file delivered at output is compatible with the TRANSIENT file, and contains:
l all data for sub-structure n1
l connecting nodes between the standard sub-structures which do not belong to sub-structure nl.
Example:
X nodes added
Furthermore, codes are modified according to the nature of the problem.
At the present time, only results concerning the non-linear geometry are accessible, and the
DATA and TRANSIENT files produced by the procedure must be used together for post-
processing.
For a restart, only the TRANSIENT file is required as input data for the procedure, as the latter
works directly on the permanent data file for SUPER ELEMENT n1. On the other hand, the
RENUMBER and FUNCTION commands
Example:
SEARCH TRANSIENT p
or
ASSIGN 19 TRANp.TIT BINARY
TRANSIENT NON-LINEAR SUPER nl
----
FUNCTION i1, i2 ..., ij
----
RENUMBER
(renumbering front )
----
INITIAL CONDITION RESTART CARD c
----
RETURN
SAVE DATA TRANSIENT p
Limitations
These limitations are particular to the method. All areas considered as "standard" sub-structures
are assumed to have linear static behavior. In particular, this means that mass or inertial effects
are ignored in standard super structures. Furthermore for mechanical problems, standard super
structures are created on the assumption that the parts of the structures to which they
correspond, have a linear elastic behavior (small stresses and small displacements).
and includes nodal temperatures and gradients from a TRANSIENT SHELL THERMAL
calculation.
The command:
defines the type of non-linearity to be used. By default, the computation is linear elastic.
l In this case, the above values of ityp correspond respectively to values 3, 4, 5 or 6 of ityp defined
in the table 12.6.1.
Example:
TIME
1.5 STEP 0.5
4. STEP 0.5
ALGORITHM RIKS 2 AMPL 20 ITAN 1 IT1 1 IT2 2 IT3 1
BUCK TYPE 15 INCR 8 EIGE 2
BEHAVIOUR LAGRANGE
COND REST CARD 3
The increment number to be specified necessarily corresponds to the last increment of the Riks
method.
General Notes on the algorithm.
If all the eigenvalues obtained are negative, the algorithm employs the constant arc-length
incrementation.
The type of incrementation can be modified during a restart procedure.
For elasto-plastic materials and when using the present method, it is advisable to define ITAN
1IT1 1IT2 2.
A check parameter S (Bergan's parameter) is provided, which characterizes the global stiffness
behavior under increasing loads.
If S > 1, the structure becomes stiffer, and if S < 1, the structure becomes softer.
At a limit point type of buckling, S = 0.
The parameter S to be considered is the one printed at the end of the first iteration of each
increment.
Definition of the reference load with
TIME
t1 STEP dt1
t2 STEP dt2
....
tn STEP dtn
ALGORITHM RIKS igroup AMPL iampl
STATIC
or
TIME
ti step dti / (S=s) STORE n CREEP dtm RESTART nr
dtm If negative, the present data indicates the table (-dtm) of the creep
function defined under the command DEFINITION.
If positive, (dtm) defines the dwell time during each increment of the
group number (i). If (tflu) is the dwell time during increment j, then the
dwell time during increment (j+1) is (tflu+dtm). Reminder: Creep
analysis is activated, only if keyword CREEP is declared under the
command BEHAVIOR.
RESTART nr Storage of data useful for restart on file RESTART.DAT, every nr
times steps.
By default nr is the minimum of s values.
It is not possible to chain a first run with a restart run in the same
session.
It is necessary to perform the two computations separately.
The following command governs the print-out of the results:
PRINT COMPLETE
Example
12.6.2 Thermal
This particular transient feature allows thermal analysis (linear and nonlinear) of structures
modeled with the thermal shell elements (SHET).
The thermal transient analysis for SHET elements is activated by the command:
TIME
ti STEP dti / (S=s) STORE n RESTART m
Notes
l In stationary analysis: use ti = dti = 0.
l A transient analysis must be initiated by a stationary time step (see also 12.6.1).
The following command governs the printout of results:
PRINT (COMPLETE)
o THERMAL
o THERMAL METALLURGY
l thermal shell elements (code 240x) which are associated or not with "continuous media"
l thermal shell elements (code 240x) followed by a computation using the option spatial
SHELL SHE8
l thermal shell elements from SHET thermal shell option.
This routine creates a file including temperatures at nodes. Provided there is no abnormal
behavior, the time value for each card selected is printed.
It is possible to specify, either a list of cards, or a range of cards including the two instants tinit
and tfinal using label TIME or the last card only with the label LAST.
It is also possible to shift the time range computed during the thermal analysis by using the label
DT. For example, if the thermal results have been computed from 0 to 100 seconds, by
specifying DT 50, the instants sent to the thermomechanical analysis will be from 50 to 150
seconds. This can be useful for mechanical structures submitted to a cyclic thermal loading
saving the useless duplication of thermal computation.
Since the temporary file is lost at the end of a run, the TEMPERATURE command must be
executed before any thermo-mechanical analysis by the TRANSIENT NON-LINEAR routine.
In the mechanical analysis, the times specified in groups of data following the TIME command
need not correspond to those of the thermal analysis. If necessary, temperature fields are
interpolated linearly over a time step of the thermal analysis to obtain values required by the
mechanical analysis.
Rules used for the interpolation are:
l If temperatures are required for times less than the lowest value of time for which a card has
been calculated, or greater than the highest value of time for which a card has been
calculated, they are taken to be the temperatures at those lowest or highest values of time,
respectively.
l otherwise, the interpolation is made between values of time of the thermal analysis which
The family of thermal shell elements (code 24XX) of THERMAL ENTHALPY computation
option (Heat Transfer Reference manual - Chapter 6) is compatible with mechanical elements
(code 24XX) of SHELL spatial option.
This routine creates 2 files:
l a file including temperatures at nodes,
l a file including thermal information at shell elements.
d) Thermo-mechanical analysis for a structure partially modeled with enthalpy thermal
shell elements (type 240X) and with "continuous media with metallurgy" elements from
THERMAL METALLURGY MULTI option
The GRADIENT label is used to transfer thermal gradients in addition of mean temperatures at
model nodes.
This routine creates a file including thermal information at nodes of the model.
g) Thermo-mechanical analysis for a structure modeled with shell thermal elements
SHET followed by a spatial SHELL SHE8 computation option
To use results from a transient thermal analysis using SHET shell elements (TRANSIENT
THERMAL SHET procedure) for a thermo-mechanical analysis with SHE8 sell elements
(TRANSIENT SHELL procedure), the input data are as follows:
GRADIENT label (optional) is used if the TRANSIENT file created in thermal analysis
contains temperature gradients through the thickness of the shell besides temperatures.
This routine creates a file including thermal information at nodes of the model.
Example:
number of degrees of freedom, but with a very large number of time steps, for these had to be
very short. As a result a very large number of cards are written to the TRANSIENT file.
l 2) A static non-linear analysis with the labels PLASTIC or EULER LAGRANGIAN has
been carried out for a structure modeled by a mesh with a large number of degrees of
freedom. The amount of data in each card is greatly augmented, in particular by force or
stress history data.
Furthermore, in the second of these situations the volume of force or stress data can render
impossible the use of certain modules (POSTPROCESSING for example).
The SELECT command may be used to get around certain of these difficulties.
The OUTPUT command described in Section 3.3 may also be indirectly used to achieve
processing selectivity.
Input files:
DATA (use SEARCH to attach this file: see also notes)
TRANSIENT (automatically attached by the program: SEARCH is irrelevant).
Output files:
DATA (modified)
TRANSIENT (if the label STORE is present)
Syntax
SELECT TRANSIENT n
CARDS <list> or LAST or LAS2 (labels) COMPONENTS i1 i2 --
(NOSEQUENCE) (STORE) (PRINT) (ENVELOPE)
CARDS (numbers) (labels) COMPONENTS j1 j2 --
(STORE) (PRINT) (ENVELOPE)(TIME (t0))
RETURN
SELECT Command
TRANSIENT n Number of the TRANSIENT permanent file in which the cards will be
selected.
CARD Command
(Numbers) or Specifies which cards are to be selected: given by their numbers or by
LAST or LAS2 keywords LAST for the last card or LAS2 for the last two cards.
Several CARD commands can coexist in the same SELECT
command.
(Labels) DISPLACEMENT
VELOCITY
ACCELERATION
REACTION
FORCE
PRESSURE
VALUE
These labels specify which results are to be processed. Several of these labels can coexist in the
same CARD command.
COMPONENT i1 Specifies which components of force (or stress) are to be selected. If i2 neq
i2 0, components i1 to i2. If i2 = 0, components 1 to i1. By default (or if
i1=i2=0), the first KPC components will be selected. This selection is
useful in restoring the force data to their standard size (see further the
relation between KPC and the option).
STORE If present, creates a new TRANSIENT file which can be made permanent
by means of the SAVE command.
NOSEQUENCE If present, the original card numbers are retained on the new
TRANSIENT file. This label is only effective on the first CARD
command, in order to prevent erasing records of the new TRANSIENT
file. By default, if this label is not used, the cards written on the new
TRANSIENT file are renumbered.
PRINT If present, the selected results are printed out.
ENVELOPE If present, the extreme values of the selected results and the time at
which they occur are printed out.
TIME (t0) If present, the time of the selected cards are translated: the first card takes
for time value t0 and the next cards, the values t0+Δt with Δt equal to the
initially existing difference of time between first selected card and the next
ones (default value: t0 = 0).
Notes
l Employ SEARCH DATA before SELECT, if you previously used another SELECT command
containing a STORE command.
l If in the same SELECT you use several CARD commands in which force components are stored.
The number of such components must remain constant.
l If not, the program uses by default the number of components in the first such command.
option 1 2 3 4 5 6 7
KPC 4 3 6 6 3 6 4
Default values for number of stresses in linear elasticity analysis.
where :
Option 1 = all types of structure with d.o.f. UX, UY, RZ
Option 2 = all types of structure with d.o.f. UZ, RX, RY
Option 3 = all types of structure with d.o.f. UX, UY, UZ, RX, RY, RZ
Option 4 = all types of structure with d.o.f. UX, UY, UZ
Option 5 = thermal and electromagnetic options
Option 6 = all types of structure with d.o.f. UX, UY, UZ, RZ
Option 7 = all types of structure with d.o.f. UX, UY.
Example:
When de-refining, if met = 1 (by default), then the value assigned to the
target Gauss point is equal to the value of the Gauss point of the nearest
source mesh.
If met = 2, then the value assigned to the target Gauss point is equal to
the mean value of the Gauss points of the nearest source mesh.
PARAMETER If the point to project fall exterior to the element but within a distance
EXTERIOR EPSILON e (relative to the element size) it is consider as a surface
point.
PARAMETER If the point to project fall interior to the element but within a distance
INTERIOR EPSILON e (relative to the element size) it is consider as a surface
point.
PARAMETER Relative distance of the boundaries (default e=0)
EPSILON e
CARDS <list> List of cards to be projected and stored in the target transient file.
If no card is specified, then all cards are selected.
LAST Indicates that the last stored card is taken without specifying its number.
In this case, the label LAST must follow the label CARD in place of
the number n.
This procedure delivers a transient output file which can be used for a new analysis. However
it is impossible to recompute mean values of the FORCES file.
With the second procedure, the complete known result history for the initial mesh, contained in
a HISTORY, file is projected onto the new mesh, and a coherent TRANSIENT file for the new
mesh is then reconstituted (REBUILD label with TRANSIENT non-linear command). This
method is more cumbersome, but does allow full compliance with the material behavior law
(total coherence between nodal displacements and results for element integration points).
Interpretation of the labels is the same as for the first procedure, apart from label CARDS,
which must appear in the form indicated above (0 TO n). The projection procedure delivers a
HISTORY output file concerning the target data. A TRANSIENT file can then be
reconstituted using the TRANSIENT NON-LINEAR REBUILD procedure (Chapter 12.4).
Limitations
The procedure is limited to the following computation option elements:
l OPTION TWO-DIMENSIONAL (all sub-options except generalized plane strain option)
l OPTION THREE-DIMENSIONAL
l OPTION ELECTROMAGNETIC AV
Input data for using the "Mapping Engine Component" tool is as follows:
(NODE n1 n2 n3 In case of source and target models are not at the same position and/or
(INTERNAL) not in the same orientation, the labels NODE following both SOURCE
(layer)) and TARGET are necessary to compute the translation and/or rotation
operations of the source model to fit geometrical position of the target
one. Nodes m1, m2 and m3 of the target are then the result of the
transformation leading nodes n1, n2 and n3 from source geometry to the
position of the target one.
Remarks:
l Nodes Ids can be set by their user numbers (default) or their internal one
(INTERNAL must be specified then).
l Layer stands for LOWER, MIDDLE (default) or UPPER:
In case of mapping of a shell model onto a solid one or vice versa, the
layer the nodes of the shell model belong to must be set with respect to
their position in the solid model.
Example for a local frame defined at the lower surface:
Limitations
The procedure is limited to the following computation option:
l TWO-DIMENSIONAL (all sub-options)
l OPTION THREE-DIMENSIONAL
l OPTION SHELL SPATIAL
l OPTION SHELL SPATIAL MULTI
l OPTION THERMAL (all sub-options)
l OPTION THERMAL MULTI (all sub-options)
l TWO-DIMENSIONAL
l THREE-DIMENSIONAL
Four labels are associated with the activation/deactivation of elements. These labels complete
the normal properties describing the behavior of the material.
STATE -f specifies the number (f) of the table giving the activation/deactivation
function.
GAUSS val if val = 1, changes of state are assessed and applied at element
integration points.
if val = 0, changes of state are assessed at the centre of gravity of the
element, and applied globally to the complete element.
DEFORMED d if d = 1, changes of state are taken on the deformed geometry.
REDUCTION factor used to multiply the properties of an element when the latter is
xk deactivated.
When this value is negative (for the mechanical options only), the
Young modulus of the deactivated elements is equal to –xk, note that
the Poisson ratio is multiplied by 10-6.
Deactivated Elements
Mechanical problems
The behavior of a deactivated element is assumed to be elastic, and the Young’s modulus (E)
and Poisson’s ratio values are multiplied by reduction factor xk.
Forces are cancelled and the element is assumed to have zero mass (RHO = 0).
Thermal enthalpy
The thermal conductivity (KX,KY,KZ) of a deactivated element is multiplied by factor xk.
Heat flux is cancelled, and the element is assumed to have zero mass (ρC = 0).
Activated Elements
Mechanical problems
Thermal strain is cancelled on activation.
Stresses and strains are computed from the moment when the element is activated. The initial
strain of the element at the moment of activation is ignored.
12.9.7 Examples
Thermal Problem
Elements are activated in this example. A band with a width of two elements is initially in the
inactive state (initial conditions). A heat source is moved across the band, activating the
individual elements in the process. The heat source and the function activating the elements are
superimposed in this example. The elements are activated permanently (F = 0 outside the
activated zone). The part is allowed to cool after 2.5 seconds.
DEFINITION
2D WELDING
OPTION THERMAL ENTHALPY PLANE
TWO-DIMENSIONAL
NODES
1/
2 / 10
3 / 10 10
4 / 0 10
ELEMENTS
1001 / 1 2 3 4
2001 / 2 3
EDGES
100 / 1 2 20
200 / 2 3 20
300 / 3 4 20
400 / 4 1 20
MATERIAL PROPERTIES
E / RHO 1 KX=KY=KZ 0.05 C 0.006 STATE -10
CONSTRAINT
ELEM 2001 TO 2020 / KT 0.01
LOAD
1 / NOTHING
ELEM / QR 35 VARI -20
ELEM 2001 TO 2020 / TT 100.
…
TABLE
1 / 1 0 0 5 10
10 / FORTRAN
FUNCTION F(X)
DIMENSION X(6)
XX = X(1)
YY = X(2)
ZZ = X(3)
TT = X(4)
VT = 2*TT
AA = XX
BB = YY - VT
AA = AA * AA
BB = BB * BB
CC = AA + BB
CC = SQRT(CC)
F = 0.
IF (CC.GT.2. ) RETURN
IF (TT.GT.4.) RETURN
F = 1.
RETURN
END
20 / FORTRAN
FUNCTION F(X)
DIMENSION X(6)
XX = X(1)
YY = X(2)
ZZ = X(3)
TT = X(4)
VT = 2*TT
AA = XX
BB = YY - VT
AA = AA * AA
BB = BB * BB
CC = AA + BB
CC = SQRT(CC)
F = 0.
IF (CC.GT.3. ) RETURN
IF (TT.GT.4.) RETURN
CC = CC/2.
CC = - CC
F = EXP (CC)
RETURN
END
RETURN
RENUMBER ITERATION 30
1 INTERN
RETURN
SEARCH DATA 1
ASSI 19 TRAN1.TIT
TEMPERATURE TRANSIENT CARD
DEASSIGN 19
Mechanical Problem
The thermal analysis for the previous example is used again in this case. The elements are
consequently activated in the same way. Note that the structure should be anchored to a node
not belonging to a deactivated element, where this is possible. On the other hand, symmetry is
naturally applied to ad hoc elements, including deactivated elements.
SEARCH DATA 1
DEFINITION
THERMO-MECHANICAL COMPUTATION
OPTION TWO-DIMENSIONAL PLANE THERMO
RESTART GEOMETRY
CONSTRAINT
NODE LINE 1 4 / UX
NODE 2 / UY
MATERIAL PROPERTIES
ELEMENTS / E=-11 H = 0.2 LX=LY=LZ=15*-6 --
LIMIT=-15 MODEL=1 NU=0.3 STATE -10
LOAD
1 / NOTHING
TABLE
1 / 1 0 0 5 10
10 / FORTRAN
FUNCTION F(X)
DIMENSION X(6)
XX = X(1)
YY = X(2)
ZZ = X(3)
TT = X(4)
VT = 2*TT
AA = XX
BB = YY - VT
AA = AA * AA
BB = BB * BB
CC = AA + BB
CC = SQRT(CC)
F = 0.
IF (CC.GT.2. ) RETURN
IF (TT.GT.4.) RETURN
F = 1.
RETURN
END
11 / 1 20 200000 300 190000 500 165000 750 70000 1000 20000 1200 50
15 / 1 20 700 200 625 400 550 550 400 750 100 850 50 1000 25 1250 1
RETURN
…
CHAPTER 1 INTRODUCTION 1
1.1 Overview 1
1.1.1 Program Execution 1
1.1.2 SYSTUS 3
1.1.3 Data Structure 3
1.1.4 Preprocessing 4
1.1.5 File Management 4
1.1.6 Postprocessing 4
1.1.7 Field of Application 5
1.1.8 Substructuring 6
1.1.9 Command Compatibility with New Data 7
1.1.10 Environment Variables 10
1.2 Basic Concepts 18
1.2.1 Theoretical Basis 18
1.2.2 Discretization 18
1.2.3 Degrees of Freedom 18
1.2.4 Coordinate Systems 19
1.2.5 Boundary Conditions 19
1.2.6 Loads 20
1.2.7 Node and Element Numbering Systems 20
1.2.8 Units 20
1.3 Command Language 21
1.3.1 Overall Operation of SYSTUS 21
1.3.2 Language Characteristics 22
1.3.3 Language Structure 26
1.3.4 Automatic Number Generation 29
1.3.5 Node and Element Groups 30
1.3.6 Comments 31
1.3.7 Batch and Interactive Processing 32
1.4 Notation 35
1.4.1 Boxes 35
1.4.2 Labels and Numbers 36
1.4.3 Prefixes 36
1.4.4 Choice of Labels 37
1.4.5 Coordinate System 38
1.4.6 Node and Element Groups 40
1.5 Groups of Nodes or Elements 41
1.5.1 Creation of a Group 42
CONTENTS
INDEX i
INDEX i
INDEX i
INDEX i
EXPORT 97
INDEX EXTRACT
FILE
322, 333
300
FRACTURE 300
FUNCTION 299
(VOL.1) GEOMETRY 118
GOTO 54
GROUP 42
A HISTORY 64
HYDROSTATIC 299
ASCII 72 IMPORT 97
Automatic Number Generation 29 INCREMENTATION 301
INPUT 83
INPUT/OUTPUT 80
B LEVEL 299
LEVELSETS 195
Batch Mode 32 LINEARISE 266
Boundary Conditions 19 LIST 331
LOADS 146
MASS 169
C MASTER 293
MATERIAL PROPERTIES 132
Command MEDIA 136, 319
$FORMAT 123 MISES 299
ADD 241 MODE 63
ASEMBLE 264 MODIFY 251
ASSIGN 56 MULTI 310
ASSOCIATE 273, 280, 290-292, 295 NODES 214
AXES 122, 238 OFFSET 259
CHECK 324 ONE-DIMENSIONAL 204
CONSTRAINTS 138 OPTION 114, 301
CONTACTOR 283, 290 OUTPUT 80
COORDINATES 120, 238 PAGING 63
COORDINATES CARTESIAN 120 PARALLELEPIPED 202
COORDINATES CYLINDRICAL 120 POINTS 54
COORDINATES POLAR 120 RBE 293
COORDINATES SPHERICAL 121 READ 57
CORRECT 239, 262 RECTANGLE 200
Crack 300 RELATION 161
CRACK 260 RELATION ELIMINATION 164
CRACKS 197 RELEASE 158
CRITERION 274, 281, 290-291, 295 REMESH 298
DAMPING 172 RENUMBER 243, 261
DEASSIGN 57 RESTART 244, 250
DEFAULT 134 RETRO 309
DEFINITION 108, 238 RETURN 108
DEFINITION RESTART 110 REWIND 58
DELETE 77 SAVE 72
EDGE 215 SEARCH 75
ELEMENTS 125, 215 SECTION 270
END 62 SECTOR 207
EQUILIBRIUM ENERGY 337 SINGULARITY 215
EQUILIBRIUM FORCE 334 SLAVE 293
EQUILIBRIUM MASS 336 SLIDING 281, 290
EQUILIBRIUM STIFFNESS 335 SOURCE 274, 290-291
ERROR 54 STOP 58
N V
Node and Element Numbering Systems 20 Vectorial Functions of Several Variables 189
O
Operating System Commands 64
P
PAMSYS 101
Periodic Structures Subjected to Fourier
Decomposed Loads 117
Postprocessing 4
Prefixes 36
Preprocessing 4
Programmed Meshes
Specification of Trajectories 314
Types of Element Generated 312
R
RBE Elements 292
REFERENTIEL
DEGRES 38
GRADES 38
RADIANS 38
TANGENTE 38
S
Simplified Fortran 185
Singularities 220
Slide Plane 277
Substructuring 6
T
Thermal Analysis in Linear Transient 171
Thermal Analysis in Non-Linear Transient 174
Thermoelastic Analysis 116
SHELL 565
INDEX SOLID
SOLVE HARMONIC
565
464
Computed Results
(VOL.2) Axisymmetric or Translation Shell Option
Beams Option
401
377
Grids Option 361
Plane Truss Option 382
A Plate Option 394
SHED axisymmetric harmonic shell option 553
Applied Loads Shell Axisymmetric Fluid Harmonic Option 493
Axisymmetric or Translation Shell Option 399 Shell Fluid Option 498
Beams Option 373 Shell Harmonic Option 475
Grids Option 359 Shell Option 431
Plane Beam Option 351 Shell SHE8 Option 545
Plane Truss Option 380 SHET Thermal Shell Option 558
Plate Option 393 Three Dimensional Harmonic Option 485
SHED Axisymmetric Harmonic Shell Option 551 Three Dimensional Option 462
Shell Axisymmetric Fluid Harmonic Option 492 Truss Option 388
Shell Fluid Option 498 Two Dimensional Option 450
Shell Harmonic Option 473 Constraints
Shell Option 422 Axisymmetric or Translation Shell Option 398
Shell SHE8 Option 541 Beams Option 371
SHET Thermal Shell Option 557 Grids Option 358
Three Dimensional Harmonic Option 483 Plane Beam Option 350
Three Dimensional Option 459 Plane Truss Option 379
Truss Option 386 Plate Option 392
Two Dimensional Option 447 SHED Axisymmetric Harmonic Shell Option 550
Axisymmetric of Translation Shell Option Shell Axisymmetric Fluid Harmonic Option 492
Applied Loads 399 Shell Fluid Option 497
Computed Results 401 Shell Harmonic Option 472
Constraints 398 Shell Option 421
Material Properties 396 Shell SHE8 Option 538
SHET Thermal Shell Option 557
B Three Dimensional Harmonic Option 482
Three Dimensional Option 458
Truss Option 384
Beam/Solid Transition Element 510 Two Dimensional Option 446
Beams Option 363 Contacts 582
Applied Loads 373
Computed Results 377 Continuous Three Dimensional Frames 363
Constraints 350, 371 Coupled Analyses 488
Geometry of Section 366 Criterion
List of Elements 364 Hill 417
Material Properties 365 Tsai 418
C E
Coaxial Cylinders 500 Eccentricity
Command Shell option 411
BEAM 565 Electromagnetic Analyses 487
COMPATIBILITY 565
DEFINITION MEDIUM 609
MATERIAL PROPERTIES 515
MIXING 565
F SC3D 601
SL2D 592, 601
SPHERE 597
Family TORUS 597
2003-2004 405
2006-2008 407
Material Properties
2203-2204 405 Axisymmetric or Translation Shell Option 396
2403-2404 405, 407 Beams Option 365
Grids Option 356
Fluid-Structure Interaction 489 Plane Beam Option 348
Plane Truss Option 379
G Plate Option 391
SHED Axisymmetric Harmonic Shell Option 549
Generalized Contact 597 Shell Axisymmetric Fluid Harmonic Option 491
Grids Option 355 Shell Fluid Option 497
Applied Loads 359 Shell Harmonic Option 467
Computed Results 361 Shell Option 410
Constraints 358 Shell SHE8 Option 535
Material Properties 356 SHET Thermal Shell Option 556
Three Dimensional Harmonic Option 479
Three Dimensional Option 454
I Truss Option 384
Two Dimensional Option 440
Interpretation of Computed Results 341 Mixing Option 559
L O
Label Option
FAST 584 Plane Beam 347
MULTI 560 Plane Truss 378
PENALTY 584 Plate 389
List of Elements SHE8 Spatial Shell 533
Beams Option 364 SHED Axisymmetric Harmonic Shell 548
Plane Truss Option 379 Shell 402
SHED Axisymmetric Harmonic Shell Option 549 Shell Axisymmetric Fluid Harmonic 490
Shell Axisymmetric Fluid Harmonic Option 491 Shell Fluid 495
Shell Fluid Option 496 Shell Harmonic 465
Shell Option Option 403 SHET Thermal Shell Option 555
Shell SHE8 Option 535 Three Dimensional 451
SHET Thermal Shell Option 556 Three Dimensional Harmonic 477
Three Dimensional Harmonic Option 479 Two Dimensional 435
Three Dimensional Option 454
Truss Option 384
Two Dimensional Option 438
P
Plane Beam Option
M Applied Loads 351
Constraints 350
Macro Elements Material Properties 348
CONE 597 Plane Truss Option
CYLINDER 597 Applied Loads 380
GAP 583 Computed Results 382
INTERFACE 587 Constraints 379
NIT3D 605 List of Elements 379
PLANE 589, 597 Material Properties 379
RECTANGLE 597
INITIAL
INDEX CONDITION
S 776, 869, 879, 903, 937, 941
INITIAL CONDITIONS RESTART 905
(VOL.3) INTEGRATION 793
LEVEL CROSSING 795
LIMITS 677
A LINK 659, 749
LOAD 632, 735, 759
MAXIMA 795
Automatic Link Generator 658
METHOD 845, 900
MODAL ANALYSIS 770
C MODE 674, 677, 680, 682-683, 687,
698, 705, 786, 805, 826
Combination between Non-Linearities and MODIFICATION 830
Algorithms 913 MODIFY 727
Command MOMENT 793
ACCELEROGRAM 845 PARTICIPATION 758, 761
ACCELEROGRAM TRANSIENT 846 PARTICIPATION MULTI 762
ALGORITHM 775, 894, 930 PRECISION 898
ALGORITME IMPLICIT 868 PRINT 779, 942
ASSEMBLE DYNAMIC 731, 734 PRINT COMPLETE 938
AXES 831 PRINT VALUES 676
BANDWITH 795 PROJECT 953
BEHAVIOUR 900, 929, 941 PSD 791
CALCULATE 632, 656, 792-794 PSEUDO MODAL 833
CALCULATE DYNAMIC 728 REGIME 941
CARD 780, 792-793, 950 RELATION 883
CHECK 659 RELEASE 659
CIRCULAR FREQUENCY 786, 826 RESIDUAL TERMS 762
COMBINATION 778 RESPONSE 825, 831
COMBINE 805 RESULTS 779
COMPLEX MODES 710 RUNEND 911
COMPONENTS 726 SELECT 676, 771, 779, 785, 834, 842, 949
CONNECTION 658, 727, 749 SENSITIVITY 830
CONTACT 632 SEQUENCE 691
DAMPING754, 772, 785, 806, 825, 842, 845, 910 SIL 910
DEFECT TRANSIENT 937 SOLVE 613, 616, 668, 877
DIMENSION 941 SOLVE BUCKLING 695
DIMENSION MEMORY 929 SOLVE FORCE TRANSIENT 885
DIRECTION 806, 814 SOLVE FORCES MODAL 689, 707
DISPLACEMENT 813 SOLVE HARMONIC 622
DURATION 806 SPECTRAL ANALYSIS 805
EFFECTIVE MASS 852 SPECTRAL ANALYSIS MULTIPLE 812
ELEMENT 710, 834 SPECTRUM 806, 813
FATIGUE 795 STATIC 812, 933
FREQUENCY 786, 826 STIFFNESS 727
FREQUENCY INITIAL 845 STOCHASTIC 794
FRONT 699 STOP 772
FUNCTION 771, 779, 785, 791, 878, 903 STRAIN GENERALISED 902
GAP 632 SUBSTRUCTURE 653
GENERATION 845 SUBSTRUCTURE DYNAMIC 725-726
HARMONIC RESPONSE MODAL 785 SUMMATION 695
SUPER ELEMENT 650
SUPERELEMENT DYNAMIC 725
SUPPRESS 676
TABLE 711, 772, 779, 786, 792, 807, 834, 878, G
903
TEMPERATURE 624, 944 Generation of Spectral Response 844
TEMPERATURE TRANSIENT 944
TEMPERATURE TRANSIENT SHELL 929, 947 H
TIME 777, 941
TITLE 771, 779 Harmonic Loads - Direct Method 823
TRANSFER 786
TRANSIENT 864, 877-878
Harmonic Response 783
TRANSIENT NON LINEAR 892, 921
UNIDIRECTIONAL 773 L
VECTOR 680, 687, 699, 706, 711, 835
Complex Eigenfrequencies and Eigenmodes Label
(Non-Symmetrical Matrices 709 ABSOLUTE 814, 898
ACCELERATION 692
Complex Eigenfrequencies and Eigenmodes ALL 710
(Symmetrical Matrices) 704 AMPLITUDE 855
Computation of Effective Masses and ARCHIVATE 892
Effective Inertia 852 BEAMS 692
Coupled Friction 917 BERGAN 895
BFGS 894
Cylindrical Stop 773 BINARY 921
BLOCK 687, 730
D BUCKLING 674
CARD 624, 855, 869, 879, 904, 956
Decoupled Friction 917 CARDS 885
Detail of Time-Related Scheme 891 CIRCULAR FREQUENCY 677, 680, 687
COMPACT 870, 880, 906
DYNAMIC 687, 698, 793-794, 812, 825
COMPLEX 885
Dynamic Responses 665 CONDENSATION 698
CONSTRAINT REMOVE 905
E CONTACT 896
DATA 954
Eigenfrequencies and Eigenmodes 665 DEFORMED CARD 953
Enumeration of Eigenvalues (Sturm DESACTIVATE 910
DIAGONAL 675, 682, 687
Sequences) 694 DIRECTION 761
Equivalent Modal Damping 841 DISPLA 770
DISPLACEMENT 899
F DT tshift 944
DTLAST 906
File DURATION dt 905
FMDAT.DAT 895 EDIT 791
HISTORY 920 ENERGY 899
MASD 723 EQUIVALENT STRESSES 692
RESU 722 EULER 901
RIGD 723 EXTRACT 729, 885, 893, 905
SPECnf.DAT 848 FACTOR 680
TRANSIENT 920 FILE 732
VECD 723 FILTER 885
FLUID 771, 779
FLUID STRUCTURE 698
FORCE 706, 759, 770, 898
FORCES 675
FREQUENCY 677, 680, 684, 687
N Z
Non-linear Dynamic Analysis 888 Zero Level Substructure 650
Non-Linear Static Analysis 923
O
Option
Axisymmetrical Harmonic 802
Fluid-Structure 802
P
Participation Factors 756
Plane Frictional Stop 774
R
Real Eigenfrequencies and Eigenmodes with
Fluid-Structure Interaction 697
Response Power Spectral Density 790
S
Seismic Response - Spectral Method 800
Sensitivity Analysis 666, 829
Simultaneous Modification of Stiffness and
Mass Matrices 671
Stochastic Dynamics 790
Substructure Modes 714
Substructuring 925
T
Thermoelastic Analysis 623
Transfer Functions 783
Transient Dynamic 854
Transient Response (Modal Analysis) 767
Transient Shell 929
U
Unidirectional Frictional Stop 773
Unidirectional Stop 772
Utility Modules 666
MEDIUM 1117
INDEX MERGE
METHOD
1004
1324
MINMAX 1056
(VOL.4) MODE
NODES
1035
1013, 1025, 1117
NUMBER 1084
A OPEN 1081
POINTS 1013
POSTPROCESSING 1021
Anisotropic Linear Elasticity 1151 PRINT 1033, 1035
Average for Thick Shells 988 PRINT CONSTRAINTS 1037
PRINT POINT 1034
C REARRANGE 1039
RETURN 992, 1028, 1066
Code SECTION 1048
COMMON PAGE 1289 SELECT 1066
KDEBT 1289 SELECT ELEMENTS 1048
LCODE 1289 SET 1113
NCODE 1289 SHELL 1057
SOLVE 978
Command
SPLIT 1054
AVERAGE 981, 988
STRESS 1042
AVERAGE TRANSIENT 991
STRESS SECTION 1048
AVERAGING POST-PROCESSING 1099
THREE-DIMENSIONAL 1062-1063
AXES 1026
TRANSF 1027
CARD 991, 1116
TRANSIENT 1035
CLOSE 1081
UNITS 1048
COMBINE 994, 1028, 1063
USER 1048
COMBINE FOURIER 999
VISUALIZATION 1101
COMPONENTS 992
CORRECTION 981 Criterion
CROSS SECTION 1011 Drucker-Prager 1162, 1173
CURVE AVERAGE 1124 Stassi 1163, 1174
CURVE FUNCTION 1121 Tsai-Hill 1163, 1174
CURVE POST PROCESSING 1112 User 1164
CURVE SATURATION 1122 User-programmed 1174
CURVE TABULATION 1123 Von Mises 1161, 1172
CURVE TRANSIENT 1116 Cross Sections 1009
DAMPING GYROSCOPIC 1340
DEASSIGN33 1338 D
DISTRIBUTION 1062, 1064
EDGE 1014 Dataset 1085
ELEMENT 981, 1117
ELEMENTS 989, 991, 1014, 1024
EXTRACT CODES 1289 E
FORCE 1040
FORMAT 1113, 1117 Elasticity 1147
GENERATE 1062, 1064 Elastoplasticity 1159
GRID 1014
HEADER 1014, 1026
LINE 1014
I
LOAD 1048, 1063, 1066
LOCAL 1066 Isotropic Linear Elasticity 1149
MATERIALS 1025
MATRIX 1026
TUBE 1022
L UNBALANCE 1339
UPPER 989
Label VARIABLE 989, 1002
ALL 1040 VMISE 1023
AMPLITUDE 1001 Large Displacements and Large Strains 1247
ASYNCHRONOUS 1339
Local Axes
AXES 989
One-Dimensional Elements 1134
AXISYMMETRIC 1023
Three-Dimensional Elements 1143
BEAM 1024
Two-Dimensional Elements 1136
CARD 1035
COARSE 1066 Location of Gaussian Points for Calculation
COMPONENTS 1034 of Stress 983
COMPOSITE 1059
CRITERION HILL 1059 M
DATASET 1084, 1100
DELTA 991 Merging of Cards 1007
DISPLACEMENT 991
ENERGY 1040 Merging of Posts 1006
EXTRACT 989, 1046 Merging of Results 1003
EXTRAPOLATION 1099 Merging of Transient Files 1008
FINE 1066 Model
FORCE COMPOSITE 1024 Chaboche 1216
GYROSCOPIC 1338-1339 Generalized Maxwell Hardening (Linear
HILL 1060 Viscoelasticity 1227
INPLA 1023 Model=1, Perfect Plasticity 1172
INPUT 1100 Model=10, Ductile Damage Model with
INREV 1023 Isotropic Strain Hardening - Leblond
INTEGRAL 991 – Perrin Model 1195
INTRI 1023 Model=11, Plasticity with combined Strain-
LAST 991, 1035
Hardening - isotropic strain
LAYER 989, 1046, 1059
hardening and PRAGER kinematic
LOWER 989
strain hardening 1184
MASS 1040
Model=12, Plasticity with combined strain
MISE 1022
MODE 1035 hardening and RAMBERG-OSGOOD
MULTIPLY 1022 strain hardening law 1184
NORM 1119 Model=13, Ductile Damage Model with
ONLY 989 Combined Isotropic and Kinematic
OUTPUT 1100 Strain Hardening - Leblond – Perrin
PART 1100 Model 1196
PHASE 1001 Model=14 Gurson - Tvergaard Model 1197
PRESSURE 1099 Model=2, Plasticity with Isotropic Strain-
PRINCIPAL 1022 Hardening 1175
PRINT 1084 Model=2, Plasticity with Prager Kinematic
SHELL 989 Strain-Hardening 1179
SHELL COMPOSITE 989 Model=21 Global Perfect Elastoplasticity
SPATIAL 1012 Model Applicable to Thin Shell
STRAIN 1046 Elements - Simplified Model 1197
STRESS 991 Model=22 Global Perfect Elastoplasticity
SUM 1118 Model Applicable to Thin Shell
THIN SHELL 1023 Elements 1198
THREE DIMENSIONAL 1022
THRESHOLD 991, 1099
TRANSF 1024
TSAI 1059
N
Non-Linear Elasticity 1154
Non-Unified Elasto-Viscoplastic Models 1210
P
Postprocessing 1018
enables manufacturers to deliver smarter and connected products, to predict product performance and to
anticipate maintenance needs.
ESI is a French company listed in compartment B of NYSE Euronext Paris. Present in more than 40 countries,
and addressing every major industrial sector, ESI Group employs about 1200 high-level specialists around the
world and reported annual sales of €141 million in 2016. For more information, please visit www.esi-group.com.
All PAM- and SYS- product names as well as other products belonging to ESI's portfolio are tradenames or trademarks of ESI Group, unless specifically mentioned.
All other trademarks are the property of their respective owners - Specifications are subject to change without prior notice.
www.esi-group.com