24小时热门版块排行榜    

查看: 1338  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 天航空航 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-28 05:13 by winsaint
[基金申请] 基金系统什么内容也没有 30+3 winsaint 2026-08-27 6/300 2026-08-28 01:46 by alongwaytogo
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +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
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见