Skip to content

Commit 36ee755

Browse files
committed
Added new rod outputs to inputs.rst
1 parent ebb5f92 commit 36ee755

File tree

2 files changed

+19
-11
lines changed

2 files changed

+19
-11
lines changed

docs/inputs.rst

Lines changed: 9 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -429,7 +429,15 @@ outputs are wanted. Eight output properties are currently possible:
429429
- c – internal damping force at each segment
430430
- s – strain of each segment
431431
- d – rate of strain of each segment
432-
432+
- f – net node force
433+
- W – weight at each node
434+
- B – buoyancy force at each node
435+
- P – dynamic pressure at each node
436+
- X – transverse drag force at each node
437+
- Y – tangential drag force at each node
438+
- A – transverse inertia force at each node
439+
- a – axial fluid inertia force at each node
440+
- b – bottom contact force
433441
For example, outputting node positions and segment tensions could be achieved by writing “pt” for
434442
this last column. These outputs will go to a dedicated output file for each rod. For sending
435443
values to the global output file, use the Outputs section instead.

source/Rod.cpp

Lines changed: 10 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -280,13 +280,13 @@ Rod::openoutput()
280280
*outfile << "Node" << i << "Box \t Node" << i << "Boy \t Node"
281281
<< i << "Boz \t ";
282282
}
283-
// output Dynamic Pressure
283+
// output transverse drag
284284
if (channels.find("X") != string::npos) {
285285
for (unsigned int i = 0; i <= N; i++)
286286
*outfile << "Node" << i << "Dpx \t Node" << i << "Dpy \t Node"
287287
<< i << "Dpz \t ";
288288
}
289-
// output Tangential drag
289+
// output tangential drag
290290
if (channels.find("Y") != string::npos) {
291291
for (unsigned int i = 0; i <= N; i++)
292292
*outfile << "Node" << i << "Dqx \t Node" << i << "Dqy \t Node"
@@ -305,7 +305,7 @@ Rod::openoutput()
305305
<< i << "AqZ \t ";
306306
}
307307
// output dynamic pressure
308-
if (channels.find("Pd") != string::npos) {
308+
if (channels.find("P") != string::npos) {
309309
for (unsigned int i = 0; i <= N; i++)
310310
*outfile << "Node" << i << "Pdx \t Node" << i << "Pdy \t Node"
311311
<< i << "Pdz \t ";
@@ -331,12 +331,12 @@ Rod::openoutput()
331331
for (unsigned int i = 0; i <= 3 * N + 2; i++)
332332
*outfile << "(m) \t";
333333
}
334-
// output velocities?
334+
// output velocities
335335
if (channels.find("v") != string::npos) {
336336
for (unsigned int i = 0; i <= 3 * N + 2; i++)
337337
*outfile << "(m/s) \t";
338338
}
339-
// output net node force?
339+
// output net node force
340340
if (channels.find("f") != string::npos) {
341341
for (unsigned int i = 0; i <= 3 * N + 2; i++)
342342
*outfile << "(N) \t";
@@ -351,13 +351,13 @@ Rod::openoutput()
351351
for (unsigned int i = 0; i <= 3 * N + 2; i++)
352352
*outfile << "(N) \t";
353353
}
354-
// output Dynamic Pressure
355-
if (channels.find("Y") != string::npos) {
354+
// output transverse drag
355+
if (channels.find("X") != string::npos) {
356356
for (unsigned int i = 0; i <= 3 * N + 2; i++)
357-
*outfile << "(Pa) \t";
357+
*outfile << "(N) \t";
358358
}
359-
// output Tangential drag
360-
if (channels.find("X") != string::npos) {
359+
// output tangential drag
360+
if (channels.find("Y") != string::npos) {
361361
for (unsigned int i = 0; i <= 3 * N + 2; i++)
362362
*outfile << "(N) \t";
363363
}

0 commit comments

Comments
 (0)