如何根据位移函数画出圆管的振动模态？

圆管示意图：

位移函数形式是切比雪夫多项式

现在已经在Matlab中用 [V，D] = eig（K，M）求出了特征频率和特征向量，特征向量就是的值，但位移函数的自变量有三个(r,z,θ)，r的范围是[0.5,1]，z的范围是[0,1]，θ的范围是[0,2*pi]，这样把自变量的取值带入位移函数后u，v，w分别都是N1×N2×N3的数组，这该如何根据三个位移分量u，v，w画出圆管的振动模态图