论文的直径是0.5 高度是是1 按照这个刚体的惯性矩公式
我带进去计算和论文的惯性矩差很远
请问是我哪里错了 
还有就是他的升力系数中的pitchaxis不是固定是z轴吗 cofr在圆心，我也不明白为什么他设置在这2点
还有按着这个图设置弹簧，那个长度不是应该超出计算域了吗
付上这篇论文2015 Tawekal_Jessica master VIV.zip (3.8 MB)
Hi. 我看了下，
- 他的 MomentOfInertia ，我觉得是计算上有问题。MomentOfInertial 是用来计算转矩的，他的论文只有一个自由度，所以这个数据可能是任意的。
- pitchaxis 和 CoR 在这里没有使用，他论文也不需要计算 moment 系数，因此 CoR 可能是任意的
- anchor 和计算域无关，计算弹簧的伸长量是 mag(restraintPosition - anchor) - restLength
大佬你好、任意选取pitchaxis和cor不会影响升阻力系数的计算吗，还有随意选择anchor也不会影响圆柱在y上的振动吗… 还有大佬如果我想实现2个圆柱都这样单自由度的抖动，是只能用重叠网格吗，还是可以像受迫振动那样直接在dymnamicdict里面设置muitisoliodmovesholve那样直接在里面设置两个圆柱，我发现识别不了2个圆柱，如果要用重叠网格有没资料 
不会影响升阻力的计算，anchor 也不会影响，但要注意初始长度与 (refAttachmentPt -anchor) 一致。
两个圆柱单自由度振动，可以采用这种方式。
你说的那种在 dynamicDict 里面设置的， 是比较新的版本，也可以使用。
大神你好，我不是很明白你下面那个意思，你的意思是我可以用老版本的pimpledymfoam直接在dynamicDict设置两个圆柱都这样六自由度吗，请问要怎么设置才能2个圆柱都这样六自由度网格控制。
dynamicFvMesh dynamicMotionSolverFvMesh;
motionSolverLibs ("libsixDoFRigidBodyMotion.so");
solver sixDoFRigidBodyMotion;
sixDoFRigidBodyMotionCoeffs
{
patches (CYLINDER);
innerDistance 1;
outerDistance 5;
centreOfMass (0 0 0);
momentOfInertia (5.984 5.984 11.632);
mass 11.961;
rho rhoInf;
rhoInf 1000;
report on;
reportToFile on;
accelerationRelaxation 0.8;
solver
{
type Newmark;
}
constraints
{
yLine
{
sixDoFRigidBodyMotionConstraint line;
centreOfRotation (0 0 0);
direction (0 1 0);
}
rotation
{
sixDoFRigidBodyMotionConstraint orientation;
}
}
restraints
{
translationLinearSpring
{
sixDoFRigidBodyMotionRestraint linearSpring;
anchor (0 -0.0445 0);
refAttachmentPt (0 -0.0445 0);
stiffness 400;
damping 0.12;
restLength 0;
}
}
value uniform (0 0 0);
}
}
还是只能用新版本的overpimpledymfoam，那个我不会用，他好像还要两套网格，能帮我修改下这个dynamicDict给我参考下吗，
你好，我最近也在研究六自由度，想请教一下这里的rhoInf是指物体的密度吗？我在foam-extend自带的wingMotion算例里面看到它设置的rhoInf是1，这是不是太小了？
不是 应该是流体密度，没有设置小…因为那在空气中，用的应该是空气的密度 可以加个QQ交流下…
对的，rhoInf 是流体的密度
多谢大佬，我还有一个问题没想明白，像floatingobject这样的两相流算例中流体的密度应该选什么呢，我看自带算例里面也是1