Transient Results and Response Spectra: Sofistik 2020
Transient Results and Response Spectra: Sofistik 2020
Transient Results and Response Spectra: Sofistik 2020
SOFiSTiK | 2020
DYNR
Transient Results and Response Spectra
SOFiSTiK AG
This manual is protected by copyright laws. No part of it may be translated, copied or reproduced, in any form or by
any means, without written permission from SOFiSTiK AG. SOFiSTiK reserves the right to modify or to release
new editions of this manual.
The manual and the program have been thoroughly checked for errors. However, SOFiSTiK does not claim that
either one is completely error free. Errors and omissions are corrected as soon as they are detected.
The user of the program is solely responsible for the applications. We strongly encourage the user to test the
correctness of all calculations at least by random sampling.
Front Cover
Project: Queensferry Crossing | Photo: Bastian Kratzke
Contents | DYNR
Contents
Contents i
3 Examples 3-1
3.1 Dynamic Load Factors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-1
3.2 Story Response Spectrum . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-3
SOFiSTiK 2020 i
DYNR | Contents
ii SOFiSTiK 2020
Task Description | DYNR
1 Task Description
The program DYNR enables the graphical presentation of transient results from the programs
DYNA and HYDRA. It plots time-controlled or frequency-controlled results and, as an additional
option, calculates response spectra for a given loading or acceleration history.
Using the input record HIST during a dynamic analysis with the program DYNA one can store
in the database the time variations of selected degrees of freedom or internal forces and mo-
ments. By the same record HIST, the program HYDRA stores in the database the unsteady
temperature or potential variation. DYNR can present these curves graphically in any combi-
nation.
Furthermore, it is capable of using either the already stored acceleration variations of individual
points or any newly defined function in order to determine the response spectra.
Each loading (e.g. earthquake, impact at corner A, impact at corner B, etc.) is stored in DYNA
by a certain number.
DYNR can now plot either the given HIST curves or the response spectra in the various coor-
dinate directions under consideration of the various damping values.
2 Input Description
Record Items
SIZE DINA S W H MARG FORM TMAX
SCHH H1 H2 H3 H4 H5
RESP TYPE TMIN TMAX TDIF
D1 D2 D3 D4 D5
C1 C2 C3 C4 C5
REF1 RC1 REF2 RC2 REF3 RC3
FUNC T F DT NO TITL
EXCI LCS NO COMP F TITL
SPEC TYPE
HIST LCS TYPE
NO0 C0 NO1 C1 NO2 C2
NO3 C3 NO4 C4 NO5 C5
NO6 C6 NO7 C7 NO8 C8
NO9 C9 YMIN YMAX
TMIN TMAX TDIM YTOL LIM1 LIM2
TEFF CLIM
XTYP NX TITL
ECHO OPT VAL
The records HEAD, END and PAGE are described in the general manual SOFiSTiK: "‘FEA /
STRUCTURAL Installation and Basics"’.
The vertical scale is selected with S. In addition, the value range can be arbitrarily changed by
HIST.
More explanations can be found in the general manual SOFiSTiK: ’FEA / STRUCTURAL In-
stallation and Basics’.
A value of 0. suppresses the corresponding designation. The values remain valid until a new
definition of SIZE or SCHH is given.
The single mass oscillator, which is subjected to the various excitations, must be described via
a frequency range with fixed subdivisions.
In the case of excitation through an acceleration, TYPE AA can be used to plot the total accel-
eration (ground+oscillator). This is usually the desired value. For TYPE one can also enter P+,
S+, V+, A+, AA+ as well as P-, S-, V-, A-, AA-. Then only the positive or negative components
are evaluated.
The evaluation of the oscillation equations is performed only for the times defined by the exci-
tation function. The duration of the evaluation is also determined by the highest time value of
the excitation value. A somewhat longer time range with exciting function 0.0 should thus be
defined in the case of high frequency excitation.
Up to five different damping values can be examined simultaneously. The damping values D1
through D5 can be used arbitrarily according to the desired marking. Different line types and
colour values can be defined for the curves. The colour values and the line types are added
here to the values defined inthe record EXCI minus 1001 (e.g. 2001 ”+” 3002 - 1001 = 4002).
Therefore it is reasonable to input one record with line types (1001 through 1010) and the
others with colours (e.g. 1001 through 9001).
Up to three external reference curves may be shown. These functions have to be defined via
FUNC and must be defined after an external load function for the response analysis.
Any functions that can be selected either as acceleration or loading of the single mass oscillator
or as comparison curves to HIST can be defined in DYNR. Since the mass of the single mass
oscillator is 1.0, no need of the differentiation exists there between acceleration and force.
The values of all excitation functions get stored in the database. They are maintained until
a new definition takes place. If FUNC values are input after the records HIST or SPEC, all
functions are automatically deleted.
NO and TITL are input by the first function value only. If no time value is entered, the last time
increment is added to the last time value. In this case for regular time intervals the input of time
values can be omitted for all higher values.
EXCI activates the variation curves that are stored in the database. All curves selected with
EXCI are plotted simultaneously.
If the excitation variation were defined in DYNR by the record FUNC, such a function must be
selected with LCS 0. The first defined function FUNC is used, if EXCI is not input.
When TYPE FREQ is input, the values TMIN, TMAX and TDIF of the RESP input are inter-
preted as frequencies too.
If LOGT or LOGF are defined, a logarithmic scale is selected. TDIF defines then a subdivision
factor relative to 1.0 for the number of intermediate values (e.g. TDIF 0.05 = 20 intermediate
values).
Each record SPEC or HIST produces a drawing with the current parameters.
With the record HIST one can request up to ten transient curves per picture.
The displayed range of values is determined automatically based on the found values and
the scale data given in the record SIZE. YMIN and YMAX or TMIN and TMAX can be used
to modify this range. If you specify the values for LIM1 and/or LIM2, these limit values will
be included in the display, if the current maximum values pass YTOL times the specified limit
values.
P PT Spring forces
SPX SPY SPZ Sum of spring forces or total drag force or coefficient
SPXX SPYY SPZZ Sum of spring moments or total drag moment or coeffi-
cient
TRUS CABL Truss bar force / cable force
N VY VZ Beam forces axial and shear
MT MY MZ Beam moments torsion and bending
SIG TAU SIGV Stresses in beams
MXX MYY MXY Shell moments
VXX VYY Shell shear forces
NXX NYY NXY Shell axial forces
H PN Potential, Pressure
QH QBH QHPI Flow quantities nodes, regions, pipes
VQUA VBRI Flow velocities QUAD, BRIC-Elements
U V W Flow velocities in monitor points
TKP TDN Turbulence parameter k and eps in monitor points
T Temperature
QT QBT QTPI Heat flow nodes, regions, links
QQUA QBRI Heat flow densities QUAD, BRIC-Elements
It is also possible to show similar values in the same picture together. This will be controlled
with an additional decimal at the number used as a shift. Thus for a specified TYPE U-X a
selection of 1001.0 is the displacement of node 1001 in X, while 1001.2 is that in Z direction
etc.
Beside the variations that are stored in the database, one can also define and plot comparison
curves with the record FUNC. The negative value of the function number must be given to one
of the NOi-entries for this purpose. In order to plot the functions only, one must enter the value
of 0 for LCS.
With LIM1 and LIM2 it is possible to show a constant limit or reference value with a horizontal
line and a printed value. The definition of literals MEAN or SDEV will show the mean values or
values shifted by the standard deviation. The definition of literals MAX, MIN or MAMI will show
the maximum and/or minimum values.
Each record SPEC or HIST produces a drawing with the current parameters.
ECHO
3 Examples
The input files which are explained here are to be found in the installation directory SOFiSTiK
in the subdirectory dyna.dat\english. Alternatively you will find these examples via the TEDDY
menu HELP > EXAMPLES sorted by program name and language.
For additional help like tutorials, tutorial movies and practical examples please refer to the
SOFiSTiK Infoportal (www.sofistik.com/Infoportal).
PROG DYNR
HEAD RESPONSE SPECTRA OF TRAPEZOIDAL FUNCTIONS
HEAD MANUAL DYNR
RESP P 0 6 0.10 0.0 0.20 1.00
FUNC 0.00 0. TITL 'TRAPEZOIDAL LOAD 1'
0.25 1. 0.025
0.75 1. 0.025
1.00 0. 0.025
SPEC
FUNC 0.00 0. TITL 'TRAPEZOIDAL LOAD 2'
0.33 1. 0.025
0.67 1. 0.025
1.00 0. 0.025
SPEC
FUNC 0.00 0. TITL 'TRIANGULAR LOAD'
0.50 1. 0.025
1.00 0. 0.025
SPEC
END
TRAPEZOIDAL LOAD 1
Damping 0.000
Damping 0.200
Damping 1.000
Force [kN]
1.8
1.6
1.4
1.2
1.0
0.8
0.6
0.4
0.2
0.0 Period
1.0 2.0 3.0 4.0 5.0 6.0 [sec]
TRAPEZOIDAL LOAD 2
Damping 0.000
Damping 0.200
Damping 1.000
Force [kN]
1.8
1.6
1.4
1.2
1.0
0.8
0.6
0.4
0.2
0.0 Period
1.0 2.0 3.0 4.0 5.0 6.0 [sec]
Damping 0.000
Damping 0.200
Damping 1.000
1.4
1.2
1.0
0.8
0.6
0.4
0.2
0.0 Period
1.0 2.0 3.0 4.0 5.0 6.0 [sec]
First an oscillation analysis of the entire structure with the program DYNA is necessary. This is
done e.g. by the following input:
PROG DYNA
HEAD THREE-STORY FRAME BASE POINT EXCITATION
! MAXIMUM VALUES U3 = 10.2 CM ; V = 1360 KN ; M = 2040 KNM
! 0.5 CM ; V = 92 KN ; M = 138 KNM
! 0.2 CM ; V = 56 KN ; M = 84 KNM
! NATURAL FREQUENCIES 1.501 / 3.185 / 4.862 (HERTZ)
EIGE 3 REST
STEP 200 0.06
LC 51
FUNC 0.0000 0.0000; 0.1000 -0.5765; 0.2000 0.1357
FUNC 0.3000 -0.5040; 0.4000 0.5614; 0.5000 0.6055
. . .
ACCE 0 0.0 3.0
HIST U-Y 4 LCST 51
HIST A-Y 3 4 1 LCST 51
HIST MYA 1 3 1 LCST 51
END
The time variation of the displacements as well as the accelerations of the upper node (No. 4)
is now available in the database. Then DYNR plots both these variations and calculates the
story response spectrum with the following input:
PROG DYNR
HEAD RESPONSE SPECTRUM OF ROOF STOREY
HIST 51 U-Y 4 1001
HIST 51 A-Y 4 2001
HIST 51 MYA 1 1001 2 2001 3 3001
RESP A 0 1.5 0.02 0.001 0.02
EXCI 51 4 Y 1001
EXCI 51 3 Y 2001
SPEC
END
Displacement uY 4
Displacement uY [mm]
100.000
50.000
0.000 Time
2.0 4.0 6.0 8.0 10.0 12.0 [sec]
−50.000
−100.000
Acceleration AY 4
Acceleration AY [m/sec2]
10.00
5.00
0.00 Time
2.0 4.0 6.0 8.0 10.0 12.0 [sec]
−5.00
−10.00
4000.00
3000.00
2000.00
1000.00
0.00 Time
2.0 4.0 6.0 8.0 10.0 12.0 [sec]
−1000.00
−2000.00
−3000.00
−4000.00
LOADCASE 51 NODE 4 Y
LOADCASE 51 NODE 3 Y
Damping 0.001
Damping 0.020
Acceleration [m/sec2]
200.00
150.00
100.00
50.00
0.00 Period
0.2 0.4 0.6 0.8 1.0 1.2 1.4 [sec]