向高手求助动网格模拟平板阶跃响应的UDF问题

ranjhzhl

物理问题：为了计算一块平板在阶跃响应下的气动力

已知： 1.平板处于有位移状态下的CFD网格坐标,数据按J=1,JN...I=1,IN;...END;END的格式存入mode1.dat文件中;
       2.平板的初始状态只能用离散点描述,但确保Fluent导入的cas网格文件中物面每点均在model1.dat中有点与之对应；
　　　 3.cas网格文件中重叠面的网格(若是自编程，常认为是ZETA=1的面)在mode1.dat中也有点与之对应;
       4.平板运动规律是：从初始位移状态随着fluent求解器的非定常推进计算，平板回归到水平位置;

本人拟定的思路：
       1.用Gridgen生成平板处于平衡位置(水平状态)时的3D网格，分好边界类型并将Cas文件导入Fluent；
       2.在求解器推进计算的初始瞬间(以Time为标志)，读入mode1.dat文件内的数据；
       3.在3D网格中将mode1.dat内x/y值均相同的点处，模态z坐标＋3D网格系的z坐标,以获得阶跃运动的初始状态;
       4.将非定常推进计算过程中，x/y值不变，z值逐渐回复到水平位置的平衡状态；
       5.网格运动采用DEFINE_GRID_MOTION这个函数。

根据上述想法，我编写了一个UDF文件，但总是运行失败,目前还有如下问题没搞明白，觉得比较致命：

       1.3D网格在自编程序中有严格的I,J,K方向及编号，用Tecplot显示即可，但在Fluent里面却不知怎么对应；
       2.Gridgen生成的网格导入Fluent后，除了物面及重叠面上有与mode1.dat中x/y坐标相同的点外，空间区域的点x/y不一定相同；
       3.点击mesh motion后，网格没反应，似乎是读网格就出错了。

       4.请帮我看看我的UDF程序错在哪里？该怎么改或者说怎么做是解决该问题的最好方法？

我为这个问题已经花了尽两个月的时间了，到目前为止还没搞定，特请高手们指点一下迷津，非常感谢，我的qq是61410210,Email是ranjhzhl@163.com