| 查看: 1275 | 回复: 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
» 猜你喜欢
估计是周四
已经有3人回复
-
基金啊基金
已经有3人回复
-
科研孤儿太难了
已经有15人回复
-
让我中一个面上吧！
已经有10人回复
-
言之凿凿，明天官网会通知放榜？
已经有17人回复
-
时间戳又变了
已经有17人回复
-
今天系统多次维护，明天很可能放榜！
已经有11人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
放榜前的不淡定
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有限元钢架模型怎么画
已经有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
40