与外部数据交互: table命令与table FISH

1年前浏览990

1 引言

table本质上是一个(x,y)值的集合，其中这些值对可以按任何顺序排列，并通过索引进行查找。经常用到的一些场景是通过table储存一个复杂的边坡面，或者通过table指定一个本构模型中的参数。table的操作是双向的，也就是说，既可以从一个table中读取数值，也可以给table插入新值。本笔记简单描述了table的用法。

2 table命令

table命令能使用于所有的软件。这是目前版本中包含的命令。

对应地，这些命令也可以在FISH中使用，事实上在FISH中使用table命令的情形更多一些。

3 命令解释

下面简要介绍table的使用方法。定义一个FISH函数granulometry，在这个函数中首先创建一个表'mwu', 取得这个表的指针exptab，然后赋值6个所谓的(x,y)值对。

def granulometry

table.delete('mwu')

global exptab = table.create('mwu')

table(exptab,0.00236) = 0.01

table(exptab,0.005) = 0.152

table(exptab,0.01) = 0.902

table(exptab,0.014) = 0.989

table(exptab,0.02) = 0.991

table(exptab, 0.028) = 1.0

end

@granulometry

在这段程序中，我使用了 table.delete('mwu') 这句代码首先删除掉表，然后再建立新表，如果不使用这句代码，第二次运行这个FISH函数就会出错，因为表'mwu'已经存在，这一点和数据库语言SQL的用法一样。在这里我只是为了演示函数的功能，如果在代码开始使用model new 则不存在这个问题，因为model new清除了内存中的所有东西。

(1) 使用 table list命令可以列出table中x,y的最大值和最小值，但是这些值好像不能直接提取出来，还没试验python模块，不知道使用python能不能把最大值和最小值提取出来。

(2) 可以根据索引值取出table中的x,y 值，但有一个问题不知道如何解决，即如何获取索引值，比如一个数据集合中有1000个数据，想要获取其中的某个值，如何才能知道这个值的索引值，python语言在循环时可以自动建立索引值idx。

[x1 = table.x(exptab,6)]

[y1 = table.y(exptab,6)]

print @x1, @y1

(3) table具有插值功能，这意味着我们可以在x的取值范围内给出任意一个x值，从而可求出对应的y值。例如，

[y = table(exptab, 0.025)]

print @y

结果为y=0.996625, table使用了线性插值，这在数据点比较密集的情况下能够满足一定的精度。

(4) 可以在一个FISH内创建多个table, 例如：

fish define yield_surface

m = 1.28128

mmc = 1.2

s1 = table.create('mc')

table.x(s1,1) = 0.0

table.y(s1,1) = 0.0

table.x(s1,2) = 80

table.y(s1,2) = 80*mmc

s2 = table.create('ellipse')

loop local i (1,2001)

x = 0.05*(i-1)

table.x(s2,i) = x

table.y(s2,i) = m*50.0*math.sqrt(1.0 - (x/50.0-1.0)*(x/50.0-1.0))

endloop

end

@yield_surface

3 应用举例

PFC中有一个生成Particle-Size Distribution Curve(PSD)的例子，为PSD计算提供了思路，根据上面所述的方法可以计算出Effective Size D10, D30和D60, 从而可以计算出Uniformity Coefficient (Cu)和Coefficient of Gradation (Cc)。不过应当注意，真实的PSD曲线筛分尺寸是从大到小排列的，与演示例子正好相反，而且颗粒尺寸横坐标使用的是对数坐标轴log(d)。