| 查看: 1195 | 回复: 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
» 猜你喜欢
08工科275分求调剂
已经有13人回复
-
材料工程322分
已经有6人回复
-
一志愿北京科技大学085601材料工程英一数二初试总分335求调剂
已经有7人回复
-
274求调剂
已经有5人回复
-
【求调剂】新能源材料本科，一志愿211，初试321
已经有3人回复
-
材料考研调剂
已经有3人回复
-
材料调剂
已经有12人回复
-
英一数一408，总分284，二战真诚求调剂
已经有14人回复
-
085410 一志愿211 22408分数359求调剂
已经有4人回复
-
271求调剂
已经有19人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有限元钢架模型怎么画
已经有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
