CFD中文网

@李东岳 再次感谢李老师回复，我后续按照您的指导试一下，另外再查阅一下其他资料学习学习

@李东岳 感谢李老师回复。目前边界是静止的，弹簧一端连接静止边界，一端连接运动物块。从了解到的资料看，目前OF中的弹簧应该也只支持一端连接动物体，一端连接静物体。因为OF通过坐标确定连接点，我想知道如果连接点坐标在边界上，那么这个弹簧是否真的连在了边界上，能否和边界传递力？感觉也想不到什么能够验证的办法。

求助各位，目前of的动网格技术可不可以实现下图功能：在计算域顶端和计算域内的物块之间连接弹簧，同时给整个计算域施加简谐运动？
9dbaecc8-35ff-49cb-a874-3db65bdce331-1653113557(1).png
目前学习of五个月有余，从我的了解来看，of自带动网格库中的弹簧只支持一端连接固定点，一端连接运动的物体。想问问大家有没有做过类似的东西，能提供些思路，不动底层代码能够解决这个问题，不胜感激。

求助各位。我的OF算例计算时间步和输出时间步设置的都比较小（0.002s），导致计算结果有4000个时程文件，请问导入paraview后处理的时候，可以间隔几个文件导入一个么，4000帧的话处理起来太费时费空间了。

想请教各位，我对OF后处理中力系数forceCoeffs中的几个概念不太能理解，想请教一下各位。

forceCoeffs { type forceCoeffs; libs ( "libforces.so" ); writeControl timeStep; writeInterval 1; patches ( "floatingObject" ); rho rhoInf; log true; rhoInf 1000; liftDir (0 0 1); //升力方向 dragDir (0 0 1); //拖曳力方向 CofR (5 5 2); //用于计算力矩的点 pitchAxis (0 0 1); //pitch方向,围绕转动的轴 magUInf 0; //无限远处来流处的速度 lRef 2; //参考长度 Aref 8; //参考面积 }

首先这里对liftDir升力和dragDir拖曳力的定义是怎样的 ?举个例子，如果模拟的是一个物体在水中匀速运动，升力可否理解为浮力？那这里的拖曳力又是指什么呢？
还是这个例子，pitchAxis解释为围绕转动的轴，这又是指什么？另外如果在静水中运动，那magUInf可以设置为0吗？lRef和Aref是相对谁的参考呢，计算域吗？

@tens 感谢您的回复，我试着把这个头文件加入到.C文件中，但是对编译结果貌似并没有作用，报错是一样的。已经把我觉得可能能用到的头文件都放进去了:135:

#include "mylinearSpring.H" #include "addToRunTimeSelectionTable.H" #include "sixDoFRigidBodyMotion.H" // Foam classes #include "OFstream.H" #include "IFstream.H" #include "Vector.H" // std classes #include <stdio.h> #include <fstream> #include <iostream> #include <string.h> #include <stdlib.h> #include "pointPatchFields.H" #include "addToRunTimeSelectionTable.H" #include "Time.H" #include "polyMesh.H" #include "fixedValuePointPatchField.H" #include "fvCFD.H" #include "fvMesh.H" #include "volFields.H"

求助各位老师，我试着在sixDoFRigidBodyMotion库中建立一个新的mylinearSpring文件夹（类似于原linearSpring），意图改动其中的锚点anchor的位置。先试着让它z坐标加1，在编译时出错，错误信息如下。

sixDoFRigidBodyMotion/restraints/mylinearSpring/mylinearSpring.C:error: assignment of read-only location ‘((const Foam::sixDoFRigidBodyMotionRestraints::mylinearSpring*)this)->Foam::sixDoFRigidBodyMotionRestraints::mylinearSpring::anchor_.Foam::Vector<Cmpt>::z<double>()’ anchor_.z()=anchor_.z()+1;

提示在给只读位置幅值，是不是意味着在这里anchor是不能变动的量，那要怎样才能改动它呢？

@bestucan 感谢您的回复，文件不能编译的问题我目前已经解决。新的库编译好后我试着用新定义的弹簧mylinearspring连接我案例的两个物体，跑了一下case发现结果不尽人意，可能这个弹簧只能像报告中那样用在两个都是被动运动的浮体上，我这里有个主动运动的，就用不了了。可能还要后续对库进行另外的改动，鉴于现在刚学openfoam不久，之前没学过c++，所以源文件还读不懂，要先去补补c++的东西。再次感谢您的帮助:xinxin:

@axy-0 QQ图片20220216155209.png

@李东岳 感谢李老师的回复。刚刚我又试着修改了一下（代码在第一张图），但是这次修改过后运行结果仍然存在有问题，程序没有运行起来，请您指点。
![3699c109-7ea4-4bd2-a3d8-8143185f9510-_157$9_AFAWJ1{)8W35YL3.png
709379c4-f503-43a1-99b2-40762da25237-8L9F0U3P()586BOLA@{2NAP.png

@李东岳 您好，我刚刚解决了之前权限的问题，可是这次键入interFoam命令之后又出现了另一种错误，看提示好像是wmake出错？请问这又是什么原因，这个算例本身应该是没问题的。
![)K435ZJCL}TWILHXJP`O42.png