Input to OUTMOD
1. General Information
The postprocessing module OUTMOD
has two main purposes:

Generate result printout from the
INPMOD
,STAMOD
andDYNMOD
modules. 
Prepare a plot file (
IFNPLO
) for later use by the plot module (PLOMOD
). Note that this functionality is deprecated.
Locations for result presentation are identified by the line identifier
and segment and element numbers specified as input to INPMOD
for all single riser systems (i.e. SA, SB, SC, SD and AR systems).
The user chooses the amount of printout by giving the appropriate
options as input data to OUTMOD
. As for the other modules, input data
to OUTMOD
are organised in groups. Some groups consist of the data
group identifier only, whilst other groups have additional input lines.
The first part of the input data is data groups selecting data to be
printed. Then a command must be given to start the printing. Then a new
set of data groups selecting data may be given etc, see the figure
below.
2. Data Group A: OUTMOD Identification and Control Data
2.1. OUTMOD identification text
If you want an identification text to be printed on the front page of
the OUTMOD
printout, the following data group may be given as the
first data group:
2.2. The PRINT command
When OUTMOD
is used for print generation, the following data group
identifier must be given subsequent to the specifications:
After printing, all specifications will be deleted. Subsequent to the print command, a new sequence of specifications and a new print command may be given. This is useful if you want to repeat one or more specifications with different parameters.
2.3. Plot generation
This functionality is deprecated!
If you want to have plots from STAMOD
or DYNMOD
, you have to run
OUTMOD
first to produce the plot file, IFNPLO
, which is the only
file the plot module PLOMOD
reads when plotting from the above
mentioned modules.
It is not possible to generate plots of all data groups. This is marked
by the word Plot
or NoPlot
in the right part of the data group
identifier frame:
WFMOtion TIME SERies Plot
can be plotted, while the following can not:
INFIrr CONTRol INFOrmation NoPlot
When you want to use OUTMOD
to build up file IFNPLO
, the data group
identifiers and print options are exactly the same as when you use
OUTMOD
for normal printout, except for the following:

IFNPLO
must be initialized by one of the data group identifiers described inInitialization of the plot file IFNPLO
(below). 
Instead of the
PRINT
command, thePLOT
command must be given (see The PLOT command). TheOUTMOD
print file will contain the same information as if thePRINT
command was given, i.e. thePLOT
command may be considered as an extension to thePRINT
command.
Normally, one specification gives one plot which may consist of up to three graphs. In some cases, e.g. when one specification gives results for six degrees of freedoms, two plots are produced per specification.
2.3.1. Initialization of the plot file IFNPLO
When OUTMOD
is to be used for plot generation, one of the following
initialization commands must be given after specification of OUTMOD
identification text (if given), but before any specification described
in Sections Data Group B: Output from STAMOD, or Data Group C: Output from DYNMOD:
NEW PLOT FILE
or
APPEnd PLOT FILE
If the command NEW PLOT FILE
is given, OUTMOD
writes the plot arrays
from the beginning of the file, i.e. the previous contents of the file,
if any, are overwritten. It is, however, possible to append new plots
after already existing plots on file IFNPLO
. This is achieved by
giving the command APPEND PLOT FILE
. A check is carried out to ensure
that the file already contains plots.
A maximum of 100 plots may be stored in one IFNPLO
file.
If more than one initialization command is given throughout an OUTMOD
run, they are simply ignored.
2.3.2. The PLOT command
When OUTMOD
is used for plot generation, the following command is
given instead of the PRINT
command (see The PRINT command):
PLOT
If a plot command is given for a data group that cannot be plotted, a
warning message is issued on OUTMOD
print file and execution continues
with the next data group specified.
There may be more than one PLOT
within one run of OUTMOD
, following
the same rules as for the PRINT
command. PRINT
and PLOT
commands
may be mixed within one run.
If the PLOT
command is given, and neither NEW PLOT FILE
nor
APPEND PLOT FILE
has been given, the program will terminate with an
error message.
2.4. Communication with the STARTIMES programs
The following command has been included for communication with the
STARTIMES
programs for statistical analysis of time series
STARtimes FILE NoPlot
This command specifies that time series of a selected response quantity
shall be written to a file in STARTIMES
format (i.e. to a file
readable by the STARTIMES
programs). The STARTIMES FILE
command can
be used in connection with the following data group identifiers:

WAVE ELEVATION

WFMOTION TIME SERIES

LFMOTION TIME SERIES

TOMOTION TIME SERIES

DYNDISP TIME SERIES

TOTDISP TIME SERIES

DYNFORC TIME SERIES

DYNCURV TIME SERIES

SUPPF TIME SERIES

ELMANGLE TIME SERIES

TOTFORC TIME SERIES

DISTANCE TIME SERIES

CALCURV TIME SERIES

STRESS TIME SERIES

STROKE TIME SERIES
A response quantity is written to the STARTIMES
file by giving
STARTIMES FILE
immediately before the PRINT
or PLOT
command.
The name of the STARTIMES
file is <prefix>_outmod.ts
, and it is
stored in the current working directory. A description of this file is
found in Description of STARTIMES file.
3. Data Group B: Output from STAMOD
Description of result presentation from static analyses is given in the following.
3.1. Results from static fixed parameter analysis
Displacement and force data from static fixed parameter analysis are
established by the STAMOD
module and stored on file IFNSTA
.
Specifying this output after a parameter variation run will produce the results of the last parameter variation step.
3.1.1. Static dimension information
If you want dimension parameters, such as no of load steps, no of nodes etc, to be printed, give
STATic DIMENsion PARAmeters NoPlot
3.1.2. System information
If you want information about the connection between the local line,
segment and element number given as input to INPMOD
and the global
FEM
element/nodal numbers generated by STAMOD
, give
STATic SYSTem INFOrmation NoPlot
A more detailed description is given on the STAMOD
print file.
3.1.3. Coordinates of final static configuration
Print options, one input line
ICONF LINEID IPROJ

CONF: integer
: Configuration switch
ICONF=1
: Initial configuration (catenary configuration) 
ICONF=2
: Final configuration (Results from FEM or CATFEM analysis)


LINEID: character(8)
: Line identifier for which coordinates are wanted. You may specifyALL
to include all lines in the system.
Note that specifying a specific line gives a 2Dplot, while specifying
ALL
gives a 3Dplot 
LINEID=0
: Plot of 2D geometry of all lines


IPROJ: integer
: Projection code
dummy if
LINEID=ALL

IPROJ=1
: Output of xy coordinates 
IPROJ=2
: Output of xz coordinates 
IPROJ=3
: Output of yz coordinates

3.1.4. Axial forces from catenary analysis
Note that no moments are included in the catenary analysis.
3.1.5. Forces from static fixed parameter analysis
Forces are printed as force, bending and torsional moments.
Data group identifier, one input line
FINAl STATic FORCes Plot
Pipe wall force, axial:
\(\mathrm {T_W=T_e+p_iA_ip_eA_e[+m_iv_i^2]}\)
(In cases with high pressure(s) it may be important to include the radial stress when material strain is to be evaluated)
This is identical with the flange force in case of a double seal (at \(\mathrm {r=ri}\) and \(\mathrm {r=re}\))

\(\mathrm {T_F=T_W}\)
In the case of an inner seal only:

\(\mathrm {T_F=T_e+p_iA_ip_eA_i[+m_iv_i^2]}\)
Any other sealing radius:

\(\mathrm {T_F=T_e+(p_ip_e)A_s[+m_iv_i^2]}\)
Where:  \(\mathrm {A_s=\pi r_s^2}\)  \(\mathrm {r_s=}\) sealing radius
rs = sealing radius
\(\mathrm {m_iv_i^2}\) is an additional term for cases with internal fluid flow.
Print options, one input line
LINEID IDOF1 IDOF2 IDOF3 IDOF4

LINEID: character(8)
: Line identifier for which forces are wanted. You may specifyALL
to include all lines in the system. 
IDOF1: integer
: Degree of freedom for first figure
IDOF1=0
: Not included 
IDOF1=1
: Axial force 
IDOF1=2
: Torsional moment 
IDOF1=3
: Bending moment about local yaxis 
IDOF1=4
: Bending moment about local zaxis 
IDOF1=5
: Pipe wall force, incl hydrostatic pressures 
IDOF1=6
: Shear force in local ydirection 
IDOF1=7
: Shear force in local zdirection


IDOF2: integer
: Degree of freedom for second figure
Interpretation as for
IDOF1


IDOF3: integer
: Degree of freedom for third figure
Interpretation as for
IDOF1


IDOF4: integer
: Degree of freedom for fourth figure
Interpretation as for
IDOF1

No of figures in one plot may vary from 13 depending on the number of response quantities specified (e.g. IDOFi).
Note that the print part of this option always will produce results for all stored degrees of freedom, i.e. axial force, torsional moment and bending moments about local y and zaxes. The parameters are used to specify the dof’s to be plotted.
3.1.6. Stress from static analysis
Output options, one input line
LINEID IDOF

LINEID: character(8)
: Line identifier for which stresses are wanted. You may specifyALL
to include all lines in the system. 
The following parameter is used to specify the dof to be considered

IDOF: integer
: Stress component
IDOF=1
: Axial stress 
IDOF=2
: Torsional stress 
IDOF=3
: Bending stress 
IDOF=4
: Axial + bending stress 
IDOF=5
: Shear stress 
IDOF=6
: Shear + torsional stress 
IDOF=7
: Equivalent stress 
IDOF=8
: Hoop stress 
IDOF=9
: Radial stress


Effect of internal/external pressure and fluid velocity are included
Specification of point for stress calculation, one input line
IMAX THETA INEX

IMAX: integer, default: 1
: Stress location option
IMAX=1
: Maximum stresses in cross section estimated 
IMAX=0
: Stresses calculated at location specified byTHETA
andINEX


THETA: real, default: 0
: Angle (in degrees) from local yaxis for stress calculation.
Dummy for
IMAX=1


INEX: integer, default: 2
: Location code
Dummy for
IMAX=1

INEX=1
: Inner wall 
INEX=2
: Outer wall

For IMAX=1
, the maximum stresses of type IDOF
in the cross section
are estimated. The equivalent stress (von Mises) is supposed to be
maximum where the bending stress is maximum or minimum.
3.2. Output from static parameter variation analysis
Displacement and force data from static parameter variation analysis are
established by the STAMOD
module and stored on file IFNSTA
. Result
presentation from static parameter variation analysis is described in
the following.
3.2.1. System geometry from parameter variation analysis
Line specification, one input line
LINEID IOTYP IPV1 NVP

LINEID: character(8)
: Line number for which geometry are wanted. You may specifyALL
to include all lines in the system.
ALL
gives a 3D plot of all lines. 
LINEID = 0
gives a 2D plot of all lines


IOTYP: integer
: Degree of freedom specification
Dummy if
ILINE = ALL

IOTYP=1
: xy coordinates 
IOTYP=2
: xz coordinates 
IOTYP=3
: yz coordinates


IPV1: integer
: First parameter variation step to be included 
NVP: integer
: No of parameter variation steps to be included
The first plot to appear will be for step no NSTEP+IPV1
where NSTEP
is total number of load steps used in the static analysis with fixed
parameters.
Negative value of IPV1
is possible, which allows for plotting of
static configuration at all load steps in static analysis with fixed
parameters.
It is also possible to plot static configurations from 1st load step to last successful solution when static analysis fails, which can be very useful for detection of possible instability problems.
3.2.2. Displacement of selected nodes from parameters variation analysis
Output code, one input line
IPV1 NPV IDOF1 IDOF2 IDOF3 NNODC

IPV1: integer
: First parameter variation step to be included 
NPV: integer
: No of parameter load steps to be included. (A large number includes the remaining steps) 
IDOF1: integer
:
IDOF1=1
: Translation in xdirection 
IDOF1=2
: Translation in ydirection 
IDOF1=3
: Translation in zdirection


IDOF2: integer
:
Interpretation as for
IDOF1


IDOF3: integer
:
Interpretation as for
IDOF1


NNODC: integer
: No. of input lines used for node specification
No of figures on each plot may vary from 1 to 3, depending on IDOFi
The first plot to appear will be for step no NSTEP+IPV1
where NSTEP
is the total number of load steps in the static analysis with fixed
parameters.
Negative value of IPV1
is allowed (see System geometry from parameter variation analysis).).
Node specification, NNODC input lines
LINEID ISEG INODE

LINEID: character(8)
: Line identifier.
You may specify
ALL
to include all lines


ISEG: integer/character
: Segment number.
You may specify
ALL
to include all segments. 
ENDS
includes the end segments on the line


INODE: integer/character
: Node number.
ALL
includes all nodes 
ENDS
includes end nodes on the above specified segment

3.2.3. Forces on selected elements from parameter variation analysis
Output options, one input line
IPV1 NPV IDOF1 IDOF2 IDOF3 NNELC

IPV1: integer
: First parameter variation step to be included 
NPV: integer
: No of parameter load steps to be included. (A large number includes the remaining steps) 
IDOF1: integer
: Degree of freedom specification for Figure 1
IDOF1=0
: No output 
IDOF1=1
: Axial force 
IDOF1=2
: Torsional moment 
IDOF1=3
: Bending moment about local yaxis 
IDOF1=4
: Bending moment about local zaxis


IDOF2: integer
:
Interpretation as for
IDOF1


IDOF3: integer
:
Interpretation as for
IDOF1


NNELC: integer
: No. of input lines used for element specification
The first plot to appear will be for step no NSTEP+IPV1
where NSTEP
is the total number of load steps used in the static analysis with fixed
parameters.
Negative value of IPV1
is allowed (see System geometry from parameter variation analysis).
No of figures in one plot may vary from 1 to 3, depending on the number
of response quantities specified (e.g. IDOFi
).
Element specification, NNELC input lines
LINEID ISEG IELM

LINEID: character(8)
: Line identifier.
You may specify
ALL
to include all lines


ISEG: integer/character
: Segment number.
You may specify
ALL
to include all segments. 
ENDS
includes the end segments on the line


IELM: integer/character
: Element number.
ALL
includes all elements 
ENDS
includes end elements on the above specified segment

4. Data Group C: Output from DYNMOD
4.1. Results from irregular wave analysis
Results from the irregular wave analysis consists of:  sampled Fourier
components of waves stored on file IFNIRR
at global origin, x=y=z=0
 motion of the support vessel, stored on file IFNIRR
 motion
transfer functions for the support vessel
4.1.2. Sampled Fourier components
Output parameters, one input line
ICOMP IDIR ISEC IW1 NW IJP

ICOMP: integer
: Component code
ICOMP=1
: Wind sea 
ICOMP=2
: Swell


IDIR: integer
: Direction no wanted 
ISEC: integer
: Sequence no wanted (dummy) 
IW1: integer
: Number of the first frequency for which Fourier components are wanted 
NW: integer
: No of frequencies for which Fourier components are wanted 
IJP: integer, default: 1
: Jump parameter
Fourier components are printed for frequencies no IW1
, IW1+IJP
,
IW1
+2x`IJP`, …, IW1
+(NW
1)x`IJP`
The components are printed/plotted as amplitude and phase angle (degrees)
4.1.3. Wave elevation
Output parameters, one input line
ICOMP IDIR ISEC IT1 NTS XP1 XP2

ICOMP: integer
: Component code
ICOMP=1
: Wind sea 
ICOMP=2
: Swell


IDIR: integer
: Direction no wanted 
ISEC: integer
: Sequence no wanted (dummy) 
IT1: integer
: First time step included 
NTS: integer
: Number of time steps included 
XP1: real, default: 0
: Global xcoordinate for wave elevation 
XP2: real, default: 0
: Global ycoordinate for wave elevation
A Fourier transformation of the wave spectrum is performed. Maximum
number of time steps will be (NWIMAX1)*2. Use the option
IFNIRR CONTROL INFORMATION
(see Control information).
In case of longcrested sea one direction is applied. In case of shortcrested sea, 11 directions are used and mean wave direction is no. 6. The other directions are spread around the mean direction in the interval \(\mathrm {[75^{\circ},75^{\circ}]}\) in intervals of \(\mathrm {15^{\circ}}\).
4.1.4. Wave frequency motion time series
Output options, one input line
IOP IMOT IDERIV ISEQ1 NSEQ IT1 NTS ITJMP IVES

IOP: integer
: Code for type of output
IOP=1
: Time series 
IOP=2
: Time series statistics 
IOP=3
: Spectral analysis


IMOT: integer
: Direction
IMOT=1
: Displacement in global xdirection 
IMOT=2
: Displacement in global ydirection 
IMOT=3
: Displacement in global zdirection 
IMOT=4
: Rotation about xaxis 
IMOT=5
: Rotation about yaxis 
IMOT=6
: Rotation about zaxis


IDERIV: integer
: Code for derivative of response
IDERIV=0
: Analyse original series 
IDERIV=1
: Analyse 1st derivative 
IDERIV=2
: Analyse 2nd derivative


ISEQ1: integer
: First sequence to be included (dummy) 
NSEQ: integer
: No of sequence to be included (dummy) 
IT1: integer
: First time step of each sequence to be included 
NTS: integer
: No of time steps of each sequence to be included 
ITJMP: integer, default: 1
: Jump parameter
Time step nos.
IT1, IT1+ITJMP, IT1+2xITJMP,…, IT1+(NTS1)xITJMP
are included


IVES: integer, default: 1
: Vessel number reference in case of multivessel systems.The vessels are numbered from 1 toNVES
Note that IMOT
refers to the global coordinate system, not the vessel
coordinate system.
Transformation of wave frequency motion time series, one input line
ITRANS XP YP ZP

ITRANS: integer, default: 0
: Transformation code
ITRANS=0
: No transformation, motions of vessel reference point 
ITRANS=1
: Transformation gives motionIMOT
(see previous input line) of point defined byXP
,YP
andZP


XP: real, default: 0
: Xcoordinate in global system, relative to the vessel reference point 
YP: real, default: 0
: Ycoordinate in global system, relative to the vessel reference point 
ZP: real, default: 0
: Zcoordinate in global system, relative to the vessel reference point
If ITRANS=0
, XP
, YP
and ZP
are dummy parameters
Options for the output distribution functions of the high frequency motion time series statistics, one input line
This input line is given only if IOP=2
.
NCL XCMIN XCMAX

NCL: integer
: No of classes in the output distribution functions (i.e. no of points on the abscissa axis)
0<NCL<41


XCMIN: real
: Range of argument values for output distribution functions isXCMIN*sx(1)  XCMAX*sx(1)
in whichsx(1)
is the standard deviation ofx
estimated from the first sequence. 
XCMAX: real
:
4.1.5. Low frequency motion time series
Output options, one input line
IOP IMOT IDERIV ISEQ1 NSEQ IT1 NTS ITJMP IVES

IOP: integer
: Code for type of output
IOP=1
: Time series 
IOP=2
: Time series statistics 
IOP=3
: Spectral analysis


IMOT: integer
: Direction code
Legal values:

IMOT=1
: Surge 
IMOT=2
: Sway 
IMOT=6
: Yaw



IDERIV: integer
: Code for derivative of response
IDERIV=0
: Analyse original series 
IDERIV=1
: Analyse 1st derivative 
IDERIV=2
: Analyse 2nd derivative


ISEQ1: integer
: First sequence to be included (dummy) 
NSEQ: integer
: No of sequence to be included (dummy) 
IT1: integer
: First time step of each sequence to be included 
NTS: integer
: No of time steps of each sequence to be included 
ITJMP: integer, default: 1
: Jump parameter
Time step nos.
IT1, IT1+ITJMP, IT1+2xITJMP,…, IT1+(NTS1)xITJMP
are included


IVES: integer, default: 1
: Vessel number reference in case of multivessel systems.The vessels are numbered from 1 toNVES
Options for the output distribution functions of the low frequency motion time series statistics, one input line
This input line is given only if IOP=2
.
4.1.6. Total motion time series
Output options, one input line
IOP IMOT IDERIV ISEQ1 NSEQ IT1 NTS ITJMP IVES

IOP: integer
: Code for type of output
IOP=1
: Time series 
IOP=2
: Time series statistics 
IOP=3
: Spectral analysis


IMOT: integer
: Direction code
IMOT=1
: Displacement in global xdirection 
IMOT=2
: Displacement in global ydirection 
IMOT=6
: Rotation about zaxis


IDERIV: integer
: Code for derivative of response
IDERIV=0
: Analyse original series 
IDERIV=1
: Analyse 1st derivative 
IDERIV=2
: Analyse 2nd derivative


SEQ1: integer
: First sequence to be included (dummy) 
NSEQ: integer
: No of sequence to be included (dummy) 
IT1: integer
: First time step of each sequence to be included 
NTS: integer
: No of time steps of each sequence to be included 
ITJMP: integer, default: 1
: Jump parameter
Time step nos.
IT1, IT1+ITJMP, IT1+2xITJMP,…, IT1+(NTS1)xITJMP
are included


IVES: integer, default: 1
: Vessel number reference in case of multivessel systems.The vessels are numbered from 1 toNVES
Options for the output distribution functions of the time series statistics of total motion, one input line
This input line is given only if IOP=2
.
4.1.7. Vessel motion transfer functions
group identifier, one input line
WFTRansfer FUNCtion DOF Plot
DOF
means degree of freedom, and may be XG
, YG
, ZG
, XGROT
,
YGROT
or ZGROT
.
Output options, one input line
IOP IDIR1 NDIR ITRAN IVES

IOP: integer
: Code for type of output
IOP=1
: Complex form (real, imaginary) 
IOP=2
: Real form (amplitude ratio, phase (degrees)) 
IOP=3
: Real form (amplitude ratio, phase (radians))


IDIR1: integer
: First direction to be included 
INDIR: integer
: No of directions to be included 
ITRAN: integer
: Code for transformation
ITRAN=0
: No transformation 
ITRAN=1
: Transformation of origin motion to point (XV1
,XV2
,XV3
), see next input line. 
Dummy if degree of freedom is
XGROT
,YGROT
orZGROT


IVES: integer, defaul: 1
: Vessel number
The coordinates of the point on the vessel for which the vessel motion transfer functions are wanted, one input line
If ITRAN=0
, or the degree of freedom is XGROT
, YGROT
or ZGROT
,
this input line is skipped.
XV1 XV2 XV3

XV1: real
: Xcoordinate of the point 
XV2: real
: Ycoordinate of the point 
XV3: real
: Zcoordinate of the point
The coordinates are referred to the global coordinate system, relative to the vessel reference point.
The transfer functions for different degrees of freedom may be given
without the PRINT
or PLOT
statement between.
4.2. Results from time domain dynamic analysis
4.2.1. Storage information
Print options, one input line
IDNOD IFNOD ICNOD

IDNOD: integer, default: 1
: Switch for printing of nodes for which displacements are stored
IDNOD=0
: No print 
IDNOD=1
: The nodes, for which displacements are stored, are printed


IFNOD: integer, default: 1
: Switch for printing of elements for which force data are stored
IFNOD=0
: No print 
IFNOD=1
: The nodes, for which force data are stored, are printed


ICNOD: integer, default: 1
: Switch for printing of elements for which curvature data are stored
ICNOD=0
: No print 
ICNOD=1
: The elements, for which curvature data are stored, are printed

4.2.2. Snapshot plot from time domain analysis
This option will create pictures of the dynamic configuration at several time steps.
Plot options
IPROJ IT1 NTS NLIC IJUMP

IPROJ: integer
: Project in code
IPROJ=1
: xz coordinates 
IPROJ=2
: yz coordinates 
IPROJ=3
: xy coordinates


IT1: integer
: First stored time step to be included 
NTS: integer/character
: No of stored time steps to be included.
You may specify
REST
to include the remaining time steps


NLIC: integer
: No. of input lines to describe line specification 
IJUMP: integer, default: 1
: Plot everyIJUMP
stored time step
4.2.3. System snapshot plot from time domain analysis
This option is an extension to the option DYNAMIC SNAPSHOT PLOT
. You
are able to plot the wave particle motion, the vessel motion and the
riser motion in one plot.
Plot options, one input line
IPROJ IT1 NTS IJUMP NLIC NPVESP NPWAPO IVES XCGVES YCGVES ZCGVES

IPROJ: integer
: Projection code
IPROJ=1
: XZ coordinates 
IPROJ=2
: YZ coordinates 
IPROJ=3
: XY coordinates


IT1: integer
: First stored time step to be included 
NTS: integer
: No of stored time steps to be included. You may specify REST to include the remaining time steps 
IJUMP: integer
: Include everyIJUMP
stored time steps 
NLIC: integer
: No. of input lines to describe line specification
NLIC=0
: No riser snapshot plot


NPVESP: integer
: No of coordinates to describe the vessel
NPVESP=0
: No vessel snapshot plot


NPWAPO: integer
: No of coordinates to describe the wave particle motion
NPWAPO=0
: No wave particle snapshot plot


IVES: integer, default: 1
: Vessel number 
XCGVES: real
: Static X coordinate of the vessel 
YCGVES: real
: Static Y coordinate of the vessel 
ZCGVES: real
: Static Z coordinate of the vessel
Line specification, NLIC input lines
LINEID

LINEID: integer/character(8)
: Line identifier to be plotted. You may specifyALL
to include all lines in the system
The lines are plotted only if at least the end node coordinates are stored.
Vessel description, NPVESP input lines. The specified points are connected by one line to illustrate a part of the vessel contour
IPV XVT YVT ZVT

IPV: integer
: Coordinate no. 
XVT: real
: Vessels Xcoordinate in global system referred from vessel origin \(\mathrm {[L]}\) 
YVT: real
: Vessels Ycoordinate \(\mathrm {[L]}\) 
ZVT: real
: Vessels Zcoordinate \(\mathrm {[L]}\)
Wave particle description, NPWAPO input lines
IPW XPW YPW ZPW

IPW: integer
: Coordinate no.
If
IPW<0
, then the intermediate coordinates between the previous coordinate specification and this one are automatically calculated. The intermediate coordinates are equally spaced on a straight line


XPW: real
: Xcoordinate of the wave particle \(\mathrm {[L]}\) 
YPW: real
: Ycoordinate of the wave particle \(\mathrm {[L]}\) 
ZPW: real
: Zcoordinate of the wave particle \(\mathrm {[L]}\)
The wave particle coordinates are given in the global coordinate system
in calm water, i.e. (0.,0.,0.) is wave at global origin. Specifying
ZPW
\(\mathrm {\equiv}\) 0. for all points will create a
plot of the wave surface elevation.
4.2.4. Dynamic displacement time series from time domain analysis
Results include only the dynamic time dependant displacements (static values are not included).
Output options, one input line
IOP IDOF IT1 NTS NNODC Plot

IOP: integer
: Code for type of output
IOP=1
: Time series 
IOP=2
: Time series statistics 
IOP=3
: Spectral analysis


IDOF: integer
: Code for degree of freedom
Rotational degrees of freedom are only to be presented from linearized dynamic analysis.

IDOF=1
: Translation in xdirection 
IDOF=2
: Translation in ydirection 
IDOF=3
: Translation in zdirection 
IDOF=4
: Rotation about xaxis 
IDOF=5
: Rotation about yaxis 
IDOF=6
: Rotation about zaxis


IT1: integer
: First stored time step to be included 
NTS: integer
: Number of stored time steps to be included (fromIT1
).
A large number includes the remaining time steps


NNODC: integer
: No. of input lines used for node specification
For IOP=3
an FFT analysis is carried out. If NTS
is not an integer
power of 2, a reduced time series will be analysed. In order to get an
effective analysis, IT1
and NTS
should be selected so that 
\(\mathrm {IT1=NT2^M+1}\)  \(\mathrm {NTS=2^M}\)
Where \(\mathrm {NT}\) is the total number of stored time steps and \(\mathrm {M}\) is the largest integer so that \(\mathrm {NTS<=NT}\). Normally it is preferable to omit the first part of the time series due to transients
Node specification, NNODC input lines
LINEID ISEG INODE

LINEID: character(8)
: Line identifier.
You may specify
ALL
to include all lines


ISEG: integer/character
: Segment number.
You may specify
ALL
to include all segments. 
ENDS
includes the end segments on the line


INODE: integer/character
: Node number.
ALL
includes all nodes 
ENDS
includes end nodes on the above specified segment

Displacements are not necessarily stored for all nodes, see data group File storage of displacement response for storage information. If the user specifies nodes for which displacements are not stored, these nodes are ignored.
The data group Storage information may be used to obtain an overview of the stored data.
Options for the output distribution functions of the displacement time series statistics, one input line
This input line is given only if IOP=2
.
NCL XCMIN XCMAX

NCL: integer
: No of classes in the output distribution functions (i.e. no of points on the abscissa axis)
0<NCL<41


XCMIN: real
: Range of argument values for output distribution functions isXCMIN*sx(1)  XCMAX*sx(1)
in whichsx(1)
is the standard deviation ofx
estimated from the first sequence 
XCMAX: real
: See above
4.2.5. Dynamic resulting force time series from time domain analysis
The results include only the dynamic time dependent force. Static values are not included.
Output options, one input line
IOP IDOF IT1 NTS NNELC

IOP: integer
: Code for type of output
IOP=1
: Time series 
IOP=2
: Time series statistics 
IOP=3
: Spectral analysis


IDOF: integer
: Code for degree of freedom
IDOF=1
: Axial force 
IDOF=2
: Torsional moment 
IDOF=3
: Bending moment about local yaxis, end 1 
IDOF=4
: Bending moment about local yaxis, end 2 
IDOF=5
: Bending moment about local zaxis, end 1 
IDOF=6
: Bending moment about local zaxis, end 2 
IDOF=7
: Shear force in local ydirection, end 1 
IDOF=8
: Shear force in local ydirection, end 2 
IDOF=9
: Shear force in local zdirection, end 1 
IDOF=10
: Shear force in local zdirection, end 2


IT1: integer
: First stored time step to be included 
NTS: integer
: Number of stored time steps to be included (fromIT1
).
A large number includes the remaining time steps.


NNELC: integer
: No. of input lines used for element specification
For IOP=3
an FFT analysis is carried out. If NTS
is not an integer
power of 2, a reduced time series will be analysed. In order to get an
effective analysis, IT1
and NTS
should be selected so that 
\(\mathrm {IT1=NT2^M+1}\)  \(\mathrm {NTS=2^M}\)
Where \(\mathrm {NT}\) is the total number of stored time steps and \(\mathrm {M}\) is the largest integer so that \(\mathrm {NTS<=NT}\). Normally it is preferable to omit the first part of the time series due to transients
Element specification, NNELC input lines
LINEID ISEG IELM

LINEID: character(8)
: Line identifier.
You may specify
ALL
to include all lines


ISEG: integer/character
: Segment number.
You may specify
ALL
to include all segments. 
ENDS
includes the end segments on the line


IELM: integer/character
: Element number.
ALL
includes all Elements 
ENDS
includes end elements on the above specified segment

Forces are not necessarily stored for all elements, see data group File storage for internal forces for storage information. If the user specifies elements for which forces are not stored these elements are ignored.
The data group Storage information may be used to obtain an overview of the stored data.
Options for the output distribution functions of the force time series statistics, one input line
This input line is given only if IOP=2
.
NCL XCMIN XCMAX

NCL: integer
: No of classes in the output distribution functions (i.e. no of points on the abscissa axis)
0<NCL<41


XCMIN: real
: Range of argument values for output distribution functions isXCMIN*sx(1)  XCMAX*sx(1)
in whichsx(1)
is the standard deviation ofx
estimated from the first sequence 
XCMAX: real
:
4.2.6. Curvature time series from time domain analysis
Results include only the dynamic time dependant curvature (static values are not included)
See also data group Curvature time series calculated from dynamic nodal displacements.
Output options, one input line
IOP IDOF IT1 NTS NNELC

IOP: integer
: Code for type of output
IOP=1
: Time series 
IOP=2
: Time series statistics 
IOP=3
: Spectral analysis


IDOF: integer
: Code for degree of freedom
IDOF=1
: Curvature about local yaxis, end 1 
IDOF=2
: Curvature about local yaxis, end 2 
IDOF=3
: Curvature about local zaxis, end 1 
IDOF=4
: Curvature about local zaxis, end 2


IT1: integer
: First stored time step to be included 
NTS: integer
: Number of stored time steps to be included (fromIT1
).
A large number includes the remaining time steps.


NNELC: integer
: No. of input lines used for element specification
For IOP=3
an FFT analysis is carried out. If NTS
is not an integer
power of 2, a reduced time series will be analysed. In order to get an
effective analysis, IT1
and NTS
should be selected so that 
\(\mathrm {IT1=NT2^M+1}\)  \(\mathrm {NTS=2^M}\)
Where \(\mathrm {NT}\) is the total number of stored time steps and \(\mathrm {M}\) is the largest integer so that \(\mathrm {NTS<=NT}\). Normally it is preferable to omit the first part of the time series due to transients
Options for the output distribution functions of the curvature time series statistics, one input line
This input line is given only if IOP=2
.
4.2.7. Curvature time series calculated from dynamic nodal displacements
See also Curvature time series from time domain analysis for curvature component time series.
This option gives absolute value of curvature in 3D space at a specified node. Calculation of curvature is based on the interpolating polynomial through the positions of 3 adjacent nodes in the same line. Curvature can therefore only be calculated if displacement time series are stored for the specified node and two neighbouring nodes (see data group File storage of displacement response for storage information). The data group Storage information may be used to obtain an overview of the stored data.
Calculation of curvature at line ends is omitted.
Data group identifier, one input line
CALCurv TIME SERIes Plot
Total curvature calculated from the selected node and the two neighbouring nodes.
Output options, one input line
IOP IT1 NTS NNODC

IOP: integer
: Code for type of output
IOP=1
: Time series 
IOP=2
: Time series statistics 
IOP=3
: Spectral analysis


IT1: integer
: First stored time step to be included 
NTS: integer
: Number of stored time steps to be included (fromIT1
).
A large number includes the remaining time steps.


NNODC: integer
: No. of input lines used for element specification
For IOP=3
an FFT analysis is carried out. If NTS
is not an integer
power of 2, a reduced time series will be analysed. In order to get an
effective analysis, IT1
and NTS
should be selected so that 
\(\mathrm {IT1=NT2^M+1}\)  \(\mathrm {NTS=2^M}\)
Where \(\mathrm {NT}\) is the total number of stored time steps and \(\mathrm {M}\) is the largest integer so that \(\mathrm {NTS<=NT}\). Normally it is preferable to omit the first part of the time series due to transients.
Options for the output distribution functions of the curvature time series statistics, one input line
This input line is given only if IOP=2
.
4.2.8. Displacement envelope curves
Envelope curves of displacements from time domain analysis are presented as:  Minimum, static and maximum x, y and z displacements for regular analysis  Mean, static and mean + standard deviation for irregular analysis.
Static values are identified as dashed lines while the others are solid.
Print options, one input line
LINEID IPDOF1 IPDOF2 IPDOF3

LINEID: character(8)
: Line identifier for which displacements are wanted.
You may specify
ALL
to include all lines in the system. 
The print part of this option will always produce results for all stored degrees of freedom, i.e. x, y and zdisplacements. The following parameters are used to specify the dof’s to be plotted


IPDOF1: integer
: Degree of freedom for first figure
IPDOF1=0
: Not included 
IPDOF1=1
: xdisplacement 
IPDOF1=2
: ydisplacement 
IPDOF1=3
: zdisplacement


IPDOF2: integer
: Degree of freedom for second figure.
Interpretation as for
IPDOF1


IPDOF3: integer
: Degree of freedom for third figure.
Interpretation as for
IPDOF1

Each figure is presented on separate plot.
4.2.9. Force envelope curves
Envelope curves of forces from time domain analysis are presented as:

Minimum, static and maximum axial force torsional moment or bending moments for regular analysis

Mean, static and mean + standard deviation for irregular analysis
Static values are identified as dashed lines while the others are solid.
Print options, one input line
LINEID IDOF1 IDOF2 IDOF3

LINEID: character(8)
: Line identifier for which forces are wanted.
You may specify
ALL
to include all lines in the system. 
The print part of this option will always produce results for all stored degrees of freedom, i.e. axial force, torsional moment and bending moments about local y and zaxes. The following parameters are used to specify the dof’s to be plotted


IDOF1: integer
: Degree of freedom for first figure.
IDOF1=0
: Not included 
IDOF1=1
: Axial force 
IDOF1=2
: Torsional moment 
IDOF1=3
: Bending moment about local yaxis 
IDOF1=4
: Bending moment about local zaxis 
IDOF1=5
: Pipe wall force, incl. hydrostatic pressures
Pipe wall force is only avaivable for PLOT


IDOF1=6
: Shear force in local ydirection 
IDOF1=7
: Shear force in local zdirection


IPDOF2: integer
: Degree of freedom for second figure.
Interpretation as for
IPDOF1


IPDOF3: integer
: Degree of freedom for third figure.
Interpretation as for
IPDOF1

Each figure is presented on separate plot.
4.2.10. Curvature envelope curves
Envelope curves of curvatures from time domain analysis are presented as:  Minimum, static and maximum values of curvatures for a regular analysis  Mean, static and mean + standard deviation for irregular analysis
Static results are dashed, while the others are solid.
Print options, one input line
LINEID IDOF1 IDOF2 IDOF3

LINEID: character(8)
: Line identifier for which curvatures are wanted.
You may specify
ALL
to include all lines in the system. 
The print part of this option will always produce results for all stored degrees of freedom, i.e. local y and zcurvatures and resulting curvature. The following parameters are used to specify the dof’s to be plotted


IPDOF1: integer
: Degree of freedom for first figure
IDOF1=0
: Not included 
IDOF1=1
: Curvature about local yaxis 
IDOF1=2
: Curvature about local zaxis 
IDOF1=3
: Resulting curvature 
Resulting curvature is taken as the vector sum of the curvatures about local y and zaxis and will therefore always be positive


IPDOF2: integer
: Degree of freedom for second figure.
Interpretation as for
IPDOF1


IPDOF3: integer
: Degree of freedom for third figure.
Interpretation as for
IPDOF1

Each figure is presented on separate plot.
4.2.11. Support forces
Forces in both ends of specified lines are analyzed and presented in the global coordinate system. Forces due to static and dynamic loads are included. Forces due to hydrostatic pressures are not included, i.e. the axial component is the effective tension.
Output options, one input line
IOP IDOF IT1 NTS NLINC

IOP: integer
: Code for type of output
IOP=1
: Time series 
IOP=2
: Time series statistics 
IOP=3
: Spectral analysis


IDOF: integer
: Code for degree of freedom
IDOF=1
: Global xdirection 
IDOF=2
: Global ydirection 
IDOF=3
: Global zdirection


IT1: integer
: First stored time step to be included 
NTS: integer
: Number of stored time steps to be included 
NLINC: integer
: Number of input lines used for line specifications
For IOP=3
an FFT analysis is carried out. If NTS
is not an integer
power of 2, a reduced time series will be analysed. In order to get an
effective analysis, IT1
and NTS
should be selected so that 
\(\mathrm {IT1=NT2^M+1}\)  \(\mathrm {NTS=2^M}\)
Where \(\mathrm {NT}\) is the total number of stored time steps and \(\mathrm {M}\) is the largest integer so that \(\mathrm {NTS<=NT}\). Normally it is preferable to omit the first part of the time series due to transients.
Element specification, NLINC input lines
LINEID

LINEID: character(8), default: 0
: Line number. You may specifyALL
to include all lines
Options for the output distribution functions of the force time series statistics, one input line
This input line is given only if IOP=2
.
NCL XCMIN XCMAX

NCL: integer
: No of classes in the output distribution functions (i.e. no of points on the abscissa axis)
0<NCL<41


XCMIN: real
: Range of argument values for output distribution functions isXCMIN*sx(1)  XCMAX*sx(1)
in whichsx(1)
is the standard deviation ofx
estimated from the first sequence 
XCMAX: real, default: 0
:
4.2.12. Element angle time series from time domain analysis
Output options, one input line
IOP IT1 NTS NNELC

IOP: integer
: Code for type of output
IOP=1
: Time series 
IOP=2
: Time series statistics 
IOP=3
: Spectral analysis


IT1: integer
: First stored time step to be included 
NTS: integer
: Number of stored time steps to be included (fromIT1
).
A large number includes the remaining time steps


NNELC: integer
: No. of pairs of input lines used for element specification
Two of the subsequent input lines (Code for element specification and either
Global or vessel axis and element
specification or Element pair specification) given NNELC
times.
Code for element specification
IRELCO

IRELCO: integer
: Code for type of output
IRELCO=0
: Angle between fixed global axis and one specified element 
IRELCO=1
: Angle between support vessel coordinate axis and one specified element 
IRELCO=2
: Angle between two elements

Global or vessel axis and element specification
This input line is given only for IRELCO=0
or 1
.
IAXIS IVES LINEID ISEG IELM HEAD

IAXIS: integer
: Code for axis
IAXIS=1
: xaxis 
IAXIS=2
: yaxis 
IAXIS=3
: zaxis


IVES: integer, default: 1
: Vessel number ifIRELCO=1
else dummy 
LINEID: character(8)
: Line identifier 
ISEG: integer
: Segment number 
IELM: integer
: Element number 
HEAD: integer
: Vessel heading in final static position ifIRECLCO=1
, else dummy \(\mathrm {[deg]}\)
The angle output is between 0 and 180 degrees. If the element direction (from end 1 to end 2) is along the specified axis, the relative angle is 0. Otherwise, if the element direction is along the negative axis direction, the angle is 180 degrees. The element direction is calculated as the direction along the secant from local end no 1 to local end no 2 (i.e. local element xaxis).
Element pair specification
This input line is given only for IRELCO=2
.
LINEID1 ISEG1 IELM1 LINEID2 ISEG2 IELM2

LINEID1: character(8)
: Specification of first element 
ISEG1: integer
: 
IELM1: integer
: 
LINEID2: character(8)
: Specification of second element 
ISEG2: integer
: 
IELM2: integer
:
The angle output is between 0 and 180 degrees. If the element direction (from end 1 to end 2) is along the specified axis, the relative angle is 0. Otherwise, if the element direction is along the negative axis direction, the angle is 180 degrees. The element direction is calculated as the direction along the secant from local end no 1 to local end no 2 (i.e. local element xaxis).
4.2.13. Total displacement time series from time domain analysis
Results include the total dynamic displacements (static values are included)
Output options, one input line
IOP IDOF IT1 NTS NNODC

IOP: integer
: Code for type of output
IOP=1
: Time series 
IOP=2
: Time series statistics 
IOP=3
: Spectral analysis


IDOF: integer
: Code for degree of freedom
IDOF=1
: Translation in xdirection 
IDOF=2
: Translation in ydirection 
IDOF=3
: Translation in zdirection


IT1: integer
: First stored time step to be included 
NTS: integer
: Number of stored time steps to be included (from IT1).
A large number includes the remaining time steps.


NNODC: integer
: No of input lines used for node specification
For IOP=3
an FFT analysis is carried out. If NTS
is not an integer
power of 2, a reduced time series will be analysed. In order to get an
effective analysis, IT1
and NTS
should be selected so that 
\(\mathrm {IT1=NT2^M+1}\)  \(\mathrm {NTS=2^M}\)
Where \(\mathrm {NT}\) is the total number of stored time steps and \(\mathrm {M}\) is the largest integer so that \(\mathrm {NTS<=NT}\). Normally it is preferable to omit the first part of the time series due to transients.
Options for the output distribution functions of the displacement time series statistics, one input line
This input line is given only if IOP=2
.
4.2.14. Total resulting force time series from time domain analysis
The result force includes both the dynamic time dependent force and the static force.
Output options, one input line
IOP IDOF IT1 NTS NNELC

IOP: integer
: Code for type of output
IOP=1
: Time series 
IOP=2
: Time series statistics 
IOP=3
: Spectral analysis


IDOF: integer
: Code for degree of freedom
IDOF=1
: Axial force 
IDOF=2
: Torsional moment 
IDOF=3
: Bending moment about local yaxis, end 1 
IDOF=4
: Bending moment about local yaxis, end 2 
IDOF=5
: Bending moment about local zaxis, end 1 
IDOF=6
: Bending moment about local zaxis, end 2 
IDOF=7
: Shear force in local ydirection, end 1
Nonlinear dynamic analysis only in present version


IDOF=8
: Shear force in local ydirection, end 2
Nonlinear dynamic analysis only in present version


IDOF=9
: Shear force in local zdirection, end 1
Nonlinear dynamic analysis only in present version


IDOF=10
: Shear force in local zdirection, end 2
Nonlinear dynamic analysis only in present version


IDOF=11
: Axial wall force


IT1: integer
: First stored time step to be included 
NTS: integer
: Number of stored time steps to be included (fromIT1
).
A large number includes the remaining time steps


NNELC: integer
: No of input lines used for element specification
For IOP=3
an FFT analysis is carried out. If NTS
is not an integer
power of 2, a reduced time series will be analysed. In order to get an
effective analysis, IT1
and NTS
should be selected so that 
\(\mathrm {IT1=NT2^M+1}\)  \(\mathrm {NTS=2^M}\)
Where \(\mathrm {NT}\) is the total number of stored time steps and \(\mathrm {M}\) is the largest integer so that \(\mathrm {NTS<=NT}\). Normally it is preferable to omit the first part of the time series due to transients.
Options for output distribution functions. Given only if IOP=2
This input line is given only if IOP=2
.
4.2.15. Distance time series calculated from the time domain analyses
This option is mainly to be used in order to perform a check of collision risk between two risers, between a riser and the vessel or between a riser and a fixed structure. The minimum distance is calculated for only a part of the riser. All elements within the specified segments are searched to find this minimum distance at each time step.
Note that the distances are absolute, they are always positive values. The program cannot identify a line crossing situation.
Output options, one input line
IOP IT1 NTS IDITYP IMETH IVES XCGVES YCGVES ZCGVES

IOP: integer
: Code for type of output
IOP=1
: Time series 
IOP=2
: Time series statistics 
IOP=3
: Spectral analysis


IT1: integer
: First stored time step to be included 
NTS: integer
: No of stored time steps to be included (fromIT1
).
A large number includes the remaining time steps


IDITYP: integer
: Type of distance to be calculated
IDITYP=1
: Distance between specified segments on lines 
IDITYP=2
: Distance between specified segments on a line and a globally fixed line 
IDITYP=3
: Distance between specified segments on a line and a line fixed on the vessel


IMETH: integer, default: 1
: Method option
IMETH=1
: Distance between elements are calculated 
IMETH=2
: Distance between nodes are calculated


IVES: integer, default: 1
: Vessel number in case of multivessel analysis andIDITYP = 3

XCGVES: real, default: 0
: Static X coordinate of the vessel in case ofIDITYP = 3

YCGVES: real, default: 0
: Static Y coordinate of the vessel in case ofIDITYP = 3

ZCGVES: real, default: 0
: Static Z coordinate of the vessel in case ofIDITYP = 3
With the distance, we here mean the minimum distance. All elements within the specified segment(s) are scanned for each time step in order to find the one with the minimum distance.
Method 1 is more accurate, but more time consuming than method 2.
Specification of segments on lines which the minimum distance should be calculated from, one input line
LINEID NSEG ISEG1 ISEG2 . . ISEG(NSEG)

LINEID: character(8)
: Line identifier 
NSEG: integer/character
: No of segments for which the minimum distances are to be calculated from
You may specify
ALL
in order to include all segments


ISEG: integer
: The included segment numbers
Searching through all elements may cause rather large computation time.
Specified segments to which the minimum distance are calculated, to be given only if IDITYP=1. One input line
LINEID NSEG ISEG1 ISEG2 ... ISEG(NSEG)

LINEID: character(8)
: Line identifier 
NSEG: integer/character
: No of segments for which the minimum distances are to be calculated to
You may specify
ALL
in order to include all segments


ISEGj: integer
: The included segment numbers
Searching through all elements may cause rather large computation time.
Specification of a line in the global coordinate system to which the minimum distance are to be calculated, to be given only if IDITYP=2. One input line
XG1 YG1 ZG1 XG2 YG2 ZG2

XG1: real
: Global xcoordinate, end 1 
YG1: real
: Global ycoordinate, end 1 
ZG1: real
: Global zcoordinate, end 1 
XG2: real
: Global xcoordinate, end 2 
YG2: real
: Global ycoordinate, end 2 
ZG2: real
: Global zcoordinate, end 2
Specification of a line in the global coordinate system relative to the vessel reference point to which the minimum distance are to be calculated, to be given only if IDITYP=3
XV1 YV1 ZV1 XV2 YV2 ZV2

XV1: real
: Vessel xcoordinate, end 1 
YV1: real
: Vessel ycoordinate, end 1 
ZV1: real
: Vessel zcoordinate, end 1 
XV2: real
: Vessel xcoordinate, end 2 
YV2: real
: Vessel ycoordinate, end 2 
ZV2: real
: Vessel zcoordinate, end 2
Options for the output distribution functions of the distance time series statistics, one input line
This input line is given only if IOP=2
.
NCL XCMIN XCMAX

NCL: integer
: No of classes in the output distribution functions (i.e. no of points on the abscissa axis)
0<NCL<41


XCMIN: real
: Range of argument values for output distribution functions isXCMIN*sx(1)  XCMAX*sx(1)
in whichsx(1)
is the standard deviation ofx
estimated from the first sequence 
XCMAX: real
:
4.2.16. Generate snapshot file from time domain analysis (special option)
This is a special option specified and commissioned by Norsk Hydro, for generation of input files for an animation program used by Norsk Hydro.
Nodes coordinates, element forces and curvatures from dynamic analysis
are written to the following files:  SNAPSNxx.DAT
 Node coordinates
 SNAPFOxx.DAT
 Element forces  SNAPCUxx.DAT
 Element curvatures
Element forces and/or curvatures will only be written for lines for which the storage coincide with the storage of node displacements.
Print options, one input line
IT1 NTS IJUMP NLIC NPVESD IVES LFORCE LCURV IASCII XCGVES YCGVES ZCGVES

IT1: integer
: First stored time step to be included 
NTS: integer
: Number of stored time steps to be included.
You may specify
REST
to include the remaining time step


IJUMP: integer
: Include every"IJUMP"
stored time step 
NLIC: integer
: No. of input lines to describe the line specification
NLIC=0
: No riser snapshot


NPVESD: integer
: No of coordinates to describe the vessel
NPVESD=0
: No vessel snapshot


IVES: integer, default: 1
: Vessel number in case of multivessel analysis 
LFORCE: integer, default: 0
: Control parameter
LFORCE=0
: Element forces are not written to file 
LFORCE=1
: Element forces are written to file


LCURV: integer, default: 0
: Control parameter
LCURV=0
: Element curvatures are not written to file 
LCURV=1
: Element curvatures are written to file


IASCII: integer, default: 0
: Control parameter
IASCII=0
: Unformatted snapshot files 
IASCII=1
: Formatted snapshot files


XCGVES: real
: Static X coordinate of vessel CG 
YCGVES: real
: Static Y coordinate of vessel CG 
ZCGVES: real
: Static Z coordinate of vessel CG
Line specification, NLIC input lines
LINEID

LINEID: character(8)
: Line identifier to be written to file. You may specifyALL
to include all the lines in the system
The lines are written only if at least the displacements of the end nodes are stored, see data group File storage of displacement response for storage information.
Vessel description, NPVESD input lines. The specified points are connected by one line to illustrate a part of the vessel contour
IPV XVT YVT ZVT

IPV: integer
: Coordinate number 
XVT: real
: Vessel’s Xcoordinate in global system, relative to the vessel reference point \(\mathrm {[L]}\) 
YVT: real
: Vessel’s Ycoordinate \(\mathrm {[L]}\) 
ZVT: real
: Vessel’s Zcoordinate \(\mathrm {[L]}\)
The vessel points are in global system, but they are relative to the vessel reference point (the attachment point).
4.2.17. Stress time series calculated from the time domain analysis
This option allows for calculation of stresses in circular metallic homogeneous risers.
The stress time series are calculated based on the stored force time
series from DYNMOD
and the component properties specified in INPMOD
.
Stresses may only be calculated for CRS1
and CRS0
components.
Stress time series are calculated for specified points on the tube circumference.
Output options, one input line
IOP IDOF IT1 NTS ISUBST NNELC

IOP: integer
: Code for type of output
IOP = 1
: Time series 
IOP = 2
: Time series statistics 
IOP = 3
: Spectral analysis 
IOP = 1
in present version


IDOF: integer
: Stress components type 1
IDOF = 1/11
: Axial stress at end 1/2 
IDOF = 2/12
: Torsional stress at end 1/2 
IDOF = 3/13
: Bending stress at end 1/2 
IDOF = 4/14
: Axial + bending stress at end 1/2 
IDOF = 5/15
: Shear stress at end 1/2 
IDOF = 6/16
: Shear stress + torsional stress at end 1/2 
IDOF = 7/17
: Equivalent stress at end 1/2 
IDOF = 8/18
: Hoop stress at end 1/2 
IDOF = 9/19
: Radial stress at end 1/2 
IDOF = 21/22
: External pressure at end 1/2 
IDOF = 23/24
: Internal pressure at end 1/2


IT1: integer
: First stored time step to be included 
NTS: integer
: Number of stored time steps to be included (fromIT1
).
A large number includes the remaining time steps


ISUBST: integer, default: 0
: Code for subtracting the static stress contributions
ISUBST = 0
: Total stresses calculated 
ISUBST = 1
: Static stress is subtracted


NNELC: integer
: Number of lines used for element specification
Point for stress calculation, one input line
THETA INEX IOPPRE

THETA: real, default: 0.0
: Angle from local yaxis for stress calculation \(\mathrm {[Deg]}\) 
INEX: integer, default: 2
: Stress location switch
INEX = 1
: Inner wall 
INEX = 2
: Outer wall


IOPPRE: integer, default: 1
: Code for updating inner and outer pressure values.
IOPPRE = 1
: Static inner and outer pressure used.
Outer pressure is calculated as hydrostatic pressure from MWL.


IOPPRE = 2
: Updated inner and outer pressure used.
Outer pressure is calculated as hydrostatic pressure from MWL.


IOPPRE < 0
: Wall forces calculated using outer area given by the pipe diameter or the alternative cross section diameter.
Corresponds to evenly distributed shear forces between buoyancy material and pipe.

Warning: This option is under development!


Nonlinear time domain analysis only.
In the present version, the external pressure is calculated as a hydrostatic pressure from the MWL. The external pressure is updated for all structural elements.
The internal pressure is updated for all elements that are part of a Main Riser Line.
Element specification, NNELC input lines
LINEID ISEG IELM

LINEID: character(8)
: Line identifier.
You may specify
ALL
to include all lines


ISEG: integer/character
: Segment number.
You may specify
ALL
to include all segments. 
ENDS
includes the end segments on the line


IELM: integer/character
: Element number.
ALL
includes all elements, and 
ENDS
includes end elements on the above specified segment

Stresses may only be calculated for elements for which forces are stored, see data group File storage for internal forces for storage information. If the user specifies elements for which forces are not stored, these elements are ignored.
The data group Storage information may be used to obtain an overview of the stored data.
Options for the output distribution functions of the stress time series statistics, one input line
This input line is given only if IOP=2
.
NCL XCMIN XCMAX

NCL: integer
: No of classes in the output distribution functions (i.e. no of points on the abscissa axis)
0<NCL<41


XCMIN: real
: Range of argument values for output distribution functions isXCMIN*sx(1)  XCMAX*sx(1)
in whichsx(1)
is the standard deviation ofx
estimated from the first sequence. 
XCMAX: real
:
4.2.18. Stress envelope curves
This option allows for calculation of stress envelopes from the element
forces stored in DYNMOD
, see data group File storage for internal forces for
storage information.
The data group Storage information may be used to obtain an overview of the stored data.
Print options, one input line
LINEID IDOF1 IDOF2 IDOF3

LINEID: character(8)
: Line identifier
ILINE = ALL
: Stresses in all lines calculated


IDOF1:integer
: Stress component type 1
IDOF1 = 1
: Axial stress 
IDOF1 = 2
: Torsional stress 
IDOF1 = 3
: Bending stress 
IDOF1 = 4
: Axial + bending stress 
IDOF1 = 5
: Shear stress 
IDOF1 = 6
: Shear + torsional stress 
IDOF1 = 7
: Equivalent stress 
IDOF1 = 8
: Hoop stress 
IDOF1 = 9
: Radial stress


IDOF2:integer
: Stress component type 2
See
IDOF1

Dummy at present


IDOF3:integer
: Stress component type 3
See
IDOF1

Dummy at present

Stress calculations options, one input line
TSTA TEND IOP DUR

TSTA: real, default: 0
: Start time in stress time series \(\mathrm {[T]}\) 
TEND: real, default: 0
: End time in stress time series \(\mathrm {[T]}\)
TEND = 0.0
: Until last time step used


IOP: integer, default: 0
: Code for envelope type
IOP = 1
: Min and max values presented 
IOP = 2
: Maximum range 
IOP = 3
: Standard deviations 
IOP = 4
: Estimated extreme values (not yet implemented)


DUR: real, default: 10800
: Duration used in extreme value estimation \(\mathrm {[T]}\)
Dummy parameter in present version

Stress calculation location, one input line
NPCS IOPPR THETA INEX IOPPRE

NPCS: integer, default: See below
: Number of points around the crosssection
= 0
: max stresses estimated


IOPPR: integer, default: 0
: Print option
IOPPR = 0
: Print maximum stresses only 
IOPPR > 0
: Print stresses at allNPRCS
points


THETA: real, default: 0
: Angle for stress calculation \(\mathrm {[Deg]}\)
Dummy for
NPCS>1


INEX: integer, default: 2
: Stress loction switch
INEX = 1
: Inner wall 
INEX = 2
: Outer wall


IOPPRE: integer, default: 1
: Code for updating inner and outer pressure values.
IOPPRE = 1
: Static inner and outer pressure used. 
IOPPRE = 2
: Updated inner and outer pressure used. 
Outer pressure calculated as hydrostatic pressure from MWL.

IOPPRE < 0
: Wall forces calculated using outer area given by the pipe diameter or the alternative cross section diameter.
Corresponds to evenly distributed shear forces between buoyancy material and pipe.

Warning: This option is under development!


Nonlinear time domain analysis only.
The default value for NPCS
is dependent on the value specified above
for IOP
: Default is 0
for IOP = 1
, otherwise it is 4
.
In the present version, the external pressure is calculated as a hydrostatic pressure from the MWL. The external pressure is updated for all structural elements.
The internal pressure is updated for all elements that are part of a Main Riser Line.
Stress calculation parameters, one input line
IOPSTR ASTI WSTI DIASTI THSTI EMOD

IOPSTR: integer, default: 0
: Option for stress calculation
IOPSTR=0
: Stresses calculated from bending moment (recommended) 
IOPSTR=1
: Stresses calculated from curvatures


ASTI: real, default: 0
: Alternative cross sectional area \(\mathrm {[L^2]}\) 
WSTI: real, default: 0
: Alternative cross section modulus \(\mathrm {[L^3]}\) 
DIASTI: real, default: 0
: Alternative cross section diameter \(\mathrm {[L]}\) 
THSTI: real, default: 0
: Alternative cross section wall thickness \(\mathrm {[L]}\) 
EMOD: real, default: 0
: Modulus of elasticity \(\mathrm {[F/L^2]}\)
Bending stresses are calculated from curvature, diameter and
EMOD
ifIOPSTR=1
andEMOD>0

\(\mathrm {WST=\frac{2}{EMOD\times DIAST}}\)

The default values of 0 for ASTI, WSTI, DIASTI, THSTI
and EMOD
are
interpreted as no change from the crosssectional properties given in
INPMOD
.
4.2.19. Riser stroke time series from time domain analysis
The riser stroke is calculated for the supernode specified in DYNMOD
from the motions of the vessel and the vertical displacement of
specified supernode.
This option is not of interest if the terminal point of the riser is vertically fixed to the vessel.
Option to calculate the riser stroke time series, one input line
IOP IMOT IDERIV IT1 NTS

IOP: integer
: Code for type of output
IOP = 1
: Time series 
IOP = 2
: Time series statistics 
IOP = 3
: Spectral analysis


IMOT: integer
:
IMOT = 1
: Stroke 
IMOT = 2
: Platform heave motion only 
IMOT = 3
: Risers upper end heave motion only


IDERIV: integer
:
IDERIV = 0
: Original 
IDERIV = 1
: First derivative 
IDERIV = 2
: Second derivative


IT1: integer
: First stored time steps to be included 
NTS: integer
: Number of stored time steps to be included
4.2.20. Code check curves
This option allows for code check of the response.
Main output options, one input line
LINEID IOPCOD IOP IDIST DUR PROB

LINEID: character(8)
: Line identifier
LINEID = ALL
: All lines checked


IOPCOD: integer, default: 1
: Option for type of code check
IOPCOD = 1
: titanium code check


IOP: integer, default: 2
: Option for using maximum or estimated extreme values
IOP = 1
: Maximum values from stress time series used 
IOP = 2
: Estimated extreme values used


IDIST: integer, default: 2
: Distribution type used in extreme value estimation
IDIST = 1
: Rayleigh distribution used 
IDIST = 2
: Three parameter Weibull used 
Dummy for
IOP = 1


DUR: real, default: 10800
: Duration used in extreme value estimation \(\mathrm {[T]}\)
Dummy for
IOP = 1


PROB: real, default: 0
: Probability level used in extreme value estimation
PROB = 0.0
: Expected maximum value used 
Dummy for
IOP = 1

Time range and crosssection points, one input line
TSTA TEND NPCS IOPPR

TSTA: real, default: 0
: Start time in stress time series \(\mathrm {[T]}\) 
TEND: real, default: 0
: End time in stress time series \(\mathrm {[T]}\)
TEND = 0.0
: Until last time step used


NPCS: integer >= 0, default: see below
: Number of points around the crosssection 
IOPPR: integer, default: 0
: Print option
The default value for NPCS
is dependent on the value specified above
for IOP
:
Default is 0
for IOP = 1
, otherwise it is 4
.
Static load step and load factors, one input line
ISTEPF GAMF GAMC GAME GAMR

ISTEPF: integer, default: 0
: Static step number for functional loads
ISTEPF = 0
: Final static load step is used


GAMF: real, default: 1
: Load factor for functional loads 
GAMC: real, default: 1
: Load effect factor for condition 
GAME: real, default: 1
: Load factor for environmental loads 
GAMR: real, default: 1
: Resistance factor
Stress calculation parameters, one input line
SMYS EMOD NU F0 SMYSB TADD

SMYS: real > 0
: Specified minimum yield stress \(\mathrm {[F/L^2]}\) 
EMOD: real > 0
: Modulus of elasticity \(\mathrm {[F/L^2]}\) 
NU: real, default: 0.3
: Poisson’s ratio 
F0: real, default: 0.005
: Initial ovality
\(\mathrm {=(D_{max}D_{min})/D}\)


SMYSB: real, default: SMYS
: Specified minimum stress used in axial capacity \(\mathrm {[F/L^2]}\) 
TADD: real, default: 0
: Additional torsion moment \(\mathrm {[FL]}\)
Typical values of SMYS
and EMOD
for steel are in the order of
\(\mathrm {[SMYS=220.0E3kN/m^2]}and\) if the units m
and
kN
were chosen in INPMOD
.
Crosssection parameters, one input line
ASTI WSTI DIASTI THSTI

ASTI: real, default: 0
: Alternative cross sectional area \(\mathrm {[L^2]}\) 
WSTI: real, default: 0
: Alternative cross section modulus \(\mathrm {[L^3]}\) 
DIASTI: real, default: 0
: Alternative cross section diameter \(\mathrm {[L]}\) 
THSSTI: real, default: 0
: Alternative cross section wall thickness \(\mathrm {[L]}\)
The default values of 0 are interpreted as no change from the
crosssectional properties given in INPMOD
4.2.21. Time domain fatigue damage
This option allows for calculation of fatigue damage calculation from axial and bending stresses in circular metallic homogeneous risers using a specified SN curve and rainflow cycle counting.
The calculated fatigue damage is per year of the specified environmental conditions.
The fatigue damage is calculated based on the stored force time series
from DYNMOD
(see data group File storage for internal forces for storage
information) and the component properties specified in INPMOD
.
Stresses may only be calculated for CRS1
and CRS0
components.
The fatigue damage is calculated for a specified number of points on the tube circumference.
Control data, one input line
NSECT NPCS IOPPR TBEG TEND IOPSTR FATID

NSECT: integer
: Number of riser cross sections to be considered
NSECT = 0
: All cross section where forces are available is included in the analysis


NPCS: integer
: Number of points in the cross section where fatigue is calculated 
IOPPR: integer
: Print option for fatigue results
IOPPR = 0
: Print results only for most critical point in cross section 
IOPPR > 0
: Print results for allNPCS
points


TBEG: real
: Beginning of stored stress time series for fatigue calculation Number \(\mathrm {[T]}\) 
TEND: real
: End of stored stress time series for fatigue calculation \(\mathrm {[T]}\)
Default is the last stored time step


IOPSTR: integer, default: 0
: Option for stress calculation
IOPSTR=0
: Bending stresses calculated from bending moment (recommended) 
IOPSTR=1
: Bending stresses calculated from curvatures. EMOD and DIAST must be given


FATID: character(16)
: Identifier for fatigue calculation. Used in result presentation only
The remaining of the time series is used if TEND
is less or equal to
TBEG
(Default is full time series).
Crosssectional data, one input line
DSCFA DSCFY DSCFZ ASI WSTI DIAST EMOD CFRS LFRS TEFF

DSCFA: real, default: 1
: Default stress concentration factor for axial force contribution 
DSCFY: real, default: DSCFA
: Default stress concentration factor for bending about the local Y axis 
DSCFZ: real, default: DSCFA
: Default stress concentration factor for bending about the local Z axis 
ASI: real, default: See below
: Optional crosssectional area \(\mathrm {[L^{2}]}\) 
WSTI: real, default: See below
: Optional section modulus \(\mathrm {[L^{3}]}\). Dummy if stresses are are calculated from curvature (IOPSTR = 1
) 
DIAST: real, default: See below
: Cross section diameter. Used to calculate tresses from curvature ifIOPSTR = 1
. Otherwise not used. 
EMOD: real, default: See below
: Modules of elasticity. Used to calculate tresses from curvature ifIOPSTR = 1
. Otherwise not used. 
CFRS: real, default: 0
: Constant correction coefficient, friction stress \(\mathrm {[FL^{2}]}\) 
LFRS: real, default: 0
: Linear correction coefficient \(\mathrm {[L^{2}]}\) 
TEFF: real, default: See below
: Effective thickness used together with the reference thicknessTREF
given below in thickness correction \(\mathrm {[L]}\)
The crosssectional area, modulus and thickness defined for each cross
section in INPMOD
are used as defaults for ASI
, WSTI
and TEFF
.
Stress range correction due to friction is given as:
\(\mathrm {\Delta \sigma _f=CFRS+LFRS\times T_{avg}}\)
\(\mathrm {T_{avg}}\) is the static value of the tension. The friction stress correction is added after Rainflow counting of the stress time series due to axial force and bending
\(\mathrm {\sigma _{tot}(t)=\sigma _{axial}(t)\times SCFA+\sigma _{Ybending}(t)\times SCFY+\sigma _{Zbending}(t)\times SCFY}\)
The units of the friction correction coefficients must be consistent
with the Selection of unit system physical constant
in INPMOD
.
SN curve data, two input lines
Fatigue capacity curve description
NOSL LIMIND FATLIM RFACT TREF KEXP

NOSL: integer ⇐ 5, default: 1
: Number of straight lines defining the SN curve 
LIMIND: integer, default: 0
: Fatigue limit indicator
LIMIND < 0
: Fatigue limit in terms of stress cycles is specified 
LIMIND = 0
: No fatigue limit for present curve 
LIMIND > 0
: Fatigue limit in terms of stress range is specified


FATLIM: real, default: 0
: Fatigue limit, interpretation dependent onLIMIND
. See Example 1: 2 segments and fatigue limit.
LIMIND < 0
: Base 10 logarithm of number of stress cycles for which the SN curve becomes horizontal 
LIMIND = 0
:FATLIM
is dummy 
LIMIND > 0
: Stress range level for which the SN curve becomes horizontal \(\mathrm {[S]}\) 
See
RFACT
below


RFACT: real, default: 1
: Factor between the stress unit \(\mathrm {[S]}\) used to define the SN curve and the force and length units \(\mathrm {[F]}\) and \(\mathrm {[L]}\) chosen inINPMOD

\(\mathrm {S\times RFACT=\frac{F}{L^2}}\)


TREF: real, default: 0
: Reference thickness for thickness correction \(\mathrm {[L]}\). IfTREF = 0
the thickness correction will be omitted. 
KEXP: real, default: 0
: Exponent for thickness correction
If \(\mathrm {kN}\) and \(\mathrm {m}\) were chosen
as force and length units while the SN curve is given in
\(\mathrm {MPa}\), RFACT
should be set to 0.001
.
If the SI units \(\mathrm {N}\) and \(\mathrm {m}\)
were chosen for force and length and the SN curve is in
\(\mathrm {MPa}\), RFACT
should be set to 1.0E6
.
Fatigue capacity curve constants
RM1 RC1 RMi RNCi ...

RM1: real
: Slope of the SN curve. First curve segment forNOSL>1
, total curve forNOSL=1
. (log cycles / log stress) 
RC1: real
: Constant defining the SN curve. First segment or total curve 
RMi: real
: Slope of curve segment i, i=2, …,NOSL

RNCi: real
: Transition point between curve segment (i1), and i, i=2,…,NOSL
 (log cycles)
See Frequency domain fatigue damage for details of the fatigue curve specification.
For a single slope SN curve, log cycles as a function of log stress:
\(\mathrm {logN=RC1RM1\times log\Delta S}\)
Where:

\(\mathrm {N}\): Number of cycles to failure

\(\mathrm {\Delta S}\): Stress range
or log stress as a function of log cycles:
\(\mathrm {log\Delta S=\frac{RC1}{RM1}\frac{logN}{RM1}}\)
Cross section specification, NSECT input lines
LINEID ISEG IEL IEND SCFA SCFY SCFZ

LINEID: character(8)
: Line identifier 
ISEG: integer >= 0
: Segment number on line
= 0
: All segments in specified line


IEL: integer >= 0
: Local element number on specified segment
= 0
: All elements in specified segment


IEND: integer
:
IEND = 0
: Cross sections at both ends checked 
IEND = 1
: Cross section at end with smallest node number checked 
IEND = 2
: Cross section at end with largest node number checked


SCFA: real, default: DSCFA
: Stress concentration factor for axial force contribution 
SCFY: real, default: SCFA
: Stress concentration factor for bending about local Y axis 
SCFZ: real, default: SCFA
: Stress concentration factor for bending about local Z axis
Time domain forces for the specified elements must be stored in DYNMOD
, see data group File storage for internal forces for storage information.
The data group Storage information may be used to obtain an overview of the stored data.
If several specifications match an element, the first specification will be used.
4.2.22. Time domain longterm data
This option allows for calculation of transfer function modulus and, in
the future, also distribution parameters for the stresses from axial and
bending force in circular metallic homogeneous risers. The results are
intended to be processes in a longterm analysis like in LONFLX
and
LOSSTA
.
The results are calculated based on the stored force time series from
DYNMOD
(see data group File storage for internal forces for storage information)
and the component properties specified in INPMOD
. Stresses may only be
calculated for CRS1
and CRS0
components.
The transfer functions are calculated for a specified number of points on the tube circumference.
Control data, one input line
NSECT NPCS TBEG TEND

NSECT: integer, default: 0
: Number of riser cross sections to be considered
NSECT = 0
: All cross section where forces are available is included in the analysis


NPCS: integer, default: 16
: Number of points in the cross section where fatigue is calculated 
TBEG: real, default: 0
: Beginning of stored stress time series for fatigue calculation Number \(\mathrm {[T]}\) 
TEND: real, default: 0
: End of stored stress time series for fatigue calculation \(\mathrm {[T]}\)
Default is the last stored time step.

The remaining of the time series is used if TEND
is less or equal to
TBEG
(Default is full time series).
Calculation control data, one input line
MXFRQ FLOW FHIG IDIST

MXFRQ: integer
: Maximum number frequencies in the output of transfer functions 
FLOW: real
: Lower frequency limit in the printing 
FHIG: real
: Upper frequency limit in the printing 
IDIST: integer
: Distribution type (Future use)
The actual number of frequencies in the output will usually be somewhat
less than MXFRQ
because the printing is going in integer steps over
the calculated Fourier components. The intermediate points is used for
smoothing of the transfer function
Cross sectional data, one input line
DSCFA DSCFY DSCFZ ASI WSTI

DSCFA: real, default: 1
: Default stress concentration factor for axial force contribution 
DSCFY: real, default: DSCFA
: Default stress concentration factor for bending about Y axis 
DSCFZ: real, default: DSCFA
: Default stress concentration factor for bending about Z axis 
ASI: real, default: 0
: Optional cross sectional area 
WSTI: real, default: 0
: Optional section modulus
The cross sectional area and modulus defined in INPMOD
is used by
default.
Cross section specification NSECT input lines
ILIN ISEG IEL IEND

ILIN: integer
: Line number 
ISEG: integer
: Segment number on line 
IEL: integer
: Local element number on specified segment 
IEND: integer
:
IEND = 1
: Cross section at end with smallest node number checked 
IEND = 2
: Cross section at end with largest node number checked

Time domain forces for the specified elements must be stored in DYNMOD
, see data group File storage for internal forces for storage information.
The data group Storage information may be used to obtain an overview of the stored data.
5. Data Group D: Output from FREMOD
5.1. Frequency domain layer damage
This data group may be used to calculate wear and fatigue of tendons in a nonbonded flexible pipe cross section.
5.1.2. Control data, one input line
NLAYER NSECT

NLAYER: integer
: Number of layers to be considered 
NSECT: integer
: Number of riser cross sections to be considered
5.1.3. Layer data, 2 × NLAYER input lines
Axial stress and friction per unit (pressure/axial force/curvature)
IDLAY ALFA1 ALFA2 ALFA3 ALFA4 ALFA5

IDLAY: integer
: Unique identification number for the layer data 
ALFA1: real
: Axial stress in helix per unit pressure (difference) 
ALFA2: real
: Axial stress in helix per unit axial force 
ALFA3: real
: Axial stress in helix per unit pipe curvature \(\mathrm {[1/bendingradius]}\) 
ALFA4: real
: Friction stress per unit pressure (difference) 
ALFA5: real
: Friction stress per unit axial force
5.1.4. Wear and geometrical data; 1 data string:
BETA1 BETA2 THICK WSAFE SIGUL SIGLI

BETA1: real
: Wear factor per unit curvature and pressure 
BETA2: real
: Wear factor per unit curvature and axial force 
THICK: real
: Thickness of layer 
WSAFE: real ⇐ 1
: Safety factor for wear 
SIGUL: real
: Ultimate stress 
SIGLI: real
: Limit stress
5.1.5. Cross section specification, NSECT input lines
LINEID ISEG IEL IEND IDLAY1 ... IDLAYn

LINEID: character(8)
: Line identifier (dummy forIEL < 0
) 
ISEG: integer
: Segment number (dummy forIEL < 0
) 
IEL: integer
: Element number
IEL > 0
: local element number 
IEL < 0
: global element number


IEND: integer
:
IEND = 1
: Cross section at end with smallest node number checked 
IEND = 2
: Cross section at end with largest node number checked


IDLAY1: integer
: First layer to be checked 
IDLAYn: integer
: Last layer to be checked
Frequency domain results for the specified element/ends must be stored
on the FREMOD
result file ifnfre
.
5.2. Frequency domain fatigue damage
5.2.2. Control data, one input line
NOFC NSECT IRES

NOFC: integer
: Number of SN curves 
NSECT: integer
: Number of riser cross sections to be considered 
IRES: integer
: Response print option
IRES>0
: print of total fatigue damage only 
IRES<0
: print of fatigue contributions

5.2.3. SN data, 2 × NOFC input lines
Fatigue capacity curve description
ISNC NOSL LIMIND FATLIM RFACT

ISNC: integer
: SN curve number  must be given in ascending order 
NOSL: integer, default: 1
: Number of straight lines defining the SN curve 
LIMIND: integer, default: 0
: Fatigue limit indicator
LIMIND < 0
: Fatigue limit in terms of stress cycles is specified 
LIMIND = 0
: No fatigue limit for present curve 
LIMIND > 0
: Fatigue limit in terms of stress range is specified


FATLIM: real, default: 0
: Fatigue limit, interpretation dependent on LIMIND. See Example 1: 2 segments and fatigue limit.
LIMIND < 0
: Base 10 logarithm of number of stress cycles for which the SN curve becomes horizontal 
LIMIND = 0
:FATLIM
is dummy 
LIMIND > 0
: Stress range level for which the SN curve becomes horizontal


RFACT: real, default: 1
: Factor between the stress unit \(\mathrm {[S]}\) used to define the SN curve and the force and length units \(\mathrm {[F]}\) and \(\mathrm {[L]}\) chosen inINPMOD

\(\mathrm {S\times RFACT=\frac{F}{L^2}}\)

If \(\mathrm {kN}\) and \(\mathrm {m}\) were chosen
as force and length units while the SN curve is given in
\(\mathrm {MPa}\), RFACT
should be set to 0.001
.
If the SI units \(\mathrm {N}\) and \(\mathrm {m}\)
were chosen for force and length and the SN curve is in
\(\mathrm {MPa}\), RFACT
should be set to 1.0E6
.
Fatigue capacity curve constants
RM1 RC1 RMi RNCi

RM1: real
: Slope of the SN curve. First curve segment forNOSL>1
, total curve forNOSL=1

RC1: real
: Constant defining the SN curve. First segment or total curve 
RMi: real
: Slope of curve segmenti, i=2, …, NOSL

RNCi: real
: Transition point between curve segment(i1)
, andi
,i=2, …, NOSL
(log cycles)
Explanation of the input parameters in input lines SN data, 2 × NOFC input lines
(above). All parameters are found in figure below.
Example 1: 2 segments and fatigue limit. Note that FATLIM can be alternatively specified
Example 2: 3 segments and not fatigue limit. Illustration of input data for fatigue capacity curve definition
The SN curves defined by input parameters are always assumed to relate \(\mathrm {\Delta S}\) (stress range) to number of cycles before failure.
A straightlined SN curve in loglog scale is in general defined as
\(\mathrm {N=C\times \Delta S^m}\)
or
\(\mathrm {logN=logC+m\times log\Delta S}\)
Where:

\(\mathrm {N}\): Number of cycles to failure

\(\mathrm {\Delta S}\): Stress range
The two input parameters used to define the SN curves are directly found in the equation above, namely

\(\mathrm {RC=logC\quad }\) (always positive)

\(\mathrm {RM=m\quad }\) (always negative)
If the user has an SN curve without having these parameters explicitly defined, they can be calculated as follows:
Using the two points A and B on the figure to define the straight line, we have
\(\mathrm {logN=\frac{logN_2logN_1}{log\Delta S_2log\Delta S_1}\times log\Delta S\frac{logN_2logN_1}{log\Delta S_2log\Delta S_1}\times log\Delta S_1+logN_1}\)
Hence:

\(\mathrm {RM= \frac{logN_2logN_1}{log\Delta S_2log\Delta S_1}}\) (always negative)

\(\mathrm {RM\times log\Delta S_1+logN_1}\) (always positive)
The relation between these parameters specified for different unit systems is easily found from the equations above.
5.2.4. Cross section specification, NSECT input lines
LINEID ISEG IEL IEND SCF IFAT1 ... IFATn

LINEID: character(8)
: Line identifier 
ISEG: integer
: Segment number in line 
IEL: integer
: Local element number in specified segment 
IEND: integer
:
IEND = 1
: Cross section at end with smallest node number checked 
IEND = 2
: Cross section at end with largest node number checked


SCF: real, default: 1
: Stress concentration factor 
IFAT1: integer
: First SN curve to be checked 
.

.

.

IFATn: integer
: Last SN curve to be checked
Frequency domain results for the specified element/ends must be stored
in FREMOD
5.3. Frequency domain force results
5.3.2. Specification of number of sections, one input line
NSECT

NSECT: integer
: Number of sections to be specified
5.3.3. Section specification, NSECT input lines
LINEID ISEG IELM IEND

LINEID: character(8)
: Line identifier 
ISEG: integer
: Segment number 
IELM: integer
: Element number 
IEND: integer
: Element end (1 or 2)
Both local (LINEID SEG ELML
) numbering and global element numbering
(ELMG
) can be given
This option is valid for linear bending stiffness only, i.e. IEJ=1
in
INPMOD
.
Example:
parameter/numbering  lineid  iseg  ielm  iend 

local: 
1 
2 
4 
1 
Eqv.global: 
0 
0 
25 
1 
6. Description of STARTIMES File
6.1. Description of STARTIMES file generated by OUTMOD
6.1.1. General comments
A time series generated by OUTMOD
is identified by a time series
number and a version number. Response type is identified by the time
series number, see description below. The selected node number or
element number is identified by the version number.
For one response type versions are numbered 1,2,3,4……N
according
to:

version
1
for 1st selected element/node, 
version
2
for 2nd selected element/node, 
…

version
N
for last selected element/node
Time series number and version number are printed to OUTMOD
result
file for each selected response. In addition, FEM
element/node number
is included in identification text for each time series stored on
STARTIMES
file.
6.1.2. Example
Output of time series of dynamic axial force for element 1, 33, 45
is
specified in OUTMOD
. Identifiers to generated time series are:

Element 1: 40.01 (time series number 40, version number 1)

Element 33: 40.02 (time series number 40, version number 2)

Element 45: 40.03 (time series number 40, version number 3)
6.2. Description of time series numbers
Time series number  Contents 

WF motion time series 

1 
HF surge 
2 
HF sway 
3 
HF heave 
4 
HF roll 
5 
HF pitch 
6 
HF yaw 
LF motion time series 

7 
LF surge 
8 
LF sway 
9 
LF yaw 
WF and LF motion time series 

10 
HF + LF surge 
11 
HF + LF sway 
12 
HF + LF yaw 
Wave elevation time series 

13 
Wave elevation 
Wave kinematics 
(Not implemented in 
14 
Water particle velocity: xdirection 
15 
Water particle velocity: ydirection 
16 
Water particle velocity: zdirection 
17 
Water particle acceleration: xdirection 
18 
Water particle acceleration: ydirection 
19 
Water particle acceleration: zdirection 
Dyndisp time series 

20 
Dynamic displacement: global xdirection 
21 
Dynamic displacement: global ydirection 
22 
Dynamic displacement: global zdirection 
23 
Dynamic rotation about: xdirection 
24 
Dynamic rotation about: ydirection 
25 
Dynamic rotation about: zdirection 
Calcurv time series 

26 
Total curvature calculated from nodal coordinates 
Element angle time series 

27 
Element angle \([\mathrm {deg}]\) 
IRELCO = 1: Angle between global zaxis and one element 

IRELCO = 2: Angle between support vessel axis and one element 

IRELCO = 3: Angle between two elements 

Distance time series 

28 
Distance 
IDITYP = 1: Distance between specified segments 

IDITYP = 2: Distance between specified segments on a line and a globally fixed line 

IDITYP = 3: Distance between specified segments on a line and a line fixed on the vessel 

Stroke time series 

29 
Stroke 
IMOT = 1: Stroke 

IMOT = 2: Platform heave motion only 

IMOT = 3: Riser heave motion only 

Support force time series 

30 
Support force component: global xdirection 
31 
Support force component: global ydirection 
32 
Support force component: global zdirection 
Dynforce time series 

40 
Dynamic axial force 
41 
Dynamic Torsional moment 
42 
Dynamic bending moment about local yaxis: End 1 
43 
Dynamic bending moment about local yaxis: End 2 
44 
Dynamic bending moment about local zaxis: End 1 
45 
Dynamic bending moment about local zaxis: End 2 
46 
Dynamic shear force in local ydirection: End 1 
47 
Dynamic shear force in local ydirection: End 2 
48 
Dynamic shear force in local zdirection: End 1 
49 
Dynamic shear force in local zdirection: End 2 
Dyncurv time series 

50 
Dynamic curvature about local yaxis: End 1 
51 
Dynamic curvature about local yaxis: End 2 
51 
Dynamic curvature about local zaxis: End 1 
53 
Dynamic curvature about local zaxis: End 2 
Totforce time series 

54 
Axial force 
55 
Torsional moment 
56 
Bending moment about local yaxis: End 1 
57 
Bending moment about local yaxis: End 2 
58 
Bending moment about local zaxis: End 1 
59 
Bending moment about local zaxis: End 2 
60 
Shear force in local ydirection: End 1 
61 
Shear force in local ydirection: End 2 
62 
Shear force in local zdirection: End 1 
63 
Shear force in local zdirection: End 2 
64 
Axial wall force 
Totdisp time series 

65 
Total displacements in global xdirection 
66 
Total displacements in global ydirection 
67 
Total displacements in global zdirection 
Stress time series 

75 
Axial + Bending stress: End 1 
76 
Axial + Bending stress: End 2 
77 
Torsional stress 
78 
Equivalent stress: End 1 
79 
Equivalent stress: End 2 