小白求助——floatingobject算例浮块突然失去浮体性质

[bestucan]

我抽空看了一下，然后，没看懂...

报错提示

you are probably trying to solv3e for a field with a default boundary condition.

你可以从这个 搜索链接 上看cfd online上关于该报错的搜索结果。我没做过这个，不太了解。

[AXY 0]

@bestucan 不好意思，可能我的表述有所欠缺。我试着换一种问法，请问一个动网格文件dynamicMeshDict里可以计算计算域中两处网格的运动吗（我目前涉及到的只有计算域中一个东西运动）？就像上面贴出的代码，那是我的动网格文件里的，代码上面的部分是控制水中物体上浮，下面的是控制水面以上的物体运动。我分别运行是没有问题的，但是合起来就报错了，我怀疑是放到一起导致的。不知这次有没有将问题描述清楚，求解答，感谢老师。

[李东岳]

看一下你cellDisplacement的边界条件。这个需要时calculated，你可能给了边界条件

[bestucan]

@axy-0 我说的看不懂不是你的描述看不懂，是我看不懂这个领域的问题和代码

不过你描述的更详细，我不懂也能找到一些相关资料。

看这个的第七页底部的 motionSolverLibs 配置。两个浮体求解器。你的只有一个。但估计是求解方法不一样。

还有这个，应该是用一个求解器求解两个浮体。

看哪个对你有用。

[AXY 0]

@李东岳 感谢李老师回复，我零文件下的pointdisplacement代码如下，水中物体的边界类型我一开始设置为calculated，就出现了上面的报错提示，报错后我设置为了固定数值fixedvalue，就可以运行了，但与此带来的就是水中物体失去浮体性质。

dimensions      [0 1 0 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{    

    stationaryWalls
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }

    atmosphere
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }

    floatingObject
    {
        type         	calculated;   			
        value           uniform (0 0 0);
    }
  
	Object
    {
        type            oscillatingDisplacement;
		omega           6.28;
		amplitude       (0 0 0.5);
		value           uniform (0 0 0);
    }

}

[AXY 0]

@bestucan 感谢这位老师的帮助，今天有点小忙，先马一下您的资料

[李东岳]

因为你设置为固定值了，所以就失去浮动特征了。你要设置成calculated

[AXY 0]

@bestucan 您好，现在我解决了之前的问题，两个物体可以保持自身的运动状态了。所以又想继续在两个物体之间添加一个弹簧，因为我现在还看不懂c++的代码，因此只能依照您给的双附体的那份资料依葫芦画瓢，在修改sixdofrigidbodymotion代码后进行编译，在编译过程中出现了点问题，但也没有报error，请您指点以下问题所在。

[AXY 0]

@李东岳 感谢老师回复，现在我解决了之前的问题，在进一步往下做。

[bestucan]

myLinearSpring.C 无法编译。

想把它整的能编译就比较麻烦了， 这不是 C++ 而是 makefile 的知识。本来就一个 makefile 的话也不算特别麻烦。 OpenFOAM 自己封装了个 wmake，把 makefile 分放到 Make 文件夹里的 file 和 option 了。

依赖的源文件和库都在那里写着。myLinearSpring.C.dep 就是编译源文件依赖的文件。要么是 myLinearSpring.C 的位置不对，要么是编译路径没配置对。

参考这个，不过这个是 makefile 不是 wmake

[AXY 0]

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