确定外倾角和厚度翼型- xy坐标系

153(30天)
亲爱的所有,
我有翼型坐标,一个封闭多边形的形式(xy点下面)。
我需要一个函数来确定最大可能内切圆多边形中的每个点(顶点)。
ps:内切圆顶点点被定义为最大可能的内圈画多边形顶点的切线。
机翼下面给出包含33分,我需要33镌刻圈……
否则需要确定的外倾角和厚度分布在下面的评论中提到的机翼。
任何帮助都是高度赞赏。
提前谢谢,
K Vijay Anand
%机翼的坐标
xy = [1.00000 0.94908 0.89816 0.79638 0.69491 0.59377 0.49322 0.39328 0.29389 0.19489 0.14562 0.09650 0.07203 0.04764 0.02338 0.01143 0.00000 0.01330 0.02596 0.05095 0.07573 0.10048 0.14992 0.19932 0.29821 0.39752 0.49722 0.59742 0.69804 0.79883 0.89957 0.94979 1.00000;
0.00000 0.02835 0.05669 0.11138 0.15658 0.19180 0.20853 0.20678 0.18805 0.15733 0.13473 0.10764 0.09134 0.07255 0.04976 0.03287 0.00000 -0.02457 -0.02966 -0.02934 -0.02254 -0.01473 0.00237 0.02098 0.05519 0.07642 0.08566 0.07942 0.06019 0.03596 0.01324 0.00637 - 0.00000);
情节(xy (1:), xy(2:),“啊——”,“线宽”,2);网格;轴相等;
5个评论

登录置评。

答案(2)

为副总经理
为副总经理 2021年9月30日
这并不是一个很好的方法,但我更比polyshape用于滥用图像处理工具()等。泰森多边形法图的提到让我认为距离地图将是一个不错的起点。
profileline = [1.00000 0.94908 0.89816 0.79638 0.69491 0.59377 0.49322 0.39328 0.29389 0.19489 0.14562 0.09650 0.07203 0.04764 0.02338 0.01143 0.00000 0.01330 0.02596 0.05095 0.07573 0.10048 0.14992 0.19932 0.29821 0.39752 0.49722 0.59742 0.69804 0.79883 0.89957 0.94979 1.00000;
0.00000 0.02835 0.05669 0.11138 0.15658 0.19180 0.20853 0.20678 0.18805 0.15733 0.13473 0.10764 0.09134 0.07255 0.04976 0.03287 0.00000 -0.02457 -0.02966 -0.02934 -0.02254 -0.01473 0.00237 0.02098 0.05519 0.07642 0.08566 0.07942 0.06019 0.03596 0.01324 0.00637 - 0.00000);
stepsize = 0.0001;%有效决议
xrange = (-0.1 - 1.1);
yrange = (-0.1 - 0.25);
%使用一个虚拟的图像显示生成一个图像的概要文件
x = xrange (1): stepsize: xrange (2);
y = yrange (1): stepsize: yrange (2);
h =图像(yrange, xrange(元素个数(y),元素个数(x)));
L = images.roi.Polygon (gca);
l位置= profileline。”;
掩码= ~ createMask (L);
%找到camberline坐标,从山脊线距离的地图
%,这只是一个不完整的解决方案由于前缘角
深度贴图= bwdist(面具);
[~,idx] = max(深度贴图,[],1);
camberline1 = [x;y (idx)];
%使用初始估计找到平分这个概要文件
%,这样可以估计前沿
[~,断点]= max (camberline1 (2:));
断点=圆(断点/ 2);%避免肤浅的斜坡
[~,idx] = max(深度贴图(:,1:断点),[],2);
camberline2 = [x (idx);y);
%情节估计演示
次要情节(2,1,1)
情节(profileline (1:), profileline (2:));持有;网格
情节(camberline1 (1:), camberline1 (2:)“:”);
情节(camberline2 (1:), camberline2 (2:)”乙:“);
xlim xrange ()
ylim (yrange)
%清理和合并估计
:camberline1 (:, camberline1 (2) = = yrange(1)) =南;
:camberline2 (:, camberline2 (1) = = xrange(1)) =南;
:camberline2 (:, camberline2 (1) = = x(断点))=南;
camberline = [camberline2 camberline1(:,断点+ 1:结束);
%的阴谋合并估计
次要情节(2,1,2)
情节(profileline (1:), profileline (2:));持有;网格
情节(camberline (1:), camberline (2:)“b——”);
xlim xrange ()
ylim (yrange)
如前所述,前缘区域是有问题的。我怀疑要么估计是正确的在这一地区,但第二个看上去略微更合理,我想至少证明一个分段的解决方案可能是一个选项。它会花费我更多的时间使我的大脑更健壮的方法。
5个评论
Vijay Anand
Vijay Anand 2021年10月7日
编辑:Vijay Anand 2021年10月7日
亲爱的先生。 @DGM 和其他人,
一个小小的突破……
我之前认为,这些圈不需要上圈……
弧线的定义说,应该正常厚度曲面线(点)…厚度也等于在顶部和底部的弧线。
请参阅图1、2、3 图1
1
虚线竖线是厚度线绘制verticaly在轴。
倾斜的线厚度局部曲面线绘制正常点
圆圈是与曲率半径和中心厚度/ 2点
图2
图2
图3
图3
从前缘附近的数据(图2),很明显,圆圈更大。
请给你的有价值的反馈在这个…
不过,我与算法来确定弧线和厚度分布。
ps:这些数字我从已知曲面生成和厚度分布方程(NACA 4312 Aifoil)和坐标点
提前谢谢,
K Vijay Anand

登录置评。


布鲁诺陈德良
布鲁诺陈德良 2021年10月4日
如果你能忍受一些有限的精度,你开启离散化室内为黑白图像,寻找骨架
每个点的骨架(圆心)半径最小的距离边界(不是很难计算)。

类别

找到更多的在机翼的工具帮助中心文件交换

s manbetx 845


释放

R2020a

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!