AQWA Reference Manual
AQWA Reference Manual
AQWA Reference Manual
ANSYS, Inc.
Southpointe
275 Technology Drive
Canonsburg, PA 15317
ansysinfo@ansys.com
http://www.ansys.com
(T) 724-746-3304
(F) 724-514-9494
Release 15.0
November 2013
ANSYS, Inc. is
certified to ISO
9001:2008.
Disclaimer Notice
THIS ANSYS SOFTWARE PRODUCT AND PROGRAM DOCUMENTATION INCLUDE TRADE SECRETS AND ARE CONFIDENTIAL AND PROPRIETARY PRODUCTS OF ANSYS, INC., ITS SUBSIDIARIES, OR LICENSORS. The software products
and documentation are furnished by ANSYS, Inc., its subsidiaries, or affiliates under a software license agreement
that contains provisions concerning non-disclosure, copying, length and nature of use, compliance with exporting
laws, warranties, disclaimers, limitations of liability, and remedies, and other provisions. The software products
and documentation may be used, disclosed, transferred, or copied only in accordance with the terms and conditions
of that software license agreement.
ANSYS, Inc. is certified to ISO 9001:2008.
Third-Party Software
See the legal information in the product help files for the complete Legal Notice for ANSYS proprietary software
and third-party software. If you are unable to access the Legal Notice, please contact ANSYS, Inc.
Published in the U.S.A.
Table of Contents
1. Introduction ............................................................................................................................................ 1
1.1. Conventions ..................................................................................................................................... 1
1.1.1. Axes Systems ........................................................................................................................... 1
1.1.2. Wave/Wind/Current Direction .................................................................................................. 2
1.1.3. Phase Angle ............................................................................................................................. 3
1.2. The Structure Definition and Analysis Position ................................................................................... 3
1.3. Units ................................................................................................................................................ 4
1.4. Aqwa Files ........................................................................................................................................ 5
2. Analysis Stages ....................................................................................................................................... 9
2.1. General Description .......................................................................................................................... 9
2.2. The Restart Stages .......................................................................................................................... 10
2.3. Restart Stage 1 ............................................................................................................................... 11
2.4. Restart Stage 2 ............................................................................................................................... 14
2.5. Restart Stage 3 ............................................................................................................................... 18
2.6. Restart Stage 4 ............................................................................................................................... 18
2.7. Restart Stage 5 ............................................................................................................................... 22
2.8. Restart Stage 6 ............................................................................................................................... 23
3. Data Preparation for the Aqwa Suite .................................................................................................... 25
3.1. Features Common to All Data Categories ......................................................................................... 26
3.1.1. Optional User Identifier .......................................................................................................... 27
3.1.2. Compulsory Data Category Keyword ...................................................................................... 27
3.1.3. Compulsory 'END' Statement .................................................................................................. 27
3.1.4. Format Requirements ............................................................................................................. 27
3.2. Classification of Data Categories 1 TO 18 ......................................................................................... 28
3.2.1. Stopping and Starting the Program During an Analysis .......................................................... 28
3.2.2. Data Categories 1 to 5 - Geometric Definition and Static Environment ..................................... 28
3.2.3. Data Categories 6 to 8 - The Radiation/Diffraction Analysis Parameters ................................... 29
3.2.4. Data Categories 9 to 18 - Definition of Analysis-Dependent Parameters ................................... 29
3.3. Default Values Assumed by the Program ......................................................................................... 29
3.3.1. Omission of Data Categories ................................................................................................... 30
3.3.2. Omission of Data Records within a Data Category ................................................................... 30
3.3.3. Omission of Fields on a Data Record ....................................................................................... 30
3.3.4. Default Values ........................................................................................................................ 30
3.4. The Data Category Series for One or More Structures ....................................................................... 31
3.4.1. Data Category Series - Definition ............................................................................................ 31
3.4.2. The Data Category Series Terminator - FINI .............................................................................. 31
3.4.3. Omission of Data Categories within a Data Category Series ..................................................... 31
4. The Preliminary Data (Data Category 0) ............................................................................................... 33
4.1. The JOB Data Record ....................................................................................................................... 33
4.2. The TITLE Data Record ..................................................................................................................... 34
4.3. The OPTIONS Data Record ............................................................................................................... 34
4.4. The RESTART Data Record ............................................................................................................... 35
5. Node Number and Coordinates (Data Category 1) ............................................................................... 37
5.1. General Description ........................................................................................................................ 37
5.2. Data Category Header ..................................................................................................................... 37
5.3. The Coordinate Data Record ............................................................................................................ 38
5.4.The Coordinate Data Record with Rotational Node Generation ......................................................... 38
5.5. The STRC Data Record - Coordinate Structure Association ................................................................ 40
5.6. The Coordinate Data Record with Offset .......................................................................................... 41
5.7. The Coordinate Data Record with Translation ................................................................................... 41
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
iii
Reference Manual
5.8. The Coordinate Data Record with Mirror Node Generation ............................................................... 42
5.9. The NOD5 Data Record - 5-digit node numbers ............................................................................... 43
6. Element Topology - ELM* (Data Category 2) ........................................................................................ 45
6.1. General Description ........................................................................................................................ 45
6.2. Data Category Header ..................................................................................................................... 46
6.3. Element Topology Data Record ....................................................................................................... 46
6.3.1. The QPPL Element .................................................................................................................. 48
6.3.2. The TPPL Element ................................................................................................................... 48
6.3.3. The TUBE Element .................................................................................................................. 49
6.3.4. The STUB Element .................................................................................................................. 49
6.3.5. The DISC Element ................................................................................................................... 50
6.3.6. The PMAS Element ................................................................................................................. 50
6.3.7. The PBOY Element .................................................................................................................. 50
6.3.8. The FPNT Element .................................................................................................................. 50
6.3.9. Examples of Element Specifications ........................................................................................ 51
6.4. SYMX and SYMY Data Records - X and Y Symmetry .......................................................................... 52
6.5.The HYDI Data Record - Hydrodynamic Interaction .......................................................................... 53
6.6. The RMXS/RMYS Data Records- Remove Symmetry .......................................................................... 54
6.7. The MSTR Data Record - Move Structure .......................................................................................... 54
6.8. The FIXD Data Record - Fix Structure ................................................................................................ 55
6.9. The VLID Data Record - Suppression of Standing Waves. ................................................................... 55
6.10. The ASYM Data Record - Axisymmetric Structure Generation ......................................................... 56
6.11.The ILID Data Record - Suppression of Irregular Frequencies ........................................................... 58
6.12. The ZLWL Data Record - Waterline Height ...................................................................................... 59
6.13. The SEAG Data Records - Creation of Wave Grid Pressures .............................................................. 60
7. Material Properties - MATE (Data Category 3) ...................................................................................... 63
7.1. General Description ........................................................................................................................ 63
7.2. Data Category Header ..................................................................................................................... 63
7.3. Material Property Data Record ........................................................................................................ 63
7.4.The STRC Data Record - Material Structure Association ..................................................................... 64
8. Geometric Properties - GEOM (Data Category 4) ................................................................................. 65
8.1. General Description ........................................................................................................................ 65
8.2. Data Category Header ..................................................................................................................... 65
8.3. Geometric Property Data Record ..................................................................................................... 65
8.4. The STRC Data Record - Geometry Structure Association .................................................................. 67
9. Constant Parameters - GLOB (Data Category 5) .................................................................................... 69
9.1. General Description ........................................................................................................................ 69
9.2. Data Category Header ..................................................................................................................... 69
9.3. The DPTH Data Record (Optional) - Water Depth .............................................................................. 69
9.4. The DENS Data Record (Optional) - Water Density ............................................................................ 69
9.5. The ACCG Data Record (Optional) - Acceleration Due to Gravity ....................................................... 70
10. Frequency and Directions Table - FDR* (Data Category 6) ................................................................. 71
10.1. General Description ...................................................................................................................... 71
10.2. Data Category Header ................................................................................................................... 71
10.3.The FREQ/PERD/HRTZ Data Record - Frequencies/Periods at which the Hydrodynamic Parameters
are Computed ...................................................................................................................................... 72
10.4. The DIRN Data Record - Directions at which the Hydrodynamic Parameters are Computed ............. 74
10.5. The MVEF Data Record - Move Existing Frequency Parameters ....................................................... 75
10.6.The DELF Data Record - Delete Frequency Parameters .................................................................... 77
10.7.The CSTR Data Record - Copying Existing Hydrodynamic Parameters for a Specific Structure Number ....................................................................................................................................................... 78
10.8. The FILE Data Record - Copying Existing Hydrodynamic Parameters from an External File ............... 78
iv
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Reference Manual
10.9. The CPYF/CPYP/CPYH Data Records - Copy Frequency Parameters ................................................. 79
10.10. The CPYS Data Record - Copy Stiffness Matrix .............................................................................. 81
10.11. The CPDB Data Record - Copy Existing Hydrodynamic Data Base .................................................. 81
10.12. The FWDS Data Record - Define Forward Speed ........................................................................... 82
10.13. Database Import ......................................................................................................................... 82
11. Wave Frequency Dependent Parameters and Stiffness Matrix - WFS* (Data Category 7) .................. 85
11.1. General Description ...................................................................................................................... 85
11.2. Data Category Header ................................................................................................................... 85
11.3. The LSTF Data Record - Linear Hydrostatic Stiffness Matrix ............................................................. 86
11.4. The ZCGE Data Record - Z Coordinate of the Center of Gravity at Equilibrium .................................. 87
11.5. The BFEQ Data Record - Buoyancy Force at Equilibrium .................................................................. 88
11.6. The FREQ/PERD/HRTZ Data Record - Frequencies/Periods at which the Parameters are Defined ...... 88
11.7.The WAMS Data Record - Wave Frequency Added Mass Matrix and the WDMP Data Record - Wave
Frequency Damping Matrix ................................................................................................................... 89
11.8.The WDGA Data Record - Wave Frequency Diagonal Added Mass and the WDGD Data Record Wave Frequency Diagonal Damping ..................................................................................................... 90
11.9. The AAMS/ADMP Data Record - Additional Added Mass/Damping Matrices ................................... 91
11.10.The FIAM/FIDP Data Record - Frequency Independent Additional Added Mass/Damping
Matrices ............................................................................................................................................... 92
11.11.The FIDA/FIDD Data Record - Frequency Independent Additional Diagonal Added Mass/Damping ....................................................................................................................................................... 92
11.12. The DIRN Data Record - Directions at which the Parameters are Defined ....................................... 93
11.13.The TDIF/RDIF/TFKV/RFKV/TRAO/RRAO Data Records - Wave Frequency Diffraction Forces and
Responses ............................................................................................................................................ 94
11.14. The GMXX/GMYY Data Record - User-Specified Metacentric Height .............................................. 95
11.15. The SSTR Data Record - Submerged Structure .............................................................................. 95
11.16. The ASTF Data Record - Additional Hydrostatic Stiffness Matrix .................................................... 96
11.17. The SSTF/SPOS/SFRC Data Records - Additional Structural Stiffness .............................................. 97
11.17.1. The SSTF Data Record - Additional Structural Stiffness Matrix ............................................... 97
11.17.2. The SPOS Data Record - Equilibrium Position ....................................................................... 98
11.17.3. The SFRC Data Record - Force at Equilibrium Position .......................................................... 98
12. Drift Force Coefficients - DRC* (Data Category 8) ............................................................................. 101
12.1. General Description .................................................................................................................... 101
12.2. Data Category Header ................................................................................................................. 102
12.3.The FREQ/PERD/HRTZ Data Records - Frequencies/Periods at which the Drift Coefficients are
defined ............................................................................................................................................... 102
12.4. The DRFX/DRFY/DRRZ Data Records - Drift Force Coefficients ....................................................... 103
12.5. The CQTF Data Record - Import of QTF Database .......................................................................... 104
13. Drift Motion Parameters - DRM* (Data Category 9) .......................................................................... 107
13.1. General Description .................................................................................................................... 107
13.2. Data Category Header ................................................................................................................. 107
13.3.The DGAM Data Record - Diagonal Added Mass Matrix ................................................................. 108
13.4. The DGDP Data Record - Diagonal Damping Matrix ...................................................................... 108
13.5.The LFAD Data Record - Use Lowest Wave Frequency Added Mass and Damping for Drift Frequency
Analysis .............................................................................................................................................. 109
13.6. The YRDP Data Record - Yaw-Rate Drag Parameters ...................................................................... 109
13.7.The FIDA/FIDD Data Record - Frequency Independent Additional Diagonal Added Mass and
Damping ........................................................................................................................................... 110
13.8.The FIAM/FIDP Data Record - Frequency Independent Additional Added Mass/Damping
Matrices ............................................................................................................................................. 111
14. Hull Drag Coefficients and Thruster Forces - HLD* (Data Category 10) ............................................ 113
14.1. General Description .................................................................................................................... 113
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Reference Manual
14.2. Data Category Header ................................................................................................................. 113
14.3.The CUFX/CUFY/CURZ Data Record - Current Force Coefficients and the WIFX/WIFY/WIRZ Data Record
- Wind Force Coefficients ..................................................................................................................... 113
14.4. The THRS Data Record - Thruster Forces ....................................................................................... 115
14.5.The NLRD/BOFF/BASE Data Records - Non Linear Roll Damping .................................................... 116
14.6. The MDIN/MDSV Data Records - Morison Hull Drag ...................................................................... 117
14.6.1. The MDIN Data Record ........................................................................................................ 117
14.6.2. The MDSV Data Record ....................................................................................................... 118
14.7. The IUFC/RUFC Data Records - External Forces ............................................................................. 118
14.8. The DIRN Data Record - Directions at which the Drag Coefficients are defined .............................. 119
14.9. SYMX and SYMY Data Records - X AND Y Symmetry ..................................................................... 120
14.10. The DPOS Data Record - Position for Drag calculation ................................................................. 120
14.11. The DDEP Data Record - Depth for current measurement ........................................................... 121
15. Current/Wind Parameters - ENVR (Data Category 11) ...................................................................... 123
15.1. General Description .................................................................................................................... 123
15.2. Data Category Header ................................................................................................................. 123
15.3. CURR/WIND Data Records - Uniform Current/Wind Velocity .......................................................... 123
15.4. CPRF Data Record - Profiled Current Velocity ................................................................................ 124
15.5. CDRN Data Record - Current Direction ......................................................................................... 124
15.6. The TOWS Data Record - Tether Tow Speed .................................................................................. 125
16. Motion Constraints on Structures - CONS (Data Category 12) .......................................................... 127
16.1. General Description .................................................................................................................... 127
16.2. Data Category Header ................................................................................................................. 127
16.3. The DACF Data Record - Deactivate Freedom at the Center of Gravity ........................................... 128
16.4. The DCON Data Record - Define Constraint Position ..................................................................... 129
16.5. The KCON/CCON/FCON Data Records - Define Articulation Stiffness, Damping and Friction ........... 132
16.5.1. KCON and CCON Data Records ............................................................................................ 132
16.5.2.The FCON Data Record ........................................................................................................ 133
17. Wind/Wave Spectrum Definition - SPEC (Data Category 13) ............................................................ 135
17.1. General Description .................................................................................................................... 135
17.2. Data Category Header ................................................................................................................. 135
17.3. Wind Spectra Definition .............................................................................................................. 135
17.4. The HRTZ Data Record - Change Units of Frequency to Hertz ........................................................ 137
17.5. The RADS Data Record - Change Units of Frequency to Radians/Second ....................................... 138
17.6. The SPDN Data Record - Wave Spectral Direction ......................................................................... 138
17.7. The SEED Data Record - Wave Spectral Seed ................................................................................. 139
17.8. The CURR / WIND Data Records - Current and Wind Speed and Direction ...................................... 139
17.9. The PSMZ Data Record - Pierson-Moskowitz Spectrum ................................................................. 140
17.10. The GAUS Data Record - Gaussian Spectrum .............................................................................. 141
17.11. The JONS Data Record - JONSWAP Spectrum ............................................................................. 141
17.12. The JONH Data Record - JONSWAP Spectrum ............................................................................ 142
17.13.The UDEF Data Record - User-Defined Wave Spectrum and the FINI Data Record - User-Defined
Spectrum Separator ............................................................................................................................ 143
17.14. The NSPL Data Record - Number of Spectral Lines ...................................................................... 143
17.15. The NRST Data Record - Number of Spectral Rasters ................................................................... 144
17.16. Input of a Time History of Wind Speed and Direction .................................................................. 144
17.17. Input of a Time History of External Forces on Structures ............................................................. 145
17.18. The XSWL/GATP Data Records - Input of Cross Swell ................................................................... 146
17.18.1. The GATP Data Record - Use of peak period ....................................................................... 147
17.18.2. Where Cross-Swell is Implemented .................................................................................... 148
17.18.3. Where Cross-Swell is not Implemented .............................................................................. 148
17.19. The IWHT Data Record - Import of Wave Height Time History ...................................................... 149
vi
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Reference Manual
17.20. The SSDN/UDDS/SSWT Data Records - User-Defined Spread Seas Spectrum ............................... 152
17.21. The NAME Data Record - Spectrum Name .................................................................................. 154
17.22. The NODR Data Record - No Drift Forces .................................................................................... 154
17.23. The SPGR Data Record - Spectral Group Keyword ....................................................................... 154
17.24. The SGNM Data Record - Spectral Group Name .......................................................................... 156
17.25. The TOWS Data Record - Tether Tow Speed ................................................................................ 156
18. Regular Wave (Aqwa-Naut) - WAVE (Data Category 13N) ................................................................. 157
18.1. General Description .................................................................................................................... 157
18.2. Data Category Header ................................................................................................................. 157
18.3. The WAMP Data Record - The Regular Wave Amplitude ................................................................ 157
18.4. The PERD Data Record - Regular Wave Period ............................................................................... 157
18.5. The WDRN & WVDN Data Records - Wave Direction ...................................................................... 157
18.6. The AIRY Data Record - Linear Wave ............................................................................................. 158
18.7.The WRMP Data Record - Wave Ramp ........................................................................................... 158
19. Mooring Lines and Attachment Points - MOOR (Data Category 14) ................................................. 161
19.1. General Description .................................................................................................................... 161
19.2. Data Category Header ................................................................................................................. 161
19.3. The LINE/WNCH/FORC Data Records - Linear Cables ..................................................................... 162
19.3.1. The LINE Data Record - A Conventional linear elastic Cable .................................................. 163
19.3.2.The WNCH Data Record - A Winch Adjusted to Constant Tension .......................................... 164
19.3.3. The FORC Data Record - Constant Force .............................................................................. 164
19.4. The POLY Data Record - Polynomial Nonlinear Properties ............................................................. 165
19.5. The NLIN Data Record - Nonlinear Cables ..................................................................................... 166
19.5.1. Parameters Applied to a Polynomial (POLY) Mooring Line .................................................... 166
19.5.2. Parameters Applied to a Composite Catenary (COMP/ECAT) Mooring Line ........................... 167
19.5.3. Parameters Applied to a Steel Wire (SWIR) Nonlinearity ....................................................... 168
19.6. The COMP/ECAT Data Records - Composite Catenary Mooring Line .............................................. 168
19.6.1. The COMP Data Record ....................................................................................................... 168
19.6.2.The ECAT Data Record ......................................................................................................... 170
19.7. The ECAX Data Record - Nonlinear Axial Stiffness Polynomial ....................................................... 171
19.8. The BUOY/CLMP Data Records - Intermediate Buoys and Clump Weights ...................................... 172
19.9. Cable Dynamics - Additional Data Requirements ......................................................................... 173
19.9.1. Introduction to Cable Dynamics .......................................................................................... 173
19.9.2. The NLID Data Record - Non-linear Dynamic Cable .............................................................. 174
19.9.3. The ECAH Data Record - Elastic Catenary Hydrodynamic Properties ..................................... 175
19.9.4. The ECAB Data Record - Elastic Catenary Bending Stiffness .................................................. 175
19.9.5. The NCEL Data Record - Number of Cable Elements ............................................................ 176
19.9.6. The DYNM/DOFF Data Records - Cable Dynamics ON/OFF Switch ........................................ 176
19.10. The LBRK Data Record - Mooring Line Break ............................................................................... 177
19.11. The PULY Data Record - Linear Cables ........................................................................................ 178
19.12. The DWT0/LNDW/DWAL Data Records - Linear Drum Winches ................................................... 180
19.12.1. The DWT0 Data Record ..................................................................................................... 180
19.12.2. The LNDW Data Record ..................................................................................................... 181
19.12.3. The DWAL Data Record - Drum Winch Additional Length ................................................... 183
19.13. The FEND/FLIN Data Records - Fenders ...................................................................................... 183
19.13.1. The FEND Data Record - Fender properties ........................................................................ 183
19.13.2. The FLIN Data Record - Fender connections ....................................................................... 184
19.14. The LE2D/*RNG Data Records - 2-Dimensional Load Extension Database .................................... 185
19.14.1.The LE2D Data Record ....................................................................................................... 185
19.14.2. The ZRNG/XRNG/HRNG/VRNG Data Records ...................................................................... 187
19.15.The PRIC Data Record - Print Initial Condition of Mooring Lines ................................................... 189
19.16. The FINI Data Record - Mooring Configuration Separator ............................................................ 189
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
vii
Reference Manual
19.17.The FILE Data Record - Mooring Description ............................................................................... 189
19.18. The SWIR Data Record - Steel Wire Nonlinear Properties ............................................................. 190
19.19. Tether Additional Data Requirements ........................................................................................ 191
19.19.1.The TELM Data Record - Tether Element ............................................................................. 191
19.19.2. The TSPA/TSPV Data Records - Tether Anchor and Vessel Springs ........................................ 192
19.19.3. The TEIG Data Record - Tether Eigensolution ...................................................................... 192
19.19.4. The TFAT Data Record - Tether Fatigue Parameters (Aqwa-Drift only) .................................. 193
19.19.5. The TPSH Data Record - Tether Peak Stress Hours (Aqwa-Drift only) .................................... 194
19.19.6. The TSLK Data Record - Tether Printing when Slack ............................................................ 194
19.19.7. The TEGR Data Record - Tether Group Factor ...................................................................... 195
19.19.8.The TCAP Data Record - Tether End Cap Areas .................................................................... 195
19.19.9. The TIFL Data Record - Tether Internal Fluid Properties ....................................................... 196
19.19.10. The TIMP Data Record - Tether Impact Parameters ........................................................... 196
19.19.11. The TLOW Data Record - Tether Lower Stop Position ........................................................ 197
19.19.12. The TETH Data Record - Tether Vessel and Anchor/Trailing End Position ............................ 197
19.19.13. The TLAC/TROC/TLAV/TROV Data Records -Tether Constraint Data Records ...................... 198
20. Starting Conditions Definition - STRT (Data Category 15) ................................................................ 201
20.1. General Description .................................................................................................................... 201
20.2. Data Category Header ................................................................................................................. 202
20.3. The POS* Data Record - Starting Positions .................................................................................... 202
20.4. The REA* Data Record - Initial Articulation Reactions .................................................................... 203
21. Starting Conditions Definition (Aqwa-Drift) - STRT (Data Category 15D) ........................................ 205
21.1. General Description .................................................................................................................... 205
21.1.1. Analysis Type - Drift only ..................................................................................................... 205
21.1.2. Analysis Type - Drift + Wave Frequency (WFRQ on JOB data record) ..................................... 206
21.2. Data Category Header ................................................................................................................. 207
21.3. The POS* Data Record - Starting Positions .................................................................................... 207
21.4. The VEL* Data Record - Starting Velocities .................................................................................... 208
21.5. The SLP* Data Record - Slow Starting Positions ............................................................................ 208
21.6. The SLV* Data Record - Slow Starting Velocities ............................................................................ 209
22. Starting Conditions Definition (Aqwa-Naut) - STRT (Data Category 15N) ........................................ 211
22.1. General Description .................................................................................................................... 211
22.2. Data Category Header ................................................................................................................. 212
22.3. The POS* Data Record - Starting Positions .................................................................................... 212
22.4. The VEL* Data Record- Starting Velocities ..................................................................................... 212
23. Time Integration Parameters - TINT (Data Category 16) ................................................................... 215
23.1. General Description .................................................................................................................... 215
23.2. Data Category Header ................................................................................................................. 215
23.3. The TIME Data Record - Time Integration Parameters .................................................................... 215
23.4. The HOTS Data Record - For Hot-Start Run ................................................................................... 216
24. Aqwa-Librium Iteration Parameters - LMTS (Data Category 16B) .................................................... 219
24.1. General Description .................................................................................................................... 219
24.2. Data Category Header ................................................................................................................. 219
24.3. The MXNI - Maximum Number of Iterations ................................................................................. 219
24.4. The MMVE Data Record - Maximum Movement Per Iteration ........................................................ 219
24.5.The MERR Data Record - Maximum Error Allowable for Equilibrium ............................................... 220
24.6. The STRP Data Record - Output Stability Report ........................................................................... 220
25. Change Geometric/Mass Characteristics - GMCH (Data Category 16L) ............................................ 223
25.1. General Description .................................................................................................................... 223
25.2. Data Category Header ................................................................................................................. 224
25.3.The SCAL Data Record - Length Scale Factor ................................................................................. 224
25.4.The MASS Data Record - Mass Scaling Parameter .......................................................................... 224
viii
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Reference Manual
25.5. The INER Data Record - New Inertia Values ................................................................................... 225
25.6. The NCOG Data Record - New Center of Gravity ........................................................................... 225
25.7. The REFP Data Record - New Hydrodynamic Reference Point ........................................................ 225
26. Hydrodynamic Parameters for Non-Diffracting Elements - HYDC (Data Category 17) .................... 227
26.1. General Description .................................................................................................................... 227
26.2. Data Category Header ................................................................................................................. 227
26.3. The SC1/ Data Record - Scale Factor for Model Test Simulation ..................................................... 228
26.4. The DRGM/ADMM Data Record - Drag/Added Mass Multiplying Factor ......................................... 229
26.5. The SLMM Data Record - Slam Multiplying Factor ......................................................................... 230
27. Additional File Output - PROP (Data Category 18) ........................................................................... 231
27.1. General Description .................................................................................................................... 231
27.2. Data Category Header ................................................................................................................. 232
27.3.The NODE Data Record - Nodal Position for Listing file Output ...................................................... 232
27.4. The ALLM Data Record - All Motions ............................................................................................ 233
27.5. The PCGP Data Record - Print Center of Gravity Parameters .......................................................... 233
27.6. The PREV Data Record - Print Every nth Timestep ......................................................................... 234
27.7. The PRNT/NOPR Data Record - Print/No Print Parameters Values ................................................... 235
27.8. The PTEN Data Record - Print Cable Tensions ................................................................................ 237
27.9. The ZRON/ZROF Data Records - Print Z Coordinate Relative to Wave Surface ON/OFF ................... 237
27.10.The ZRWS Data Record - Print Z Coordinate Relative to Wave Surface .......................................... 238
27.11. The PPRV Data Record - Print POS Every nth Timestep ................................................................ 238
27.12. The GREV Data Record - Graphics Output Every nth Timestep ..................................................... 238
27.13. The PRMD Data Record - Print Mooring Drag .............................................................................. 239
27.14.The PMST Data Record - Print Mooring Section Tensions ............................................................. 239
27.15. The SSPC Data Record - Print Sub-Spectrum Response ............................................................... 240
27.16. Tether Additional Data Requirements ........................................................................................ 240
27.16.1. The TPRV Data Record - Tether Printing Interval ................................................................. 241
27.16.2. The TGRV Data Record - Tether Graphics/Statistics Interval ................................................. 241
27.16.3. The TSTS/TSTF Data Record - Tether Start/Finish Timesteps for Statistics ............................. 241
28. Element and Nodal Loads - ENLD (Data Category 21) ...................................................................... 243
28.1. General Description .................................................................................................................... 243
28.2. Data Category Header ................................................................................................................. 243
28.3. The ISEL and LSEL Data Records - Element/Nodal Load Record Selection ...................................... 243
28.4. The RISR Data Record - Nodal Load Output for a Riser Structure .................................................... 244
28.5. The TUBE Local Axes .................................................................................................................... 244
29. Options for Use in Running Aqwa Programs ..................................................................................... 247
29.1. Command Line Options .............................................................................................................. 247
29.2. Job Type Options ........................................................................................................................ 247
29.3. Program Options for Use in Aqwa Program Suite .......................................................................... 248
29.3.1. Printing of the Expanded Input Data List for each Data Category ......................................... 248
29.3.2. Printing Options for Output of Calculated Parameters ......................................................... 250
29.3.3. Administration and Calculation Options for the Aqwa Suite ................................................. 252
30. External Force Calculation ................................................................................................................ 259
30.1. Example of Fortran Subroutine user_force ................................................................................... 259
30.2. Conversion of Euler Angles to LSA ............................................................................................... 261
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
ix
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 1: Introduction
The Aqwa suite is a set of advanced hydrodynamic analysis programs. This document is the Reference
Manual for the Aqwa suite of programs.
The programs in the suite are:
Aqwa-Line
Aqwa-Librium
Aqwa-Fer
Aqwa-Drift
Aqwa-Naut
The Reference Manual defines the input data format for the above programs. Most of the input information is applicable to more than one program. When the information relates to only one program, this
is made clear in the text.
The manual also contains details of certain topics common to all the programs.
The Aqwa Reference Manual should be used in conjunction with the appropriate individual Program
User Manual when running any of the programs in the Aqwa suite.
1.1. Conventions
The following conventions are adopted in Aqwa programs. Users should keep them in mind during input
data preparation and results interpretation.
Local System Axes (LSA)x,y,z with origin through the body's center of gravity
The Local System Axis (LSA) has its origin at the CoG of the vessel, with X, Y and Z axes parallel to the
FRA (see 1 above) when the vessel is in its initial definition position. With conventional modelling, X is
along the length of the vessel, Y along the beam to port, and Z in the direction of the cross product of
X and Y. This axis system moves with the vessel.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Introduction
Figure 1.1: Axes Systems
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Introduction
3. Motions of the structure are often output in the Aqwa suite in the form of translations and rotations
about the X, Y and Z axis and are termed surge, sway, heave, roll, pitch, and yaw. These motions may
be successive 'LARGE' rotations about the X, Y and Z axes of the FRA axes system in the definition position
(Aqwa-Drift/Librium/Naut) or 'small' rotations in arbitrary order about the analysis position in the local
LSA axes system (Aqwa-Line/Drift/Fer). It can be seen that interpretation of the results can be made extremely difficult by an unsuitable choice of the definition position of the structure.
For example, if the structure is a ship or barge, conventional terminology for motion along, and rotation about, the longitudinal center line, is surge and roll. However, if the longitudinal center line
is defined parallel to the FRA Y-axis, then rotational motion about this axis will be termed pitch, and
translational motion along the axis, sway.
For other structures, e.g. semi-submersibles, this may not be relevant. The user must take due note
of the terms associated with the motions about the axes and is recommended to define all ship/barge
shaped structures with the longitudinal axis parallel to the FRA X-axis.
1.3. Units
The units used in Aqwa are decided by the input values for the water density and gravitational acceleration. For example, if metre, Newton are to be used as the units for the length and force, users should
use 1025 for the water density and 9.806 for the gravitational acceleration.
In the output, the unit for the rotational motions is in degrees (although they are originally calculated
in radians), while the rotational terms in the stiffness and damping matrices are output in radians.
The user is free to choose any system of units for the data, with the proviso that the system must be
consistent. This means that the unit of mass must be consistent with the units of length and force
already selected.
Examples of consistent sets of units are:
SI units Force in newtons, length in metres, mass in kilograms, time in seconds, acceleration in meters/sec2
Imperial units Force in poundals, length in feet, mass in pounds, time in seconds, acceleration in feet/sec2,
or Force in pounds, length in feet, mass in slugs, time in seconds, acceleration in feet/sec2
For any other set of units, the consistent unit of mass will be a multiple of the basic unit of mass because
it is a derived unit.
The consistent unit of mass is obtained by dividing the unit of force by the acceleration due to gravity,
which itself has units of length divided by time squared. A change in the unit of length, for example,
from feet to inches or metres to millimeters, requires a corresponding change in the unit of mass used
for calculating the density. A list of sets of consistent units is given below.
_____________________________________________________________________________
|
|
|
|
|
|
|
| Unit of | Unit of | Typical |
| Consistent | Density (mass/unit volume) |
| force
| length | value of| g | unit of
| -------------------------- |
|
|
| E for
|
| mass
|
Steel
| Sea Water
|
|
|
| Steel
|
|
|
|
|
|-----------------------------------------------------------------------------|
| Newton | metre
| 2.1*E11 | 9.81| 1.0 Kg
| 7850
| 1025
|
|
|
|
|
|
|
|
|
| Newton | cm
| 2.1*E07 | 981| 100 Kg
| 7.85*E-5
| 1.025*E-5
|
|
|
|
|
|
|
|
|
| Newton | mm
| 2.1*E05 | 9810| 1000 Kg
| 7.85*E-9
| 1.025*E-9
|
|-----------------------------------------------------------------------------|
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Aqwa Files
| Kilopond| metre
| 2.14*E10| 9.81| 9.81 Kg
| 800
| 104.5
|
|
|
|
|
|
|
|
|
| Kilopond| cm
| 2.14*E6 | 981| 981 Kg
| 8.00*E-06 | 1.045*E-06
|
|
|
|
|
|
|
|
|
| Kilopond| mm
| 2.14*E4 | 9810| 9810 Kg
| 8.00*E-10 | 1.045*E-10
|
|-----------------------------------------------------------------------------|
| KNewton | metre
| 2.1*E8 | 9.81| 1000 Kg
| 7.85
| 1.025
|
|
|
|
|
|
|
|
|
| KNewton | cm
| 2.1*E4 | 981| 1.0*E5 Kg | 7.85*E-08 | 1.025*E-08
|
|
|
|
|
|
|
|
|
| KNewton | mm
| 2.1*E2 | 9810| 1.0*E6 Kg | 7.85*E-12 | 1.025*E-12
|
|-----------------------------------------------------------------------------|
| Tonne(f)| metre
| 2.14*E7 | 9.81| 9.81*E3 Kg | 0.800
| 0.1045
|
|
|
|
|
|
|
|
|
| Tonne(f)| cm
| 2.14*E3 | 981| 9.81*E5 Kg | 8.0*E-09
| 1.045*E-09
|
|
|
|
|
|
|
|
|
| Tonne(f)| mm
| 2.14*E1 | 9810| 9.81*E6 Kg | 8.0*E-13
| 1.045*E-13
|
|-----------------------------------------------------------------------------|
| Poundal | Foot
| 1.39*E11| 32.2| 1.0 lb
| 491
| 64.1
|
|
|
|
|
|
|
|
|
| Poundal | Inch
| 9.66*E8 | 386 | 12 lbs
| 2.37*E-2
| 3.095*E-3
|
|-----------------------------------------------------------------------------|
| Pound(f)| Foot
| 4.32*E9 | 32.2| 32.2 lbs
| 15.2
| 1.985
|
|
|
|
|
|
|
|
|
| Pound(f)| Inch
| 3.0*E7 | 386 | 386 lbs
| 7.35*E-4
| 9.597*E-5
|
|-----------------------------------------------------------------------------|
| Kip
| Foot
| 4.32*E6 | 32.2| 3.22*E4 lbs| 1.52*E-2
| 1.985*E-3
|
|
|
|
|
|
|
|
|
| Kip
| Inch
| 3.0*E4 | 386 | 3.86*E5 lbs| 7.35*E-7
| 9.597*E-8
|
|-----------------------------------------------------------------------------|
| Ton(f) | Foot
| 1.93*E6 | 32.2| 7.21*E4 lbs| 6.81*E-3
| 8.892*E-4
|
|
|
|
|
|
|
|
|
| Ton(f) | Inch
| 1.34*E4 | 386 | 8.66*E5 lbs| 3.28*E-7
| 4.283*E-8
|
-----------------------------------------------------------------------------
Note
1 Kip = 1000 pounds force
1 Kilopond = 1 Kilogram force
all times are in seconds
assumed specific gravity of steel = 7.85
assumed specific gravity of sea water = 1.025
Input Files
DAT -- ASCII file for model definition and analysis parameters. Used by all Aqwa programs.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Introduction
XFT -- ASCII file defining a time history of external force on a structure or structures in six degrees of
freedom in local axis system. Used for time domain analysis (optional).
WVT -- ASCII file defining a time history of wind velocity and direction. Used for time domain analysis
(optional).
WHT -- ASCII file defining a time history of water surface elevation. Used for time domain analysis (optional).
LIN -- ASCII file defining ship offsets. Used by AGS Mesh Generator to define hull shape.
MSD -- ASCII file defining the mass distribution of a vessel. Used by AGS for shear force and bending
moment calculation.
SFM -- ASCII file defining the mass distribution of a vessel. Used by AGS for splitting force calculation.
EQP -- Binary file containing the equilibrium positions of structures. Created by Aqwa-Librium and used
(optional) by Fer, Drift or Naut; see RDEP option (Administration and Calculation Options for the Aqwa
Suite (p. 252)).
Output Files
LIS -- ASCII file containing model definition/analysis parameters and the analysis results.
MES -- ASCII file containing messages issued during an Aqwa analysis.
QTF -- ASCII file containing fully populated matrix of Quadratic Transfer Functions.
HYD -- Binary file containing the hydrodynamic results calculated in Aqwa-Line. Can be used for further
Aqwa analysis.
RES -- Binary file containing the model definition/analysis parameters and the hydrodynamic results calculated in Aqwa-Line. Can be used for further Aqwa analysis or structure visualisation etc in AGS.
EQP -- Binary file containing the equilibrium positions of structures. Created by Aqwa-Librium and used
(optional) by Fer, Drift or Naut; see RDEP option (Administration and Calculation Options for the Aqwa
Suite (p. 252)).
ENL -- Binary file containing Morison element/nodal loading. Only created for tether analysis, or tube
elements at analysis stage 6.
POS -- Binary file containing structures' positions at each time step. Used by AGS for generating animation.
PLT -- Binary file containing Aqwa analysis results. Used by AGS for plotting graphs.
POT -- Binary file containing potentials. Used by AGS or Aqwa-WAVE for element pressure calculation. This
file is only created by Aqwa-Line when LDOP option is on.
USS -- Binary file containing source strengths. Used by Aqwa-WAVE for Morison force calculation. This file
is only created by Aqwa-Line when LDOP option is on.
SEQ -- Binary file containing the animation of structure motion. Created and used by AGS.
TAB -- ASCII file containing the statistics table from Aqwa-Drift tether analysis.
PAC -- Binary file containing pressures at element centroids. Used by AGS for post-processing involving
pressures.
6
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Aqwa Files
VAC -- Binary file containing fluid velocities at element centroids. Used by AGS for wave contour plotting.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Drag Linearization
There are several types of nonlinear drag in Aqwa, most of which can be linearized and included in the
linear programs Aqwa-Line and Aqwa-Fer. TUBE, DISC and STUB elements have drag calculated according
to Morison's equation and can be included in both Aqwa-Line and Aqwa-Fer. Hull drag is modeled with
drag force coefficients input in Data Category 10. As with most other Stage 4 data, this is not used in
Aqwa-Line but can be included in Aqwa-Fer.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Analysis Stages
Aqwa-Line
The linearization is based on a spectrum that must be input in Data Category13. There can be only one
spectrum; i.e. a spectral group cannot be used for linearization. The spectral direction is ignored and the
same spectrum is used with each of the wave directions specified in Data Category 6.
In the linearization process the drag force term CD(1/2)v2 is replaced by CD(1/2)vVRMS where
VRMS is the rms velocity in the specified spectrum and is a factor calculated so that the total energy
dissipated is the same. Because the values of and VRMS depend on the spectrum, the revised RAOs
are strictly only applicable in that spectrum.
Aqwa-Line must be run for stages 1 to 5, either in one run or two.
If there are multiple structures they cannot be connected to each other by an additional stiffness
matrix. Connection to the ground ("structure zero") is allowed.
New RAOs including the linearized drag are output to the .LIS file, and are also available for
plotting in the AGS. Note that the new RAOs are NOT written to the database and are NOT used
for calculation of the 2nd order drift coefficients.
Aqwa-Fer
In Aqwa-Fer, the linearized drag will be applied in all spectra in all spectral groups. The same linearization
process is utilized as in Aqwa-Line.
10
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Restart Stage 1
11
Analysis Stages
No Axial Motion - Towed tethers are not considered to move in the axial direction or rotate about the
axis of the tether; i.e., displacements of the tether are 2 translations and 2 rotations at each node. These
displacements are considered as small motions from the tether axis (TLA q.v.)
Note
Although current in the axial direction will produce stabilizing effects, if the tether spring
at the ends are very soft, large rotations (>30 degrees) may be produced, which will invalidate the analysis. The program also takes full account of the change in encounter frequency,
due to the component of the current in the direction of the waves.
Axial Tension - Both the wall and effective tensions in a towed tether are assumed to be zero, and hence
the bending stiffness is purely structural. The tether responses, especially in the fundamental mode, may
be inaccurate if this tension is significant.
Note
This also means that the tether may not be analyzed, if any point moves to a depth where
the effective tension is significant, i.e. for upending.
Small Motions - It is assumed that the lateral and rotational motions of the tether from the defined tether
axis are small. This means that the program is unsuitable for large rotations about the Y or Z axis, e.g. for
upending. However, full account is taken of the phase shift of the waves, due to movement in the direction
of the wave/wave spectrum.
Mass/Stiffness - The mass/stiffness ratio of any element must not be too small. Very short elements inherently have small mass/stiffness ratios. This gives rise to very high frequencies. These high frequencies may
cause stability problems and roundoff errors in the programs. A general rule is that natural periods of less
than 1/100th second are not allowed. These periods are output from the eigenvalue analysis.
Very short elements should therefore be modelled with a value of Young's modulus reduced so that
no periods less than 1/100th second are present. The user can check that the bending of short elements
is still small, using the graphical output.
Timestep - The timestep must be small enough to resolve the response motion of the tether. This includes
any transients that may be present either initially or, more importantly, throughout the analysis. Although
a good rule of thumb is that the timestep should be 1/10th of the period of any response, the best
method of checking the timestep is to re-run a short simulation with half the timestep and compare the
bending moments or stresses for both runs. These should be approximately the same for both runs.
Timesteps of 0.25 seconds are typically used.
For towed tethers, the local axis (TLA) must be defined parallel to, and in the same direction as, the
X axis of the fixed reference axes (FRA) i.e. XY in the water plane and Z vertical. The X axis coincides
with the zero current wave direction. The nodes of the tether increase with positive X. The last node
of the tether, at zero FRA displacement, lies at the TLA origin. For installed tethers, the TLA is parallel
to the FRA, when the tether is vertical. In general, the TLA X axis goes from the anchor node to the
attachment node, the Y axis is in the plane of the XY FRA, and the Z axis follows the right hand rule.
The TLA origin is at the anchor node.
12
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Restart Stage 1
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
13
Analysis Stages
14
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Restart Stage 2
Diffraction Forces
Froude-Krylov Forces
Response Motions (or RAOs)
Data Category 8 - No Input or Drift Coefficients
Usually, not all of the above data are required for a particular mode of analysis, in which case, the user
may simply omit the data which are not applicable.
Note
Although the Radiation/Diffraction parameters calculated by Aqwa-Line can be transferred
to other programs in the Aqwa suite, this is not mandatory. This means that, if the backing
file produced by an Aqwa-Line run is not available (i.e. Aqwa-Line has not been run previously)
or if the user wishes to input values from a source other than Aqwa-Line, he may do so in
these data categories.
The following sections show the required data input for the available modes of analysis:
15
Analysis Stages
Note that the frequencies input in Aqwa-Line runs must differ from those in previous runs or they will
be automatically rejected. As all parameters are defined for a unique range of directions, these directions
may not be redefined General Description (p. 71).
16
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Restart Stage 2
Data Category 8 - No input required
Input for Aqwa-Drift/Fer/Librium/Naut with results from a previous AqwaLine run and a source other than Aqwa-Line
The new user is advised to ignore this facility
If the user wishes to APPEND to or CHANGE the parameters calculated by a previous Aqwa-Line run for
the current analysis, this is achieved by using the backing file from a previous Aqwa-Line run (i.e.
automatically read), together with specifying the parameters to be appended or changed.
To APPEND to the parameters calculated in a previous run, additional frequencies which differ from
those existing may be input in Data Category 6, together with values of the appropriate frequency dependent parameters in Data Categories 7 and 8, at these additional frequencies. Note that, as all parameters are defined for a unique range of directions, these directions may not be re-defined (see Frequency and Directions Table - FDR* (Data Category 6) (p. 71)).
To CHANGE the parameters calculated in a previous run, these parameters are simply input in Data
Categories 7 and 8 and, depending on the type of input, (see Wave Frequency Dependent Parameters
and Stiffness Matrix - WFS* (Data Category 7) (p. 85) and Drift Force Coefficients - DRC* (Data Category
8) (p. 101)), the parameters will be either overwritten with the input values or become the sum of input
values and original values.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
17
Analysis Stages
Data Category 6 - No input required
Data Category 7 - Input Appended or Changed Data
Data Category 8 - No Input or Drift Coefficients
Aqwa-Line
Data Category 9 - Drift Added Mass and Damping (only used for scaling, Data Category 16)
18
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Restart Stage 4
Data Category 10 - No input required
Data Category 11 - No input required
Data Category 12 - No input required
Data Category 13 - No input required except:
Spectrum Information (PSMZ/JONS/JONH/UDEF), only used when linearizing Morison drag using the LDRG
option
Data Category 14 - No input required
Data Category 15 - No input required
Data Category 16 - Geometrical Changes
Length Scaling of Results
Mass Scaling of Results
Change of Body's Center of Gravity
Change of Body's Inertia
Definition of New Hydrodynamic Reference Point
Data Category 17 - No input required
Data Category 18 - Printing Options for Nodal RAOs
Aqwa-Librium
Data Category 9 - Drift Added Mass and Damping
Data Category 10 - Hull Drag
Current Drag Coefficients
Wind Drag Coefficients
Thruster Forces
Data Category 11 - Current Velocity Profile and Direction
Data Category 12 - Elimination of Degrees of Freedom
Articulations between Bodies
Articulations between Global Points and Bodies
Data Category 13 - Wave Spectra Details
Data Category 14 - Mooring Lines
Data Category 15 - Motion Analysis Starting Positions
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
19
Analysis Stages
Data Category 16 - Iteration Limits
Data Category 17 - No input required
Data Category 18 - Printing Options for Nodal Positions
Aqwa-Fer
Data Category 9 - Drift Added Mass and Damping
Data Category 10 - Hull Drag
Current Drag Coefficients
Wind Drag Coefficients
Thruster Forces
Data Category 11 - No input required
Data Category 12 - Elimination of Degrees of Freedom
Articulations between Bodies
Articulations between Global Points and Bodies
Data Category 13 - Wave Spectra Details
Data Category 14 - Hawser/Mooring Lines
Data Category 15 - Motion Analysis Starting Positions
Data Category 16 - No input required
Data Category 17 - No input required
Data Category 18 - Printing Options for Nodal Positions
Aqwa-Naut
Data Category 9 - No input required
Data Category 10 - Hull Drag
Current Drag Coefficients
Wind Drag Coefficients
Thruster Forces
Data Category 11 - Current Velocity, Direction and Profile
Wind Velocity and Direction
Data Category 12 - Elimination of Degrees of Freedom
20
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Restart Stage 4
Articulations between Bodies
Articulations between Global Points and Bodies
Data Category 13N - Regular Wave Properties - default analysis
Data Category 13 - Wave Spectrum Details - irregular wave analysis
Data Category 14 - Hawser/Mooring Lines
Data Category 15 - Motion Analysis Starting Positions
Data Category 16 - Time Integration Parameters
Data Category 17 - Hydrodynamic Scaling Factors for Morison Elements
Data Category 18 - Printing Options for Nodal positions
Aqwa-Drift
Data Category 9 - Drift Added Mass and Damping
Yaw Rate Drag Coefficient
Data Category 10 - Hull Drag
Current Drag Coefficients
Wind Drag Coefficients
Thruster Forces
Data Category 11 - Current/Wind Velocity and Direction (only when no spectrum)
Data Category 12 - Elimination of Degrees of Freedom
Articulations between Bodies
Articulations between Global Points and Bodies
Data Category 13 - Wave Spectrum Details
Data Category 14 - Hawser/Mooring lines
Data Category 15 - Motion Analysis Starting Positions for both Wave and Drift Frequency Motions
Data Category 16 - Time Integration Parameters
Data Category 17 - No input required
Data Category 18 - Printing Options for Nodal Positions
21
Analysis Stages
Data Category 9
No additional input required
Data Category 10
No additional input required
Data Category 11 (for towed tethers)
Current speed and direction
Tow speed and direction
Data Category 12
No additional input required
Data Category 13 (for towed tethers)
Tow speed and direction
Data Category 14
Description of tether elements and boundary conditions
Fatigue and extreme value information (for towed tethers in Aqwa-Drift only)
Data Category 15 (for towed tethers)
Initial position of the vessel, which must be the same as the coordinate of the center of gravity
Data Category 16
No additional input required
Data Category 17 (for towed tethers)
Slam coefficient multiplier if required
Data Category 18
Frequency of tether information output to both the listing file and the plot file
Start and finish times for the statistical/fatigue analysis (for towed tethers in Aqwa-Drift only)
22
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Restart Stage 6
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
23
24
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
_(1) Optional User Identifier (A2) (see Optional User Identifier (p. 27))
|
|
_(2) Compulsory Data Category Keyword (A4) (see Compulsory Data Category Keyword (p.
|
|
|
|
2
5|
11 |
- --- -- ---- ---- ---------------------|X|XXX| |XXXX|COOR|XXXXXXXXXXXXXXXXXXXXX )
- --- -- ---- ---- ---------------------|
.
)
.
|
.
)
.
|
.
) etc.
---------------------------------------|
numerical input.....
) to 80
---------------------------------------|
numerical input.....
) columns
- --- ---------------------------------|X|END|
numerical input.....
)
- --- ---------------------------------|
|
|
|_(3) Compulsory END Code (see Compulsory 'END' Statement (p. 27))
The 'X' characters in columns 1-4, 7-10, etc. indicate that no data is to be input in these columns.
Note that the data category keyword for Data Category 1 is 'COOR' (COORdinate positions). The
keyword will be different for other data categories.
The precise format for every data record in every data category is given in the following sections:
The Preliminary Data (Data Category 0) (p. 33)
Node Number and Coordinates (Data Category 1) (p. 37)
Element Topology - ELM* (Data Category 2) (p. 45)
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
25
26
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
27
28
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
29
Warning
THE DEFAULT VALUES SUPPLIED ARE IN S.I. UNITS
Therefore, where Aqwa returns non-zero default values (for physical quantities), only users employing
S.I. units may make use of the default facility.
In any event, all data either input by the user or assumed by the program are automatically sent to the
listing file, unless the user specifically requests otherwise.
30
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
31
32
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) The 4-letter code is for the convenience of the user and is not used by the program. If the user
wishes to have the same title for several runs, the 4-letter character identifier (which should always be
unique) enables the user to distinguish between these runs. This is particularly relevant when using the
graphics, as this code is output on every plot.
(2) An abbreviation of the program name must be input to specify the overall data input format to be
expected by the program. If left blank, or the incorrect name is input, the program will output an error
message and abort after the preliminary data category has been read. The accepted abbreviations for
the program names are as follows:
Aqwa-Drift = DRIF
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
33
34
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) Up to 14 options may be input on each options data record. Each data record must have the data
record keyword OPTIONS in columns 1-7 and the last option on the last OPTIONS data record must be
END as shown above. In practice, only a few options are used for any one run. It is therefore extremely
unusual to input more than one OPTIONS data record.
The list of options common to all programs, and those which are have special applicability for a certain
program are listed for each program in Options for Use in Running Aqwa Programs (p. 247).
1
9 12
21
------- - --- --- ------ ----------|RESTART|X|
|
|XXXXXX|FILENAME
------- - --- --- ------ ----------|
|
|
|
|
|
|
|_(3)Optional. Name of database to be copied (without extension)
|
|
|
|
|
|_(2)Finish Stage, i.e. Finish at the End of this Stage (I1)
|
|
|
|_(1)Start stage, i.e. Start at the Beginning of this Stage (I1)
|
|_Compulsory Data Record Keyword
(1) The start stage will depend on which stages have been previously run. Stages in the program may
not be omitted. This means that the finish stage of the previous run must either be the previous stage
or a later stage in the program analysis.
When errors are detected, the user may start the analysis again at any of the preceding stages before
the error occurred. If the previous stage was correctly performed, the user instructs the program to
execute the next stage, i.e. simply to continue the analysis.
Note that as no analysis is performed within Stage 3 (except for Aqwa-Line) a run starting and finishing
at Stage 3 may be used, in conjunction with the PRDL option, to print out the contents of the restart
file on the listing file. As the restart file is a binary file and cannot be printed out directly, this is very
useful if there is some question as to the contents.
See Analysis Stages (p. 9) for a description and list of restart stages for each program.
(3) The program will automatically copy any relevant files with this name to use for the current run.
This may be .RES, .EQP and/or .POS, depending on the program which is being run. For example, using
the name ABTEST1 in the Aqwa-Fer data file AFTEST2.DAT will result in:
ABTEST1.RES being copied to AFTEST2.RES
ABTEST1.EQP being copied to AFTEST2.EQP.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
35
36
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
37
(1) This is the node number (N1) whose coordinate position is given by the values X, Y, Z (4).
(2) Leave blank (as in first data record shown above) if only 1 node is to be generated. A total of N2
nodes are generated automatically if this number (N2) is greater than 1. The Nth node is generated:
with a node number N1+ (N-1)N3
with coordinates X+(N-1)DX, Y+(N-1)DY, Z+(N-1)DZ
(3) Leave blank (as in first data record shown above) if only one node is to be generated. If this field is
left blank when generating nodes, then a value of 1 is assumed for this Nodal Increment (N3).
(4) Coordinate Position (X,Y,Z) associated with node number N1 (1).
(5) Leave blank (as in first data record shown above) if only 1 node is to be generated. These 3 values
are referred to as DX,DY,DZ in (2).
38
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) This single letter code indicates that a set of nodes is to be generated by rotating a set of nodes
already input.
(2) This is the starting node of the set to be copied (input sequence = NS1). The finish node of the set
is the last node input so far (input sequence = NF1).
(3) This is the number of node sets to be generated (NSETS). The total number of nodes automatically
generated in addition to those already input is given by:
Total additional nodes = NSETS (NF1-NS1+1)
If this field is left blank, or zero is input, NSETS=1 is assumed.
(4) This node increment number (N3) is the number by which the node numbers already input are to
be incremented. In other words, if the nodes previously input are numbered NP1, NP2, NP3 etc., and
the new node numbers are NN1, NN2, NN3 etc., these new node numbers are given by:
NN1 = NP1 + N3; NN2 = NP2 + N3; NN3 = NP3 + N3; etc.
If no increment value is present, or zero is input, the program will issue a warning and default to a
value of 100.
(5) These three values are the coordinates of the point about which the rotation of the nodes already
input takes place. This field may be left blank, in which case the point is assumed to be 0.0, 0.0, 0.0, i.e.
the FRA origin.
(6) These three values (RX,RY,RZ) are the components of the angular vector of rotation, i.e. the angle
and axis about which the rotation takes place is given by
Angle = modulus (RX,RY,RZ)
Axis = vector (RX,RY,RZ) / Modulus (RX,RY,RZ)
Although this means that the user has the facility to generate nodes rotated about any axis, the most
common usage is as follows (using 30 degrees as an example)
Rotation about X-axis RX = 30.0 RY = 0.0 RZ = 0.0
Rotation about Y-axis RX = 0.0 RY = 30.0 RZ = 0.0
Rotation about Z-axis RX = 0.0 RY = 0.0 RZ = 30.0
EXAMPLE
Generation of nodes for six, 50m tall, 8m diameter cylindrical legs of an axisymmetric semi-submersible
at a radius of 75 metres
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
39
Note
If the STRC data record is not used all user node numbers must be unique.
If the STRC data record is used in the COOR data category, then it may also be used in the
MATE and GEOM data categories if the use of the same property number is required from different structures.
If the STRC data record is not used in the COOR Data Category, it cannot be used in the
MATE(materials) and GEOM(geometry) Data Categories.
40
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(3) The structure number must start from 1 on the 1st data record and increment by 1 for each structure
data record input. The number of ELM* data categories must correspond to the number of structure
data records.
(1) Note that upper case 'O' will define an absolute offset and lower case 'o' will define a relative offset.
All node coordinates input after this data record will be generated offset by these 3 values.
For example, if the nodes for 3 identical legs along the X axis, 30m apart, with the 1st leg at X = 5.0
are defined, then the 2 absolute X offset data records input just before the node coordinate for legs 2
and 3 would be 35.0,65.0. The 2 relative X offsets would be 30.0,30.0, as the default offset is zero.
Note
If a STRC data record (The STRC Data Record - Coordinate Structure Association (p. 40)) is
input, the OFFSET will automatically be re-set to zero.
41
(1) This single letter code indicates that new sets of nodes are to be generated by translating a set of
nodes already input.
(2) This is the starting node of the set (input sequence = NS1). The finish node of the set is the last
node input so far ( input sequence = NF1).
(3) This is the number of node sets to be generated (NSETS). The total number of nodes automatically
generated in addition to those already input is given by:
Total additional nodes = NSETS (NF1-NS1+1)
If this field is left blank, or zero is input, NSETS = 1 is assumed.
(4) This node increment number (N3) is the number by which the node numbers already input are to
be incremented. In other words, if the node set previously input are numbered NP1, NP2, NP3, etc., and
the new node numbers are NN1, NN2, NN3, etc., these new node numbers are given by:
NN1 = NP1 + N3; NN2 = NP2 + N3; NN3 = NP3 + N3
If no increment value is present, or zero is input, the program will issue a warning and default to a
value of 100.
(5) This node coordinate increment number (DX,DY,DZ) is the VALUE by which the node set coordinates
are to be incremented, i.e. if the node set previously input are values XP1, XP2, XP3 etc and the new X
coordinates XN1, XN2, XN3, etc, these new values are given by:
XN1 = XP1 + DX; XN2 = XP2 + DX; XN3 = XP3 + DX
Similarly for Y and Z.
42
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) This single letter code indicates that a set of nodes is to be generated by reflecting a set of nodes
already input.
(2) This is the starting node of the set to be copied (input sequence = NS1). The finish node of the set
is the last node input so far ( input sequence = NF1). The number of nodes generated will always be
(NF1-NS1+1).
(3) This node increment number (N3) is the number by which the node numbers already input are to
be incremented. In other words, if the nodes previously input are numbered NP1, NP2, NP3, etc., and
the new node numbers are NN1, NN2, NN3, etc., these new node numbers are given by:
NN1 = NP1 + N3; NN2 = NP2 + N3; NN3 = NP3 + N3; etc.
If no increment value is present, or zero is input, the program will issue a warning and default to a
value of 100.
(4) These are the co-efficients A, B, C, D which define the plane of reflection. The plane is defined by
the equation:
Ax + By + Cz = D
When the NOD5 data record is used the format for all following coordinate data records are changed
slightly as shown below. The field for node numbers becomes 5 columns, while the field for the number
of generated nodes is reduced to 4 columns. All the other fields are unchanged.
2
5 7
12
16
21
31
51
61
71
- --- -- ----- ---- ----- --------- -- -- --------- --------- ------|X|
| |
4|
|
| 54.6
| ... |
|
|
...
- --- -- ----- ---- ----- --------- -- -- --------- --------- ------|X|END| |
23|
3|
10| 54.6
| ... |
4.5
|
|
...
- --- -- ----- ---- ----- --------- -- -- --------- --------- ------| |
|
|
|
|
| |
|
|
|
|_(4)3 Coordinates for Start Node (3F10.0)
| |
|
|
|
| |
|
|
|_(3)Increment for Node Number(I5)
| |
|
|
| |
|
|_(2)Number of Nodes to be Generated(I4)
| |
|
| |
|_(1)Starting Node Number(I5)
| |
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
43
44
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Note
At this stage of the modeling (i.e. input of Data Category 2) the program does not know how
many structures the user wishes to define. Therefore The Data Category Series Terminator FINI (p. 31) must be used to indicate that no more structures will be input. However, if the
maximum number of structures is input, the FINI data record must be omitted, as termination
of the Data Category series is then mandatory.
The 32-bit version of Aqwa can accept up to 18000 elements explicitly defined in Data Category 2. Of
these, 12000 may be diffracting elements. These limits are raised to 36000/24000 elements with one
symmetry data record, or 72000/48000 elements with two symmetry data records. These limits are for
all the elements in a complete model, not for each individual structure.
The 64-bit version of Aqwa can accept up to 40000 elements explicitly defined in Data Category 2. Of
these, 30000 may be diffracting elements. These limits are raised to 80000/60000 elements with one
symmetry data record, or 160000/120000 elements with two symmetry data records. These limits are
for all the elements in a complete model, not for each individual structure.
Aqwa can accept up to 50 structures, each including diffracting or non-diffracting elements, but the
total number of elements on all structures must not exceed the limits above. It is not possible for all
these structures to interact hydrodynamically with each other, and the limits on the number of structures
are given below. If there are: S sets of interacting structures, with In structures in each (n = 1, S), D additional diffracting structures, and ND additional non-diffracting structures, then interaction limits are
defined by the following rules:
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
45
=
+ +
Rule 2
Rule 3
+ +
46
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| (3)Columns 21-80 free format data for
|
|
generating NODE,MATERIAL,and GEOMETRIC
|
|
GROUP NUMBERS. Format(60A1)
|
|
|
|_(2)Element Group Number(I5)
|
|_(4) DIFF identifier (diffracting plates only)
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) The element type (always four characters) provides the classification for a particular element, i.e. the
number of nodes, and whether material and geometric properties are required. Valid codes for element
types are QPPL, TPPL, TUBE, PMAS, PBOY, FPNT, DISC and STUB as shown below:
Description
No of Nodes
Material Property
Geometric Property
QPPL
Quadrilateral Panel
none
none
TPPL
Triangular Panel
none
none
TUBE
Tube
density
geometry
PMAS
Point Mass
mass
inertia
PBOY
Point Buoyancy
displaced mass
none
FPNT
Field Point
none
none
STUB
Slender Tube
mass, inertia
geometry
DISC
Circular Disc
none
geometry
(2) The element group number associated with each element is used to divide the elements defining
the structure into groups. Groups can be used for plotting and to identify special sets of elements (see
The ILID Data Record - Suppression of Irregular Frequencies (p. 58) and The VLID Data Record - Suppression of Standing Waves. (p. 55)), but the group number may be left blank if desired.
(3) Free Format data generation is achieved by specifying several bracketed sets of Topological variables
in columns 21-80. The number of bracketed sets is given by
Number of Sets = 1 + (Number of Nodes)
+ 1 if a Material Group is required
+ 1 if a Geometry Group is required
In general we have a format of
(N)(N1,N2,N3)
- -------|
|
first
node
(N1,N2,N3)
-------|
|
second node
or parameter
(N1,N2,N3).....
-------|
|
third node
or parameter
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
47
Note
The nodes defining TPPL and QPPL elements must be ordered in a counterclockwise direction
from the perspective of an observer external to the structure.
In an Aqwa-Line data file, QPPL and TPPL elements which are below the still water line in the
Aqwa-Line analysis position must be denoted as diffracting elements by entering the identifier
DIFF in columns 12 to 15.
For structures which cross the waterline, the top row of diffracting plate elements must have
their top edges aligned with the still water line (i.e. diffracting plates must not cross the still
water line).
48
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Tubular element with uniform circular cross-section and constant wall thickness.
Forces on this element are calculated using Morisons equation.
Open TUBE elements are considered to have open ends. The water surface inside the TUBE is at the
same level as the surface outside. The transverse added mass is based on the outside diameter, but the
axial added mass only uses the cross-sectional area of the TUBE material.
TUBE elements have an axial drag coefficient of 0.016 that cannot be changed.
Slender tube element. The STUB element differs from the TUBE element in the following respects:
STUB elements permit tubes of non-circular cross-section to be modeled, by allowing the tube properties
(diameter, drag coefficient, and added mass coefficient) to be specified in two directions at right angles.
Longer lengths of tube can be input, as the program automatically subdivides STUB elements into sections
of shorter length for integration purposes.
STUB elements should only be employed if the mean diameter is small compared to the length.
The local axis of the element runs from Node 1 to Node 2 (origin at Node 1). The perpendicular from
Node 3 to the local X-axis, together with the local X-axis itself, defines the XZ-plane, and hence the Zaxis. The local Y-axis forms a right-handed set.
If any Z parameter for a STUB is omitted or set to zero it will default to the equivalent Y value.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
49
A circular disc element with no thickness and no mass. The DISC element has a drag coefficient and
added mass coefficient in its normal direction. To add mass, the user can define a PMAS element and
attach it to the DISC element.
The DISC element has two nodes; the first defines the center of the element and the second gives the
normal direction of the element.
The force on a DISC element has two components: added mass force and drag force. As DISC elements
have no thickness, the Froude-Krylov force and hydrostatic force are zero. This is different from a TUBE
end on which the Froude-Krylov force and hydrostatic force are non-zero and therefore are calculated
in Aqwa. The application point of the force is at the centroid of the DISC element when fully submerged,
and at the force center computed by the program if partially submerged. The direction of the force on
a DISC is parallel to the normal direction of the DISC.
50
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
node 2
node 3
node 4
Element 1
i=1
11+0(1+0) = 11
71+0(1+0) = 71
72+0(1+0) = 72
12+0(1+0) = 12
Element 2
i=2
11+1(1+0) = 12
71+1(1+0) = 72
72+1(1+0) = 73
12+1(1+0) = 13
Element 3
i=3
11+2(1+0) = 13
71+2(1+0) = 73
72+2(1+0) = 74
12+2(1+0) = 14
etc.
i.e. 30 elements are generated each with 4 node numbers where
Element 1 nodes are 11,71,72,12
Element 2 nodes are 12,72,73,13
Element 3 nodes are 13,73,74,14
etc.
Note
N3 is zero in all cases, which is very common.
Example 2
For the TUBE element, the data
(10)(11,3)(71,3)(2)(1)
will generate the elements as follows:
node 1
node 2
material group
geometry group
Element 1
i=1
11+0(3+0) = 11
71+0(3+0) = 71
2+0(0+0) = 2
1+0(0+0) = 1
Element 2
i=2
11+1(3+0) = 14
71+1(3+0) = 74
2+1(0+0) = 2
1+1(0+0) = 1
Element 3
i=3
11+2(3+0) = 17
71+2(3+0) = 77
2+1(0+0) = 2
1+2(0+0) = 1
etc.
i.e. 10 elements are generated
Element 1: nodes are 11,71; Material Group = 2; Geometry Group = 1
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
51
Note
N2 (as well as N3) is zero for the material/geometric group number increments, which is very
common. This means that all 10 elements will have the same material and geometric properties.
Example 3
For the DISC element, the data
(1)(11)(71)(1)
will generate one DISC element as indicated by the number in the first pair of brackets. The centroid
position of the disc is defined by the first node number (node 11) and the normal direction of the disc
is decided by the vector from the first node (11) to the second node (71). It should be noted that
whether the normal vector defined as from node 11 to 71 or from node 71 to 11 has no effect on the
results. The last number in the DISC data record is the geometry property group number for this DISC
which is to be defined in Geometric Properties - GEOM (Data Category 4) (p. 65) (DISC has no material
properties).
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Note
The Symmetry data record only applies to TPPL and QPPL elements. All other element types
are are unaffected by the introduction of these data records and must be fully described as
the physical geometry dictates.
5 7
11
16
21
- --- -- ---- ----- ----- -------X|XXX|02|HYDI|XXXXX|
1|XXXXX
- --- -- ---- ----- ----- -------| |
|
| |
|
| |
|
| |
|_(3)Structure Number.
| |
| |_(2)Mandatory data record keyword.
|
|_(1)Optional user data category identifier
(3) All structures between the one specified on the HYDI data record and the current one will be considered to be interacting. For example, if an HYDI data record specifying structure 2 is input for structure
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
53
Note
Structures for which interaction effects are to be calculated must be defined with consecutive
structure numbers.
5 7
- --- -- ---- -----X|XXX| |RMXS|XXXXX
- --- -- ---- -----| |
| |
| |_(2)Symmetry Specification
|
|_(1)Optional user data category identifier
These data records are most useful when structures have been created using the mesh generator, which
automatically creates a symmetric structure using the SYMX data record. When the full structure has
been created it can then be put into its correct position using the MSTR data record.
5 7
11
16
21
- --- -- ---- ----- ----- -------X|
|02|MSTR|XXXXX|XXXXX| (NNN) (X,Y,RZ)
- --- -- ---- ----- ----- --------
54
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(4) The new position of the structure is defined by giving a node number on the structure, the position
to which it is to be moved and the yaw angle (degrees). The vessel is yawed about the specified node.
Only X and Y positions and yaw can be given as the Z position is defined by the ZCGE or ZLWL data
records and Aqwa-Line does not permit the model to be rotated in roll or pitch.
Note
All the output from Aqwa-Line (RAOs, hydrodynamic force coefficients etc.) is in the FRA.
This is still so if a structure is rotated using the MSTR data record, so, for example, motions
about the X-axis may no longer correspond to roll of the vessel.
5 7
11
16
21
- --- -- ---- ----- ----- -------X|
|02|FIXD|XXXXX|XXXXX|
- --- -- ---- ----- ----- -------|
| |
|
| |
|
| |_(3)Mandatory data record keyword.
|
|
|
|_(2)Optional user data category identifier.
|
|_(1)Compulsory END on Last Data Record in Data Category(A3)
55
2
5 7
11
16
21
- --- -- ---- ----- ----- ----------------------------|X|
| |VLID|
|
|(DAMP=???,GAP=???)
- --- -- ---- ----- ----- ----------------------------| |
|
|
|
| |
|
|
|_(2)Mandatory parameters
| |
|
|
| |
|
|_(1)Group number
| |
|
| |
|_Compulsory Data Record Keyword (A4)
| |
| |_Optional User Identifier (A2)
|
|_Compulsory END on last data record in Data Category (A3)
(1) This is the group number that will be used for the lid elements. If there is more than one structure
with a lid a different group must be used for each lid. If a group number is not specified a default value
of 999 will be used.
(2) Lid Parameters
DAMP: damping factor for the lid, typically between 0.0 and 0.2; 0 will give no effect, 0.2 will result
in heavy damping of surface elevation at the elements.
GAP: characteristic length for the lid. This will typically be the gap between two adjacent vessels or
the smallest dimension of a moonpool. It is not used to define the size of the lid itself.
56
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
5 7
16
21
- --- -- ---- - ---- ----- ------------------------X|XXX|01|ASYM|X|DIFF|
1|(32)(111)(217)XXXXXXXXXXX
- --- -- ---- - ---- ----- ------------------------| |
|
|
|
|
|
| |
|
|
|
|
|
| |
|
|
|
|
|_(7)Finish user node no.
| |
|
|
|
|
(MAX # nodes total = 100)
| |
|
|
|
|
| |
|
|
|
|_(6)Start user node no.
| |
|
|
|
| |
|
|
|_(5)No. of Axi-symmetric rotations 8,16,32,64
| |
|
|
| |
|
|_(4)Element Group Number
| |
|
| |
|
| |
|_(3)Diffracting - leave blank for non diffracting.
| |
| |_(2)Mandatory data record keyword.
|
|_(1)Optional User Identifier
(3) The above example generates diffracting elements. Leave this field blank for non-diffracting elements
or elements above the water line.
(4) The element group that will be used for all elements generated by this data record. If 0 or blank is
input then the group number will start from 1 and for each set of elements generated along the axis
the group number will increment by 1. i.e. groups 1-32 will be generated if 32 is the number of axisymmetric rotations.
(5) No. of Axi-symmetric rotations to produce a full 360 degree description of the structure. This is restricted to 8,16,32 and 64.
(6/7) The start/finish user node numbers of the line of nodes to be rotated. This defines the 'profile'
line. e.g.:
if the start node is 111
the finish node is 217
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
57
2
5 7
11
16
21
- --- -- ---- ----- ----- ----------------------------|X|
| |ILID| AUTO|
|(LID_SIZE=????,START_NODE=NNNNN)
- --- -- ---- ----- ----- ----------------------------| |
|
|
|
|
| |
|
|
|
|_(3)Optional parameters
| |
|
|
|
| |
|
|
|_(2)Group number
| |
|
|
| |
|
|_(1)Automatic lid generation
| |
|
| |
|_Compulsory Data Record Keyword (A4)
| |
| |_Optional User Identifier (A2)
|
|_Compulsory END on last data record in Data Category (A3)
58
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
User-Defined Lid
2
5 7
11
16
21
- --- -- ---- ----- ----- ----------------------------|X|
| |ILID|
|
|
- --- -- ---- ----- ----- ----------------------------| |
|
|
| |
|
|_(1)Group number
| |
|
| |
|_Compulsory Data Record Keyword (A4)
| |
| |_Optional User Identifier (A2)
|
|_Compulsory END on last data record in Data Category (A3)
(1) The lid elements must be put into a single group with this number. If there is more than one structure
with a lid a different group must be used for each structure. If a group number is not specified a value
of 999 will be assumed.
The user must then define a number of elements that should:
be within the vessel
have their normals pointing upwards (i.e. the nodes must be listed anti-clockwise when looking down
on the elements)
be at the still water surface
Except for the requirement that the centroids should be below the water surface, the usual modeling
checks are still applied, and the TOTAL number of elements must still be fewer than the maximum
permitted.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
59
(4) The position of the waterline is defined in the same axes as used to define the nodes in Data Category
1.
NOTES.
If a ZLWL data record is used, the ZCGE data record in Data Category 7 is no longer required. An error
message will be issued if a ZCGE data record is present and is not consistent with the ZLWL data record.
If the ZCGE and ZLWL data records are both absent, the definition position based on the nodes in
Data Category 1 will be used.
This position is NOT USED as the starting position for an analysis using Aqwa-Librium, Fer, Drift or
Naut. This must be specified in Data Category 15 or by using the RDEP option (see Administration
and Calculation Options for the Aqwa Suite (p. 252)).
60
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Y-direction
41
26
81
51
41
26
(2) The last four parameters give the X and Y limits of the grid over which the wave pressures will be
calculated. These parameters cannot be set, and should be omitted, when Aqwa is run as a batch program
or from the AGS.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
61
62
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) Material group number - This number is referred to in the Element Topology Data Record (Data
Category 2) and is the index to the values given in columns 21-40 (2).
(2) Material property values. Material properties may, optionally, be divided into sets for each structure
(see The STRC Data Record - Coordinate Structure Association (p. 40)). Any element whose material
group number corresponds to that given in columns 16-20 (1) has these values associated with it. Details
of the properties appropriate to each element type are given below.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
63
Property 1
Property 2
Property 3
T/QPPL
None
None
None
TUBE
Density
None
None
STUB
Mass/unit
length
Y axis inertia/unit
length
Z axis inertia/unit
length
PMAS
Mass
None
None
TELM
Density
Young's Modulus
None
DISC
None
None
None
(3) The structure number must start from 1 on the 1st data record and increment by 1 for each structure
data record input. The number of ELM* data categories must correspond to the number of structure
data records.
64
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Continuation Data Record - Only required for Morison elements. It may be omitted if the appropriate
default values are required (see (4) and (5) below)
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
65
(1) The Element Type (always four characters) provides the classification for a particular element. This
must correspond to the element type which references the geometric group number on the Element
Topology Data Record in Data Category 2. Valid codes for element types in this data category are TUBE,
STUB, DISC and PMAS.
Note
QPPL, TPPL and PBOY do not have any geometric properties.
(2) Geometric group number - This number is referred to in the Element Topology Data Record (Data
Category 2) and is the index to the values given in columns 21-70 (2).
(3) Geometric parameters - Any element on any structure whose geometric group number corresponds
to that given in columns 16-20 (2) has up to 6 of these values associated with it. A summary of the
properties appropriate to each element is shown below. (Units are consistent with the physical quantity
defined by each parameter.)
(4)/(5) Viscous drag/added mass coefficient associated with elements having hydrodynamic properties
in the Morison regime. If the continuation data record is omitted, or if either field is left blank, or if zero
value is input, default values are used for the appropriate element type. These values are shown below
in Columns 7 and 8.
Element
Type
TUBE
Diameter
Thickness
Continuation
0.
Sealed
Cut
1
Cut 2
0.75
1.0
Cd Z
Added Mass
Z
1.
Open
STUB
Diameter
Y
Diameter
Z
Area
Cd
Y
Added Mass
Y
PMAS
Ixx
Ixy
Ixz
Iyy
Iyz
DISC
Diameter
1.0
Note
Within the program, the Cd and Ca values for a DISC will always be halved; therefore the
calculated forces are only for the side of the DISC which is subject to the hydrodynamic
66
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(3) The structure number must start from 1 on the 1st data record and increment by 1 for each structure
data record input. The number of ELM* data categories must correspond to the number of structure
data records.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
67
68
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Note: If data records denoted (Optional) are omitted, the program will produce the appropriate default
value shown for each parameter.
(1) Note that the water depth is fundamental to the calculation of all wave properties.
69
70
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
71
(1) The data record keyword defines the units used for the values in columns 21-80, as follows.
FREQ - frequencies defined in radians/sec
HRTZ - frequencies defined in hertz
PERD - periods defined in seconds
(2)-(3) These are the frequency/period numbers associated with the values in columns 21-80. These
numbers must be used in subsequent data categories when referring to the frequency values where
they are used to cross check the data input. The initial number refers to the first value input (columns
21-30) and the terminal number refers to the last value input. If the initial number is zero or blank the
program will generate frequencies automatically - see below.
------------------------------------------------| If the Frequency/Period
| Initial |Terminal |
|Values Correspond to Numbers | Number | Number |
|-----------------------------+---------+---------|
|
1,2,3,4,5 and 6
|
1
|
6
|
|-----------------------------+---------+---------|
|
4,5,6, and 7
|
4
|
7
|
|-----------------------------+---------+---------|
|
3 only
|
3
|
3
|
|
or |
3
|
|
-------------------------------------------------
72
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Frequency Generation
Automatic generation of frequencies is specified by setting the 1st frequency number to 0 (or blank).
Aqwa will then generate frequencies numbered from the previous frequency to the final frequency
number on the data record. The generated frequencies will be equally spaced between the first and
last frequencies on the data record.
This data record will generate 20 frequencies numbered from 1 to 20, equally spaced between 0.1 and
1.8 hertz.
- --- -- ---- ----- ----- --------- --------| |
| |HRTZ|
0|
20|
0.1|
1.8|
- --- -- ---- ----- ----- --------- ---------
These data records will generate 5 frequencies numbered from 1 to 5, equally spaced between 0.1 and
0.5 rad/s, followed by 10 frequencies numbered from 6 to 15, equally spaced between 0.55 and 1.0
rad/s.
- --- -- ---- ----- ----- --------- --------| |
| |FREQ|
0|
5|
0.10|
0.50|
| |
| |FREQ|
0|
15|
0.55|
1.00|
- --- -- ---- ----- ----- --------- ---------
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
73
(1)-(2) These are the direction numbers associated with the values in columns 21-80. These numbers
must be used in subsequent data categories when referring to the direction values where they are used
to cross check the data input. The first number refers to the first value input (columns 21-30) and the
second number refers to the last value input. If the first number is zero or blank the program will generate directions automatically - see below.
(3)-(8) These are the values of direction at which the hydrodynamic parameters (Data Categories 7 and
8) are to be calculated. They are also the default directions at which the hull drag coefficients (Data
Category 10) are to be input. These values are associated with the numbers defined in columns 11-20
(see parameters (3)-(8) in The FREQ/PERD/HRTZ Data Record - Frequencies/Periods at which the Hydrodynamic Parameters are Computed (p. 72)).
74
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
EXAMPLES
These two data records define 9 directions, numbered from 1 to 9, with values from -180 to 0.
- --- -- ---- ----- ----- --------- --------- --------- --------- ---------- ---------| |
| |DIRN|
1|
5|
-180.0|
-157.5|
-135.0|
-112.5|
-90.0|
|
| |
| |DIRN|
6|
9|
-67.5|
-45.0|
-22.5|
0.0|
|
|
- --- -- ---- ----- ----- --------- --------- --------- --------- ---------- ----------
Direction Generation
Automatic generation of directions is specified by setting the 1st direction number to 0 (or blank). Aqwa
will then generate directions numbered from the previous direction to the final direction number on
the data record. The generated directions will be equally spaced between the first and last directions
on the data record.
This data record will generate 21 directions numbered from 1 to 21, equally spaced between 0.0 and
180.
- --- -- ---- ----- ----- --------- --------| |
| |DIRN|
0|
21|
0.0|
180.0|
- --- -- ---- ----- ----- --------- ---------
These data records will generate 5 directions numbered from 1 to 5, equally spaced between 0 and
40, followed bya single direction at 45, followed by 5 more directions numbered from 7 to 11, equally
spaced between 50 and 90.
- --- -- ---- ----- ----- --------- --------| |
| |DIRN|
0|
5|
0.0|
40.0|
| |
| |DIRN|
6|
6|
45.0|
|
| |
| |DIRN|
0|
11|
50.0|
90.0|
- --- -- ---- ----- ----- --------- ---------
75
(1) The value of the frequency and frequency-dependent parameters associated with this number are
moved from their existing position to the position specified in columns 16-20, leaving a position in the
range of frequencies whose value of frequency is undefined. An additional frequency may subsequently
be input using a FREQ data record to specify the value of the frequency at the position corresponding
to this frequency number.
(2) This is the number of the position to which the value of frequency and frequency- dependent
parameters (associated with the number specified in columns 11-15) are to be moved.
76
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
The parameters indicated as 'UN' (undefined) are then calculated (Aqwa-Line only) or expected as input
in Data Categories 7 and 8. Note that if a frequency of 1.0 was required to be input, a FREQ data record
with frequency number 5 and this frequency value is valid WITHOUT a preceding MVEF data record.
Note
The DELF data record relates only to the structure indicated by the data category keyword.
2
5 7
11
16
- --- -- ---- ----- ----|X|
| |DELF|XXXXX|
|
- --- -- ---- ----- ----| |
|
|
| |
|
|
| |
|
|_Number of Frequency and Associated
| |
|
Parameters which are to be Deleted(I5)
| |
|
| |
|_Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on Last Data Record in Data Category(A3)
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
77
Note
Redefinition of a different range of directions will cause an ERROR, if Data Category 10 is
used to input direction-dependent parameters based on the ORIGINAL range of directions
which became undefined. A warning message is therefore output if the only existing frequency
is deleted.
See Database Import (p. 82) for examples showing the import of data from an existing Aqwa-Line
database.
78
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) This may be any of the numbers assigned to Aqwa scratch files. See Section 9 entitled Running the
Program in each of the Aqwa program manuals, for the FORTRAN units which are valid on computer
installations.
(2) This is a file name, usually a *.HYD file from a previous Aqwa-Line analysis, whose database is to be
copied.
Note
It is preferable to use option (2) in conjunction with the CSTR data record (The CSTR Data
Record - Copying Existing Hydrodynamic Parameters for a Specific Structure Number (p. 78))
and the CPDB data record (The CPDB Data Record - Copy Existing Hydrodynamic Data
Base (p. 81)). In this case, option (1), the file unit number, is not required.
(1) The value of the frequency-dependent parameters associated with this number within the structure
and backing file specified on the preceding CSTR and FILE data record, are to be copied to the position
specified in columns 16-20, within the structure indicated by the data category keyword. These parameters are referred to as those within the source file. Note that if there is no preceding CSTR or FILE
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
79
80
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
10.11. The CPDB Data Record - Copy Existing Hydrodynamic Data Base
This data record is used when copying the hydrodynamic data from a backing file (*.HYD) from a previous
run in order to duplicate the data base for the structure indicated by the data category keyword. The
preceding CSTR and/or FILE data records define the structure and the source file from which the data
base is to be copied. See also The CPYF/CPYP/CPYH Data Records - Copy Frequency Parameters (p. 79),
Note 2.
See Database Import (p. 82) for examples showing the import of data from an existing Aqwa-Line
database.
2
5 7
- --- -- ---|X|
| |CPDB|
- --- -- ---| |
|
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
81
(1) Speed must be defined in units consistent with those used for gravity and density in Data Category
5.
AL_RUN1.HYD
END
06
FDR2
06FILE
06CSTR
1
06CPDB
AL_RUN3.HYD
END
82
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Database Import
06
FDR1
06FILE
06CSTR
1
06CPDB
AL_RUN1.HYD
END
06
FDR2
06FILE
06CSTR
1
06CPDB
AL_RUN2.HYD
END
06
FDR5
06FILE
06CSTR
1
06CPDB
AL_RUN3.HYD
END
06
FINI
Note
A FINI data record is required. This is because the program is expecting to read 5 FDR* data
records and there are only three.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
83
84
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
85
Wave Frequency Dependent Parameters and Stiffness Matrix - WFS* (Data Category
7)
|
|
|
|_Compulsory Data Category Keyword(A4)
|
|_Optional User Identifier(A2)
Note
LSTF data records must be used to define the whole hydrostatic stiffness matrix for a particular
structure. If any LSTF data records are used, undefined rows will be re-set to zero.
The Z coordinate of the center of gravity and the buoyancy force at equilibrium must always
be defined when using LSTF data records to specify the stiffness matrix.
See the notes at the end of this section about the effects of the LSTF data record on the calculations.
2
5 7
11
16
21
31
41
- --- -- ---- ----- ----- --------- --------- -----|X|
| |LSTF|XXXXX|
|
|
|
...
- --- -- ---- ----- ----- --------- --------- -----| |
|
|
|
|
|
| |
|
|
|---------------------| |
|
|
|
| |
|
|
|_(2)-(7)6 Stiffness Values
| |
|
|
(6E10.0)(Units for freedoms
| |
|
|
1-3 = force/length,
| |
|
|
4-6 = force*length/RADIAN)
| |
|
|
| |
|
|_(1)Row Number of Stiffness Matrix(I5)
| |
|
| |
|_Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on Last data record in data category(A3)
(1) This number indicates which row of the stiffness matrix the values input in columns 21-80 relate to.
(2)-(7) These are the values which replace the row (1) in the hydrostatic stiffness matrix.
Note
Stiffness related forces acting on the structure
The linear hydrostatic stiffness matrix relates to the hydrostatic forces contributing to the
equations of static equilibrium of a structure. Specifically, the net linear hydrostatic forces
F(s), acting at the center of gravity of a structure, when the structure is at an arbitrary
position X, is given by:
F(s) = K ( X(e) - X ) + B(e)
86
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
87
Wave Frequency Dependent Parameters and Stiffness Matrix - WFS* (Data Category
7)
|_Compulsory END on Last data record in data category(A3)
Note
1. An alternative method of specifying the free-floating equilibrium position is to use the ZLWL
data record in Data Category 2 (The ZLWL Data Record - Waterline Height (p. 59)).
2. If the ZCGE and ZLWL data records are both absent, the definition position based on the
nodes in Data Category 1 will be used.
3. This position is not USED as the starting position for an analysis using Aqwa-Librium, Fer, Drift
or Naut. This must be specified in Data Category 15 or by using the RDEP option (Administration
and Calculation Options for the Aqwa Suite (p. 252)).
88
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
The WAMS Data Record - Wave Frequency Added Mass Matrix and the WDMP Data
Record - Wave Frequency Damping Matrix
| |
|
|
(HRTZ=hertz)
| |
|
|
| |
|
|_(1)Frequency/Period Number(I5)
| |
|
| |
|
| |
|
| |
|_Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on last data record in data category(A3)
(1) This is the number of the frequency/period which corresponds to the value in columns 21-30. This
number must correspond to one of the frequency numbers specified in Data Category 6.
(2) The frequency/period value should be the same as that specified in Data Category 6 and is used
only as a check that the frequency number(1) has been input correctly. It does not redefine the frequency
value. If the value is not the same as in Data Category 6, an error will occur.
11.7. The WAMS Data Record - Wave Frequency Added Mass Matrix and
the WDMP Data Record - Wave Frequency Damping Matrix
These data records may be used to input the added mass matrix (WAMS) and linear damping matrix
(WDMP) at the frequency/period specified on the preceding FREQ/PERD/HRTZ data record. If the matrix
has been read from backing file, the WAMS or WDMP data record will replace (i.e. overwrite) the existing
values within that matrix. (NB see Note 1 at the end of this section, when using Aqwa-Line.)
2
5 7
11
16
21
31
41
- --- -- ---- ----- ----- --------- --------- -----|X|
| |WAMS|
|
|
|
|
...
- --- -- ---- ----- ----- --------- --------- -----|X|
| |WDMP|
|
|
|
|
...
- --- -- ---- ----- ----- --------- --------- -----| |
|
|
|
|
|
|
| |
|
|
|
|---------------------| |
|
|
|
|
| |
|
|
|
|_(3)-(8) 6 Mass/Inertia or Linear Damping Values(6E10.0)
| |
|
|
|
| |
|
|
|_(2) Row Number of Added Mass or Damping Matrix(I5)
| |
|
|
| |
|
|_(9),(10) Structure Number(I5)
| |
|
| |
|_(1) Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on last data record in data category(A3)
(1) This data record keyword indicates whether the values (3) and row number (2) refer to the added
mass (WAMS) or damping matrix (WDMP).
(2) This number indicates to which row of the added mass/damping matrix the values in columns 2180 correspond.
(3)-(8) These are the values which replace those in the row (2) of the existing added mass/damping
matrix.
(9) If the values input in columns 21-80 are to be used for convolution, the structure number must be
input here.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
89
Wave Frequency Dependent Parameters and Stiffness Matrix - WFS* (Data Category
7)
(10) If the analysis includes hydrodynamically interacting structures, this is the number of the structure
which interacts with that defined on the data category keyword (WFS*). The structure number must be
input even if the co-efficients are those giving the effect of the structure on itself (i.e. structure 1 in
data category WFS1). Both complete added mass and damping matrices should be input; if part matrices
are input the results will probably be inaccurate or the program may fail. When using Aqwa-Line, see
note 1 below.
Note
1. Effect of the WAMS/WDMP Data Record when Running Aqwa-Line
Input of the added mass/damping matrix is redundant information when running AqwaLine, as the program normally calculates these matrices. Therefore, use of one or more
WAMS/WDMP data records will instruct Aqwa-Line not to perform a radiation/diffraction
analysis at this frequency. The user must then specify all other parameters which would
otherwise be calculated at this frequency, i.e. the linear damping/added mass matrix,
Froude Krylov and diffraction forces.
2. Convolution
When the added mass and damping have been input by the user, it is not possible to use
convolution to calculate radiation forces in Aqwa-Drift or Naut.
11.8. The WDGA Data Record - Wave Frequency Diagonal Added Mass
and the WDGD Data Record - Wave Frequency Diagonal Damping
These data records represent a more convenient manner to input the added mass matrix (WDGA) and
linear damping matrix (WDGD) at the frequency/period specified on the preceding FREQ/PERD/HRTZ
data record, when there are no off diagonal terms coupling the degrees of freedom of motion. They
are used instead of the data records WAMS and WDMP respectively, which are used only when these
coupling terms exist. If the added mass or damping matrix has been read from backing file, the WDGA
or WDGD data record will replace (i.e. overwrite) the existing values within that matrix. See also the
note at the end of this section, when using Aqwa-Line.
2
5 7
11
21
31
41
- --- -- ---- --------- --------- --------- -----|X|
| |WDGA|XXXXXXXXX|
|
|
...
- --- -- ---- --------- --------- --------- -----|X|
| |WDGD|XXXXXXXXX|
|
|
...
- --- -- ---- --------- --------- --------- -----| |
|
|
|
|
| |
|
|---------------------| |
|
|
| |
|
|_(2)-(7)6 Mass/Inertia or Linear Damping Values(6E10.0)
| |
|
| |
|
| |
|_(1)Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on last data record in data category(A3)
(1) This data record keyword indicates whether the values (3) refer to the added mass matrix (WDGA)
or to the damping matrix (WDGD)
90
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Note
Effect of the WDGA/WDGD Data Record when Using Aqwa-Line
Input of the wave frequency added mass/damping matrix is redundant information when
running Aqwa-Line as the program normally calculates these matrices. Therefore, use of a
WDGA/WDGD data record will instruct Aqwa-Line not to perform a radiation/diffraction
analysis at this frequency. The user must then specify all other parameters which would
otherwise be calculated at this frequency, i.e. the linear damping/added mass matrix, Froude
Krylov and diffraction forces.
Note
If the convolution method is to be used in a time domain analysis, the AAMS and ADMP data
records which define the frequency dependent additional added mass and damping matrices
should be replaced by the frequency independent additional added mass and damping
defined by FIAM, FIDP (for matrices) ( The FIAM/FIDP Data Record - Frequency Independent
Additional Added Mass/Damping Matrices (p. 92)) or FIDA, FIDD (for diagonal terms) (The
FIDA/FIDD Data Record - Frequency Independent Additional Diagonal Added Mass/Damping (p. 92)).
2
5 7
11
16
21
31
41
- --- -- ---- ----- ----- --------- --------- -----|X|
| |AAMS|XXXXX|
|
|
|
...
- --- -- ---- ----- ----- --------- --------- -----|X|
| |ADMP|XXXXX|
|
|
|
...
- --- -- ---- ----- ----- --------- --------- -----| |
|
|
|
|
|
| |
|
|
|---------------------| |
|
|
|
| |
|
|
|_(3)-(8)6 Added Mass/Damping Values(6F10.0)
| |
|
|
| |
|
|_(2)Row Number of Added Mass/Damping Matrix(I5)
| |
|
| |
|_(1)Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on last data record in data category(A3)
(1) This data record keyword indicates whether the values (3) and row number (2) refer to the added
mass (AAMS) or damping (ADMP) matrix
(2) This number indicates the row of the added mass/damping matrix to which the values in columns
21-80 correspond.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
91
Wave Frequency Dependent Parameters and Stiffness Matrix - WFS* (Data Category
7)
(3)-(8) These are the values which add TO those in the row (2) of the existing added mass/ damping
matrix or add TO any subsequent values input or calculated by Aqwa-Line.
(1) This data record keyword indicates whether the values refer to the added mass (FIAM) or damping
(FIDP)
(2) This number indicates the row of the added mass/damping matrix to which the values in columns
21-80 correspond.
(3)-(8) These are the values which add TO those in the existing added mass/ damping matrix or add TO
any subsequent added mass/damping values input or calculated by Aqwa-Line.
92
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
The DIRN Data Record - Directions at which the Parameters are Defined
- --- -- ---- ----- ----- --------- --------- -----| |
|
|
|
|
| |
|
|---------------------| |
|
|
| |
|
|_(2)-(7)6 Added Mass/Damping Values(6F10.0)
| |
|
| |
|
| |
|
| |
|_(1)Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on last data record in data category(A3)
(1) This data record keyword indicates whether the values refer to the added mass (FIDA) or damping
(FIDD).
(2)-(7) These are the values which add TO the diagonal terms of the existing added mass/ damping
matrix or add TO any subsequent diagonal added mass/damping values input or calculated by AqwaLine.
Note
Only one FIDA data record is needed for defining the diagonal terms of an additional added
mass matrix. If more than one FIDA data record is used for a structure, only those values in
the last FIDA data record will be used. This also applies to the FIDD data record.
11.12. The DIRN Data Record - Directions at which the Parameters are
Defined
This data record is used to specify the direction to which the following diffraction forces, Froude Krylov
forces or response amplitude operators, relate. The direction specified is operative until another DIRN
data record is input. (NB See also the note at the end of this section when using Aqwa-Line.
2
5 7
11
16
21
- --- -- ---- ----- ----- --------|X|
| |DIRN|XXXXX|
|
|
- --- -- ---- ----- ----- --------| |
|
|
|
| |
|
|
|
| |
|
|
|
| |
|
|
|_(2)Values of Direction (DEGREES)(6F10.0)
| |
|
|
| |
|
|_(1)Direction Number(I5)
| |
|
| |
|_Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on last data record in data category(A3)
(1) This is the number of the direction which corresponds to the values input in columns 21-30. This
number must be one of the direction numbers specified in Data Category 6.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
93
Wave Frequency Dependent Parameters and Stiffness Matrix - WFS* (Data Category
7)
(2) The direction value should be the same as that specified in Data Category 6 and is used only as a
check that the direction number (1) has been input correctly. It does not redefine the direction value.
If the value is not the same, an error will occur.
Note
Effect of the DIRN Data Record when Running Aqwa-Line
Input of the Froude Krylov forces, diffraction forces or RAOs is redundant, when running
Aqwa-Line, as the program normally calculates these parameters. The DIRN data record is a
precursor to inputting values of these parameters. Therefore, use of a DIRN data record will
instruct Aqwa-Line not to perform a radiation/diffraction analysis at this frequency. The user
must then specify all the other parameters which would otherwise have been calculated by
the analysis at this frequency, i.e. the linear damping/added mass matrix, Froude Krylov and
diffraction forces.
11.13. The TDIF/RDIF/TFKV/RFKV/TRAO/RRAO Data Records - Wave Frequency Diffraction Forces and Responses
These data records are used to input the vectors of diffraction forces, Froude Krylov forces and RAOs
at the frequency/period specified on the preceding FREQ/PERD/HRTZ data record, and at the direction
specified on the preceding DIRN data record. If the vectors have been read from backing file, the values
input on each data record will replace (i.e. overwrite) the existing values.
2
5 7
11
21
31
41
51
61
71
- --- -- ---- --------- ----------- ----------- ----------- ---------- ----------- ----------|X|
| |TDIF|XXXXXXXXX|
|
|
|
|
|
|
- --- -- ---- --------- ----------- ----------- ----------- ---------- ----------- ----------|X|
| |RDIF|XXXXXXXXX|
|
|
|
|
|
|
- --- -- ---- --------- ----------- ----------- ----------- ---------- ----------- ----------|X|
| |TFKV|XXXXXXXXX|
|
|
|
|
|
|
- --- -- ---- --------- ----------- ----------- ----------- ---------- ----------- ----------|X|
| |RFKV|XXXXXXXXX|
|
|
|
|
|
|
- --- -- ---- --------- ----------- ----------- ----------- ---------- ----------- ----------|X|
| |TRAO|XXXXXXXXX|
|
|
|
|
|
|
- --- -- ---- --------- ----------- ----------- ----------- ---------- ----------- ----------|X|
| |RRAO|XXXXXXXXX|
|
|
|
|
|
|
- --- -- ---- --------- ----------- ----------- ----------- ----------- ----------- ----------| |
|
|
|
|
|
|
|
| |
|
(2)Surge(X) (3)Surge(X) (4)Sway(Y)
(5)Sway(Y) (6)Heave(Z) (7)Heave(Z)
| |
|
Roll(RX)
Roll(RX)
Pitch(RY)
Pitch(RY)
Yaw(RZ)
Yaw(RZ)
| |
|
Amplitude
Phase
Amplitude
Phase
Amplitude
Phase
| |
|
| |
|
| |
|_(1)Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on last data record in data category(A3)
(1) This data record keyword indicates whether values in columns 21-80 are diffraction forces, Froude
Krylov forces or Response Amplitude Operators (RAOs) and whether the values relate to the translational
or rotational degrees of freedom, i.e.:
TDIF - Translation Diffraction Forces
RDIF - Rotational Diffraction Forces
TFKV - Translational Froude Krylov Forces
94
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) Aqwa firstly calculates the hydrostatic stiffness matrix based only on the cut water plane and displaced
volume properties. It then adjusts the second moments of area IXX, IYY and recalculates its associated
properties, PHI (principal axis), GMX/GMY, BMX/BMY etc. to give the required GM values. The associated
additional hydrostatic stiffness is calculated automatically and stored in the hydrodynamic database.
Note that if the GM value input is less than that based on the geometry alone, the resulting additional
stiffness will be negative. This would be the case if ballast tanks were being modeled, making the
structure less stable, statically.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
95
Wave Frequency Dependent Parameters and Stiffness Matrix - WFS* (Data Category
7)
- --- -- ---X|
| |SSTR|
- --- -- ---| |
|
| |
|_(1)Submerged structure(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on Last
Data Record in Data Category.(A3)
Note
If the analysis includes stage 2 any additional stiffness already in the database will be re-set
to zero and must be re-defined.
See the notes at the end of this section about the effects of this data record when using different programs.
2
5 7
11
16
21
31
41
- --- -- ---- ----- ----- --------- --------- -----|X|
| |ASTF|XXXXX|
|
|
|
...
- --- -- ---- ----- ----- --------- --------- -----| |
|
|
|
|
|
| |
|
|
|---------------------| |
|
|
|
| |
|
|
|_(2)-(7)6 Stiffness Values
| |
|
|
(6E10.0)(Units for freedoms
| |
|
|
1-3 = force/length,
| |
|
|
4-6 = force*length/RADIAN)
| |
|
|
| |
|
|_(1)Row Number of Stiffness Matrix(I5)
| |
|
| |
|_Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on Last data record in data category(A3)
(1) This number indicates which row of the stiffness matrix the values input in columns 21-80 relate to.
(2)-(7) These are the values which add TO the row (1) in the hydrostatic stiffness matrix.
Stiffness related forces acting on the structure
The linear stiffness matrix relates to the hydrostatic forces contributing to the equations of static
equilibrium of a structure. Specifically, the net linear hydrostatic forces F(s), acting at the center of
gravity of a structure, when the structure is at an arbitrary position X, is given by:
F(s) = K ( X(e) - X ) + B(e)
where
K = stiffness matrix
X(e) = equilibrium position
96
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
97
Wave Frequency Dependent Parameters and Stiffness Matrix - WFS* (Data Category
7)
| |
|
|
| |
|
|_(2)2nd Structure Number (I5)
| |
|
| |
|_(1)Compulsory Data Record Keyword(A4)
| |
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on last data record in data category(A3)
(2) The data on this data record defines one row of a stiffness matrix coupling this structure to the one
indicated on the WFS* data record in the data category keyword (referred to here as Structure#1). If
this is zero or omitted, Structure#2 is set to be the same as Structure#1. Otherwise this number must
not be less than Structure#1.
(3) This number indicates which row of the stiffness matrix the values input in columns 21-80 relate to.
(4)-(9) These are the values for the row (3) in the 6x6 additional structure stiffness sub-matrix. The units
for angles are radians.
98
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
99
100
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
101
102
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) The data record keyword indicates the freedom of the structure, as defined in the Fixed Reference
Axis system (FRA), to which the drift coefficients apply, where:
DRFX - X translational freedom of motion
DRFY - Y translational freedom of motion
DRFZ - Z translational freedom of motion
DRRX - X rotational freedom of motion
DRRY - Y rotational freedom of motion
DRRZ - Z rotational freedom of motion, i.e. XY plane rotation
(2)-(3) These are the directions, as defined in Data Category 6, to which the values of the drift coefficients
in columns 21-80 apply, e.g.
----------------------------------------------------|If the Drift Coefficient(s) | Initial | Terminal |
|Correspond(s) to Direction(s)| Direction | Direction |
|-----------------------------+-----------+-----------|
|
1,2,3,4,5 and 6
|
1
|
6
|
|-----------------------------+-----------+-----------|
|
4,5,6, and 7
|
4
|
7
|
|-----------------------------+-----------+-----------|
|
3
|
3
|
3
|
-----------------------------------------------------
(4)-(9) These are the values of the drift coefficients which are defined at the directions specified by (2)
and (3) and at the frequency/period specified on the preceding FREQ/PERD/HRTZ data record. Note
that when less than 6 coefficients are input on one data record, the values in the extra fields on the
data record are ignored.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
103
(1) QTFs are copied from this structure in the QTF file to the structure number indicated by the data
category keyword (DRC* data record).
(2) The format of the input file is the same as the format of the .QTF file written by Aqwa when the
AQTF option is used. The .QTF file is ASCII with the first line containing the version number of the file
itself and optional user title. At present all files must start with version 1.0. A simple example is shown
here and the format is explained in detail below.
DATA FILE
30.00000
45.00000
90.00000
0.4188791
0.5235988
1.5438E-02 9.2001E+05 3.7491E+00 4.0415E+06 2.9768E-01
0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
1.4039E-02 2.9887E+06 -2.8746E+00 -5.9962E+06 -1.4445E-01
1.5872E-02 1.6689E+06 1.0933E+00 3.1541E+07 8.5249E-01
4.9260E-02 -1.6417E+05 4.2471E-01 1.0464E+07 5.7609E-01
2.6780E-02 4.5498E+05 1.8771E+00 1.3861E+07 -3.7482E-01
7.1545E-02 -2.1790E+06 1.2171E+00 -1.3207E+07 -7.2684E-01
1.0648E-02 1.5811E+06 3.7608E-01 9.6265E+06 3.0844E-01
The format descriptions below make use of FORTRAN style codes. The meaning of these is:
nX
n spaces
nIw
nFw
104
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
105
106
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
107
(1) These 6 values specify the Drift Added Mass of the structure indicated by the Data Category Keyword,
at each of the 6 degrees of freedom of motion.
Typical values of drift added mass for a ship structure in the surge, sway direction are 10 to 100 per
cent of the mass of displaced water. The rotational inertia in yaw is typically the same as that of the
displaced water (i.e. the same as a solid mass of the same density and shape as that of the displaced
water).
(1) These 6 values specify the Drift Damping of the structure indicated by the Data Category Keyword,
at each of the 6 degrees of freedom of motion.
108
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Note
In the present version of the program, only the surge, sway (X and Y translation) and yaw
(Z rotation) need to be specified. This also means that there are no off-diagonal terms in the
damping matrix.
13.5. The LFAD Data Record - Use Lowest Wave Frequency Added Mass
and Damping for Drift Frequency Analysis
The default low frequency added mass is the added mass calculated in Aqwa-Line at the lowest wave
frequency.
The default low frequency damping is zero.
2
5 7
11
- --- -- ---|X|
| |LFAD|
- --- -- ---| |
|
| |
|
| |
|
| |
|
| |
|
| |
|_Optional Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on last data record in data categoy(A3)
Note
When this data record is on, the DGAM and DGDP data records should be omitted, and the
program will use the added mass and damping calculated in Aqwa-Line at lowest wave frequency
for drift frequency analysis.
When the convolution method is used, i.e. CONV option is on, in a time domain analysis, Aqwa
automatically calculates the low frequency (asymptotic) added mass and this data record is
ignored. Any additional added mass or damping should be frequency-independent, specified
using the FIDA, FIDD, FIAM or FIDP data records ( The FIDA/FIDD Data Record - Frequency Independent Additional Diagonal Added Mass/Damping (p. 92) and The FIAM/FIDP Data Record
- Frequency Independent Additional Added Mass/Damping Matrices (p. 92)).
109
(1)-(2) These node numbers specify the positions defining the 'length' on the geometric center of the
structure along which the drag force is integrated.
For ship structures these nodes should be positioned where the longitudinal centerline cuts the waterplane at the bow and stern.
For a structure whose plan section (viewed along the negative Z axis) is more square than a ship, the
positions of the nodes (1) and (2) should be at the ends of a horizontal diagonal across the structure.
Note
When the center of gravity is vertically above or below the geometric center of the structure,
the forces in the surge and sway directions are zero.
(3) The yaw-rate drag coefficient is the force per unit length per unit velocity squared where the velocity
is that of a side-on current (i.e. a current at right angles to the direction along which its 'length' is
defined). A function of this force, together with the local current velocity is integrated along the 'length'
of the structure to give the total drag force.
The yaw-rate drag coefficient will depend on the cross sectional shape of the structure which generally
varies along the 'length'. As the coefficient is assumed constant it therefore represents an average value.
Note
For structures where the cross sectional area varies greatly along its length, the choice of
node positions and the yaw-rate drag coefficient becomes more complicated and reference
should be made to the literature.
110
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) This data record keyword indicates whether the values refer to the added mass (FIDA) or damping
(FIDD).
(2)-(7) These are the diagonal values in the 6 X 6 wave drift damping matrix.
(1) This data record keyword indicates whether the values refer to the added mass (FIAM) or damping
(FIDP).
(2) This number indicates the row of the added mass/damping matrix to which the values in columns
21-80 correspond.
(3)-(8) These are the values which add to those in the existing frequency independent added mass or
damping matrix.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
111
112
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 14: Hull Drag Coefficients and Thruster Forces - HLD* (Data
Category 10)
When entering HLD data categories, the * indicates the structure number; for example, enter HLD1 for
Structure 1, HLD2 for Structure 2, .... HL10 for Structure 10.
113
Hull Drag Coefficients and Thruster Forces - HLD* (Data Category 10)
|X|
| |CUFX|
|
|
|
|
|
...
- --- -- ---- ----- ----- --------- --------- --------- --------|X|
| |CUFY|
|
|
|
|
|
...
- --- -- ---- ----- ----- --------- --------- --------- --------|X|
| |CURZ|
|
|
|
|
|
...
- --- -- ---- ----- ----- --------- --------- --------- --------|X|
| |WIFX|
|
|
|
|
|
...
- --- -- ---- ----- ----- --------- --------- --------- --------|X|
| |WIFY|
|
|
|
|
|
...
- --- -- ---- ----- ----- --------- --------- --------- --------|X|
| |WIRZ|
|
|
|
|
|
...
- --- -- ---- ----- ----- --------- --------- --------- --------| |
|
|
|
|
|
|
|
| |
|
|
|
|
|
|
|
| |
|
|
|
|------------------------------------| |
|
|
|
|
| |
|
|
|
|_(4)-(9)Up to 6 Force Coefficients(6E10.0)
| |
|
|
|
| |
|
|
|_(3)Terminal Direction Number(I5)
| |
|
|
| |
|
|_(2)Initial Direction Number(I5)
| |
|
| |
|_(1)Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on last data record in data category(A3)
(1) The data record keyword indicates the freedom of the structure, as defined in the Fixed Reference
Axis system (FRA), to which the force coefficients apply where
CUFX/WIFX - X translational freedom of motion (SURGE)
CUFY/WIFY - Y translational freedom of motion (SWAY)
CURZ/WIRZ - Z rotational freedom of motion (YAW), i.e. rotation in the XY plane
Similar conventions apply to forces or moments in heave (CUFZ/WIFZ), roll (CURX,WIRX), and pitch
(CURY,WIRY). CUFZ may be utilized to simulate the sinkage effect that can be induced due to forward
speed or high current. The RX and RY parameters can be included to account for the distance between
the physical center of load application and the position of the center of gravity where the actual load
will be applied in the numerical model.
(2)-(3) These are the directions, as defined using the DIRN data records (either within this data category
or in the Frequencies and Directions data category), to which the values of the force coefficients in
columns 21-80 apply, e.g.
----------------------------------------------------|If the Force Coefficient(s) | Initial | Terminal |
|Correspond(s) to Direction(s)| Direction | Direction |
|-----------------------------+-----------+-----------|
|
1,2,3,4,5 and 6
|
1
|
6
|
|-----------------------------+-----------+-----------|
|
4,5,6, and 7
|
4
|
7
|
|-----------------------------+-----------+-----------|
|
3
|
3
|
3
|
-----------------------------------------------------
114
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Note
The current and wind force coefficients are defined as the force or moment per unit velocity
squared. The moment is about the center of gravity of the structure.
These forces are a function of the relative velocity between the structure and current/wind.
This means that the current/wind coefficient should still be input even when there is no
current/wind present, as the relative velocity is generally non-zero for a dynamic analysis,
even when there is no current/wind present.
(1) This node number and its corresponding position (defined in Data Category 1) specifies the position
on the structure at which the Thruster Force (2)- (4) acts.
(2)-(4) These are the three components of force in the X,Y, and Z directions (see note below) exerted
by the thruster at the position of the node (1) specified.
Note
The magnitude of the thruster forces is assumed to be constant throughout the analysis.
The three components of force define the direction relative to the structure which is also
assumed constant. This means that their direction relative to the Fixed Reference Axis system
(FRA) will change with the structure position; for example, if the thruster force is defined
with a component in the X direction only (i.e. zero for parameters (3) and (4)) and the
structure at some stage of the analysis is yawed (Z rotation) by 90 degrees, then the direction
of the force will act in the Y direction of the Fixed Reference Axes (FRA).
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
115
Hull Drag Coefficients and Thruster Forces - HLD* (Data Category 10)
2
5 7
11
16
21
31
41
- --- -- ---- ----- ----- --------- --------- --------|X|
| |NLRD|
XXXXXXXXXXXXXX
- --- -- ---- ----- ----- --------- --------- --------|X|
| |BOFF|
|
|
|
- --- -- ---- ----- ----- --------- --------- --------|X|
| |BASE|
|
|XXXXXXXXX|
- --- -- ---- ----- ----- --------- --------- --------| |
|
|
|
|
| |
|
|
|
|
| |
|
|
|
|_(3)Bilge Radius (only for BOFF data record) (F10.0)
| |
|
|
|
| |
|
|
|_(2)Node Number Two(I5)
| |
|
|
| |
|
|
| |
|
|_(1)Node Number One(I5)
| |
|
| |
|
| |
|_Compulsory Data Record Keyword(A4)
| |
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on last data record in data category(A3)
Note
No parameters are required for the NLRD data record
(1)-(2) On the BOFF data record, these two nodes define the lateral and vertical position of the two
bilges. They should be at opposite ends of a line that is perpendicular to the line defined by the nodes
on the BASE data record.
On the BASE data record, these two nodes define the longitudinal extent of the bilges (two bilges are
assumed to have the same longitudinal extent). They should normally be on the centerline of the vessel.
Z is unimportant but is normally near the baseline.
(3) This value is the bilge radius (assuming the two bilges have the same radius), and is only for BOFF
data record.
116
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
_
_
|Vx.mod(Vx)|
|Vy.mod(Vy)|
|Vz.mod(Vz)|
|Rx.mod(Rx)|
|Ry.mod(Ry)|
|Rz.mod(Rz)|
(1) This data record keyword indicates that the drag force on this structure is to be calculated as Morison
hull drag using the coefficients specified.
(2)-(3) Row number and column number of the coefficients. These entries as interpreted as follows:
(a) Row and column number zero or omitted: The 6 values are assigned to the lead diagonal.
(b) Row only specified: The 6 values are assigned to the specified row.
(c) Column only specified: The 6 values are assigned to the specified column.
(d) Row and column specified: The value of C is assigned to the specified row/column. The coefficient
C must be in columns 21 - 30.
The default fluid velocity U is relative steady fluid velocity for translational freedoms, i.e:
U(translational) = U(uniform + current profile) - U(structure)
U(rotational) = -U(structure)
These default values can be changed by using the MDSV data record (see below).
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
117
Hull Drag Coefficients and Thruster Forces - HLD* (Data Category 10)
(1) This data record may (optionally) be used in conjunction with the MDIN data record to change the
value of the relative fluid velocity used in the calculation of the Morison hull drag. It specifies that the
structure velocity only is to be used for translational freedoms. For example,
MDSV
2
5 7
11
16
21
31
41
51
61
- --- -- ---- ----- ----- --------- --------- --------- --------- --|X|
| |IUFC|
|
|
|
|
|
|
...
- --- -- ---- ----- ----- --------- --------- --------- --------- --|X|
| |RUFC|
|
|
|
|
|
|
...
- --- -- ---- ----- ----- --------- --------- --------- --------- --| |
|
|
|
|
|
|
|
| |
|
|
|
|-----------------------------------------| |
|
|
|
|
| |
|
|
|
|_(3)-(8)Up to 6 integer (IUFC)
| |
|
|
|
or real (RUFC) parameters
| |
|
|
|
| |
|
|
|_(2)Last parameter number (I5)
| |
|
|
| |
|
|_(1)1st parameter number(I5)
| |
|
| |
|_Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on last data record in data category (A3)
These data records can be repeated until all the required parameters are defined.
118
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
The DIRN Data Record - Directions at which the Drag Coefficients are defined
14.8. The DIRN Data Record - Directions at which the Drag Coefficients
are defined
If this data record is not present the directions used for the drag coefficients are those defined on the
DIRN data records in Data Category 6.
2
5 7
11
16
21
31
41
- --- -- ---- ----- ----- --------- --------- -----|X|
| |DIRN|
|
|
|
|
...
- --- -- ---- ----- ----- --------- --------- -----| |
|
|
|
|
|
|
| |
|
|
|
|---------------------| |
|
|
|
|
| |
|
|
|
|_(3)-(8)Up to 6 Values of Direction(6F10.0)(DEGREES)
| |
|
|
|
| |
|
|
|
| |
|
|
|_(2)Terminal Direction Number(I5)
| |
|
|
| |
|
|_(1)Initial Direction Number(I5)
| |
|
| |
|
| |
|_Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on last data record in data category (A3)
(1)-(2) These are the direction numbers associated with the values in columns 21-80. The first number
refers to the first value input (columns 21-30) and the second number refers to the last value input.
(3)-(8) These are the values of direction at which the following hull drag coefficients are to be defined.
These values are associated with the numbers defined in columns 11-20. The following rules govern
input of directions:
(a) Direction numbers must be unique.
(b) There must be no gaps in direction numbers; for example, if five directions are input then these
must be numbered 1 through 5.
(c) Directions must be distinct and unique.
(d) If no symmetry is defined, direction values must be input for the directions from -180 to 180.
If SYMX is specified, direction values must be input for the directions from -180 to 0 or from 0 to 180.
If SYMY is specified, direction values must be input for the directions from -90 to +90 or from +90 to
-90.
If both SYMX and SYMY are specified, direction values must be input for one 90 quadrant.
(e) When all data records in this data category have been input, ascending numbers must correspond
to ascending direction values.
(f ) Although more than one direction data record can be used, the total number of directions cannot
exceed 41. 41 directions can be input when there is no symmetry, 21 directions can be input when
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
119
Hull Drag Coefficients and Thruster Forces - HLD* (Data Category 10)
there is one plane of symmetry, and 11 directions when there are two symmetry planes. Symmetry for
hull drag coefficients is defined using the SYMX and SYMY data records in Data Category 10.
(g) The directions defined here for drag coefficients are distinct from those used for the wave directions.
They do not have to be the same in number or value as those in Data Category 6.
Note
This data category cannot specify a greater degree of symmetry than is defined for the
structure in Data Category 2. In other words, SYMX/SYMY cannot be used here if the structure
does not have the corresponding symmetry.
2
5 7
- --- -- ---|X|
| |SYMX|
- --- -- ---|X|
| |SYMY|
- --- -- ---| |
|
| |
|_(1)Symmetry Specification(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on Last
Data Record in Data Category.(A3)
1) The same restrictions on directions apply as for the diffraction analysis directions in Data Category
6.
If no symmetry is defined, coefficients must be input for the directions from -180 to 180.
If SYMX is specified, coefficients must be input for the directions from -180 to 0 or from 0 to -180.
If SYMY is specified, coefficients must be input for the directions from -90 to +90 or from +90 to -90.
If both SYMX and SYMY are specified, coefficients must be input for one 90 quadrant.
120
61
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) This is a node on the structure defined by the HLD* data category (Hull Drag Coefficients and Thruster
Forces - HLD* (Data Category 10) (p. 113)). When calculating the hull drag force on the structure, the
current at this point will be used.
Note
The drag force is still applied at the center of gravity.
The node moves with the structure.
The current used is the total of uniform + profiled current.
If the node is zero or omitted, the current at the center of gravity of the structure will be used.
Both current profile and constant current extend above the water surface, so if the node is
above the water surface there may still be a drag force. See the CPRF data record (CPRF Data
Record - Profiled Current Velocity (p. 124)).
61
(1) This is a vertical position in the FRA. When calculating the hull drag force on the structure defined
by the HLD* data category (Hull Drag Coefficients and Thruster Forces - HLD* (Data Category 10) (p. 113)),
the current at this point will be used.
Note
The drag force is still applied at the center of gravity.
This is a fixed height in the FRA.
The current used is the total of uniform + profiled current.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
121
122
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) Enter CURR or WIND in columns 7-10 as appropriate. Note that both data records may be input if
required.
(2) Current/Wind speed is a scalar quantity and must always be a positive value. The current is uniform
from the sea bed to the water surface and the wind is uniform from the water surface upwards. Their
components in the directions of the Fixed Reference Axes are given by:
X - Component = (Speed)( COS (Direction))
Y - Component = (Speed)( SIN (Direction))
(3) The direction of the current/wind, if left blank or zero is input, is in the positive X direction of the
Fixed Reference axes.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
123
(1) The Z position at which the current is defined is with reference to the Fixed Reference Axes whose
origin is at the water surface. These values will therefore always be negative and must also appear in
order i.e. from the sea bed, at Z = - (Water Depth), up to the water surface at Z = 0.0.
(2) Current/Wind speed is a scalar quantity and must always be a positive value. The current components
at the Z position in columns 11-20 (1) in the directions of the Fixed Reference Axes are given by
X-component = (Current Speed) * cos(Current Direction)
Y-component = (Current Speed) * sin(Current Direction)
(3) This is the direction of the current at the Z position in columns 11- 20 (1). If this field is left blank
or zero is input, it will be the positive X direction of the Fixed Reference axes. As it is common to define
a current profile with this value the same on all data records, the user may wish to input a single CDRN
(CDRN Data Record - Current Direction (p. 124)) instead of repeating the current direction on each data
record.
Note
The maximum number of CPRF data records which may be used to define a current profile is
25.
The current profile remains constant below the lowest Z position or above the highest; it does
not drop to zero outside the defined range.
If a wave height time-history is imported using the IWHT card it may be necessary to define
the current profile with negative velocities. This is because the WHT file must contain a positive
constant current in order for the waves to be calculated correctly. The total current velocity
will be the sum of the value in the WHT file and the values defined here.
124
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) This value of the current direction will apply to all CPRF (current profile) data records where the
direction is given as zero or left blank. It enables the user to avoid repeating the direction of the current
on each or any of the CPRF data records.
Below is an example showing the effect of the CDRN data record.
AS INTERPRETED BY AQWA
Depth
-800.0
-700.0
-600.0
-500.0
-400.0
-300.0
-200.0
-100.0
0.0
Speed
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.7
0.9
Direction
180.0
180.0
180.0
180.0
45.0
45.0
45.0
45.0
45.0
2
5 7
11
21
31
41
- --- -- ---- ---------- ---------- ---------- ---------|X|XXX| |TOWS|XXXXXXXXXX|XXXXXXXXXX|
|
|
- --- -- ---- ---------- ---------- ---------- ---------|
|
|
|
|
|
|
|
|
|
|
|_(2) Tow Direction (E10.0)
|
|
|
|
|
|_(1) Tow Speed (E10.0)
|
|
|
|_Compulsory Data Record Keyword (A4)
|
|_Optional User Identifier (A2)
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
125
126
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
127
2
5 7
11
16
- --- -- ---- ----- ----|X|
| |DACF|
|
|
- --- -- ---- ----- ----| |
|
|
|
| |
|
|
|_(2)Freedom Number(I5)
| |
|
|
| |
|
|_(1)Structure Number(I5)
| |
|
| |
|_Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on last data record in data category(A3)
(1) The structure number must correspond to one of the structures defined in Data Category 2; in other
words, if 1 is input, then this will correspond to the structure defined in Data Category ELM1. If 2 is
input, then this will correspond to the structure defined in Data Category ELM2, and so on.
(2) Freedom numbers must be from 1 through 6; in other words, 1, 2 and 3 correspond to X, Y, and Z
respectively for the translational freedoms and numbers 4, 5 and 6 correspond to X, Y, and Z respectively
for the rotational freedoms.
Note
Axis Systems Associated with the Deactivated Freedoms - if a freedom is deactivated
then no motion of the center of gravity will occur in that freedom, where a freedom is defined
as motion about an axis system parallel to the fixed reference axes, not the local structure
axes.
Example
If translational and rotational motion in the X,Y and Z directions of the structure (as originally defined
in Data Category 2) is surge, sway and heave (translational) and roll, pitch and yaw (rotational) respectively, then the following situations are possible:
(i) If the X translational freedom is deactivated, and if the structure is analyzed in a position where its
local axes are parallel to the Fixed Reference Axes then surge motion will be eliminated and there will
be no restrictions on the sway motion. However, if the structure is yawed by 90 degrees then the local
sway motion will be eliminated and there will be no restrictions on the local surge motion.
128
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
2
5 7
11
16
21
26
31
36
41
46
51
- --- -- ---- ----- ----- ----- ----- ----- ----- ----- ----- ----X|
| |DCON|
|
|
|
|XXXXX|
|
|
|XXXXX|
- --- -- ---- ----- ----- ----- ----- ----- ----- ----- ----- ----| |
|
|
|
|
|
|
|
|
| |
|
|
|
|----|
|-----| |
|
|
|
|
|
|
| |
|
|
|
|
|
|_(7)-(9)Node Numbers(2I5)
| |
|
|
|
|
|
| |
|
|
|
|
|
| |
|
|
|
|
|_(6)Second Structure Number(I5)
| |
|
|
|
|
| |
|
|
|
|
| |
|
|
|
|_(3)-(5)Node Numbers(2I5)
| |
|
|
|
| |
|
|
|_(2)First Structure Number(I5)
| |
|
|
| |
|
|_(1)Number of Locked Rotational Freedoms at the Constraint(I5)
| |
|
| |
|
| |
|_Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on last data record in data category(A3)
Constraint Definition
(1) Four types of constraint are valid in the program which are coded as 0, 1, 2 and 3, representing the
number of rotational freedoms which are locked at the constraint. Therefore, rotational reactions or
moments transmitted through the joint at the constraint, from the first structure to the second structure,
and vice versa, also correspond to these numbers. The four types of joint are as follows:
0
Universal
Hinged
Locked
Transmitting a moment in all three freedoms and not free to rotate at all.
This type of constraint enables the user to find the reactions between two
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
129
Joint
Mandatory Input
For Structure 1
For Structure 2
B/Socket
(3)
(7)
Universal
(3) (4)
(7) (8)
Hinged
(3) (4)
(7) (8)
Locked
(3)
(7)
(6) This is the number of the second structure on which the constraint is defined otherwise as in (2)
except that '0' structure number is illegal. See rules at the end of this section.
(7)-(9) As for the first node. See (3)-(5).
Universal Joint
The axis defined on each structure by the second node will give two axes
which will always be at right angles to each other. The joint will only allow
relative motion about these two axes.
Hinged Joint
The axis defined on each structure by the second node will give two axes
which will always be coincident. The joint will only allow relative motion
about this coincident axis.
Locked Joint
130
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Universal Joint
Hinged Joint
Locked Joint
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
131
Note
132
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1), (2), (3), (4) Coulomb friction always uses a local axis system for calculation with the X-axis being the
axis of the instantaneous relative rotational velocity of the two parts of the articulation. For a hinge this
will always be the axis of the hinge; for a ball joint it will vary as the two structures move relative to
one another. The frictional moment is given by
+ +
+ + +
where:
= 0 if the relative rotational velocity is less than 0.001 rad/s, 1 otherwise.
k1 - k4 are coefficients. Note that these are not conventional dimensionless friction coefficents, as used
in the equation F = R. These coefficients are factors to be applied to the appropriate forces to give
frictional moments, and they must include effects of the bearing diameter etc.
k1, k2 and k3 must not be negative. k1 and k3 have dimensions of length, and the maximum value allowed is 0.025g / 9.81 where g is the acceleration due to gravity. k2 is non-dimensional and has a
maximum value of 0.025.
This moment is transformed into articulation, structure or global axes as appropriate before output.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
133
134
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
135
2
5 7
11
21
31
41
- --- -- ---- --------- --------- --------- --------|X|
| |OCIN|
|
-- -- -- ---- --------- --------- --------- --------|X|
| |APIR|
|
-- -- -- ---- --------- --------- --------- --------|X|
| |NPDW|
|
-- -- -- ---- --------- --------- --------- --------|X|
| |ISOW|
|
-- -- -- ---- --------- --------- --------- --------|X|
| |UDWD| XXXXX |cf
|cs
|I(z)
|
-- -- -- ---- --------- --------- --------- --------|X|
| |WIND| XXXXX |Speed Uz |Direction| Ref Ht |
-- -- -- ---- --------- --------- --------- --------| |
|
| |
|_Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on last data record in data category(A3)
The UDWD data record should be followed by UDWS data records defining dimensionless frequencies,
f, and spectral ordinates, U(f ), (as defined below). A user defined wind spectrum can consist of up to
50 UDWS data records.
2
-|X|
-|X|
-|X|
-|X|
--
5 7
11
21
31
-- -- ---- --------- --------- --------| |UDWS| XXXXX |f
|U(f)
|
-- -- ---- --------- --------- --------| |UDWS| XXXXX |f
|U(f)
|
-- -- ---- --------- --------- --------| |UDWS| XXXXX |f
|U(f)
|
-- -- ---- --------- --------- --------| |UDWS| XXXXX |f
|U(f)
|
-- -- ---- --------- --------- --------| |
|
| |
|_Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on last data record in data category(A3)
Note
The OCIN/APIR/NPDW/ISOW/UDWD/UDWS data records must be input before the WIND data
record in each wind spectrum definition.
136
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
This data record acts as a switch which changes the units for all following frequencies input to hertz in
wave spectra definitions.
If omitted, the program will expect all frequencies to be input in radians/second.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
137
This data record acts as a switch which changes the units for all following frequencies input to radians/second. Note that this data record is needed only if a HRTZ data record has been input, as the default
units for frequency are radians/second in wave spectra definitions.
(1) This is only required for a spreading sea. This value defines the power N of the wave spreading
function which is in the form of COS**N. This space should be left blank if N=2.
(2) This value defines the total spreading angle for the case when (1) is left blank (i.e. N=2). If N>2, this
space should be left blank as the total spreading angle will be set to 180 degrees.
For a non-spreading sea, both (1) and (2) should be left blank.
(3) This value is the direction of waves within a wave spectrum (for a spreading sea, this is the direction
in which the wave spectra values are the highest). All the following specifications of wave spectra take
this value until another SPDN data record is input. The next group of wave spectra will then take the
new value as the wave direction. It is therefore a mandatory requirement that this data record is first
data record of the wave spectrum data category, as the wave direction will be undefined until this data
record is input.
138
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
The CURR / WIND Data Records - Current and Wind Speed and Direction
The introduction of a SPDN data record thus enables the user to alter the value of the wave direction
for a complete group of wave spectra.
Note
Although there is no limit to the number of SPDN data records that may be input, a wave
spectrum can have only one direction. Hence the number of SPDN data records can never
exceed the number of wave spectra.
Wave spreading is not available in Aqwa-Naut at present.
The results output for plotting in the AGS are for the main direction only.
Within a spectral group, the default values of seed for each spectrum will be:
1. If there is no SEED data record or the SEED number is set to 0:
SEED = 1+(I-1) * 1,000,000
where I represents the i-th spectrum.
2. If a SEED data record is used for some but not all of the spectra:
SEED = N_SEED + J * 1,000,000
where N_SEED is the number of SEED data records defined in this group,
J ranges from 1 to the number of spectra which do not have a SEED data record.
17.8. The CURR / WIND Data Records - Current and Wind Speed and Direction
These data records are optional and their omission indicates that no current or wind is present. A value
of zero is then assumed for the wind and current speeds associated with each wave spectrum. In general both data records will be input before each wave spectrum.
2
5 7
11
21
31
41
- --- -- ---- ---------- --------- --------- --------|X|
| |CURR|XXXXXXXXXX|
|
|
|
- --- -- ---- ---------- --------- --------- --------|X|
| |WIND|XXXXXXXXXX|
|
|
|
- --- -- ---- ---------- --------- --------- --------| |
|
|
|
|
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
139
(1) The Data Record Keyword indicates whether parameters (2) and (3) apply to current or wind.
(2)-(3) These values are magnitude and direction of the current/wind speed for a wave spectrum. All
following specifications of wave spectra (q.v.) take these values until another CURR/WIND data record
is input.
(4) The reference height defines the height where wind speed is measured and is only required for wind
spectrum definition. If a non-zero height is entered, a wind spectrum will be used even if not specifically
requested. The default spectrum is Ochi and Shin. See also Wind Spectra Definition (p. 135).
Although there is no limit to the number of CURR/WIND data records that may be input, a wave spectrum
can only be associated with one current magnitude and direction and one wind magnitude and direction.
Hence the number of functional CURR or WIND data records can never exceed the number of wave
spectra.
(1)-(2) Start/Finish Frequency - The lowest/highest frequency at which the spectrum is defined. The
program will assume that the frequencies are in radians/sec, unless a HRTZ data record (see The HRTZ
Data Record - Change Units of Frequency to Hertz (p. 137)) has been used to change this to Hertz
(cycles/sec).
If these fields are left blank, defaults will be assumed as follows:
Start frequency = Peak frequency * 0.58
Finish frequency is evaluated numerically so that approximately 99% of the spectral energy is included.
140
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1)-(2) Start/Finish Frequency - The lowest/highest frequency at which the spectrum is defined. Aqwa
assumes that there is no wave frequency energy outside this range. If these fields are left blank, defaults
will be assumed as follows:
Start frequency: s = p- 3, subject to a max of 100 and min of 0.001
Finish frequency: f = p+ 3, subject to a max of 100 and min of 0.001
If f - s < 0.001 then s = 0.1, f = 6.0
(3) Standard deviation, . The minimum allowed value of is 0.08p
(5) Peak frequency, p
141
(1)-(2) Start/Finish Frequency - The lowest/highest frequency at which the spectrum is defined. The
program will assume that the frequencies are in radians/sec, unless a HRTZ data record (see The HRTZ
Data Record - Change Units of Frequency to Hertz (p. 137)) has been used to change this to Hertz
(cycles/sec).
If these fields are left blank, defaults will be assumed as follows:
Start frequency = Peak frequency * (0.58+(-1.0)*0.05/19.0)
Finish frequency is evaluated numerically so that approximately 99% of the spectral energy is included.
(3)-(5) Gamma, Alpha and Peak Frequency.
The JONSWAP wave spectrum can be used to describe a wave system where there is an imbalance of
energy flow (i.e. sea not fully developed). This is nearly always the case when there is a high wind speed.
It may be considered as having a higher peak spectral value than the Pierson-Moskowitz spectrum (The
PSMZ Data Record - Pierson-Moskowitz Spectrum (p. 140)) but is narrower away from the peak in order
to maintain the energy balance.
Parameterization of the classic form of the JONSWAP spectrum (with parameters of fetch and wind
speed) was undertaken by Houmb and Overvik (BOSS Trondheim 1976,Vol 1). These parameters are
used by Aqwa. These empirical parameters are termed gamma (3), alpha (4) and peak frequency (5)
(the frequency at which the spectral energy is a maximum).
142
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
17.13.The UDEF Data Record - User-Defined Wave Spectrum and the FINI
Data Record - User-Defined Spectrum Separator
This facility may be used to input any spectrum the user wishes. It is normally employed for input of
non-deterministic spectra such as tank spectra, recorded full-scale spectra, or simply where the formulated
spectrum is not yet available in Aqwa.
2
5 7
11
21
31
- --- -- ---- ---------- --------- --------|X|
| |UDEF|XXXXXXXXXX|
|
|
- --- -- ---- ---------- --------- --------| |
|
|
|
| |
|
|
|
| |
|
| (2)Spectral Ordinate(F10.0)
| |
|
|
| |
|
|
| |
|
(1)Frequency Value(F10.0) [*]
| |
|
| |
|_Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on last data record in data category(A3)
* See The HRTZ Data Record - Change Units of Frequency to Hertz (p. 137) and The RADS Data Record
- Change Units of Frequency to Radians/Second (p. 138).
(1) The value of frequency at which the value of the spectral ordinate (2) is given.
These frequencies must appear in ascending order (lowest first). The maximum number of frequencies
is 50.
(2) Value of the spectral ordinate at the frequency given in columns 21-30.
If more than one spectrum is input, then the FINI data record must be used to separate each set of
UDEF data records defining a spectrum. The format is shown below:
2
5 7
- --- -- ---|X|
| |FINI|
- --- -- ---| |
|
| |
|_Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on last data record in data category(A3)
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
143
(1) The number of spectral lines is the number of spectral ordinates required to define the spectrum in
order to achieve an accuracy comparable to the formulation of the analysis in which it is used. This is
necessarily a general definition as the program uses this number in several different ways.
This number of spectral lines will be used for ALL subsequent spectra in all groups, unless a particular
spectrum type requires a different number. For example, a spectrum imported using the IWHT data
record will always use 200 lines, irrespective of any value input using the NSPL data record.
(1) The number of spectral rasters is the number of spectral frequencies required to achieve sufficiently
small frequency differences to accurately define the spectrum at peak of response at the small frequency
values. These low frequencies are normally associated with slow drift motions.
144
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Note
As long as this file exists, the program will automatically read it in and use the wind speed and
direction defined in the file at each time step during the analysis. There is no data record
needed for implementing this feature.
The time defined in the WVT file does not need to match the time step defined in the corresponding DAT file as the program will interpolate the wind speed and direction when necessary,
using a cubic spline interpolation technique. When modeling periods of constant wind velocity
adequate data points must be provided to satisfy the interpolation method.
There is no limit on the length of the WVT file.
Note
As long as this file exists, the program will automatically read it in and apply the forces defined
in the file to the structures at each time step during the analysis. There is no data record needed
for implementing this feature.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
145
2
5 7
11
21
31
41
51
61
- --- -- ---- ----- ----- -------- -------- -------- -------- -------|X|
| |XSWL|X
|XXXXX|XXXXXXXX|XXXXXXXX|
|
|
- --- -- ---- ----- ----- -------- -------- -------- -------- -------| |
|
|
|
|
|
| |
|
|
|
|
|
146
71
---------|
|
---------|
|
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
The parameters used to define the spectra are the same as on the JONH, JONS and PSMZ data records,
except that the start and finish frequencies are calculated by the program.
----------------------------------------------------------------|Spectrum Type
| Parameter 1 | Parameter 2 | Parameter 3
|
----------------------------------------------------------------|JONSWAP (JONH)
|
gamma
|
Hs
| peak frequency |
----------------------------------------------------------------|JONSWAP (JONS)
|
gamma
|
alpha
| peak frequency |
----------------------------------------------------------------|Pierson Moskowitz |
Hs
|
Tz
| not used
|
----------------------------------------------------------------|Gaussian
|
sigma
|
Hs
| peak frequency |
-----------------------------------------------------------------
where the user may specify , Hs, fp (as p or Tp) and the direction of the spectrum, and:
Hs is significant wave height
fp is peak frequency in rad/sec
is standard deviation, always in rad/sec
Note
Hydrodynamic interaction between the normal and the cross swell spectrum in the case of
coupled Drift Force Coefficients or QTFs has been ignored; i.e the spectra are considered
totally separately and for totally linear systems will give the same results as the sum of the
results for each spectrum analyzed separately. Note that results such as significant motions
will be the square root of the sum of the squares.
2
5 7
11
21
- --- -- ---- -------- -------- -------- -------|X|
| |GATP|XXXXXXXX|XXXXXXXX|XXX
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
147
Note
All of the above and other places where wave spectra are used but cannot use cross swell
will cause a fatal error. This failure mechanism will be implemented in as many places as
possible where the wave spectra are used, but may not include every occurrence.
148
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
17.19. The IWHT Data Record - Import of Wave Height Time History
A time history series of wave elevation may be imported into Aqwa-Drift, in order to reproduce model
test wave conditions as accurately as possible. This facility uses the IWHT data record and an external
data file. See the section below for a discussion on the accuracy of the reproduction.
Import of the time series will also generate a user-defined spectrum, using a Fast Fourier Transform,
whose frequency range is based on a JONSWAP fit of the wave elevation spectral density. In AqwaLibrium, Aqwa-Fer and Aqwa-Naut this spectrum will be used in the same way as a normal user-defined
spectrum.
Note
As the phases of the spectral wavelets are allocated randomly, the input wave elevation
time-history will not be reproduced.
The number of spectral lines of 200 is mandatory for this facility.
Note
Current is ignored when calculating the phase wave forces on the structure and the wave
kinematics for Morison elements.
This data record is used to define the file from which the wave height time history is to be copied.
2
5 7
11
16
21
- --- -- ---- ----- ----- --------|X|
| |IWHT|XXXXX|XXXXX|
- --- -- ---- ----- ----- --------| |
|
|
| |
|
|_(1)Filename from which the data is to be copied
| |
|
| |
|_Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on Last data record in data category(A3)
(1) This is a file name with a .WHT extension. Up to 5 WHT files may be imported in a spectral group.
Imported wave elevations are treated as user-defined spectra, so they must be listed consecutively
within a spectral group and the last IWHT data record must be followed by a FINI data record. The
format is shown below:
2
5 7
- --- -- ---|X|
| |FINI|
- --- -- ---| |
|
| |
|_Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on last data record in data category(A3)
149
Note
The X_REF and Y_REF values are used in the calculation of the phase of the wave and are the
position where the wave elevation was measured. For example (in SI units), if the direction of
the wave is zero degrees (i.e. along the positive X-axis in the FRA) then values of X_REF/Y_REF
of 100.0/0.0 will indicate that the wave elevation was measured 100 metres downstream of the
0,0 wave reference point. Omission of these data will default the reference point to 0,0. i.e. the
wave elevation will be calculated using the origin of the FRA as the point at which the wave
elevation will be reproduced.
The Spectrum Name will be used for graphs and tables where appropriate throughout the
program.
Current speed and direction are needed for calculation of the wavelengths of the wavelets
used to reproduce the wave elevation. If omitted it is assumed that there is no current. Special
treatment is necessary for a profiled current. See CPRF Data Record - Profiled Current Velocity (p. 124) for additional information.
The duration of the time history in the file should be at least 7200s. This duration is necessary
in order to give sufficient resolution of low frequency resonant responses. If the file contains
less data than this, the data will be extended automatically up to 7200s, using a process of
mirroring and copying.
The maximum number of timesteps in the .WHT file is 150000.
Comments (starting with * in Column 1) may be added anywhere in the file.
150
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
onl
Accuracy of Wave Surface Modeling
In Aqwa-Drift the wave elevation time-history will be reproduced exactly, within the frequency range
of the fitted spectrum and subject to the limitations of roundoff error. This is achieved by multiplying
each of the spectral wavelets (as in standard Aqwa) by a different Low Frequency Perturbation (LFP)
Function . i.e.:
Wave elevation =Sigma(j=1,N) { a(j) cos(-w(j)t+k(j)x+(j)) * LFP(j)(t) }
Where:
N is the number of spectral lines (normal Aqwa N=NSPL)
j is the wavelet number
t is time
w is frequency (as normally output in the .LIS file)
is phase (as normally output in the .LIS file)
a(j) is the amplitude
k is the wave number
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
151
2
5 7
11
16
21
31
41
51
61
71
80
- --- -- ---- ----- ----- ---------- ---------- ---------- ---------- ---------- ---------|X|XXX| |SSDN|
|
|XXXXXXXXXX|
|
|
|
|
|
- --- -- ---- ----- ----- ---------- ---------- ---------- ---------- ---------- ---------|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|_Up to 5 Directions in Ascending Order (Degrees) (F10.0)
|
|
|
|
|
|
|
|_Finish Direction Number (Maximum = 41) (I5)
|
|
|
|
|
|_Start Direction Number (I5)
|
|
|
|_Compulsory Data Record Keyword(A4)
|
|_Optional User Identifier(A2)
This may then optionally be followed by specification of the individual weighting factors for each direction. This uses the SSWT data record and is described in detail below.
This is then followed by specification of the wave spectral density ordinates. The spectrum in each direction can be defined by ordinates at up to 50 frequencies: i.e. there can be up to 450 (41/5 * 50) UDDS
data records. Note that the same frequencies must be used for each spectrum.
2
5 7
11
16
21
31
41
51
61
71
80
- --- -- ---- ----- ----- ---------- ---------- ---------- ---------- ---------- ---------|X|
| |UDDS|XXXXX|XXXXX|
|
|
|
|
|
|
- --- -- ---- ----- ----- ---------- ---------- ---------- ---------- ---------- ---------| |
|
|
|
|
|
|
|
| |
|
|
|
|
|
|
|
| |
|
|
|_(3)Spectral ordinate values for each direction (F10.0)
| |
|
|
| |
|
|_(2)Frequency of Spectral Ordinates (F10.0)
| |
|
| |
|_(1)Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on last data record in data category(A3)
152
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
NOTES
The results output for plotting in the AGS are for the main direction only.
If there is no SSWT data record the weighting factors i.e. the contribution of each spectrum to the total
force (Aqwa-Librium) or motion (Aqwa-Fer) are calculated by a simple trapezoidal integral and will add
up to unity. The weighting factors are calculated as follows:
Wt(i) = 0.5(Dn(i+1)-Dn(i-1))/(Total direction range)
In the case of the 1st direction, the weighting function is:
Wt(1) = 0.5(Dn(2)-Dn(1))/(Total direction range)
and for the last direction of N total directions:
Wt(N) = 0.5(Dn(N)-Dn(N-1))/(Total direction range)
In general, it is up to the user to input the spectral ordinates factored so the integral of the spreading
function divided by the spread range is unity.
As an example, if a carpet spectrum is defined to model a cos2 spreading function, then the spectral
ordinates must be factored by 2.0, as the integral of cos2 divided by the direction range is 0.5.
Alternatively the user may exercise total control over the contribution of each spectrum at each direction
by specifying the individual weighting factors. This information must be input directly after the SSDN
input and has the following format.
2
5 7
11
16
21
31
41
51
61
71
80
- --- -- ---- ----- ----- ---------- ---------- ---------- ---------- ---------- ---------|X|
| |SSWT|XXXXX|XXXXX|XXXXXXXXXX|
|
|
|
|
|
- --- -- ---- ----- ----- ---------- ---------- ---------- ---------- ---------- ---------|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|_(1)Spectral weights for each direction (F10.0)
|
|
|
|_Compulsory Data Record Keyword(A4)
|
|_Optional User Identifier(A2)
For each direction, a spectral weight of 1.0 will give the same results as a single long-crested wave. This
method of input can therefore be used to describe a long crested wave with 2 simultaneous cross swell
spectra with 3 direction entries. Note that if columns 31-80 are left blank then the default weighting
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
153
Note
It is important to note that the drift force coefficients or QTFS are not the same as those
used by the COSN spreading functions. Hence the low frequency response in Aqwa-Fer will
not be the same. For the COSN the QTFs used for each direction are evaluated by integrating
the QTFs over the 'effective' spread range specified. This is not the case with user defined
directional spread and the QTFs are simply interpolated at the user defined directions.
(1) NODR specifies that the 2nd order drift forces are not to be calculated for the spectrum being
defined. If more than one spectrum is input then a NODR data record is required for every spectrum
for which drift forces are to be omitted.
(2) Where a User Defined Spread Seas Spectra is defined, this item specifies a direction for which 2nd
order drift forces are to be omitted. For all other spectrum types this data is not required.
154
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
LIMITATIONS
1. Program applicability:
Aqwa-Line - the SPGR card is recognized, but there can be only one spectral group containing only
one spectrum.
Aqwa-Librium - up to 20 spectral groups.
Aqwa-Fer - up to 20 spectral groups.
Aqwa-Drift - only one spectral group is allowed.
Aqwa-Naut - only one spectral group is allowed.
2. The overall maximum total number of spectra in all groups is 20.
3. In the following sections the term sub-spectrum is used.
A user-defined 2D (carpet) spectrum is a single spectrum, but it may consist of up to 41 sub-spectra.
A spectrum with a spreading function is modeled in Aqwa with 7 sub-spectra, at Gauss integration
points.
A spectrum with cross-swell contains 2 sub-spectra.
4. A spectral group may contain up to 20 spectra, subject to the overall limit above. Therefore, if a group
contains 20 spectra, there can be only one group. The spectra can be any of the available types, i.e.
JONSWAP, Pierson-Moskowitz, Gaussian, user-defined or imported using the IWHT data record.
5. There may be up to 5 IWHT data records, but these must all be in the same spectral group. Within a
group they must be listed consecutively and the list terminated with a FINI data record.
6. A spectral group may contain up to 41 sub-spectra in total.
7. A spectral group may contain only one definition of wind and/or current.
8. If there is no wind or current defined in a spectral group, the values will be taken from the last WIND or
CURR data record read. This is consistent with behavior in earlier versions of Aqwa.
An example of a spectral group definition is shown below. This group contains a wind spectrum, a
current, a JONSWAP spectrum and a user-defined spectrum. It also makes use of both spectrum and
group names.
SPGR
SGNM
NPDW
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
155
12.000
126.00
13.0
1.000
30.00
20.0
Spectrum 1 - JONSWAP
0.2
1.5
2.0
1.00000
-135.0
Spectrum 2 - simple user-defined spectrum
0.100
0.000
0.300
2.000
0.700
2.000
0.900
0.000
0.6
*
SPGR
2
5 7
11
21
31
41
- --- -- ---- ---------- ---------- ---------- ---------|X|XXX| |TOWS|XXXXXXXXXX|XXXXXXXXXX|
|
|
- --- -- ---- ---------- ---------- ---------- ---------|
|
|
|
|
|
|
|
|
|
|
|_(2) Tow Direction (E10.0)
|
|
|
|
|
|_(1) Tow Speed (E10.0)
|
|
|
|_Compulsory Data Record Keyword (A4)
|
|_Optional User Identifier (A2)
156
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
157
(1) The wave ramp is introduced to reduce the transient motion of the structure at the beginning of a
time domain analysis. When this data record is defined, the wave ramp will take effect from t = 0.0 to
t = tw, during which a wave ramp factor f (0.0 < f < 1.0) will be calculated and then multiplied to the
wave frequency forces.
The wave ramp factor f is decided by
158
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
159
160
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
161
(1) The mooring line code indicates which type of linear mooring line the user wishes to model. Note
that the term 'linear' is used generally to denote that the tension is either constant or linearly proportional to the extension. The cable codes available at present are:
LINE - a conventional linear elastic cable
WNCH - a winch adjusted to constant tension
FORC - a constant force
Note that these mooring lines are assumed to have no mass and are therefore represented geometrically
by a straight line. See below for further details.
(2) This is the number of the structure to which the cable is attached and must correspond to one of
the structures defined within Data Category 2. If 1 is input, this will correspond to the structure defined
in Data Category ELM1. If '2' is input, this will correspond to the structure defined in Data Category
ELM2 etc. Note that Structure 0 (in other words, a fixed node) is an ILLEGAL structure and will produce
an error (see (4)).
(3) This is the node number whose position is the attachment point of the end of the mooring line on
the structure specified (1). The position of this node on the structure (2) must have been defined in
Data Category 1.
162
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
163
164
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
2
5 7
11
31
41
51
61
- --- -- ---- ----------------- --------- --------- --------- ----|X|
| |POLY|XXXXXXXXXXXXXXXXX|
|
|
|
...
- --- -- ---- ----------------- --------- --------- --------------| |
|
|
|
|
|
| |
|
|--------------------------------| |
|
|
| |
|
|
| |
|
(1)-(5) 5 Coefficients of the Polynomial (5E10.0)
| |
|
| |
|_Compulsory Data Record Keyword (A4)
| |
| |_Optional User Identifier (A2)
|
|_Compulsory END on last data record in Data Category (A3)
(1)-(5) These values represent the coefficients of the polynomial which defines the force in the mooring
as a function of extension or compression.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
165
Note
When used to define the stiffness properties of a fender the force and displacement used
are compressive instead of tensile. For a fender with linear stiffness P1 will still be positive.
2
5 7
11
16
21
26
31
41
51
61
71
- --- -- ---- ----- ----- ----- ----- --------- --------- --------- --------- --------|X|
| |NLIN|
|
|
|
|
|
|
|
|
|
- --- -- ---- ----- ----- ----- ----- --------- --------- --------- --------- --------| |
|
|
|
|
|
|
|
|
|
| |
|
|
|
|
|
|
|
|
|
| |
|
|
|
|
|
|
|
| (9) Paying-out Friction
| |
|
|
|
|
|
|
|
|
Factor Fp (F10.0)
| |
|
|
|
|
|
|
|
|
| |
|
|
|
|
|
|
| (8) Winding-in Friction
| |
|
|
|
|
|
|
|
Factor Fw (F10.0)
| |
|
|
|
|
|
|
|
| |
|
|
|
|
|
|
|
| |
|
|
|
|
|
| (7) Unstretched Length (F10.0)
| |
|
|
|
|
|
|
| |
|
|
|
|
|
(6) Winch Tension Ts (F10.0)
| |
|
|
|
|
|
| |
|
|
|
|
|_(5) Node Number (I5)
| |
|
|
|
|
| |
|
|
|
|_(4) Structure Number (I5)
| |
|
|
|
| |
|
|
|_(3) Node Number (I5)
| |
|
|
| |
|
|_(2) Structure Number (I5)
| |
|
| |
|_(1) Nonlinear Mooring Line Code (A4)
| |
| |_Optional User Identifier (A2)
|
|_Compulsory END on last data record in Data Category (A3)
166
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Fp
Description
Calculation
>0.0
>0.0
Mean tension
(Tmin+Tmax)/2
900t
>0.0
< 0.0
Winding-in tension
Tmin
700t
< 0.0
>0.0
Paying-out tension
Tmax
1100t
< 0.0
< 0.0
Winch tension
1000t
For Aqwa-Librium and Aqwa-Fer, this tension value will not change unless the line goes slack. (Range
is less than the initial length specified.)
Aqwa-Drift and Aqwa-Naut will vary the tension according to whether the range (distance between the
anchor and vessel attachment point) is increasing or decreasing. If the range is less than the initial
length specified the line becomes slack and the tension is zero.
The initial paid-out length is initialized to give the specified tension in the initial position specified in
Data Category 15 at the beginning of the time history.
* If the range is increasing and the tension is at Tmax then the winch is paying out and the tension remains
at Tmax.
* If the range is decreasing and the tension is at Tmin then the winch is winding in and the tension remains at Tmin as long as the range is greater than the initial length specified.
* At all other times the line acts as a normal NLIN(POLY) line.
167
2
5 7
11
16
21
31
41
51
- --- -- ---- ----- ----- ----- ----- --------- --------- ----------|X|
| |COMP|
3|
10|11112|
2|
|
|
|
- --- -- ---- ----- ----- ----- ----- --------- --------- ----------| |
|
|
|
|
|
|
|
|_(7) Sea bed slope (F10.0)
| |
|
|
|
|
|
|
|
| |
|
|
|
|
|
|
|_(6) Z Maximum Value (F10.0)
| |
|
|
|
|
|
|
| |
|
|
|
|
|
|_(5) Z Minimum Value (F10.0) Must be > zero
| |
|
|
|
|
|
| |
|
|
|
|
|_(4) Number of lines in the cable
| |
|
|
|
|
composite (I5) (NO DEFAULT) (Max=10)
| |
|
|
|
|
| |
|
|
|
|_(3) Warning flags and Symmetry forcing (I5)
| |
|
|
|
(Def = all warnings issued)
| |
|
|
|
| |
|
|
|_(2) Number of X coordinates (I5) (Def =40,Max = 40) or 600 divided by (1), whichever is les
| |
|
|
| |
|
|_(1) Number of Z coordinates (I5) (Def =15, Max = 20) or 600 divided by (2), whichever is less
| |
|
168
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) The number of Z values to be used in creating the 2-dimension load extension database for this
composite mooring line
(2) The number of X values to be used in creating the 2-dimension load/extension database for this
composite mooring line.
For composite catenaries linking two structures (instead of sea bed and structure), it is strongly recommended that the full database size (600 points) be used. The default numbers, when these fields are
left blank, for Z and X are 15 and 40 respectively. Also the Z becomes the angular coordinate and the
X the radial coordinate, and the program will ignore user specifications for the Z range and use the
default values (-90 to + 90 degs). The sea bed slope will also be ignored for catenaries between structures
as they are not allowed to touch the sea bed.
(3) This parameter indicates whether warnings should be issued when the position of the attachment
point of the mooring line relative to the anchor point is outside the range of the database created for
this mooring line characteristic. Warnings about the degree of symmetry are not issued, as the stiffness
matrix is automatically symmetric for composite cables. See The LE2D Data Record (p. 185) for further
details on stiffness matrix asymmetry.
This parameter should be thought of as five separate flags, which indicate the following:
Flag
Column
Meaning
21
0 = Warnings are issued when the X position exceeds the range specified
22
0 = Warnings are issued when the X position is below the range specified
23
0 = Warnings are issued when the Z position exceeds the range specified
24
0 = Warnings are issued when the Z position is below the range specified
25
0 = Symmetry Forcing
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
169
170
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) Mass per unit length of the section of the composite mooring line.
(2) The equivalent cross sectional area of the mooring line. It is often more convenient, especially with
wire lines, to specify this parameter so the buoyancy of the line may be calculated and subtracted from
the structural weight to give the 'weight in water'. This parameter may also be specified as zero if the
mass per unit length (1) is input as the mass of the line less the mass of the displaced water per unit
length (this does not apply to the cases when cable dynamic analysis is required, for which a non-zero
equivalent cross section area must be defined).
(3) The stiffness of the line is specified in terms of EA, where E is Youngs modulus and A is the cross
sectional area of the line. The default value is chosen to give a typical value based on the mass/unit
length. Clearly this may be in error if the mass per unit length specified (1) includes buoyancy effects.
(4) The maximum tension is the highest value of tension that should used in the database created for
this composite mooring line.
(5) The length of the section of the composite mooring line. Note that values less than 1 are not accepted.
N.B. For a composite mooring line containing more than one ECAT data record, the definition of ECATs
should start from the anchor point. If a composite mooring line links two structures the ECAT data records
can start from either end, but the start must correspond to the second structure on the NLIN data record.
171
(2)-(4) These coefficients define a nonlinear variation in axial stiffness. The stiffness is calculated using
the formula
EA() = EA(const) + k1.+ k2.2 + k3.3
(0 tmax)
( > tmax)
where
EA() = EA as a function of strain
EA(const) = the EA value on the ECAT data record (The ECAT Data Record (p. 170))
k1, k2, k3 = coefficients input on the ECAX data record (The ECAX Data Record - Nonlinear Axial Stiffness
Polynomial (p. 171))
= linear strain L/L
tmax = strain at Tmax
Tmax = maximum tension specified on the preceding ECAT
Note
The stiffness must increase monotonically over a tension range from zero to the maximum
tension specified on the preceding ECAT data record (Tmax).
For tension greater than Tmax, the stiffness is assumed to increase linearly with strain, with
slope equal to the slope at Tmax.
172
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(2) This is the structural mass of the buoy or clump weight. This must be smaller than the mass of displaced water (3) for a buoy, or larger for a clump weight. This can be positive, zero or negative.
(3) The mass of water displaced, i.e. the buoyancy/gravity. This can be positive, zero or negative.
(4) Total (constant) added mass, i.e. not added mass coefficient. Applicable to Cable Dynamics only.
(5) The drag force (cable dynamics only) will be in the direction of the relative velocity of the fluid, VR.
The magnitude of the force is given by
FD = 0.5 (CDA )(VR )(|VR|)
where CDA = Drag coefficient projected area.
173
174
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
2
5 7
11
16
21
26
31
41
51
61
71
81
- --- -- ---- ----- ----- ----- ----- ---------- ---------- ---------- ---------- ---------|X|
| |ECAH|XXXXX|XXXXX|XXXXX|XXXXX|
|XXXXXXXXXX|
|
|
|
- --- -- ---- ----- ----- ----- ----- ---------- ---------- ---------- ---------- ---------| |
|
|
|
|
|
|
| |
|
|
|
|
|
|_(6)Inline Drag Coefficie
| |
|
|
|
|
|
| |
|
|
|
|
|_(5) Equivalent Diameter for Drag D
| |
|
|
|
|
| |
|
|
|
|_(4) Transverse Drag Coefficient Cd (F10.0)
| |
|
|
|
| |
|
|
|_(3) Leave Blank (for future use)
| |
|
|
| |
|
|_(2) Added Mass Coefficient Ca (F10.0)
| |
|
| |
|_(1) Data Record Keyword(A4)
| |
| |_Optional User Identifier (A2)
|
|_Compulsory END on last data record in Data Category (A3)
(2) Added mass is calculated by RHO*Ca*A per unit length in which RHO is the water density and A is
the equivalent cross section area of the ECAT The COMP/ECAT Data Records - Composite Catenary
Mooring Line (p. 168). In other words, the added mass is equal to the displaced mass of water multiplied
by Ca. For cable dynamic analysis, the equivalent cross section area A must not be omitted in the ECAT
definition. The default is 1.0.
(4) Transverse drag force is calculated by 0.5*RHO*Cd*V2*De per unit length where V is the relative
transverse velocity. The default is 1.0.
(5) Equivalent diameter for drag. This allows the drag to be based on a different diameter from the
added mass. The default is
(6) Inline drag force is calculated by 0.5*RHO*Cx*V2*De per unit length where V is the relative inline
velocity. The default is 0.025.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
175
(2) For a prismatic section the bending stiffness is the product of the elastic modulus and the 2nd moment of area. For a chain or laid hawser an equivalent value must be defined.
Note
It is not possible to input a very high stiffness to create a straight rod. The maximum bending
stiffness is related to the axial stiffness as follows.
EImax = (EA)(L2/30)
This check is carried out for the individual elements into which the cable is split, so the length
(L) is the length of the element. This is not available to the user but is approximately the
total length of the line divided by the total number of elements. See the NCEL data record.
2
5 7
11
16
21
26
31
41
51
61
71
- --- -- ---- ----- ----- ----- ----- --------- --------- --------- --------- --------|X|
| |NCEL|
|
|
|
|
|
|
|
|
|
- --- -- ---- ----- ----- ----- ----- --------- --------- --------- --------- --------| |
|
|
| |
|
|
| |
|
|_(2) Number of elements required for each line (I5)
| |
|
| |
|_(1) Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier (A2)
|
|_Compulsory END on last data record in Data Category (A3)
Each NLID data record that follows this specification will have the required number of elements.
176
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
2
7
11
16
21
31
41
51
- ---- ---- ---- ---- --------- --------- --------- --------|X|
|LBRK|
|
|XXXXXXXXX|
|
|
|
- ---- ---- ---- ---- --------- --------- --------- --------|
|
|
|
|
|
|
|
|
|
|
|_(4) Breaking tension at 2nd structure (AQWA-DRIFT
|
|
|
|
|
|
|
|
|
|_(3) Breaking tension at 1st structure (AQWA-DRIFT/NAUT)
|
|
|
|
|
|
|
|_(2) Time at which to break the mooring line (AQWA-DRIFT/NAUT)
|
|
|
|
|
|_(1) Mooring line number to break (I5)
|
|
|
|_Spectrum number (AQWA-LIBRIUM/FER) (max 20)
|
|_LBRK Data Record Keyword
(1) The mooring line number is taken from the order in which the lines are defined.
(2) The time has two meanings.
If no tension is defined the line will break at this time.
If a tension is also defined, the line will break the 1st time the tension exceeds the specified value, after
this time.
(3) The mooring line will break when the tension at the 1st structure on the LINE/NLIN/NLID data record
is greater than this value, if the time is greater than the specified breaking time (2).
(4) The mooring line will break when the tension at the 2nd structure on the LINE/NLIN/NLID data record
is greater than this value, if the time is greater than the specified breaking time (2).
Note
Logic for breaking lines.
not
not
not
not
TIME
given
given
given
given
t1
t1
t1
or
or
or
or
0
0
0
0
TENSION 1
not given or 0
F1
not given or 0
F1
not given or 0
F1
F1
TENSION 2
not given or
not given or
F2
F2
not given or
not given or
F2
0
0
0
0
RESULT
Always broken
Breaks when tension
Breaks when tension
Breaks when tension
Breaks at time=t1
Breaks when tension
Breaks when tension
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
1 > F1
2 > F2
1 > F1 or tension 2 > F2
177
(1) The PULY data record indicates that a pulley is positioned on a line. A maximum of 2 pulleys is allowed
for each pulley set (i.e. two pulleys per LINE), with the first PULY data record following a LINE data record,
i.e. applying to the most recently defined LINE. In the case of a 2 pulley set the second PULY data record
will follow the first PULY data record. A PULY has the effect of intersecting a LINE and will effectively
extend the LINE from the LINE's node to the second PULY node.
(2) This is the structure number to which the pulley is attached and must correspond to one of the
structures defined in Data Category 2. If '1' is input, this will correspond to the structure defined in Data
Category ELM1. If '2' is input, this will correspond to the structure defined in Data Category ELM2 etc.
Note that Structure '0' (i.e. a fixed node) is an ILLEGAL structure and will produce an error (see (4)).
(3) This is the node number whose position is the attachment point of the pulley on the structure specified (1). The position of this node on the structure (2) must have been defined in Data Category 1.
Items (2) & (3) must appear on the previous LINE data record (see examples).
(4)-(5) This structure number (4) and its corresponding node number (5) define the attachment position
of the other end of the mooring line. The position of this node on the structure (4) must have been
defined in Data Category 1.
178
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(6) This value represents the angle at which the pulley lies about the Z axis (FRA). The default angle (0)
places the pulley in the orientation shown in the diagram.
(7) The radius is used to indicate the radius of the pulley, but is not used for the purpose of the calculation. The radius of the pulley will be displayed in the Aqwa Graphical Supervisor.
(8)-(9) There are two possible formulations for the friction. In each case the friction of the pulley is
represented by T2/T1, where T2 is the larger tension and T1 the smaller. T2/T1 should be defined for the
situation where the line turns through 180 around the pulley.
T2/T1 in columns 51-60 represents bearing friction of a pulley rotating around an axle
T2/T1 in columns 61-70 represents sliding between the rope and the pulley, for example if the pulley
represents a fairlead.
If values are entered in both fields the 2nd definition overrides the 1st.
Pulley with bearing friction: a friction factor is calculated such that = (T2-T1) / (T2+T1). The friction
is then varied depending on how far around the pulley the line passes.
Pulley with sliding friction: a friction factor is calculated such that (T2/T1) = exp(). The friction is
then varied depending on how far around the pulley the line passes.
In each case T2/T1 must be in the range 1 T2/T1 2.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
179
Examples of Connectivity
The following three cases demonstrate the required input of connectivity for systems including pulleys.
CASE 1 Is a standard LINE connected between a single structure and earth. The typical connectivity input
for this system would be:
LINE
101
201
...
CASE 2 Is a LINE connected between a single structure and earth. A pulley is inserted on the line and
extends the line to node 102. The typical connectivity input for this system would be:
LINE
PULY
1
1
101
101
0
1
201
102
...
...
CASE 3 Is a LINE connected between a two structures and earth, with a 2 pulley set on the LINE. The
typical connectivity input for this system would be:
LINE
PULY
PULY
1
1
1
101
101
102
0
1
2
201
102
201
...
...
...
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
2
5 7
11
16
21
26
31
41
51
61
- --- -- ---- ----- ----- ----- ----- --------- --------- --------X|
| |DWT0|
|
|
|
|
|
|
- --- -- ---- ----- ----- ----- ----- --------- --------- --------|
|
|
|
|
|
|
|
|
|
|
|_ (2) Starting Time of the Winching Action
|
|
|
|_(1) Compulsory Data Record Keyword(A4)
|
|_Optional User Identifier (A2)
(2) The line will will be treated as a standard linear elastic line until this time is reached or exceeded at
any particular timestep.
2
5 7
11
16
21
26
31
41
51
61
71
81
- --- -- ---- ----- ----- ----- ----- --------- --------- --------- --------- --------X|
| |LNDW|
|
|
|
|
|
|
|
|
|
- --- -- ---- ----- ----- ----- ----- --------- --------- --------- --------- --------| |
|
|
|
|
|
|
|
|
|
|
| |
|
|
|
|
|
|
|
|
|
|_(10)T_Max. Maximum Tensi
| |
|
|
|
|
|
|
|
|
|
| |
|
|
|
|
|
|
|
|
|_(9) Speed. Non-zero positive
| |
|
|
|
|
|
|
|
|
or negative winch speed (F10.0)
| |
|
|
|
|
|
|
|
|
| |
|
|
|
|
|
|
|
|_(8) L_End. Final Length (F10.0)
| |
|
|
|
|
|
|
|
| |
|
|
|
|
|
|
|_(7) L_Beg. Initial Length (F10.0)
| |
|
|
|
|
|
|
| |
|
|
|
|
|
|_(6) K_Beg. Initial Stiffness (E10.0)
| |
|
|
|
|
|
| |
|
|
|
|
|_(5) 2nd Node Number (I5)
| |
|
|
|
|
| |
|
|
|
|_(4) 2nd Structure Number (I5)
| |
|
|
|
| |
|
|
|_(3) 1st Node Number (I5)
| |
|
|
| |
|
|_(2) 1st Structure Number (I5)
| |
|
| |
|_(1) Compulsory Data Record Keyword (A4)
| |
| |_Optional User Identifier (A2)
|
|_Compulsory END on last data record in Data Category (A3)
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
181
CAVEATS/ASSUMPTIONS
Note
The exact length of the line at any time will depend on the previous motions which have been
encountered by the structure connected to this line. This in turn means that the length of the
line has 'memory'. The implication of this is that in situations where initial or specific positions
are used in Aqwa, the line length cannot be determined and will be assumed to be the initial
length. An example of this is the hot start. A warning message to this effect will be issued in
these cases.
The resolution of switching on and off the drum winch can only be the same as the time step.
This means that the winch drum can only be switched at the beginning or the end of a time
step and not in the middle. In order to conserve energy/momentum in the equations of motion,
the length of the line can only be changed in steps of (time step)(speed of the winch). The
182
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
2
5 7
11
16
21
26
31
41
51
61
- --- -- ---- ----- ----- ----- ----- --------- --------- --------X|
| |DWAL|
|
|
|
|
|
|
- --- -- ---- ----- ----- ----- ----- --------- --------- --------|
|
|
|
|
|
|
|
|
|
|
|_ (2) Additional length required.
|
|
|
|_(1) Compulsory Data Record Keyword(A4)
|
|_Optional User Identifier (A2)
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
183
(1) The size of the fender is the distance between the contact planes at which the fender just touches
both of them.
(2) This is the friction coefficient . The friction force is given by F = R, where R is the normal reaction.
The maximum value of is 0.5. See note below.
(3) Material (or structural) damping coefficient . Damping is modeled as linear material damping, where
the damping coefficient is ()(stiffness). Damping is only applied in the direction perpendicular to the
contact points.
Note
Fender friction works best in situations where the friction force is smaller than other forces
in the same direction. Friction will slow down relative motion between two structures, but
is not suitable for keeping them fixed together - there is no "stiction". When the relative velocity changes sign the friction force must also change sign, but to avoid an instantaneous
change in force (and therefore an instantaneous change in acceleration) a smoothing function
is applied. This means that when the relative velocity is very small the friction force is also
small, and the structures can move relative to each other.
When an FLIN data record is input, the program assumes that the nonlinear properties correspond to
those specified by the most recently input nonlinear properties (POLY and FEND) data records. Failure
to input these properties before a FLIN data record means that they are undefined and will result in an
error. Note that one nonlinear properties data record may apply to several FLIN data records. This not
184
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
185
(1) The number of values of Z to be input on the ZRNG data record and hence the number of values
on each XRNG, HRNG and VRNG data records.
(2) The number of values XRNG, HRNG and VRNG data records making up the load extension curves.
This data record together with (1) form a database which is represented by a matrix of values, e.g. for
the default size of the database Z(5), X(12,5), H(12,5), V(12,5) would be created.
(Further details may be found in The 2 Dimensional Load Extension Characteristics (p. 187).)
(3) This parameter indicates whether warnings should be issued when the position of the attachment
point of the mooring line relative to the anchor point is outside of the range specified in the ZNRG and
XRNG data records. In addition, warnings about the degree of symmetry of the resulting stiffness matrix
can also be issued, omitted or the matrix can be forced to be symmetrical. This is explained further
below.
The parameter should be thought of as five separate flags which indicate the following:
If flag 1 = 0 or is blank, Warnings are issued when the X position exceeds the range specified
If flag 2 = 0 or is blank, Warnings are issued when the X position is below the range specified
If flag 3 = 0 or is blank, Warnings are issued when the Z position exceeds the range specified
If flag 4 = 0 or is blank, Warnings are issued when the Z position is below the range specified
If flag 5 = 0 or is blank, Warnings are issued when the degree of asymmetry of the stiffness matrix
formed from the database is considered unacceptable.
If any flag value is non-zero then warnings are not issued for that flags conditions.
For the stiffness asymmetry, if the value is greater than unity, symmetry of the stiffness matrix will be
imposed. Clearly in this case no warning asymmetry warning message will be issued. Note that in time
history programs (Aqwa-Naut/Drift) the stiffness matrix is not used except for output of information to
the user.
DH
-DX
DH
-DZ
H
X
DV
-DX
DV
-DZ
where
H = the horizontal force
V = the vertical force
D/DX, D/DZ = differential operators
186
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1)-(5) ZRNG-Values on this data record represent values of Z in the Fixed Reference axis at which the
load extension characteristics are defined on the XRNG, HRNG, VRNG data records. This value is the
expected range of the attachment points in the analysis. These values will normally be the Z distance
between the attachment point on the vessel and the anchor plus or minus the expected amplitude of
motion of the NLIN line using this composite, i.e. if the sea bed and anchor are at -100 and the Z position
of the attachment point when the vessel is in equilibrium is at -10, ZMIN and ZMAX would be 80 and
100 respectively assuming a maximum amplitude of motion of 10.
XRNG-Values on this data record represent values of X (horizontal distance in the fixed reference axis
FRA) between the attachment and anchor position axis at which the horizontal and vertical forces are
defined on the HRNG, VRNG data records.
HRNG/VRNG-Values on these data records represent values of the horizontal/vertical forces at the previously specified X and Z values.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
187
VERTICAL
FORCE
IN
COLUMN 1
^
|
|
|
*
|
*
|
*
|
*
|
*
|*
*
|
|____________________>
X-VALUES IN COLUMN 1
VERTICAL
FORCE
IN
COLUMN 2
^
|
|
|
*
|
*
|
*
|
*
|
*
|*
*
|
|____________________>
X-VALUES IN COLUMN 2
VERTICAL
FORCE
IN
COULMN 3
^
|
|
|
*
|
*
|
*
|
*
|
*
| * *
|
|____________________>
X-VALUES IN COLUMN 3
If the above graphs represent a complete set of values to define the 2-D load extension characteristics
of a mooring line then the ZRNG data record will contain three values (there is only 1 ZRNG data record
in each set). If the number of values of X, H and V are 7 (as shown above) then the complete set of data
record required to define the 2-D load/extension curve is as follows:
Values of Z, X, H, V
1-Data Record
ZRNG . . . . . . . xxx
xxx
xxx
xxx
xxx
xxx
xxx
xxx
xxx
1-Data Record
188
LE2D
(3 values)
7 . . .
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
19.15. The PRIC Data Record - Print Initial Condition of Mooring Lines
This data record causes details of the mooring lines to be printed in the initial position at the start of
the analysis.
2
5 7
- --- -- ---- --X|
| |PRIC|
- --- -- ---- --|
|
|
|_Compulsory Data Record Keyword (A4)
2
5 7
- --- -- ---- -X|
| |FINI|
- --- -- ---- -| |
|
| |
|
| |
|_Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier (A2)
|
|_Compulsory END on last data record in data category (A3)
2
5 7
11
21
- --- -- ---- ---------- ------------------X|
| |FILE|XXXXXXXXXX|
- --- -- ---- ---------- ------------------|
|
|
|
|
|_(1)Name of the file containing mooring data
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
189
(1) The name can include the path of the file. If no path is specified, the local directory will be assumed.
2
5 7
31
41
- --- -- ---- -------------------- --------- --------|X|
| |SWIR|XXXXXXXXXXXXXXXXXXXX|
|
|
- --- -- ---- -------------------- --------- --------| |
|
|
|
| |
|
|
|
| |
|
|
|_(2) Asymptotic Offset (F10.0)
| |
|
|
| |
|
|_(1) Asymptotic Stiffness (E10.0)
| |
|
| |
|_Compulsory Data Record Keyword (A4)
| |
| |_Optional User Identifier (A2)
|
|_Compulsory END on last data record in Data Category (A3)
(1)-(2) These fields contain the values of the two constants in the equation defining the tension of the
line as a function of the extension (see below). Values must be specified for both fields.
This facility enables the user to input the physical properties (constants defining the tension/extension
curve) of a steel wire mooring line. Note that this data record does not define any mooring lines having
these properties. This information must be supplied on following NLIN data records.
The mooring line properties defined on the SWIR data record will apply to all NLIN data records that
follow until another nonlinear property data record is input. (Note that POLY is also a nonlinear property
data record.)
Tension in a steel wire mooring line is given by:
T = k (x - d (tanh(x/d)))
where
x = extension of mooring line
k = asymptotic stiffness (constant)
d = asymptotic offset (constant)
The names of the constants k and d arise from the fact that, at large values of extension, tanh(x/d) tends
to unity and the equation tends to the asymptotic form:
T = k (x - d)
190
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
2
5 7
11
16
21
26
- --- -- ---- ----- ----- ----- ----|X|
| |TELM|
|
|
|
|
- --- -- ---- ----- ----- ----- ----| |
|
|
|
|
|
| |
|
|
|
|
|
| |
|
|
|
|
|
| |
|
|
|
|
|
| |
|
|
|
|
|_(4) Geometric Group Number (I5)
| |
|
|
|
|
| |
|
|
|
|_(3) Material Group Number (I5)
| |
|
|
|
| |
|
|
|_(2) Node Number 2 (I5)
| |
|
|
| |
|
|_(1) Node Number 1 (I5)
| |
|
| |
|_Compulsory Data Record Keyword (A4)
| |
| |_Optional User Identifier (A2)
|
|_Compulsory END on last data record in data category (A3)
(1)-(2) These are the nodes input in Data Category 1 and define the length of the tether element only.
The first element is considered to have Node 1 attached to the anchor, for installed tethers. It is the
trailing node for towed tethers.
(3) The material group number (input in Data Category 3) for this element. There are two parameters
input for the material properties of tether elements. These are density and Young's Modulus of elasticity,
i.e.
2
5 7
16
21
31
- --- -- --------- ----- ---------- ---------|X|
| |XXXXXXXXX|
|
|
|
- --- -- --------- ----- ---------- ---------|
|
|
|
|
|
|
|
|
|
|
|_(3) Young's Modulus (E10.0)
|
|
|
|
|
|
|
|
|_(2) Material Density or Mass (E10.0)
|
|
|
|
|
|_(1) Material Group Number (I5)
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
191
(4) The Geometric Group for this element. Geometric properties for tether elements are the same as for
TUBE elements, except that tether elements cannot be free flooding or have end cuts, i.e. they have
diameter, wall thickness, drag and added mass coefficients specified.
19.19.2. The TSPA/TSPV Data Records - Tether Anchor and Vessel Springs
Only one TSPA and one TSPV data record may be input for each tether.
2
5 7
11
31
41
51
- --- -- ---- ----- ... ---------- ---------- ---------|X|XXX| |TSPA|
1|
|
|
|
|
- --- -- ---- ----- ... ---------- ---------- ---------|X|XXX| |TSPV|
1|
|
|
|
|
- --- -- ---- ----- ... ---------- ---------- ---------|
|
|
|
|
|
|
|
|
|----------------------|
|
|
|
|
|
|
|_(2) 3 Spring stiffnesses (3E10.0)
|
|
|
|
|
|
|
|
|_(1) Non-tension element flag (I5)
|
|
|
|_Compulsory Data Record Keyword (A4)
|
|_Optional User Identifier (A2)
(1) For installed tethers, this should be left blank. For towed tethers, a '1' should be entered in Column
15.
(2) The values of the stiffnesses of the springs at the anchor end should be specified on the TSPA data
record. The stiffnesses of the springs at the vessel end should be specified on the TSPV data record.
For installed tethers, the spring stiffnesses are the inline/vertical stiffness and the two rotational stiffnesses
at the ends of the tether. Default values of 1.0E15 are used, if this data record is omitted. A default
value of 1.0E15 for the inline stiffness is used if the 1st field is left blank or a negative or zero values is
input. For the rotational fields, any value may be entered (except negative values which will be set to
zero).
For towed tethers, the stiffnesses are assumed to represent soft mooring line stiffness and are the three
stiffnesses in the translational directions. Note that the higher the stiffnesses input here, the smaller
the time steps will need to be in Data Category 16. This data record should always be present for towed
tethers.
2
5 7
11
16
- --- -- ---- -----
192
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) The number of modes to be output for the pre-processing eigensolution. The total number of modes
available:
= total number of degrees of freedom
= number of nodes * number of degrees of freedom/node
= number of nodes * 4
19.19.4. The TFAT Data Record - Tether Fatigue Parameters (Aqwa-Drift only)
This data record may be omitted, in which case the default values shown below will be used.
2
5 7
11
31
41
51
61
- --- -- ---...
---------- ---------- ---------- ---------|X|XXX| |TFAT|
|
|
|
|
|
- --- -- ---...
---------- ---------- ---------- ---------|
|
|
|
|
|
|
|
|--------------------------------|
|
|
|
|
|_(1) 4 Fatigue parameters( 4E10.0)
|
|
|
|_Compulsory Data Record Keyword (A4)
|
|_Optional User Identifier (A2)
=
=
where
R(i) = stress range (computed from rainflow count of time history stresses)
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
193
19.19.5. The TPSH Data Record - Tether Peak Stress Hours (Aqwa-Drift only)
This data record may be omitted, in which case the default value of 3 hours will be used.
2
5 7
11
16
- --- -- ---- ----|X|XXX| |TPSH|
|
- --- -- ---- ----|
|
|
|
|
|
|
|
|
|
|
|_(1) Number of Hours (I5)
|
|
|
|_Compulsory Data Record Keyword (A4)
|
|_Optional User Identifier (A2)
(1) The number of hours for which the expected peak stress is calculated
The extreme values of stress are based on the assumption that stress has a Rayleigh distribution. The
peak stress is given by
where
n = the number of cycles/hour, based on the stress time history calculation of the mean of the number
of positive and negative peaks
2
5 7
11
21
31
- --- -- ---- ---------- ---------- ---------|X|XXX| |TSLK|XXXXXXXXXX|XXXXXXXXXX|
|
- --- -- ---- ---------- ---------- ---------|
|
|
|
|
|
|
|
|_(1) Time (F10.0)
|
|
|
|_Compulsory Data Record Keyword (A4)
|
|_Optional User Identifier (A2)
194
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
2
5 7
11
- --- -- ---- ----|X|XXX| |TEGR|
|
- --- -- ---- ----|
|
|
|
|
|
|
|
|_(1) Tether group factor (F10.0)
|
|
|
|_Compulsory Data Record Keyword (A4)
|
|_Optional User Identifier (A2)
2
5 7
11
21
31
41
- --- -- ---- ---------- ---------- ---------- ---------|X|XXX| |TCAP|XXXXXXXXXX|XXXXXXXXXX|
|
|
- --- -- ---- ---------- ---------- ---------- ---------|
|
|
|
|
|
|
|
|
|
|
|_(2) Vessel Cap Area (E10.0)
|
|
|
|
|
|_(1) Anchor Cap Area (E10.0)
|
|
|
|_Compulsory Data Record Keyword (A4)
|
|_Optional User Identifier (A2)
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
195
2
5 7
11
21
31
41
- --- -- ---- ---------- ---------- ---------- ---------|X|XXX| |TIFL|XXXXXXXXXX|XXXXXXXXXX|
|
|
- --- -- ---- ---------- ---------- ---------- ---------|
|
|
|
|
|
|
|
|
|
|
|_(2) Density (E10.0)
|
|
|
|
|
|_(1) Pressure (E10.0)
|
|
|
|_Compulsory Data Record Keyword (A4)
|
|_Optional User Identifier (A2)
2
5 7
11
21
31
41
- --- -- ---- ---------- ---------- ---------- ---------|X|XXX| |TIFL|XXXXXXXXXX|XXXXXXXXXX|
|
|
- --- -- ---- ---------- ---------- ---------- ---------|
|
|
|
|
|
|
|
|
|
|
|_(2) Half Life (E10.0)
|
|
|
|
|
|_(1) Impact Factors (E10.0)
|
|
|
|_Compulsory Data Record Keyword (A4)
|
|_Optional User Identifier (A2)
196
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
2
5 7
11
21
31
- --- -- ---- ---------- ---------- ---------|X|XXX| |TLOW|XXXXXXXXXX|XXXXXXXXXX|
|
- --- -- ---- ---------- ---------- ---------|
|
|
|
|
|
|
|
|_(1) Lower Stop Position (E10.0)
|
|
|
|_Compulsory Data Record Keyword (A4)
|
|_Optional User Identifier (A2)
(1) The distance of the lower stop below the anchor. If the end of the tether is below this point, a
warning will be issued. Note that, if the lower stop distance is input as zero, the tether can never be
free hanging.
19.19.12. The TETH Data Record - Tether Vessel and Anchor/Trailing End Position
This data record should be input after the tether has been fully described i.e. all the previous data records
have been input.
For towed tethers this should be the *last* and *only* TETH data record in Data Category 14.
For installed tethers, where in general there is more than one tether, a complete tether description may
be duplicated by inputting a TETH data record immediately following another. In this case, the previous
tether will be duplicated at the positions specified by the structure/node numbers.
2
5 7
11
16
21
26
- --- -- ---- ----- ----- ----- ----|X|XXX| |TETH|
|
|
0|
|
- --- -- ---- ----- ----- ----- ----|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|_(4) Node 2 (I5)
|
|
|
|
|
|
|
|
|
|_(3) Fixed position specification (I5)
|
|
|
|
|
|
|
|_(2) Node 1 (I5)
|
|
|
|
|
|_(1) Structure Number (I5)
|
|
|
|_Compulsory Data Record Keyword (A4)
|
|_Optional User Identifier (A2)
(1) The number of the structure/vessel to which the tether is attached. This must correspond to one of
the structures defined in Data Category 2. If '1' is input, this will correspond to the structure defined in
Data Category 'ELM1'. If '2' is input, this will correspond to the structure defined in Data Category 'ELM2',
etc. Structure number '0' (i.e. a fixed node) is an illegal structure (in this position) and will produce an
error.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
197
5 7
11
21
31
---- -- ---- --------- --------- --------- ---XXXX| |TROC|XXXXXXXXX|XXXXXXXXX|XXXXXXXXX|XXXX
---- -- ---- --------- --------- --------- ---XXXX| |TROV|XXXXXXXXX|XXXXXXXXX|XXXXXXXXX|XXXX
---- -- ---- --------- --------- --------- ---XXXX| |TLAC|XXXXXXXXX|XXXXXXXXX|XXXXXXXXX|XXXX
---- -- ---- --------- --------- --------- ---XXXX| |TLAV|XXXXXXXXX|XXXXXXXXX|
0.100|XXXX
---- -- ---- --------- --------- --------- ---|
|
|
|
|
|_(3) Gap between tether and structure
|
|
|
|_(2) Mandatory Data Record Keyword
|
|_(1) Optional 2 letter user identifier
(2) Note that the rotational constraints, TROC and TROV, are rarely used, as this will cause large bending
moments at the attachment points. Weak/zero stiffness spring are normally used (see The TSPA/TSPV
Data Records - Tether Anchor and Vessel Springs (p. 192)).
(3) For a lateral constraint on the vessel, a gap can be specified, representing an opening which is wider
than the tether. It is assumed to be a frictionless circular gap in the structure, vertically below the
tether attachment point. If the total lateral movement relative to the center of the gap is greater than
the gap distance specified, the program will assume that the node at the gap is constrained laterally
by the structure.
As forces on the tether are by definition in the XY plane of the tether axis system (TLA), the reaction
on the structure must be at right angles to the TLA, i.e. for a vertical tether, the reactive force will be
in the horizontal plane of the FRA. For a sloping tether, i.e. when the TLP is offset, there will be a small
198
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
199
200
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
201
(1) The data record keyword indicates the corresponding structure number for the positions in columns
21-80; for example, enter POS1 for Structure 1, POS2 for Structure 2, .... PO10 for Structure 10.
(2)-(3) This is to indicate which spectrum and mooring configuration the structure position is defined
for, only needed for multiple spectra/mooring configuration cases.
(4)-(9) This is the position of the structure indicated by the data record keyword (1) at the start of the
analysis. The position is defined by three translations and three successive rotations. Angles are in degrees.
202
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) The data record keyword indicates the corresponding articulation number for the positions in columns
21-80; for example, enter REA1 for Articulation 1, REA2 for Articulation 2, .... RE10 for Articulation 10.
(2)-(3) This is to indicate which spectrum and mooring configuration the reaction is defined for, only
needed for multiple spectra/mooring configuration cases. If omitted the same initial reactions will be
used for all spectra and /mooring configurations.
(4)-(9) These are the initial reactions applied by the articulation to the 1st structure on the relevant
DCON data record (The DCON Data Record - Define Constraint Position (p. 129)) in Data Category 12.
The reactions are three forces and three moments. The reactions are applied in the global axis system.
The LAAR (Administration and Calculation Options for the Aqwa Suite (p. 252)) and LSAR (Administration
and Calculation Options for the Aqwa Suite (p. 252)) options, if used, apply only to output reactions.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
203
204
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
205
Initial Reactions
It is not possible to specify the initial articulation reactions in Aqwa-Drift. If the positions are correct
the articulation reactions will reach equilibrium in the 1st time step.
21.1.2. Analysis Type - Drift + Wave Frequency (WFRQ on JOB data record)
Starting Conditions - Defaults
The default starting slow position is the position of the structure as originally defined by the user in
Data Categories 1 to 4. The default starting slow velocity is zero in all 6 degrees of freedom.
The fast position and velocity are calculated by the program and added to the slow position and velocity
to form the total position and velocity. If any two structures are connected by an articulation this calculation is omitted for those structures, the fast position and velocity are zero and the total position and
velocity are the same as the slow position and velocity.
Starting Conditions with the RDEP option, but NONE in Data Category 15
If the RDEP option (Administration and Calculation Options for the Aqwa Suite (p. 252)) is used with a
filename given on the RESTART data record (e.g FILE01), then the position in the file FILE01.EQP will
override the default slow position. However, if the slow position is defined in Data Category 15 this will
in turn override the data in the .EQP file. The default starting slow velocity is zero in all 6 degrees of
freedom.
The fast position and velocity are calculated by the program and added to the slow position and velocity
to form the total position and velocity. If any two structures are connected by an articulation this calcu-
206
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Initial Reactions
It is not possible to specify the initial articulation reactions in Aqwa-Drift. If the positions are correct
the articulation reactions will reach equilibrium in the 1st time step.
207
(1) The data record keyword indicates the corresponding structure number for the positions in columns
21-80; for example, enter POS1 for Structure 1, POS2 for Structure 2, .... PO10 for Structure 10.
(2)-(7) This is the total position of the structure indicated by the data record keyword (1) at the start of
the time-history simulation. Angles are in degrees.
(1) The data record keyword indicates the corresponding structure number for the positions in columns
21-80; for example, enter VEL1 for Structure 1, VEL2 for Structure 2, .... VE10 for Structure 10.
(2)-(7) This is the total velocity of the structure indicated by the data record keyword (1) at the start of
the time-history simulation.
208
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
2
5 7
11
21
31
41
51
- --- -- ---- --------- --------- --------- --------- |X|
| |SLP1|XXXXXXXXX|
|
|
|...
- --- -- ---- --------- --------- --------- --------- | |
|
|
|
|
|
| |
|
|---------------------------| |
|
|
| |
|
|
| |
|
|_(2)-(7)6 Starting Slow Positions (default as
| |
|
defined in POS* data record in the same data category)(6F10.0)
| |
|
| |
|_(1)Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on Last Data Record in Data Category(A3)
(1) The data record keyword indicates the corresponding structure number for the positions in columns
21-80; for example, enter SLP1 for Structure 1, SLP2 for Structure 2, .... SP10 for Structure 10.
(2)-(7) This is the slow position of the structure indicated by the data record keyword (1) at the start of
the time-history simulation.
(1) The data record keyword indicates the corresponding structure number for the positions in columns
21-80; for example, enter SLV1 for Structure 1, SLV2 for Structure 2, .... SV10 for Structure 10.
(2)-(7) This is the slow velocity of the structure indicated by the data record keyword (1) at the start of
the time-history simulation.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
209
210
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Initial Reactions
It is not possible to specify the initial articulation reactions in Aqwa-Naut. If the positions are correct
the articulation reactions will reach equilibrium in the 1st time step.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
211
(1) The data record keyword indicates the corresponding structure number for the positions in columns
21-80; for example, enter POS1 for Structure 1, POS2 for Structure 2, .... PO10 for Structure 10.
(2)-(7) This is the position of the structure indicated by the data record keyword (1) at the start of the
time-history simulation. Angles are in degrees.
212
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) The data record keyword indicates the corresponding structure number for the positions in columns
21-80; for example, enter VEL1 for Structure 1, VEL2 for Structure 2, .... VE10 for Structure 10.
(2)-(7) This is the velocity of the structure indicated by the data record keyword (1) at the start of the
time-history simulation.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
213
214
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) This number (NT) governs the length of real time simulated, i.e. simulation time is given by (NT-1)
DT (see (2)). The maximum number of time-steps is 1,000,001.
(2) This important parameter (DT) governs the accuracy of the integration of the equations of motion.
This value must be small enough to enable an accurate representation of the highest frequency present
in the motion of a structure. Failure to do so will at best give an inaccurate simulation, and at worst
will cause divergence of the integration scheme and the program will abort.
For pure sinusoidal motion, this value should not be greater than 1/10th of the period of that motion.
For non-sinusoidal motion, the user should consider how accurately a set of discrete points (whose in-
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
215
Number of TimeSteps
Value of Time-Step
(Seconds)
Simulation Time
(Seconds)
2000
9,995
400
0.5
199.5
200
0.40
79.6
200
0.75
149.25
Aqwa-Naut
216
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
2
5 7
11
21
31
- --- -- ---- ----------- --------- --------|X|
| |HOTS|
|
|
|
- --- -- ---- ----------- --------- --------|
|
|
|_(1)Time Step Number at which Hot-Start run beings (I10)
|
|_Compulsory Data Record Keyword (A4)
(1) This is the time step number in a previous run at which the user wants the hot-start run to begin.
This is only needed when the start time is not defined in a subsequent TIME record; see (3) below.
When the HOTS data record is defined, the parameters for the TIME data record are now as follows:
2
5 7
11
21
31
- --- -- ---- ----------- --------- --------|X|
| |TIME|
|
|
|
- --- -- ---- ----------- --------- --------| |
|
|
|
_ (3) Start Time of the Hot-Start Run
| |
|
|
|
| |
|
|
|_(2) Value of Time-step (F10.0)
| |
|
|
| |
|
|_(1)Number of Time-steps for The Hot-Start Run(I10)
| |
|
| |
|
| |
|
| |
|_Compulsory Data Record Keyword (A4)
| |
| |_Optional User Identifier (A2)
|
|_Compulsory END on last data record in Data Category (A3)
(1) The number of time steps for this hot-start run only.
(2) This defines the time step size for this hot-start run.
(3) The start time of the hot-start run. This is only needed when the starting time-step number is not
defined on the HOTS card.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
217
218
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) If NONE is entered for the data category keyword, the program assumes default values.
(1) If the user does not specify the maximum number of iterations, then the default value of 100 is used.
219
(2) These are the maximum movements allowed per iteration, for each degree of freedom. Rotations
are in degrees. The default values are shown below.
Default values are: 2.00, 2.00, 0.50, 0.573, 0.573, 1.432
(G=9.81)
To take account of different unit systems Aqwa multiplies the default translation values only by G/9.81,
where G is the value of gravity input in Data Category 5. These defaults will be over-written by any
values input on the MMVE data record, but Aqwa does not apply any additional factors; i.e. the value
on the data record is the value that is used.
24.5. The MERR Data Record - Maximum Error Allowable for Equilibrium
2
5 7
11
21
31
41
51
- --- -- ---- ----- ----- --------- --------- --------- --------- --------- --------|X|
| |MERR|
|XXXXX|
|
|
|
|
|
|
- --- -- ---- ----- ----- --------- --------- --------- --------- --------- --------| |
|
|
|
| |
|
|
|_(2)Maximum Errors for System to
| |
|
|
be Considered in Equilibrium (6F10.0)
| |
|
|
| |
|
|_Structure Number (I5)
| |
|
| |
|
| |
|_Compulsory Data Record Keyword (A4)
| |
| |_Optional User Identifier (A2)
|
|_Compulsory END on last data record in Data Category (A3)
(2) These are the maximum errors allowed in the final equilibrium position, for each degree of freedom.
I.e. if the calculated movement is less than this, the structure is assumed to be in equilibrium. Rotations
are in degrees. The default values are shown below.
Default values are: 0.02, 0.02, 0.02, 0.057, 0.057, 0.143
(G=9.81)
To take account of different unit systems Aqwa multiplies the default translation values only by G/9.81,
where G is the value of gravity input in Data Category 5. These defaults will be over-written by any
values input on the MERR data record, but Aqwa does not apply any additional factors; i.e. the value
on the data record is the value that is used.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
2
5 7
11
16
21
31
41
51
61
71
- --- -- ---- ----- ----- --------- --------- --------- --------- --------- --------|X|
| |STRP|
|
|
|
|
|
|
|
|
- --- -- ---- ----- ----- --------- --------- --------- --------- --------- --------| |
|
|
|
|
|
|
| |
|
|
|
|
|
|_(5) The Angular Increment
| |
|
|
|
|
|
of the Hinge about
| |
|
|
|
|
|
Global Z (default=1 deg)
| |
|
|
|
|
|
| |
|
|
|
|
|_(4)The Angular Increment about
| |
|
|
|
|
the Hinge (default=1 deg)
| |
|
|
|
|
| |
|
|
|
|_(3) The Definition of Hinge Axis Vector (3F10.0)
| |
|
|
|
| |
|
|
|_(2) Number of Rotations of the Hinge (default=1) (I5)
| |
|
|
| |
|
|_(1)Number of Rotations about the Hinge (default=30) (I5)
| |
|
| |
|_Compulsory Data Record Keyword (A4)
| |
| |_Optional User Identifier (A2)
|
|_Compulsory END on last data record in Data Category (A3)
(1) This is the number of angular positions (typically roll angles) for which the forces are calculated.
(2) It is possible to calculate stability about a number of axes. Typically rotation about the X-axis will
give roll stability. The stability axis (the hinge) can itself be rotated about the global Z axis, so that pitch
stability, or stability for any intermediate position, can be obtained. This parameter is the number of
positions of the hinge. If only 1 hinge vector is required, the number of rotations of the hinge and angular increment of the hinge can be omitted.
(3) This is the initial direction of the hinge axis defined as a vector in the FRA. For example, for typical
roll stability this would be (1, 0, 0).
(4) This is the increment between the angular positions defined in (1).
(5) This is the increment between the hinge positions defined in (2).
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
221
222
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
223
Note
The results obtained for the new reference point are not written to the database backing
files.
224
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) The new inertia values are input in the following order:
Ixx, Ixy, Ixz, Iyy, Iyz, Izz
(1) The new coordinates of the center of gravity are given with respect to the original definition position
of the structure.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
225
226
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Note
The default value for the factor multiplying the slam coefficient is zero, as there are no general formulae documented in the literature (see The SLMM Data Record - Slam Multiplying
Factor (p. 230)).
The scale factor is commonly used where the effects of Morison drag on TUBE elements is considered
important (e.g. simulating tests at model scale). The multiplying factor is used for parametric studies
where the effect of the hydrodynamic coefficients is of particular interest.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
227
26.3. The SC1/ Data Record - Scale Factor for Model Test Simulation
If this data record is used the RNDD option (Administration and Calculation Options for the Aqwa
Suite (p. 252)) is also needed in the preliminary data category.
2
5 7
11
- --- -- ---- ----- ----|X|
| |SC1/|
|XXXXX|
- --- -- ---- ----- ----| |
|
|
| |
|
|_(1)Scale Factor (I5) (default implies constant drag coefficients)
| |
|
| |
|
| |
|_Compulsory Data Record Keyword (A4)
| |
| |_Optional User Identifier (A2)
|
|_Compulsory END on last data record in Data Category (A3)
(1) This scale factor applies to all elements on all structures but in the present version of the program
only affects Morison TUBE elements on those structures.
If this data record is omitted the program assumes that the Reynolds Number is sufficiently large for
the drag coefficient to be considered constant. If this is not a reasonable assumption then this data
record should be input with a value of unity.
The user should read the following sections before using this facility.
228
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) The structure number must correspond to one of the structures defined in Data Category 2. If '1' is
input, this will correspond to the structure defined in Data Category ELM1. If '2' is input, this will correspond to the structure defined in Data Category ELM2, etc (see note 1).
(2) This value is the multiplying factor for all drag/added mass coefficients for non-diffracting elements
on the structure specified (1).
Note
The multiplying factors relate to drag/added mass coefficients for the structure specified on
this data record only. The program does NOT multiply the values in the Geometric Properties
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
229
(1) The structure number must correspond to one of the structures defined in Data Category 2; in other
words, if '1' is input then this will correspond to the structure defined in Data Category ELM1. If '2' is
input then this will correspond to the structure defined in Data Category ELM2, etc. Note that this
multiplying factor relates to slam coefficient for the structure specified on this data record only.
(2) This value is the multiplying factor for all slam coefficients for MORISON elements on the structure
specified (1).
230
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Note
This data category is optional. Enter NONE if no additional information is required.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
231
27.3. The NODE Data Record - Nodal Position for Listing file Output
Maximum of the total number of NODE data records that may be input is 35.
2
5 7
11
16
21
26
- --- -- ---- ----- ----- ----- ----|X|
| |NODE|
|
|
|
|
- --- -- ---- ----- ----- ----- ----| |
|
|
|
|
|
| |
|
|
|
|
|_(4) Node Number (I5)
| |
|
|
|
|
| |
|
|
|
|_(3) Structure Number (I5)
| |
|
|
|
| |
|
|
|_(2) Node Number (I5)
| |
|
|
| |
|
|_(1) Structure Number (I5)
| |
|
| |
|_Compulsory Data Record Keyword (A4)
| |
| |_Optional User Identifier (A2)
|
|_Compulsory END on last data record in Data Category (A3)
(1) The structure number must correspond to one of the structures defined in Data Category 2. If '1' is
input then this will correspond to the structure defined in Data Category ELM1. If '2' is input then this
will correspond to the structure defined in Data Category ELM2, etc. As the NODE data record is a request
for output of the position/motion of the node number specified (2), structure number '0' (i.e. a fixed
node) is illegal in this field and will produce an error, since zero structure indicates a fixed node.
(2) This is the node number whose position/motions are requested during the analysis. The position of
this node on the structure (1) must be defined in Data Category 1.
Note that these motions are with respect to the position defined by (3) and (4).
(3)-(4) This structure number and its corresponding node number (4) define the reference point for the
positions/motions defined by parameters (1) and (2).
Both these fields may be left blank, in which case the program will assume that the output at the position defined by (1) and (2) is with respect to the origin of the Fixed Reference Axes, i.e. the ABSOLUTE
values.
If '0' is input as the structure number (3), together with a node number (4), the program will recognize
that this node number (4) references the fixed position as defined in Data Category 1 in the Fixed Reference Axis System (FRA). Note that a non-zero structure number (3) must be followed by a valid node
number (4) on that structure.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
The ALLM data record enables the user to expand or limit the information output associated with the
positions on the NODE data record. The default information, at these positions specified (where the
output format is freedom dependent), is displacement only in the X, Y, and Z translational freedoms.
Input of this data record, which is optional, has the following effect:
ALLM expands the output associated with the positions specified on all NODE data records to include
velocity and acceleration as well as displacement. When this information is not available this data record
has no function and is ignored.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
233
(1) The structure number must correspond to one of the structures defined in Data Category 2. If '1' is
input, this will correspond to the structure defined in Data Category ELM1. If '2' is input this will correspond to the structure defined in Data Category ELM2 etc.
(2) Not yet implemented.
(1) The timestep increment determines how often the full printout of positions and forces is printed to
the output listing file. For example, if the timestep increment was five then the printout would contain
information on each structure at timesteps 1,6,11,16 etc. This facility can be used to limit the size of
the printout to a manageable size for long simulations.
Note that recording of positions and forces in the graphics file is independent of the parameter on this
data record. See the GREV data record for more information.
234
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) PRNT causes the specified parameter to be added to the list of output. NOPR causes it to be removed
from the list.
(2) The structure number must correspond to one of the structures defined in Data Category 2. If '1' is
input, this will correspond to the structure defined in Data Category ELM1. If '2' is input, this will correspond to the structure defined in Data Category ELM2, etc.
(3) The parameter number refers to a type of output as shown in the list below. If the parameter number
is omitted all the parameters for the specified structure will be output (PRNT) or omitted (NOPR).
(4) For articulation reactions, the 3rd number on the data record is the number of the articulation (150) for which reactions are required.
The FROUDE-KRYLOV and DIFFRACTION forces are sub-divided differently for Aqwa-Drift and Naut, as
shown in the table below.
---------------------------------------------------------------------------------| PARAMETER (#)
|
AQWA-DRIFT
|
AQWA-NAUT
|
|---------------------------------------------------------|------------------------|
| DIFFRACTION (16) | Diffraction and Froude-Krylov force | Diffraction force on
|
|
| on diffracting panels
| diffracting panels
|
|
|
|
|
| FROUDE-KRYLOV (20)| Froude-Krylov force on Morison
| Froude-Krylov force on |
|
| elements only
| all elements
|
|
|
|
|
| WAVE INERTIA (23) | Diffraction force on
| Diffraction force on
|
|
| Morison elements
| Morison elements
|
----------------------------------------------------------------------------------
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
235
236
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
The ZRON/ZROF Data Records - Print Z Coordinate Relative to Wave Surface ON/OFF
|3O. ERROR PER TIMESTEP
|
|3O. ERROR PER TIMESTEP
|
|31. TOTAL REACTION FORCE |
|31. TOTAL REACTION FORCE |
|34. EXTERNAL FORCE
|
|
|
|35. RADIATION FORCE
|
|
|
|36. FLUID MOMENTUM
|
|
|
|38. FLUID GYROSCOPIC FORCE|
|
|
|39. ADD STRUCT STIFF FORCE|
|
|
|47. ARTICULATION REACTION |
SEE (3) ABOVE
|
|5O. TOTAL FORCE
|5O. TOTAL FORCE
|5O. TOTAL FORCE
|
------------------------------------------------------------------------------
This is the full list of output parameters for Aqwa-Librium, Drift and Naut. Further output can be requested
using additional data records in Data Category 18.
(1) The structure number indicates that the tensions in all cables, catenaries and hawsers attached to
this defined structure are to be printed in the listing file (how often is governed by the PREV data record)
and written to backing file for plotting in the AGS.
Note
The ZRON and ZROF data records act as switches during the input of the NODE data records
and have the function of switching the Z relative to the wave surface option on and off.
The ZRWS data record is global, and may be input anywhere in Data Category 18, these data
records therefore replace the existing ZRWS data record when only some nodes with Z relative
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
237
27.10. The ZRWS Data Record - Print Z Coordinate Relative to Wave Surface
This data record is used to output the relative vertical distance from the undiffracted wave surface to
a node defined in a NODE data record. All nodes defined in NODE data records in Data Category 18
will be affected except those for relative distances between two nodes.
2 5
7
11
- --- -- ---|X|
| |ZRWS|
- --- -- ---| |
|
| |
|
| |
|
| |
|_Data Record Keyword (A4)
| |
| |_Optional User Identifier (A2)
|
|_Compulsory END on last data record in Data Category (A3)
27.11. The PPRV Data Record - Print POS Every nth Timestep
2
5 7
11
- --- -- ---- ----|X|
| |PPRV|
|
- --- -- ---- ----| |
|
|
| |
|
|_(1)Timestep increment(I5)
| |
|
| |
|_Compulsory Data Record Keyword(A4)
| |
| |_Optional User Identifier(A2)
|
|_Compulsory END on last data record in Data Category (A3)
(1) The timestep increment determines how often the program will print positions and velocities etc
to the output *.POS file. For example, if the timestep increment was five then the printout would
contain information on each structure at timesteps 1,6,11,16 etc. Note that this data record only affects
the output in the *.POS file.
This facility can be used to limit the number of position/velocity output within the given period of time,
which is particularly useful in creating an animation (sequence file) without having to have one picture
for each time step.
27.12. The GREV Data Record - Graphics Output Every nth Timestep
2
5 7
11
- --- -- ---- ----|X|
| |GREV|
|
- --- -- ---- ----| |
|
|
238
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) The timestep increment determines how often the full graphics and plotting results are output. For
example, if the timestep increment was five then the backing files would contain information on each
structure at timesteps 1,6,11,16 etc.
This facility can be used to limit the size of the printout to a manageable size for long simulations.
In order to limit the size of the files the default maximum number of timesteps on the .PLT file is
10,000. If the GREV data record is not used the program will calculate a timestep increment such that
the number of timesteps does not exceed 10,000. For example, if a run has 20,000 timesteps and it is
required to plot every step, a data record "GREV 1" is required. If this is not included, the .PLT file will
contain results for every 2nd timestep.
The absolute maximum number of timesteps on the .PLT file is 100,000.
239
2
5 7
11
- --- -- ---- ----|X|
| |PMST|
|
- --- -- ---- ----| |
|
| |
|
| |
|
| |
|_Compulsory Data Record Keyword (A4)
| |
| |_Optional User Identifier (A2)
|
|_Compulsory END on last data record in Data Category (A3)
(1) If this number is '0', blank or a number greater than the number of spectra in the group, the main
sub-spectrum number of each wave spectral group is used.
240
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(1) Enter a non-zero integer 'N' where listing file time history output for tethers is required every 'N'
timestep.
(1) Enter a non-zero integer 'N' where graphics output and statistics post-processing for tethers is required
every 'N' timesteps.
27.16.3. The TSTS/TSTF Data Record - Tether Start/Finish Timesteps for Statistics
Input for Aqwa-Drift only.
If these data records are not input tether statistics post-processing will be on all records specified by
the TGRV Data Record.
2
5 7
11
- --- -- ---- ----|X|
| |TSTS|
|
- --- -- ---- ----|X|
| |TSTF|
|
- --- -- ---- ----| |
|
|
| |
|
|_(1) Start/Finish Timestep (I5)
| |
|
| |
|_Compulsory Data Record Keyword (A4)
| |
| |_Optional User Identifier (A2)
|
|_Compulsory END on last data record in data category (A3)
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
241
242
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 28: Element and Nodal Loads - ENLD (Data Category 21)
28.1. General Description
This data category is used to input the request for output of loads on Morison elements. At present,
this is only available for TUBE elements in Aqwa-Drift and Aqwa.
The loads are output in the form of loads at the nodes joining each element. The user should input
"RESTART 6 6" on the RESTART data record in Data Category 0. The *.RES file and *.POS file from a preceeding Drift or Naut run (stages 4 to 5) should be copied before this run can be performed.
There are two basic forms of nodal load output. The first is for space frames where all elements are assumed encastre-encastre and more than two elements can be joined at a single node. The second is
for riser-type structures, where riser geometry is assumed and, by definition, only two tube elements
can join at a node. In addition, there are rules concerning the element description, which the user must
be aware of in order to obtain the correct riser-type format (see The RISR Data Record - Nodal Load
Output for a Riser Structure (p. 244)).
The user should note that the post-processing Stage 6 graphics output overwrites any previous graphics
output. This means that if the user wishes to keep the output from Stage 5, then Stage 6 post-processing
must be carried out as a separate run.
28.3.The ISEL and LSEL Data Records - Element/Nodal Load Record Selection
These two data records perform the same function, however the LSEL data record allows input of timestep numbers with more than 5 digits.
2
5 7
11
16
- --- -- ---- ----- ----|X|
| |ISEL|
|
|
- --- -- ---- ----- ----| |
|
|
|
| |
|
|
|
| |
|
|
|
| |
|
|
|_(2) Terminal Record for output (I5)
| |
|
|
| |
|
|_(1) Initial Record for output (I5)
| |
|
| |
|_Compulsory Data Record Keyword (A4)
| |
| |_Optional User Identifier (A2)
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
243
(1) The initial record is the first record (starting time step number) that the user wishes to include in
the output and statistical post-processing of the nodal loads.
(2) The terminal record is the last record (finishing time step number) that the user wishes to include
in the output and statistical post-processing of the nodal loads. Thus the total number of records for
post-pressing is Terminal-Initial+1.
28.4. The RISR Data Record - Nodal Load Output for a Riser Structure
The riser data record specifies that inter-element forces are to be calculated assuming that a riser-type
structure has been defined.
When this card is input, the program makes the following assumptions:
Only tube elements are used to describe the riser in Data Category 2.
The first node of the first tube describing the riser is the node at the seabed end. The first node of all the
other elements must be the second node of the previous element.
If there is a connection to a fixed point this must be at the seabed end of the riser.
2
5 7
11
- --- -- ---- ----|X|
| |RISR|
|
- --- -- ---- ----| |
|
|
| |
|
|
| |
|
|
| |
|
|_(1) Structure Number (I5)
| |
|
| |
|_Compulsory Data Record Keyword (A4)
| |
| |_Optional User Identifier (A2)
|
|_Compulsory END on last data record in Data Category (A3)
(1) The structure number must correspond to one of the structures defined in Data Category 2. If '1' is
input, this will correspond to the structure defined in Data Category ELM1. If '2' is input, this will correspond to the structure defined in Data Category ELM2, etc.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
245
246
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
would run the data file ALTEST.DAT using the option NOWIND.
At present there are only two command line options available.
/NOWIND - No Windows
This option automatically closes all windows that may be opened during a run. This means that no response is required from the user at the end of a run.
/STD - Use Command File
This option instructs Aqwa to accept commands from an Aqwa command file. For compatibility with
previous versions of Aqwa this option (only) may be entered without the leading forward slash.
247
29.3.1. Printing of the Expanded Input Data List for each Data Category
By default, the printing of the expanded input data list for each data category is output to the listing
file at the end of the three data record image input Stages (1, 2, 4). Default output is therefore
For a run of Stages 1 to 5 - Output of expanded data list for Data Categories 1 to 18
For a run of Stages 2 to 5 - Output of expanded data list for Data Categories 6 to 18
For a run of Stages 3 to 5 - No output
For a run of Stages 4 to 5 - Output of expanded data list for Data Categories 9 to 18
For a run of Stage 5 - No output
248
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
249
250
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
251
2.9763E-01
5.3036E-02
6.3355E-01
1.0078E+00
-2.1092E+04
-2.4030E-01
-1.9725E+04
-7.8328E+03
Header Record#1:
Columns 1-10 - Reserved for Version Header (AQTF-1.0 :)
Columns 11-80 - Run Title
For each structure:
Logical Record#1:
Columns 1-2 - Structure Number.
Columns 3-4 - Number of Directions for this structure.
Columns 5-6 - Number of Frequencies for this structure.
Columns 9-80 - 6 Directions (degrees) PER LINE in field widths of 12.
NB:If more than 6 directions are input then columns 9-80 are used on the next line.
Logical Record#2:
252
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
253
254
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
(29.1)
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
255
Note
The RDDB option is only needed if the hydrodynamics file from the previous Aqwa-Line
run has been accidentally deleted.
256
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
257
258
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
259
260
TIME
TIMESTEP
STAGE(Int)
USER_FORCE(.....,TIME=0.0,TIMESTEP=1,STAGE=1
USER_FORCE(.....,TIME=1.0,TIMESTEP=1,STAGE=2
USER_FORCE(.....,TIME=1.0,TIMESTEP=2,STAGE=1
USER_FORCE(.....,TIME=2.0,TIMESTEP=2,STAGE=2
...)
...)
...)
...)
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Euler angles
Unit vectors in
LSA
261
Overall motion
The columns of q.p.r give the direction cosines of the LSA x,y,z axes in the FRA, e.g.
Note that the inverse of qpr is the same as its transpose, so use qprT to convert from FRA to LSA.
If the position of the CoG in the FRA is CG, then the position of x in the FRA is given by
262
Release 15.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.