用tecplot导出圆柱绕流中的表面平均压力系数

需要软件：tecplot15，matlab

需要语音讲解请转到我录的b站视频，ID：带一司机，还上传了其他前处理（画网格），后处理视频等等，后续会继续更新，欢迎关注：）

https://www.bilibili.com/video/BV1Mr4y1w71c/

1. 导入算好的云图，这里是双圆柱绕流的算例，Re=2e6，使用DES湍流模型，经过一段时间计算之后，平均压力场如下：

这里平均压力场（pMean）是用OpenFOAM的后处理工具导出的，这里不讨论其他CFD软件。

1.5 计算压力系数

点击DATA=>Alter=>Specify Equation，计算表面压强系数，公式是

在窗口上输入：

{Cp}={pMean}*2/1.225/44/44

数值根据自己的算例进行修正。

2. 转到3D的视角

3. 在左侧的工具栏里，勾选Slice，然后点击detail，主要看划红线的位置，Slice location选择对应你要切的平面，我这里选择Z平面。Slice through 那里要选Surface Zones，不然不能与圆柱表面相切。（选Volume Zones的话会截取整一个Z平面）

4. 关掉Details，在左边工具栏点Zone Style，把除了你想要截的那个圆柱都隐藏掉。如图，只剩下上游的圆柱。

5. 由于我的算例是2D拉伸的，不需要截特定的Z方向的位置（都一样），如果是3D算例，可以在3.那个Details窗口进行设置，或者手动拖拽。现在，Slice算是截好了，我们将这一串数据提取出来：

点击上方Data=>Extract=>Current Slice，再点Extract确定。这里不好截图。

6. 提取数据点。点击DATA=>Spreadsheet，在上方Zone选择刚刚截好的Slice: Z=0的Zone:

7. 用Ctrl选择X、Y、Cp三列，复制到txt文件中。

8. 用Matlab打开进行处理，也可用其他软件，这里主要介绍方法。将刚刚的txt文件以数值矩阵的格式导入到工作区：

9. 对于圆柱（圆形），只要有了xy坐标就可以改成角度，这里用matlab自带的atan2函数，编辑函数：

function fun(b,x0,y0)
a(:,1)=atan2(b(:,2)-y0,b(:,1)-x0);
a=a*180/3.14;
a(:,2)=b(:,3);
plot(a(:,1),a(:,2),'d')
end

其中，b为刚刚复制下来的3列的数据，x0为圆心x坐标，y0为圆心y坐标。因为我要算两个圆柱，所以在这里圆心位置也是必要的输入。如果是其他外形，则需要更复杂的处理手法，这里暂时不谈。

输入fun(Cp,0,0)，对应上游的圆柱，得到下图：

同样方法得到下游圆柱表面压强分布系数图：

结束。