Skip to content

Introducing SubDyn DirCos =-1 and sanitization of input files #95

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
ebranlard merged 3 commits into from
Apr 14, 2023
Merged

Introducing SubDyn DirCos =-1 and sanitization of input files #95

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
5a2db9b
SD: update of input files for COSMID and cross-section comments
ebranlard Apr 14, 2023
97acecc
SD: adding test for DirCos
ebranlard Apr 14, 2023
665a584
SD: adding arbitrary section test (only one with NDIV=1), but not par…
ebranlard Apr 14, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
50 changes: 25 additions & 25 deletions ...-codes/openfast/5MW_OC3Mnpl_DLL_WTurb_WavesIrr/NRELOffshrBsline5MW_OC3Monopile_SubDyn.dat
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,29 +1,29 @@
----------- SubDyn MultiMember Support Structure Input File ------------
----------- SubDyn MultiMember Support Structure Input File ---------------------------
OC3 Monopile configuration (pile only).
-------------------------- SIMULATION CONTROL ---------------------------------
-------------------------- SIMULATION CONTROL -----------------------------------------
False Echo - Echo input data to "<rootname>.SD.ech" (flag)
"DEFAULT" SDdeltaT - Local Integration Step. If "default", the glue-code integration step will be used.
3 IntMethod - Integration Method [1/2/3/4 = RK4/AB4/ABM4/AM2].
True SttcSolve - Solve dynamics about static equilibrium point
False GuyanLoadCorrection - Include extra moment from lever arm at interface and rotate FEM for floating.
-------------------- FEA and CRAIG-BAMPTON PARAMETERS---------------------------
-------------------- FEA and CRAIG-BAMPTON PARAMETERS ---------------------------------
3 FEMMod - FEM switch: element model in the FEM. [1= Euler-Bernoulli(E-B); 2=Tapered E-B (unavailable); 3= 2-node Timoshenko; 4= 2-node tapered Timoshenko (unavailable)]
3 NDiv - Number of sub-elements per member
True CBMod - [T/F] If True perform C-B reduction, else full FEM dofs will be retained. If True, select Nmodes to retain in C-B reduced system.
0 Nmodes - Number of internal modes to retain (ignored if CBMod=False). If Nmodes=0 --> Guyan Reduction.
1 JDampings - Damping Ratios for each retained mode (% of critical) If Nmodes>0, list Nmodes structural damping ratios for each retained mode (% of critical), or a single damping ratio to be applied to all retained modes. (last entered value will be used for all remaining modes).
2 GuyanDampMod - Guyan damping {0=none, 1=Rayleigh Damping, 2=user specified 6x6 matrix}
0.000, 0.000 RayleighDamp - Mass and stiffness proportional damping coefficients (Rayleigh Damping) [only if GuyanDampMod=1]
0.000, 0.000 RayleighDamp - Mass and stiffness proportional damping coefficients (Rayleigh Damping) [only if GuyanDampMod=1]
6 GuyanDampSize - Guyan damping matrix (6x6) [only if GuyanDampMod=2]
0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00
0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00
0.0000e+00 0.0000e+00 1.4625e+06 0.0000e+00 0.0000e+00 0.0000e+00
0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00
0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00
0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00
---- STRUCTURE JOINTS: joints connect structure members (~Hydrodyn Input File)---
---- STRUCTURE JOINTS: joints connect structure members (~Hydrodyn Input File) --------
4 NJoints - Number of joints (-)
JointID JointXss JointYss JointZss JointType JointDirX JointDirY JointDirZ JointStiff
JointID JointXss JointYss JointZss JointType JointDirX JointDirY JointDirZ JointStiff ![Coordinates of Member joints in SS-Coordinate System][JointType={1:cantilever, 2:universal joint, 3:revolute joint, 4:spherical joint}]
(-) (m) (m) (m) (-) (-) (-) (-) (Nm/rad)
1 0.00000 0.00000 -20.00010 1 0.0 0.0 0.0 0.0
2 0.00000 0.00000 -10.00000 1 0.0 0.0 0.0 0.0
Expand All @@ -39,47 +39,47 @@ RJointID RctTDXss RctTDYss RctTDZss RctRDXss RctRDYss RctRDZss
IJointID ItfTDXss ItfTDYss ItfTDZss ItfRDXss ItfRDYss ItfRDZss ![Global Coordinate System]
(-) (flag) (flag) (flag) (flag) (flag) (flag)
4 1 1 1 1 1 1
----------------------------------- MEMBERS --------------------------------------
3 NMembers - Number of frame members
MemberID MJointID1 MJointID2 MPropSetID1 MPropSetID2 MType COSMID
----------------------------------- MEMBERS -------------------------------------------
3 NMembers - Number of members (-)
MemberID MJointID1 MJointID2 MPropSetID1 MPropSetID2 MType COSMID ![MType={1:beam circ., 2:cable, 3:rigid, 4:beam arb.}. COMSID={-1:none}]
(-) (-) (-) (-) (-) (-) (-)
1 1 2 1 1 1
2 2 3 1 1 1
3 3 4 1 1 1
------------------ MEMBER X-SECTION PROPERTY data 1/2 [isotropic material for now: use this table for circular-tubular elements] ------------------------
4 NPropSets - Number of structurally unique x-sections (i.e. how many groups of X-sectional properties are utilized throughout all of the members)
1 1 2 1 1 1 -1
2 2 3 1 1 1 -1
3 3 4 1 1 1 -1
------------------ CIRCULAR BEAM CROSS-SECTION PROPERTIES -----------------------------
4 NPropSets - Number of structurally unique cross-sections
PropSetID YoungE ShearG MatDens XsecD XsecT
(-) (N/m2) (N/m2) (kg/m3) (m) (m)
1 2.10000e+11 8.08000e+10 8500.00 6.000000 0.060000
2 2.10000e+11 8.08000e+10 8500.00 6.000000 0.027000
3 2.10000e+11 8.08000e+10 8500.00 3.870000 0.019000
4 2.10000e+14 8.08000e+14 1.00 1.000000 0.010000
------------------ MEMBER X-SECTION PROPERTY data 2/2 [isotropic material for now: use this table if any section other than circular, however provide COSM(i,j) below] ------------------------
0 NXPropSets - Number of structurally unique non-circular x-sections (if 0 the following table is ignored)
----------------- ARBITRARY BEAM CROSS-SECTION PROPERTIES -----------------------------
0 NXPropSets - Number of structurally unique non-circular cross-sections (if 0 the following table is ignored)
PropSetID YoungE ShearG MatDens XsecA XsecAsx XsecAsy XsecJxx XsecJyy XsecJ0
(-) (N/m2) (N/m2) (kg/m3) (m2) (m2) (m2) (m4) (m4) (m4)
-------------------------- CABLE PROPERTIES -------------------------------------
-------------------------- CABLE PROPERTIES -------------------------------------------
0 NCablePropSets - Number of cable cable properties
PropSetID EA MatDens T0 CtrlChannel
(-) (N) (kg/m) (N) (-)
----------------------- RIGID LINK PROPERTIES ------------------------------------
(-) (N) (kg/m) (N) (-)
----------------------- RIGID LINK PROPERTIES -----------------------------------------
0 NRigidPropSets - Number of rigid link properties
PropSetID MatDens
(-) (kg/m)
---------------------- MEMBER COSINE MATRICES COSM(i,j) ------------------------
(-) (kg/m)
---------------------- MEMBER COSINE MATRICES COSM(i,j) -------------------------------
0 NCOSMs - Number of unique cosine matrices (i.e., of unique member alignments including principal axis rotations); ignored if NXPropSets=0 or 9999 in any element below
COSMID COSM11 COSM12 COSM13 COSM21 COSM22 COSM23 COSM31 COSM32 COSM33
(-) (-) (-) (-) (-) (-) (-) (-) (-) (-)
------------------------ JOINT ADDITIONAL CONCENTRATED MASSES--------------------------
0 NCmass - Number of joints with concentrated masses; Global Coordinate System
CMJointID JMass JMXX JMYY JMZZ JMXY JMXZ JMYZ MCGX MCGY MCGZ
CMJointID JMass JMXX JMYY JMZZ JMXY JMXZ JMYZ MCGX MCGY MCGZ
(-) (kg) (kg*m^2) (kg*m^2) (kg*m^2) (kg*m^2) (kg*m^2) (kg*m^2) (m) (m) (m)
---------------------------- OUTPUT: SUMMARY & OUTFILE ------------------------------
True SumPrint - Output a Summary File (flag).It contains: matrices K,M and C-B reduced M_BB, M-BM, K_BB, K_MM(OMG^2), PHI_R, PHI_L. It can also contain COSMs if requested.
---------------------------- OUTPUT: SUMMARY & OUTFILE --------------------------------
True SumPrint - Output a Summary File (flag)
0 OutCBModes - Output Guyan and Craig-Bampton modes {0: No output, 1: JSON output}, (flag)
0 OutFEMModes - Output first 30 FEM modes {0: No output, 1: JSON output} (flag)
False OutCOSM - Output cosine matrices with the selected output member forces (flag)
False OutAll - [T/F] Output all members' end forces
False OutAll - [T/F] Output all members' end forces
2 OutSwtch - [1/2/3] Output requested channels to: 1=<rootname>.SD.out; 2=<rootname>.out (generated by FAST); 3=both files.
True TabDelim - Generate a tab-delimited output in the <rootname>.SD.out file
1 OutDec - Decimation of output in the <rootname>.SD.out file
Expand Down
Loading