24小时热门版块排行榜

返回列表

类人之子

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 5
帖子: 4
在线: 1.5小时
虫号: 3910283
注册: 2015-06-02
专业: 河流海岸动力学与泥沙研究

[求助] 【fortran求助】如何将三角形网格节点上数据赋值给四边形网格节点已有1人参与

现有同一区域范围内两种网格（三角形网格和四边形网格），有三角形网格“三角形网格编号.txt”、“三角形网格各节点坐标.txt”、“三角形网格各节点分流速.txt”三个文件，四边形网格有文件“四边形网格坐标.txt”（三角形网格数32465个，每个三角形6个节点（三个顶点和三条边中点），因相邻三角形共用三个节点，共有节点数65302个，各文件见附件），为简化程序，四边形网格只给了1000个节点及其坐标。现需要将三角形网格上的X、Y方向上的分流速赋值到四边形节点坐标上。三角形网格与四边形网格形式见图片。
   赋值计算如下：M点分流速Um=（sABM*Uc+sACM*Ub+sBCM*Ua）/sABC，同理M点分流速Vm=（sABM*Vc+sACM*Vb+sBCM*Va）/sABC。也就是说三个顶点的值乘以对面小三角形面积占大三角形面积权重的和值。（U、V为节点分流速；sABC等为三角形面积）
   需要说明的是：三角形相邻网格不一定是连续编号，而且每个三角形网格的节点编号不是连续的。
程序的总体思路是：对于每一个四边形网格节点，搜索这个节点在哪一个编号的三角形内（也有可能这个节点刚好落在三角形某条边上，如果是这种情况，则直接根据直线插值），然后再根据上面的公式计算其分流速U和V。
请各位大神帮帮忙，谢谢！
  最后生成结果如下形式：
四边形节点编号          X坐标          Y坐标       分流速U       分流速V
   1                            X1                Y1             U1                V1
   2                            X2                Y2             U2                V2
   3                            X3                Y3             U3                V3
。
。
。
1000                      X1000          Y1000       U1000       V1000

三、四变形.jpg