| 查看: 1219 | 回复: 2 | ||
[求助]
【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  000000.jpg |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 三角形网格编号.txt
- 附件 2 : 三角形网格各节点分流速.txt
- 附件 3 : 三角形网格各节点坐标.txt
- 附件 4 : 四边形网格坐标.txt
2015-06-05 17:11:39, 1.29 M
2015-06-05 17:11:40, 1.23 M
2015-06-05 17:11:42, 2.11 M
2015-06-05 17:11:42, 25.08 K
» 猜你喜欢
今年审到国自然15份，谈谈感受
已经有15人回复
-
国自然上会要求
已经有5人回复
-
面上本子正文33页，违规吗？会被低分嘛？
已经有7人回复
-
上海大学实验技术岗位非升即走
已经有8人回复
-
考博自荐
已经有6人回复
-
青C资助名额大幅增加！
已经有16人回复
-
重磅！青年科学基金项目（C类）资助增幅预计超过50%
已经有10人回复
-
我在等一个没有答案的答案
已经有3人回复
-
半夜喝咖啡
已经有5人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有限元钢架模型怎么画
已经有6人回复
- 关于网格数量与节点数量的讨论
已经有11人回复
- 在用差分法插值时边界怎么处理
已经有6人回复
- fluent 仿真温度出现负值
已经有9人回复
- 一模型用gambit能划三角形网格但为什么划不了四边形网格啊？？
已经有11人回复
- 求助：如何让fluent在服务器多个节点上进行运算
已经有5人回复
- 请教Ansys 车轮 模型车轮原地转动，如何做到？
已经有11人回复
- ABAQUS里面约束参考点问题
已经有5人回复
- QUICK格式应该如何使用？
已经有8人回复
- VOF模型中，体积分数云图中液体和气体之间总存在一个边界层，如何去掉？
已经有14人回复
- 使用AdNDP方法以及ELF/LOL、多中心键级研究多中心键
已经有14人回复
- Fluent 结构化网格与非结构化网格
已经有7人回复
- 集群上面的节点需要定时间的重启吗？
已经有18人回复
 |
2楼2015-06-08 09:07:11
3楼2015-11-16 12:38:00
20