我是用的这个

https://blog.csdn.net/liuxiabing150/article/details/51920182/

没错吧。。

《Fluid Mechanics》 Frank M.White P239

《工程流体力学》 李玉柱

弹性力也是这么定义的

@bestucan

最后找到了解决办法，是边界条件的原因。本来上下圆柱连接的上下两个边界用的symmetry边界条件，改成cylic边界之后就能算了。。。具体原因还没研究出来

@fireztw

https://openfoamwiki.net/index.php/Parameter_Definitions_-_dynamicMotionSolverFvMesh

https://forum.cfdac.com/t/topic/660/2

这两个网站可以参考一下~~

@赵海盛 看起来没有问题的。。如果是自激振荡就是caculated，如果是强迫就要oscillatingdisplacement了。。。看边界条件是一个三维的圆形流场吗？。。

二维还是三维呢。。边界条件问题吧一般是

CPU型号： Intel(R) Xeon(R) Gold 6248R
系统：ubuntu 18.04
版本：openfoam-8

48 86.08
32 93.05
16 125.13
8 153.34
4 287.14
2 524.41
1 1013.26

正常了。。。

CPU型号： Intel(R) Xeon(R) Gold 6248R
系统：ubuntu 18.04
版本：openfoam-8

48 98.12
32 96.14
16 128.01
8 155.18
4 304.95
2 542.82
1 1068.49

没关超线程，物理核48 感觉不合理啊哈哈 明天再弄一次

尝试使用sixDoFRigidBodyMotion计算一个弹簧连接的三维圆柱在来流下的振动情况，但是结果发散。

在算到发散的前几步之前，库朗数开始急剧增加。

Courant Number mean: 0.097781 max: 12712.2

然后后处理发现本来应该随着刚体一起运动的近壁网格发生了变形，图如下（红色是库朗数大的地方，同样是网格变形的地方）

二维圆柱动的算例之前已经跑通过了，所以我先以为是constraints约束不够，后来加了几个约束后还是发散。 然后考虑到是不是icem没有关联好，cylinder边界重新关联了一下还是没有解决。

现在不知道问题出在哪里：有可能是感觉是innerDistance 0.1没有起到作用，但是使用slice切片观察Co发现圆柱中部的库朗数没有边界处这么大。所以不知道是不是约束没有设置合理。

[0_1607610864136_pointDisplacement](正在上传 100%)
附件中是相关的文件，不知道有没有人能帮忙解答一下。。网格太大了传不上来 ，，，
dynamicDictandpointdisplacement.zip

目前在用icoFoam算圆柱绕流，算例和 http://www.dyfluid.com/icoFoam.html 中设置基本相同。

现在想生成无量纲化平均应变张量$S_{ij}=\dfrac12\dfrac{k}{\epsilon}\left(\dfrac{\partial \bar{u_i}}{\partial \bar{x_j}}+\dfrac{\partial \bar{u_j}}{\partial \bar{x_i}}\right)$

OpenFOAM利用fieldAverage可以直接生成平均速度与雷诺应力$R_{ij}$。湍动能$k$即$tr(R_{ij})$

平均速度梯度可以后处理“postProcess -func “grad(UMean)””得到。

那耗散率在openfoam里面可以怎么求呢？$\epsilon=\nu \overline{\dfrac{\partial u_{i}^{\prime}}{\partial x_{k}} \dfrac{\partial u_{i}^{\prime}}{\partial x_{k}}}$ 看定义是可以通过改编fieldAverage求。但是试了一下没有成功。。请问一下能通过什么方法求得呢。