24小时热门版块排行榜    

查看: 1344  |  回复: 6

天航空航

新虫 (初入文坛)

[求助] matlab插值

大家好,我用fluent得到了圆柱面的压力分布,知道每个节点坐标及压力值。然后想在patran里进行结构分析,但是两个软件划分的网格不一致,需要使用matlab进行差值。知道了fluent节点的(x,y,z,p)和patran里节点的(x,y,z),怎么插值好点。我是使用了griddata进行差值,结果出现很多NaN。很头疼,着急啊。而且我是用3000个节点来差值40000个节点,感觉好难啊!求各位帮忙,谢谢
A=xlsread('c:\Users\yu\Desktop\fluent.xls');
B=xlsread('c:\Users\yu\Desktop\patran.xls');
x=A(:,1);y=A(:,2);z=A(:,3);v=A(:,4);
xx=B(:,1);yy=B(:,2);zz=B(:,3);
vv=griddata(x,y,z,v,xx,yy,zz)
回复此楼
阿斯顿所得税大声道
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bucheron

金虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
1592203609: 金币+1, 谢谢回帖 2013-12-06 21:29:53
怀疑是patran 网格坐标问题, 尝试 在matlab 中 用 meshgrid 生成你的patran 坐标网格节点,然后重新 用 griddata
Born to suffer !Born to burn!
2楼2013-12-05 23:47:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

言西早幻城

新虫 (正式写手)

内容已删除
阅读开阔眼界,交流助力成功
3楼2013-12-06 07:46:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunshine819

金虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
1592203609: 金币+1, 谢谢回帖 2013-12-06 21:30:04
一般CFD的网格密度会远远大于CSD的,楼主的情况恰好相反,这样的问题定义是否合理楼主考虑一下吧。
日光之下皆覆辙 月光之下皆旧梦
4楼2013-12-06 11:14:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天航空航

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by sunshine819 at 2013-12-06 11:14:07
一般CFD的网格密度会远远大于CSD的,楼主的情况恰好相反,这样的问题定义是否合理楼主考虑一下吧。

假如我fluent计算飞行器表面的压强,要对整个庞大计算域进行网格划分,飞行器表面网格虽然细化了,但patran只对飞行器划分网格,表面网格比我多很正常啊
阿斯顿所得税大声道
5楼2013-12-06 16:24:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天航空航

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by bucheron at 2013-12-05 23:47:12
怀疑是patran 网格坐标问题, 尝试 在matlab 中 用 meshgrid 生成你的patran 坐标网格节点,然后重新 用 griddata

我的节点全部在圆柱面上,用meshgrid的话,得到的节点还在这个面上吗
阿斯顿所得税大声道
6楼2013-12-06 16:25:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xingangahu

木虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
出现NaN 可能是由于matlab中不支持外插值造成的 插值务必保证加密网格[xmin1 xmax1]只能在初始网格[xmin xmax]内部
7楼2013-12-07 09:09:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 天航空航 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 有没有仍没收到信息的 +3 德尚中行 2026-08-27 3/150 2026-08-28 15:24 by rongshuxia
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +4 gdd2018 2026-08-28 9/450 2026-08-28 14:22 by jklily
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见