Tecplot后处理扩展工具:自定义函数

释放双眼,带上耳机,听听看~!

最近回家练字多一些,所以断更了几天。写得一手好字也是妥妥的加分项,可惜这次买的字帖差点意思,是那种凹陷设计的,有一说一不是很舒服写起来,还是原来那种每一页附带1张白色透明纸的用起来舒服。

官方帮助文档链接点击跳转

结合最近一些群友遇到的问题,以及我自己在利用Tecplot后处理时遇到的问题,参数不全。就比如在处理CFX计算结果时,要看速度云图却只有X、Y、Z三个分量值,所以这个时候我们还需要学习一些Specify Equations(变量),比如常见的Pressure、Temperature等都是已经编辑好的公式供我们使用,没有的要自己在Tecplot中扩展!
Tecplot后处理扩展工具:自定义函数

依次点击软件顶部的Data-Alter-Specify Equations进入公式编辑器,顶部的空白区域用于输入公式,左下角的Zones to Alter用于选择需要参与计算的区域,比如我想计算叶片上的速度分布,那么在Zone区域中选择关于叶片的边界即可,通常这里默认即可!

另外可以看到每一个边界前面都有一个数字,这个值代表了这个边界,所以在输入公式时只需要前面这个数字作为代码即可!同样我们点击Data Set Info按钮也可以看到现有的变量前面都有数字的,后面会演示用途!
Tecplot后处理扩展工具:自定义函数
Tecplot后处理扩展工具:自定义函数

那么首先来编辑一个刚需的求解速度场的函数:Velocity

因为速度是矢量,不仅有大小还有方向,所以三维中有3个X、Y、Z三个方向的速度,这个值在默认函数中已经有了,可以看上图中的47、48、49就是,公式如下:

{velocity}=sqrt(v47**2+v48**2+v49**2)

其中“velocity”是你自己定义的一个变量名而已,可以是V1或者其他词组,目测不能是中文

sqrt()是开根号的意思,v47**2的意思是计算v42也就是X方向速度的平方值,然后将3个方向上的速度平方值之和开根号就能计算出合速度。

注意格式,变量要用{}包裹,47、48、49这三个值来自Tecplot,对应变量前的ID,可能每个人每个案例这个值都不一样,所以要搭配自己的案例设置。

接下来对比一下CFD-POST和Tecplot做出的压力云图和速度云图的对比:(左侧为Tecplot右侧为CFD-POST)
Tecplot后处理扩展工具:自定义函数Tecplot后处理扩展工具:自定义函数

Tecplot的公式功能还有很多用途,我们参照帮助文档中的说明来学习一波:
Tecplot后处理扩展工具:自定义函数

其中+、-、*、/的用法自不必说,加减乘除,**指平方,上面我们已经在实际案例中用到了

包括一些三角函数的应用,其中角度都需要指定为弧度rad才可以
Tecplot后处理扩展工具:自定义函数条件表达式我们应该很少用到,如果是什么那么结果取什么,可以用于流场中筛选数据,格式为:IF(P,T,F)

如果谓词表达式 P 为真(具有非零值),则返回表达式 T 的值,否则返回表达式 F 的值。无论 P 的值如何,T 和 F 都被完全计算;没有“短路”。

配合一些等于、不等于、大于等于以及布尔运算等功能会使得IF条件筛选更强大

例如:{A} = IF(X==2.0, 3, 4)

意思是:如果 X 等于 2.0,则A赋值为 3,否则赋值为 4
Tecplot后处理扩展工具:自定义函数

求导也可以,看着就头疼,用不到,略
Tecplot后处理扩展工具:自定义函数

下面是官方的公式举例,已通过谷歌翻译,凑活看吧
1、在以下等式中,V1(数据文件中定义的第一个变量)被 V1 现有值的两倍半替换:

V1 = 2.5*V1

2、以下等式将名为 Density 的变量的值设置为 205。如果名为 Density 的变量不存在,则会创建一个新变量。

{Density} = 205

3、在下一个等式中,Y(分配给 Y 轴的变量)的值被 X(分配给 X 轴的变量)的值平方的负数替换:

Y = -X**2

以下等式将 V3 的值替换为 V2 的值,四舍五入到最接近的整数。如果数据集中当前只有两个变量,则创建一个新变量。

V3 = round(V2)

4、在以下等式中,数据集中第四个变量的值被第三个变量的值的对数(以 10 为底)替换。

V4 = ALOG10(V3)

5、假设第三个和第四个变量是速度的 X 和 Y 分量,并且目前总共有五个变量。以下示例创建了一个新变量 (V6),它是速度分量的大小。

V6 = (V3*V3+V4*V4)**0.5

或者

V6 = sqrt(V3**2+V4**2)

6、您还可以使用以下等式完成上述操作(假设您已经定义了活动框架的矢量分量):

{Mag} = sqrt(U*U + V*V)

7、以下等式将名为 diff 的变量的值设置为从现有深度值中减去名为 depth 的变量的截断值:

{diff} = {depth} – trunc({depth})

8、在下一个等式中,假设 C 和 S 都已定义,则 C(轮廓变量)设置为 S(散布大小变量)的绝对值:

C = abs(S)

8、在以下示例中,创建了一个新变量(假设数据文件中最初只存在七个变量)。 V8(新变量)的值是根据现有变量的函数计算得出的:

V8 = SQRT((V1*V1+V2*V2+V3*V3)/(287.0*V4*V6))

9、上述操作可以通过两个更简单的步骤执行,如下所示:

V8 = V1*V1+V2*V2+V3*V3
V8 = SQRT(V8/(287.0*V4*V6))

10、以下等式将名为 TIME 的任何低于 5.0 的变量值替换为 5.0。换句话说,TIME 的值被替换为 TIME 的当前值和 5.0 中的最大值:

{TIME} = max({TIME},5)

11、以下等式创建变量 V4,其在 X<0 的点处具有 X 值;在其他点,它的值为零(这不会影响 X 的任何值):

V4 = min(X,0)

12、另一个使用内部函数的示例如下所示:

V8 = 55.0*SIN(V3*3.14/180.0) + ALOG(V4**3/(v1+1.0))

13、您还可以在方程式中引用 I、J 和 K 指数。例如,如果您想使用值消隐来切出一个区域的一部分,您可以创建一个新变量,它是 I 和 J 指数的函数(用于 IJ 有序数据)。然后,通过使用值消隐,您可以删除值消隐变量的值小于或等于值消隐截止值的某些单元格。

14、下面是一个计算值消隐变量的示例,该变量在 I=10 到 30 的单元格块中为零,在其他单元格中等于 1:

V3 = min((max(I,30)-min(I,10)-20),1)

15、以下等式将 Y 的所有值替换为 Y 的当前值与区域 1 中的 Y 值之间的差值。(如果区域 1 用于数据更改,则 Y 的新值将在该区域中为零。)

Y = Y – Y[1]

16、以下等式将 V3 的值(在 IJ 有序区域中)替换为四个相邻数据点的 V3 值的平均值:

V3 = (V3(i+1,j)+V3(i-1,j)+V3(i,j+1)+V3(i,j-1))/4

17、以下等式将名为 TEMP 的变量的值设置为在四个位置测量的名为 T 的变量的值的乘积:在区域 1 中,在当前数据点之前的两个索引值处,在当前区域中,在绝对索引处 3 ,在当前数据点的区域 4 中,在当前数据点的当前区域中。

{TEMP} = {T}[1](i-2) * {T}(3) * {T}[4] * {T}

18、指数可以与区域规格相结合。首先列出区域,然后列出索引偏移量。例如:

V3 = V3 – V3[1](i+1)
Y = Y[1] – Y[2](1) + Y(1,j+3) + Y

19、通过字母代码引用变量:

V3 = I + J
V4 = cos(X) * cos(Y) * cos(Z)
{距离} = sqrt(U*U + V*V + W*W)
{温度} = min(B,1)

20、为给定变量指定区域编号:

V3 = V3 – V3[1]
X = ( X[1] + X[2] + X[3] ) / 3
{TempAdj} = {Temp}[7] – {Adj}
V8 = V1[19] – 2*C[21] + {R/T}[18]


Tecplot后处理扩展工具:自定义函数

免责声明

泵小丫收集资源均来自网络

仅用于内部交流之用,如有需要请支持正版

声明:本站所有文章,如无特殊说明或标注,均为本站原创发布。任何个人或组织,在未征得本站同意时,禁止复制、盗用、采集、发布本站内容到任何网站、书籍等各类媒体平台。如若本站内容侵犯了原著者的合法权益,可联系我们进行处理,邮箱admin@7b3.cn。

给TA充电
共{{data.count}}人
人已充电
仿真模拟

Fluent meshing中如何查看网格质量分布

2022-7-11 22:56:32

仿真模拟

Cfturbo导入TurboGrid的三种方法

2022-10-5 14:41:58

2 条回复 A文章作者 M管理员
  1. 子丘

    以前就会求个速度,了解了

个人中心
购物车
优惠劵
今日签到
有新私信 私信列表
搜索