openFoam涡激振动
-
我在配置圆柱涡激振动算例的时候，运算了3步就算崩了，提示浮点数例外，我进去看了前几步的运算结果，发现网格变形并没有按照设定的那样沿着y方向改变，而是扭曲了。
mover { type motionSolver; libs ("libfvMeshMovers.so" "libsixDoFRigidBodyMotion.so"); motionSolver sixDoFRigidBodyMotion; sixDoFRigidBodyMotionCoeffs { patches ( CYLINDER ); innerDistance 5; outerDistance 25; mass 500; centreOfMass ( 0 0 0 ); momentOfInertia ( 1e-08 1e-08 1e-08 ); orientation ( 1 0 0 0 1 0 0 0 1 ); g ( 0 0 0 ); rho rhoInf; rhoInf 1000; report on; solver { type Newmark; gamma 0.5; beta 0.25; } accelerationRelaxation 0.4; constraints { yline { sixDoFRigidBodyMotionConstraint line; direction ( 0 1 0 ); } } restraints { ySpring { sixDoFRigidBodyMotionRestraint linearSpring; anchor ( 0 0 0 ); refAttachmentPt ( 0 0 0 ); stiffness 1000; damping 0; restLength 0; } } }这是我配置的动网格文件，请各位大佬帮我看一下，到底是哪儿出现问题了，为什么会出现这种情况呢，谢谢大家。
-
@李东岳 李老师，我还想请教您几个问题。我们在做圆柱涡激振动的时候，每一个时间步都会求解结构运动方程，意味着每一个时间步都需要对圆柱表面进行积分，来计算圆柱的受力。
我现在想复现下面这篇文献的内容，意味着圆柱表面存在不同的边界，即一个固定壁面边界，一个射流口边界，那么在计算的时候，程序默认是对这两个边界都积分呢？还是只对壁面边界进行积分。
第二张图是我们输出力系数的一个代码，里面可以定义监测的边界具体是哪个，这个代码会影响前面涡激振动力的计算么（意思就是，如果我这儿设置了只监测壁面边界，那么第一个问题就只对壁面边界进行积分）
第三个问题就是，如果我想把这两个边界分开，分别积分，我应该去哪儿修改，是dynamicMeshDict文件里面么？
-
@Joker 在 openFoam涡激振动 中说：
@李东岳 李老师，我还想请教您几个问题。我们在做圆柱涡激振动的时候，每一个时间步都会求解结构运动方程，意味着每一个时间步都需要对圆柱表面进行积分，来计算圆柱的受力。
我现在想复现下面这篇文献的内容，意味着圆柱表面存在不同的边界，即一个固定壁面边界，一个射流口边界，那么在计算的时候，程序默认是对这两个边界都积分呢？还是只对壁面边界进行积分。
第二张图是我们输出力系数的一个代码，里面可以定义监测的边界具体是哪个，这个代码会影响前面涡激振动力的计算么（意思就是，如果我这儿设置了只监测壁面边界，那么第一个问题就只对壁面边界进行积分）
第三个问题就是，如果我想把这两个边界分开，分别积分，我应该去哪儿修改，是dynamicMeshDict文件里面么？
朋友，这篇文章能给个标题吗？最近也打算做做这块内容，想看看这个文献
