24小时热门版块排行榜    

查看: 1099  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 天航空航 的主题更新
信息提示
请填处理意见